考虑钻头转速影响的新三维钻速方程
2015-04-07刘军波韦红术赵景芳张日鹏
刘军波， 韦红术， 赵景芳， 张 辉， 张日鹏
(1.中海油田服务股份有限公司钻井事业部，河北三河 065201；2.中海石油(中国)有限公司深圳分公司，广东深圳 518067；3.中国石油大学(北京)石油工程学院，北京 102249)
考虑钻头转速影响的新三维钻速方程
刘军波1， 韦红术2， 赵景芳1， 张 辉3， 张日鹏3
(1.中海油田服务股份有限公司钻井事业部，河北三河 065201；2.中海石油(中国)有限公司深圳分公司，广东深圳 518067；3.中国石油大学(北京)石油工程学院，北京 102249)
为了更全面地考虑影响机械钻速的因素、提高预测机械钻速和井眼轨迹的准确性，深入研究了涉及钻头和地层相互作用的三维钻速方程。考虑钻头转速的影响，重新定义了钻头各向异性指数和地层各向异性指数；利用地层各向异性指数和钻头各向异性指数的定义以及坐标变换，建立了考虑钻头转速影响、钻头与地层相互作用的新三维钻速方程。利用新三维钻速方程，实例分析了钻头转速对机械钻速的影响，随着钻头转速的增加，井底平面Xd方向、Yd方向及井底法向Zd方向的机械钻速都随之提高，但井底3个方向的机械钻速随钻头转速增加而提高的幅度并不相同，井底法向Zd方向机械钻速的提高幅度要大于井底Xd和Yd方向机械钻速的提高幅度。研究表明，钻井过程中，利用新三维钻速模型可以预测钻头转速变化对机械钻速和井眼轨迹的影响，也可以通过调整钻头转速、钻头各向异性指数、钻头机械力等参数对井眼轨迹进行控制。
钻头转速 钻速方程 各向异性指数
对钻头和地层相互作用情况的研究已有二三十年的历史，许多学者都建立了相关模型。其中，白家祉等人[1]建立的“多元幂积侧向切削模型”考虑的因素比较全面，但涉及到大量参数反演，需要大量的现场资料或实验资料，因而现场应用价值不高；高德利、李子丰和闫铁等人[2-4]的三维钻速方程都只是建立了力与机械钻速的关系，没有考虑钻头转速对机械钻速以及转速对钻头和地层各向异性指数的影响。而按文献[5]对钻头各向异性指数的定义，不同钻头转速下测量出来的钻头各向异性指数是不同的；同理可知，按以前的地层各向异性指数定义，不同钻头转速下测量出来的地层各向异性指数也是不同的。但是，钻头各向异性指数与地层各向异性指数是钻头和地层本身的性质参数，在自身没有发生改变的情况下，不会因为外界测量参数的改变而改变。为此，笔者重新定义了钻头各向异性指数和地层各向异性指数，建立了考虑钻头转速影响的新三维钻速方程，并利用该方程进行了实例分析，分析了其实用价值。
1 各向异性指数的定义
1.1 地层的各向异性
地层的各向异性程度可以用地层各向异性指数来表征。假设钻头为各向同性，所钻地层为正交各向异性地层，定义钻井效率和地层各向异性指数为：
(1)
(2)
式中：Ddip，Dstr和Dn分别为地层倾向、走向及法向的钻井效率；vdip，vstr及vn分别为钻头沿地层倾向、走向及法向的机械钻速，m/h；Fdip，Fstr和Fn分别为钻头沿地层倾向、走向及法向的净作用力，kN；n为钻头转速，r/min；λ1，λ2和λ3分别为地层倾向、走向和法向的转速指数；Ir1和Ir2分别为地层倾向和走向的各向异性指数。
1.2 钻头的各向异性
钻头的各向异性是指钻头在不同方向上钻进能力的差异。一般而言，无论是牙轮钻头还是PDC钻头，其结构都有利于向前(即沿轴向)钻进，钻头沿轴向的钻进效率优于沿侧向的钻进效率，即存在钻进效率的各向异性。钻头的各向异性可用钻头各向异性指数Ib来表征。假设各向异性钻头在各向同性地层里钻进，可定义钻井效率和钻头各向异性指数为：
(3)
式中：Da和Dl分别为钻头轴向和侧向的钻井效率；Fa和Fl分别为钻头轴向和侧向的净作用力，kN；va和vl分别为钻头轴向和侧向的机械钻速，m/h；λ4和λ5分别为钻头侧向和轴向的转速指数；Ib为钻头的各向异性指数[6-9]。
2 考虑钻头转速新三维钻速方程的建立
2.1 坐标系的建立
建立了大地坐标系、地层坐标系、井底坐标系、钻头坐标系等4个坐标系，均取钻头处为坐标原点O。
2.1.1 大地坐标系(O-XYZ)
建立的大地坐标系(O—XYZ)中，X在水平面内并指向正东方向，Y在水平面内并指向正北方向，Z垂直于XOY平面并指向上方，如图1所示。
2.1.2 地层坐标系(O-XfYfZf)
建立的地层坐标系(O—XfYfZf)中，Xf沿地层层面指向上倾方向，Yf沿地层走向并符合右手法则，Zf垂直于地层层面并指向外法线方向，如图2所示。图2中，β为地层倾角，(°)。
2.1.3 井底坐标系(O-XdYdZd)
建立的井底坐标系(O-XdYdZd)中，Xd指向井眼低边，Zd垂直于井底平面并指向外法线方向，Yd垂直于XdOZd平面并符合右手法则，如图3所示。图3中，α为井斜角，(°)。
2.1.4 钻头坐标系(O-XbYbZb)
建立的钻头坐标系(O-XbYbZb)中，Zb沿钻柱变形后的切线方向指向上方，Xb轴垂直于Zb轴指向井眼低边，Yb轴垂直于XbOZb平面并符合右手法则，如图4所示。
2.2 坐标系之间的转换关系
上述4个坐标系是可以相互转换的。由地层坐标系到大地坐标系的转换关系为：
(4)
(5)
式中：φf为地层下倾方位角，(°)。
由井底坐标系到大地坐标系的变换关系为：
(6)
(7)
式中:φ为井斜方位角，(°)。
由于地层坐标系和井底坐标系是不断变化的，直接求地层坐标系到井底坐标系的转换矩阵比较困难，因此要先求出地层坐标系和井底坐标系到恒定大地坐标系的转换矩阵，再据此求得到地层坐标系到井底坐标系的变换关系。由地层坐标系到井底坐标系的变换关系为：
(8)
其中D=B-1A
(9)
(10)
(11)
Δφ=φ-φf
(12)
由钻头坐标系到井底坐标系的变换关系为：
(13)
(14)
式中：θx，θy分别为钻头转角在XdOZd平面和YdOZd平面上的投影角，rad。
钻头坐标系和井底坐标系的相对位置关系如图5所示。θx符号约定：正对Yd正方向，如果Zd轴正向旋转到Zb轴正向所转过的角度为逆时针，则θx取正，反之取负。θy符号约定：正对Xd正方向，如果Yd轴正向旋转到Yb轴正向所转过的角度为逆时针，则θy取正，反之取负。
2.3 方程的推导
设vxd，vyd和vzd分别为井底平面上钻头机械钻速沿Xd，Yd和Zd方向的分量，Fxd，Fyd和Fzd分别为井底平面上钻头机械力沿Xd，Yd和Zd方向的分量，由地层各向异性和钻井效率的定义，它们之间的关系可表示为：
(15)
(16)
(17)
净作用力Fdip，Fstr和Fn可以通过井底平面上钻头机械力Fxd，Fyd和Fzd转换而来。由于得出式(15)的前提是钻头为各向同性，因此，需要对钻头的机械力进行等效化处理：首先，将钻头机械力转换到钻头坐标系中；然后，假设钻头为各向同性，即如果钻头侧向和轴向的受力大小相等，那么钻头在侧向和轴向的机械钻速也相同，进一步可以表述为“此时钻头的侧向和轴向切削能力均为Dnnλ5”；在钻头各向同性的假设下，钻头机械力的等效力在钻头坐标系下可以表示为：
(18)
式中：F1,F2和F3分别为钻头各向同性条件下钻头机械力的等效力在Xb,Yb和Zb方向的分量，kN。
将钻头机械力的等效力从钻头坐标系转换到地层坐标系，转换式为：
(20)
(21)
(22)
(23)
将式(19)代入式(15)，得：
(24)
将地层坐标系中的机械钻速转换到井底坐标系中，可得：
(25)
式(25)即为考虑钻头转速影响的新三维钻速方程。
针对某种钻头，式(25)中的Ib，λ4和λ5可以通过试验测得，具体方法可以参考文献[5]和修正的杨格钻速方程之五点法钻速试验求解系数，也可以通过三维钻速方程与BHA力学分析反演求取。Ir1,Ir2及λ1,λ2,λ3可以通过岩心地面试验测得[8]，也可以通过该三维钻速方程与BHA力学分析反演求取；Ir1,Ir2还可以通过分析随钻声波测井资料求取[10-12]。
3 钻头转速影响机械钻速实例分析
由式(25)可得，钻头转速、转速指数、钻头机械力、钻头各向异性指数、地层各向异性指数等参数的变化都对机械钻速有影响。钻头机械力、钻头各向异性指数和地层各向异性指数对机械钻速的影响在以前的三维钻速方程中已经分析过，这里仅分析钻头转速对机械钻速的影响规律。
由式(25)可知机械钻速是由钻头转速的幂函数参与组成的多项式，为了定量描述机械钻速随钻头转速的变化规律，选取实例如下：井斜角α=0°，地层倾角β=30°，井斜方位角和地层下倾方位角之差Δφ=30°。则地层坐标系与井底坐标系的变换矩阵D为：
(26)
其逆矩阵为：
(27)
为减少计算量，对矩阵F和F-1进行简化。一般情况下，钻头转角小于0.4°可视为小量[1]，则可以作如下近似处理：
(28)
其中 θx,θy∈[0,0.007]rad
(29)
将式(28)带入式(14)，并忽略二阶小量，则钻头坐标系与井底坐标系的变换矩阵F为：
(30)
其逆矩阵为：
(31)
取钻头转角θx=0.003 rad,θy=0.005 rad，则矩阵F为：
(32)
其逆矩阵为：
(33)
一般而言，钻头的轴向切削能力高于侧向，沿着地层法线方向钻进要比沿着地层倾向或地层走向钻进容易，所以：0
从图6可以看出，随着钻头转速从40 r/min增至180 r/min，井底平面Xd方向的机械钻速从2.3 m/h提高至6.1 m/h，井底平面Yd方向的机械钻速从4.2 m/h提高至12.5 m/h，井底法向Zd方向的机械钻速从8.4 m/h提高至23.4 m/h。通过实例分析可知：井底3个方向的机械钻速随钻头转速增加而提高的幅度并不相同，井底法向Zd方向机械钻速的提高幅度要大于井底Xd和Yd方向机械钻速的提高幅度。
4 结论与认识
1) 建立的新三维钻速方程考虑了钻头转速对机械钻速的影响，这样更符合现场实际情况。对钻头和地层的各向异性指数重新进行了定义，考虑钻头转速因素，避免了因钻头转速改变而导致测出的各向异性指数不同。
2) 井底3个方向的机械钻速随钻头转速增加而提高的幅度并不相同，井底法向Zd方向机械钻速的提高幅度要大于井底Xd和Yd方向机械钻速的提高幅度。
3) 新三维钻速方程与BHA力学分析相结合，作为正演模型可以预测钻头转速变化对机械钻速和井眼轨迹的影响，进一步可以通过调节钻头转速、钻头机械力、钻头各向异性指数等参数来控制井眼轨迹；作为反演模型可以进行钻头和地层各向异性参数反演。
References
[1] 白家祉，苏义脑.井斜控制理论与实践[M].北京：石油工业出版社，1990：3-25. Bai Jiazhi,Su Yinao.Theory and practice of deviation control[M].Beijing:Petroleum Industry Press,1990:3-25.
[2] 高德利，刘希圣，黄荣樽.钻头与地层相互作用的三维宏观分析[J].石油大学学报:自然科学版，1989，13(1)：23-31. Gao Deli,Liu Xisheng,Huang Rongzun.Three dimensional macro analysis of rock bit interaction[J].Journal of the University of Petroleum，China：Edition of Natural Science,1989,13(1):23-31.
[3] 闫铁，张建群.钻头与地层相互作用的三维理论分析[J].天然气工业，1991，11(6)：40-46. Yan Tie,Zhang Jianqun.Analysis of three dimensional theory of rock bit interaction[J].Natural Gas Industry,1991,11(6):40-46.
[4] 李子丰，李敬元.三维钻速方程及其反演模型[J].西部探矿工程，1994，6(2)：8-11.Li Zifeng,Li Jingyuan.Three dimensional drilling rate equation and the inversion model[J].West-China Exploration Engineering,1994,6(2):8-11.
[5] 高德利，聂翠平，刘希圣.钻头各向异性钻井特性的理论分析和试验[J].石油大学学报:自然科学版，1991，15(6)：28-31. Gao Deli,Nie Cuiping,Liu Xisheng.Theoretical analysis and test of bit anisotropic drilling characteristics[J].Journal of the University of Petroleum,China:Edition of Natural Science,1991,15(6):28-31.
[6] 高德利.钻头和地层各向异性钻井特性的一种表达方法[J].石油学报，1994，15(2)：126-132. Gao Deli.An expression method of bit and formation of anisotropic drilling characteristics[J].Acta Petrolei Sinica,1994,15(2):126-132.
[7] 高德利，刘希圣，徐秉业.井眼轨迹控制[M].东营：石油大学出版社，1994：130-190. Gao Deli,Liu Xisheng,Xu Bingye.Prediction and control of wellbore trajectory[M].Dongying:Petroleum University Press,1994:130-190.
[8] 高德利.地层各向异性的评估方法[J].石油学报，1993，14(2)：96-101. Gao Deli.Method for the evaluation of formation anisotropy[J].Acta Petrolei Sinica,1993,14(2):96-101.
[9] 高德利，刘希圣.典型地层的各向异性钻井特性[J].石油大学学报:自然科学版，1990，14(5)：1-8. Gao Deli,Liu Xisheng.Anisotropic drilling characteristics of typical formation[J].Journal of the University of Petroleum，China:Edition of Natural Science,1990,14(5):1-8.
[10] 潘起峰，高德利.用声波法评价地层可钻性各向异性的实验研究[J].岩石力学与工程学报，2006，25(1)：162-167. Pan Qifeng,Gao Deli.Experimental study on evaluation of formation drillability anisotropy by acoustic wave velocity[J].Chinese Journal of Rock Mechanics and Engineering,2006,25(1):162-167.
[11] 周祖辉，尹宏锦.岩石可钻性各向异性的微钻头试验研究[J].华东石油学院学报:自然科学版，1986，10(4)：22-24. Zhou Zuhui,Yin Hongjin.Study of micro drill test of rock drillability anisotropy[J].Journal of East China Petroleum Institute：Edition of Natural Science,1986,10(4):22-24.
[12] 王瑞甲，乔文孝，鞠晓东.利用随钻正交偶极子声波测井评价地层各向异性的数值研究[J].地球物理学报，2012，55(11)：3871-3873. Wang Ruijia,Qiao Wenxiao,Ju Xiaodong.Numerical study of formation anisotropy evaluation using cross dipole acoustic LWD[J].Chinese Journal of Geophysics,2012,55(11):3871-3873.
[编辑 令文学]
A New 3D ROP Equation Considering the Rotary Speed of Bit
Liu Junbo1， Wei Hongshu2, Zhao Jingfang1, Zhang Hui3， Zhang Ripeng3
(1.DrillingDivisionofChinaOilfieldServicesLimited,Sanhe,Hebei,065201,China;2.CNOOC(Shenzhen)Company,Shenzhen,Guangdong,518067,China;3.CollegeofPetroleumEngineering,ChinaUniversityofPetroleum(Beijing),Beijing,102249,China)
In order to fully consider the factors affecting the ROP and improve the accuracy of predicting the ROP and wellbore trajectory,the 3D ROP equation involving the interaction between bit and formation was studied. Considering the effect of rotary speed of the bit,the bit anisotropic index and the formation anisotropic index were redefined. Based on the definition of both anisotropic indexes and the coordinate transformation, a new 3D ROP equation considering the rotary speed of bit and the bit-formation interaction was established. Using the new 3D ROP equation,a practical case was taken to analyze the effect of rotary speed on ROP.It was found that all ROPs inXddirection,Yddirection and normalZddirection increased with the increase of rotary speed,but to different extents.ROP in normalZddirection increased more than that inXdandYddirections.The results showed that the new model could be used to predict the change of ROP and wellbore trajectory with the change of rotary speed in drilling operations,and that the wellbore trajectory could be controlled by adjusting the rotary speed,bit anisotropic index,bit mechanical force and other parameters.
rotary speed of the bit;ROP equation;anisotropic index
2014-01-17；改回日期：2014-11-30。
刘军波(1971—)，男，湖北武汉人，1994年毕业于石油大学(华东)钻井工程专业， 2011年获中国石油大学(华东)石油与天然气工程专业工程硕士学位，工程师，现主要从事钻井工艺与技术方面的研究工作。
国家自然科学基金项目“气体钻井井眼轨迹偏斜机理与控制理论研究”(编号：51174220)、中国石油科技创新基金项目“页岩气水平井钻井优化设计研究”(编号：2011D-5006-0310)资助。
◀钻井完井▶
10.11911/syztjs.201501009
TE21
A
1001-0890(2015)01-0052-06
联系方式：(022)66916393，liujb13@cosl.com.cn。
