井筒环空气液流动力学模型及环空流型过渡准则

井筒环空气液流动力学模型及环空流型过渡准则

0环空气液两相流型的流动方程在大部分石油开采中，井筒环周围的空中有两种气相流动：井筒环上有溢流，井筒环上有两相流。井筒环境中的欠平衡也包括井筒环的空间概念。当自喷井的产量高时，油袋环被用作气生物和素的形状。井筒气液两相流的研究方法可以分为两类:经验模型法和机理模型法。经验模型法中的环空水力学模型多是在圆管经验公式的基础上采用水力学直径校正的方法建立起来的[7,8,9,10,11,12,13,14],误差较大。此外,这种方法不能体现环空和圆管在结构上的差异性,通常忽略流型对参数分布的影响,将两相流看成拟单相流,受限于实验条件而不能广泛应用。机理模型法则是考虑流体的形态及性质、环空的结构及流速的影响,针对每种流型分别建立流动方程。20世纪70年代以来,已有较多学者对流型转换模型进行了研究,但关于持液率及压降预测的模型较少。油管膜、套管膜双层液膜的存在,使得单元体的受力发生了变化,对流型过渡、持液率、压降等的预测影响较大,在研究环空气液两相流动规律时需要加以考虑。本文基于段塞流的动态变化特征,考虑油管液膜、套管液膜、倾角及气芯中液滴的影响,对井筒环空中段塞流、环状流流型下的水动力学模型及各流型间的过渡准则进行研究,并与前人的气-煤油实验数据进行对比验证。1环空两相流中的液膜井筒垂直环空中上升气液两相流流型的分类与圆管中类似,可以分为5种流型(见图1):泡状流、分散泡状流、段塞流、搅拌流及环状流,其中段塞流及环状流与圆管中流型差别最大。环空两相流中的段塞流存在两层液膜,一层是与套管壁(或井壁)接触的套管膜,另一层是与油管壁接触的油管膜。此外,根据CaetanoEF的实验观察,Taylor泡不再对称,由于在背部方向存在一个连接油管膜与套管膜的通道,Taylor泡后存在高紊流区域。环空两相流中的环状流发生在很高的气体流速下,气芯中的气相速度非常高,可能含有液滴,气芯周围是很薄的液膜。由于环空的结构,也存在两种液膜,一种是与油管壁接触的内膜,一种是与套管壁接触的外膜,外膜比内膜厚。2处理工艺-膜区域不存在时段塞流演变规律2003年,ZhangHQ等基于“段塞流体力学”理论建立了适用于各个流型的统一的圆管水动力学模型(本文称Zhang模型),认为各个流型都可以由段塞流过渡而来:当液塞区域不存在时,段塞流演变为环状流;当液膜区域不存在时,段塞流演变为泡状流或分散泡状流。本文在段塞流体力学理论基础上,考虑套管膜及与油管膜两层液膜的影响、气液间的传质作用及气芯中液滴的作用,以液膜区域为控制单元建立环空段塞流水动力学模型。2.1液相质量守恒方程将套管液膜与油管液膜区域作为控制单元(见图2),考虑Taylor泡区域存在液滴,假设液滴速度与液塞运移速度相等,并且控制单元中流体不可压缩。对于稳定的段塞流,从液膜底部边界流入的液相质量与从液膜顶部边界流出的液相质量相等,即:同理,进入液膜的气相质量等于离开液膜的气相质量,即:由(1)式、(2)式可得:若Taylor泡通过某个横截面所用的时间为ΔtTB,则ΔtTB内液膜中流动的液相体积为:若液塞通过相同横截面所用时间为ΔtLS,则ΔtLS内液膜中流动的液相体积为:在时间ΔtTB+ΔtLS内,段塞单元通过相同横截面时流动的液相体积为:由质量守恒定律可得:由(4)~(7)式可得段塞单元中液相质量守恒方程:同理,段塞单元中气相质量守恒方程为:此外,Taylor泡中液滴体积分数可以表示为:2.2环空段塞流动量恒方程对套管膜进行受力分析(见图3),根据动量守恒定律可得:(11)式可变换为:同理,对于油管膜:同理,对于Taylor泡:其中将(12)式、(13)式分别与(14)式联立,可得环空段塞流动量守恒方程:2.3界面接触本构模型①液塞运移速度。NicklinDJ提出液塞运移速度可以表示为混合速度的函数:(17)式中,CS是稳定段塞的最大速度与平均速度之比,随着流动条件的变化而变化。ZhangHQ等认为层流时该值等于2.0,紊流时该值等于1.3。在过渡区域(Re在2000~4000):SantosOLA等建立了流体为钻井液时Taylor泡的上升速度模型:其中②界面摩擦因子。段塞流中液膜与Taylor泡接触的界面会产生摩擦力,AndritsosN等建立了计算界面摩擦因子的模型,本文采用该模型计算套管膜、油管膜与Taylor泡的界面摩擦因子。对于套管膜与Taylor泡接触面:其中vSGt=5(ρG0ρG)0.5对于油管膜与Taylor泡接触面:③壁面剪切力。环空段塞流动量守恒方程中的壁面剪切力计算方法如下:其中fLFc=CRe-nLFcfLFd=CRe-nLFd采用水力学当量直径方法,套管膜、油管膜与Taylor泡界面的雷诺数可以表示为:其中ALFc=HLFcAcALFd=HLFdAc④管壁及界面的周长。当气体流速较高或者倾斜角较大时,管壁处液膜厚度就会增加,导致管壁与液膜界面周长的增大,可以采用Grolman校正公式计算管壁润湿分数:⑤套管膜与油管膜持液率之比。采用CaetanoEF建立的环空环状流持液率模型及套管膜、油管膜厚度比模型计算套管膜与油管膜持液率之比:由于δc、δdp与套管、油管直径相比小很多,则(34)式可简化为:其中⑥压力梯度。整个段塞单元的平均压力梯度可以表示为:液膜区及液塞区的压力梯度均包含重力项、加速度项及摩擦阻力项。摩擦阻力项均包括3部分:套管膜阻力、油管膜阻力及Taylor泡阻力。3环境空广环的水动力学模型环空中的环状流(见图4)可以等效为高速流动的气芯被套管膜、油管膜包围,同时气芯中含有一些小液滴。3.1液滴3部分,质量守恒方程由于液相包含套管膜、油管膜及气芯中液滴3部分,其质量守恒方程为:由于气相仅存在于气芯中,其质量守恒方程为:假设气芯中气液均匀分布,则:3.2环空环形水流动能恒方程(15)式、(16)式等号左边第1项均表示液塞区与液膜区的动量交换,当段塞流过渡到环状流时,不再存在液塞区,即液塞区与液膜区的动量交换变为零,也就是(15)式、(16)式等号左边第1项变为零,即环空环状流下动量守恒方程为:3.3环探气模型的计算周长、水力学直径、横截面积、套管膜与油管膜持液率之比等参数的计算方法参照环空段塞流。①界面摩擦因子。本文采用AmbrosiniW等修正的Asali模型计算套管膜、油管膜与气芯的界面摩擦因子。套管膜与气芯的界面摩擦因子为:油管膜与气芯的界面摩擦因子为:其中WeG=ρGvc2dR/σReG=ρGvcdR/μG②管壁摩擦因子。CaetanoEF基于连续性方程、动量方程及Fanning方程推导了环空中管壁摩擦因子方程。对于层流:对于紊流:其中4循环流转移基准基于不同流型的形成机理,在前人研究成果的基础上建立环空流型过渡准则。4.1采用barnea模型当vSG大于0.1m/s时,采用ZhangHQ建立的分散泡状流到段塞流的过渡模型,并采用水力学当量直径的方法进行校正:其中Tsm=2fSρSvS2/3当vSG小于0.1m/s时,采用Barnea模型:4.2环时空含气率的影响当环空中气体流速较低时,离散的气泡不会发生碰撞及合并,气泡直线上升。当气体流速较高时,气泡开始变大,达到临界尺寸时气泡开始变形并以Z型路线运动,然后开始碰撞、合并,形成球形泡,与段塞流中的Taylor泡类似。这时就会发生向段塞流的过渡。在纯水中,TaitelY等认为当含气率达到25%时,泡状流会向段塞流发生过渡,当气泡间的距离是气泡半径的一半时就会发生急剧的合并。CaetanoEF通过实验研究认为,当环空中含气率为20%时,会发生泡状流到段塞流的过渡。所以,泡状流到段塞流的过渡准则为:4.3向非面者不液塞和液膜间的涤纶配置当环空中的气体流速很高时会发生段塞流到环状流的过渡。液相以液膜的形式沿着管壁向上运动,同时在气芯中还可能以小液滴的形式存在。当段塞流中的液膜变得无限长时,就会发生段塞流向环状流的过渡,使得液塞与液膜间的动量交换变为零。基于ZhangHQ建立的段塞流到环状流的过渡模型得到环空段塞流到环状流的过渡准则:5环空流型及与实验结果对比CaetanoEF在外径76.2mm、内径42.2mm的环空中,对不同温度、压力、表观速度下的气-煤油两相流流型、持液率及压降梯度进行了实验研究。根据CaetanoEF的实验条件,利用本文建立的环空段塞流、环状流水动力学模型及经过当量直径校正后的Zhang模型,对不同液相表观速度、气相表观速度下的流型、持液率及压降梯度进行预测,并与实验结果进行对比。模型求解的具体步骤(以环空段塞流为例)如下:①输入基本参数,先求解vT、lS和FE,估算HLS,给定lF的初值;②基于辅助参数计算方法,求解质量守恒方程,得到VLFc、VLFd、Vc、HLFc、HLFd和HLc;③计算SLFc、SLFd、SIc、SId、fLFc、fLFd、fIc和fId;④计算摩阻和HLS;⑤利用动量守恒方程,求解lF;⑥判断lF是否满足迭代精度,若满足,求解结束,若不满足,返回①重新开始。表1、表2分别为段塞流、环状流下采用本文模型和校正后的Zhang模型预测的流型及与实验结果的对比,可以看出:采用本文建立的模型可以比较准确地预测环空流型,且对段塞流的预测比对环状流的预测更准确。表3—表6分别为段塞流、环状流下采用本文模型和校正后的Zhang模型预测的持液率、压力梯度及与实验结果的对比。可以看出:采用本文模型预测的环空段塞流下持液率、压力梯度与实验结果相比误差均在15%以内,预测的环空环状流下持液率、压力梯度与实验结果相比误差均在25%以内;采用Zhang模型预测的环空段塞流及环状流下持液率、压力梯度与实验结果相比平均误差在50%左右,且在液相相对速度较大、气相相对速度较小时误差较小,在液相相对速度较小、气相相对速度较大时误差很大。因此,采用本文建立的环空段塞流、环状流水力学模型可以比较准确地预测持液率、压力梯度,且准确性优于采用当量直径方法校正后的Zhang模型。6环空段塞流、环形扫码模型分析了井筒环空气液两相流流型的分类及特点,并与圆管气液两相流流型进行了对比。基于段塞流体力学理论,以液膜区域为控制单元,考虑环空段塞流下油管膜与套管膜的影响及气芯中液滴的影响,建立了适合于直井、斜井的井筒环空段塞流、环状流水动力学模型,推导了质量守恒方程和动量守恒方程,并给出了辅助参数计算方法。考虑环空中气液两相流流动特点和环空结构对流型过渡的影响,建立了环空流型过渡准则。利用已发表文献中实验数据对本文建立的环空段塞流、环状流水动力学模型进行了验证,并与经过当量直径校正后的圆管模型进行了对比。结果表明:与实验结果相比,本文模型对流型的预测比较准确,段塞流模型对持液率及压降梯度的预测误差均在15%以内,环状流模型对持液率及压降梯度的预测误差均在25%以内,而校正后的圆管模型对段塞流及环状流下持液率、压降梯度的平均预测误差约为50%。酶活剂参数HLFc——套管膜持液率,%;vLFc——套管膜流动速度,m/s;vT——液塞运移速度,m/s;HLFd——油管膜持液率,%;vLFd——油管膜流动速度,m/s;HLc——Taylor泡持液率,%;vc——Taylor泡流动速度,m/s;HLS——液塞持液率,%;vS——液塞流动速度,m/s;ΔtTB,ΔtLS——Taylor泡、液塞通过某个横截面所用的时间,s;ρLF——液膜密度,kg/m3;Ac——井筒环空横截面积,m2;lU——段塞单元长度,m;lF——液膜区长度,m;lS——液塞区长度,m;vTB——Taylor泡的上升速度,m/s;VLF——ΔtTB内液膜中流动的液相体积,m3;VLS——ΔtLS内液膜中流动的液相体积,m3;VSU——ΔtTB+ΔtLS内液膜中流动的液相体积,m3;vSL——液相表观速度,m/s;vSG——气相表观速度,m/s;FE——Taylor泡区域中液滴体积分数,%;p1,p2——段塞单元入口端、出口端压力,Pa;τIc——套管膜与Taylor泡接触界面的剪切力,Pa;SIc——套管膜与Taylor泡接触界面的湿周,m;τLFc——套管壁与套管膜接触界面的剪切力,Pa;SLFc——套管壁湿周,m;τId——油管膜与Taylor泡接触界面的剪切力,Pa;SId——油管膜与Taylor泡接触界面的湿周,m;τLFd——油管壁与油管膜接触界面的剪切力,Pa;SLFd——油管壁湿周,m;g——重力加速度,9.8m/s2;θ——井筒倾斜角,(°);ρc——Taylor泡区域当量密度,kg/m3;ρG——Taylor泡区域气相密度,kg/m3;ρL——Taylor泡区域液相密度,kg/m3;CS——稳定段塞的最大速度与平均速度之比;Re——雷诺数;dc——气泡直径,m;C1——考虑环空结构的影响系数;K——钻井液稠度系数,Pa·sn;C2——考虑钻井液非牛顿性质的影响系数;Reb——气泡雷诺数;n——钻井液流动特性指数;dco——环空外径,m;dci——环空内径,m;fc——Taylor泡的界面摩擦因子;fIc——套管膜与Taylor泡接触面的摩擦因子;fId——油管膜与Taylor泡接触面的摩擦因子;δc,δdp——套管膜和油管膜厚度,m;ρG0——大气压下的气体密度,kg/m3;fLFc——套管壁与套管膜接触面的摩擦因子;fLFd——油管壁与油管膜接触面的摩擦因子;ReLFc——套管膜与Taylor泡界面的雷诺数;C——局部变量;dLFc——套管膜与Taylor泡界面直径,m;μLF——液膜黏度,m