基于支持向量回归的抽油机井检泵周期预测研究.pdf

2023 年第 8 期计算机与数字工程收稿日期：2023年2月8日，修回日期：2023年3月20日基金项目：国家自然科学基金项目（编号：61801517）；中央高校基本科研业务专项经费（编号：19CX02029A，19CX02027A）资助。作者简介：邓杰，男，硕士研究生，研究方向：油田大数据分析、深度学习。刘新平，男，博士研究生，副教授，研究方向：智能感知与计算。杨鹏磊，男，硕士研究生，研究方向：数据驱动、预测控制。1引言目前，在我国机械采油工艺中，有杆泵采油方式在原油开采中占据着十分重要的地位。有效降低油井检泵率，提高油井利用率是采油工程控制成本的一个有效途径。研究人员在抽油机井检泵预测方面做了大量的研究工作，分析了影响油井检泵周期的因素并提出了预防措施15。影响抽油机井检泵因素较多。其中，文献 67 从抽油杆寿命方面进行了深入的研究，文献 8 研究了抽油机杆悬点载荷对检泵周期的影响，文献 9 开展了抽油机游动凡尔罩断裂机理及防治对策研究，文献 10 分析了抽汲参数调整对油井检泵率的影响。由于抽油机井生产参数较多（实际检测的生产参数20余种），各类参数与抽油机井的检泵因素有着直接关联关系，但是各类生产参数之间以及各类检泵因素之间存在耦合影响和非线性因素，难以判断各类生产参数对抽油机检泵周期的影响权重。同时，基于机理分析的预测模型需要与生产实际情况相结合，准确率和有效性有待进一步提高。因此，结合现场经验，从数据挖掘的角度出发，利用油田生产数据和油井作业数据，建立抽油机井检泵周期与各生产基于支持向量回归的抽油机井检泵周期预测研究邓杰刘新平杨鹏磊（中国石油大学（华东）计算机科学与技术学院青岛266580）摘要为了降低抽油机井检泵频率，合理安排检泵作业，减少生产作业成本，开展了抽油机井检泵周期预测方法研究。通过对抽油机井生产数据和作业数据的预处理，生成有效数据集，采用灰度关联方法优选主要影响参数，选取支持向量回归方法建立了抽油机井主要生产参数与检泵周期的关系模型，基于生产参数预测检泵周期，获得了较好的效果，证明了方法的有效性和可行性，同时，为进一步指导油田作业提供了实验基础和理论依据。关键词抽油机井；检泵周期；支持向量回归；数据预处理中图分类号TE355DOI：10.3969/j.issn.1672-9722.2023.08.038Research on Pump Detection Period Predicting Based onSupport Vector RegressionDENG JieLIU XinpingYANG Penglei（College of Computer Science and Technology，China University of Petroleum（East China），Qingdao266580）AbstractIn order to reduce the frequency of pump detection period，reasonably arrange the operation of pump detection andreduce the cost of production，the research on the prediction method of pump detection period in pumping well is carried out.Through the preprocessing of production data and operation data of pumping well，a valid data set is generated，the gray correlationmethod is used to select the main influencing parameters，the main production parameters of pumping well and pumping detectionperiod the relationship model are established by support vector regression method，through pump detection period predicting on production parameters，the good effect is obtained，and the validity and feasibility of this method are proved.At the same time，the experimental basis and theoretical basis are improved for further guidance of oilfield operation.Key Wordspumped well，pump detection period，SVR，data preprocessingClass NumberTE355总第 406期2023 年第 8 期计算机与数字工程Computer&Digital EngineeringVol.51No.81893第 51 卷要参数的关系模型预测抽油机井检泵周期，对有效降低检泵率具有重要的研究意义。常用的数据挖掘技术主要涉及机器学习和深度学习11，如线性回归、支持向量机（SVM）12、卷积神经网络13、循环神经网络14等算法。其中，支持向量回归（SVR）1516是支持向量机在回归领域的应用，它在解决非线性和高维问题中表现出许多特有的优势，被广泛应用于各种预测问题。本文采用SVR算法开展抽油机井检泵周期的研究，通过对抽油机井生产数据的预处理，基于灰度关联方法优选主要生产参数，建立抽油机井检泵周期与主要生产参数的关系模型。2数据预处理2.1数据采集与抽油机井检泵因素相关的数据来源于抽油机井生产数据和检泵作业数据。主要包括生产时间、日产液量、日产油量、含水率、上行电流、下行电流、油压、套压、流压、泵径、泵深、冲程、冲次、动液面、累产油量、累产水量、最大载荷、最小载荷、检泵作业时间等。由于抽油机井原始生产数据存在多源、异构、缺失、重复、不准确以及各类生产参数量纲不同等情况。需要对数据预处理后开展实验研究，有效提高模型的准确性。2.2数据预处理1）数据清洗统计分析各类生产参数的特征，掌握各类数据的范围，有效去除各类奇异数据和错误数据，提高模型的有效性和准确性。抽油机井的生产参数较多，以日产液量数据统计为例，如图1所示，某一区域的抽油机井单井日产液量集中在1545（t/d）范围内。因此，当日产液量超出这个范围后，需要分析产液量降低或是升高的因素。从而统计正常生产情况下的抽油机井的编号及数量。4003002001000井口数量0102030405060日产液量/（t/d）图1日产液量数据分布图基于此方法，统计生产参数正常情况下的油井的型号和数量，清洗掉奇异数据，生成有效数据体，作为训练和预测数据，提高预测模型的准确性和实用性。同时，基于程序判断去除重复数据，若某一属性缺失值较多，则直接删除该属性，若缺失值较少，则进行填充。对于数值型变量，采用中位数进行填充，对于类别型数据直接填充为“空”。2）去量纲抽油机井生产参数较多，数据单位不统一，如：体积单位（/m3），压力单位（/Mpa），时间单位（/d）等。为了便于不同单位或量级的指标能够进行比较和加权。将数据统一变换为无单位（统一单位）的数据集，采用最大-最小标准化方法消除量纲：x*=x-xminxmax-xmin（1）式（1）中：x为数据清洗后的样本数据，xmin，xmax为样本数据中最大值和最小值，x*表示归一化后的样本数据。3主要参数提取各类生产参数对抽油机井检泵周期影响权重不同。如果将全部生产参数作为自变量会增加计算负担，提高模型的复杂程度，降低准确率。因此，采用灰度关联方法优选与检泵周期关联性较大的参数做为自变量。通过灰度关联分析算法，可得到特征值（子序列）与周期（母序列）关联度大小的排序，关联度越高则表示两个因素变化的趋势具有强一致性，即同步变化程度越高。子序列的各个指标与母序列的关联系数计算方法如下：i(k)=minimink|x0(k)-xi(k)|+p maximaxk|x0(k)-xi(k)|x0(k)-xi(k)|+p maximaxk|x0(k)-xi(k)|（2）式（2）中，为分辨系数，01，若越小，关联系数间差异越大，区分能力越强，通常取0.5；x0(k)和xi(k)分别表示母序列第k个数和子序列第i个特征值的第k个数；i(k)则表示第i个特征的第k个值与母序列第k个值的关联系数。通过关联系数来计算关联度再进行最后的排序。其中关联度得计算如下：ri=1nk=1ni(k)（3）式（3）中，ri为第i个特征与母序列的关联度大小，n为样本数量。基于计算结果，对所有ri进行排序即可得到关联度排序。通过使用灰度关联分邓杰等：基于支持向量回归的抽油机井检泵周期预测研究18942023 年第 8 期计算机与数字工程析算法选择的主要参数为日产液量、冲程、冲次、泵径、泵深、沉没度、含水率、最大载荷、最小载荷。其所选参数与因变量之间关联度如表1所示。表1特征参数关联度表参数关联度参数关联度冲程0.78冲次0.76泵径0.69泵深0.63沉没度0.73含水率0.68最大载荷0.74最小载荷0.61日产液量0.64驱替方式三元驱4预测模型基于SVR的油井检泵周期预测模型的函数为F(x)=w(x)+b（4）式中，w为权值向量，(x)为非线性映射函数，b为偏置向量。由于SVR存在容忍偏差，于是SVR问题可形式化为minwb12|w|2+Ci=1mLoss（5）其中，C为正则化常数，Loss为损失函数：Loss=0if|f(xi)-yi|f(xi)-yi|-otherwise（6）为了确保大部分数据参与模型训练，引入松弛变量i和i*，则式（5）优化为minwb12|w|2+Ci=1m（i+i*）s.t.f(xi)-yi+iyi-f(xi)+i*i0i*0i=12m.（7）由拉格朗日乘子法可得拉格朗日函数：L=12|w|2+Ci=1m（i+i*）-i=1muii-i=1mui*i*+i=1mai(f(xi)-yi-i)+i=1mai*(yi-f(xi)-i*)（8）在KKT条件下，拉格朗日的对偶形式为maxaa*i=1myi(ai*-ai)-(ai*+ai)-12i=1mj=1m(ai*-ai)(aj*-aj)xTixjs.t.i=1m(ai*-ai)=00aiai*C.（9）将其对偶形式求解获得回归函数为F(x)=w(x)+b=i=1m(a*i-ai)(xi)(x)+b=i=1m(a*i-ai)K(xix)+b（10）其中，K(xixj)=(xi)T(xj)为核函数。在SVR周期预测模型中，核函数K（x，x）的类型对模型的性能影响较大，可以通过比较不同核函数的性能情况来选择最佳的和函数类型。5实验结果将原始数据经过清洗、归一化以及参数优选后开展周期预测建模实验。对数据样本进行随机划分，训练集占70%，测试集占30%。如果经过生产参数优选，控制算法参与遍历优化后仍不能达到预期准确率（预期准确率为85%以上），重新开展数据预处理和主要参数提取步骤，过程如图2所示。图2模型构建流程图实验分别测试了在不同核函数下以及不同的算法下的模型准确率。常见的核函数类型有线性核函数、多项式核函数、径向基RBF核函数和Sigmoid核函数。当选取不同的核函数时，模型准确率如表2所示。表2核函数的选取kernelC=1.0，gamma=0.1，Epsilon=3Linear76.55%Poly80.19%Rbf84.19%Sigmoid82.35%选择SVR模型的核函数为Rbf，遍历SVR模型参数 C的范围为（0.01，100）此时参数 gamma设置为默认值，即：1/n_features，其中 n_features 为输入特征数量，从图 3可知当 C为 9时准确率最高，为86.55%。确定 C 值后，再遍历参数 gamma 的范围为（0.01，10），从图4可知gamma为1.8时结果最优。确定各项参数后，仿真结果如图5所示，模型1895第 51 卷预测结果准确率高达90.76%。86858483预测准确率020406080100参数C图3遍历参数C实验图91908988878685预测准确率0246810参数gamma图4遍历参数gamma实验图1300120011001000900800700600500周期/天020406080100数据量真实值预测值图5SVR预测示意图基于实验数据，测试了多元线性回归算法，BP神经网络算法，并将模型预测结果与SVR算法相比较。实验结果如表3所示。表3三种算法预测误差对比情况预测算法SVR多元线性回归BP神经网络平均绝对误差值/天4.1319.2836.064相对误差平均值/%0.7422.3911.163模型的平均绝对误差值：MAE(x)=1mi=1m|h(xi)-yi|（11）模型的平均相对误差值：MRE(x)=1mi=1m|h(xi)-yi|yi（12）其中 m 为数据量，h(xi)和 yi分别为预测值和真实值。基于SVR的检泵周期预测结果准确率率为达到了预期目标，与多元线性回归算法和BP神经网络算法相比较，获得了较好的预测结果。6结语本文从数据挖掘的角度出发，采用灰度关联算法，优选了抽油机井主要生产参数，基于SVR算法开展了抽油机井检泵周期预测的实验研究，建立了检泵周期与生产参数之间的关系模型，获得了较好的实验结果。将大数据分析方法与油田措施业务相结合，充分挖掘数据价值的可行性和有效性，为智慧油田建设提供了新的研究方向，拓展了解决工程技术难题的思路和方法。同时，在实验过程中，由于数据噪声的影响，人工检泵质量不同，多因素检泵等因素共存，对检泵周期的预测和分析影响极大。参 考 文 献1王汉雄.油井检泵原因分析及延长检泵周期的措施J.化学工程与装备，2020（11）：110-111.WANG Hanxiong.Cause analysis of oil well pump inspection and measures to prolong pump inspection period J.Chemical Engineering&Equipment，2020（11）：110-111.2任军，刘金森，张慧敏.影响油井检泵周期的原因分析及预防措施的研究 J.中国石油和化工标准与质量，2014，34（12）：42-43.REN Jun，LIU Jinsen，ZHANG Huimin.Analysis of thecauses and preventive measures affecting the period of oilwell inspection pump J.China Petroleum and ChemicalStandard and Quality，2014，34（12）：42-43.3王平.油井检泵周期影响因素分析及预防措施探讨J.石化技术，2018，25（5）：153.WANG Ping.Factors and preventive measures of oil wellinspection period J.Petrochemical Industry Technology，2018，25（5）：153.4刘刚.油井检泵原因及延长检泵周期的措施 J.化学工程与装备，2021（01）：115-116.LIU Gang.Reasons for oil well pump inspection and measures to prolong pump inspection period J.Chemical Engineering&Equipment，2021（01）：115-116.5赵妤蕾.影响油井检泵周期原因及提升泵效措施研究J.石油知识，2020（01）：54-55.邓杰等：基于支持向量回归的抽油机井检泵周期预测研究18962023 年第 8 期计算机与数字工程ZHAO Yulei.Study on the reasons affecting the period ofoil well inspection pump and the measures to improvepump efficiencyJ.Petroleum Knowledge，2020（01）：54-55.6Ulmanu V V，Ghofrani R.Fatigue life prediction methodfor sucker rods based on local concept J.Erdoel ErdgasKohle，2001，117（4）：189-195.7李德君，王伟，庞斌，等.煤层气井用抽油杆腐蚀疲劳寿命的影响因素J.材料热处理学报，2017，38（3）：121-127.LI Dejun，WANG Wei，PANG Bin，et al.Influence factors on corrosion fatigue life of sucker rods used in coalbedmethane wellJ.Transactions of Materials and HeatTreatment，2017，38（3）：121-127.8朱治国.抽油机最优悬点载荷利用率的判断方法研究J.石油石化节能，2018，8（11）：4-5，8.ZHU Zhiguo.Study on the Optimal Suspension Load Utilization of Pumping Unit J.Energy Conservation in Petroleum&Petrochemical Industry，2018，8（11）：4-5，8.9李亚峰，张浩毅.抽油泵游动凡尔罩断裂机理及防治对策研究 J.内蒙古石油化工，2010，36（16）：26-28.LI Yafeng，ZHANG Haoyi.Study on the Fracture mechanism and prevention countermeasures of the Moving ValJar of oil pump J.Inner Mongolia Petrochemical Industry，2010，36（16）：26-28.10张世超.浅析抽油机井的抽汲参数的调整与优化 J.化学工程与装备，2016（05）：111-112.ZHANG Shichao.The adjustment and optimization ofswabbing parameters of pumping Wells are analyzed J.ChemicalEngineering&Equipment，2016（05）：111-112.11LeCun Y，Bengio Y，Hinton G.Deep Learning J.Nature，2015，521（7553）：436-444.12Md A S，Sultana N，Hossain S M Z.Bayesian optimization algorithm based support vector regression analysisfor estimation of shear capacity of FRP reinforced concrete membersJ.Applied Soft Computing Journal，2021，105（20）：107281.13 Krizhevsky A，Sutskever I，Geoffrey H E.ImageNetClassification with Deep Convolutional Neural NetworksJ.Communications of the ACM，2017，60（6）：84-90.14Zaremba W，Sutskever I，Vinyals O.Recurrent NeuralNetwork Regularization J.Eprint Arxiv，2014.15Sabzekar M，Hasheminejad S M H.Robust regression using support vector regressionsJ.Chaos，Solitons andFractals，2021，144（2）：110738.16Liu Q，Liu W，Mei J，et al.A New Support Vector Regression Model for Equipment Health Diagnosis withSmall Sample Data Missing and Its ApplicationJ.Shock and Vibration，2021（2）：1-11.19刘胜久，李天瑞，珠杰.Zipf定律与网络信息计量学J.中文信息学报，2015，29（04）：89-94.LIU Shengjiu，LI Tianrui，ZHU Jie.Zipf Law and Network Information MetrologyJ.Journal of Chinese Information，2015，29（04）：89-94.20马剑锋.基于Web结构挖掘中HITS算法研究 J.计算机产品与流通，2020（02）：252.MA Jianfeng.Research on HITS algorithm in web-basedstructure miningJ.Computer Products and Circulation，2020（02）：252.21李开拓，彭慧，周晓锋，等.基于多种相关性度量的特征选择方法研究 J.小型微型计算机系统，2017，38（04）：696-700.LI Kaituo，PENG Hui，ZHOU Xiaofeng，et al.Researchon feature selection method based on multiple correlationmeasuresJ.Miniature Microcomputer System，2017，38（04）：696-700.（上接第1892页）1897