CN102011576A - 基于示功图的有杆抽油***故障递阶诊断方法 - Google Patents
基于示功图的有杆抽油***故障递阶诊断方法 Download PDFInfo
- Publication number
- CN102011576A CN102011576A CN2010105572348A CN201010557234A CN102011576A CN 102011576 A CN102011576 A CN 102011576A CN 2010105572348 A CN2010105572348 A CN 2010105572348A CN 201010557234 A CN201010557234 A CN 201010557234A CN 102011576 A CN102011576 A CN 102011576A
- Authority
- CN
- China
- Prior art keywords
- fault
- sample
- distribution
- interval
- load
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005086 pumping Methods 0.000 title claims abstract description 44
- 238000010586 diagram Methods 0.000 title abstract description 5
- 238000009826 distribution Methods 0.000 claims description 66
- 238000012360 testing method Methods 0.000 claims description 52
- 238000006073 displacement reaction Methods 0.000 claims description 24
- 238000002405 diagnostic procedure Methods 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 9
- 230000005856 abnormality Effects 0.000 claims description 7
- 230000007812 deficiency Effects 0.000 claims description 7
- 238000003556 assay Methods 0.000 claims description 6
- 239000007789 gas Substances 0.000 claims description 6
- 210000002445 nipple Anatomy 0.000 claims description 6
- 239000007921 spray Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000012353 t test Methods 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 abstract description 15
- 239000003129 oil well Substances 0.000 description 14
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000012423 maintenance Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公布了一种基于示功图的有杆抽油***故障递阶诊断法,分为故障分辨和故障识别两个阶段,即首先在故障分辨阶段根据正常样本的统计规律把示功图分为故障类和非故障类;然后采用基于示功图的有杆抽油***故障识别的搜索树方法,对故障样本进行故障类型的识别。它提高了故障诊断***的故障诊断性能,且不需要建立和求解有杆抽油***力学模型,也不存在训练集问题,并能反映出有杆抽油***的自身的基本特征。
Description
技术领域
本发明涉及一种基于示功图的有杆抽油***故障诊断方法,尤其涉及一种基于示功图的有杆抽油***故障递阶诊断方法,属于有杆抽油***故障诊断方法的技术领域。
背景技术
有杆抽油***工作环境恶劣,发生故障的概率很高,不及时诊断出来,就会造成能源的浪费,影响生产,给企业、国家造成损失。
目前有杆抽油***的故障诊断主要是先根据地面示功图,通过有杆抽油***力学模型得到泵功图,再以正常样本和尽可能多类型故障样本的泵功图作为训练集,训练出相应模型后再对各样本进行故障诊断。这类方法主要存在两个方面的问题:一是有杆抽油***是一个机电液耦合的复杂非线性***,边界条件和阻尼系数很难准确确定,进而很难精确建立和求解该***的力学模型;二是有杆抽油***具有样本量大、平稳状态多、类别不均衡的特点,使训练集中缺少很多类型的故障样本,尤其是严重故障的样本,而严重故障的准确、实时诊断对于有杆抽油***意义重大。由于有杆抽油***的复杂性,至今还没有一个完全有效的手段来对井下故障进行实时、准确的故障诊断。
反映抽油***悬点载荷随其位移变化规律的图形称为地面示功图或光杆示功图,简称示功图,它是示功仪在抽油机一个抽吸周期内测取的封闭曲线。地面示功图是抽油井采油现场采集的第一手资料,基于地面示功图来实现有杆抽油***的实时诊断对解除油井故障、保证油井正常生产或提高油井产量等就显得非常重要而且具有现实意义。
每组示功图数据包括载荷和位移,共216对数据。油田一个小时采集一次示功图,一天就有24组数据,一个星期有168组数据,一个月(30天计算)有720组数据。所以,油井示功图数据多、样本量大。并且,对于油田来讲,绝大多数都是平稳的工作状态,包括两种情况:一是油井正常,大部分示功图均正常,故障样本和故障类型较少;二是有的油井以某种故障,如含气、出砂、供液不足、轻微碰挂、振动影响等某类对生产影响不太大的故障为主,出现其它故障的概率较低。以上这两种情况不需要立即停产进行检修,可以待下次修井时或在方便和必要时处理,也就是说大量的油井示功图数据都是不需要立即进行故障处理的。而一旦出现凡尔失灵、抽油杆断脱、泵卡死等严重故障,均需立即停产,及时采取检泵和修井措施,以减少损失。油井修井措施,特别是对油井大修或更换设备会对油井的生产状况和功图产生影响,检修前后的生产状况和功图往往会有较大差别,检修后的训练集和样本特征往往需要重新形成。另外,现有的采油工业中抽油机数量大、分布广,这更增加了井下状况的复杂性,使各口井均有自身特点,可能出现的故障类型也不固定。
发明内容
技术问题:本发明目的是针对有杆抽油***故障诊断目前存在的缺陷和油井示功图样本量大、平稳状态多、类别不均衡(尤其缺少严重故障的样本)的特点,以正常(或平稳状态)样本地面示功图作为切入点,提供一种基于示功图的有杆抽油***故障递阶诊断法,即首先在故障分辨阶段根据正常(或平稳状态)样本的统计规律把示功图分为故障类和非故障类;然后在故障诊断阶段针对故障样本进行详细的故障诊断识别。由于大部分的示功图都是正常(或平稳)的,先进行故障分辨就节省了大量的时间,从而提高油田管理的效率。另外,故障分辨的训练样本来自于正常(或平稳状态)样本,样本丰富且能反映出有杆抽油***的自身的基本特征。
本发明为实现上述目的,采用如下技术方案:
本发明基于示功图的有杆抽油***故障递阶诊断方法包括如下步骤:
1)故障分辨阶段,故障分辨阶段又分为两个阶段:
11)训练阶段:在训练阶段由正常或平稳样本提取15个示功图上的几何特征量,根据这些样本的特征量得到相关统计信息:
111)人工选取n个正常或平稳样本,进行15个特征量的提取,其中n≥50;
112)对提取的特征量进行异常值检验,如果异常,则剔除该样本,并重新进行异常值检验,直至没有异常值;
113)对异常值检验后的样本进行分布规律检验;
114)根据分布规律检验结果计算正常区域和故障区域;
12)分辨阶段:根据测试样本本身是否存在打结点和测试样本是否有特征量落入故障区域来判断是否故障:
121)判断测试样本本身是否存在打结点,如果存在,则为故障样本,跳转步骤2);否则跳转步骤122);
122)提取测试样本的15个特征量;
123)测试样本是否有特征量落入故障区域,如果有,则为故障样本,跳转步骤2);否则跳转步骤124);
124)判定测试样本为正常样本,不需要进行故障识别,结束;
2)故障识别阶段:采用基于示功图的有杆抽油***故障识别的搜索树方法,对故障样本进行故障类型的识别。
优选地,所述15个特征量为下死点(E点)载荷、上死点(F点)位移和载荷、固定凡尔打开点(B点)的位移和载荷、游动凡尔打开点(D点)的位移和载荷、面积、最大载荷、最小载荷、载荷总的平均变化量,EB、BF、FD、DE段载荷的平均变化量。
优选地,步骤112)所述剔除异常样本的方法采用t检验准则剔除异常数据法:
21)在n个观测值中找出与平均值相比误差最大的正常样本的特征量值作为可疑值xk,其中观测值为正常样本的一个特征量;
22)对不包括可疑值xk在内的n-1个观测值,计算平均值:
以及标准差估计值:
23)确定风险率αT,αT=0.001,从t分布表上查出AT(αT;n-2),再计算得到:
(3)
24)检查xk,如果有
成立,则xk为异常值,应予剔除,跳转步骤25);否则xk不是异常值,不能剔除,并终止检查;
25)若xk为异常值,将它剔除后,对余下的n-1个观测值重复上述步骤,直到不再有异常值。
优选地,步骤113)所述对样本进行分布规律检验为在对样本分布进行假设基础上,采用χ2拟合优度检验法检验假设分布与实际分布是否相符合:
31)提出假设H0:X服从某种备选的分布类型;选取6个常用的连续性随机分布(分别为:指数分布、均匀分布、威布尔分布、正态分布、瑞利分布和伽马分布)为备选的分布类型,分别进行假设检验,所述分布类型的分布函数中含有r个未知数,r为自然数;
32)将实数轴均分为k个不相交的区间(a0,a1],(a1,a2],...,(ak-1,ak],其中:当n≤200时,k分别取5、7、...13;当n>200时,k分别取11、13、...21;a0,ak按式(5)循环计算:
a0=xmin-(xmax-xmin)×j×1%,j=1,2,L 20
(5)
ak=xmax+(xmax-xmin)×j×1%,j=1,2,L 20
33)计算数据落入各区间的频数ni,i=1,2,3...k;
34)在H0成立的条件下,计算X落入各区间的概率pi:
pi=P(ai-1<X≤ai) (6)
进而得到理论频数npi(i=1,2,Λ,k);
优选地,步骤114)所述特征量正常区域和故障区域计算方法如下:
41)对各特征量选取统一的概率值α,α为样本落入正常区域区间[b0,a0]的概率,α=99.99%,并且样本落入故障区域区间(-∞,b0)和故障区域区间(a0,+∞)的概率相同,则
P{-∞<X<b0}+P{b0≤X≤a0}+P{a0<X<+∞}=1 (7)
P{-∞<X<b0}=P{a0<X<+∞} (8)
所以
42)根据式(9)、(10)及该特征量的分布类型与参数,计算出a0、b0,即可得到正常区域区间[b0,a0]与故障区域区间(-∞,b0)∪(a0,+∞);对于载荷平均变化量,有意义的故障区域为(a0,+∞);
43)如果某特征量的分布规律检验结果是符合多个分布,则分别按多个分布规律计算正常区域区间,并取它们的并集为该特征量正常区域区间。
45)令b1=b2·λ(a2-b2),a1=a2·λ(a2-b2),λ为设定的阈值,λ=0.2,其中b2、a2分别为样本的最小值与最大值。区间[b1,a1]与步骤41)至44)计算得到的正常区域区间的并集为该特征量最终的正常区域区间[b,a],则该特征量最终的故障区域区间(-∞,b)∪(a,+∞)。
优选地,步骤121)所述判断测试样本本身是否存在打结点的方法为:
51)计算测试样本位移的最大和最小值;
52)将测试样本示功图上的位移均分为J份,J=10000,j=1,2,L J-1;
53)计算第j条直线xj=xmin+(xmax-xmin)×j/J与示功图的交点集合U{H1,H2};
54)用插值法求出H1对应的载荷yj1;
55)用插值法求出H2对应的载荷yj2;
56)若|yj1-yj2|/yj1<ε,ε为设定的阈值,ε=0.001,则存在打结点,结束;
57)如果j<J,则更新j=j+1,跳转步骤53),否则判定测试样本不存在打结点,结束。
优选地,步骤2)所述搜索树分为四大类:面积变化类、打结类、振动类、碰挂类,其中面积变化类又分为三类:面积增大类、面积极小类、面积偏小或正常类。面积类故障搜索树用于识别的故障类型包括:固定阀漏失、抽油杆中上部断脱、游动阀失灵、油管底部漏失、连抽带喷、抽油杆底部断脱、气锁、泵卡死、严重供液不足、油管漏失、游动阀漏失、游动阀关闭迟缓、柱塞脱出工作筒、供液不足、气体影响、惯性载荷大;打结类故障搜索树用于识别的故障类型包括:下碰泵、二级振动、柱塞脱出工作筒;振动类故障搜索树用于识别的故障类型包括:出砂、振动过大;碰挂类故障搜索树用于识别的故障类型包括:上碰挂。
本发明方法中,步骤23)的αT、步骤36)的αχ、步骤41)的α、步骤45)的λ、步骤56)的ε、图5的符号说明中ε1、ε2、ε3、ξ均为设定的阈值,步骤32)至34)的k和步骤52)的J均为设定的初值,这些阈值和初值主要根据精度要求及计算经验确定,本发明公开了其优选的数值。
有益效果:本发明基于示功图的有杆抽油***故障递阶诊断方法在保持计算简单的前提下,提高了分类的精度,尤其是针对小样本类标的异常事件,从而提高了故障诊断***的故障诊断性能;它不需要建立和求解有杆抽油***力学模型,也不存在训练集问题;先进行故障分辨可节省大量时间,提高油田管理的效率。实例表明故障分辨能从大量的数据中筛选出故障样本,分辨的正确率比较高;训练样本的统计规律能反映油井实时的生产状况。通过实例和多种典型故障示功图,包括严重故障,验证了该方法能准确识别出有杆抽油***的单个故障和组合故障。
附图说明
图1:基于示功图的有杆抽油***故障递阶诊断方法的流程图
图2:训练阶段流程图;
图3:正常区域和故障区域概率示意图;
图4:分辨阶段流程图;
图5:面积类故障搜索树;
图5的符号说明如下:
T1-1:面积不大于ε1b面积(设面积的正常区域为[b面积,a面积]),ε1为小于1的小值,ε1为设定的阈值,ε1=0.25
T1-2:面积小于a面积且大于ε1b面积
T1-3:面积不小于a面积
T2-1:(最大载荷-最小载荷)/(Z最大载荷-Z最小载荷)小于ε2,ε2为小于1的小值,ε2为设定的阈值,ε2=0.01
(设最大载荷的正常区域为[b最大载荷,a最大载荷],最小载荷的正常区域为[b最 小载荷,a最小载荷],Z最大载荷=0.5(a最大载荷-b最大载荷),Z最小载荷=0.5(a最小载荷-b最小载荷))
T2-2:(最大载荷-最小载荷)/(Z最大载荷-Z最小载荷)不小于ε2
T3-1:最小载荷不小于ξa最小载荷,ξ大于1,ξ为设定的阈值,ξ=1.2
T3-2:最小载荷不大于ε3b最小载荷,ε3为小于1的小值,ε3为设定的阈值,ε3=0.8
T3-3:最小载荷小于ξa最小载荷且大于ε3b最小载荷
T3-4:B点不正常
T3-5:B点正常
T4-1:FD是凸函数
T4-2:FD不是凸函数
T5-1:B点或F点位置异常,D点正常
T5-2:D点位置异常、B点正常
T5-3:B、D点位置异常
T6-1:示功图是平行四边形
T6-2:示功图不是平行四边形
T7-1:F点不异常
T7-2:F点异常
T8-1:EB是凹函数
T8-2:EB不是凹函数
T9-1:BF是凹函数
T9-2:BF不是凹函数
T10-1:FD是凹函数
T10-1:FD是凹函数
T11-1:FD不是圆弧
T11-2:FD是圆弧
T12-1:示功图是平行四边形
T12-2:示功图不是平行四边形
T13-1:F点载荷不小于bF点载荷(设F点载荷的正常区域为[bF点载荷,aF点载荷])
T13-2:F点载荷小于bF点载荷
T14-1:EB不是凹函数
T14-2:EB是凹函数
T15-1:FD是凹函数
T15-2:FD不是凹函数
T16-1:FD是圆弧
T16-2:FD不是圆弧
T17-1:FD是凹函数
T17-2:FD不是凹函数
T18-1:FD是圆弧
T18-2:FD不是圆弧
T19-1:BF是凹函数
T19-2:BF不是凹函数
T20-1:FD是凹函数
T20-2:FD不是凹函数
T21-1:FD是圆弧
T21-2:FD不是圆弧
T22-1:FD是凹函数
T22-2:FD不是凹函数
T23-1:FD是圆弧
T23-2:FD不是圆弧
图6:振动类故障搜索树;
图6的符号说明如下:
ZT1-1:5个载荷变化量均进入故障区域
ZT1-2:5个载荷变化量未全部进入故障区域
ZT2-1:BF间有振动
ZT2-2:BF间没有振动
ZT2-3:BF间有振动
ZT2-4:BF间没有振动
图7:打结类故障搜索树;
图7的符号说明如下:
DT1:打结点的位移小于总位移的1/4
DT2:打结点的位移大于总位移的3/4并且F点的载荷大于载荷中值
DT3:打结点的位移大于总位移的3/4并且F点的载荷小于载荷中值
图8:碰挂类故障搜索树;
图8符号说明如下:
P1:最大载荷大于a最大载荷且最大载荷点在F点附近(设最大载荷的正常区域为[b最大载荷,a最大载荷])
P2:不满足P1规则。
图5至图8符号说明中相关名词定义如下:
由于抽油杆的力学特性及油田的现场测试仪器的精密性等条件的限制,所述的凹函数、圆弧、平行四边形等并不是数学意义上严格的定义,本发明采用的定义如下:
定义1由地面示功图上E、B、F、D四点坐标计算出的四点连线的斜率,如果EB和FD、BF和DE的斜率差在一个极小值范围内,则称该示功图是平行四边形;
定义2设f(x)为由多点拟合而成的曲线,将f(x)的两端点连一直线,曲线中的点距离直线的距离最远的五个点都在下侧,则称f(x)是凸函数;
定义3设f(x)为由多点拟合而成的曲线,将f(x)的两端点连一直线,曲线中的点距离直线的距离最远的五个点都在上侧,则称f(x)是凹函数;
定义4设f(x)为由多点拟合而成的曲线,将f(x)的两端点I、J连一直线,曲线中的点距离该直线距离最远的点记为H,IH和HJ均为凸函数,则称f(x)是圆弧;
定义5由地面示功图上B、F间的点拟合成一条直线,若B、F间距离该直线最远的三个点,位移是依次递增,而距离是依次递减的,则称BF间有振动。
具体实施方式
本发明针对有杆抽油***的特点及其故障诊断目前存在的问题,提供一种基于示功图的有杆抽油***故障递阶诊断法,分为故障分辨和故障识别两个阶段。即首先在故障分辨阶段根据正常样本的统计规律把示功图分为故障类和非故障类;然后将统计理论与搜索树相结合,针对故障样本进行详细的故障类型识别。它不需要建立和求解有杆抽油***力学模型,也不存在训练集问题。
所提供的有杆抽油***故障分辨方法,分为训练和分辨两个阶段。在训练阶段,t检验剔除异常数据后,通过χ2拟合优度检验确定样本的随机分布形式及其参数,然后,计算出正常区域和故障区域;在分辨阶段,根据测试样本是否有特征量落入故障区域来判断其分类。故障分辨能从大量的数据中筛选出故障样本,分辨的正确率比较高;训练样本的统计规律能反映油井实时的生产状况。
所提供的有杆抽油***故障识别的搜索树法,在故障分辨的基础上,将正常样本的统计信息与搜索树相结合,建立起基于示功图的有杆抽油***故障识别的搜索树,对故障样本进行故障类型的详细识别。该方法能准确识别出有杆抽油***的单个故障和组合故障。
如图1,本发明基于示功图的有杆抽油***故障递阶诊断方法包括如下步骤:
1)故障分辨阶段,故障分辨阶段又分为两个阶段:
11)训练阶段:在训练阶段由正常或平稳样本提取15个示功图上的几何特征量,根据这些样本的特征量得到相关统计信息,如图2:
111)人工选取n个正常或平稳样本,进行15个特征量的提取,其中n≥50;
112)对提取的特征量进行异常值检验,如果异常,则剔除该样本,并重新进行异常值检验,直至没有异常值;
113)对异常值检验后的样本进行分布规律检验;
114)根据分布规律检验结果计算正常区域和故障区域;
12)分辨阶段:根据测试样本本身是否存在打结点和测试样本是否有特征量落入故障区域来判断是否故障,如图4:
121)判断测试样本本身是否存在打结点,如果存在,则为故障样本,跳转步骤2);否则跳转步骤122);
122)提取测试样本的15个特征量;
123)测试样本是否有特征量落入故障区域,如果有,则为故障样本,跳转步骤2);否则跳转步骤124);
124)判定测试样本为正常样本,不需要进行故障识别,结束;
2)故障识别阶段:采用基于示功图的有杆抽油***故障识别的搜索树方法,对故障样本进行故障类型的识别。
优选地,步骤11)、步骤11)、步骤122)所述15个特征量为下死点(E点)载荷、上死点(F点)位移和载荷、固定凡尔打开点(B点)的位移和载荷、游动凡尔打开点(D点)的位移和载荷、面积、最大载荷、最小载荷、载荷总的平均变化量,EB、BF、FD、DE段载荷的平均变化量。
优选地,步骤112)所述剔除异常样本的方法采用t检验准则剔除异常数据法:
21)在n个观测值中找出与平均值相比误差最大的正常样本的特征量值作为可疑值xk,其中观测值为正常样本的一个特征量;
22)对不包括可疑值xk在内的n-1个观测值,计算平均值:
以及标准差估计值:
23)确定风险率αT,αT=0.001,从t分布表上查出AT(αT;n-2),再计算得到:
(3)
24)检查xk,如果有
成立,则xk为异常值,应予剔除,跳转步骤25);否则xk不是异常值,不能剔除,并终止检查;
25)若xk为异常值,将它剔除后,对余下的n-1个观测值重复上述步骤,直到不再有异常值。
优选地,步骤113)所述对样本进行分布规律检验为在对样本分布进行假设基础上,采用χ2拟合优度检验法检验假设分布与实际分布是否相符合:
31)提出假设H0:X服从某种备选的分布类型;选取6个常用的连续性随机分布(分别为:指数分布、均匀分布、威布尔分布、正态分布、瑞利分布和伽马分布)为备选的分布类型,分别进行假设检验,所述分布类型的分布函数中含有r个未知数,r为自然数;
32)将实数轴均分为k个不相交的区间(a0,a1],(a1,a2],...,(ak-1,ak],其中:当n≤200时,k分别取5、7、...13;当n>200时,k分别取11、13、...21;a0,ak按式(5)循环计算:
a0=xmin-(xmax-xmin)×j×1%,j=1,2,L 20
(5)
ak=xmax+(xmax-xmin)×j×1%,j=1,2,L 20
33)计算数据落入各区间的频数ni,i=1,2,3...k;
34)在H0成立的条件下,计算X落入各区间的概率pi:
pi=P(ai-1<X≤ai) (6)
进而得到理论频数npi(i=1,2,Λ,k);
37)若则拒绝H0,否则,可接受H0。
优选地,步骤114)所述特征量正常区域和故障区域计算方法如下,如图3:
41)对各特征量选取统一的概率值α,α为样本落入正常区域区间[b0,a0]的概率,α=99.99%,并且样本落入故障区域区间(-∞,b0)和故障区域区间(a0,+∞)的概率相同,则
P{-∞<X<b0}+P{b0≤X≤a0}+P{a0<X<+∞}=1 (7)
P{-∞<X<b0}=P{a0<X<+∞} (8)
所以
42)根据式(9)、(10)及该特征量的分布类型与参数,计算出a0、b0,即可得到正常区域区间[b0,a0]与故障区域区间(-∞,b0)∪(a0,+∞);对于载荷平均变化量,有意义的故障区域为(a0,+∞);
43)如果某特征量的分布规律检验结果是符合多个分布,则分别按多个分布规律计算正常区域区间,并取它们的并集为该特征量正常区域区间。
根据切比雪夫(Chebyshev)不等式:
或
对任意的特征量分布规律,若正常区域区间取为[μ-3σ,μ+3σ],则样本落入正常区域区间的概率α≥88.89%;若正常区域区间取为[μ-4σ,μ+4σ],则α≥93.75%;若正常区域区间取为[μ-5σ,μ+5σ],则α≥96%。在此情况下,考虑到故障的特征量一般偏离正常区域区间较远,本发明取正常区域区间为[μ-5σ,μ+5σ]。
45)令b1=b2·λ(a2-b2),a1=a2·λ(a2-b2),λ为设定的阈值,λ=0.2,其中b2、a2分别为样本的最小值与最大值。区间[b1,a1]与步骤41)至44)计算得到的正常区域区间的并集为该特征量最终的正常区域区间[b,a],则该特征量最终的故障区域区间(-∞,b)∪(a,+∞)。
优选地,步骤121)所述判断测试样本本身是否存在打结点的方法为:
51)计算测试样本位移的最大和最小值;
52)将测试样本示功图上的位移均分为J份,J=10000,j=1,2,L J-1;
53)计算第j条直线xj=xmin+(xmax-xmin)×j/J与示功图的交点集合U{H1,H2};
54)用插值法求出H1对应的载荷yj1;
55)用插值法求出H2对应的载荷yj2;
56)若|yj1-yj2|/yj1<ε,ε为设定的阈值,ε=0.001,则存在打结点,结束;
57)如果j<J,则更新j=j+1,跳转步骤53),否则判定测试样本不存在打结点,结束。
优选地,如图5-图8,步骤2)所述搜索树分为四大类:面积变化类、打结类、振动类、碰挂类,其中面积变化类又分为三类:面积增大类、面积极小类、面积偏小或正常类。面积类故障搜索树用于识别的故障类型包括:固定阀漏失、抽油杆中上部断脱、游动阀失灵、油管底部漏失、连抽带喷、抽油杆底部断脱、气锁、泵卡死、严重供液不足、油管漏失、游动阀漏失、游动阀关闭迟缓、柱塞脱出工作筒、供液不足、气体影响、惯性载荷大;打结类故障搜索树用于识别的故障类型包括:下碰泵、二级振动、柱塞脱出工作筒;振动类故障搜索树用于识别的故障类型包括:出砂、振动过大;碰挂类故障搜索树用于识别的故障类型包括:上碰挂。
由于油稠、油井结蜡、管线堵摩阻大三种故障的示功图特征类似,在搜索树中未将其进一步分类,需根据油井实时生产状况进行判断(可参见相应的故障机理分析)。游动阀失灵、油管底部漏失、连抽带喷、抽油杆底部断脱四种故障的示功图也特征类似,且这四种故障均是极其严重的故障,均需立即停止生产进行处理,进一步区分意义不大,因而在搜索树中也未将其进一步分类。
Claims (7)
1.一种基于示功图的有杆抽油***故障递阶诊断方法,其特征在于包括如下步骤:
1)故障分辨阶段,故障分辨阶段又分为两个阶段:
11)训练阶段:在训练阶段由正常或平稳样本提取15个示功图上的几何特征量,根据这些样本的特征量得到相关统计信息:
111)人工选取n个正常或平稳样本,进行15个特征量的提取,其中n≥50;
112)对提取的特征量进行异常值检验,如果异常,则剔除该样本,并重新进行异常值检验,直至没有异常值;
113)对异常值检验后的样本进行分布规律检验;
114)根据分布规律检验结果计算正常区域和故障区域;
12)分辨阶段:根据测试样本本身是否存在打结点和测试样本是否有特征量落入故障区域来判断是否故障:
121)判断测试样本本身是否存在打结点,如果存在,则为故障样本,跳转步骤2);否则跳转步骤122);
122)提取测试样本的15个特征量;
123)测试样本是否有特征量落入故障区域,如果有,则为故障样本,跳转步骤2);否则跳转步骤124);
124)判定测试样本为正常样本,不需要进行故障识别,结束;
2)故障识别阶段:采用基于示功图的有杆抽油***故障识别的搜索树方法,对故障样本进行故障类型的识别。
2.根据权力要求1所述的基于示功图的有杆抽油***故障递阶诊断方法,其特征在于所述15个特征量为下死点(E点)载荷、上死点(F点)位移和载荷、固定凡尔打开点(B点)的位移和载荷、游动凡尔打开点(D点)的位移和载荷、面积、最大载荷、最小载荷、载荷总的平均变化量,EB、BF、FD、DE段载荷的平均变化量。
3.根据权力要求1所述的基于示功图的有杆抽油***故障递阶诊断方法,其特征在于步骤112)所述剔除异常样本的方法采用t检验准则剔除异常数据法:
21)在n个观测值中找出与平均值相比误差最大的正常样本的特征量值作为可疑值xk,其中观测值为正常样本的一个特征量;
22)对不包括可疑值xk在内的n-1个观测值,计算平均值:
以及标准差估计值:
23)确定风险率αT,αT=0.001,从t分布表上查出AT(αT;n-2),再计算得到:
(3)
24)检查xk,如果有
成立,则xk为异常值,应予剔除,跳转步骤25);否则xk不是异常值,不能剔除,并终止检查;
25)若xk为异常值,将它剔除后,对余下的n-1个观测值重复上述步骤,直到不再有异常值。
4.根据权力要求1所述的基于示功图的有杆抽油***故障递阶诊断方法,其特征在于步骤113)所述对样本进行分布规律检验为在对样本分布进行假设基础上,采用χ2拟合优度检验法检验假设分布与实际分布是否相符合:
31)提出H0:X服从某种备选的分布类型;选取6个常用的连续性随机分布:指数分布、均匀分布、威布尔分布、正态分布、瑞利分布和伽马分布为备选的分布类型,分别进行假设检验,所述分布类型的分布函数中含有r个未知数,r为自然数;
32)将实数轴均分为k个不相交的区间(a0,a1],(a1,a2],...,(ak-1,ak],其中:当n≤200时,k分别取5、7、...13;当n>200时,k分别取11、13、...21;a0,ak按式(5)循环计算:
a0=xmin-(xmax-xmin)×j×1%,j=1,2,L 20
(5)
ak=xmax+(xmax-xmin)×j×1%,j=1,2,L 20
33)计算数据落入各区间的频数ni,i=1,2,3...k;
34)在H0成立的条件下,计算X落入各区间的概率pi:
pi=P(ai-1<X≤ai) (6)
进而得到理论频数npi(i=1,2,Λ,k);
37)若则拒绝H0,否则,可接受H0。
5.根据权力要求1所述的基于示功图的有杆抽油***故障递阶诊断方法,其特征在于步骤114)所述特征量正常区域和故障区域计算方法如下:
41)对各特征量选取统一的概率值α,α为样本落入正常区域区间[b0,a0]的概率,α=99.99%,并且样本落入故障区域区间(-∞,b0)和故障区域区间(a0,+∞)的概率相同,则
P{-∞<X<b0}+P{b0≤X≤a0}+P{a0<X<+∞}=1 (7)
P{-∞<X<b0}=P{a0<X<+∞} (8)
所以
(10)
42)根据式(9)、(10)及该特征量的分布类型与参数,计算出a0、b0,即可得到正常区域区间[b0,a0]与故障区域区间(-∞,b0)∪(a0,+∞);对于载荷平均变化量,有意义的故障区域为(a0,+∞);
43)如果某特征量的分布规律检验结果是符合多个分布,则分别按多个分布规律计算正常区域区间,并取它们的并集为该特征量正常区域区间;
45)令b1=b2-λ(a2-b2),a1=a2+λ(a2-b2),λ为设定的阈值,λ=0.2,其中b2、a2分别为异常值检验后的训练样本的最小值与最大值;区间[b1,a1]与步骤41)至44)计算得到的正常区域区间的并集为该特征量最终的正常区域区间[b,a],则该特征量最终的故障区域区间(-∞,b)∪(a,+∞)。
6.根据权力要求1所述的基于示功图的有杆抽油***故障递阶诊断方法,其特征在于步骤121)所述判断测试样本本身是否存在打结点的方法为:
51)计算测试样本位移的最大和最小值;
52)将测试样本示功图上的位移均分为J份,J=10000,j=1,2,L J-1;
53)计算第j条直线xj=xmin+(xmax-xmin)×j/J与示功图的交点集合U{H1,H2};
54)用插值法求出H1对应的载荷yj1;
55)用插值法求出H2对应的载荷yj2;
56)若|yj1-yj2|/yj1<ε,ε为设定的阈值,ε=0.001,则存在打结点,结束;
57)如果j<J,则更新j=j+1,跳转步骤53),否则判定测试样本不存在打结点,结束。
7.根据权力要求1所述的基于示功图的有杆抽油***故障递阶诊断方法,其特征在于步骤2)所述搜索树分为四大类:面积变化类、打结类、振动类、碰挂类,其中面积变化类又分为三类:面积增大类、面积极小类、面积偏小或正常类。面积类故障搜索树用于识别的故障类型包括:固定阀漏失、抽油杆中上部断脱、游动阀失灵、油管底部漏失、连抽带喷、抽油杆底部断脱、气锁、泵卡死、严重供液不足、油管漏失、游动阀漏失、游动阀关闭迟缓、柱塞脱出工作筒、供液不足、气体影响、惯性载荷大;打结类故障搜索树用于识别的故障类型包括:下碰泵、二级振动、柱塞脱出工作筒;振动类故障搜索树用于识别的故障类型包括:出砂、振动过大;碰挂类故障搜索树用于识别的故障类型包括:上碰挂。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010557234 CN102011576B (zh) | 2010-11-24 | 2010-11-24 | 基于示功图的有杆抽油***故障递阶诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010557234 CN102011576B (zh) | 2010-11-24 | 2010-11-24 | 基于示功图的有杆抽油***故障递阶诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102011576A true CN102011576A (zh) | 2011-04-13 |
CN102011576B CN102011576B (zh) | 2013-09-25 |
Family
ID=43841885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010557234 Expired - Fee Related CN102011576B (zh) | 2010-11-24 | 2010-11-24 | 基于示功图的有杆抽油***故障递阶诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102011576B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184414A (zh) * | 2011-05-16 | 2011-09-14 | 中国石油天然气股份有限公司 | 泵示功图的识别和判断方法及其*** |
CN103541723A (zh) * | 2013-11-12 | 2014-01-29 | 丁涛 | 基于地面示功图面积变化的抽油机井实时工况诊断方法 |
CN103806893A (zh) * | 2012-11-14 | 2014-05-21 | 中国石油天然气股份有限公司 | 智能检测方法、装置及检测*** |
CN104110251A (zh) * | 2014-06-24 | 2014-10-22 | 安徽多杰电气有限公司 | 一种基于art2的抽油机示功图识别方法 |
CN104948436A (zh) * | 2014-03-25 | 2015-09-30 | 布里斯托尔D/B/A远程自动化解决方案公司 | 用于确定井下泵的产量的方法和装置 |
CN105257277A (zh) * | 2015-05-15 | 2016-01-20 | 渤海大学 | 基于多变量灰色模型的有杆泵抽油井井下故障预测方法 |
CN105672988A (zh) * | 2015-12-30 | 2016-06-15 | 中国石油天然气股份有限公司 | 抽油机示功图诊断***和诊断方法 |
CN106089184A (zh) * | 2016-07-18 | 2016-11-09 | 中国石油天然气股份有限公司 | 一种井下抽油泵工况的诊断方法和装置 |
CN106121622A (zh) * | 2016-07-27 | 2016-11-16 | 渤海大学 | 一种基于示功图的有杆泵抽油井的多故障诊断方法 |
CN106250692A (zh) * | 2016-07-29 | 2016-12-21 | 广东工业大学 | 一种心率数据质量评估方法 |
CN106682815A (zh) * | 2016-11-30 | 2017-05-17 | 西安石油大学 | 一种抽油井监控***及监控方法 |
WO2020019681A1 (zh) * | 2018-07-25 | 2020-01-30 | 北京国双科技有限公司 | 采油设备故障诊断方法及装置 |
CN110821478A (zh) * | 2018-08-13 | 2020-02-21 | 中国石油天然气股份有限公司 | 油井抽油泵漏失检测方法及装置 |
CN111472726A (zh) * | 2020-03-27 | 2020-07-31 | 上海复泉工程技术有限公司 | 产量检测智能化抽油机 |
CN112031748A (zh) * | 2020-09-14 | 2020-12-04 | 南京富岛信息工程有限公司 | 一种基于示功图特征的抽油机井异常工况诊断方法 |
CN113445993A (zh) * | 2021-07-06 | 2021-09-28 | 燕山大学 | 一种基于离散***的抽油机井泵功图转化模型 |
CN113898334A (zh) * | 2021-10-14 | 2022-01-07 | 辽宁弘毅科技有限公司 | 一种抽油机井多功能综合测试仪参数智能分析方法及*** |
CN114263458A (zh) * | 2021-12-27 | 2022-04-01 | 西安健尚智能科技有限公司 | 一种油井工况全感知智能诊断自动处理的方法和*** |
CN114444620A (zh) * | 2022-04-08 | 2022-05-06 | 中国石油大学(华东) | 一种基于生成式对抗神经网络的示功图故障诊断方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6748341B2 (en) * | 2002-04-12 | 2004-06-08 | George E. Crowder, Jr. | Method and device for machinery diagnostics and prognostics |
CN101285463A (zh) * | 2008-05-30 | 2008-10-15 | 华中科技大学 | 一种有杆泵抽油机抽油控制方法及其*** |
CN101660402A (zh) * | 2009-09-11 | 2010-03-03 | 河海大学 | 基于物理意义的抽油井示功图凡尔开闭点提取方法 |
CN101819436A (zh) * | 2010-04-20 | 2010-09-01 | 山东卓尔电气有限公司 | 一种油田抽油机运行的监控*** |
-
2010
- 2010-11-24 CN CN 201010557234 patent/CN102011576B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6748341B2 (en) * | 2002-04-12 | 2004-06-08 | George E. Crowder, Jr. | Method and device for machinery diagnostics and prognostics |
CN101285463A (zh) * | 2008-05-30 | 2008-10-15 | 华中科技大学 | 一种有杆泵抽油机抽油控制方法及其*** |
CN101660402A (zh) * | 2009-09-11 | 2010-03-03 | 河海大学 | 基于物理意义的抽油井示功图凡尔开闭点提取方法 |
CN101819436A (zh) * | 2010-04-20 | 2010-09-01 | 山东卓尔电气有限公司 | 一种油田抽油机运行的监控*** |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184414B (zh) * | 2011-05-16 | 2013-03-13 | 中国石油天然气股份有限公司 | 泵示功图的识别和判断方法及其*** |
CN102184414A (zh) * | 2011-05-16 | 2011-09-14 | 中国石油天然气股份有限公司 | 泵示功图的识别和判断方法及其*** |
CN103806893A (zh) * | 2012-11-14 | 2014-05-21 | 中国石油天然气股份有限公司 | 智能检测方法、装置及检测*** |
CN103806893B (zh) * | 2012-11-14 | 2016-07-06 | 中国石油天然气股份有限公司 | 智能检测方法、装置及检测*** |
CN103541723A (zh) * | 2013-11-12 | 2014-01-29 | 丁涛 | 基于地面示功图面积变化的抽油机井实时工况诊断方法 |
CN103541723B (zh) * | 2013-11-12 | 2016-03-30 | 丁涛 | 基于地面示功图面积变化的抽油机井实时工况诊断方法 |
CN104948436A (zh) * | 2014-03-25 | 2015-09-30 | 布里斯托尔D/B/A远程自动化解决方案公司 | 用于确定井下泵的产量的方法和装置 |
CN104948436B (zh) * | 2014-03-25 | 2018-11-13 | 布里斯托尔D/B/A远程自动化解决方案公司 | 用于确定井下泵的产量的方法和装置 |
CN104110251A (zh) * | 2014-06-24 | 2014-10-22 | 安徽多杰电气有限公司 | 一种基于art2的抽油机示功图识别方法 |
CN105257277B (zh) * | 2015-05-15 | 2018-04-06 | 渤海大学 | 基于多变量灰色模型的有杆泵抽油井井下故障预测方法 |
CN105257277A (zh) * | 2015-05-15 | 2016-01-20 | 渤海大学 | 基于多变量灰色模型的有杆泵抽油井井下故障预测方法 |
CN105672988B (zh) * | 2015-12-30 | 2019-09-06 | 中国石油天然气股份有限公司 | 抽油机示功图诊断***和诊断方法 |
CN105672988A (zh) * | 2015-12-30 | 2016-06-15 | 中国石油天然气股份有限公司 | 抽油机示功图诊断***和诊断方法 |
CN106089184A (zh) * | 2016-07-18 | 2016-11-09 | 中国石油天然气股份有限公司 | 一种井下抽油泵工况的诊断方法和装置 |
CN106089184B (zh) * | 2016-07-18 | 2019-12-10 | 中国石油天然气股份有限公司 | 一种井下抽油泵工况的诊断方法和装置 |
CN106121622A (zh) * | 2016-07-27 | 2016-11-16 | 渤海大学 | 一种基于示功图的有杆泵抽油井的多故障诊断方法 |
CN106121622B (zh) * | 2016-07-27 | 2019-10-25 | 渤海大学 | 一种基于示功图的有杆泵抽油井的多故障诊断方法 |
CN106250692A (zh) * | 2016-07-29 | 2016-12-21 | 广东工业大学 | 一种心率数据质量评估方法 |
CN106682815A (zh) * | 2016-11-30 | 2017-05-17 | 西安石油大学 | 一种抽油井监控***及监控方法 |
WO2020019681A1 (zh) * | 2018-07-25 | 2020-01-30 | 北京国双科技有限公司 | 采油设备故障诊断方法及装置 |
CN110821478A (zh) * | 2018-08-13 | 2020-02-21 | 中国石油天然气股份有限公司 | 油井抽油泵漏失检测方法及装置 |
CN110821478B (zh) * | 2018-08-13 | 2022-07-05 | 中国石油天然气股份有限公司 | 油井抽油泵漏失检测方法及装置 |
CN111472726B (zh) * | 2020-03-27 | 2022-03-25 | 上海复泉工程技术有限公司 | 产量检测智能化抽油机 |
CN111472726A (zh) * | 2020-03-27 | 2020-07-31 | 上海复泉工程技术有限公司 | 产量检测智能化抽油机 |
CN112031748A (zh) * | 2020-09-14 | 2020-12-04 | 南京富岛信息工程有限公司 | 一种基于示功图特征的抽油机井异常工况诊断方法 |
CN112031748B (zh) * | 2020-09-14 | 2023-09-01 | 南京富岛信息工程有限公司 | 一种基于示功图特征的抽油机井异常工况诊断方法 |
CN113445993A (zh) * | 2021-07-06 | 2021-09-28 | 燕山大学 | 一种基于离散***的抽油机井泵功图转化模型 |
CN113898334A (zh) * | 2021-10-14 | 2022-01-07 | 辽宁弘毅科技有限公司 | 一种抽油机井多功能综合测试仪参数智能分析方法及*** |
CN113898334B (zh) * | 2021-10-14 | 2023-09-05 | 辽宁弘毅科技有限公司 | 一种抽油机井多功能综合测试仪参数智能分析方法及*** |
CN114263458A (zh) * | 2021-12-27 | 2022-04-01 | 西安健尚智能科技有限公司 | 一种油井工况全感知智能诊断自动处理的方法和*** |
CN114263458B (zh) * | 2021-12-27 | 2023-04-07 | 西安健尚智能科技有限公司 | 一种油井工况全感知智能诊断自动处理的方法和*** |
CN114444620A (zh) * | 2022-04-08 | 2022-05-06 | 中国石油大学(华东) | 一种基于生成式对抗神经网络的示功图故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102011576B (zh) | 2013-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102011576B (zh) | 基于示功图的有杆抽油***故障递阶诊断方法 | |
CN104712542B (zh) | 一种基于物联网的往复压缩机敏感特征提取与故障诊断方法 | |
CN101251564B (zh) | 一种利用可拓学与粗糙集理论相结合的电力变压器故障诊断方法 | |
CN201943686U (zh) | 油井故障实时预警*** | |
CN103754718B (zh) | 电梯运行安全监测***及其方法 | |
CN102043900B (zh) | 基于示功图的有杆抽油***故障预测方法 | |
CN101196185B (zh) | 自动变速箱油泵总成的试验方法 | |
CN102094626A (zh) | 油井故障实时预警方法和*** | |
CN111441864A (zh) | 发动机健康诊断方法及发动机诊断*** | |
CN108022058A (zh) | 一种风力机状态可靠性评估方法 | |
CN102445941B (zh) | 列控***车载设备互联互通测试结果自动判定及分析方法 | |
CN105930629A (zh) | 一种基于海量运行数据的在线故障诊断方法 | |
CN109886433A (zh) | 智能识别城市燃气管道缺陷的方法 | |
WO2015101856A2 (en) | A method and system for fault detection during a plunger lift operation in hydrocarbon wells | |
CN103389701A (zh) | 基于分布式数据模型的厂级过程故障检测与诊断方法 | |
CN105388884A (zh) | 一种基于由数据驱动的识别算法检测供热管网泄漏故障的报警***及方法 | |
CN110486350B (zh) | 电液伺服阀故障诊断方法及装置、存储介质和电子设备 | |
CN111144433A (zh) | 基于svm模型油井工况智能诊断分析方法及装置 | |
CN116561670B (zh) | 一种金属屋面健康状态识别及报警方法 | |
CN105784340A (zh) | 基于混合智能技术的气阀故障诊断方法 | |
CN105760672B (zh) | 一种机械设备故障诊断方法 | |
CN109271314A (zh) | 一种软件测试人员绩效评估***及其使用方法 | |
CN107784148B (zh) | 一种集输管道失效率评估方法及装置 | |
CN108346111A (zh) | 集输管网泄漏风险评估方法及装置 | |
CN110805549A (zh) | 一种注聚泵故障诊断*** |
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: 20130925 Termination date: 20161124 |