直井环空多相流井底流动型流动构型模型.docxVIP

直井环空多相流井底流动型流动构型模型 1 多相流压降和确定 在过去的几十年里，多相流规律的研究发生了变化。自18世纪以来，有许多计算模型可以计算管道、环流和静态变换的参数。美国学者提出的Poettman模型的不足之处是忽略了气体的摩擦力和滑脱造成的影响;以后有学者先后提出了水平管中多相流动压降计算和流动型态确定的相关方法,这些方法在充分考虑两相共存时,存在相间的速度差,即滑脱现象,同时提出了两相流动型态的确定方法和HK模式。HK模式从整体看与实际情况吻合好,偏差较小。然而,他们都没有考虑井筒温度和多相流环空岩屑的影响。西南石油学院提出的多相流模型,是利用室内实验方法模拟多相介质流动,推导出了岩屑的最低携带速度,增加了固相对多相流动的影响,但对流动状态的多相流井底压力、循环摩阻以及加速度摩阻没有进行理论和实际工程的验证,而且室内实验受设备条件限制,其理论计算结果存在一定的误差。通过井底压力测量仪器实测数据分析可见,原有的气液两相流计算方法误差在10%左右,最高达到13%,因此有必要研究和建立新的井底压力计算模型。 2 流动常态模型 井筒多相流模型主要有4种:①均相流动模型;②分相流动模型;③漂移流动模型;④流动型态模型。本研究采用半经验方法,在前人研究的基础上考虑了岩屑固相对井底静液柱压力和循环摩阻的影响,并对气液流动加速度压降进行了分析,使本模型更接近欠平衡钻井的实际情况,提高了井底压力模型的计算精度。 2.1 水基钻井液气流系统 根据原有的钻井过程气侵时的两相流模型和欠平衡理论,应用经验公式建立了流钻欠平衡钻井气侵的多相流模型。在流钻欠平衡钻井过程中,当钻遇气层时,由于地层压力高于井底压力,气体就会进入井筒内,随着钻井液一起上返并产生溢流。对于水基钻井液,气体溶解量很小,可以忽略不计。为了描述整个多相流体系需要11个变量,即:气体、岩屑和钻井液的密度;气体、岩屑和钻井液的体积分数;气体、岩屑和钻井液的流速以及压力和温度。在整个流动过程中,可以假设钻井液自身各相之间的整体能量保持不变,这样就可以对与剩余变量相关的10个方程进行理论求解。 2.1.1 岩屑质量守恒方程 钻井液的质量守恒方程为 ??t(ρmξ)+??z(ρmυmξ)=0(1) 气体的质量守恒方程为 ??t(ρgλ)+??z(ρgυgλ)=Qg(2) 岩屑的质量守恒方程为 ??t(ρsδ)+??z(ρsυsδ)=Qs(3) 式中υm为钻井液上返速度,m/s;υg为气体上返速度,m/s;υs为岩屑上返速度,m/s;ρm为钻井液密度,g/cm3;ρg为气体密度,g/cm3;ρs为岩石密度,g/cm3;ξ为钻井液体积系数;δ为气体体积系数;λ为岩屑体积系数;Qs为岩屑体积量,L/s;Qg为随钻产气量,L/s。 2.1.2 井底压力mpa 在模型中的气液两相动量方程中考虑固相颗粒的影响,则流体总的动量守恒方程可以表达如下: ??t(ρmυmξ+ρgυgλ+ρsυsδ)+??z(ρmυ2mξ+ρgυ2gλ+ρsυ2sδ)+(?p?z)+(?p?z)f+(ρmξ+ρgλ+ρsδ)g=0(4) 式中p为井底压力,MPa;g为重力加速度,m/s2。 2.1.3 g0岩屑下降速度m/s 根据文献中的公式,气体随钻井液上升的速度可以表达如下 υg=C0(δυg0+ξνm)(5) 式中C0为经验系数,由流道尺寸决定;υg0为气体在静止液中上升速度,m/s。 岩屑随钻井液的上升速度为 υs=υm-υs0(6) 式中υs0为岩屑下降速度,m/s;υs0是采用Gray推导出来的临界速度约等式。 多相流总的流动速度为 υz=υm+υg+υs(7) 式中υz为混合物流动速度,m/s。 2.1.4 井筒内流体温度t 气体状态方程为 ρg=pΜgΖRΤ(8) 约束方程为 δ+λ+ξ=1(9) 根据文献的井筒内流体的温度变化方程可得 Τ=Τ0+h/A(10) 式中T为井底温度,℃;T0为地面温度,℃;A为井筒过流面积,mm2;h为井深,m;Mg为气体质量,kg;Z为气体压缩因子;R为气体常数。 2.2 降计算的推导 欠平衡钻井井筒环空多相流井底压力梯度公式的推导方法,可根据Beggs-Brill推导的以流动型态模型为基础的压降计算方法进行推导,并增加岩屑固相的影响。 假设外界没有对气液固多相混合物作功,气液固多相混合物也没有向外界作功,根据稳定流动的多相流机械能守恒方程式的微分形式,得出单位质量的气液固多相混合物的计算公式 dpdz=(?p?z)1+(?p?z)2+(?p?z)3(11) 即总压力梯度是多相流摩擦压力梯度、重力压力梯度和加速压力梯度三者之和。 2.2.1 钻井液变形参数影响 在流钻欠平衡钻井过程中,新钻出的岩屑要进入井筒内,并从井底上返至泥浆池中。原为气液两相的体系中,开始掺入固体颗粒