amer模型在气井临界取液中的应用.docxVIP

amer模型在气井临界取液中的应用 0 临界携带液流量模型 如果气体是连续的，液体是分散的。液体以颗粒的形式被气体带到地面。然而，如果气体的速度降低，传输的能力将降低。如果降低到一定程度，井筒中的液体将继续流入井，这样液体将集中在井底，形成积液。为保证气井不积液,气井产量必须大于临界携液流量,因此,气井的临界携液流量是非常重要的一个数据。1969年,Turner等人建立了液滴模型,并将理论结果上调20%来提高与现场的接近程度,该模型得到广泛应用。但是,近年来国内外研究者提出了多种新的临界携液流量模型[4~9],这些模型可以被认为是对Turner液滴模型进行的修正或改进。本文在邵明望等人得到的拽力系数CD和雷诺数之间的全域拟合关联式的基础上,根据Turner液滴理论推导出临界携液流量模型。通过实例分析,计算的临界携液流量更符合实际生产情况。 1 临界流速公式 1969年,Turner针对大于1 400 m3/m3的高气液比,通过对管壁液膜移动模型和高速气流携带液滴模型的比较,认为高速气流携带液滴模型更适合于气井的积液问题。他通过在假设高速气流携带的液滴是圆球形的前提下,推导出临界流速计算公式,并对公式进行修正后得出气井临界流量公式。近年来,国内外许多研究者[4~9]在Turner液滴模型的基础上,提出了多种新的计算模型。Coleman认为对于低压气井没有必要对Turner模型中的临界流速进行20%的修正,而是认为对临界流速公式前的系数进行修正后预测结果比较理想;Nosseir以及刘广峰等人认为雷诺数取值范围不一样时,所取的拽力系数CD不一样,推导出雷诺数在(3×105,5×106)范围内模型;Li Min等人认为液滴在高速气流作用下会由圆球形变成一个椭球形,且椭球形时所取拽力系数CD近似等于1,推导出新的计算模型;王毅忠等人根据流体力学的最新成果,推导出气井携液过程中的液滴形状是以球帽形为主的临界流速公式,并对公式前的系数进行25%的修 2 临界携带液流速公式 根据Turner液滴理论推导出的气井临界携液流速公式为: 其中:a=(4g/3CD)0.25 式中υ——临界流速,m/s; ρl、ρg、ρm——液相、气相和气液混合物的密度,kg/m3; σ——气液表面张力,N/m。 换算成标况下的气井临界携液流量公式: 油管面积: 式中qc——临界流量,104m3/d; T——温度,K; dc——管内径,m; 井筒中天然气的摩尔质量: 式中Mg——天然气的摩尔质量,g/mol; γg——气体的相对密度,无因次。 将式(5)、(6)代入式(2)中简化得气井临界携液流速公式: 将式(4)、(7)代入式(2)中简化得气井临界携液流量公式: 根据实验数据绘制出的绕流阻力系数曲线图,在不同的雷诺数范围,会引起液滴形状呈现不同的状态,且会影响液滴的拽力系数CD值的计算。邵明望根据实验数据利用非线性拟合方法得到一种拟合程度较高的关系式: 式中Re——雷诺数,无因次。 根据文献可知,该公式是全域拟合关联式,公式的拟合程度较高,误差较小,因此,笔者在计算临界携液流量时利用这个公式计算拽力系数值。根据气井的实际参数先计算出流体的雷诺数,再依据域拟合关联式计算拽力系数CD值,再根据临界携液流量模型,计算气井的临界流量,来判断气井是否积液。 3 预测结果分析 在相同条件下,利用本文中总结的国内外研究者推导出的气井临界携液流量模型和文中推导的临界携液流量模型对气井进行了连续临界携液流量的对比计算,并预测气井的生产状态,实例对比分析结果见表2。气井的部分生产数据见表2,一般天然气的相对密度0.58~0.62,计算时天然气的相对密度取0.6,气井的井口温度为322 K,气井的井口压力根据气井的实测数据,气的临界温度为193 K,临界压力为4.6MPa,采气管柱的内径为0.062m。其中天然气偏差系数Z采用D-A-K方法计算,同时为了计算方便,对于水的密度和表面张力,可采用以下经验数据,气—水间的界面张力取0.06 N/m,水密度为1 074 kg/m3。 从表1可以看出模型1和模型3的计算结果偏大,预测结果的误差较大,而模型5的计算结果偏小,但是与实际生产有误差,文中推导的气井临界携液流量模型与模型4的计算结果相近,但文中推导模型的预测结果与实际生产较一致。文中推导模型是根据气井的实际气井参数计算出雷诺数,根据拽力系数CD和雷诺数之间的拟合的关联式(式(8)),计算出颗粒的拽力系数CD值,再利用临界携液流量模型计算出的临界流量,因此,推导模型的预测结果与实际生产较一致。 假设气井产量等于其临界携液流量,对角线为基准线,那么正常生产气井(圆形符号)应处于对角线上方,接近积液的气井(菱形符号)应分布在对角线附近,而积液井(三角形符号)应在下方。从图1中可以说明临界携液