基于Matlab的控压钻井压力控制建模方法

韩 煦， 韩玉安

(1.中国矿业大学(北京)机电与信息工程学院，北京 100083；2.中国石油集团大庆钻探工程公司钻井工程技术研究院，黑龙江大庆 163413)

基于Matlab的控压钻井压力控制建模方法

韩 煦1， 韩玉安2

(1.中国矿业大学(北京)机电与信息工程学院，北京 100083；2.中国石油集团大庆钻探工程公司钻井工程技术研究院，黑龙江大庆 163413)

为了提高控压钻井的控压精度，研究了基于Matlab的控压钻井压力控制建模方法。分析了节流阀的节流控压原理，根据伯努利方程建立了表征节流阀物理特性的方程;利用采集的节流阀节流控压试验数据，通过Matlab软件进行线性回归，得到了可计算节流阀控压的四参数公式；根据试验选用的节流阀开度与当量直径的特性曲线，采用Matlab软件进行拟合，得到了节流阀开度与节流面积的关系方程；根据节流阀开关试验得到的节流阀开度与开关时间数据，利用Matlab软件拟合出了开度与开关时间的关系方程，将上述各方程相结合得到了节流控压的理论模型。利用Matlab软件的半实物仿真插件将模型转化为通用实时代码，并采用模块化建模方法，编制了控压钻井自动控制软件，形成了控压模型。循环模拟系统验证显示，控制系统的响应时间小于10 s，控制精度达到了±0.2 MPa，满足控压钻井精细控压的需求，这证明基于Matlab的控压钻井压力控制建模方法是可行的。

控压钻井；控压方法；井口回压；控制模型

控压钻井过程中，需要根据采集的井口压力数据实时自动控制井口回压，进而动态控制井眼压力剖面使之在设定的压力范围内[1-3]。井口回压的自动控制是通过自动调节专用自动节流管汇的节流阀开度来实现的。控压钻井要实现井口回压的精确动态控制，就需要研究节流阀的节流特性，建立一套能够精确控制节流阀开度的方法和模型，并将模型嵌入控制软件中。

节流阀的压力控制是基于经典控制理论的单一输入和输出系统。根据节流阀节流特性，建立线性微分方程和基于拉普拉斯变换传递函数关系，形成理论推导的控压模型。理论计算结果与实际应用数据还存在一定误差，虽然通过工程试验可以不断去完善，使控压精度满足要求，但推导过程繁琐，控压精度验证也无法通过仿真实现。为了更加便捷和准确地建立节流控压模型，笔者在分析节流阀节流控压原理的基础上，应用伯努利方程，并结合选用节流阀的物理特性和大量的试验数据，应用Matlab软件的微积分算法、逻辑算法、图形绘制、数据处理、实时代码生成和仿真等便利软件工具[4-5]，建立了一套控压钻井系统控制模型。

1 节流压力控制模型的建立和验证

1.1 模型的建立

目前，控压钻井系统普遍采用节流阀实现节流控压，液控节流阀的控制原理(如图1所示)为：通过液控部分的三位四通电磁阀换位切换，控制液压油正反循环进出液压马达，驱动液压马达正反转，带动蜗轮蜗杆正反转，推拉节流阀阀芯往复运动，调节节流阀的开度，进而实现节流压力的控制。节流阀的节流控压是典型的流体力学问题，根据流体力学理论，可以建立一套节流控压的理论模型。

图1 液控节流阀控制原理示意Fig.1 Sketch of control principle of the hydraulic control throttle valve

节流压力的大小主要与流经节流阀流体的流量、流体密度和节流阀的节流面积有关[6]。井口返出的钻井液流经节流阀，产生的节流压力可应用伯努利方程进行计算，计算公式为[7]：

(1)

式中：Δp为节流阀两端进出口压力差，Pa；μ为流量系数；Q为流经节流阀流体的流量，m3/s；ρ为节流阀循环流体密度，kg/m3；A为节流阀节流面积，m2。

式(1)为理想状态下的计算公式，工程应用中由于流体的密度、流量和节流面积都处在动态变化中，并且由于钻井液中含有大量岩屑，其节流口的形状也在发生动态变化。所以，要针对具体的节流特性，对伯努利方程进行相应的修正，修正后的节流压力计算公式为：

(2)

式中：R为节流压差与密度、流量和节流面积线性关系的系数；a，b，c分别为密度、流量和节流面积的幂指数。

为便于Matlab软件处理，对式(2)两边同时取对数，可得：

lgΔp=lgR+algρ+blgQ+clgA

(3)

参数Δp、Q、ρ和A可以通过测量得到，这样就可以通过采集大量数据，进行线性回归后得到R、a、b和c等4个参数的值。为此，设计了如图2所示的节流循环控压试验系统，通过调节和改变流量、节流阀节流面积(开度)和钻井液密度等参数，采集了上万组流体密度、流量、节流面积和节流压力数据。

图2 节流循环控压试验系统示意Fig.2 Sketch of the throttling-circulating pressure control testing system

应用Matlab软件中的xlsread命令将数据代入Matlab数值计算拟合工具箱中，通过线性回归计算出未知量R，a，b和c的值，将其代入式(3)进行整理，可得出节流阀控压的四参数公式：

(4)

式(4)的通用性受到节流循环控压试验系统本身和所使用循环介质的限制，但是可以通过不断积累各种工况下的数据，应用该方法对R，a，b和c的值进行优化，进而不断提高该公式的通用性。

在控压钻井过程中，井口压力的控制是通过调整节流面积即调整节流阀开度来实现的，所以要确定开度与节流面积的关系。节流阀生产厂商给出了节流阀开度与当量直径的关系曲线(如图3所示)，可以通过当量直径计算出节流阀的节流面积，从而确定开度与节流面积的关系。利用Matlab软件的拟合工具箱cftool命令对有限的点进行拟合，得到节流阀开度与节流面积的关系曲线，如图4所示。

图3 选定的节流阀开度与当量直径的关系曲线Fig.3 Relationship between the equivalent diameter and the opening size of selected throttle valve

图4 节流阀节流面积与开度的拟合曲线Fig.4 Fitting curve of throttling area and opening size of throttle valve

从图4可以看出，随着节流阀开度不断加大，节流面积也不断增大，开度达80%以后节流阀的节流面积不变，开度小于80%时节流面积与开度之间近似呈线性关系，可应用Matlab软件中的polyfit命令进行数据拟合，得出关系式：

(5)

式中：k为节流阀开度。

由节流阀控制原理可知，它通过控制时间的脉冲来控制三位四通换向阀的给电时间，从而控制液压马达带动蜗轮蜗杆的旋转时间，进而调节节流阀的开度，因此需要建立时间与开度的关系。通过节流阀开关试验，使节流阀开度在30%～80%(节流阀特性曲线最佳范围)间调节，获得大量的节流阀开度k和开关时间t的对应数据，利用Matlab软件拟合工具中polyval显示数据的散点分布曲线，设定置信度为0.03(根据Matlab软件统计显示能够设置的最小值，不然不稳定点就会增多)，从而去掉试验过程中受外界因素影响的不稳定点(Matlab软件具有筛出不稳定数据的功能)。然后对这些稳定的数据应用Matlab拟合工具箱对节流阀开度与控制时间进行拟合，得出公式：

k=7.134t+0.071

(6)

式中：t为节流阀开关时间，s。

整理式(4)—式(6)，得到节流阀控压的理论模型，可表示为：

(7)

1.2 模型的验证

使用图2所示的节流循环控压试验系统，变换各种参数重新采集了2 500组数据，对式(7)的可信度进行校验分析。应用Matlab软件的rcoplot命令对式(7)进行验证，图5为验证公式的残差图(图5中，绿色部分表示数据在置信区间，红色线部分表示超出了置信区间)。图5中(-0.2，0.2)区间以外的点表示超出期望值的数据，通过2 500组后续测量的数据分析，模型的拟合置信度为97.36%，拟合的曲线与测量数据散点相差不大，说明式(7)准确可行，可以嵌入到控制模型中。

图5 控压理论模型验证残差图Fig.5 Residual plot of control formula verification

2 控压钻井系统控制模型的建立

2.1 控制方式

Matlab软件中的simulink采用模块组合方式来建模[8]，可以准确地创建动态系统的仿真模型，特别是对控压钻井这种复杂的非线性系统，它提供了一些成熟的系统模块(如非线性控制设计模块集、实时工作空间库、实时工作窗口目标库等)，其中real-time work半实物仿真是用simulink设计出来的控制器，直接去控制实际被控对象，通过半实物仿真过程来观察控制效果，如果效果不理想，则可以直接在simulink上调整控制器[9]的结构参数，直至获得满意的控制结果为止。

2.2 控制模型的建立

图6所示为控压钻井系统控制流程。在该流程中，监测单元采集节流控制装置的信号，数据接收及数模转换模块接收来自监测单元采集的流量、压力和密度等数据，经数模转换等处理后将数据传输到计算处理模块，计算得出控制量即开关节流阀的时间t，并将其传输给发送模块，经发送模块的发送，通过执行单元控制节流阀，实现节流控压，整套系统连续地采集数据，不断调整控制参数t，实现动态连续闭环控压。

图6 控压钻井系统控制流程Fig.6 Control flow chart of managed pressure drilling system

基于以上控压思路，应用Matlab软件中半实物仿真工具箱RTW建立了控制程序，形成了控压钻井系统控制模型，如图7所示。

图7 控压钻井系统控制模型Fig.7 Control model of managed pressure drilling system

控制程序的控制过程涉及5个控制参数，即流体流量、流体密度、节流阀当前开度、井口真实回压及计算的井口回压，数据通过采集板卡研华PCI-1710HG进行模拟信号的接收(分别来自Q7、F9、K9、P6支路和程序计算值)和A/D转换，建立数据接收模块并连接到板卡相应端口上接收数据，建立一个增益模块接收电压信号数据进行放大还原，同时添加一个选择开关模块筛选采集过程中受外界因素影响而产生的没有物理意义的负值。选择开关模块中间的是控制信号，控制信号与设定信号进行对比，决定上、下2个端口哪个能通过。

利用Matlab软件中的S-Function模块建立传递函数，将理论模型嵌入到控制模型中。采集到的压力、流量和密度数据，通过S-Function模块的计算，可求出设定的井口回压值对应的节流阀开度。计算出的节流阀开度与节流阀真实开度进行对比，系统按照软件程序计算结果自动控制节流开度，使真实值与计算值之间不断反馈修正，从而实现程序内部的二级反馈控制，使实时测量的井口回压越来越接近设定的井口回压，进而实现精细控压钻井。

3 控制模型的验证

为验证控制模型的准确性和有效性，应用压力循环系统(钻井泵、钻井液罐、节流控压装置和配套管线)、数据采集和传输系统对控制模型进行了地面模拟试验，0.5～4.0 MPa压力变化区间节流压力跟随设定压力曲线的控制模型控压效果如图8所示。

图8 控制模型控压效果曲线Fig.8 Control pressure effect curve of control model

从图8可以看出，随着设定压力的变化，节流阀开度连续调整，节流压力曲线跟随设定压力的变化较好，在各压力区间压降可以快速稳定地跟随目标压力，系统响应时间(节流压降曲线跟随设定压降曲线到重合的时间)不到10.0 s，同时控制精度[10]小于0.2 MPa(即节流压降曲线与设定压力曲线重合后波动幅度小于0.2 MPa)，模型控制精度较好。

4 结论与建议

1) 应用Matlab软件强大的数据处理功能，建立了选定节流阀的压力控制理论模型，建模方法便捷，准确性较高，置信度达到了97%以上。

2) 应用Matlab软件先进的控制算法以及控制策略，使控制系统的响应时间小于10.0 s，控制精度达到了±0.2 MPa。

3) 地面试验表明，自动控制软件已经达到要求，但还需要通过现场试验不断积累数据，并进行回归计算，使理论模型和控压模型更加完善，以满足精细控压钻井的需求。

References

[1] 王洋，韩来聚，步玉环，等.控压钻井井口回压控制技术仿真研究[J].石油机械，2014，42(2)：10-13. WANG Yang,HAN Laiju,BU Yuhuan,et al.The simulation of wellhead back pressure control technology for managed pressure drilling[J].China Petroleum Machinery,2014,42(2):10-13.

[2] 张兴全，周英操，刘伟，等.控压欠平衡钻井井口回压控制技术[J].天然气工业，2013，33(10)：75-79. ZHANG Xingquan,ZHOU Yingcao,LIU Wei,et al.Wellhead backpressure control in under-balanced and managed pressure drilling[J].Natural Gas Industry,2013,33(10):75-79.

[3] 苏勤，赵向阳.一种精细控压钻井流动模型的研究与应用[J].石油钻探技术，2013，41(1)：8-13. SU Qin,ZHAO Xiangyang.The research and application of a fine MPD flow model[J].Petroleum Drilling Techniques,2013,41(1):8-13.

[4] 张豪阳，李二欠，吕德瑾.基于MATLAB与AMESim的液压系统仿真特点[J].煤矿机械，2016，37(7)：163-165. ZHANG Haoyang，LI Erqian，LYU Dejin.Simulation characteristics of hydraulic system based on MATLAB and AMESim[J].Coal Mine Machinery,2016,37(7):163-165.

[5] 杨贵恒，文武松，张颖超，等.基于RTW的单项半桥逆变器SPWM控制算法实现[J].电源世界，2015(4)：21-24，32. YANG Guiheng，WEN Wusong，ZHANG Yingchao，et al.Realization of the SPWM algorithm for single-phase half-bridge inverter based on RTW[J].The World of Power Supply，2015(4):21-24,32.

[6] 周英操，杨雄文，方世良，等.精细控压钻井系统研制与现场试验[J].石油钻探技术，2011，39(4)：7-12. ZHOU Yingcao，YANG Xiongwen，FANG Shiliang，et al.Development and field test of PCDS-Ⅰprecise managed pressure drilling system[J].Petroleum Drilling Techniques，2011，39(4)：7-12.

[7] 刘琛，唐亮.控压钻井自动节流控制系统的研制与应用[J].新疆石油天然气，2013，9(4)：78-83. LIU Chen,TANG Liang.The research and application of MPD automatic throttle control system[J].Xinjiang Oil & Gas,2013,9(4):78-83.

[8] 袁亚锐，徐莉萍，任德志，等.金刚石压机位置控制数学模型的建立及Simulink仿真[J].液压与气动，2005(2)：55-57. YUAN Yarui，XU Liping，REN Dezhi，et al.Mathematical model for diamond press position control and its Simulink simulation[J].Chinese Hydraulics & Pneumatics,2005(2):55-57.

[9] 史勇，刘小勇，高翔，等.电液比列加载系统的模糊PID控制策略的研究[J].流体传动与控制，2014(4)：24-27. SHI Yong，LIU Xiaoyong，GAO Xiang，et al.Research on fuzzy-PID control strategy of electro-hydraulic loading system[J].Fluid Power Transmission and Control，2014(4):24-27.

[10] 金业权，孙泽秋，刘刚.控压钻井钻井液动节流压力控制系统仿真分析与试验研究[J].石油钻探技术，2013，41(2)：109-113. JIN Yequan，SUN Zeqiu，LIU Gang.Simulation analysis and experimental study of managed pressure drilling hydraulic throttle pressure control system[J].Petroleum Drilling Techniques，2013，41(2)：109-113.

[编辑 令文学]

A Pressure Control Modeling Approach for Managed Pressure Drilling Using Matlab

HAN Xu1,HAN Yu’an2

(1.CollegeofMechanicalElectronic&InformationEngineering,ChinaUniversityofMiningandTechnology(Beijing),Beijing，100083,China;2.DrillingEngineeringTechnologyResearchInstitute,CNPCDaqingDrilling&ExplorationEngineeringCorporation,Daqing,Heilongjiang，163413,China)

In this paper,a pressure control modeling approach using Matlab was studied to improve the pressure control accuracy of managed pressure drilling.Firstly,throttling control pressure principle of throttle valve was analyzed，the equation representing physical characteristics of throttle valve was established on the basis of Bernoulli’s equation.To do so,it was necessary to acquire throttling control pressure experiment data of throttle valve and to develop a four-parameter equation that could be used to calculate the pressure control of throttle valve.This was done by means of linear regressions using Matlab software.Second,the relation equation between the opening and throttling areas of the throttle valve was established by using Matlab software which perform goodness of fit analysis based on the characteristic curve of opening and equivalent diameter of the throttle valve selected in this experiment.Third,the opening and switching time data of the throttle valve was collected from its on-off experiment,and the relation equation between the opening and switching times was developed by using Matlab software fitting,and the theoretical model of throttling control pressure was finally established by combining above-mentioned equations.And fourth,the model was translated into the universal real time code by using the semi-physical simulation plug-in of Matlab software and the automatic control software for managed pressure drilling was prepared by adopting the modularized modeling method.In this way,the pressure control model was constructed.It was verified by a circulation analog system that this model could meet the fine pressure control requirements of managed pressure drilling with the response time of control system less than 10 s and the automatic control pressure accuracy ±0.2 MPa.It is proved that this Matlab based pressure control modeling approach for managed pressure drilling is viable.

managed pressure drilling;pressure control approach;wellhead back pressure;control model

2016-12-12；改回日期：2017-05-10。

韩煦(1996—)，男，黑龙江大庆人，测控技术与仪器专业在读本科生。E-mail:hhyyaa2011@sina.com。

10.11911/syztjs.201703012

TE249

A

1001-0890(2017)03-0067-05