CN111335812A - 一种钻井工具面角自校正闭环控制方法 - Google Patents
一种钻井工具面角自校正闭环控制方法 Download PDFInfo
- Publication number
- CN111335812A CN111335812A CN202010138385.3A CN202010138385A CN111335812A CN 111335812 A CN111335812 A CN 111335812A CN 202010138385 A CN202010138385 A CN 202010138385A CN 111335812 A CN111335812 A CN 111335812A
- Authority
- CN
- China
- Prior art keywords
- identification parameter
- parameter vector
- establishing
- face angle
- objective function
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005553 drilling Methods 0.000 title claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 239000013598 vector Substances 0.000 claims description 40
- 230000008859 change Effects 0.000 claims description 8
- 238000013178 mathematical model Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
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
- E21B7/00—Special methods or apparatus for drilling
- E21B7/04—Directional drilling
-
- 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
- E21B44/00—Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions
Abstract
Description
技术领域
本发明涉及钻采机械领域,尤其涉及一种钻井工具面角自校正闭环控制方法。
背景技术
在定向井钻井作业过程中,要保证实钻的井眼轨迹贴近设计的井眼轨道,现场工程技术人员,需要有丰富经验。在滑动模式下,首先要将钻柱提离井底并停止旋转,然后定向司钻才能通过摆工具面使井下马达弯头的朝向指向井眼轨迹设计所需的方向。摆工具面需要考虑钻头扭矩、钻柱反扭矩以及钻柱和井壁之间的摩阻,所以司钻要在地面上缓缓旋转钻柱,然后参考MWD测量的工具面变化。在较深井眼里,钻柱会吸收大量的扭矩,地面可能旋转几圈后井底的工具才转一圈。当合适的工具面被确认后,司钻会启动井下马达开始沿着目标方向滑动钻进。由于钻头切入地层产生的反扭矩会使工具面左右摆动,所以有时需反复调整工具面的方向。显然,这种靠司钻经验的摆工具面方式会影响钻井进度,增加钻井成本,易发生工具面漂移、托压问题。
比如,国家专利文献CN109098659A,公开了“一种滑动钻井工具面调整方法”,该发明步骤为:a)钻井参数采集;b)根据井眼轨迹要求,确定工具面中值φ及工具面波动范围δφ;c)完成防托压定向工具面调整的参数调试,拟合出扭矩-角速度曲线;d)根据扭矩-角速度曲线,调试井口工具旋转钻柱的能量参数,系统自动选择钻柱旋转的步骤,进行防托压定向钻井;e)持续更新扭矩-角速度曲线,结合钻井参数曲线系列,优化钻井参数。该发明在拟合扭矩-角速度曲线的具体步骤中,需要人工监测定向工具面变化情况,确认工具面顺时针转动或逆时针转动,计算工具面转动角速度。然而通过人为方式进行钻井工具面调整,误差大,易发生工具面漂移、托压情况。
发明内容
本发明的目的在于提供一种钻井工具面角自校正闭环控制方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明采用的技术方案如下:
一种钻井工具面角自校正闭环控制方法,包括步骤:
S1)采集输入数据及输出数据,输入数据为电机转速,输出数据为实时采集的工具面角;
S2)根据输入数据及所述输出数据建立CARMA模型;
S4)设置增益矩阵L(t)和协方差矩阵P(t),进行参数估计,获得满足要求的辨识参数;
S5)建立控制目标函数,根据控制目标函数以及辨识参数对电机进行控制,得到所需的工具面角,返回步骤S1)。
由于钻具系统不能够建立精确的数学模型,所以本发明采用自适应辨识模型进行辨识控制,CARMA模型为受控自回归滑动平均模型,根据系统的输入(电机转速)和系统的输出(实时采集的工具面角)实时变化,实时地辨识对象模型。
进一步的,步骤S2)中CARMA模型为A、B和C均为单位滞后算子z-1的常系数时不变多项式,z-d表示纯时延,d表示滞后量,A=1+a1z-1+…+anaz-na,na为输出的阶次,B=b1z-1+…+bnbz-nb,b0≠0,nb为输入的阶次,C=c0+c1z-1+…+cncz-nc,nc为噪声的阶次,ξ(t)为标准正态分布的白噪声序列。
S31)将CARMA模型转换成差分方程:
y(t)+a1y(t-1)+...+anay(t-na)=b1u(t-d-1)+...bnbu(t-d-nb)+c0ξ(t)+c1ξ(t-1)+...+cncξ(t-nc);
θ=(a1,...,ana,b1,...,bnb,c0,c1,...,cnc)T。
本发明在辨识过程中考虑了噪声模型的参数,噪声为具有有理谱密度的零均值平稳序列,即采用了白噪声序列,适用于CARMA模型。
进一步的,步骤S4)中设置增益矩阵L(t)和协方差矩阵P(t),进行参数估计,包括步骤:
S42)计算t时刻需要进行辨识的辨识参数向量
S44)设置递推条件,包括步骤:
S443)按递推时间顺序依次获得最近n次的辨识参数向量的相对变化量结果:ωj、ωj-1、…、ωj-n,ωj表示当前进行递推的辨识参数向量相对变化量,ωj-1表示上一次进行递推时的辨识参数向量相对变化量;
本发明进行参数估计过程中采用了递推最小二乘法,利用递推最小二乘法具有无偏性、有效性以及一致性的特点,根据最新得到的钻井输入及输出数据,在原有参数估计值的基础上,对预测误差乘上增益矩阵L(t),不同程度地修正辨识参数向量中的每个辨识参数,进行动态参数调整,这样不断地递推更新,最后得到最优辨识参数。递推最小二乘法为在线算法,能够在线跟踪被估计参数,而一般最小二乘递推估计算法缺少在线跟踪被估计参数的功能,本发明相对于一般最小二乘递推估计算法,实用性更强。另外,本发明通过设置遗忘因子t0,增加了新测的数据的作用,相对地减少了旧的历史数据的影响,使得参数估计过程随着数据的累积,更多地依靠当前数据,从而遗忘旧的历史数据,减少旧的历史数据对建模的影响,更好地适应参数的变化。
对辨识参数向量和协方差矩阵P(t)需要设置初值。还有另外一种方法:先取最初N组数据,包括N个观测向量和对应的输出值y,再对参数进行成批处理,以此作为初值,再从N+1时刻进行递推计算。而本发明将辨识参数向量的初值设置为0,p0为充分大的实数,这样的设置使得本发明相对于另外一种方法,递推过程计算量小,便于应用。
进一步的,步骤S5)中建立控制目标函数,根据所述控制目标函数以及所述辨识参数对电机进行控制,包括步骤:
F=1+f1z-1+...+fd-1z-d+1,G=g0+g1z-1+...+gn-1z-n+1;
S53)将到t时刻为止的所有输入输出观测数据记为{Y,U}={y(t),y(t-1),...,u(t),u(t-1),...},基于{Y,U}对t+d时刻实际输出工具面角y(t+d)进行预测;
S54)获得使预测误差方差最小的d步最优预测值
S55)建立控制目标函数Jmin=Emin{|y(t+d)-yr|2},使实际输出工具面角y(t+d)与希望输出工具面角yr之间的方差最小;
本发明采用了最小方差自校正控制解决定向井钻作业工程中对工具面角进行控制的问题,由于钻井作业存在纯时延,当前的控制作用要滞后于采样周期才能影响输出。因此,要使输出方差最小,就需要提前步对输出量做出预测,然后,根据所得的预测值来设计所需的控制。然后通过连续不断的预测和控制,就能保证稳态输出方差最小。本发明引入了方程,进行最优d步预测,获得使预测误差方差最小的d步最优预测值yr=y*(t+d|t),最后获得电机控制量u(t),通过电机控制量u(t)对电机进行控制,从而获得需要的工具面角。
本发明的有益效果是:本发明采用CARMA模型,根据系统的输入(电机转速)和系统的输出(实时采集的工具面角)实时变化,实时地辨识对象模型,解决了建立钻具系统数学模型变量不确定性及模型形式不确定性的问题;并通过递推最小二乘法在线进行参数动态调整,获得最优参数,针对被控对象模型受随机扰动,致使模型参数变化,进行实时的辨识出模型参数;本发明引入的C/A方程,能够将预测方程简化,以便进行求解;通过d步最优控制,避免了控制量输出过大,系统失衡。
附图说明
图1是实施例一获得最优辨识参数的流程示意图。
图2是实施例一控制框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
实施例一,一种钻井工具面角自校正闭环控制方法,包括步骤:
S1)采集输入数据及输出数据,输入数据为电机转速,输出数据为实时采集的工具面角;
S2)根据输入数据及所述输出数据建立CARMA模型;CARMA模型为A、B和C均为单位滞后算子z-1的常系数时不变多项式,z-d表示纯时延,d表示滞后量,A=1+a1z-1+…+anaz-na,na为输出的阶次,B=b1z-1+…+bnbz-nb,b0≠0,nb为输入的阶次,C=c0+c1z-1+…+cncz-nc,nc为噪声的阶次,ξ(t)为标准正态分布的白噪声序列。
S31)将CARMA模型转换成差分方程:
y(t)+a1y(t-1)+...+anay(t-na)=b1u(t-d-1)+...bnbu(t-d-nb)+c0ξ(t)+c1ξ(t-1)+...+cncξ(t-nc);
θ=(a1,...,ana,b1,...,bnb,c0,c1,...,cnc)T。
S4)设置增益矩阵L(t)和协方差矩阵P(t),进行参数估计,如图1所示,包括步骤:
S42)计算需要进行辨识的辨识参数向量
S44)设置递推条件,包括步骤:
S443)按递推时间顺序依次获得最近n次的辨识参数向量的相对变化量结果:ωj、ωj-1、…、ωj-n,ωj表示当前进行递推的辨识参数向量相对变化量,ωj-1表示上一次进行递推时的辨识参数向量相对变化量;
获得满足要求的辨识参数。
S5)建立控制目标函数,根据控制目标函数以及辨识参数对电机进行控制,钻具系统的控制框图如图2所示,包括步骤:
F=1+f1z-1+...+fd-1z-d+1,G=g0+g1z-1+...+gn-1z-n+1;
S53)将到t时刻为止的所有输入输出观测数据记为{Y,U}={y(t),y(t-1),...,u(t),u(t-1),...},基于{Y,U}对t+d时刻实际输出工具面角y(t+d)进行预测;
S54)获得使预测误差方差最小的d步最优预测值
S55)建立控制目标函数Jmin=Emin{|y(t+d)-yr|2},使实际输出工具面角y(t+d)与希望输出工具面角yr之间的方差最小;
得到所需的工具面角,返回步骤S1)。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明提供了一种基于最小方差自校正的工具面闭环控制算法,通过本发明能够自动判断井下工具面的偏移情况并且发送指令至顶驱,实现滑动钻进过程对工具面的实时监控和自动控制。
以上仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (6)
4.根据权利要求1或3所述的一种钻井工具面角自校正闭环控制方法,其特征在于,步骤S4)中设置增益矩阵L(t)和协方差矩阵P(t),进行参数估计,包括步骤:
S42)计算t时刻需要进行辨识的辨识参数向量
S44)设置递推条件,包括步骤:
S443)按递推时间顺序依次获得最近n次的辨识参数向量的相对变化量结果:ωj、ωj-1、…、ωj-n,ωj表示当前进行递推的辨识参数向量相对变化量,ωj-1表示上一次进行递推时的辨识参数向量相对变化量;
6.根据权利要求3或5所述的一种钻井工具面角自校正闭环控制方法,其特征在于,步骤S5)中建立控制目标函数,根据所述控制目标函数以及所述辨识参数对电机进行控制,包括步骤:
F=1+f1z-1+...+fd-1z-d+1,G=g0+g1z-1+...+gn-1z-n+1;
S53)将到t时刻为止的所有输入输出观测数据记为{Y,U}={y(t),y(t-1),...,u(t),u(t-1),...},基于{Y,U}对t+d时刻实际输出工具面角y(t+d)进行预测;
S54)获得使预测误差方差最小的d步最优预测值
S55)建立控制目标函数Jmin=Emin{|y(t+d)-yr|2},使实际输出工具面角y(t+d)与希望输出工具面角yr之间的方差最小;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010138385.3A CN111335812B (zh) | 2020-03-03 | 2020-03-03 | 一种钻井工具面角自校正闭环控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010138385.3A CN111335812B (zh) | 2020-03-03 | 2020-03-03 | 一种钻井工具面角自校正闭环控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111335812A true CN111335812A (zh) | 2020-06-26 |
CN111335812B CN111335812B (zh) | 2021-08-31 |
Family
ID=71179882
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010138385.3A Expired - Fee Related CN111335812B (zh) | 2020-03-03 | 2020-03-03 | 一种钻井工具面角自校正闭环控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111335812B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100101860A1 (en) * | 2008-10-29 | 2010-04-29 | Baker Hughes Incorporated | Phase Estimation From Rotating Sensors To Get a Toolface |
CN105102762A (zh) * | 2013-03-29 | 2015-11-25 | 普拉德研究及开发股份有限公司 | 钻井工具面的闭环控制 |
CN106522924A (zh) * | 2016-11-15 | 2017-03-22 | 北京恒泰万博石油技术股份有限公司 | 一种随钻测量中方位角的获取方法 |
CN107515001A (zh) * | 2017-09-05 | 2017-12-26 | 中国石油大学(华东) | 旋转导向稳定平台重力工具面角动态测量方法及装置 |
US20200063546A1 (en) * | 2017-08-10 | 2020-02-27 | Motive Drilling Technologies, Inc. | Apparatus and methods for uninterrupted drilling |
-
2020
- 2020-03-03 CN CN202010138385.3A patent/CN111335812B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100101860A1 (en) * | 2008-10-29 | 2010-04-29 | Baker Hughes Incorporated | Phase Estimation From Rotating Sensors To Get a Toolface |
CN105102762A (zh) * | 2013-03-29 | 2015-11-25 | 普拉德研究及开发股份有限公司 | 钻井工具面的闭环控制 |
CN106522924A (zh) * | 2016-11-15 | 2017-03-22 | 北京恒泰万博石油技术股份有限公司 | 一种随钻测量中方位角的获取方法 |
US20200063546A1 (en) * | 2017-08-10 | 2020-02-27 | Motive Drilling Technologies, Inc. | Apparatus and methods for uninterrupted drilling |
CN107515001A (zh) * | 2017-09-05 | 2017-12-26 | 中国石油大学(华东) | 旋转导向稳定平台重力工具面角动态测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111335812B (zh) | 2021-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11421520B2 (en) | Drilling parameter optimization for automated well planning, drilling and guidance systems | |
EP2592224B1 (en) | Methods and systems for drilling | |
CN105283629B (zh) | 钻柱在定向钻井的滑动模式期间的优化的旋转 | |
US10370902B2 (en) | Downhole steering control apparatus and methods | |
US10612989B2 (en) | Method and apparatus for automated drilling rig sheave friction calibration | |
US10767462B2 (en) | Method and system for performing automated drilling of a wellbore | |
US11118440B2 (en) | Method, system and computer-readable medium for automatically controlling a drilling operation | |
US11585203B2 (en) | Closed loop control of drilling curvature | |
US11408268B2 (en) | Autonomous directional drilling directional tendency estimation | |
Zhang et al. | Overview of rotary steerable system and its control methods | |
CN107201877B (zh) | 一种旋转导向钻井的闭环控制方法及系统 | |
US20230046043A1 (en) | Adaptive trajectory control for automated directional drilling | |
CN111335812B (zh) | 一种钻井工具面角自校正闭环控制方法 | |
CN115142791A (zh) | 旋转导向系统自动曲率控制方法和系统、计算机设备 | |
CN114856540A (zh) | 一种基于在线学习的水平井机械钻速随钻预测方法 | |
CN116498294A (zh) | 一种弯外壳螺杆钻具组合智能定向闭环调控方法及系统 | |
CN116607929A (zh) | 一种煤田钻探深孔钻探中的钻孔自动纠偏方法 | |
Spencer et al. | Estimation and control for efficient autonomous drilling through layered materials | |
US11802473B2 (en) | Methods, systems, and media for controlling a toolface of a downhole tool | |
US11585205B2 (en) | Methods, systems, and media for controlling a toolface of a downhole tool | |
CN111335876A (zh) | 一种石油钻井井眼轨迹自适应跟踪预测控制方法 | |
CN114856428B (zh) | 一种滑动导向钻井前馈控制方法、前馈工控机及系统 | |
WO2023114938A1 (en) | Drilling with miso, toolface wagging, return to neutral, and spindle reaction time improvements |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210831 |