白田增,任勇强,杨春华,吴 德,郑 博,傅小勇
渤海钻探工程有限公司井下作业分公司,河北 任丘
连续油管具有成本低廉、作业方便,可靠性高等特点,被广泛应用于油气田勘探开发作业各个领域。在油气开发作业中,井下液体内外压及油管失稳弯曲等因素严重影响连续油管作业的安全性与可靠性[1][2]。笔者基于实际连续油管相关研究,对连续油管下入作业中的力学特性进行研究。
利用任意井深坐标参数建立井眼轴线轨迹几何模型并绘制井眼轨迹和计算连续油管下入过程的载荷,通过数值积分方法[3][4],求得各点坐标的坐标参数,得出井眼轴坐标参数。
已知任意井的井深坐标参数(井深、井斜、方位角),绘制其3维井眼轨迹。将对应井深坐标参数的3维井眼轨迹曲线进行拟合,对整体井眼轨迹拟合曲线进行等步长划分,共划分出n个节点,对节点依次编号,取其中第i个节点作为分析对象,其下一个节点编号对应为i+1。基于3维井眼轨迹曲线建立直角坐标系(图1),坐标原点为井眼轨迹起点即井口位置。取井眼轨迹上任意微元段为研究对象,长度为dL,与z轴夹角为井斜角1,该微元段在xoy平面的投影长度记为dH,其与y轴夹角为方位角φ。
根据几何关系,得到如下关系式:
Figure 1.The micro-elemental model of well trajectory图1.井眼轨迹微元模型
由于坐标原点为井眼轨迹起点,因此井口位置边界条件为:x0=0,y0=0,z0=0,H0=0。任意井深L处的坐标参数关系式为:
基于上述关系式,对于已知井斜角、方位角等井眼轨迹参数,可以得出任意井的任意点坐标参数,为后续连续油管下入分析建立数值基础。
连续油管的压力计算可用于连续油管的强度校核和管内流体计算,也可判断内外压力对连续油管强度的影响,指导连续油管安全性能评估。依照上述节点的坐标参数,对于任意连续油管节点i处内外压的差分方程其关系式为[5]:
式中:p1、p2为应节点处连续油管的内、外流体的压力,Pa;ρ1、ρ2为对应节点处内、外流体的密度,kg/m3;υ1、υ2为对应节点处内、外流体的流速,m/s;λ1,i、λ2,i为相应节点处内、外流体的无因次黏滞摩阻系数,1;Δli为连续油管第i单元的长度Δli=li+1-li,m;r1为连续油管内径,mm;r2为油管等效环空水力半径,m。
通过求解连续油管轴向力,可以计算出任意节点处的无因次轴力,得到连续油管在任意节点处的无因次正压力,用以判断油管在下入过程中是否会发生屈曲。
基于微元体物理方程,可以得到任意节点处的等效轴向力方程[6]:
式中:Feτ为对应节点处的等效轴向力,N;qe为连续油管节点单元的单位长度重量,N;fi为节点单元与井壁或套管壁之间的库仑摩擦系数,1;f1、f2为对应节点的内、外壁黏滞摩阻力,N/m。
真实轴力的计算公式为:
式中:Fτ为对应节点处的真实轴力,N;A1、A2分别为连续油管任意节点处的内圆和外圆的截面积,mm2。
通过计算连续油管任意一点弯矩,可以判断连续油管是否出现弯矩过大,指导连续油管安全评估。
计算任意节点处弯矩的关系式为[7]:
式中:k0,i为井眼轨迹对应节点曲率,1/m;E为连续油管对应材料的弹性模量,Pa;I为截面的惯性矩,mm4。ir为轴线与套管内壁轴线的距离,mm;b为无因次轴力,1;M为弯距,N∙mm。
1)根据钻机作业下的井深坐标参数,建立井眼轨迹几何模型,对井斜角、方位角函数简单的情况直接积分计算;当函数较为复杂时,可以采用复化辛普森求积公式的数值积分方法计算井眼轴坐标参数。
2)对于已知井眼轴坐标参数的情况下,由环空压力、轴向力、弯矩等力学特性计算公式,可以求得任意井对应节点处的力学特性,为连续油管作业提供力学分析基础。(2018ZD21Y-03)。