基于曲波变换和核稀疏的抽油井半监督故障诊断方法技术

本发明专利技术涉及一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法，步骤为：获取多个示功图数据作为训练样本；将多个示功图转化为井下泵功图，再将每个泵功图转化为灰度图像；对每个泵功图进行曲波变换得到系数矩阵；将全部有标签的泵功图的特征向量作为字典，对每个未标记泵功图特征向量求稀疏系数；利用稀疏系数计算每个无标签泵功图的虚拟标签；将训练样本中所有泵功图的特征向量作为字典；对每一个待诊断的测试样本计算其特征向量求得稀疏系数；利用稀疏系数计算待诊断样本的虚拟标签，断故障类型。本发明专利技术能精确的描述出泵功图的特征，基于核方法的半监督稀疏表达分类器不仅可以有效的利用未标记数据的信息，而且对有标记的样本数量要求不高。

Semi supervised fault diagnosis method for pumping well based on Qu Bo transform and kernel sparse

The invention relates to a semi supervised method for fault diagnosis of oil pumping, Qu Bo transform and kernel based on sparse steps: obtaining a plurality of dynamometer data as the training sample; a plurality of dynagrams into downhole pump dynamograph, then each pump dynamograph into gray image; for each pump dynamograph of Qu Bo transform coefficient matrix; all the feature vectors of pump dynamograph label as a dictionary, for each unlabeled pump dynamograph feature vector for sparse coefficient; virtual label label free pump dynamograph is calculated for each use of sparse coefficient vectors; all pump dynamograph in the training sample as a dictionary for each; a diagnostic test sample to calculate the feature vector to obtain sparse coefficient; virtual label to calculate diagnostic samples using sparse coefficient, fault type. The invention can accurately describe the characteristics of pump dynamograph, kernel method for semi supervised sparse representation classifier can not only effectively using unlabeled data based on the information, but not high requirements on the number of labeled samples.

【技术实现步骤摘要】
基于曲波变换和核稀疏的抽油井半监督故障诊断方法
本专利技术涉及一种抽油井故障诊断技术，具体地说是一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法。
技术介绍
在实际的石油开采过程中，常用的有杆泵抽油井系统所处的生产环境较为恶劣，故障发生率高，经常无法正常工作。各种故障会导致诸如产油量减少，油井停产甚至设备损坏等严重后果。在传统的石油生产中，一般依靠有经验的技术人员人为的对采集到的井上示功图或者经转化而得到的井下泵功图进行分析，来判断系统当前是否正常工作。这种方法效率低，且受主观影响大，无法满足企业的生产需求。因此，利用计算机自动实现有杆泵抽油井系统的故障诊断是一种非常有意义的技术。有杆泵抽油井系统的故障诊断可以看成是一个模式识别问题，通常分为两个步骤：特征提取和分类。常用的针对井上示功图或者井下泵功图的特征提取方法主要有网格法，链码法，傅里叶级数法等；分类一般采用经典的支持向量机和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；5)将全部有标签的泵功图的特征向量作为字典D，D＝{Vi}，i＝1，…，l；对每个未标记泵功图特征向量Vj,j＝l,…,l+u，求稀疏系数6)利用稀疏系数，计算每个无标签泵功图的虚拟标签j＝l,…,l+u；7)将训练样本中所有泵功图的特征向量作为字典D'，D'＝{Vi}，i＝1,…,l+u；对每一个待诊断的测试样本，计算其特征向量Vtest，并求得其对应的稀疏系数8)利用稀疏系数计算待诊断样本的虚拟标签ytest，并据此判断故障类型。步骤4)中，根据第一个尺度的系数矩阵ci1计算单幅泵功图的特征向量Vi为：对由曲波变换得到的尺度系数矩阵中的第一个尺度的系数矩阵ci1，采用非重叠的方式按照2×2的大小进行二进剖分，分别计算每个剖分子块的灰度均值，得到一个特征矩阵c′i1，将矩阵c′i1按照如下公式进行归一化：将归一化后的矩阵c′i1的每一列级联成一个向量，即得到单幅泵功图的特征向量Vi。步骤6)利用稀疏系数计算每个无标签泵功图的虚拟标签j＝l,…,l+u为：利用稀疏系数和全部有标签的泵功图的标签集合YD＝{yi},i＝1,...,l直接计算每个无标签泵功图的虚拟标签j＝l,…,l+u。步骤8)利用稀疏系数计算待诊断样本的标签ytest为：利用稀疏系数全部有标签泵功图的标签集合YD＝{yi},i＝1,...,l，以及全部无标签泵功图的虚拟标签集合i＝l,...,l+u直接计算待诊断样本的虚拟标签标签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；5)将全部有标签的泵功图的特征向量作为字典D，D＝{Vi}，i＝1，…，l；对每个未标记泵功图特征向量Vj，j＝l,…,l+u，求稀疏系数6)利用稀疏系数，计算每个无标签泵功图的虚拟标签j＝l,…,l+u。7)将训练样本中所有泵功图的特征向量作为字典D'，D'＝{Vi}，i＝1,…,l+u；对每一个待诊断的测试样本，其特征向量Vtest，求得其稀疏系数8)利用稀疏系数计算待诊断样本的标签ytest。步骤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杆柱末端(泵)处，载荷为：根据所求得的泵位移和本文档来自技高网...

【技术保护点】
一种基于曲波变换和核稀疏的抽油井半监督故障诊断方法，其特征在于包括以下步骤：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；5)将全部有标签的泵功图的特征向量作为字典D，D＝{Vi}，i＝1，…，l；对每个未标记泵功图特征向量Vj,j＝l,…,l+u，求稀疏系数

【技术特征摘要】
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；5)将全部有标签的泵功图的特征向量作为字典D，D＝{Vi}，i＝1，…，l；对每个未标记泵功图特征向量Vj,j＝l,…,l+u，求稀疏系数6)利用稀疏系数，计算每个无标签泵功图的虚拟标签j＝l,…,l+u；7)将训练样本中所有泵功图的特征向量作为字典D'，D'＝{Vi}，i＝1,…,l+u；对每一个待诊断的测试样本，计算其特征向量Vtest，并求得其对应的稀疏系数8)利用稀疏系数计算待诊断样本的虚拟标签ytest，并据此判断故障类型。2.按权利要求1所述的基于曲波变换和核稀疏的抽油井半监督故障诊断方法，其特征在于：步骤4)中，根据第一个尺度的系数矩阵ci1计算单...

【专利技术属性】
技术研发人员：高宪文，王明顺，张遨，张平，魏晶亮，郑博元，陈星宇，宋乐，
申请(专利权)人：东北大学，
类型：发明
国别省市：辽宁,21