1. 引言
井下油管柱是油气从井底顺利到达井口的唯一通道,保证油管柱在各种外载荷下的强度安全性是实现油气正常开采的基本条件之一。一体化油管柱下入井筒后可以完成酸压、生产等多种井下作业,避免管柱更换所需的多次上提、下放,节省了施工成本,提升了作业效率,特别适合深井和超深井,已广泛应用在高压气田的开发中 [1] 。由于一体化油管柱经历了不同的作业过程,承受的载荷变化复杂,其安全性取决于最恶劣工况,需要对一体化管柱从下入井筒到生产整个作业过程进行详细的力学分析,以获得管柱的最危险工况和危险位置,从而为管柱的合理设计和安全应用提供理论指导 [2] 。当前,许多学者对井下管柱的屈曲行为进行了大量研究,建立了预测管柱发生正弦或螺旋屈曲等永久变形的临界条件,可用于判断管柱的平衡状态,为稳定性分析提供依据 [3] 。因此,拟在准确描述井眼轨迹和管柱结构的基础上,建立考虑众多因素的管柱力学计算模型,分析一体化管柱在各工况下的受力和变形特点,提出油管柱的力学计算方法和封隔器作用力计算步骤,从强度和稳定性两方面综合评价一体化油管柱的安全性,为其设计和安全应用提供科学依据。
2. 力学模型的建立及求解
2.1. 三维弯曲井眼中管柱的几何描述
管柱在井筒中常因外力发生失稳,出现正弦或螺旋的屈曲变形 [4] ,导致其轴线与井眼轨迹不再重合或平行。如图1所示,在法平面
上,管柱轴线上C点相对于井眼轴线上的O点,有一偏离主法线
的偏转角
。假定弯曲失稳后管柱始终与井壁保持连续接触,则C点在以O点为圆心、半径为r的圆上,即:
(1)
或
(2)
式中:rw为井眼半径,m;ruo为管柱外半径,m;
为井眼轴线上O点的矢径,m;
为偏转角,(˚);
为主法线单位矢量,1;
为次法线单位矢量,1。
![](//html.hanspub.org/file/23061x19_hanspub.png)
Figure 1. The geometric description of downhole tubing string
图1. 井下管柱的几何描述
2.2. 管柱受力分析
管柱微元ds在井筒中的受力如图2所示。在单位长度上分布的外力矢量有: ① 管柱自重
,N/m (其中,q为管柱重量,N/m;
为竖直方向的单位矢量,1)。 ② 弯曲管柱与井壁之间的接触正压力
,N/m。 ③ 流体黏滞摩阻力
,N/m (其中,fui、fuo分别为管柱内、外的流体流动摩阻力,N/m;
为井眼轴线的切向单位矢量,1)。 ④ 管柱与井壁之间的轴向摩擦力
,N/m (其中,f1为管柱与井壁之间的轴向摩擦系数,1)。 ⑤ 环向摩擦力
,N/m (其中,f2为管柱与井壁之间的环向摩擦系数,1)。
管柱微元体内、外侧壁上的流体压力为pi(s)、po(s),Pa;其静力等效作用相当于分别在s和s + ds截面处作用有一对轴向压缩载荷(pi(s)Ai、po(s + ds)Ai)、拉伸载荷(po(s)Ao、po(s + ds)Ao)和一个分布载荷(
、
)。其中,Ai、Ao分别为管柱的内、外圆截面积,m2。
管柱微元在以上外力及内、外流体压力的作用下处于平衡状态,有平衡方程:
(3)
![](//html.hanspub.org/file/23061x30_hanspub.png)
Figure 2. The force diagram of micro-element
图2. 微元体受力示意图
化简,得:
(4)
引入等效力:
(5)
则
(6)
将单位长度上的外力代入式(6),并沿
投影,得管柱轴向力的平衡微分方程:
(7)
式中:
为单位管柱所受外力,N/m;
为等效力,N;
为等效力沿
的投影分量,N;
、
分别为微元两端面的内力矢量,N;fuw为井壁摩阻,N/m;Ah为环空面积,m2;
为井斜角,(˚);ρi和ρo分别为管柱内、外流体的密度,kg/m3。
2.3. 轴力计算
将总长为L的油管柱划分为n段,并从底端开始排列序号:0,1,…,n。根据方程(7),利用差分方法可以确定j+1节点处的轴力为:
(8)
j节点处的真实轴力为:
(9)
节点j处管柱受到的正压力按下式计算:
(10)
其中:
(11)
式中:
、
为第j个节点和第j + 1个节点的轴力,N;Δsj为管柱微元长度,m;qe,j为管柱的等效线重,N/m;ψ为方位角,(˚)。
系数nj随无量纲轴力βj的取值不同,井下管柱处于不同的平衡状态:当βj < l时,管柱处于稳定状态;当1 < βj < 1.449时,管柱处于正弦屈曲状态;当βj > 1.449时,管柱处于螺旋屈曲状态 [5] 。
(12)
其中,
式中:Ej为节点j处管柱的弹性模量,Pa;Ij为节点j处的极惯性矩,m4;rj为节点j处的半径,m。
2.4. 边界条件
利用式(8)自下而上迭代计算管柱上各点的轴向力时,需首先确定管柱底端及封隔器处的边界条件。
2.4.1. 力边界条件
封隔器坐封时的受力情况如图3所示。封隔器随管柱下到预定位置后,投球封闭管柱下端,第一次加压使封隔器的卡瓦张开锚定在套管内壁上(图3(a)),有:
(13)
在管柱下入过程中,管柱端部不封闭,则
,有:
(14)
式中:Fe,0为一次憋压封隔器受力,N;Ao、Ai分别为与封隔器相接的作业管柱的内、外横截面积,m2;Apo、Api分别为封隔器密封筒的内、外横截面积,m2;
为封隔器上部受到的真实轴力,N;po、pi分别为管柱内、外流体压力,Pa。
当封隔器坐封后,油管内压力变化(二次憋压)引起的管柱变形已经受到限制,封隔器组(包括卡瓦、水力锚、胶筒等)将对封隔器密封筒产生作用力(图3(b)),则
(15)
在酸压、生产过程中,封隔器处的受力如图3(c)所示,有:
(16)
式中:Fe,1为二次憋压封隔器受力,N;Fp为封隔器作用力,N。
由式(11)~(14)可以看出,在封隔器坐封前及一次憋压过程中,等效力为0,管柱不会发生屈曲,在生产和酸压过程中,等效轴力考虑了“虚构力”
,当内压远大于外压时,在管柱下部可能发生屈曲变形。另外,由于Fp未知,需根据封隔器处的位移边界条件来确定。
![](//html.hanspub.org/file/23061x58_hanspub.png)
Figure 3. The setting force of the packer
图3. 封隔器坐封时的受力
2.4.2. 位移边界条件
封隔器坐封后,由于二次憋压、酸压或生产过程中温度、压力的变化,封隔器上方的管柱变形发生变化,引起封隔器处管柱的位移。对于带插管或只带卡瓦的封隔器管柱,管柱只能上移,但不能下移;对于不能移动的封隔器,该处的管柱被固定。这些被限制的位移都将转化为各工况下封隔器的作用力。因此,以坐封状态为基准计算管柱的变形变化量,记为:
(17)
式中:L0为坐封后初始状态的管柱长度,m;L1(s)为坐封后酸压、生产等工况下轴力产生的轴向变形,m;L2(s)为管柱弯曲产生的轴向变形,m;L3(s)为管柱内外压产生的轴向变形,m;L4(s)为温度产生的轴向变形,m。
封隔器对油管的作用力Fp按下面步骤进行迭代计算:
1) 假设封隔器作用力Fp为0,根据公式
计算轴向变形变化总量,此时
不等于0。
2) 由于台阶存在,这个变形变化量
实际上是不可能发生的(图4),令
。
![](//html.hanspub.org/file/23061x64_hanspub.png)
Figure 4. The deformation diagram of the pipe string above the packer
图4. 封隔器上方管柱的变形示意图
3) 重新计算
。
4) 若
不等于0,修正
为
。
5) 重复步骤3)、4),直至
,输出Fp。
3. 典型算例计算及分析
3.1. 已知条件
B101井的三维井眼轨迹如图5所示。一体化油管柱规格为∅88.9 mm × 6.45 mm,材料和钢级分别为G3、125,线重为13.7 kg/m,屈服强度为862 MPa。油管柱下深5020 m,套管内、外径分别为152.5 mm和177.8 mm,封隔器坐封深度4966 m。不同作业的压力和温度参数如表1和表2所示,温度和压力随垂深的增加按线性增大。
![](//html.hanspub.org/file/23061x70_hanspub.png)
Figure 5. 3D well traiectory of B101 well
图5. B101井的三维井眼轨迹
![](Images/Table_Tmp.jpg)
Table 1. The calculation results of wellbore temperature
表1. 井筒温度计算结果
![](Images/Table_Tmp.jpg)
Table 2. Calculation results of wellbore pressure
表2. 井筒压力计算结果
3.2. 载荷、变形及应力分析
油管的轴向应力和等效应力(第四强度理论)分布如图6所示,油管柱在井口处的载荷及不同作业过程中的变形、强度校核结果如表3~5所示。
由图6(a)可见,在下入井筒、一次坐封、二次憋压、酸压增产、生产过程中,油管柱的轴向应力沿着井深方向基本呈线性递减变化,到达一定深度以后,该点处的管柱轴向应力降为0,则该点以上管柱受拉,以下部分受压。由图6(b)可见,按第四强度理论计算的等效应力为三向应力,受周向、径向应力及急剧变化的狗腿度(约井下4800 m)影响,生产过程中呈曲线分布。油管柱的最大轴向应力和等效应力均出现在井口处,管柱在井口处的安全系数最小,是受力的最危险位置。
![](//html.hanspub.org/file/23061x71_hanspub.png)
Figure 6. The axial and equivalent stress distributions of tubing string in different operation processes
图6. 不同操作过程中油管柱的轴向、等效应力分布
![](Images/Table_Tmp.jpg)
Table 3. The axial deformations of tubing string
表3. 管柱轴向变形
表4. 管柱载荷
![](Images/Table_Tmp.jpg)
Table 5. The maximum stresses and minimum safety factors of tubing string
表5. 管柱的最大应力及最小安全系数
由表3、表4可看出,油管内的压力显著影响管柱的变形和应力大小。由于压力增大导致管柱的鼓胀变形增加,管柱轴向缩短,封隔器对上部管柱产生较大的向下拉伸力,增大了管柱在井口处(危险位置)应力,使最小安全系数下降,因此,具有较大作业压力的酸压是一体化油管柱的最恶劣工况,需重点对该工况下的管柱安全性进行分析。目前投入使用的一体化油管柱在不同作业过程中的最小抗屈服安全系数高于1.75 (表5),无正弦和螺旋变形发生,可以满足酸压、生产联作的安全要求。
4. 结论
基于三维弯曲井眼中管柱的几何描述,建立了考虑井身结构、管柱几何、物理特性和作业参数的管柱力学模型,根据一体化油管柱的工作特点,提出了封隔器作用力的详细计算过程,形成一体化油管柱的力学计算方法,确定了B101井在不同工况下的油管柱变形和应力,进行了强度和稳定性评价。研究表明,井口处的管柱安全系数最低,酸压过程相对于管柱下入和生产工况,管内的作用压力较高,使得油管柱的最小安全系数从下入工况时的2.51大幅下降至1.75,基本满足强度要求;同时,油管柱在各工况下无正弦和螺旋等永久变形发生,预示油管柱是稳定的。通过该研究,为一体化油管柱的设计和应用提供了科学依据,有利于保障井下管柱的安全。
基金项目
国家科技重大专项(2016ZX05056004);国家自然科学基金(61572084);湖北省教育厅科学研究计划项目(B2016032)。
NOTES
*通信作者。