-
底部钻具组合(BHA)是井眼轨迹控制的核心部件,其造斜率的准确预测和控制是复杂结构井工程的关键技术之一[1-10]。国内外学者对导向钻具造斜率的预测方面做了大量的研究,其预测方法主要经历了几何法[11]、力学法[12-13]、轨迹法[14]、极限造斜率法[15]的过程。
几何法假设BHA处于小变形状态或钻具组合长度较短的情况下,将钻头、下稳定器、上切点当作同一圆弧上的3个点,求解该圆弧的曲率作为BHA的造斜率,然而该方法并未考虑钻具刚度及钻井参数[16-17]。力学法采用静力学计算导向钻具组合钻头处的机械合力,将合力作用线与井眼轴线的夹角作为导向工具造斜率的预测结果,但该方法并未考虑钻头切削性能和地层各向异性[18-21]。轨迹法通过计算钻头侧向力和钻压联合作用下的合位移,确定合位移方向与钻头轴线的夹角,阐明了钻进过程中的运动趋势,但是钻进趋势角与造斜率之间缺少明确的定量关系,不利于造斜率的准确评估。极限曲率法采用了最小势能原理,将钻头侧向力为0或趋近于0时的井眼曲率作为导向工具造斜率的评价标准,但是未充分考虑钻头切削性能和地层各向异性的影响。总之,目前造斜率的计算方法仍存在不足,加之参数的不确定性,使得导向工具造斜率的理论计算值与实际值存在较大的误差。
本文在前人研究的基础上,综合考虑工具结构、钻井参数、轨迹参数及钻头-地层相互作用等影响,建立了基于零侧向钻速准则的滑动钻进造斜率计算方法,并从“杠杆效应”和“钟摆效应”的角度揭示了造斜率的影响机理,为导向钻具组合的结构设计、造斜率预测、钻具优选及钻井参数优化等提供了分析手段。
-
在定向钻进过程中,导向工具受到钻压、扭矩、支反力、自重等载荷作用。将BHA分解成多个梁的组合,针对每一段梁,考虑纵向和横向载荷的作用来计算梁的变形,建立BHA的三维力学模型。由于钻头处在井斜平面、方位平面内都会存在侧向力的作用,在钻压、侧向力的作用下,井下管柱所受多种载荷的大小、方向也会进行相应变化,钻出的井眼轨迹为三维曲线,因此可将三维力学问题转化为两个二维力学问题进行求解,控制方程为
$$ \frac{{{{\text{d}}^4}{u_k}}}{{{\text{d}}{s^4}}} + \frac{{\text{d}}}{{{\text{d}}s}}\left ( {\frac{F}{{EI}}\frac{{{\text{d}}{u_k}}}{{{\text{d}}s}}} \right) = {K_{\text{B}}}{q_k} \quad\quad (k = 1,2) $$ (1) 式中,k=1代表井斜平面,k=2代表方位平面;
${u_k}$ 为管柱在井斜平面或方位平面上的横向位移,m;s为管柱任意一点的弧长,m;EI为相应的抗弯刚度,N · m2;F为轴向力,N;KB为管柱的浮力系数,KB=1−ρd/ρs;ρd为钻井液密度,g/cm3;ρs为管柱的材料密度,g/cm3;qk为管柱线重投影到井斜平面或方位平面的分量,N/m;由于井斜平面与方位平面相互垂直,管柱线重qk在方位平面上的投影始终为0,即q2为0。本文利用加权余量法求解了导向钻具组合的变形问题,需要构造BHA的变形试挠度函数作为控制方程的近似解,即可令BHA第j跨的试挠度函数为
$$ {u_{k,j}}\left ( s \right) = \sum\limits_{i = 1}^5 {{C_{k,i,j}}{s^{5 - i}}\quad\quad {\left ( {j = {\text{1,2}}, \cdots ,n{\text{ + 1}}} \right)} } $$ (2) 式中,Ck,i, j为井斜平面和方位平面的待定系数,无因次;n表示导向钻具单元节点个数,节点包括BHA的结构弯角、稳定器、管柱与井壁的接触点等。
假定底部钻具的下部边界为铰支约束,钻具上端躺在井眼底边,即
$$ \left\{ \begin{gathered} {u_{k,n + 1}}({L_{k,n + 1}}) = ({1 \mathord{\left/ {\vphantom {1 {{\kappa _k}}}} \right. } {{\kappa _k}}} + {\lambda _k} \cdot {r_{n + 1}})\cos {\psi _{k,n + 1}} - {1 \mathord{\left/ {\vphantom {1 {{\kappa _k}}}} \right. } {{\kappa _k}}} \\ {{u'}_{k,n + 1}}({L_{k,n + 1}}) = - (1 + {\lambda _k} \cdot {r_{n + 1}} \cdot {\kappa _k})\sin {\psi _{k,n + 1}} \\ {{u''}_{k,n + 1}}({L_{k,n + 1}}) = - (1 + {\lambda _k} \cdot {r_{n + 1}} \cdot {\kappa _k})\cos {\psi _{k,n + 1}} \cdot {\kappa _k} \end{gathered} \right. $$ (3) 式中,
$ {u'_k} $ 和$ {u''_k} $ 分别表示$ {u_k} $ 的1阶导数和2阶导数;Lk,n+1表示钻具上切段的长度在井斜平面和方位平面的投影,m;λk为无因次系数(如果将Lk,n+1投影至井斜平面,λ1=1;如果投影至方位平面,需要考虑投影所产生的居中效应,此时λ2=0);rn+1表示上切点处的井眼视半径,m;ψk,n+1表示钻头至钻具上切点的狗腿角在井斜平面和方位平面上的投影,rad;κk表示井眼曲率在井斜和方位平面上的分量,rad/30 m。将式(2)代入到边界条件式(3)中可得到关于Ck,i, j的线性方程,将式(2)代入到控制方程式(1)并利用子域法可得到关于Ck,i, j的线性方程。由于钻具上切段长度是待定的,因此需要采用迭代法求解所有的未知量。
钻头处(j=1)的轴向力F1和支反力Nsk,1可表示为
$$ \left\{ \begin{gathered} {F_1} = {W_{{\text{ob}}}} \\ {N_{{\text{s}}k,1}} = 6E{I_1}{C_{k,2,1}} + {W_{{\text{ob}}}}{C_{k,4,1}} \\ \end{gathered} \right. $$ (4) 其余各跨节点(j=2,···,n+1)的轴向力Fj和支反力Nsk,j表示为
$$ \left\{ \begin{gathered} {F_j} = {F_{j - 1}} + {K_{\rm{B}}} \cdot \sum\limits_{j = 2}^{N + 1} {{q_{{\text{1}},j - 1}} \cdot {L_{j - 1}} \cdot \cos {\beta _{j - 1}}} \\ {N_{{\text{s}}k,j}} = - E{I_j}(24{L_{k,j}}{C_{k,1,j}} + 6{C_{k,2,j}}) + 6E{I_{j{\text{ + 1}}}}{C_{k,2,j + 1}} \\ \end{gathered} \right. $$ (5) 式中,Wob为钻压,一般无论对于井斜平面还是方位平面,钻头处的轴向力F1等于钻压,N; Nsk,1为钻头处的支反力,可表示为钻头处井斜平面和方位平面的侧向力,N;L j−1为管柱第j−1跨的弧长,m;Lk, j为管柱在井斜平面和方位平面第j跨的弧长,m;βj−1为管柱第j−1跨的井斜角,rad。
-
在定向造斜过程中,一般采用滑动钻进方式进行导向钻进。然而在实际工况中,接触点产生较高的摩阻会使得BHA难以连续滑动,甚至会大大降低管柱的轴向力传递效率,因此需要合理计算滑动钻进时管柱的摩阻分布。依据式(4)、式(5),可得到三维井眼中各跨段管柱的轴向载荷和摩阻计算公式。
$$ {F_{{\text{r}},j}} = {F_j} \pm {f_j} $$ (6) $$ {f_j} = \mu \sqrt {N_{{\text{s}}1,j}^2 + N_{{\text{s}}2,j}^2} $$ (7) 式中,Fr, j为第j跨管柱的实际轴向载荷,N;fj为第j跨管柱的摩阻,N;μ为管柱的轴向摩阻系数。
-
在钻头与地层相互作用的影响下,钻头的位移方向并非与其机械合力方向一致。为了得到更加精确的造斜率计算模型,本文考虑BHA实际钻进过程的影响,即侧向钻速为0,建立滑动钻进造斜率的计算方法。
根据式(4),进一步考虑钻头与地层相互作用[18],建立井斜平面和方位平面上的钻速方程。
$$ \left[ {\begin{array}{*{20}{c}} {{R_{\text{a}}}} \\ {{R_k}} \end{array}} \right]{\text{ = }}{D_{\text{n}}}{{\boldsymbol{R}}^{\text{T}}}\left[ {\begin{array}{*{20}{c}} 1&{} \\ {}&{{I_{{\text{r}}k}}} \end{array}} \right]{\boldsymbol{R}}{{\boldsymbol{S}}^{\text{T}}}\left[ {\begin{array}{*{20}{c}} 1&{} \\ {}&{{I_{\text{b}}}} \end{array}} \right]{\boldsymbol{S}}\left[ {\begin{array}{*{20}{c}} {{F_1}} \\ {{N_{{\text{s}}k,1}}} \end{array}} \right] $$ (8) 式中,Ra、Rk分别为沿井眼轴线方向、井斜以及方位平面侧向方向上的钻速分量,m/h;Dn为标准(法向)钻井效率;R为井底坐标与地层坐标的转换矩阵;S为井底坐标与钻头坐标的转换矩阵;Irk为地层各向异性指数;Ib为钻头的各向异性指数。
由于式(3)中的井眼曲率κk并未得到解决,因此第1次迭代所得到钻头在井斜平面和方位平面处的侧向力Nsk,1中含有井眼曲率κk。为了确定式(3)中的井眼曲率κk,需要在第1次迭代的基础上,基于零侧向钻速准则,对式(8)再次进行牛顿迭代求解,可构造目标函数为
$$ {R_k}\left ( {{\kappa _k}} \right) \to 0 $$ (9) 考虑滑动钻进过程中工具面的影响,结构弯角的存在会导致BHA的作用方向线与管柱轴线构成一个工具面。在滑动钻进时,由于钻头侧向力的变化,导致工具造斜能力无法用某一特定的工具面去进行表述,因此需要通过调整装置角从0°到 360°的变化范围来表征结构弯角的方位,如图1所示。
图1中x、y、z轴相互正交并构成了井斜α平面和方位φ平面,交点为钻头;N表示正北方向;Φc为装置方位角,rad;Azc为导向钻具投影在水平面上的方位角,rad;ω为高边模式中,工具面与井底平面的交线所转过的装置角,rad;θ为导向钻具的结构弯角值,rad。
根据图1,结构弯角的方位表示为
$$ {\overline \theta _k} = \left\{ \begin{gathered} \theta \sin \omega \begin{array}{*{20}{c}} {}&{(k = 1)} \end{array} \\ \theta \cos \omega \begin{array}{*{20}{c}} {}&{(k = 2)} \end{array} \\ \end{gathered} \right. $$ (10) 式中,
$ {\overline \theta _k} $ 为结构弯角在井斜平面和方位平面的投影值,rad。根据式(9),利用牛顿迭代法进行求解,并引入式(10)中结构弯角方位变化的影响,将侧向钻速为0时井眼曲率作为滑动钻进的造斜率,计算公式为
$$ {\kappa _{\text{L}}} = \sqrt {\kappa _{{\text{1}},\omega }^{\text{2}}{\text{ + }}\kappa _{{\text{2}},\omega }^{\text{2}}} $$ (11) 式中,κL为滑动钻进在某一工具面下的理论造斜率,κ1,ω和κ2,ω分别为BHA在某一工具面对应的造斜率在井斜平面和方位平面上的投影,rad/30 m。
-
根据钻速方程式(8)可知,钻头各向异性指数Ib和地层各向异性指数Irk的确定是准确预测工具造斜率的关键所在。一般地,Ib和Irk可以通过实钻资料进行反演计算。当Ib=1时,钻头为各向同性;当0<Ib<1时,钻头的轴向切削能力强,并且当Ib越小时,钻头的轴向切削能力越强;当Ib>1时,钻头的侧向切削能力强。然而,即使在均质的各向同性地层条件下,钻头实际的运动轨迹也并非沿着机械合力的方向。对于Irk而言,在考虑地层层面各向差异和Ib的影响下,会使得Irk确定更加困难。因此,本文假设地层为各向同性的均质地层(Irk=1),进一步简化了式(8)计算的复杂性,并引入折算系数对式(11)进行修正,得到实际造斜率的计算公式为
$$ {\kappa _{\text{A}}} = \lambda {\kappa _{\text{L}}} $$ (12) 式中,
$ {\kappa _{A} } $ 为滑动钻进时真实造斜率,rad/30 m;λ为折算系数,无因次。引入折算系数,可将地层各向异性、实钻过程等参数引入到造斜率预测模型中,λ的取值范围大致为0.3~0.9。
Prediction and analysis of build-up rate during sliding drilling
-
摘要: 滑动钻进工况下造斜率的准确预测是钻井工程面临的难点问题。首先利用加权余量法,建立了导向钻具组合的三维力学模型,得到了滑动钻进过程下钻头机械载荷及导向钻具组合的摩阻分布;其次考虑滑动钻进过程中工具的造斜特性,建立了零侧向钻速准则下的滑动钻进造斜率预测模型及折算系数反演方法;最后利用一口页岩气水平井钻井数据,反演计算了该井造斜率折算系数,验证了滑动钻进造斜率预测模型的精度,并分析了钻井参数、轨迹参数、钻头切削特性、工具结构等对工具造斜率和摩阻的影响规律。研究结果表明:滑动导向模式下,通过适当增加工具结构弯角、降低钻压、提高钻头侧向切削能力等,有利于发挥工具的造斜性能,提高滑动钻进效率,为导向钻井参数优化、工具结构优化设计及井眼轨迹控制等提供重要依据。Abstract: The drilling engineering is facing difficulties in accurately predicting the build-up rate during the sliding drilling operation. First, this research developed the three-dimensional mechanical model of the steering bottomhole assembly (BHA), based on the weighted residual method, and calculated the mechanical load of the bit and friction distribution in the steering BHA during sliding. Moreover, given the build-up characteristics of the tool during sliding, the sliding build-up rate prediction model with a zero side penetration rate and the inversion method of the conversion coefficient were developed. At last, the build-up rate conversion coefficient of a shale gas horizontal well was inverted from its drilling data and this application confirms the precision of the proposed sliding build-up rate prediction model. Besides, the effects of drilling parameters, trajectory parameters, cutting characteristics of the bit, and tool structure on the build-up rate and friction of the tool are investigated. It is shown that during the sliding steering, it is in favor of improving the build-up performance of the tool and the sliding drilling efficiency to properly increase the structural bend angle of the tool, lower the weight on bit, and enhance the side cutting performance of the bit. The findings of this research provide vital references for optimizing the steering drilling parameters, tool structure design, and well trajectory control.
-
Key words:
- sliding drilling /
- steering drilling tool /
- build-up rate /
- friction /
- structural bend angle /
- trajectory control
-
-
[1] 高德利. 复杂结构井优化设计与钻完井控制技术[M]. 东营: 中国石油大学出版社, 2011. GAO Deli. Optimized design and control techniques for drilling & completion of complex-structure wells[M]. Dongying: China University of Petroleum Press, 2011. [2] 张鹏程. 钻井工程井斜产生原因及其控制措施[J]. 西部探矿工程, 2019, 31(7):69-70. doi: 10.3969/j.issn.1004-5716.2019.07.026 ZHANG Pengcheng. Causes and control measures of well deviation in drillingengineering[J]. West-China Exploration Engineering, 2019, 31(7): 69-70. doi: 10.3969/j.issn.1004-5716.2019.07.026 [3] 高德利, 狄勤丰, 张武辇. 南海西江大位移井定向控制技术研究[J]. 石油钻采工艺, 2004, 26(2):1-4. doi: 10.3969/j.issn.1000-7393.2004.02.001 GAO Deli, DI Qinfeng, ZHANG Wunian. Study on directioanl control technology for Xijiang extended - reach wells in South China sea[J]. Oil Drilling & Production Technology, 2004, 26(2): 1-4. doi: 10.3969/j.issn.1000-7393.2004.02.001 [4] 高德利. 典型导向钻具组合的力学分析方法[J]. 石油大学学报(自然科学版), 1993, 17(6):35-40. GAO Deli. A method for modeling static behavior of downhole motor assembly with bent housing and stabilizers[J]. Journal of China University of Petroleum (Edition of Natural Science), 1993, 17(6): 35-40. [5] 狄勤丰. 滑动式导向钻具组合复合钻井时导向力计算分析[J]. 石油钻采工艺, 2000, 22(1):14-16. doi: 10.3969/j.issn.1000-7393.2000.01.004 DI Qinfeng. Calculation and analysis of the steering force of the compounding slide steering assembly[J]. Oil Drilling & Production Technology, 2000, 22(1): 14-16. doi: 10.3969/j.issn.1000-7393.2000.01.004 [6] 秦永和, 付胜利, 高德利. 大位移井摩阻扭矩力学分析新模型[J]. 天然气工业, 2006, 26(11):77-79. doi: 10.3321/j.issn:1000-0976.2006.11.023 QIN Yonghe, FU Shengli, GAO Deli. A new model for analyzing torque and drag in extended reach wells[J]. Natural Gas Industry, 2006, 26(11): 77-79. doi: 10.3321/j.issn:1000-0976.2006.11.023 [7] 祝效华, 刘清友, 李红岩, 等. 井眼轨迹预测理论及方法研究[J]. 天然气工业, 2004, 24(4):38-40. doi: 10.3321/j.issn:1000-0976.2004.04.014 ZHU Xiaohua, LIU Qingyou, LI Hongyan, et al. Study on theory and method of borehole trajectory prediction[J]. Natural Gas Industry, 2004, 24(4): 38-40. doi: 10.3321/j.issn:1000-0976.2004.04.014 [8] BISCARO E, D’ALESSANDRO J D, MORENO A, et al. New rotary steerable drilling system delivers extensive formation evaluation for high build rate wells[C]//Paper presented at the SPE Western Regional Meeting, April 2015, Garden Grove, California, USA: SPE-174016-MS. [9] AL DABYAH A A, ZUBAIRI B A, AL UMAIRIN I S, et al. Proper planning supported by latest technology leads to drilling the longest 8.5-in. Section in an extended-reach well: Case study[C]//Paper presented at the SPE Annual Technical Conference and Exhibition, September 2016, Dubai, UAE: SPE-181363-MS. [10] KENUPP R, LOURENCO A, SOLTVEDT S, et al. The longest horizontal section ever drilled in an extended-reach well in Brazil[C]//Paper presented at the OTC Brasil, October 2017, Rio de Janeiro, Brazil: OTC-28091-MS. [11] HASSEN B R, MACDONALD A J. Field comparison of medium and long radius horizontal wells drilled in the same reservior[C]//Paper presented at the IADC/SPE Drilling Conference, February 1990, Houston, Texas: SPE-19986-MS. [12] BRADLEY W B. Factors affecting the control of borehole angle in straight and directional wells[J]. Journal of Petroleum Technology, 1975, 27(6): 679-688. doi: 10.2118/5070-PA [13] 王恒, 管志川, 史玉才, 等. 柔性短节对推靠式旋转导向底部钻具组合造斜能力的影响分析[J]. 钻采工艺, 2018, 41(6):19-22. doi: 10.3969/J.ISSN.1006-768X.2018.06.06 WANG Heng, GUAN Zhichuan, SHI Yucai, et al. Effect of flex sub on build-up performance of push the bit RSBHA[J]. Drilling & Production Technology, 2018, 41(6): 19-22. doi: 10.3969/J.ISSN.1006-768X.2018.06.06 [14] 夏焱, 管志川, 史玉才. 底部钻具组合参数设计的评价方法[J]. 石油大学学报(自然科学版), 2005, 29(4):48-51,63. doi: 10.3321/j.issn:1000-5870.2005.04.011 XIA Yan, GUAN Zhichuan, SHI Yucai. Evaluation method of parameters design for Bottom hole assembly[J]. Journal of China University of Petroleum (Edition of Natural Science), 2005, 29(4): 48-51,63. doi: 10.3321/j.issn:1000-5870.2005.04.011 [15] 苏义脑. 极限曲率法及其应用[J]. 石油学报, 1997, 18(3):112-116. doi: 10.7623/syxb199703017 SU Yi’nao. A method of limiting curvature and its applications[J]. Acta Petrolei Sinica, 1997, 18(3): 112-116. doi: 10.7623/syxb199703017 [16] 狄勤丰, 彭国荣. 下部钻具组合造斜能力的预测方法研究[J]. 石油钻采工艺, 2000, 22(2):17-19. doi: 10.3969/j.issn.1000-7393.2000.02.004 DI Qinfeng, PENG Guorong. A new method for predicting the deflecting capacity of Bottom hole assembly[J]. Oil Drilling & Production Technology, 2000, 22(2): 17-19. doi: 10.3969/j.issn.1000-7393.2000.02.004 [17] 杨君, 何元君, 郝荣明. 适应单弯滑动导向钻进的新型动力钻具优化设计[J]. 石油钻探技术, 2003, 31(3):40-42. doi: 10.3969/j.issn.1001-0890.2003.03.016 YANG Jun, HE Yuanjun, HAO Rongming. Optimum design of new style mud motor applied to steering drilling[J]. Petroleum Drilling Techniques, 2003, 31(3): 40-42. doi: 10.3969/j.issn.1001-0890.2003.03.016 [18] 高德利. 井眼轨迹控制[M]. 东营: 石油大学出版社, 1994. GAO Deli. Hole trajectory control[M]. Dongying: Petroleum University Press, 1994 [19] 郑德帅, 张华卫, 朱晔, 等. 岩石力学性质对井眼轨迹走向的影响机理[J]. 石油钻探技术, 2014, 42(3):45-49. doi: 10.3969/j.issn.10010890.2014.03.009 ZHENG Deshuai, ZHANG Huawei, ZHU Ye, et al. Affecting mechanism of rock mechanical property on well trajectory[J]. Petroleum Drilling Techniques, 2014, 42(3): 45-49. doi: 10.3969/j.issn.10010890.2014.03.009 [20] 王家骏, 邹德永, 姜鑫, 等. PDC钻头各向异性指数的计算方法研究[J]. 石油机械, 2014, 42(9):27-30. doi: 10.3969/j.issn.1001-4578.2014.09.007 WANG Jiajun, ZOU Deyong, JIANG Xin, et al. Calculation method of anisotropy index of PDC bit[J]. China Petroleum Machinery, 2014, 42(9): 27-30. doi: 10.3969/j.issn.1001-4578.2014.09.007 [21] 陈军海, 王怡, 韩艳浓, 等. 一种利用岩石强度刻划试验确定岩石PDC钻头可钻性级值的方法[J]. 中国石油大学学报(自然科学版), 2017, 41(2):94-99. doi: 10.3969/j.issn.1673-5005.2017.02.011 CHEN Junhai, WANG Yi, HAN Yannong, et al. A new method for assessing rock’s drillabity for PDC bits via a scratch testing[J]. Journal of China University of Petroleum (Edition of Natural Science), 2017, 41(2): 94-99. doi: 10.3969/j.issn.1673-5005.2017.02.011