CN113486599B - 一种抽油机有效冲程的计算方法 - Google Patents
一种抽油机有效冲程的计算方法 Download PDFInfo
- Publication number
- CN113486599B CN113486599B CN202110927362.5A CN202110927362A CN113486599B CN 113486599 B CN113486599 B CN 113486599B CN 202110927362 A CN202110927362 A CN 202110927362A CN 113486599 B CN113486599 B CN 113486599B
- Authority
- CN
- China
- Prior art keywords
- indicator diagram
- data
- regression model
- map
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005086 pumping Methods 0.000 title claims abstract description 45
- 238000010586 diagram Methods 0.000 claims abstract description 116
- 238000012360 testing method Methods 0.000 claims abstract description 34
- 238000012549 training Methods 0.000 claims abstract description 21
- 239000003129 oil well Substances 0.000 claims abstract description 20
- 230000002159 abnormal effect Effects 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 238000005065 mining Methods 0.000 claims abstract description 4
- 238000013528 artificial neural network Methods 0.000 claims description 38
- 238000006073 displacement reaction Methods 0.000 claims description 29
- 239000000725 suspension Substances 0.000 claims description 27
- 238000013507 mapping Methods 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 16
- 238000012795 verification Methods 0.000 claims description 13
- 238000010606 normalization Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 5
- 238000011176 pooling Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims 1
- 239000007787 solid Substances 0.000 abstract 1
- 239000003921 oil Substances 0.000 description 24
- 238000004364 calculation method Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 238000012935 Averaging Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000013135 deep learning Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Medical Informatics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种抽油机有效冲程的计算方法,首先获取油井现场采集到的抽油机示功图历史数据,并对其进行异常数据过滤;然后将有效工况数据归一化并绘制为闭合的示功图图像;接着为示功图数据集标定与有效冲程相关的四个凡尔开闭点;再以带有标签的数据集训练凡尔开闭点回归模型,并根据测试结果挖掘困难样本,进而调整网络结构以提高回归模型性能;最后将现场的实时数据输入到训练完成的点回归模型中得到预测的凡尔开闭点,从而推算出实际的抽油机有效冲程。该方法能够准确计算抽油机的有效冲程,从而为功图量油提供坚实的技术支撑。
Description
技术领域
本发明涉及基于油田量油技术领域,具体而言,是通过深度学习技术,建立数据驱动的回归模型,从而预测抽油机井有效冲程的方法。
背景技术
在油田开采过程中,油井产量计量是油田日常生产管理的一项基础且重要的工作,目前国内外传统的计量方式是将各井口的产油输送到计量站进行集中计量,这一计量方式存在安装设备繁多、工艺流程冗杂、实时性低等问题,难以满足智慧油田的发展需求。在这种需求下,功图量油成为替代目前计量站形式量油的主要手段。
功图量油是通过对油井示功图的实时采集和处理分析来计算油井产液量的,而柱塞有效冲程的计算直接影响功图量油的准确性。柱塞有效冲程由抽油泵固定凡尔和游动凡尔的开启、关闭四个关键点所决定。四个凡尔开闭点的求解,目前工业中普遍采用五点平均法,该方法通过计算功图中每个采样点的平均曲率变化值,并分别将功图上、下冲程中曲率变化最大的两个采样点作为凡尔开闭点。五点平均法在工程应用中具有计算简单和实时性好的优点,但同时也存在着对噪声干扰敏感和对功图完整性要求高的缺点。在油田实际开采过程中,噪声影响和极端故障导致的功图不规整、毛刺多和曲率变化波动大的情况时有发生,此时五点平均法会出现计算误差大,可靠性降低,甚至难以使用的情况。
近年来,随着数据驱动技术的快速发展,将深度学习等人工智能技术应用于油田生产开采过程,有望为解决上述问题提供另一种新的思路,为准确定位凡尔开闭点从而计算有效冲程和油井产量提供技术支撑。
发明内容
针对上述问题,本发明提供一种抽油机有效冲程的计算方法,通过示功图的历史采样数据训练卷积神经网络回归模型,对抽油机的示功图凡尔开闭点进行实时检测,从而预测示功图的有效冲程,具有计算流程简单、实时性高和泛化能力强等特点。
本发明提出一种抽油机有效冲程的计算方法,具有如下步骤:
(1)获取油井现场采集到的抽油机示功图历史数据并对其进行异常数据过滤处理,得到有效工况数据集D:
D={di|i=1,2,...,N}
di={(xij,yij)|j=1,2,...,m}
其中di表示第i组抽油机示功仪采样数据,共N组;其中每组数据包含m对抽油机悬点位移和载荷采样值(xij,yij),xij表示第j个悬点位移采样值,yij表示第j个悬点载荷采样值;
(2)将有效工况数据归一化并绘制为闭合的示功图图像,得到有效示功图数据集P:
P={pi|i=1,2,...,K}
其中pi表示由第i组抽油机示功仪采样数据di绘制成的示功图样本,共K个;
(3)由工艺专家为示功图数据集P标定与有效冲程相关的四个凡尔开闭点作为数据集P的标签L:
L={li|i=1,2,...,K}
li={aiA,biA,aiB,biB,aiC,biC,aiD,biD}
其中li表示第i组抽油机示功图样本pi的标签,共K组;其中每组标签包含4个凡尔开闭点在示功图中的像素坐标值,分别为固定凡尔开启点A的坐标(aiA,biA)、固定凡尔关闭点B的坐标(aiB,biB)、游动凡尔开启点C的坐标(aiC,biC)和游动凡尔关闭点D的坐标(aiD,biD);
(4)将示功图数据集P按照留出法划分为训练集Ptrain、验证集Pval和测试集Ptest,对应的标签L也被划分为训练集标签Ltrain、验证集Lval和测试集Ltest;
(5)根据凡尔开闭点识别任务的特性,选取神经网络回归模型M;
(6)利用训练集Ptrain,结合对应的训练集标签Ltrain,对神经网络回归模型M进行示功图到凡尔开闭点位置的映射关系的学习;
(7)根据神经网络回归模型M训练过程中,在验证集Pval上的损失函数的最小值选取模型的权重参数;
(8)通过测试集Ptest评估(7)中已确定网络权重参数的回归模型M的泛化能力,并根据测试结果挖掘困难样本,针对困难样本调整神经网络回归模型M的网络结构以达到提高回归模型性能的目的;
(9)重复步骤(6)至步骤(8)直至回归模型M的性能达到工程应用要求,获得神经网络回归模型M;
(10)获取油井现场采集到的实时悬点位移及载荷并对其进行异常数据过滤和归一化处理,将结果绘制成示功图输入至神经网络回归模型M,根据模型M输出的凡尔开闭点像素坐标和归一化相关参数反推出实际有效冲程。
进一步地,所述获取油井现场采集到的抽油机示功图历史数据并对其进行异常数据过滤处理具有如下步骤:
(1-1)将超出示功仪位移及载荷量程或空值数据定义为异常数据,统计一组示功图据中异常数据的占比αerror:
其中Nerroe表示异常数据点个数,Nall表示采样总数;
(1-2)若αerror大于设定阈值则丢弃当前样本;
(1-3)若αerror小于设定阈值则采用如下策略对异常采样点perror(xerror,yerror)进行修正:
其中(xpre,ypre)表示(xerror,yerror)的前一个有效采样点,(xafter,yafter)表示(xerror,yerror)的后一个有效采样点。
进一步的,所述有效工况数据归一化并将其绘制为示功图图像具有如下步骤:
(2-1)设定每台抽油机设备的示功图坐标系以其生产参数的最大载荷ymax为纵轴最大量程,以最大冲程xmax为横轴最大量程,对示功图原始二维数据进行坐标归一化,将原始数据转换到[0,1]的范围内,归一化结果如下:
dnorm_i={(xnorm_ij,ynorm_ij)|j=1,2,...,m}
其中:
(2-2)根据归一化结果映射到w×h的网格中并进行向下取整,映射结果dmap_i如下:
dmap_i={(xmap_ij,ymap_ij)|j=1,2,...,m}
其中:
(xmap_ij,ymap_ij)=(xnorm_ij,ynorm_ij)×(w,h)
(2-3)根据(2-2)中的映射结果绘制示功图图像,绘制时以归一化并映射后的悬点位移xmap_ij为图像的横坐标,以归一化后的悬点载荷ymap_ij为图像的纵坐标,将(xmap_i1,ymap_i1)与其后的一个点相连,并以此类推,将数据点两两相连至最后一点(xmap_im,ymap_im),再将(xmap_im,ymap_im)与(xmap_i1,ymap_i1)相连,形成示功图封闭曲线,示功图大小为w×h。
进一步地,所述针对困难样本调整神经网络回归模型M的网络结构包括以下措施:
(8-1)将神经网络回归模型M的高层特征通过双线性插值法上采样与模型M的低层特征融合;
(8-2)对(8-1)中融合了神经网络回归模型低层特征和高层特征的最终特征图进行卷积和池化处理并展开,用于关键点回归任务,使神经网络回归模型M在学习示功图到凡尔开闭点的映射关系时,既能考虑神经网络低层的高分辨率特征即示功图的局部信息,还能考虑到神经网络高层的语义特征也即示功图的全局信息。
进一步的,所述根据回归模型输出的凡尔开闭点像素坐标和归一化相关参数反推出实际有效冲程具有如下步骤:
(10-1)将实时的示功图有效数据dr归一化为dnorm_r,其中:
dr={(xrj,yrj)|j=1,2,...,m}
dnorm_r={(xnorm_rj,ynorm_rj)|j=1,2,...,m}
两者间的转换关系如下:
(10-2)将归一化后的数据映射至w×h大小的网格内并绘制成示功图后输入神经网络回归模型,则示功图像素点(arj,brj)和归一化后的示功图悬点及位移数据(xnorm_rj,ynorm_rj)的转换关系如下:
(arj,brj)=(xnorm_rj,ynorm_rj)×(w,h)
(10-3)结合(10-1)和(10-2)可得示功图像素点(arj,brj)和油井现场示功图悬点及位移数据(xrj,yrj)的映射关系如下:
(10-4)神经网络回归模型输出的四个凡尔开闭点像素坐标如下:
lr={arA,brA,arB,brB,arC,brC,arD,brD}
结合(10-3)中建立的映射关系反推出四个凡尔开闭点实际悬点位移及载荷物理坐标如下:
lreal_r={x1,y1,x2,y2,x3,y3,x4,y4}
其中(x1,y1)为固定凡尔的开启点A的悬点及位移、(x2,y2)为固定凡尔关闭点B的悬点及位移,(x3,y3)为游动凡尔开启点C的位移及载荷,(x4,y4)为游动凡尔关闭点D的位移及集合;
(10-5)根据(10-4)中的四个凡尔开闭点实际悬点位移及载荷物理坐标可以计算出示功图的有效冲程Spe即为固定凡尔开闭点的位移差和游动凡尔开闭点的位移差中的较小值:
Spe=min{x2-x1,x3-x4}。
有益效果
本发明公开了一种抽油机有效冲程的计算方法,该方法充分利用抽油机历史示功图数据,建立基于卷积神经网络的回归模型,实时预测示功图的四个凡尔开闭点,准确计算抽油机井的有效冲程,具有计算流程简单、实时性高和泛化能力强等特点,从而为功图量油提供坚实技术支撑。
附图说明
图1为本发明的抽油机有效冲程计算方法的系统流程框图;
图2为本发明具体算例数据预处理后的示功图;
图3为本发明具体算例的示功图样本标签可视化图;
图4为本发明具体算例回归模型的损失函数训练结果图;
图5为本发明具体算例测试集困难样本示功图;
图6为本发明具体算例调整网络结构后的回归模型损失函数训练结果图;
图7为本发明具体算例调整网络结构后的回归模型在困难样本上的预测结果图;
图8为本发明具体算例油井现场若干典型工况的示功图;
图9为本发明具体算例最终回归模型对典型工况示功图的预测结果图;
图10为本发明具体算例最终回归模型对现场供液不足示功图的预测结果图。
具体实施方式
下面结合附图以及具体的算例,由具体的操作流程说明本方法在抽油机井有效冲程的实施效果。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
选取某原油开采企业某油井2020年1月1日到2020年1月25日期间8口抽油机井的工况数据,以此建立对抽油机凡尔开闭点的回归模型,从而验证本发明所提出的基于抽油机示功图的有效冲程计算方法。
本发明的总体系统流程图如图1所示,具体实施步骤如下:
(1)获取油井现场采集到的抽油机示功图历史数据并对其进行异常数据过滤处理,得到有效工况数据集D,共计6522组,每组数据包含200个采样点:
D={di|i=1,2,...,6522}
di={(xij,yij)|j=1,2,...,200}
以原始数据样本中一号油井的第一组有效数据d1为例:
d1=[(0,50.74),(0.003,52.84),(0.009,55.29),...,(0.011,46.49),(0.003,47.06),(0,48.73)]2×200
(2)将有效工况数据归一化并绘制为闭合的示功图图像,得到有效示功图数据集P:
P={pi|i=1,2,...,6522}
以原始数据样本中一号油井的第一组有效数据为例:
首先以抽油机井生产参数的最大载荷ymax=87.91kN为纵轴最大量程,以最大冲程xmax=5.981m为横轴最大量程对示功图原始二维数据进行坐标归一化,将原始数据转换到[0,1]的范围内,归一化结果如下:
dnorm_1=[(0,0.577),(0,0.601),(0.0015,0.0629),...,(0.0183,0.528),(0,0.535),(0,0.554)]2×200
然后,将归一化结果映射到256×256的网格中并进行向下取整,映射结果dmap_i如下:
dmap_1=[(0,147),(0,153),(0,161),...,(4,135),(0,136),(0,141)]2×200
再由映射结果绘制示功图图像,绘制时以归一化并映射后的悬点位移xmap_1j为图像的横坐标,以归一化后的悬点载荷ymap_ij为图像的纵坐标,采用OpenCV库绘制曲线的方法将(xmap_i1,ymap_i1)与其后的一个点相连,并以此类推,将数据点两两相连至最后一个采样点(xmap_i200,ymap_i200),再将(xmap_i200,ymap_i200)与(xmap_i1,ymap_i1)相连,形成示功图封闭曲线,示功图大小为256×256,具体示功图像如图2所示;
(3)由工艺专家为示功图数据集P标定与有效冲程相关的四个凡尔开闭点作为数据集P的标签L:
L={li|i=1,2,...,6522}
li={aiA,biA,aiB,biB,aiC,biC,aiD,biD}
其中li表示第i组抽油机示功图样本的标签,共6522组;其中每组标签包含4个凡尔开闭点在示功图中的像素坐标值,分别为固定凡尔开启点A的坐标(aiA,biA)、固定凡尔关闭点B的坐标(aiB,biB)、游动凡尔开启点C的坐标(aiC,biC)和游动凡尔关闭点D的坐标(aiD,biD),以一号油井的第一组有效数据对应的示功图为例:
li={26.52,41.56,247.09,52.85,140.93,160.47,5.78,141.54}
可视化结果如图3所示;
(4)将示功图数据集P按照留出法划分为训练集Ptrain、验证集Pval和测试集Ptest,对应的标签L也被划分为训练集标签Ltrain、验证集Lval和测试集Ltest,此处按照6:2:2的比例随机选取数据集P中的3913张样本组成训练集Ptrain,1304张样本组成验证机Pval,1305张样本组成测试集Ptest;
(5)选取ResNet50深度残差网络作为神经网络回归模型M的主干网络,将ResNet50最后一个Bottleneck的输出通过全局池化并展开,再外接全连接层输出预测的凡尔开闭点像素坐标,选取Pytorch作为深度学习网络框架。
(6)利用训练集Ptrain结合对应的训练集标签Ltrain,对回归模型M进行示功图到凡尔开闭点位置映射关系的学习,并根据验证集Pval辅助验证模型的泛化能力和调节神经网络的超参数,训练结果如图4所示;
(7)根据候选回归模型M训练过程中,在验证集Pval上的损失函数的最小值选取回归模型的权重参数;
(8)通过(4)中选取的测试集Ptest评估(7)中已确定网络权重参数的候选回归模型M的最终性能,并根据测试结果挖掘困难样本,针对困难样本调整神经网络回归模型M的网络结构以达到提高回归模型性能的目的。此处采用测试集所有样本的四个凡尔开闭点预测值与实际值的欧式距离的平均值即均方根误差作为衡量回归模型最终性能的指标,具体测试结果均方根误差为0.04435,可视化结果如图5所示。
(9)重复步骤(6)至步骤(8)直至回归模型M的性能达到工程应用要求。根据可视化的测试集预测结果,可以发现当示功图比较规整或略有缺陷时神经网络回归模型M预测结果较为准确,但当示功图存在极度不规整和频繁震荡时,神经网络回归模型M预测结果往往存在较大误差;同时,模型也存在一定的过拟合问题也即泛化能力不足。针对挖掘出的困难样本分析可知神经网络回归模型ResNet50的输出层特征具有很强的语义信息,但分辨率很低,对示功图细节的感知能力较差,在ResNet50池化层下采样的影响下,神经网络低层的高分率信息包括更多的位置、细节信息都有所丢失,因此导致最终关键点的回归缺失细节性的信息而有所偏差。此处先采用将神经网络回归模型M的主干网络ResNet50高层Bottleneck的输出特征通过1x1的卷积核和双线性插值法上采样至前一层Bottleneck的输出特征的尺寸,再将两个输出特征逐元素相叠加,然后迭代重复此操作直至ResNet50的第一个Bottleneck,最后再将叠加后的特征通过卷积层和池化层处理并展开作为最终凡尔开闭点回归的特征输入,此时叠加后的特征层既包含神经网络高层的语义信息又包含低层的轮廓及位置信息,在测试集上的均方根误差为0.04123,下降了7个百分点的同时也解决示功图极度不规整和频繁震荡时预测误差大的问题,测试可视化结果如图6所示。
(10)获取油井现场采集到的实时悬点位移及载荷并对其进行异常数据过滤和归一化处理,将结果绘制成示功图并转换成ResNet50要求的输入格式,此处选取了油井现场若干典型工况的的数据,经过以上的预处理步骤后生成的示功图如图7所示,输入至最终回归模型,回归模型的输出结果如图8所示,再根据回归模型输出的凡尔开闭点像素坐标和归一化相关参数反推出实际有效冲程,以图9示功图为例,该示功图是典型供液不足的示功图,回归模型的输出如下:
l={30,39,211,40,127,113,5,107}
结合公式以及实际生产参数如下:
xmax=5.981m
ymax=87.91kN
w=256
h=256
从而反推出实际的凡尔开闭点坐标为:
lreal={0.7,13.39,4.93,13.74,0.12,36.74,2.97,38.8}
最后根据预测的凡尔开闭的实际坐标计算有效冲程:
Spe=min{x2-x1,x3-x4}=min{4.93-0.7,2.97-0.12}=2.85m
实际现场工况的有效冲程为2.71m,误差为在可接受范围内。
由上述分析可知,本发明的抽油机有效冲程的计算方法,能够根据示功图预测抽油机的凡尔开闭点,并计算抽油机井的有效冲程。可见,通过本发明可以实现快速、高效的抽油机凡尔开闭点预测和抽油机井有效冲程计算,避免了传统五点法计算误差大,可靠性低,甚至难以使用的情况。
Claims (6)
1.一种抽油机有效冲程的计算方法,其特征在于包括以下步骤:
(1)获取油井现场采集到的抽油机示功图历史数据并对其进行异常数据过滤处理,得到有效工况数据集D:
D={di|i=1,2,...,N}
di={(xij,yij)|j=1,2,...,m}
其中di表示第i组抽油机示功仪采样数据,共N组;其中每组数据包含m对抽油机悬点位移和载荷采样值(xij,yij),xij表示第j个悬点位移采样值,yij表示第j个悬点载荷采样值;
(2)将有效工况数据归一化并绘制为闭合的示功图图像,得到有效示功图数据集P:
P={pi|i=1,2,...,N}
其中pi表示由第i组抽油机示功仪采样数据di绘制成的示功图样本,共N个;
(3)由工艺专家为示功图数据集P标定与有效冲程相关的四个凡尔开闭点作为数据集P的标签L:
L={li|i=1,2,...,N}
li={aiA,biA,aiB,biB,aiC,biC,aiD,biD}
其中li表示第i组抽油机示功图样本pi的标签,共N组;其中每组标签包含4个凡尔开闭点在示功图中的像素坐标值,分别为固定凡尔开启点A的坐标(aiA,biA)、固定凡尔关闭点B的坐标(aiB,biB)、游动凡尔开启点C的坐标(aiC,biC)和游动凡尔关闭点D的坐标(aiD,biD);
(4)将示功图数据集P划分为训练集Ptrain、验证集Pval和测试集Ptest,对应的标签L也被划分为训练集标签Ltrain、验证集标签Lval和测试集标签Ltest;
(5)根据凡尔开闭点识别任务的特性,选取神经网络回归模型M;
(6)利用训练集Ptrain,结合对应的训练集标签Ltrain,对神经网络回归模型M进行示功图到凡尔开闭点位置的映射关系学习;
(7)根据神经网络回归模型M训练过程中,在验证集Pval上的损失函数的最小值选取模型的权重参数;
(8)通过测试集Ptest评估(7)中已确定网络权重参数的回归模型M的泛化能力,并根据测试结果挖掘困难样本,针对困难样本调整神经网络回归模型M的网络结构以达到提高回归模型性能的目的;
(9)重复步骤(6)至步骤(8)直至回归模型M的性能达到工程应用要求,获得训练完成的神经网络回归模型M;
(10)获取油井现场采集到的实时悬点位移及载荷并对其进行异常数据过滤和归一化处理,将结果绘制成示功图输入至神经网络回归模型M,根据模型M输出的凡尔开闭点像素坐标和归一化相关参数反推出实际有效冲程。
2.根据权利要求1所述的一种抽油机有效冲程的计算方法,其特征在于,所述异常数据过滤处理过程具有如下步骤:
(1-1)将超出示功仪位移及载荷量程的数据或空值数据定义为异常数据,统计一组示功图据中异常数据的占比αerror:
其中Nerroe表示异常数据点个数,Nall表示采样总数;
(1-2)若αerror大于设定阈值则丢弃当前样本;
(1-3)若αerror小于设定阈值则采用如下策略对异常采样点perror(xerror,yerror)进行修正:
其中(xpre,ypre)表示(xerror,yerror)的前一个有效采样点,(xafter,yafter)表示(xerror,yerror)的后一个有效采样点。
3.根据权利要求1所述的一种抽油机有效冲程的计算方法,其特征在于,所述将有效工况数据归一化并绘制闭合的示功图图像,得到有效示功图数据集,具有如下步骤:
(2-1)设定每台抽油机设备的示功图坐标系以其生产参数的最大载荷ymax为纵轴最大量程,以最大冲程xmax为横轴最大量程,对示功图原始二维数据进行坐标归一化,将原始数据di转换到[0,1]的范围内,归一化结果dnorm_i如下:
dnorm_i={(xnorm_ij,ynorm_ij)|j=1,2,...,m}
其中:
(2-2)将(2-1)中的归一化结果映射到w×h的网格中并进行向下取整,映射结果dmap_i如下:
dmap_i={(xmap_ij,ymap_ij)|j=1,2,...,m}
其中:
(xmap_ij,ymap_ij)=(xnorm_ij,ynorm_ij)×(w,h)
(2-3)根据(2-2)中的映射结果绘制示功图图像,绘制时以归一化并映射后的悬点位移xmap_ij为图像的横坐标,以悬点载荷ymap_ij为图像的纵坐标,将(xmap_i1,ymap_i1)与其后的一个点相连,并以此类推,将数据点两两相连至最后一点(xmap_im,ymap_im),再将(xmap_im,ymap_im)与(xmap_i1,ymap_i1)相连,形成示功图封闭曲线,示功图大小为w×h。
4.根据权利要求1所述的一种抽油机有效冲程的计算方法,其特征在于,采用留出法将示功图数据集P划分为训练集Ptrain、验证集Pval和测试集Ptest,按照固定比例6:2:2将数据集静态的划分为训练集Ptrain、验证集Pval和测试集Ptest。
5.根据权利要求1所述的一种抽油机有效冲程的计算方法,其特征在于,所述针对困难样本调整神经网络回归模型M的网络结构包括以下措施:
(8-1)将神经网络回归模型M的高层特征通过双线性插值法上采样与模型M的低层特征融合;
(8-2)对(8-1)中融合了神经网络回归模型低层特征和高层特征的最终特征图进行卷积和池化处理并展开,用于关键点回归任务,使模型M在学习示功图到凡尔开闭点的映射关系时,既能考虑神经网络低层的高分辨率特征即示功图的局部信息,还能考虑到神经网络高层的语义特征也即示功图的全局信息。
6.根据权利要求1所述的一种抽油机有效冲程的计算方法,其特征在于,所述根据模型M输出的凡尔开闭点像素坐标和归一化相关参数反推出实际有效冲程,具有如下步骤:
(10-1)将实时的示功图有效数据dr归一化为dnorm_r,其中:
dr={(xrj,yrj)|j=1,2,...,m}
dnorm_r={(xnorm_rj,ynorm_rj)|j=1,2,...,m}
两者间的转换关系如下:
(10-2)将归一化后的数据映射至w×h大小的网格内并绘制成示功图后输入神经网络回归模型,则示功图像素点(arj,brj)和归一化后的示功图悬点位移及载荷数据(xnorm_rj,ynorm_rj)的转换关系如下:
(arj,brj)=(xnorm_rj,ynorm_rj)×(w,h)
(10-3)结合(10-1)和(10-2)可得示功图像素点(arj,brj)和油井现场示功图悬点位移及载荷数据(xrj,yrj)的映射关系如下:
(10-4)神经网络回归模型输出的四个凡尔开闭点像素坐标如下:
lr={arA,brA,arB,brB,arC,brC,arD,brD}
结合(10-3)中建立的映射关系反推出四个凡尔开闭点实际悬点位移及载荷物理坐标如下:
lreal_r={x1,y1,x2,y2,x3,y3,x4,y4}
(10-5)根据(10-4)中的四个凡尔开闭点实际悬点位移及载荷物理坐标计算出示功图的有效冲程Spe如下:
Spe=min{x2-x1,x3-x4}。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110927362.5A CN113486599B (zh) | 2021-08-10 | 2021-08-10 | 一种抽油机有效冲程的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110927362.5A CN113486599B (zh) | 2021-08-10 | 2021-08-10 | 一种抽油机有效冲程的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113486599A CN113486599A (zh) | 2021-10-08 |
CN113486599B true CN113486599B (zh) | 2022-07-12 |
Family
ID=77946446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110927362.5A Active CN113486599B (zh) | 2021-08-10 | 2021-08-10 | 一种抽油机有效冲程的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113486599B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116838323B (zh) * | 2023-08-04 | 2023-12-05 | 大庆石油管理局有限公司 | 一种塔架式抽油机安全防护方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101660402B (zh) * | 2009-09-11 | 2013-06-12 | 河海大学 | 基于地面示功图凡尔开闭点的抽油井状态识别方法 |
CN104504236B (zh) * | 2014-11-24 | 2017-07-21 | 贵州航天凯山石油仪器有限公司 | 泵功图计量稳定性方法 |
CN104680023B (zh) * | 2015-03-12 | 2017-08-29 | 重庆科技学院 | 基于多目标决策的抽油机参数优化方法 |
CN107762456B (zh) * | 2017-10-24 | 2018-12-25 | 中国石油天然气股份有限公司华北油田分公司 | 一种抽油机变速运行智能控制方法及控制装置 |
CN112766301B (zh) * | 2020-12-11 | 2024-04-12 | 南京富岛信息工程有限公司 | 一种采油机示功图相似性判断方法 |
-
2021
- 2021-08-10 CN CN202110927362.5A patent/CN113486599B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113486599A (zh) | 2021-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111830408B (zh) | 一种基于边缘计算和深度学习的电机故障诊断系统及方法 | |
CN107941537B (zh) | 一种机械设备健康状态评估方法 | |
CN110570396B (zh) | 一种基于深度学习的工业产品缺陷检测方法 | |
CN109508630B (zh) | 一种基于人工智能识别水尺水位的方法 | |
CN109272123B (zh) | 一种基于卷积-循环神经网络的有杆泵工况预警方法 | |
CN109635806B (zh) | 基于残差网络的电表数值识别方法 | |
CN106595551A (zh) | 基于深度学习的输电线路覆冰图像中覆冰厚度检测方法 | |
CN109636171A (zh) | 一种区域植被恢复的综合诊断与风险评价方法 | |
CN110619479B (zh) | 光伏组件的故障判别模型建模方法及故障判别方法 | |
CN109058771B (zh) | 基于样本生成和间隔Markov特征的管道异常检测方法 | |
CN112508105A (zh) | 一种采油机故障检测与检索方法 | |
CN112507479B (zh) | 一种基于流形学习和softmax的石油钻机健康状态评估方法 | |
CN113486599B (zh) | 一种抽油机有效冲程的计算方法 | |
CN111402249B (zh) | 基于深度学习的图像演化分析方法 | |
CN114639064B (zh) | 一种水位识别方法及装置 | |
CN115655717A (zh) | 一种基于深度领域泛化网络的轴承故障诊断方法 | |
CN112766301A (zh) | 一种采油机示功图相似性判断方法 | |
CN117056865B (zh) | 一种基于特征融合的机泵设备运行故障诊断方法及装置 | |
CN112686181B (zh) | 一种基于插值轴心轨迹的水轮机故障诊断方法 | |
CN110196456A (zh) | 一种基于相似年灰色关联分析的中长期降雨径流预报方法 | |
CN110320802B (zh) | 基于数据可视化的复杂系统信号时序识别方法 | |
CN115063337A (zh) | 埋地管道智能维修决策方法及装置 | |
CN114048762B (zh) | 双注意力引导的旋转机械健康评估方法 | |
CN111931574B (zh) | 一种气动调节阀鲁棒故障诊断方法 | |
CN114120043A (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 |