CN104481496B - 一种有杆泵抽油井故障诊断方法 - Google Patents
一种有杆泵抽油井故障诊断方法 Download PDFInfo
- Publication number
- CN104481496B CN104481496B CN201410473163.1A CN201410473163A CN104481496B CN 104481496 B CN104481496 B CN 104481496B CN 201410473163 A CN201410473163 A CN 201410473163A CN 104481496 B CN104481496 B CN 104481496B
- Authority
- CN
- China
- Prior art keywords
- num
- characteristic vector
- pump dynamometers
- chain code
- pump
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000003745 diagnosis Methods 0.000 title abstract description 11
- 239000013598 vector Substances 0.000 claims abstract description 155
- 239000002245 particle Substances 0.000 claims abstract description 50
- 230000003595 spectral effect Effects 0.000 claims abstract description 6
- 238000005086 pumping Methods 0.000 claims description 85
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000007599 discharging Methods 0.000 claims description 24
- 238000013480 data collection Methods 0.000 claims description 13
- 239000003129 oil well Substances 0.000 claims description 11
- 238000013139 quantization Methods 0.000 claims description 8
- 101100517651 Caenorhabditis elegans num-1 gene Proteins 0.000 claims 3
- 238000012549 training Methods 0.000 abstract description 7
- 238000010606 normalization Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract 2
- 238000013316 zoning Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 18
- 230000006870 function Effects 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/008—Monitoring of down-hole pump systems, e.g. for the detection of "pumped-off" conditions
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明涉及一种有杆泵抽油井故障诊断方法,获取已知和待诊断有杆泵抽油井地面示功图,将有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图,采用16方向链码对归一化处理后的泵示功图进行边界链码重画,将边界链码重画后的泵示功图进行区域划分,并提取各区域的特征向量,利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类,完成有杆泵抽油井故障诊断,本发明采用16方向链码的数量对各个泵示功图特征向量进行量化能够有效、细致地描述图形特征,有杆泵抽油井故障诊断是利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类,不依赖训练样本,可以提高诊断的性能。
Description
技术领域
本发明属于抽油井故障诊断领域,特别涉及一种有杆泵抽油井故障诊断方法。
背景技术
有杆泵抽油井泵示功图在油田生产中具有十分重要的作用,它是判断抽油系统井下工作状态的主要依据。泵示功图可以用来进行井下生产状况的故障诊断,也可以用来计算井口的出液量,还可以用来计算井下动液面的位置。传统生产中,一般依靠技术人员对采集的泵示功图进行人工分析,从而判断抽油系统的实际工作状况。这种人工方式工作周期长,受主观因素影响较大,越来越无法满足企业不断提高生产效率的需求。
随着计算机、通信、电子技术等领域的不断进步,利用计算机取代人力劳动方式已成为国内外油田生产的发展趋势。鉴于泵示功图在生产中的重要性,采用计算机进行计算和分析已成为该领域研究的重点和难点。利用计算机处理泵示功图,关键在于提取其显著的特征,识别越准确,计算越正确。目前,对于泵示功图的特征提取方法,主要有面积法、匹配法、网格法、矢量法等。但是这些方法都不能很好地结合实际生产过程的特点,从而使得到的特征向量不能准确、细致地反映不同过程的显著特征。
有杆泵抽油井的抽油泵一般工作在数千米的地下,生产环境恶劣,故障发生率很高,采用计算机实现诊断是很有意义的。现有的诊断方法主要以有监督学习为主,根据生产经验得到不同故障类型的训练样本集,然后由这些训练样本训练出不同故障类型的训练集模型,从而实现待诊断样本的诊断。这种有监督学习方法主要存在以下两方面问题:一是诊断的准确率很大程度上取决于所构成训练样本集的样本特征向量是否能够准确反映样本的典型特征;二是这种有监督方式主要依靠人为划分训练样本集,容易受到主观因素的影响。
发明内容
针对现有技术的不足,本发明提供一种有杆泵抽油井故障诊断方法。
本发明的技术方案是:
一种有杆泵抽油井故障诊断方法,包括以下步骤:
步骤1:获取M个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的有杆泵抽油井地面示功图;
步骤2:将M+1个有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图,并对有杆泵抽油井泵示功图进行归一化处理,得到归一化处理后的泵示功图;
步骤3:采用16方向链码对归一化处理后的泵示功图进行边界链码重画,得到边界链码重画后的泵示功图,其中,边界链码的类型由zt表示,t=1、2、…16,NUMt表示边界链码类型为zt的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数;
步骤4:将边界链码重画后的泵示功图进行区域划分,并提取各区域的特征向量,各特征向量以16方向链码形式表示;划分的区域包括:左上区、右上区、左下区和右下区;所述的特征向量包括:左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量;
步骤4.1:确定有杆泵抽油井泵示功图的上冲程和下冲程:将有杆泵抽油井泵示功图曲线的横坐标最大的点作为上冲程的结束点,则有杆泵抽油井泵示功图的上冲程为有杆泵抽油井泵示功图曲线第一个数据点到上冲程结束的数据点,有杆泵抽油井泵示功图的下冲程为有杆泵抽油井泵示功图曲线上冲程结束点的下一个数据点到最后一个数据点;
步骤4.2:提取有杆泵抽油井泵示功图各区域的泵示功图特征向量,包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量;
左上区欠缺程度特征向量表示为:
左上区快速加载程度特征向量表示为:
右上区突出程度特征向量表示为:
其中,为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值,为边界链码的类型为z3、z4和z5的边界链码重画后的泵示功图的数据点纵坐标的平均值;
右上区欠缺程度特征向量表示为:
右上区快速卸载程度特征向量表示为:
左下区突出程度特征向量表示为:
其中,为有杆泵抽油井泵示功图的下冲程数据点纵坐标的平均值,为边界链码的类型为z11、z12、z13、z14、z15和z16的边界链码重画后的泵示功图的数据点纵坐标的平均值;
左下区欠缺程度特征向量表示为:
右下区欠缺程度特征向量表示为:
右下区快速卸载程度特征向量表示为:
步骤4.3:将各区域的泵示功图特征向量进行量化,得到量化后的泵示功图特征向量;
步骤5:利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类,完成有杆泵抽油井故障诊断;
步骤5.1:设将M+1个泵示功图特征向量聚类为k类,σ为尺度参数,其中,σ∈[0,1];
步骤5.2:将聚类数k和尺度参数σ作为粒子进行聚类优化,得到最优的聚类数k和最优的尺度参数σ;
步骤5.2.1:将聚类数k和尺度参数σ作为粒子,随机产生若干个粒子,初始化粒子的速度和位置;
步骤5.2.2:M+1个量化后的泵示功图特征向量由d1、d2…dM+1表示,由高斯函数s(dp,dq)=exp(-||dp-dq||2/(2σ2))计算任意两个量化后的泵示功图特征向量的相似度,各相似度值构成相似度矩阵S∈R(M+1)×(M+1),其中,p、q=1、2…M+1;
步骤5.2.3:将每个量化后的泵示功图特征向量看作顶点V,根据量化后的泵示功图特征向量间的相似度对各顶点间的边E赋权重值W,得到一个基于量化后的泵示功图特征向量相似度的无向加权图G=(V,E);
步骤5.2.4:将相似度矩阵S的每行元素相加,即得到对应顶点的度,以所有顶点的度作为对角元素构成的对角矩阵,即度矩阵D;
步骤5.2.5:根据度矩阵、各顶点间的边的权重及单位矩阵,建立拉普拉斯矩阵Lrw=I-D-1W,其中,I为单位矩阵;
步骤5.2.6:计算拉普拉斯矩阵Lrw前k个最小的特征值对应的特征向量,k个特征向量组成特征向量矩阵U;
步骤5.2.7:将特征向量矩阵U的第g行向量(eg)g=1,2,…,M+1作为聚类数据集;
步骤5.2.8:用CritC指标函数的倒数作为各个粒子的适应度值;
步骤5.2.9:若当前粒子的适应度值小于其历史适应度值,则将当前粒子的适应度值作为其历史适应度值;
步骤5.2.10:设定粒子的适应度值终止条件,若当前粒子的适应度值小于其终止条件,则当前的聚类数k和尺度参数σ为最优的聚类数k和最优的尺度参数σ,当前的聚类数据集(eg)g=1,2,…,M+1为优化后的聚类数据集(eg)g=1,2,…,M+1,否则,更新粒子的速度和位置,返回步骤5.2.3:
步骤5.3:利用最优的聚类数k和最优的尺度参数σ进行聚类优化,得到最优的聚类中心;
步骤5.4:根据最优的聚类中心将M+1个泵示功图特征向量聚类为k类,得到最终的聚类结果:H1、H2…Hk,若待诊断故障类型的泵示功图特征向量与某知故障类型的泵示功图特征向量聚为一类,则两者为同类故障,即得到有杆泵抽油井故障诊断结果,若待诊断故障类型的泵示功图特征向量单独聚为一类,则该待诊断故障类型为新故障类型。
本发明的有益效果是:泵示功图曲线能够很好地描述曲线的轮廓特征;所提取的泵示功图曲线不同区域的典型泵示功图特征向量能够精确地描述泵示功图特征,很好地反映实际生产特点,从而提高泵示功图的识别准确率;采用16方向链码的数量对各个泵示功图特征向量进行量化能够有效、细致地描述图形特征;利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类,完成有杆泵抽油井故障诊断,不依赖训练样本,可以提高诊断的性能。
附图说明
图1为本发明具体实施方式中一种有杆泵抽油井故障诊断方法的流程图;
图2为本发明具体实施方式中的有杆泵抽油井泵示功图;
图3为本发明具体实施方式中的归一化处理后的泵示功图;
图4为本发明具体实施方式中采用的16个方向链码示意图;
图5为本发明具体实施方式中16个方向链码按逆时针顺序构成的区域边界示意图;
图6为本发明具体实施方式中采用16方向链码进行边界链码重画后的泵示功图;
图7为本发明具体实施方式中左上区欠缺程度特征链码示意图;
图8为本发明具体实施方式中左上区快速加载程度特征链码示意图;
图9为本发明具体实施方式中右上区欠缺程度特征链码示意图;
图10为本发明具体实施方式中右上区快速卸载程度特征链码示意图;
图11为本发明具体实施方式中左下区欠缺程度特征链码示意图;
图12为本发明具体实施方式中右下区欠缺程度特征链码示意图;
图13为本发明具体实施方式中右下区快速卸载程度特征链码示意图;
图14为本发明具体实施方式中对量化后的泵示功图特征向量进行聚类为k类后的示意图。
具体实施方式
下面结合附图对本发明具体实施方式加以详细的说明。
一种有杆泵抽油井故障诊断方法,如图1所示,包括以下步骤:
步骤1:获取M个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的有杆泵抽油井地面示功图。
本实施方式中,选用108个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的有杆泵抽油井地面示功图,其中,已知故障类型的有杆泵抽油井地面示功图有5种,将“正常”也看作一种故障类型,则故障类型分别为:24个“正常”类型的有杆泵抽油井地面示功图、20个“气体影响”类型的有杆泵抽油井地面示功图、24个“供液不足”类型的有杆泵抽油井地面示功图、20个“泵上碰”类型的有杆泵抽油井地面示功图、20个“泵下碰”类型的有杆泵抽油井地面示功图。
步骤2:将109个有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图,并对有杆泵抽油井泵示功图进行归一化处理,得到归一化处理后的泵示功图。
建立用来描述有杆泵抽油井抽油杆运动和应力传播的一维粘滞阻尼波动方程,采用傅里叶系数法求解该一维粘滞阻尼波动方程,得到抽油杆任意截面处的位移和载荷,进而得到有杆泵抽油井抽油泵处的位移和载荷,从而确定有杆泵抽油井泵示功图。
如图2所示,有杆泵抽油井泵示功图的每一个数据都由离散数据对(xi,yi)构成,其中:i=1,2,…,N,N为构成有杆泵抽油井泵示功图的离散数据对数;xi为有杆泵抽油井泵示功图位移数据,yi为有杆泵抽油井泵示功图载荷数据。
对有杆泵抽油井泵示功图进行归一化处理,得到归一化处理后的泵示功图。
对有杆泵抽油井泵示功图按如下归一化公式进行归一化处理:
其中,i=1,2,…,N;
xi为有杆泵抽油井泵示功图位移数据;
yi为有杆泵抽油井泵示功图载荷数据;
为归一化处理后的泵示功图的位移数据;
为归一化处理后的泵示功图的载荷数据。
归一化处理后的泵示功图如图3所示。
步骤3:采用16方向链码对归一化处理后的泵示功图进行边界链码重画,得到边界链码重画后的泵示功图,其中,边界链码的类型由zt表示,t=1、2、…16,NUMt表示边界链码类型为zt的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
如图4所示,16方向链码的角度分别为:0°、22.5°、45°、67.5°、90°、112.5°、135°、157.5°、180°、202.5°、225°、247.5°、270°、292.5°、315°、337.5°。
边界链码的类型由zt表示,t=1、2、…16。
采用的16方向链码分别为:0°(1)、22.5°(2)、45°(3)、67.5°(4)、90°(5)、112.5°(6)、135°(7)、157.5°(8)、180°(9)、202.5°(10)、225°(11)、247.5°(12)、270°(13)、292.5°(14)、315°(15)、337.5°(16)。
16方向链码在0°、90°、180°、270°的长度为1,在45°、135°、225°、315°的长度为在22.5°、67.5°、12.5°、157.5°、202.5°、247.5°、292.5°、337.5°的长度为
如图5所示,为根据16个方向链码重画泵示功图的曲线轮廓,其按照逆时针方向顺序构成的区域边界示意图。
由ωi表示第i个离散数据对到第i+1个离散数据对的角度,其计算公式如下:
那么,边界链码的类型zt可以由下式表示:
其中,表示取计算的最小值的类型,t=1、2、…16。
将归一化处理后的泵示功图上的两个相邻的数据点转换成边界链码,每一个方向的链码表示前一个数据点指向后一个数据点的角度。将归一化处理后的泵示功图上所有方向的边界链码连接成曲线,得到边界链码重画后的泵示功图。
采用16方向链码进行边界链码重画后的泵示功图如图6所示。
NUMt表示边界链码类型为zt的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
步骤4:将边界链码重画后的泵示功图进行区域划分,并提取各区域的特征向量,各特征向量以16方向链码形式表示;划分的区域包括:左上区、右上区、左下区和右下区;所述的特征向量包括:左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量。
步骤4.1:确定有杆泵抽油井泵示功图的上冲程和下冲程:将有杆泵抽油井泵示功图曲线的横坐标最大的点作为上冲程的结束点,则有杆泵抽油井泵示功图的上冲程为有杆泵抽油井泵示功图曲线第一个数据点到上冲程结束的数据点,有杆泵抽油井泵示功图的下冲程为有杆泵抽油井泵示功图曲线上冲程结束点的下一个数据点到最后一个数据点。
为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值;
为有杆泵抽油井泵示功图的下冲程数据点纵坐标的平均值。
步骤4.2:提取有杆泵抽油井泵示功图各区域的泵示功图特征向量,包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量。
由于系统加载时变得缓慢,使有杆泵抽油井泵示功图曲线的增载线后半部分曲率变化较小,造成有杆泵抽油井泵示功图的图形左上区欠缺,因此左上区欠缺程度特征向量表示为:
其中,NUM1表示边界链码类型为z1的边界链码个数,NUM2表示边界链码类型为z2的边界链码个数,NUM3表示边界链码类型为z3的边界链码个数,NUM4表示边界链码类型为z4的边界链码个数,NUM5表示边界链码类型为z5的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
左上区欠缺程度特征链码示意图如图7所示,其中,实线代表系统正常加载,虚线代表系统加载变缓。
由于系统突然急剧加载,使有杆泵抽油井泵示功图曲线的增载线前半部分曲率变化较大,因此,左上区快速加载程度特征向量表示为:
其中,NUM1表示边界链码类型为z1的边界链码个数,NUM2表示边界链码类型为z2的边界链码个数,NUM3表示边界链码类型为z3的边界链码个数,NUM4表示边界链码类型为z4的边界链码个数,NUM5表示边界链码类型为z5的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
左上区快速加载程度特征链码示意图如图8所示,其中,实线代表系统正常加载,虚线代表系统急剧加载。
在上冲程末端由于系统内部相关器件发生碰撞,使载荷突然增大,造成泵示功图曲线的上冲程尾部的曲率变化较大,因此右上区突出程度特征向量表示为:
其中,NUM3表示边界链码类型为z3的边界链码个数,NUM4表示边界链码类型为z4的边界链码个数,NUM5表示边界链码类型为z5的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数,为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值,为边界链码的类型为z3、z4和z5的边界链码重画后的泵示功图的数据点纵坐标的平均值。
由于系统卸载提前,使泵示功图曲线的上冲程线后半部分曲率变化较小,因此右上区欠缺程度特征向量表示为:
其中,NUM1表示边界链码类型为z1的边界链码个数,NUM2表示边界链码类型为z2的边界链码个数,NUM14表示边界链码类型为z14的边界链码个数,NUM15表示边界链码类型为z15的边界链码个数,NUM16表示边界链码类型为z16的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
右上区欠缺程度特征链码示意图如图9所示,其中,实线代表系统正常加载,虚线代表系统卸载变缓。
系统在上冲程结束时快速卸载,使泵示功图曲线的上死点处曲率变化较大,因此右上区快速卸载程度特征向量表示为:
其中,NUM11表示边界链码类型为z11的边界链码个数,NUM12表示边界链码类型为z12的边界链码个数,NUM13表示边界链码类型为z13的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
右上区快速卸载程度特征链码示意图如图10所示,其中,实线代表系统正常加载,虚线代表系统快速卸载。
在下冲程末端,由于系统内部相关器件发生碰撞,使载荷突然减小,造成泵示功图曲线的下冲程尾部的曲率变化较大,因此左下区突出程度特征向量表示为:
其中,NUM11表示边界链码类型为z11的边界链码个数,NUM12表示边界链码类型为z12的边界链码个数,NUM13表示边界链码类型为z13的边界链码个数,NUM14表示边界链码类型为z14的边界链码个数,NUM15表示边界链码类型为z15的边界链码个数,NUM16表示边界链码类型为z16的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数,为有杆泵抽油井泵示功图的下冲程数据点纵坐标的平均值,为边界链码的类型为z11、z12、z13、z14、z15和z16的边界链码重画后的泵示功图的数据点纵坐标的平均值。
由于系统的加载提前,使泵示功图曲线的下冲程线后半部分曲率变化较小,因此左下区欠缺程度特征向量表示为:
其中,NUM6表示边界链码类型为z6的边界链码个数,NUM7表示边界链码类型为z7的边界链码个数,NUM8表示边界链码类型为z8的边界链码个数,NUM9表示边界链码类型为z9的边界链码个数,NUM10表示边界链码类型为z9的边界链码个数,NUM10表示边界链码类型为z10的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
左下区欠缺程度特征链码示意图如图11所示,其中,实线代表系统正常加载,虚线代表系统加载提前。
由于系统在卸载时变缓,使泵示功图曲线的卸载线后半部分曲率变化较小,因此右下区欠缺程度特征向量表示为:
其中,NUM10表示边界链码类型为z9的边界链码个数,NUM10表示边界链码类型为z10的边界链码个数,NUM11表示边界链码类型为z11的边界链码个数,NUM12表示边界链码类型为z12的边界链码个数,NUM13表示边界链码类型为z13的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
右下区欠缺程度特征链码示意图如图12所示,其中,实线代表系统正常加载,虚线代表系统加载提前。
由于系统发生急剧卸载,使泵示功图曲线的卸载线的曲率变化较大,因此右下区快速卸载程度特征向量表示为:
其中,NUM9表示边界链码类型为z9的边界链码个数,NUM11表示边界链码类型为z11的边界链码个数,NUM12表示边界链码类型为z12的边界链码个数,NUM13表示边界链码类型为z13的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数。
右下区快速卸载程度特征链码示意图如图13所示,其中,实线代表系统正常加载,虚线代表系统加载提前。
步骤4.3:将各区域的泵示功图特征向量进行量化,得到量化后的泵示功图特征向量。
本发明具体实施方式中得到的一个量化后的泵示功图特征向量如表1所示:
表1 量化后的泵示功图特征向量
步骤5:利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类,完成有杆泵抽油井故障诊断。
步骤5.1:设将109个泵示功图特征向量聚类为k类,σ为尺度参数,其中,σ∈[0,1]。
步骤5.2:将聚类数k和尺度参数σ作为粒子进行聚类优化,得到最优的聚类数k和最优的尺度参数σ。
步骤5.2.1:将聚类数k和尺度参数σ作为粒子,随机产生若干个粒子,初始化粒子的速度和位置。
步骤5.2.2:109个量化后的泵示功图特征向量由d1、d2…d109表示,由高斯函数s(dp,dq)=exp(-||dp-dq||2/(2σ2))计算任意两个量化后的泵示功图特征向量的相似度,构成相似度矩阵S∈R109×109,其中,p、q=1、2…109。
步骤5.2.3:将每个量化后的泵示功图特征向量看作顶点V,根据量化后的泵示功图特征向量间的相似度对各顶点间的边E赋权重值W,得到一个基于量化后的泵示功图特征向量相似度的无向加权图G=(V,E)。
步骤5.2.4:将相似度矩阵S的每行元素相加,即得到对应顶点的度,以所有顶点的度作为对角元素构成的对角矩阵,即度矩阵D。
步骤5.2.5:根据度矩阵、各顶点间的边的权重及单位矩阵,建立拉普拉斯矩阵Lrw=I-D-1W,其中,I为单位矩阵。
步骤5.2.6:计算拉普拉斯矩阵Lrw前k个最小的特征值对应的特征向量u1、u2…uk,k个特征向量组成特征向量矩阵U∈R109×k。
步骤5.2.7:将特征向量矩阵U∈R109×k的第g行向量(eg)g=1、2...109作为聚类数据集。
步骤5.2.8:用CritC指标函数的倒数作为各个粒子的适应度值。
CritC指标函数的定义如下:
其中:
γ(Cf,O)表示聚类数据集(eg)g=1、2...109中每个聚类中心Cf和整个聚类数据集(eg)g=1、2...109的中心0之间的距离;
Cf由每个聚类中所有数据的平均值求得;
0由聚类数据集(eg)g=1、2...109中所有数据的平均值求得;
|Cf|表示γ(Cf,O)所对应的权值,等于每个聚类中数据的数量;
γ(Cf,o)表示聚类数据集(eg)g=1、2...109中每个数据点和它所在聚类的中心Cf之间的距离;
a为常数系数,本实施方式中取a=0.8;
CritC指标函数在[0,1]区间内取值。
步骤5.2.9:若当前粒子的适应度值小于其历史适应度值,则将当前粒子的适应度值作为其历史适应度值。
步骤5.2.10:设定粒子的适应度值终止条件,若当前粒子的适应度值小于其终止条件,则当前的聚类数k和尺度参数σ为最优的聚类数k和最优的尺度参数σ,当前的聚类数据集(eg)g=1、2...109为优化后的聚类数据集(eg)g=1、2...109,否则,更新粒子的速度和位置,返回步骤5.2.3。
本实施方式中,设定代表聚类数k和尺度参数σ的粒子的适应度值终止条件为0.5%,聚类数k和尺度参数σ的最优值为5和0.1237。
步骤5.3:利用最优的聚类数k和最优的尺度参数σ进行聚类优化,得到最优的聚类中心。
步骤5.3.1:将优化后的聚类数据集(eg)g=1、2...109中数据作为粒子,随机选取k=5个粒子作为初始聚类中心,初始化粒子的速度和位置。
步骤5.3.2:用Jc指标函数作为粒子的适应度值。
Jc指标函数的定义如下:
其中:φNc表示优化后的聚类数据集(eg)g=1、2...109的聚类中心,Nc=1、2、...5。
步骤5.3.3:若当前粒子的适应度值小于其历史适应度值时,则把当前粒子的适应度值作为其历史适应度值,当前粒子作为最优粒子,同时将当前粒子位置记为其历史最优位置。
步骤5.3.4:设定粒子的适应度值终止条件,若当前粒子的适应度值小于其终止条件,则当前最优粒子的聚类中心为最优的聚类中心,否则,更新粒子的速度和位置,返回步骤5.3.2。
步骤5.4:根据最优的聚类中心将109个泵示功图特征向量聚类为5类,得到最终的聚类结果:H1、H2…H5,若待诊断故障类型的泵示功图特征向量与某知故障类型的泵示功图特征向量聚为一类,则两者为同类故障,即得到有杆泵抽油井故障诊断结果,若待诊断故障类型的泵示功图特征向量单独聚为一类,则该待诊断故障类型为新故障类型。
如图14所示,为本发明具体实施方式中对量化后的泵示功图特征向量进行聚类为5类后的示意图。
本实施方式中,待诊断故障类型的量化后的泵示功图特征向量(由“☆”标记)和“正常”类型的量化后的泵示功图特征向量聚类为同一类型,因此该待诊断故障类型为“正常”。
Claims (2)
1.一种有杆泵抽油井故障诊断方法,其特征在于,包括以下步骤:
步骤1:获取M个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的有杆泵抽油井地面示功图;
步骤2:将M+1个有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图,并对有杆泵抽油井泵示功图进行归一化处理,得到归一化处理后的泵示功图;
步骤3:采用16方向链码对归一化处理后的泵示功图进行边界链码重画,得到边界链码重画后的泵示功图,其中,边界链码的类型由zt表示,t=1、2、…16,NUMt表示边界链码类型为zt的边界链码个数,NUM表示边界链码重画后的泵示功图的边界链码总数;
步骤4:将边界链码重画后的泵示功图进行区域划分,并提取各区域的特征向量,各特征向量以16方向链码形式表示;划分的区域包括:左上区、右上区、左下区和右下区;所述的特征向量包括:左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量;
步骤4.1:确定有杆泵抽油井泵示功图的上冲程和下冲程:将有杆泵抽油井泵示功图曲线的横坐标最大的点作为上冲程的结束点,则有杆泵抽油井泵示功图的上冲程为有杆泵抽油井泵示功图曲线第一个数据点到上冲程结束的数据点,有杆泵抽油井泵示功图的下冲程为有杆泵抽油井泵示功图曲线上冲程结束点的下一个数据点到最后一个数据点;
步骤4.2:提取有杆泵抽油井泵示功图各区域的泵示功图特征向量,包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量;
左上区欠缺程度特征向量表示为:
左上区快速加载程度特征向量表示为:
右上区突出程度特征向量表示为:
其中,为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值,为边界链码的类型为z3、z4和z5的边界链码重画后的泵示功图的数据点纵坐标的平均值;
右上区欠缺程度特征向量表示为:
右上区快速卸载程度特征向量表示为:
左下区突出程度特征向量表示为:
其中,为有杆泵抽油井泵示功图的下冲程数据点纵坐标的平均值,为边界链码的类型为z11、z12、z13、z14、z15和z16的边界链码重画后的泵示功图的数据点纵坐标的平均值;
左下区欠缺程度特征向量表示为:
右下区欠缺程度特征向量表示为:
右下区快速卸载程度特征向量表示为:
步骤4.3:将各区域的泵示功图特征向量进行量化,得到量化后的泵示功图特征向量;
步骤5:利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类,完成有杆泵抽油井故障诊断;
步骤5.1:设将M+1个泵示功图特征向量聚类为k类,σ为尺度参数,其中,σ∈[0,1];
步骤5.2:将聚类数k和尺度参数σ作为粒子进行聚类优化,得到最优的聚类数k和最优的尺度参数σ;
步骤5.3:利用最优的聚类数k和最优的尺度参数σ进行聚类优化,得到最优的聚类中心;
步骤5.4:根据最优的聚类中心将M+1个泵示功图特征向量聚类为k类,得到最终的聚类结果:H1、H2…Hk,若待诊断故障类型的泵示功图特征向量与某知故障类型的泵示功图特征向量聚为一类,则两者为同类故障,即得到有杆泵抽油井故障诊断结果,若待诊断故障类型的泵示功图特征向量单独聚为一类,则该待诊断故障类型为新故障类型,对该故障类型的泵示功图特征向量进行保存。
2.根据权利要求1所述的一种有杆泵抽油井故障诊断方法,其特征在于,所述的步骤5.2按以下步骤执行:
步骤5.2.1:将聚类数k和尺度参数σ作为粒子,随机产生若干个粒子,初始化粒子的速度和位置;
步骤5.2.2:M+1个量化后的泵示功图特征向量由d1、d2…dM+1表示,由高斯函数s(dp,dq)=exp(-||dp-dq||2/(2σ2))计算任意两个量化后的泵示功图特征向量的相似度,各相似度值构成相似度矩阵S∈R(M+1)×(M+1),其中,p、q=1、2…M+1;
步骤5.2.3:将每个量化后的泵示功图特征向量看作顶点V,根据量化后的泵示功图特征向量间的相似度对各顶点间的边E赋权重值W,得到一个基于量化后的泵示功图特征向量相似度的无向加权图G=(V,E);
步骤5.2.4:将相似度矩阵S的每行元素相加,即得到对应顶点的度,以所有顶点的度作为对角元素构成的对角矩阵,即度矩阵D;
步骤5.2.5:根据度矩阵、各顶点间的边的权重及单位矩阵,建立拉普拉斯矩阵Lrw=I-D-1W,其中,I为单位矩阵;
步骤5.2.6:计算拉普拉斯矩阵Lrw前k个最小的特征值对应的特征向量,k个特征向量组成特征向量矩阵U;
步骤5.2.7:将特征向量矩阵U的第g行向量(eg)g=1,2,…,M+1作为聚类数据集;
步骤5.2.8:用CritC指标函数的倒数作为各个粒子的适应度值;
步骤5.2.9:若当前粒子的适应度值小于其历史适应度值,则将当前粒子的适应度值作为其历史适应度值;
步骤5.2.10:设定粒子的适应度值终止条件,若当前粒子的适应度值小于其终止条件,则当前的聚类数k和尺度参数σ为最优的聚类数k和最优的尺度参数σ,当前的聚类数据集(eg)g=1,2,…,M+1为优化后的聚类数据集(eg)g=1,2,…,M+1,否则,更新粒子的速度和位置,返回步骤5.2.3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410473163.1A CN104481496B (zh) | 2014-09-16 | 2014-09-16 | 一种有杆泵抽油井故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410473163.1A CN104481496B (zh) | 2014-09-16 | 2014-09-16 | 一种有杆泵抽油井故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104481496A CN104481496A (zh) | 2015-04-01 |
CN104481496B true CN104481496B (zh) | 2017-02-15 |
Family
ID=52756088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410473163.1A Expired - Fee Related CN104481496B (zh) | 2014-09-16 | 2014-09-16 | 一种有杆泵抽油井故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104481496B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107165615A (zh) * | 2017-05-10 | 2017-09-15 | 东北大学 | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105257277B (zh) * | 2015-05-15 | 2018-04-06 | 渤海大学 | 基于多变量灰色模型的有杆泵抽油井井下故障预测方法 |
CN106121622B (zh) * | 2016-07-27 | 2019-10-25 | 渤海大学 | 一种基于示功图的有杆泵抽油井的多故障诊断方法 |
CN106951662B (zh) * | 2017-04-12 | 2020-03-20 | 东北大学 | 基于凡尔工作点的有杆泵抽油井井下工况诊断方法 |
CN106930751B (zh) * | 2017-04-20 | 2020-02-11 | 东北大学 | 一种有杆泵抽油井故障分离方法 |
CN106884644B (zh) * | 2017-04-26 | 2020-12-15 | 中国石油大学(华东) | 基于时序地面示功图的抽油机井实时工况诊断方法 |
CN108729902B (zh) * | 2018-05-03 | 2021-09-10 | 西安永瑞自动化有限公司 | 抽油机在线故障诊断系统及其诊断方法 |
CN109460561B (zh) * | 2018-07-25 | 2021-08-17 | 北京国双科技有限公司 | 采油设备故障诊断方法及装置 |
CN108979624B (zh) * | 2018-08-07 | 2022-03-08 | 东北大学 | 一种基于示功图矩特征的有杆抽油系统摩擦因数辨识方法 |
CN112949196B (zh) * | 2021-03-11 | 2023-05-12 | 中国石油大学(北京) | 一种基于残差神经网络的抽油机井故障诊断方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1786410A (zh) * | 2005-12-20 | 2006-06-14 | 西安理工大学 | 油井有杆抽油系统阻尼系数的辨识方法 |
CN101105119A (zh) * | 2006-07-11 | 2008-01-16 | 刘克文 | 一种抽油机平衡调节仪及抽油机平衡结构调节方法 |
CN101446283A (zh) * | 2008-12-17 | 2009-06-03 | 西安交通大学 | 一种基于进化计算诊断压缩机故障的方法 |
CN102184414A (zh) * | 2011-05-16 | 2011-09-14 | 中国石油天然气股份有限公司 | 泵示功图的识别和判断方法及其系统 |
CN103541723A (zh) * | 2013-11-12 | 2014-01-29 | 丁涛 | 基于地面示功图面积变化的抽油机井实时工况诊断方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9280517B2 (en) * | 2011-06-23 | 2016-03-08 | University Of Southern California | System and method for failure detection for artificial lift systems |
-
2014
- 2014-09-16 CN CN201410473163.1A patent/CN104481496B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1786410A (zh) * | 2005-12-20 | 2006-06-14 | 西安理工大学 | 油井有杆抽油系统阻尼系数的辨识方法 |
CN101105119A (zh) * | 2006-07-11 | 2008-01-16 | 刘克文 | 一种抽油机平衡调节仪及抽油机平衡结构调节方法 |
CN101446283A (zh) * | 2008-12-17 | 2009-06-03 | 西安交通大学 | 一种基于进化计算诊断压缩机故障的方法 |
CN102184414A (zh) * | 2011-05-16 | 2011-09-14 | 中国石油天然气股份有限公司 | 泵示功图的识别和判断方法及其系统 |
CN103541723A (zh) * | 2013-11-12 | 2014-01-29 | 丁涛 | 基于地面示功图面积变化的抽油机井实时工况诊断方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107165615A (zh) * | 2017-05-10 | 2017-09-15 | 东北大学 | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104481496A (zh) | 2015-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104481496B (zh) | 一种有杆泵抽油井故障诊断方法 | |
Li et al. | Using the curve moment and the PSO-SVM method to diagnose downhole conditions of a sucker rod pumping unit | |
CN106930751B (zh) | 一种有杆泵抽油井故障分离方法 | |
CN104110251A (zh) | 一种基于art2的抽油机示功图识别方法 | |
CN108447057B (zh) | 基于显著性和深度卷积网络的sar图像变化检测方法 | |
CN105257277B (zh) | 基于多变量灰色模型的有杆泵抽油井井下故障预测方法 | |
Zheng et al. | Sucker rod pumping diagnosis using valve working position and parameter optimal continuous hidden Markov model | |
CN105426889A (zh) | 基于pca混合特征融合的气液两相流流型识别方法 | |
CN104008551B (zh) | 一种基于可见光图像的柑橘黄龙病检测方法 | |
CN102930533A (zh) | 一种基于改进k-均值聚类的半监督高光谱影像降维方法 | |
CN106121622B (zh) | 一种基于示功图的有杆泵抽油井的多故障诊断方法 | |
CN106951662B (zh) | 基于凡尔工作点的有杆泵抽油井井下工况诊断方法 | |
CN110738271B (zh) | 一种锌浮选过程精矿品位预测方法 | |
CN104234695A (zh) | 一种基于神经网络的油井故障诊断方法 | |
CN110288257A (zh) | 一种深度超限示功图学习方法 | |
CN111639882A (zh) | 一种基于深度学习的用电风险的判定方法 | |
CN111946329A (zh) | 一种油井动液面求取方法 | |
CN107016260B (zh) | 一种基于跨平台基因表达数据的基因调控网络重建方法 | |
CN111379555A (zh) | 抽油机井的故障诊断方法及装置 | |
CN108985161A (zh) | 一种基于拉普拉斯正则化的低秩稀疏表征图像特征学习方法 | |
CN104166806B (zh) | 一种井间示踪曲线聚类分析方法及装置 | |
CN107165615B (zh) | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 | |
CN107507192A (zh) | 一种水稻剑叶夹角在体测量便携装置及在体测量方法 | |
CN109460781B (zh) | 一种基于决策树归纳学习的牛排等级划分方法 | |
Zhang et al. | Diagnosis of pumping machine working conditions based on transfer learning and ViT model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 Termination date: 20170916 |