凝析油气井沿井筒压力梯度计算方法.docVIP

中国工程热物理学会 多相流 学术会议论文 编号：046056 凝析油气井沿井筒压力梯度计算方法 李明忠’，石锋2，崔红霞。 (1．中国矿业大学力学与建筑T程学院．北京100083：2．石油大学石油工程学院，山东东营 257061；3．胜利油田胜利采油厂．山东东营257051) 电话：0546-8396444，E—mail：limzh@mail．hdpu．edu．C11 摘要 依据能量守恒定律，建立了垂直气井井筒温度分布的计算模型；借用一定的气液两相管流 理论和传热学知识并结合流体相平衡热力学闪蒸计算方法，得出了凝析气井沿井筒压力梯度的计算模 型及物性参数模型；给出了模型的求解方法；用计算实例对建立的模型进行了验证和评价。 关键词 凝析气井压力梯度多相流相态状态方程 1引言 沿井筒压力梯度的计算，是油气井举升方式选择、工艺参数设计以及工况分析的重 要基础。对凝析气井井筒压力梯度的计算，目前常用的方法(如Sukkar-Cornell方法川 等)中，没有考虑流体组分的变化或只是简单将流量、流体物性参数等进行一定校证来 近似处理，这势必会影响其计算精度。因为实际井筒内自井底至井口压力和温度变化较 大，凝析油气的相态及组成变化也必然很大。本文，依据能量守恒定律，建立了气井井 筒温度分布的计算模型：借用一定的气液两相管流理论和传热学知识并结合流体相平衡 热力学闪蒸计算方法，得出了凝析气井沿井筒压力梯度的计算模型及物性参数模型；用 计算实例对建立的模型进行了验证和评价。 2相关模型的建立及参数计算方法 2．1凝析油气混合物在井筒中的流动情况 凝析气井中流体的流态对井筒压力和温度变化较为敏感，在计算井筒压力梯度时需 考虑流体相态变化和温度、压力的相互作用关系。凝析油气混合物在井筒中的流动可分 为三种情况： (1)若计算点的压力和温度高于露点压力和温度，则该位置无凝析油析出，仅为 单相气流。这时气相的组成为井流物的总组成．气相的摩尔分数为1，液相的摩尔分数 为0。 (2)当计算点的压力和温度略低于露点压力和温度时，将有少量凝析油析出。这 时气液两相的摩尔分数将发生变化．且气体携液滴向上运动。 (3)若计算点的温度和压力大大低于露点压力和温度，凝析油析出量加大且为气  般悄相流动。 2．2气液相平衡计算 气液相平衡计算是将平衡闪蒸方法f21应用到凝析气藏：肖温度压力发生改变时。它 可准确计算出气液相的组分的含量，为物性参数的预测提供依据。 2．3井筒压力分布方程 气体从井底到井rl沿油管稳定垂直管流能量方程为 一塑：∥立+艘sin口+厶生 (1)¨’ 一磊。∥磊+昭81眦+亏■ 其中 ZRT 2。一y面一484．5，，刍 2 ㈤‘2’ V蛆k 2面PjcZT×砜q-％7．=1．25x10-9净等 ㈤ 令G．=1．25×10。9 q,c／#，G2=3484．5y，F=ZT／P，将式(2)、(3)代入 L 1)并推导得出井筒压力分布差分方程： ‰=只+?@¨聃型哮坦×甓’ ．厶垡鱼盥：垒型1 k2+‘‰】×10“ (4) 8，f『 式中 ．P一凝析气压力，MPa： p一凝析气密度，k咖3； v一凝析气流速，rids； g一重力加速度，g=9．8m／s2； 口一井身轴线与水平面的夹角，度： 厶一摩阻系数，无因次； 0一油管内径，m； 肘、M《一凝析气、空气的分子量  Z一凝析气压缩因子，无因次： R一通用气体常数，这里取0．008314； 丁一凝析气温度，K： y一凝析气的相对密度，无因次： B一气体体积系数，无因此； sc一标准状况： q一凝析气体积流量，m3／d； in、out一计算段入口、出口： h一计算点到井底的距离，m。 2．4井筒流体温度分布方程13,4．5] 井简内流体的温度分布是计算其物性参数以及流动压力降的重要基础。根据传热学 原理、井筒结构及流体性质等的综合变化情况推导得出温度分布方程： 7k，=瓦。。+gr sina／A+eA(k一“一1《r加一瓦。)一gr sin酬彳】 (5) 其中4=卺而r厕,oU,kP，t=‰喝硒n口 式中 ￡一任意处的地层温度，K； L一井筒内流体温度，K； gr一地温梯度，K／m； 女。一地层导热率，J／(m·Kh 砂一总传热系数，w／(m·K) ‰一油管半径，m： m一质量流量，kg／s f(t)--Remay无因次时间函数： 7k一井底地层温度，K C。一井筒内流体的定压比热，J／(kg‘K)。  Cp：!!!￡!!±!!旦!!! q sPg+q oPn 其中，gg、g。一凝析气和凝析油的质量流量，kg／s；pg、p0一凝析气和凝析油 的密度，kg／m3：Cg、C。一凝析气和凝析油的定压比热，J／(kg。K)。 2．5井筒内流体相对密度的计算 井筒内混合流体的相对密度16J为 ，