CN102184414A - 泵示功图的识别和判断方法及其系统 - Google Patents

Abstract

公开了一种泵示功图的识别和判断方法及其系统，所述方法包括：对泵示功图进行预处理，提取泵示功图信息；根据所述泵示功图信息对泵示功图进行归一化处理；对归一化处理后的泵示功图进行多边形逼近处理；利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别。所述系统包括泵示功图预处理模块、归一化处理模块、多边形逼近处理模块及故障识别模块。根据本发明提供的泵示功图的识别和判断方法及其系统不仅能提高深井泵示功图识别准确率，而且对油井工况判断精确。

Description

泵示功图的识别和判断方法及其系统

技术领域

本发明涉及油田采油技术领域，特别涉及一种泵示功图的识别和判断方法及其系统。

背景技术

抽油机井井下泵示功图是诊断抽油系统工作状态的依据，国内外现有油井诊断系统大部分采用这一依据。诊断的关键是井下泵示功图识别的准确性，识别越精确，诊断结论越正确。功图识别的关键在于如何提取其最具代表性的特征。目前，泵示功图特征的提取，大多采用的是专家相面法、面积法、模板匹配法网格法、矢量法等。这些方法注重提取泵示功图的特征，如几何、形状、力学等特征，普遍存在特征参数偏多的不足。所使用的判别器多为距离分类器，这类分类器实质是满足一定条件的特殊贝叶斯分类器，要求各特征分量线性无关、各特征分量的方差相同、类条件概率密度函数为正态分布。虽然，分类器在某种意义上是最优分类器，但针对采油系统工况分类问题，不能满足这一前提条件，导致识别结果误判。

发明内容

本发明的目的之一是提供一种提高泵示功图识别的准确性，实现泵工况准确判断的一种泵示功图的识别和判断方法及其系统。

根据本发明的一个方面，提供一种泵示功图的识别和判断方法包括：

对泵示功图进行预处理，提取泵示功图信息；

根据所述泵示功图信息对泵示功图进行归一化处理；

对归一化处理后的泵示功图进行多边形逼近处理；

利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别。

根据本发明的另一个方面，提供一种泵示功图的识别和判断系统包括：

泵示功图预处理模块，对泵示功图进行预处理，提取泵示功图信息；

归一化处理模块，根据提取的泵示功图信息对泵示功图进行归一化处理；

多边形逼近处理模块，对归一化处理后的泵示功图进行多边形逼近处理；

故障识别模块，将所述泵示功图的矢量链与标准故障矢量链库中矢量链比较，进行故障判别；所述标准故障矢量链库包括泵各种故障的矢量链。

根据本发明提供的泵示功图的识别和判断方法及其系统不仅能提高深井泵示功图识别准确率，而且对油井工况判断精确。

附图说明

图1是本发明实施例提供的泵示功图的识别和判断方法的流程示意图；

图2是本发明实施例提供的泵示功图的识别和判断系统的流程示意图；

图3是本发明实施例提供的泵示功图的识别和判断系统的结构框图；

图4是图3所示系统中故障识别模块的结构框图；

图5是本发明实施例提供的对一泵示功图进行多边形逼近处理结果示意图；

图6是本发明实施例提供的对已进行多边形逼近处理的泵示功图未进行校正的示意图；

图7是对图6所提供的泵示功图进行校正后的示意图；

本发明目的、功能及优点将结合实施例，参照附图做进一步说明。

具体实施方式

如图1所示，本发明实施例提供的一种泵示功图的识别和判断方法包括：

步骤S1、对泵示功图进行预处理，提取泵示功图信息。该步骤包括：

步骤S11、生成泵示功图。根据有杆泵采油系统故障诊断模型(是现有技术，不再赘述)求解得到的泵示功图是以泵位移s_p为横坐标、以对应的泵载荷F_p为纵坐标的封闭曲线，描述的是泵载荷随泵位移的变化规律。泵示功图的求解结果是以两组离散数据对的形式给出：

(θ_i，s_pi)，i＝0，1，Λ，2N     (1.1)，

(θ_i，F_pi)，i＝0，1，Λ，2N     (1.2)，

式中的θ_i为曲柄转角θ＝θ(t)在一个周期区间[0，2π]上的2N+1个等分点：

式中的N需要权衡求解精度和求解速度而定，例如，30＜N＜300。

步骤S12、对泵示功图进行函数化处理。泵示功图可以用函数F_p＝F_p(s)(1.4)来逼近，这是一个多值函数，称之为泵示功图函数。泵示功图函数可由下面的参数方程来确定：

其中函数s_p＝s_p(θ)称为泵位移-曲柄转角函数，函数F_p＝f_p(θ)称为泵载荷-曲柄转角函数，它们都是曲柄转角θ的单值函数，可以分别用离散数据对式(1.1)和(1.2)来进行Fourier级数逼近：

$s_p=s_p\left(\mathrm{\θ}\right)$

$=\frac{1}{2}{a}_0+\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}(a_k\cos \mathrm{k\θ}+b_k\sin \mathrm{k\θ})$

$\≈\frac{1}{2}{a}_0+\underset{k=1}{\overset{N}{\mathrm{\Σ}}}(a_k\cos \mathrm{k\θ}+b_k\sin \mathrm{k\θ})---\left(1.6\right),$

$F_p=f_p\left(\mathrm{\θ}\right)$

$=\frac{1}{2}{d}_0+\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}(d_k\cos \mathrm{k\θ}+e_k\sin \mathrm{k\θ})$

$\≈\frac{1}{2}{d}_0+\underset{k=1}{\overset{N}{\mathrm{\Σ}}}(d_k\cos \mathrm{k\θ}+e_k\sin \mathrm{k\θ})---(1.7),$

令

{s_p}＝(s_p0 s_p1 s_p2 Λ s_p，2N-1 s_p，2N)^T   (1.8)，

{F_p}＝(F_p0 F_p1 F_p2 Λ F_p，2N-1 F_p，2N)^T   (1.9)，

{s}＝(a₀ a₁ b₁ Λ a_N b_N)^T    (1.10)，

{d}＝(d₀ d₁ e₁ Λ d_N E_N)^T    (1.11)，

$\left[T_F\right]=\left(\begin{array}{cccccc}\frac{1}{2}& \cos {\mathrm{\θ}}_0& \sin {\mathrm{\θ}}_0& \mathrm{\Λ}& \cos N{\mathrm{\θ}}_0& \sin N{\mathrm{\θ}}_0\\ \frac{1}{2}& \cos {\mathrm{\θ}}_1& \sin {\mathrm{\θ}}_1& \mathrm{\Λ}& \cos N{\mathrm{\θ}}_1& \sin N{\mathrm{\θ}}_1\\ \frac{1}{2}& \cos {\mathrm{\θ}}_2& \sin {\mathrm{\θ}}_2& \mathrm{\Λ}& \cos N{\mathrm{\θ}}_2& \sin {\mathrm{N\θ}}_2\\ \mathrm{\Λ}& \mathrm{\Λ}& \mathrm{\Λ}& \mathrm{\Λ}& \mathrm{\Λ}& \mathrm{\Λ}\\ \frac{1}{2}& \cos {\mathrm{\θ}}_{2N-1}& \sin {\mathrm{\θ}}_{2N-1}& \mathrm{\Λ}& \cos N{\mathrm{\θ}}_{2N-1}& \sin N{\mathrm{\θ}}_{2N-1}\\ \frac{1}{2}& \cos {\mathrm{\θ}}_{2N}& \sin {\mathrm{\θ}}_{2N}& \mathrm{\Λ}& \cos N{\mathrm{\θ}}_{2N}& \sin N{\mathrm{\θ}}_{2N}\end{array}\right)---\left(1.12\right),$

则有{s_p}＝[T_F]{a}(1.13)，

{F_p}＝[T_F]{d}(1.14)，

可得：

{d}＝[T_F]^-1{F_p}(1.15)，

{a}＝[T_F]^-1{s_p}(1.16)，

以式(1.6)和(1.7)为参数方程可以确定泵示功图函数

F_p＝F_p(s)。

步骤S13、对所述函数化处理后的泵示功图进行信息提取。该步骤包括：

步骤S131、提取泵冲程与最小、最大泵载荷信息。

泵最小位移： $s_{p\min }=\underset{i=0}{\overset{2N}{\min }}{s}_{\mathrm{pi}}---\left(1.18\right),$

泵最大位移： $s_{p\min }=\underset{i=0}{\overset{2N}{\max }}{s}_{\mathrm{pi}}---\left(1.19\right),$

则泵冲程为： $S_p=s_{p_{\max }}-s_{p_{\min }}---\left(1.20\right),$

泵最小载荷为： $F_{p\min }=\underset{i=0}{\overset{2N}{\min }}{F}_{\mathrm{pi}}---\left(1.21\right),$

泵最大载荷为： $F_{p\max }=\underset{i=0}{\overset{2N}{\max }}{F}_{\mathrm{pi}}---(1.22).$

步骤S132、提取泵示功图面积与泵马力信息。

令(s_p，2N+1，F_p，2N+1)＝(s_p0，F_p0)，则泵示功图的面积为：

$A_p=\underset{i=1}{\overset{2N+1}{\mathrm{\Σ}}}\frac{1}{2}(F_{\mathrm{pi}}+F_{p,i-1})(s_{\mathrm{pi}}-s_{p,i-1})$

$=\frac{1}{2}\underset{i=1}{\overset{2N+1}{\mathrm{\Σ}}}(F_{\mathrm{pi}}+F_{p,i-1})(s_{\mathrm{pi}}-s_{p,i-1})=\frac{1}{2}\underset{i=1}{\overset{2N+1}{\mathrm{\Σ}}}(F_{p,i-1}{s}_{\mathrm{pi}}-F_{\mathrm{pi}}{s}_{p,i-1})---\left(1.2\right),$

令泵马力为 $P_p=\frac{A_{p}n}{60}---\left(1.24\right),$

步骤S133、提取泵示功图曲线上各点的形状特征；所述形状特征包括斜率、凹凸性及曲率。

其中，斜率(一阶导数)

凹凸性(二阶导数)

曲率

步骤S2、根据泵示功图信息对泵示功图进行归一化处理。

根据前面所述，泵示功图的求解结果以两组离散数据对式错误！未找到引用源。和(1.2)的形式给出，可以把它们化简，用一组离散数据对：(s_pi，F_pi)，i＝0，1，Λ，2N(1.28)来表示泵示功图，每一个离散数据对代表泵示功图上一点。其中，为了在泵示功图识别时消除泵示功图实际尺寸的影响，从而需要对泵示功图进行归一化处理，对式(1.28)进行如下变换：

${\stackrel{\‾}{s}}_{\mathrm{pi}}=\frac{s_{\mathrm{pi}}-s_{p\min }}{s_{p\max }-s_{p\min }}\×\mathrm{PPDYN}\_\mathrm{WIDTH}---\left(1.29\right),$

${\stackrel{\‾}{F}}_{\mathrm{pi}}=\frac{F_{\mathrm{pi}}-F_{p\min }}{F_{p\max }-F_{p\min }}\×\mathrm{PPDYN}\_\mathrm{HEIGHT}---(1.30),$

其中PPDYN_WIDTH、PPDYN_HEIGHT是出于绘图和识别的考虑而引入的两个放大系数。其中，0＜PPDYN_WIDTH＜500；0＜PPDYN_HEIGHT＜300；于是离散数据对式(1.28)就转换为归一化的离散数据对：

$({\stackrel{\‾}{s}}_{\mathrm{pi}},{\stackrel{\‾}{F}}_{\mathrm{pi}}),i=\mathrm{0,1},\mathrm{\Λ},2N---\left(1.31\right),$

步骤S3、对归一化处理后的泵示功图进行多边形逼近处理。

式(1.31)所表示的离散数据对对应于泵示功图上的各个点，在这些点中，有些点对泵示功图的形状起着关键作用，将这些点成为关键点；有些点则对泵示功图的形状无关紧要，将这些点称之为次要点。采用多边形逼近处理方法可找出泵示功图上的所有关键点，去除其余的次要点，并用关键点所组成的多边形来逼近泵示功图。多边形逼近处理方法包括：

若泵示功图上的两个数据点很接近，则删除其中一个；例如，当所述归一化处理后的泵示功图上的两个数据点同时满足横坐标之差小于0.05×PPDYN_WIDTH，纵坐标之差小于0.03×PPDYN_HEIGHT，则删除其中一个数据点。

若泵示功图上的几个连续数据点几乎在同一条直线上，则只保留其中第一个数据点和最后一个数据点。例如，当归一化处理后的泵示功图上的相邻数据点所构成的线段的斜率不大于0.03时(即相邻两条线段的斜率值不大于0.03)，此时，相邻两条线段可视为一条线段。则只保留第一条线段的第一个数据点和第二条线段的最后一个数据点。

将泵示功图上的所有数据点连接成封闭曲线，若相邻两条线段的夹角为接近夹角180°时，则去掉它们相交的数据点，平滑处理；例如，夹角为178°、179°等。图5是用多边形逼近法对某个示功图的处理结果。浅色线部分为计算得到的泵示功图(即未经过多边形逼近处理的泵示功图)，而黑色的多边形为经过多边形逼近处理的泵示功图，由此给出了该泵示功图多边形逼近后的数字特征(参见表1所示的多边形逼近数据表)。

表1

另外，用仪器测得的光杆示功图数据(泵示功图就是利用光杆示功图计算得到的，光杆示功图是利用油井安装的载荷、位移传感器采集的载荷、位移值绘制而成)通常有垂直误差和倾斜误差，有时这种误差会很大。根据此光杆示功图数据计算出的泵示功图数据也会存在垂直误差和倾斜误差，这将会影响泵示功图识别的准确性，有必要在泵示功图识别之前对其进行倾斜矫正。对泵示功图进行归一化处理和多边形逼近后，可在泵示功图上寻找这样的线段：

在多边形逼近处理后的泵示功图上搜索0.5×PPDYN_WIDTH＜PPDYN_WIDTH、0°＜和水平的夹角＜8°且取线段中点纵坐标为最大或最小值的线段；满足这三个条件的线段应该是上下冲程的静载线，它们应该是水平的。

根据搜索的线段对多边形逼近处理后的泵示功图进行倾斜矫正。图6所示的泵示功图经过倾斜矫正可得到如图7所示的泵示功图。

步骤S4、利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别。用多边形逼近法处理过的示功图特征明显，接下来就可以很容易地使用矢量特征法进行泵有效冲程的识别。每种故障的特征均可以用一个矢量链来表征，因此可建立一个标准故障矢量链库，这个库里包含了泵各种故障的矢量链。如图2所示，该步骤包括：

步骤S41、对多边形逼近处理后的泵示功图进行无因次化处理(标准化处理)、建立该泵示功图的矢量链；

步骤S42、将泵示功图的矢量链与标准故障矢量链库中矢量链比较，进行故障判别；标准故障矢量链库包括泵各种故障的矢量链。

如图3所示，基于图1所示的方法本发明实施例还提供一种泵示功图的识别和判断系统包括泵示功图预处理模块1、归一化处理模块2、多边形逼近处理模块3及故障识别模块4。

泵示功图预处理模块1对泵示功图进行预处理，提取泵示功图信息。归一化处理模块2根据提取的泵示功图信息对泵示功图进行归一化处理。多边形逼近处理模块3对归一化处理后的泵示功图进行多边形逼近处理。故障识别模块4利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别。

泵示功图预处理模块1包括泵示功图生成单元、函数化处理单元及信息提取单元。泵示功图生成单元生成泵示功图。函数化处理单元对生成的泵示功图进行函数化处理。信息提取单元对所述函数化处理后的泵示功图进行信息提取；所提取的信息包括泵冲程与最小、最大泵载荷信息，泵示功图面积与泵马力信息及泵示功图曲线上各点的形状特征。

多边形逼近处理模块3包括图形处理单元、倾斜矫正单元。图形处理单元当所述归一化处理后的泵示功图上的两个数据点同时满足横坐标之差小于0.05×PPDYN_WIDTH，纵坐标之差小于0.03×PPDYN_HEIGHT，则删除其中一个数据点；当所述归一化处理后的泵示功图上的相邻的数据点所构成的线段的斜率不大于0.03时，则只保留其中第一个数据点和最后一个数据点；将所述归一化处理后的泵示功图上的所有数据点连接成封闭曲线，当相邻两条线段的夹角接近夹角180°时，则去掉两条线段相交的数据点。倾斜矫正单元，在多边形逼近处理后的泵示功图上搜索0.5×PPDYN_WIDTH＜WIDTH的线段或者0°＜和水平的夹角＜8°的线段或者取线段中点纵坐标为最大或最小值的线段，根据搜索的线段对多边形逼近处理后的泵示功图进行倾斜矫正。

参见图4、故障识别模块4包括矢量链建立单元41及故障判断单元。矢量链建立单元41对多边形逼近处理后的泵示功图进行无因次化处理、建立该泵示功图的矢量链。故障判断单元42将所述泵示功图的矢量链与标准故障矢量链库中矢量链比较，进行故障判别；标准故障矢量链库包括泵各种故障的矢量链。

本发明提供的一种泵示功图的识别和判断方法及其系统，可从根本上解决了深井泵功图识别准确率低，油井工况判断不精确的问题。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (13)

1.一种泵示功图的识别和判断方法，其特征在于，包括：

对泵示功图进行预处理，提取泵示功图信息；

根据所述泵示功图信息对泵示功图进行归一化处理；

对归一化处理后的泵示功图进行多边形逼近处理；

利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别。

2.根据权利要求1所述的方法，其特征在于，所述对泵示功图进行预处理，提取泵示功图信息包括：

生成泵示功图；

对所述泵示功图进行函数化处理；

对所述函数化处理后的泵示功图进行信息提取。

3.根据权利要求2所述的泵示功图的识别和判断方法，其特征在于，所述生成泵示功图是：

根据有杆泵采油系统故障诊断模型求解得到泵示功图；所述泵示功图的求解结果以两组离散数据对的形式给出：

(θ_i，s_pi)，i＝0，1，Λ，2N，

(θ_i，F_pi)，i＝0，1，Λ，2N，

${\mathrm{\θ}}_i=\frac{2\mathrm{\π}}{2N+1}(i+0.5),i=\mathrm{0,1},\mathrm{\Λ},2N,$

式中的θ_i为曲柄转角θ＝θ(t)在一个周期区间[0，2π]上的2N+1个等分点，s_p为泵位移，F_p为泵载荷。

4.根据权利要求1所述的方法，其特征在于，所述对所述泵示功图进行函数化处理是：

$s_p=s_p\left(\mathrm{\θ}\right)$

以 $=\frac{1}{2}{a}_0+\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}(a_k\cos \mathrm{k\θ}+b_k\sin \mathrm{k\θ}),$

$\≈\frac{1}{2}{a}_0+\underset{k=1}{\overset{N}{\mathrm{\Σ}}}(a_k\cos \mathrm{k\θ}+b_k\sin \mathrm{k\θ})$

$F_p=f_p\left(\mathrm{\θ}\right)$

$=\frac{1}{2}{d}_0+\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}(d_k\cos \mathrm{k\θ}+e_k\sin \mathrm{k\θ})$

$\≈\frac{1}{2}{d}_0+\underset{k=1}{\overset{N}{\mathrm{\Σ}}}(d_k\cos \mathrm{k\θ}+e_k\sin \mathrm{k\θ})$ 为参数方程对所述泵示功图进行函数化处理，从而确定泵示功图函数F_p＝F_p(s)。

5.根据权利要求2所述的方法，其特征在于，所述对所述函数化处理后的泵示功图进行信息提取包括：

提取泵冲程与最小、最大泵载荷信息；

提取泵示功图面积与泵马力信息；

提取泵示功图曲线上各点的形状特征；所述形状特征包括斜率、凹凸性及曲率。

6.根据权利要求3所述的方法，其特征在于，所述根据所述泵示功图信息对泵示功图进行归一化处理包括：

将所述(θ_i，s_pi)，i＝0，1，Λ，2N，(θ_i，F_pi)，i＝0，1，Λ，2N所表示的泵示功图进行简化，得到由(s_pi，F_pi)，i＝0，1，Λ，2N表示的泵示功图；

将所述(s_pi，F_pi)，i＝0，1，Λ，2N进行如下处理：，

从而得归一化的离散数据对

所述PPDYN_WIDTH、PPDYN_HEIGHT是两个放大系数，所述离散数据对对应于归一化处理后的泵示功图上的各个数据点。

7.根据权利要求6所述的方法，其特征在于，所述对归一化处理后的泵示功图进行多边形逼近处理包括：

当所述归一化处理后的泵示功图上的两个数据点同时满足横坐标之差小于0.05×PPDYN_WIDTH，纵坐标之差小于0.03×PPDYN_HEIGHT时，则删除其中一个数据点；

当所述归一化处理后的泵示功图上的相邻的数据点所构成的线段的斜率不大于0.03时，则只保留其中第一个数据点和最后一个数据点；

将所述归一化处理后的泵示功图上的所有数据点连接成封闭曲线，当相邻两条线段的夹角为接近夹角180°时，则去掉两条线段相交的数据点。

8.根据权利要求7所述的方法，其特征在于，所述利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别之前还包括：

对多边形逼近处理后的泵示功图的倾斜矫正的步骤，该步骤包括：

在多边形逼近处理后的泵示功图上搜索0.5×PPDYN_WIDTH＜PPDYN_WIDTH、0°＜和水平的夹角＜8°且取线段中点纵坐标为最大或最小值的线段；

根据搜索的线段对多边形逼近处理后的泵示功图进行倾斜矫正。

9.根据权利要求1-8任一项所述的方法，其特征在于，所述利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别：

对多边形逼近处理后的泵示功图进行无因次化处理、建立该泵示功图的矢量链；

将所述泵示功图的矢量链与标准故障矢量链库中矢量链比较，进行故障判别；所述标准故障矢量链库包括泵各种故障的矢量链。

10.一种泵示功图的识别和判断系统，其特征在于，包括：

泵示功图预处理模块，对泵示功图进行预处理，提取泵示功图信息；

归一化处理模块，根据提取的泵示功图信息对泵示功图进行归一化处理；

多边形逼近处理模块，对归一化处理后的泵示功图进行多边形逼近处理；

故障识别模块，利用矢量特征法对多边形逼近处理后的泵示功图中的故障进行识别。

11.根据权利要求10所述的系统，其特征在于，所述泵示功图预处理模块包括：

泵示功图生成单元，生成泵示功图；

函数化处理单元，对生成的泵示功图进行函数化处理；

信息提取单元，对所述函数化处理后的泵示功图进行信息提取；所提取的信息包括泵冲程与最小、最大泵载荷信息，泵示功图面积与泵马力信息及泵示功图曲线上各点的形状特征。

12.根据权利要求10所述的系统，其特征在于，所述多边形逼近处理模块包括：

图形处理单元，当所述归一化处理后的泵示功图上的两个数据点同时满足横坐标之差小于0.05×PPDYN_WIDTH，纵坐标之差小于0.03×PPDYN_HEIGHT，则删除其中一个数据点；当所述归一化处理后的泵示功图上的相邻的数据点所构成的线段的斜率不大于0.03时，则只保留其中第一个数据点和最后一个数据点；将所述归一化处理后的泵示功图上的所有数据点连接成封闭曲线，当相邻两条线段的夹角为接近夹角180°时，则去掉两条线段相交的数据点；

倾斜矫正单元，在多边形逼近处理后的泵示功图上搜索0.5×PPDYN_WIDTH＜PPDYN_WIDTH、0°＜和水平的夹角＜8°且取线段中点纵坐标为最大或最小值的线段，根据搜索的线段对多边形逼近处理后的泵示功图进行倾斜矫正。

13.根据权利要求10-12任一项所述的系统，其特征在于，所述故障识别模块包括：

矢量链建立单元，对多边形逼近处理后的泵示功图进行无因次化处理、建立该泵示功图的矢量链；

故障判断单元，将所述泵示功图的矢量链与标准故障矢量链库中矢量链比较，进行故障判别；所述标准故障矢量链库包括泵各种故障的矢量链。

Images