2 气液两相管流2009

图2-1自喷井生产系统及压力损失11/8/20231

第二节气液两相管流基本概念及基本方程两相管流：占油气井系统总压降35%-90%核心问题：沿程压力变化及其影响因素11/8/20232一、气液两相管流的滑脱现象及特性参数滑脱现象气液两相上升流动时，由于气液两相间的密度差异所产生气相超越液相流动。相对于气相而言，有一部分液相被滞留于管段中。密度增加，压降损失增加11/8/20233

持液率(LiquidHoldup)

流动状态下单位长度管段内液相容积所占份额无滑脱持液率(No-slipLiquidHoldup)

管流截面上液相体积流量与气液混合物总体积流量的比值。（条件：）

关系11/8/20234流速相平均流速实际流速表观流速关系混合物流速滑脱速度11/8/20235混合物密度存在滑脱无滑脱滑脱损失

11/8/20236二、气液两相管流的流型

纯液流（p>pb）无气相，管内均质液体流体密度最大，压力梯度最大溶解气开始从油中析出，气体以小气泡分散在液相中泡流（p<pb）液相是连续相，气相是分散相液相滑脱损失严重，易水淹摩阻小，重力损失为主特点11/8/20237液相是连续相，气相是分散相气体体积变大，摩阻增加滑脱较小，总压力损失最小特点混合物继续向上流动，压力降低，小气泡合并形成大气泡，在井筒中形成一段气，一段液流动结构，气托着油向上运动。段塞流11/8/20238液相由连续相过渡为分散相，气相相反气体流量大，摩阻增加特点压力继续降低，部分大气泡从中间突破液段形成短气柱，把液体挤到环壁。液体靠中心气流的摩擦携带作用向上运移。过渡流（环流）11/8/20239气相是连续相，液相是分散相摩阻增加，重力损失最小特点压力进一步降低，中心气柱逐渐增大，壁面液膜厚度降低，液体以液滴分散于气相中。雾状流11/8/202310纯液流泡流段塞流过渡流雾状流油井生产中可能出现的流型自下而上依次为：纯油流（液流），泡流，段塞流，过渡流，雾状流。实际上，在同一口油井中，一般不会出现完整的流型变化。总结Hp11/8/2023112.水平管流气液两相流流型分层流上部气流、下部液流气液界面平滑或波状中心气流，携带液滴管壁液环流动11/8/202312大气泡沿管子顶部流动，管子下部为液流间歇流塞流段塞流大液体段塞流与几乎充满管子的高速气泡的交替流。11/8/202313气流量高、液流量低气流中夹带液滴分散流大气泡集中在管子的上半部。11/8/202314四、气液两相管流压力梯度方程及求解步骤1.压力梯度方程压降梯度=重力梯度+摩阻梯度+动能梯度15单相流多相流水平管流(θ=0)，且忽略动能？11/8/202316分析泡流段塞流过渡流雾流井筒L处p，TL=L1,L2,L3---Lnp=p1＋△p1+△p2---△pn-111/8/202317迭代计算步骤pwf，Twfp1=pwf－△p1？p1，T2p1=pwf－△pxpwfBo,Bg,Bw,Rs,μo,μg等vm,vsg,vsl判别流型计算HL(ρm)泡流段塞流过渡流雾流泡流段塞流过渡流雾流动能＋摩阻＋重力△px2△px,-

△px2<Ep1=pwf－△px2△H11/8/202318

第三节气液两相管流计算方法早期均匀流方法（总摩阻系数法）1952Poettmann—Carpenter80’s陈家琅λ'~(NRe)2经验相关式1963Duns－Ros无因次化处理NvL、Nvg、ND、NL1965Hagedorm－Brown现场实验1967Orkiszewski流型组合1973Beggs－Brill倾斜管实验1985Mukherjee－Brill改进实验条件现代机理模型1985Hasan&Kaber1990Ansari11/8/202319

Orkiszewski方法

Orkiszewski（1967）采用148口油井实测数据，对比分析了多个气液两相流模型。然后分不同流型择其优者，综合他的研究成果得出四种流型的压降计算方法。流型选用方法泡流段塞流过渡流雾状流Griffith和Wallis密度项对Griffith和Wallis公式作了修正，摩阻项用Orkiszewski方法Ros和DunsRos和Duns11/8/2023201.总压降梯度公式一般动能较小，只在雾流情况下才有意义。只考虑气体的压缩性：？11/8/202321伴随生产1m3地面脱气原油产出的油、气和水的总质量，kg/m3。总压降梯度公式？11/8/2023222、流型判别多相管流流态的影响因素共有13个，主要因素：VSL、VSG、ρL、σ1）影响流态的因素2）无因次处理11/8/202323Ros流型图版11/8/202324NGV<LS段塞流NGV>LSNGV<LM过渡流NGV>LM雾状流段塞流与过渡流界限值为:雾流与过渡界限值为:？泡流与段塞流分界25vmqG/qm1.00.1323/8段塞流泡流41/2431/227/81.9泡流与段塞流分界11/8/202326

Orkiszewski方法流型界限流型界限泡流qG/qm<0.13或0.13<qG/qm<LB段塞流qG/qm>LB，NGV<LS过渡流LM>NGV>LS雾流NGV>LM27四、混合物密度与摩阻梯度的计算a.混合物的密度1）泡流与滑脱速度有关实验表明：泡流时vs=0.244？11/8/202328p,T状态标准状态生产油气比溶解油气比11/8/20232911/8/202330b.摩阻梯度泡流中气体以小气泡分布于液体中，靠近管壁主要是液体。其摩阻压力梯度按液相计算。？紊流（NRe>2300）

层流（NRe≤2300）

11/8/2023312.段塞流a.混合物密度？液体分布系数Co由连续液相的类型及混合物速度分别选用相应的公式连续液相vm，m/sCo计算公式水<3.0481－86a水>3.0481－86b油<3.0481－86c油>3.0481－86d11/8/202332vs计算方法一11/8/202333vs计算方法二当Nb≤3000时当3000<Nb<8000时当Nb≥8000时11/8/202334b.摩阻梯度11/8/2023353.雾状流a.混合物密度雾状流一般发生在高气液比、高流速条件下，液相以小液滴形式分散在气柱中呈雾状，这种高速气流携液能力强，其滑脱速度甚小，一般可忽略不计。36b.摩阻梯度？37确定e/D

根据无因次韦伯系数选择公式。当Nw≤0.05时当Nw>0.05时液膜的相对粗糙度，取0.001－0.511/8/2023384．过渡流

用段塞流和雾流计算后内插。39qo=38m3/d，qg=2027.4m3/d，γg=0.85γg=0.65，pb=8.66MPa，pwh=2.352MPaTwh=25℃，ΔT=3℃/100m，dti=62mm。试用Orkiszewski方法计算井口压力梯度。例1-6解：（1）以井口或井底为起点(由已知压力位置定) 起始点：井口压力p1=pwh，T1=Twh，H1=0（2）选择计算区间长度:ΔH一般取50～100m

选取计算区间长度：ΔH=100m

（3）假设这一区间的压降值ΔP(由经验定)

假设深度ΔH对应的压力增量ΔP=0.6MPa11/8/202340（4）计算出区间的平均温度和平均压力Pav,Tav

Pav=

P1+ΔP/2 Tav=T1+

ΔT×ΔH/2（5）确定Pav和Tav下的物性参数

计算Rs、Bo、μo、σo、Zg、μg、Bg等（6）判断流态

计算NGV、LB、LS、LM、qG/qm等，利用表1-9判断流态（7）计算dp/dh和ΔP′

根据流态计算ρm、f等，计算dp/dh则：ΔP′=dp/dh×ΔH11/8/202341（8）比较ΔP′与ΔP,若相差超过误差限,以ΔP

代替ΔP返回到（4）重新计算到第8步。 如果，则进行下一段计算:P1=P1+ΔP,H1=H1+ΔH,T1=T1+ΔT×ΔH直至井底。11/8/202342二、倾斜（水平）管两相流计算方法Mukherjee和Brill（1985）实验装置：内径38mm的倒U形倾斜管，中部可以升降，可与水平方向在0～90°范围内变化实验介质：空气+煤油或润滑油温度：-7.8~55.6℃σ，mN/mρ，kg/m3μ，mPas煤油268172润滑油3584929根据所测得的1500个实验数据，通过回归分析，提出了倾斜管气液两相流的持液率及摩阻系数经验公式。适于垂直井、斜直井、定向井和水平井的两相管流压力计算。

11/8/202343Mukherjee和Brill（1985）对比：11/8/202344与水平方向的夹角（0～+90o）a.持液率流向向上和水平流向下流流型所有分层流其它系数值c1-0.380113-1.33082-0.516644c20.1298754.8081390.789805c3-0.1197884.1715840.551627c42.34322756.26226815.519214c50.47756860.0799510.371771c60.2886570.5048870.39395211/8/202345无因次液相粘度

无因次液相速度

无因次气相速度

11/8/202346若则为环流，否则为泡流–段塞流。条件：向上或水平（水平井段）流动分泡流–段塞流和环流，判别式为b.确定摩阻系数泡流–段塞流11/8/202347

fm为相对持液率HR和无滑脱摩阻系数fns

的函数确定步骤：(1)计算相对持液率HR＝λL/HL

(2)根据HR按下表确定摩阻系数比fR；(3)根据NRens由摩阻系数公式计算f，即为fns；(4)fm＝fR·fns

。环流HR0.100.200.300.400.500.701.00fR1.000.981.201.251.301.251.0011/8/202348三、环形空间流动的处理方法模型：采用圆管内气液两相管流压降计算方法修正环形空间流动：管径和管壁粗糙度相当粗糙度eo,ei分别为环空内管、环空外管有效粗糙度。eo应考虑油管或抽油杆接箍的局部摩阻影响。11/8/202349圆管时，Di=0，故R=D0/4，水力相当直径De=4R。所以环空的水力相当直径为：水力当量直径将圆管管径采用环空的水力半径代替。11/8/20235013-405井压力计算结果对比13-117井压力计算结果对比11/8/20235111/8/202352Hagedorn－Brown方法无滑脱滑脱压降11/8/202353

第四节油井井筒传热模型及温度计算

一、油井井筒传热模型将流体在井筒油管内流动考虑为稳定的一维问题。能量方程dq/dt？11/8/202354焓是工质在某一状态下所具有的总能量（内能U与压力势能ＰＶ之和，为一个复合状态参数。比焓梯度焦耳－汤姆逊系数11/8/202355流体-地层稳定传热油管隔热层套管地层流体环空水泥环rtirhrcorcirtorinsTfTtiTtoTinsTciTcoTh忽略油管内壁水膜及金属的热阻

？11/8/202356热流梯度方程地层内不稳定传热11/8/202357热流梯度比焓梯度能量方程井筒温降梯度方程11/8/202358松弛距离A

任意流通断面的地温（静温）按井筒内流体流动温度梯度gf，折算到流温曲线所产生的相对距离。

Te(z)Tf(z)gfA11/8/202359井筒温降计算需要确定油套环空流体和水泥环及周围地层的一系列物性参数。为此，Ramey（1962）、Satter（1965）、Shiu&Beggs（1980）、Hasan&Kabir（1990）等提出了多个井筒温度简化计算方法。11/8/202360松弛距离A

为产出流体质量流量、管径、产出流体物性和油压的函数。应用370口油气井（直井、定向井）现场测温资料进行线性回归处理得到了A的系数。

二、Shiu&Beggs温度计算方法

考虑油井在稳定生产情况下，上述物性等参数变化不大，故均视为常数，导出任意z截面的温度

11/8/202361解：

【例1-7】油层中深3000m处温度为82℃，地温