CN107165615B - 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 - Google Patents
基于曲波变换和核稀疏的抽油井半监督故障诊断方法 Download PDFInfo
- Publication number
- CN107165615B CN107165615B CN201710326671.0A CN201710326671A CN107165615B CN 107165615 B CN107165615 B CN 107165615B CN 201710326671 A CN201710326671 A CN 201710326671A CN 107165615 B CN107165615 B CN 107165615B
- Authority
- CN
- China
- Prior art keywords
- pump
- diagram
- label
- diagrams
- calculating
- 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.)
- Active
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2136—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on sparsity criteria, e.g. with an overcomplete basis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geology (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Mining & Mineral Resources (AREA)
- Bioinformatics & Computational Biology (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Complex Calculations (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明涉及一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法,步骤为:获取多个示功图数据作为训练样本;将多个示功图转化为井下泵功图,再将每个泵功图转化为灰度图像;对每个泵功图进行曲波变换得到系数矩阵;将全部有标签的泵功图的特征向量作为字典,对每个未标记泵功图特征向量求稀疏系数;利用稀疏系数计算每个无标签泵功图的虚拟标签;将训练样本中所有泵功图的特征向量作为字典;对每一个待诊断的测试样本计算其特征向量求得稀疏系数;利用稀疏系数计算待诊断样本的虚拟标签,断故障类型。本发明能精确的描述出泵功图的特征,基于核方法的半监督稀疏表达分类器不仅可以有效的利用未标记数据的信息,而且对有标记的样本数量要求不高。
Description
技术领域
本发明涉及一种抽油井故障诊断技术,具体地说是一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法。
背景技术
在实际的石油开采过程中,常用的有杆泵抽油井系统所处的生产环境较为恶劣,故障发生率高,经常无法正常工作。各种故障会导致诸如产油量减少,油井停产甚至设备损坏等严重后果。在传统的石油生产中,一般依靠有经验的技术人员人为的对采集到的井上示功图或者经转化而得到的井下泵功图进行分析,来判断系统当前是否正常工作。这种方法效率低,且受主观影响大,无法满足企业的生产需求。因此,利用计算机自动实现有杆泵抽油井系统的故障诊断是一种非常有意义的技术。
有杆泵抽油井系统的故障诊断可以看成是一个模式识别问题,通常分为两个步骤:特征提取和分类。常用的针对井上示功图或者井下泵功图的特征提取方法主要有网格法,链码法,傅里叶级数法等;分类一般采用经典的支持向量机和BP神经网络。这些经典方法往往并不能很好的结合实际生产情况。
当前的有杆泵抽油井系统的故障诊断方法大都为有监督方法,即所有的训练数据都人为的标记出故障类型。在实际生产中,有标签的数据是非常昂贵的,很难获得大量的标记数据,相对的,存在大量的未被利用的无标记的数据。如何利用好这些无标记数据的信息,是一个亟待解决的问题。
发明内容
针对现有技术中有杆泵抽油井系统故障诊断方法中存在大量的未被利用的无标记的数据、不能很好的结合实际生产情况等不足,本发明要解决的问题是提供一种可有效利用未标记数据信息的基于曲波变换和核稀疏的抽油井半监督故障诊断方法。
为解决上述技术问题,本发明采用的技术方案是:
本发明一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法,包括以下步骤:
1)通过现场示功仪获取n(n=l+u)个示功图数据作为训练样本,其中l个示功图为已知标签数据,u个示功图为无标签数据;
2)根据波动方程,利用有限差分法将n个示功图转化为井下泵功图,再将每个泵功图转化大小为m×m像素的灰度图像;
3)对每个泵功图Xi进行曲波变换,得到第i个泵功图的s个尺度的系数矩阵Ci:
Ci={cij},i=1,…,n,j=1,…,s,其中n为泵功图总数量,s=log2m-3;
4)根据第一个尺度的系数矩阵ci1计算单幅泵功图的特征向量Vi;
步骤4)中,根据第一个尺度的系数矩阵ci1计算单幅泵功图的特征向量Vi为:对由曲波变换得到的尺度系数矩阵中的第一个尺度的系数矩阵ci1,采用非重叠的方式按照2×2的大小进行二进剖分,分别计算每个剖分子块的灰度均值,得到一个特征矩阵c′i1,将矩阵c′i1按照如下公式进行归一化:
将归一化后的矩阵c′i1的每一列级联成一个向量,即得到单幅泵功图的特征向量Vi。
标签ytest中数值最大的元素所对应的序号i,即代表当前发生第i类故障。
本发明具有以下有益效果及优点:
1.本发明通过对井下泵功图进行曲波变换后得到的特征向量能更精确的描述出泵功图的特征,基于核方法的半监督稀疏表达分类器不仅可以有效的利用未标记数据的信息,而且对有标记的样本数量要求不高,可以节省大量用来标记样本的人力成本。
2.本发明通过有限差分法将井上示功图转化为井下泵功图,并对得到的井下泵功图进行曲波变换并提取特征,所得到的特征向量能更精确的描述出泵功图的特征,以便于后续诊断。
3.基于核方法的半监督稀疏表达分类器只需要计算训练样本间的稀疏关系,不需要预先测试超参数,具有很强的泛化能力,可以更好的结合实际生产情况。
附图说明
图1为本发明出基于曲波变换和核稀疏的抽油井半监督故障诊断方法流程图;
图2为本发明方法中有限差分法求解泵功图中补格法示意图;
图3为本发明方法中利用有限差分法求得的泵功图;
图4为本发明方法对泵功图进行曲波变换得到的系数矩阵可视化的结果示意图。
具体实施方式
下面结合说明书附图对本发明作进一步阐述。
如图1所示,本发明一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法包括以下步骤:
1)通过现场示功仪获取n(n=l+u)个示功图数据作为训练样本,其中l个示功图为已知标签数据,u个示功图为无标签数据;
2)根据经典波动方程,利用有限差分法将n个示功图转化为井下泵功图,再将每个泵功图转化大小为256×256像素的灰度图像;
3)对每个泵功图Xi进行曲波变换,得到第i个泵功图的s个尺度的系数矩阵Ci:
Ci={cij},i=1,…,n,j=1,…,s,其中n为泵功图总数量,s=log2256-3=5;
4)据第一个尺度的系数矩阵ci1计算单幅泵功图的特征向量Vi;
步骤2)中,以某油田为例,三级组合杆柱每根杆长为523.61m,664.32m和618.35m,杆柱的直径为0.025m,0.022m和0.019m。抽油杆密度为8456kg/m3,石油密度为843kg/m3,弹性模量E=2.1×1011Pa,粘滞阻尼系数c=0.4654,应力传播速度a=4970m/s。
根据如下波动方程:
利用有限差分法将所有井上示功图转化为井下泵功图,设驴头的下始点为坐标原点,将每根抽油杆分为k段,选取每段的步长为Δx,同理,选取Δt为时间t的步长,下角标i表示位置,j表示时间,则有:
将以上三式带入波动方程可以得到:
根据井上示功图可知地面光杆位移为u1,u2,…,uk,光杆动载荷为F1,F2,...,Fk,根据胡克定律则波动方程的边界条件为:
(1)u1,1=-u1;u1,2=-u2;…;u1,k=-uk
其中A为抽油杆截面积。
由于示功图是一个周期函数,如图2所示,采用补格法可得到边界条件:
(3)ui,k=ui,0
(4)ui+1,1=ui,k+1
步长Δx和Δt的选取,应满足如下的稳定性条件:
在杆柱末端(泵)处,位移为:
up=-ui+1
杆柱末端(泵)处,载荷为:
根据所求得的泵位移和泵载荷,即可绘制出井下泵功图。所求得的泵功图如图3所示。之后将每个泵功图转化大小为256×256像素的灰度图像。步骤3)中,利用第二代Curvelet工具包对每个泵功图Xi进行曲波变换,得到5个尺度的系数矩阵Ci,Ci={cij},i=1,…,n,j=1,2,3,4,5。泵功图及其对应的系数矩阵可视化结果如图4所示。
步骤4)中,将第一个尺度的系数矩阵ci1,采用非重叠的方式按照2×2的大小进行二进剖分,分别计算每个剖分子块的灰度均值,得到一个10×10的矩阵c′i1,将矩阵c′i1按照如下公式进行归一化:
将归一化后的矩阵c′i1的每一列级联成一个100×1的向量,得到单幅泵功图的特征向量Vi。
k(x,xc)=exp(-||x-xc||/σ)
其中,xc为核函数中心和σ为核函数的宽度参数,控制了核函数的径向控制范围。取σ=8d,d为特征向量的维度。
其中,φ(·)代表非线性映射函数,并满足φ(x1)φ(x2)=k(x1,x2),核函数向量kj=[k(Vj,V1),...,k(Vj,Vl)],K为如下形式的核函数矩阵:
其中,YD为全部有标签的泵功图的标签向量集合,YD={yi},i=1,...,l。
其中,核向量ktest=[k(Vtest,V1),...,k(Vtest,Vl)],K为如下形式的核函数矩阵:为
虚拟标签ytest中数值最大的元素所对应的序号i,即代表当前发生第i类故障。
本发明通过有限差分法将所有训练数据即井上示功图转化为井下泵功图并转为灰度图像,利用曲波变换对泵功图进行特征提取后,将其中所有的有标签的泵功图的特征向量作为字典,对每一个无标签泵功图特征向量求其相对应的稀疏表达系数,根据稀疏表达系数计算虚拟标签以补所有全训练数据的标签,最后将训练样本中所有泵功图的特征向量作为字典,对每一个待诊断的测试样本,求其特征向量所对应的稀疏系数,并计算诊断样本的标签,完成诊断。本发明所采用的曲波变换,能更精确的描述出泵功图的特征;基于核方法的稀疏表达分类器不仅可以有效的利用未标记数据以便提高故障诊断的正确率,而且对有标记的样本数量要求不高。
Claims (4)
1.一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法,其特征在于包括以下步骤:
1)通过现场示功仪获取n(n=l+u)个示功图数据作为训练样本,其中l个示功图为已知标签数据,u个示功图为无标签数据;
2)根据波动方程,利用有限差分法将n个示功图转化为井下泵功图,再将每个泵功图转化大小为m×m像素的灰度图像;
3)对每个泵功图Xi进行曲波变换,得到第i个泵功图的s个尺度的系数矩阵Ci:
Ci={cij},i=1,…,n,j=1,…,s,其中n为泵功图总数量,s=log2m-3;
4)根据第一个尺度的系数矩阵ci1计算单幅泵功图的特征向量Vi;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710326671.0A CN107165615B (zh) | 2017-05-10 | 2017-05-10 | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710326671.0A CN107165615B (zh) | 2017-05-10 | 2017-05-10 | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107165615A CN107165615A (zh) | 2017-09-15 |
CN107165615B true CN107165615B (zh) | 2020-04-24 |
Family
ID=59813916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710326671.0A Active CN107165615B (zh) | 2017-05-10 | 2017-05-10 | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107165615B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764361B (zh) * | 2018-06-01 | 2020-12-22 | 北京中油瑞飞信息技术有限责任公司 | 基于集成学习的游梁式抽油机示功图的工况识别方法 |
CN108798641B (zh) * | 2018-06-19 | 2021-06-11 | 东北大学 | 一种基于子空间迁移学习的有杆泵抽油井故障诊断方法 |
CN108979624B (zh) * | 2018-08-07 | 2022-03-08 | 东北大学 | 一种基于示功图矩特征的有杆抽油系统摩擦因数辨识方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006075337A2 (en) * | 2005-01-12 | 2006-07-20 | Oil And Natural Gas Corporation Limited | Apparatus and method for monitoring remotely located sucker rod pumps |
CN100404786C (zh) * | 2005-09-23 | 2008-07-23 | 中国石油天然气股份有限公司 | 功图法油井计量装置 |
US7823640B2 (en) * | 2007-10-23 | 2010-11-02 | Saudi Arabian Oil Company | Wellhead flowline protection and testing system with ESP speed controller and emergency isolation valve |
CN100572810C (zh) * | 2008-05-30 | 2009-12-23 | 华中科技大学 | 一种有杆泵抽油机抽油控制方法及其系统 |
CN201943686U (zh) * | 2010-12-20 | 2011-08-24 | 中国石油天然气股份有限公司 | 油井故障实时预警系统 |
CN102129573A (zh) * | 2011-03-10 | 2011-07-20 | 西安电子科技大学 | 基于字典学习和稀疏表示的sar图像分割方法 |
CN102324047B (zh) * | 2011-09-05 | 2013-06-12 | 西安电子科技大学 | 基于稀疏核编码skr的高光谱图像地物识别方法 |
CN104570991A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油天然气股份有限公司 | 基于功图的抽油机调控方法 |
CN103793694B (zh) * | 2014-02-10 | 2017-02-08 | 天津大学 | 一种基于多特征空间稀疏分类器的人脸识别方法 |
CN103899300A (zh) * | 2014-03-25 | 2014-07-02 | 中国石油天然气股份有限公司 | 一种基于示功图的二流量试井分析的方法及系统 |
CN104481496B (zh) * | 2014-09-16 | 2017-02-15 | 渤海大学 | 一种有杆泵抽油井故障诊断方法 |
CN105257277B (zh) * | 2015-05-15 | 2018-04-06 | 渤海大学 | 基于多变量灰色模型的有杆泵抽油井井下故障预测方法 |
CN106321072B (zh) * | 2015-06-15 | 2019-02-19 | 中国科学院沈阳自动化研究所 | 一种基于泵功图的抽油井故障诊断方法 |
CN105089638A (zh) * | 2015-06-26 | 2015-11-25 | 中国石油化工股份有限公司胜利油田分公司 | 一种利用抽油机泵功图在线计算油井动液面方法 |
CN106089184B (zh) * | 2016-07-18 | 2019-12-10 | 中国石油天然气股份有限公司 | 一种井下抽油泵工况的诊断方法和装置 |
CN106121622B (zh) * | 2016-07-27 | 2019-10-25 | 渤海大学 | 一种基于示功图的有杆泵抽油井的多故障诊断方法 |
-
2017
- 2017-05-10 CN CN201710326671.0A patent/CN107165615B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107165615A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107165615B (zh) | 基于曲波变换和核稀疏的抽油井半监督故障诊断方法 | |
WO2020063689A1 (zh) | 一种电力系统热负荷预测的方法和装置 | |
CN104481496B (zh) | 一种有杆泵抽油井故障诊断方法 | |
Zheng et al. | Sucker rod pumping diagnosis using valve working position and parameter optimal continuous hidden Markov model | |
CN106930751A (zh) | 一种有杆泵抽油井故障分离方法 | |
CN104361365A (zh) | 一种抽油泵运行状态识别方法及装置 | |
CN112611584A (zh) | 风力发电机组的疲劳失效检测方法、装置、设备及介质 | |
CN113095414A (zh) | 一种基于卷积神经网络和支持向量机的示功图识别方法 | |
CN114548190A (zh) | 一种基于自适应残差神经网络的风力涡轮机故障诊断方法 | |
CN115267945A (zh) | 一种基于图神经网络的雷电预警方法与系统 | |
Xin et al. | An intelligent fault diagnosis method of rotating machinery based on deep neural networks and time-frequency analysis | |
CN114021620B (zh) | 基于bp神经网络特征提取的电潜泵故障诊断方法 | |
CN117056402B (zh) | 一种基于多源信号的电机诊断方法、装置及存储介质 | |
Liu et al. | Research on image segmentation algorithm and performance of power insulator based on adaptive region growing | |
Liu et al. | Fault diagnosis of rod pump oil well based on support vector machine using preprocessed indicator diagram | |
CN112802011A (zh) | 一种基于vgg-bls的风机叶片缺陷检测方法 | |
Yang et al. | A fault identification method for electric submersible pumps based on dae-svm | |
Liu et al. | Semi-supervised deep learning recognition method for the new classes of faults in wind turbine system | |
CN108459585B (zh) | 基于稀疏局部嵌入深度卷积网络的电站风机故障诊断方法 | |
CN111523557A (zh) | 基于大数据的风电智能化故障诊断方法 | |
CN114120043A (zh) | 一种基于生产动态数据和示功图检测异常抽油井的方法 | |
Yang et al. | Convolutional autoencoder-based sensor fault classification | |
CN110348532B (zh) | 基于机器学习与残差网络构成的有杆泵工况组合监测方法 | |
CN112664185A (zh) | 一种基于示功图的抽油机井工况预测方法 | |
CN116698410B (zh) | 一种基于卷积神经网络的滚动轴承多传感器数据监控方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |