《一种有杆泵抽油井故障诊断方法.pdf》由会员分享，可在线阅读，更多相关《一种有杆泵抽油井故障诊断方法.pdf（24页珍藏版）》请在专利查询网上搜索。
1、10申请公布号43申请公布日21申请号201410473163122申请日20140916E21B44/0020060171申请人渤海大学地址121013辽宁省锦州市松山新区科技路19号72发明人李琨韩莹74专利代理机构沈阳东大知识产权代理有限公司21109代理人梁焱54发明名称一种有杆泵抽油井故障诊断方法57摘要本发明涉及一种有杆泵抽油井故障诊断方法，获取已知和待诊断有杆泵抽油井地面示功图，将有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图，采用16方向链码对归一化处理后的泵示功图进行边界链码重画，将边界链码重画后的泵示功图进行区域划分，并提取各区域的特征向量，利用基于粒子群优化的谱聚类算法。
2、对量化后的泵示功图特征向量进行聚类，完成有杆泵抽油井故障诊断，本发明采用16方向链码的数量对各个泵示功图特征向量进行量化能够有效、细致地描述图形特征，有杆泵抽油井故障诊断是利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类，不依赖训练样本，可以提高诊断的性能。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书3页说明书11页附图9页10申请公布号CN104481496A43申请公布日20150401CN104481496A1/3页21一种有杆泵抽油井故障诊断方法，其特征在于，包括以下步骤步骤1获取M个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的。
3、有杆泵抽油井地面示功图；步骤2将M1个有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图，并对有杆泵抽油井泵示功图进行归一化处理，得到归一化处理后的泵示功图；步骤3采用16方向链码对归一化处理后的泵示功图进行边界链码重画，得到边界链码重画后的泵示功图，其中，边界链码的类型由ZT表示，T1、2、16，NUMT表示边界链码类型为ZT的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数；步骤4将边界链码重画后的泵示功图进行区域划分，并提取各区域的特征向量，各特征向量以16方向链码形式表示；划分的区域包括左上区、右上区、左下区和右下区；所述的特征向量包括左上区欠缺程度特征向量、左上区快速加载。
4、程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量；步骤41确定有杆泵抽油井泵示功图的上冲程和下冲程将有杆泵抽油井泵示功图曲线的横坐标最大的点作为上冲程的结束点，则有杆泵抽油井泵示功图的上冲程为有杆泵抽油井泵示功图曲线第一个数据点到上冲程结束的数据点，有杆泵抽油井泵示功图的下冲程为有杆泵抽油井泵示功图曲线上冲程结束点的下一个数据点到最后一个数据点；步骤42提取有杆泵抽油井泵示功图各区域的泵示功图特征向量，包括左上区欠缺程度特征向量、左上区快速加载程度特征向量。
5、、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量；左上区欠缺程度特征向量表示为左上区快速加载程度特征向量表示为右上区突出程度特征向量表示为其中，为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值，为边界链码的类型为Z3、Z4和Z5的边界链码重画后的泵示功图的数据点纵坐标的平均值；右上区欠缺程度特征向量表示为权利要求书CN104481496A2/3页3右上区快速卸载程度特征向量表示为左下区突出程度特征向量表示为其中，为有杆泵抽油井泵示功图的下冲程数据点纵坐标的平均值，为。
6、边界链码的类型为Z11、Z12、Z13、Z14、Z15和Z16的边界链码重画后的泵示功图的数据点纵坐标的平均值；左下区欠缺程度特征向量表示为右下区欠缺程度特征向量表示为右下区快速卸载程度特征向量表示为步骤43将各区域的泵示功图特征向量进行量化，得到量化后的泵示功图特征向量；步骤5利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类，完成有杆泵抽油井故障诊断；步骤51设将M1个泵示功图特征向量聚类为K类，为尺度参数，其中，0，1；步骤52将聚类数K和尺度参数作为粒子进行聚类优化，得到最优的聚类数K和最优的尺度参数；步骤53利用最优的聚类数K和最优的尺度参数进行聚类优化，得到最优的聚类。
7、中心；步骤54根据最优的聚类中心将M1个泵示功图特征向量聚类为K类，得到最终的聚类结果H1、H2HK，若待诊断故障类型的泵示功图特征向量与某知故障类型的泵示功图特征向量聚为一类，则两者为同类故障，即得到有杆泵抽油井故障诊断结果，若待诊断故障类型的泵示功图特征向量单独聚为一类，则该待诊断故障类型为新故障类型，对该故障类型的泵示功图特征向量进行保存。2根据权利要求1所述的一种有杆泵抽油井故障诊断方法，其特征在于，所述的步骤52按以下步骤执行步骤521将聚类数K和尺度参数作为粒子，随机产生若干个粒子，初始化粒子的速度和位置；步骤522M1个量化后的泵示功图特征向量由D1、D2DM1表示，由高斯函数S。
8、DP，权利要求书CN104481496A3/3页4DQEXP|DPDQ|2/22计算任意两个量化后的泵示功图特征向量的相似度，各相似度值构成相似度矩阵SRM1M1，其中，P、Q1、2M1；步骤523将每个量化后的泵示功图特征向量看作顶点V，根据量化后的泵示功图特征向量间的相似度对各顶点间的边E赋权重值W，得到一个基于量化后的泵示功图特征向量相似度的无向加权图GV，E；步骤524将相似度矩阵S的每行元素相加，即得到对应顶点的度，以所有顶点的度作为对角元素构成的对角矩阵，即度矩阵D；步骤525根据度矩阵、各顶点间的边的权重及单位矩阵，建立拉普拉斯矩阵LRWID1W，其中，I为单位矩阵；步骤526计。
9、算拉普拉斯矩阵LRW前K个最小的特征值对应的特征向量，K个特征向量组成特征向量矩阵U；步骤527将特征向量矩阵U的第G行向量EGG1，2，M1作为聚类数据集；步骤528用CRITC指标函数的倒数作为各个粒子的适应度值；步骤529若当前粒子的适应度值小于其历史适应度值，则将当前粒子的适应度值作为其历史适应度值；步骤5210设定粒子的适应度值终止条件，若当前粒子的适应度值小于其终止条件，则当前的聚类数K和尺度参数为最优的聚类数K和最优的尺度参数，当前的聚类数据集EGG1，2，M1为优化后的聚类数据集EGG1，2，M1，否则，更新粒子的速度和位置，返回步骤523。权利要求书CN104481496A1。
10、/11页5一种有杆泵抽油井故障诊断方法技术领域0001本发明属于抽油井故障诊断领域，特别涉及一种有杆泵抽油井故障诊断方法。背景技术0002有杆泵抽油井泵示功图在油田生产中具有十分重要的作用，它是判断抽油系统井下工作状态的主要依据。泵示功图可以用来进行井下生产状况的故障诊断，也可以用来计算井口的出液量，还可以用来计算井下动液面的位置。传统生产中，一般依靠技术人员对采集的泵示功图进行人工分析，从而判断抽油系统的实际工作状况。这种人工方式工作周期长，受主观因素影响较大，越来越无法满足企业不断提高生产效率的需求。0003随着计算机、通信、电子技术等领域的不断进步，利用计算机取代人力劳动方式已成为国内外。
11、油田生产的发展趋势。鉴于泵示功图在生产中的重要性，采用计算机进行计算和分析已成为该领域研究的重点和难点。利用计算机处理泵示功图，关键在于提取其显著的特征，识别越准确，计算越正确。目前，对于泵示功图的特征提取方法，主要有面积法、匹配法、网格法、矢量法等。但是这些方法都不能很好地结合实际生产过程的特点，从而使得到的特征向量不能准确、细致地反映不同过程的显著特征。0004有杆泵抽油井的抽油泵一般工作在数千米的地下，生产环境恶劣，故障发生率很高，采用计算机实现诊断是很有意义的。现有的诊断方法主要以有监督学习为主，根据生产经验得到不同故障类型的训练样本集，然后由这些训练样本训练出不同故障类型的训练集模型。
12、，从而实现待诊断样本的诊断。这种有监督学习方法主要存在以下两方面问题一是诊断的准确率很大程度上取决于所构成训练样本集的样本特征向量是否能够准确反映样本的典型特征；二是这种有监督方式主要依靠人为划分训练样本集，容易受到主观因素的影响。发明内容0005针对现有技术的不足，本发明提供一种有杆泵抽油井故障诊断方法。0006本发明的技术方案是0007一种有杆泵抽油井故障诊断方法，包括以下步骤0008步骤1获取M个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的有杆泵抽油井地面示功图；0009步骤2将M1个有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图，并对有杆泵抽油井泵示功图进行归一化处理，。
13、得到归一化处理后的泵示功图；0010步骤3采用16方向链码对归一化处理后的泵示功图进行边界链码重画，得到边界链码重画后的泵示功图，其中，边界链码的类型由ZT表示，T1、2、16，NUMT表示边界链码类型为ZT的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数；0011步骤4将边界链码重画后的泵示功图进行区域划分，并提取各区域的特征向量，各特征向量以16方向链码形式表示；划分的区域包括左上区、右上区、左下区和右下区；说明书CN104481496A2/11页6所述的特征向量包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快。
14、速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量；0012步骤41确定有杆泵抽油井泵示功图的上冲程和下冲程将有杆抽油井泵示功图曲线的横坐标最大的点作为上冲程的结束点，则有杆泵抽油井泵示功图的上冲程为有杆泵抽油井泵示功图曲线第一个数据点到上冲程结束的数据点，有杆泵抽油井泵示功图的下冲程为有杆泵抽油井泵示功图曲线上冲程结束点的下一个数据点到最后一个数据点；0013步骤42提取有杆泵抽油井泵示功图各区域的泵示功图特征向量，包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区。
15、快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量；0014左上区欠缺程度特征向量表示为00150016左上区快速加载程度特征向量表示为00170018右上区突出程度特征向量表示为00190020其中，为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值，为边界链码的类型为Z3、Z4和Z5的边界链码重画后的泵示功图的数据点纵坐标的平均值；0021右上区欠缺程度特征向量表示为00220023右上区快速卸载程度特征向量表示为00240025左下区突出程度特征向量表示为00260027其中，为有杆泵抽油井泵示功图的下冲程数据点纵坐标。
16、的平均值，说明书CN104481496A3/11页7为边界链码的类型为Z11、Z12、Z13、Z14、Z15和Z16的边界链码重画后的泵示功图的数据点纵坐标的平均值；0028左下区欠缺程度特征向量表示为00290030右下区欠缺程度特征向量表示为00310032右下区快速卸载程度特征向量表示为00330034步骤43将各区域的泵示功图特征向量进行量化，得到量化后的泵示功图特征向量；0035步骤5利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类，完成有杆泵抽油井故障诊断；0036步骤51设将M1个泵示功图特征向量聚类为K类，为尺度参数，其中，0，1；0037步骤52将聚类数K和尺。
17、度参数作为粒子进行聚类优化，得到最优的聚类数K和最优的尺度参数；0038步骤521将聚类数K和尺度参数作为粒子，随机产生若干个粒子，初始化粒子的速度和位置；0039步骤522M1个量化后的泵示功图特征向量由D1、D2DM1表示，由高斯函数SDP，DQEXP|DPDQ|2/22计算任意两个量化后的泵示功图特征向量的相似度，各相似度值构成相似度矩阵SRM1M1，其中，P、Q1、2M1；0040步骤523将每个量化后的泵示功图特征向量看作顶点V，根据量化后的泵示功图特征向量间的相似度对各顶点间的边E赋权重值W，得到一个基于量化后的泵示功图特征向量相似度的无向加权图GV，E；0041步骤524将相似度。
18、矩阵S的每行元素相加，即得到对应顶点的度，以所有顶点的度作为对角元素构成的对角矩阵，即度矩阵D；0042步骤525根据度矩阵、各顶点间的边的权重及单位矩阵，建立拉普拉斯矩阵LRWID1W，其中，I为单位矩阵；0043步骤526计算拉普拉斯矩阵LRW前K个最小的特征值对应的特征向量，K个特征向量组成特征向量矩阵U；0044步骤527将特征向量矩阵U的第G行向量EGG1，2，M1作为聚类数据集；0045步骤528用CRITC指标函数的倒数作为各个粒子的适应度值；0046步骤529若当前粒子的适应度值小于其历史适应度值，则将当前粒子的适应度值作为其历史适应度值；说明书CN104481496A4/11。
19、页80047步骤5210设定粒子的适应度值终止条件，若当前粒子的适应度值小于其终止条件，则当前的聚类数K和尺度参数为最优的聚类数K和最优的尺度参数，当前的聚类数据集EGG1，2，M1为优化后的聚类数据集EGG1，2，M1，否则，更新粒子的速度和位置，返回步骤5230048步骤53利用最优的聚类数K和最优的尺度参数进行聚类优化，得到最优的聚类中心；0049步骤54根据最优的聚类中心将M1个泵示功图特征向量聚类为K类，得到最终的聚类结果H1、H2HK，若待诊断故障类型的泵示功图特征向量与某知故障类型的泵示功图特征向量聚为一类，则两者为同类故障，即得到有杆泵抽油井故障诊断结果，若待诊断故障类型的泵示。
20、功图特征向量单独聚为一类，则该待诊断故障类型为新故障类型。0050本发明的有益效果是泵示功图曲线能够很好地描述曲线的轮廓特征；所提取的泵示功图曲线不同区域的典型泵示功图特征向量能够精确地描述泵示功图特征，很好地反映实际生产特点，从而提高泵示功图的识别准确率；采用16方向链码的数量对各个泵示功图特征向量进行量化能够有效、细致地描述图形特征；利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类，完成有杆泵抽油井故障诊断，不依赖训练样本，可以提高诊断的性能。附图说明0051图1为本发明具体实施方式中一种有杆泵抽油井故障诊断方法的流程图；0052图2为本发明具体实施方式中的有杆泵抽油井泵示。
21、功图；0053图3为本发明具体实施方式中的归一化处理后的泵示功图；0054图4为本发明具体实施方式中采用的16个方向链码示意图；0055图5为本发明具体实施方式中16个方向链码按逆时针顺序构成的区域边界示意图；0056图6为本发明具体实施方式中采用16方向链码进行边界链码重画后的泵示功图；0057图7为本发明具体实施方式中左上区欠缺程度特征链码示意图；0058图8为本发明具体实施方式中左上区快速加载程度特征链码示意图；0059图9为本发明具体实施方式中右上区欠缺程度特征链码示意图；0060图10为本发明具体实施方式中右上区快速卸载程度特征链码示意图；0061图11为本发明具体实施方式中左下区欠。
22、缺程度特征链码示意图；0062图12为本发明具体实施方式中右下区欠缺程度特征链码示意图；0063图13为本发明具体实施方式中右下区快速卸载程度特征链码示意图；0064图14为本发明具体实施方式中对量化后的泵示功图特征向量进行聚类为K类后的示意图。具体实施方式0065下面结合附图对本发明具体实施方式加以详细的说明。0066一种有杆泵抽油井故障诊断方法，如图1所示，包括以下步骤说明书CN104481496A5/11页90067步骤1获取M个已知故障类型的有杆泵抽油井地面示功图和1个待诊断故障类型的有杆泵抽油井地面示功图。0068本实施方式中，选用108个已知故障类型的有杆泵抽油井地面示功图和1个待。
23、诊断故障类型的有杆泵抽油井地面示功图，其中，已知故障类型的有杆泵抽油井地面示功图有5种，将“正常”也看作一种故障类型，则故障类型分别为24个“正常”类型的有杆泵抽油井地面示功图、20个“气体影响”类型的有杆泵抽油井地面示功图、24个“供液不足”类型的有杆泵抽油井地面示功图、20个“泵上碰”类型的有杆泵抽油井地面示功图、20个“泵下碰”类型的有杆泵抽油井地面示功图。0069步骤2将109个有杆泵抽油井地面示功图转化为有杆泵抽油井泵示功图，并对有杆泵抽油井泵示功图进行归一化处理，得到归一化处理后的泵示功图。0070建立用来描述有杆泵抽油井抽油杆运动和应力传播的一维粘滞阻尼波动方程，采用傅里叶系数法。
24、求解该一维粘滞阻尼波动方程，得到抽油杆任意截面处的位移和载荷，进而得到有杆泵抽油井抽油泵处的位移和载荷，从而确定有杆泵抽油井泵示功图。0071如图2所示，有杆泵抽油井泵示功图的每一个数据都由离散数据对XI，YI构成，其中I1，2，N，N为构成有杆泵抽油井泵示功图的离散数据对数；XI为有杆泵抽油井泵示功图位移数据，YI为有杆泵抽油井泵示功图载荷数据。0072对有杆泵抽油井泵示功图进行归一化处理，得到归一化处理后的泵示功图。0073对有杆泵抽油井泵示功图按如下归一化公式进行归一化处理007400750076其中，I1，2，N；0077XI为有杆泵抽油井泵示功图位移数据；0078YI为有杆泵抽油井泵。
25、示功图载荷数据；0079为归一化处理后的泵示功图的位移数据；0080为归一化处理后的泵示功图的载荷数据。0081归一化处理后的泵示功图如图3所示。0082步骤3采用16方向链码对归一化处理后的泵示功图进行边界链码重画，得到边界链码重画后的泵示功图，其中，边界链码的类型由ZT表示，T1、2、16，NUMT表示边界链码类型为ZT的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0083如图4所示，16方向链码的角度分别为0、225、45、675、90、1125、135、1575、180、2025、225、2475、270、2925、315、3375。0084边界链码的类型由ZT表示。
26、，T1、2、16。0085采用的16方向链码分别为01、2252、453、6754、905、11256、1357、15758、1809、202510、22511、247512、说明书CN104481496A6/11页1027013、292514、31515、337516。008616方向链码在0、90、180、270的长度为1，在45、135、225、315的长度为在225、675、125、1575、2025、2475、2925、3375的长度为0087如图5所示，为根据16个方向链码重画泵示功图的曲线轮廓，其按照逆时针方向顺序构成的区域边界示意图。0088由I表示第I个离散数据对到第I1个离。
27、散数据对的角度，其计算公式如下00890090那么，边界链码的类型ZT可以由下式表示00910092其中，表示取计算的最小值的类型，T1、2、16。0093将归一化处理后的泵示功图上的两个相邻的数据点转换成边界链码，每一个方向的链码表示前一个数据点指向后一个数据点的角度。将归一化处理后的泵示功图上所有方向的边界链码连接成曲线，得到边界链码重画后的泵示功图。0094采用16方向链码进行边界链码重画后的泵示功图如图6所示。0095NUMT表示边界链码类型为ZT的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0096步骤4将边界链码重画后的泵示功图进行区域划分，并提取各区域的特征向。
28、量，各特征向量以16方向链码形式表示；划分的区域包括左上区、右上区、左下区和右下区；所述的特征向量包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量。0097步骤41确定有杆泵抽油井泵示功图的上冲程和下冲程将有杆抽油井泵示功图曲线的横坐标最大的点作为上冲程的结束点，则有杆泵抽油井泵示功图的上冲程为有杆泵抽油井泵示功图曲线第一个数据点到上冲程结束的数据点，有杆泵抽油井泵示功图的下冲程为有杆泵抽油井泵示功图曲线上冲程结束。
29、点的下一个数据点到最后一个数据点。0098为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值；0099为有杆泵抽油井泵示功图的下冲程数据点纵坐标的平均值。0100步骤42提取有杆泵抽油井泵示功图各区域的泵示功图特征向量，包括左上区欠缺程度特征向量、左上区快速加载程度特征向量、右上区突出程度特征向量、右上区欠缺程度特征向量、右上区快速卸载程度特征向量、左下区突出程度特征向量、左下区欠缺程度说明书CN104481496A7/11页11特征向量、右下区欠缺程度特征向量和右下区快速卸载程度特征向量。0101由于系统加载时变得缓慢，使有杆泵抽油井泵示功图曲线的增载线后半部分曲率变化较小，造成有杆泵抽油井。
30、泵示功图的图形左上区欠缺，因此左上区欠缺程度特征向量表示为01020103其中，NUM1表示边界链码类型为Z1的边界链码个数，NUM2表示边界链码类型为Z2的边界链码个数，NUM3表示边界链码类型为Z3的边界链码个数，NUM4表示边界链码类型为Z4的边界链码个数，NUM5表示边界链码类型为Z5的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0104左上区欠缺程度特征链码示意图如图7所示，其中，实线代表系统正常加载，虚线代表系统加载变缓。0105由于系统突然急剧加载，使有杆泵抽油井泵示功图曲线的增载线前半部分曲率变化较大，因此，左上区快速加载程度特征向量表示为01060107其。
31、中，NUM1表示边界链码类型为Z1的边界链码个数，NUM2表示边界链码类型为Z2的边界链码个数，NUM3表示边界链码类型为Z3的边界链码个数，NUM4表示边界链码类型为Z4的边界链码个数，NUM5表示边界链码类型为Z5的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0108左上区快速加载程度特征链码示意图如图8所示，其中，实线代表系统正常加载，虚线代表系统急剧加载。0109在上冲程末端由于系统内部相关器件发生碰撞，使载荷突然增大，造成泵示功图曲线的上冲程尾部的曲率变化较大，因此右上区突出程度特征向量表示为01100111其中，NUM3表示边界链码类型为Z3的边界链码个数，NU。
32、M4表示边界链码类型为Z4的边界链码个数，NUM5表示边界链码类型为Z5的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数，为有杆泵抽油井泵示功图的上冲程数据点纵坐标的平均值，为边界链码的类型为Z3、Z4和Z5的边界链码重画后的泵示功图的数据点纵坐标的平均值。0112由于系统卸载提前，使泵示功图曲线的上冲程线后半部分曲率变化较小因此右上区欠缺程度特征向量表示为0113说明书CN104481496A8/11页120114其中，NUM1表示边界链码类型为Z1的边界链码个数，NUM2表示边界链码类型为Z2的边界链码个数，NUM14表示边界链码类型为Z14的边界链码个数，NUM15表示边。
33、界链码类型为Z15的边界链码个数，NUM16表示边界链码类型为Z16的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0115右上区欠缺程度特征链码示意图如图9所示，其中，实线代表系统正常加载，虚线代表系统卸载变缓。0116系统在上冲程结束时快速卸载，使泵示功图曲线的上死点处曲率变化较大，因此右上区快速卸载程度特征向量表示为01170118其中，NUM11表示边界链码类型为Z11的边界链码个数，NUM12表示边界链码类型为Z12的边界链码个数，NUM13表示边界链码类型为Z13的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0119右上区快速卸载程度特征链码示。
34、意图如图10所示，其中，实线代表系统正常加载，虚线代表系统快速卸载。0120在下冲程末端，由于系统内部相关器件发生碰撞，使载荷突然减小，造成泵示功图曲线的下冲程尾部的曲率变化较大，因此左下区突出程度特征向量表示为01210122其中，NUM11表示边界链码类型为Z11的边界链码个数，NUM12表示边界链码类型为Z12的边界链码个数，NUM13表示边界链码类型为Z13的边界链码个数，NUM14表示边界链码类型为Z14的边界链码个数，NUM15表示边界链码类型为Z15的边界链码个数，NUM16表示边界链码类型为Z16的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数，为有杆泵抽油井泵。
35、示功图的下冲程数据点纵坐标的平均值，为边界链码的类型为Z11、Z12、Z13、Z14、Z15和Z16的边界链码重画后的泵示功图的数据点纵坐标的平均值。0123由于系统的加载提前，使泵示功图曲线的下冲程线后半部分曲率变化较小因此左下区欠缺程度特征向量表示为01240125其中，NUM6表示边界链码类型为Z6的边界链码个数，NUM7表示边界链码类型为Z7的边界链码个数，NUM8表示边界链码类型为Z8的边界链码个数，NUM9表示边界链码类型为Z9的边界链码个数，NUM10表示边界链码类型为Z9的边界链码个数，NUM10表示边界链码类型为Z10的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链。
36、码总数。0126左下区欠缺程度特征链码示意图如图11所示，其中，实线代表系统正常加载，虚说明书CN104481496A9/11页13线代表系统加载提前。0127由于系统在卸载时变缓，使泵示功图曲线的卸载线后半部分曲率变化较小，因此右下区欠缺程度特征向量表示为01280129其中，NUM10表示边界链码类型为Z9的边界链码个数，NUM10表示边界链码类型为Z10的边界链码个数，NUM11表示边界链码类型为Z11的边界链码个数，NUM12表示边界链码类型为Z12的边界链码个数，NUM13表示边界链码类型为Z13的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0130右下区欠缺程度。
37、特征链码示意图如图12所示，其中，实线代表系统正常加载，虚线代表系统加载提前。0131由于系统发生急剧卸载，使泵示功图曲线的卸载线的曲率变化较大，因此右下区快速卸载程度特征向量表示为01320133其中，NUM9表示边界链码类型为Z9的边界链码个数，NUM11表示边界链码类型为Z11的边界链码个数，NUM12表示边界链码类型为Z12的边界链码个数，NUM13表示边界链码类型为Z13的边界链码个数，NUM表示边界链码重画后的泵示功图的边界链码总数。0134右下区快速卸载程度特征链码示意图如图13所示，其中，实线代表系统正常加载，虚线代表系统加载提前。0135步骤43将各区域的泵示功图特征向量进行。
38、量化，得到量化后的泵示功图特征向量。0136本发明具体实施方式中得到的一个量化后的泵示功图特征向量如表1所示0137表1量化后的泵示功图特征向量01380139步骤5利用基于粒子群优化的谱聚类算法对量化后的泵示功图特征向量进行聚类，完成有杆泵抽油井故障诊断。0140步骤51设将109个泵示功图特征向量聚类为K类，为尺度参数，其中，0，1。0141步骤52将聚类数K和尺度参数作为粒子进行聚类优化，得到最优的聚类数K和最优的尺度参数。0142步骤521将聚类数K和尺度参数作为粒子，随机产生若干个粒子，初始化粒子的速度和位置。0143步骤522109个量化后的泵示功图特征向量由D1、D2D109表示。
39、，由高斯函数SDP，DQEXP|DPDQ|2/22计算任意两个量化后的泵示功图特征向量的相似度，说明书CN104481496A10/11页14构成相似度矩阵SR109109，其中，P、Q1、2109。0144步骤523将每个量化后的泵示功图特征向量看作顶点V，根据量化后的泵示功图特征向量间的相似度对各顶点间的边E赋权重值W，得到一个基于量化后的泵示功图特征向量相似度的无向加权图GV，E。0145步骤524将相似度矩阵S的每行元素相加，即得到对应顶点的度，以所有顶点的度作为对角元素构成的对角矩阵，即度矩阵D。0146步骤525根据度矩阵、各顶点间的边的权重及单位矩阵，建立拉普拉斯矩阵LRWID1。
40、W，其中，I为单位矩阵。0147步骤526计算拉普拉斯矩阵LRW前K个最小的特征值对应的特征向量U1、U2UK，K个特征向量组成特征向量矩阵UR109K。0148步骤527将特征向量矩阵UR109K的第G行向量EGG1、2109作为聚类数据集。0149步骤528用CRITC指标函数的倒数作为各个粒子的适应度值。0150CRITC指标函数的定义如下01510152其中01530154CF，O表示聚类数据集EGG1、2109中每个聚类中心CF和整个聚类数据集EGG1、2109的中心0之间的距离；0155CF由每个聚类中所有数据的平均值求得；01560由聚类数据集EGG1、2109中所有数据的平均值。
41、求得；0157|CF|表示CF，O所对应的权值，等于每个聚类中数据的数量；0158CF，O表示聚类数据集EGG1、2109中每个数据点和它所在聚类的中心CF之间的距离；0159A为常数系数，本实施方式中取A08；0160CRITC指标函数在0，1区间内取值。0161步骤529若当前粒子的适应度值小于其历史适应度值，则将当前粒子的适应度值作为其历史适应度值。0162步骤5210设定粒子的适应度值终止条件，若当前粒子的适应度值小于其终止条件，则当前的聚类数K和尺度参数为最优的聚类数K和最优的尺度参数，当前的聚类数据集EGG1、2109为优化后的聚类数据集EGG1、2109，否则，更新粒子的速度和位。
42、置，返回步骤523。0163本实施方式中，设定代表聚类数K和尺度参数的粒子的适应度值终止条件为05，聚类数K和尺度参数的最优值为5和01237。0164步骤53利用最优的聚类数K和最优的尺度参数进行聚类优化，得到最优的聚类中心。说明书CN104481496A11/11页150165步骤531将优化后的聚类数据集EGG1、2109中数据作为粒子，随机选取K5个粒子作为初始聚类中心，初始化粒子的速度和位置。0166步骤532用JC指标函数作为粒子的适应度值。0167JC指标函数的定义如下01680169其中NC表示优化后的聚类数据集EGG1、2109的聚类中心，NC1、2、5。0170步骤533若。
43、当前粒子的适应度值小于其历史适应度值时，则把当前粒子的适应度值作为其历史适应度值，当前粒子作为最优粒子，同时将当前粒子位置记为其历史最优位置。0171步骤534设定粒子的适应度值终止条件，若当前粒子的适应度值小于其终止条件，则当前最优粒子的聚类中心为最优的聚类中心，否则，更新粒子的速度和位置，返回步骤532。0172步骤54根据最优的聚类中心将109个泵示功图特征向量聚类为5类，得到最终的聚类结果H1、H2H5，若待诊断故障类型的泵示功图特征向量与某知故障类型的泵示功图特征向量聚为一类，则两者为同类故障，即得到有杆泵抽油井故障诊断结果，若待诊断故障类型的泵示功图特征向量单独聚为一类，则该待诊断。
44、故障类型为新故障类型。0173如图14所示，为本发明具体实施方式中对量化后的泵示功图特征向量进行聚类为5类后的示意图。0174本实施方式中，待诊断故障类型的量化后的泵示功图特征向量由“”标记和“正常”类型的量化后的泵示功图特征向量聚类为同一类型，因此该待诊断故障类型为“正常”。说明书CN104481496A1/9页16图1说明书附图CN104481496A2/9页17图2说明书附图CN104481496A3/9页18图3说明书附图CN104481496A4/9页19图4图5说明书附图CN104481496A5/9页20图6图7说明书附图CN104481496A6/9页21图8图9说明书附图CN104481496A7/9页22图10图11说明书附图CN104481496A8/9页23图12图13说明书附图CN104481496A9/9页24图14说明书附图CN104481496A。