在气田开发过程中,井底积液是最常见的问题之一。积液会导致气井产量下降,甚至会造成气井停产,临界携液流量的准确计算对于判断气井积液十分重要[1-3]。Turner模型[4]和李闽模型[5]在现场比较常用,二者均是直井携液模型,将曳力系数取为0.44。杨文明等[6]分析了天然气斜井中液滴的受力情况,假设液滴沿井筒方向运动,不与管壁发生摩擦,这在斜井的实际生产中显然是不合理的。李丽等[7]以Turner模型为基础,考虑井斜角的影响,建立了斜井携液模型,但该模型未考虑曳力系数变化对临界携液流量的影响。Turner假设气井中液滴的雷诺数(Re)范围为104<Re<2×105,但在实际生产过程中,气井雷诺数范围为104<Re<106。针对这一问题,Nosseir等[8]和刘广锋等[9]将雷诺数范围分为两段,当雷诺数范围为104<Re<2×105时,将曳力系数取为0.44,当雷诺数范围为2×105<Re<106时,将曳力系数取为0.2。根据曳力系数与雷诺数关系曲线可知[10],当雷诺数范围为104<Re<2×105时,曳力系数比较稳定,在0.44左右;当雷诺数范围为2×105<Re<106时,曳力系数先下降后上升,变化较大,此时,将曳力系数取为0.2并不合理。一些学者[11-13]考虑到曳力系数是雷诺数的函数,引入曳力系数计算模型计算气井中液滴的曳力系数,但他们所引用的曳力系数计算模型的适用范围为Re<2×105,并不能用来计算2×105<Re<106时的曳力系数。本研究首先对天然气斜井中液滴进行受力分析,认为液滴最终将沿气井管壁滑动,然后对比不同曳力系数计算模型的精度,优选出Barati模型来计算气井中液滴的曳力系数[10],最后建立了预测天然气斜井临界携液流量新方法。
在天然气斜井中,液滴受力如图 1所示。根据液滴质点受力分析可知,如果液滴在井筒中间运动,液滴会受到曳力、浮力和重力3个力作用,由于液滴所受浮力和重力是沿垂直方向的,而曳力是沿井筒方向的,故液滴在这3个力作用下无法平衡。曳力垂直分力、浮力和重力在垂直方向上达到平衡,液滴在曳力水平分力的作用下,不断地在水平方向上向井口运动,直至靠近气井管壁位置,进而沿着气井管壁发生滑动。当天然气斜井中的液滴与气井管壁发生接触后,将受到曳力、浮力、重力、支撑力和摩擦力的作用。
根据牛顿第二定律,当达到临界状态时,天然气斜井中液滴受力达到平衡,会沿着气井管壁匀速移动,此时液滴沿井筒方向的受力关系如式(1)所示。
式中:FD为曳力,N;Fb为浮力,N;Fg为重力,N;Fτ为摩擦力,N;α为井斜角,(°)。
液滴受到的曳力FD、浮力Fb与重力Fg计算公式如式(2)所示。
式中:d为液滴直径,m;ρG为气体密度,kg/m3;ρL为液滴密度,kg/m3;v为气流速,m/s;CD为曳力系数;g为重力加速度,m/s2。
液滴所受气井管壁的摩擦力遵循牛顿内摩擦定律,在管壁处流体切应力τ如式(3)所示。
式中:τ为管壁处流体切应力,Pa;Δp为液体流经长度为l的管道压力损失,Pa;r为管道半径,m;Re为液滴雷诺数。
结合式(3), 可得液滴所受摩擦力Fτ,如式(4)所示。
当气流能把直径最大的液滴带出井口时,气井可以连续携液,井筒就不会积液[14]。气井中液滴的最大直径是由韦伯数决定的,一般来说,气井中液滴的临界韦伯数为20~30。前人研究表明[15],当液滴的韦伯数超过临界值时,液滴就会发生破碎。所以,气井中液滴最大直径可由式(5)计算得出。
式中:σ为气液界面张力,N/m。
将式(2)、式(4)和式(5)代入式(1),可得临界携液流速,如式(6)所示。
临界携液流量公式如式(7)所示。
式中:qc为临界携液流量,m3/d;A为油管截面积,m2;Z为气体偏差因子;T为温度,K;p为压力,MPa。
将式(6)变换为式(8),即本研究新模型。
式中:C为修正系数,括号内为Turner模型计算结果。
理论和实验研究表明,在不可压缩流体中,光滑球体的曳力系数是雷诺数的函数。前人已经研究出很多用于计算光滑圆球曳力系数的公式,其中应用较为广泛的曳力系数计算模型如表 1所列。
为了评估表 1中各曳力系数计算模型的精确度,引入对数偏差平方和(SSLD)、对数偏差均方根(RMSLD)与相对误差和(SRE),如式(15)~式(17)所示。
采用文献[21, 22]中的实验数据将各曳力系数计算模型进行了对比,结果见表 2和图 2。
由表 2可知,一些模型具有较大的SSLD、RMSLD和SRE值(即低精度),如Brauer模型、邵明望模型和Clift模型。这是因为这些模型的适用范围为Re<2×105,作为拟合关系式,它们没有拟合2×105<Re<106的实验数据(见图 2),故不能用来预测2×105<Re<106的曳力系数。在实际生产中,气井Re一般为104<Re<106,一些学者[11-13]引入邵明望模型[16]和Brauer模型[17]计算气井中液滴的曳力系数显然是不合理的。适用范围为Re<106的3个模型中,Barati模型具有最低的SSLD、RMSLD和SRE值,精确度最高。因此,本研究临界携液流速模型中曳力系数采用Barati模型进行计算。
由式(8)可知,修正系数由井斜角、雷诺数和曳力系数决定,而曳力系数是雷诺数的函数。图 3为井斜角和雷诺数对修正系数的影响。由图 3可知,随着井斜角的增大,修正系数逐渐减小,与之对应的临界携液流量也相应减小。当104<Re<2×105时,雷诺数和曳力系数对修正系数影响不大,一些学者[11-13]考虑到曳力系数随雷诺数变化而变化,引入适用范围为Re<2×105的曳力系数计算模型来计算气井中液滴的曳力系数,实际上对临界携液流量的计算影响不大;而当2×105<Re<106时,雷诺数和曳力系数对修正系数影响很大,随着雷诺数的增加,修正系数先增大后减小,与之对应的临界携液流量也先增大后减小。因此,引入适用范围为Re<106的Barati模型计算气井中液滴的曳力系数是十分有必要的。表 3为修正系数速查表,实际运用时,可根据气井实际情况进行修正。
基于文献[23]中的37口积液气井数据,分别采用本研究模型和几种常用的斜井携液模型计算天然气斜井临界携液流速,结果见图 4。假设气井气流速等于其临界携液流速,对角线为基准线,则积液气井临界携液流速计算值应处于对角线的上方。由图 4可知,杨文明模型计算结果明显偏小[6],Belfroid模型预测准确率仅为24.3%[24],李丽模型预测准确率为45.9%[7],而本研究模型仅误判2口井,准确率为94.6%。
采用本研究模型和几种常用的斜井携液模型对文献[7]中2口临近积液气井数据进行计算,计算结果见表 4。由表 4可知,本研究模型计算结果较几种常用斜井携液模型更为准确。
由图 4和表 4可知,本研究模型计算结果具有较高的精度,与气井实际情况吻合较好,从而验证了本研究天然气斜井携液模型的可靠性,可用于对天然气斜井积液的判断。
3结论
(1) 对天然气斜井中的液滴进行受力分析,认为液滴最终沿气井管壁发生滑动。
(2) 考虑曳力系数变化对临界携液流量的影响,采用实验数据对比不同曳力系数计算模型的精确度,优选出Barati模型计算气井中液滴的曳力系数。
(3) 考虑井斜角和曳力系数,对Turner模型进行了修正,并给出了修正系数速查表。结合实例分析,本研究模型计算结果具有较高的精度,准确率为94.6%。