垂直管流中的气液两相流压力计算
气液两相流实验指导
实验三气液两相流实验气液两相流是近几十年发展起来的一门新学科，在热能、动力、化工、核能、制冷、石油、冶金、航空航天、气力输送、液力输送、叶轮机械、生物技术、电子设备冷却等领域均有重要应用，已经成为研制、设计和运转这些重要工业关键设备的必备理论知识。
通过气液两相流的实验研究，是掌握气液两相流规律的基本方法。
本实验指导书根据目前已有的科研成果和国内外有关的成就，结合热能工程专业特点，针对大型电站锅炉中的水动力问题，制定如下实验内容：①垂直上升管中气液两相流的流型和管内气液两相流的压力降；②倾斜管中气液两相流的流型和管内气液两相流的压力降；③气液两相流流经孔板的流型；④气液两相流流经文丘里管的流型；⑤水平集箱和垂直并联管的管道系统通过以上实验内容，希望能达到下列目的：①了解大型电站锅炉中的水动力特性和两相流基本现象；②能够从基本原理与动手实践的角度切实训练学生进行实验的基本能力，使学生知其然、也知其所以然；③使学生从实验设计、仪器选型、实验操作、数据提取与分析处理等各个环节能够训练出真正的实验技能，能够完成合格的实验报告。
实验1 垂直上升管中气液两相流特性实验一、实验目的：1. 在大型电站锅炉中垂直布置的锅炉水冷壁管被广泛应用，本实验将模拟其两相流现象和水动力特性；2. 通过观察垂直上升管中气液两相流的流型,进一步加深了解垂直上升管中气液两相流型的特点；3. 对垂直上升管中气液两相流的压力降有比较直观的认识，并掌握垂直上升管中气液两相流的压力降的计算方法；二、实验仪器：仪器名称型号参数范围磁力泵50CQ-50 130L/min空气压缩机V-1.2/10 1.2m3/min电磁流量计MF/E2004011100EH11 282.6 L/min涡轮气体流量计CP 32700-10 1-5L/min涡轮气体流量计CP 32700-16 5-50 L/min涡轮气体流量计CP 32700-22 50-500 L/min差压变送器1151DP4E22B3 10KPa差压变送器1151DP5E22B3 100KPa压力变送器1151GP6E22B3 300KPa三、实验原理图：11164445231298101371381 水箱2 空气压缩机3 磁力泵4 涡轮流量计 5电磁流量计 6 气液混合器7 减压阀 8 调节阀 9截止阀 10球阀 11 水集箱 12 针阀 13 过滤器四、实验任务：1.观察垂直上升管中气液两相流的流型：(1)打开系统电源，使气体、液体流量计预热2分钟；然后打开采集程序，记下采集程序上显示的气路和水路温度（根据此温度查出水和空气的密度）；(2)打开磁力泵，将主路的调节阀开度调小和旁路的调节阀开度调大,同时将垂直上升管实验段水路的球阀开启，使水缓慢地流过实验段，直到取压管内大体上充满水为止；(3)关闭磁力泵和水路的球阀,打开空气压缩机和气路的球阀,将50-500L/min涡轮流量计一路的针阀开启，调节针阀开度,使涡轮气体流量计所测得的体积流量保持在300L/min；打开磁力泵，调节主路和旁路的调节阀开度,将主路阀门开度达到最小，旁路阀门开度达到最大。
1 两相垂直管流实验
两相垂直管流实验气举井及绝大多数自喷井的油管中流动的都是油—气或油—气—水三相混合物。
对采油来说，油、气、水混合物在油管中的流动规律——多相垂直管流理论是研究自喷井、气举井生产规律的基本理论。
在许多情况下，油井生产系统的总压降大部分是用来克服混合物在油管中流动时的重力和摩擦损失。
它不仅关系到油井能否自喷，而且决定着用自喷和气举方法可能获得的最大产量。
为了掌握油井生产规律及合理地控制和调节油井工作方式，必须熟悉气—液混合物在油管中的流动规律。
在油气田开发过程中，为了充分利用天然资源和取得好的经济效果，或者要进行油气田动态分析，拟订油气田的增产及提高油气田采收率，高速度、高水平地开发油气就必须深入细致地研究地层—油管—油嘴生产衔接与协调，研究多相流在井筒中的流态变化。
使生产井的工作制度同地层变化了的情况协调起来，只有通过各个生产井的各种变化并把它们综合起来进行分析，才能为整个油气田动态分析提供准确的资料和依据，并对各个注采井提出有效的工艺措施，不断完善开发方案，改善油气田开发效果。
该实验就是研究气、液两相在垂直井筒中的流态变化及观察模拟井筒气体膨胀能参与举升液体的现象，抓住观察到的现象综合分析，并对所作的气量与液量的关系曲线作出解释。
一、实验原理在多相垂直管流中，沿井筒自下而上随着压力不断降低，气体则不断从液体中分离出来，以及压力降低气相体积流量逐渐变大。
随着液气流沿井筒上升，压力逐渐降低气体随之膨胀，不断释放出气体弹性膨胀能量，该能量要参与举升液体，膨胀能的大小与气量多少、压力变化范围有关。
该实验是研究液气两相在模拟垂直井筒中的流动变化。
也是利用气体膨胀能量来举升液体的实验，它依靠两种作用：一种是气体作用于液体上，垂直地顶推液体上升；另一种是靠气体与液体之间的摩擦作用，气体携带液体上升。
其能量来源除压能外，气体膨胀能是个很重要的方面。
因在管径不变的油管中，举升一定的油量，则单位管长上所消耗的总压头，是随着气量的不同而变化的，而只有在某一气量下，举升一定气量的液体所必须消耗的压头最小。
2-2.井筒气液两相流
• 出现雾流时，气体体积流量远大于液体体积流量。根
据气体定律，动能变化可表示为：
mvmdvm
Wt qg Ap2 P
dp
• 取 dh hk dP pk
m m
PP
•则
dp [ m g f ]dh
1
Wt qg Ap2 P
表2-2 流型界限
流动型态 泡流 段塞流 过渡流 雾流
界限
qg qt
LB
LM
75 84(vg
qL )0.75 qg
2.3.2平均密度及摩擦损失梯度的计算
(1)泡流 平均密度:
m H L L H g g (1 H g )L H g g
HL Hg 1
气相存容比(含气率)Hg ：管段中气相体积与管段容
积之比值。
液相存容比(持液率)HL ：管段中液相体积与管段容
m
LM LM
vg LS
SL
vg Ls LM LS
Mi
t
LM LM
vg LS
SL
vg Lg LM LS
Mi
图2-23 C2～NRe’曲线
(4)雾流
• 雾流混合物密度计算公式与泡流相同：
m H L L H g g (1 H g )L H g g
• 由于雾流的气液无相对运动速度，即滑脱速度接近于雾，
基本上没有滑脱:
Hg
qg qL qg
• 摩擦梯度:
f
f
g vs2g
2D
• 雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度由图查
得。
NRe
g
gvsg D g
2.4计算井筒多相管流动的Beggs-Brill方法
• Beggs-Brill方法可用于水平、垂直和任意倾斜气液两 相管流动计算。
多相管流-第三章垂直气液两相管流计算1
第三章 垂直气液两相管流压力梯度计算模型及方法垂直气液两相管流压力梯度计算模型及方法•3.13.1流动模型流动模型•3.23.2压力分布计算方法压力分布计算方法•3.33.3垂直气液两相管流压力梯度计算模型及垂直气液两相管流压力梯度计算模型及方法•3.43.4水平或倾斜管中气液两相流动计算模型水平或倾斜管中气液两相流动计算模型及方法•3.53.5水平管中气体和非牛顿液体的两相流动水平管中气体和非牛顿液体的两相流动主要内容第一节流动型态流动型态的划分方法：两类第一类划分方法：根据两相介质分布的外形划分泡状流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、环状流、雾状流第二类划分方法：按流动的数学模型或流体的分散程度划分分散流、间歇流、分离流分散流、间歇流、分离流2010-3-263垂直气液两相流流型水平气液两相流流型两种分类方法比较第一类划分方法较为直观第二类划分方法便于进行数学处理第一类划分方法•泡状流•弹状流或团状流•层状流•波状流•段塞流或冲击流•环状流•雾状流第二类划分方法•分散流•间歇流•分离流•分离流•间歇流•分离流•分散流两类划分结果的对应关系2010-3-264垂直环空两相流型第二节 压力分布计算方法第二节•由于多相管流中每相流体影响流动的物理参数（密度、粘度等）及混合物密度和流速都随压力和温度而变，沿程压力梯度并不是常数。
因此，多相管流压降需要分段计算，并要预先求得相应段的流体性质参数。
然而，这些参数又是压力和温度的函数，压力却又是计算中需要求得的未知数。
所以，多相管流通常采用迭代法进行计算。
一、常用两相流压降计算方法•早期均匀流方法（总摩阻系数法）1952 Poettmann 1952 Poettmann——Carpenter 80 80’’s 陈家琅 λ'~(N Re )2•经验相关式1963 Duns--Ros 无因次化处理 N vL 、Nvg 、N D 、N L 1965 Hagedorm--Brow 现场实验 1967 Orkiiszewski 流型组合 1973 Beggs--Brill 倾斜管实验1985 Mukherijee--Brill 改进实验条件•现代机理模型SPE20630等考虑具体流型的物理现象第二节第二节压力分布计算方法段塞流示意图环状流示意图2010-3-2682010-3-26当 单相液流， H L 、ρm 、f m 随两相流流型变化b P P ≥二、两相管流压降计算根据地面条件应用关系式计算井底流压1 输入数据油管数据：管长L 、管径D 、井斜角θ、粗糙度e 油气井产量：油气水日产量Q O 、Q SC 、Q W或Q L 、f w 、GOR P (GLR P ) QQ W = f w Q L Q O = Q L - Q wQ SC = GOR P Q O 或 Q SC = GLR P Q L 边界条件：井口压力P wh 、井口温度T wh 、地温梯度g t 考虑井温线性分布 T(Z)=T wh +g t Z 油气水相对密度γo 、γg、γw第二节第二节压力分布计算方法2 输入数据单位处理常用单位 统一单位Q —m 3/d q —m 3/s μ—Pa.s P —MPap —Pa V Pa V——m/s D —mm d mm d——m T m T——℃ T T——K 第二节第二节压力分布计算方法3 输入流体物性资料气：拟临界压力、温度 Pc ， Tc偏差系数 Zg(Pr, Tr)粘度μg油：μo ， 溶解油气比 Rs体积系数 Bo ， 油气界面张力σo水：μw ，σw ，B w第二节第二节压力分布计算方法),(223004hk p h z F k ++=第二节第二节压力分布计算方法龙格库塔数值解法•压力梯度函数F(Z，P)计算步骤(1) Z处流动温度 T(Z)=T0+g t Z(2) 计算T、P条件下的有关物性(3) 气液体积流量 q g，q L(4) 气液表观流速V sg、V sL和V m(5) 计算λL、μL、ρns、μns(6) 无因次量N Rens、N L、N gV、N LV、N gvsM(7) 计算H L、ρm(8) 判别流型，计算f m(9) 计算F(Z，P)•2、迭代计算第二节第二节压力分布计算方法误差又能提高计算速度。
垂直管流实验
垂直管流实验一、实验目的1.观察垂直井筒中出现的各种流型，掌握流型判别方法；2.验证垂直井筒多相管流压力分布计算模型；3.了解自喷及气举采油的举升原理。
二、实验原理在许多情况下，当油井的井口压力高于原油饱和压力时，井筒内流动着的是单相液体。
当自喷井的井底压力低于饱和压力时，则整个油管内部都是气-液两相流动。
油井生产系统的总压降大部分是用来克服混合物在油管中流动时的重力和摩擦损失，只有当气液两相的流速很高时（如环雾流型），才考虑动能损失。
在垂直井筒中，井底压力大部分消耗在克服液柱重力上。
在水平井水平段，重力损失也可以忽略。
所以，总压降的通式为：式中：—重力压降；—摩擦压降；—加速压降。
在流动过程中，混合物密度和摩擦力随着气-液体积比、流速及混合物流型而变化。
油井中可能出现的流型自下而上依次为：纯油流、泡流、段塞流、环流和雾流。
除某些高产量凝析气井和含水气井外，一般油井都不会出现环流和雾流。
本实验以空气和水作为实验介质，用阀门控制井筒中的气、水比例并通过仪表测取相应的流量和压力数据，同时可以从透明的有机玻璃管中观察相应的流型。
三、实验设备及材料仪器与设备：自喷井模拟器，空气压缩机，离心泵，秒表等；实验介质：空气，水。
设备的流程（如图1所示）图 1 垂直管流实验设备流程图四、实验步骤1.检查自喷井模拟器的阀门开关状态，保证所有阀门都关闭，检查稳压罐的液位（3/4液位）；2.打开空气压缩机及供气阀门；3.打开离心泵向系统供液；4.打开液路总阀，向稳压罐中供液，控制稳压罐减压阀，保证罐内压力不超过0.12MPa ；5.待液面达到罐体3/4高度，关闭液路总阀，轻轻打开气路总阀和气路旁通阀，向实验管路供气，保证气路压力不大于0.5MPa ，稳压罐压力约为0.8MPa；6.轻轻打开液路旁通阀，向系统供液，待液面上升至井口时，可以改变气液阀门的相对大小，观察井筒中出现的各种流型；7.慢慢打开液路测试阀门和气路测试阀门，然后关闭气路旁通阀和液路旁通阀，调节到所需流型，待流型稳定后开始测量；8.按下流量积算仪清零按钮，同时启动秒表计时，观察井底流压和气体浮子流量计的示数。
两相流计算方法l
蓝色流体网
液-固两相流
蓝色流体网
S = vg vs
S = vs vl
5. 压力降
混合物的两相流压力降 ∆pT 与气相、液相、固相分相压力降 ∆pg 、∆pl 、 ∆ps 之间的关系有： A. Murdock J.W.公式；林宗虎公式
蓝色流体网
蓝色流体网
2) 应用近代新技术； 3) 应用基于软测量技术的软测量方法。 2. 主要应用 相分率：
采用低能γ 射线传感器测量相分率：单能γ 射线传感器测量含气 率，双能γ 射线传感器测量相含水率； 流量：
采用文丘里流量计或文丘里流量计+互相关技术测量总流量； 根据相分率和数学模型确定油、气、水各相的流量； 压力、温度变送器
蓝色流体网
大部分旋风分离器的切割粒径50之间实际上来流中携带的绝大多数液滴直径要比切割粒径大的多实际上工业用蓝蓝蓝色流体网色流体网色流体网色流体网wwwopenfluidcn蓝色流体网蓝色流体网蓝色流体网蓝色流体网wwwopenfluidcn旋风分离器可以除去大多数液体设计中的主要问题不在于分离器能否将液体分离出来而是如何合理地处理已经被离心力甩向边壁的液大部分气液旋流分离器不采用锥体结构而采用圆筒型结构
蓝色流体网
蓝色流体网
旋风分离器可以除去大多数液体，设计中的主要问题不在于分离器能 否将液体分离出来，而是如何合理地处理已经被离心力甩向边壁的液 体相。
大部分气-液旋流分离器不采用锥体结构，而采用圆筒型结构。 分离器的容量应保证旋涡的尾部不与下部的液面接触。因此隔离板 （或消涡板或稳涡板）的作用是提供一个接触面，旋涡的末端在此面 上像陀螺一样不停地旋转，该板的目的不是破坏和干扰旋涡，而是避 免旋涡与液面接触，可以大大降低湍流和明显削弱旋涡强度。
第4章-管内气液两相流的阻力计算
平均
PF 值。一般认为
此方法适用于低质量流速工况。。
22
4.2摩擦压降计算
1 -摩擦压降计算的经验方法
分相流模型
2 奇斯霍姆对两相流动摩擦阻力压降进行了许多研究，得出了 L 与X PF C 1 2 值的关系为 L 1 2
或者
PL0
X
X
2 FG = 1 + CC + C 2
G0 Ct (
VSG D n ) VG
（9）
14
（光滑管区的Blasius公式为： 0.3164 ） Re0.25
4.2摩擦压降计算
1 -摩擦压降计算的经验方法
AG D ( G )2 ,VSG VG 则（ A D
n5 2 G 2
分相流模型
又因
所以，
m=
1 x
PF =
PF =
L V2 D 2
2 L mVm D 2
（4-1） （4-2） （4-4）
2 （mVm） L PF = [1 x( L 1)] Dቤተ መጻሕፍቲ ባይዱ2 L G
G
+
1-x
L
（x为平均干度）（4-3）
苏联50年代锅炉水循环计 算法
该式计算误差对水平均相 为（>20%）
该方法适用于双组分的气液两相流在低压时摩擦阻力计算，因其计算数据是建
立在低压的气液流动基础之上的。
为了适用于汽—水混合物的摩擦阻力压力降的计算，Martinelli-Nelson对此方 法进行了改进。
19
4.2摩擦压降计算
1 -摩擦压降计算的经验方法
分相流模型
假设汽水两相分开流动时都呈紊流状态，同时利用常压下的空气 —水混合物
气液两相管流计算
气液两相管流计算1基本要求.从能量平衡微分方程出发，与压降等经验关联式联立，得出了油气集输管线沿程温降的计算公式，并编制大型计算程序。
2 数学模型2.1 热力计算能量平衡方程假设两相之间没有温度滑移, 并不计油品的径向温度梯度，这样气液两相混合物沿管线的能量微分方程可写成:dx dv vg dxdH dx dq ++=θsin (1) q —垂直于管壁方向的热流量；H —混合物焓；v —混合物平均速度; θ—管轴线与水平面夹角；g —重力加速度。
由于混合物焓H 依赖于它自身的压力P 和温度T f , 这样dHdx 可下式来表示： dH dxH P dP dx H T dT dx C dPdx C dT dx Tf f P f J pm Pmf =+=-+()()∂∂∂∂μ (2)其中μJ 为焦耳-汤普森(Joule-Thompson)系数, 其物理意义是流体每单位压力变化引起的温度变化，C Pm 为混合物定压比热。
由(1)式和(2)式可得:dx dPdx dv v g dx dq C dx dP dxdH C dx dT JPm J pm fμθμ+--=+=)sin (11 (3) 式中单位长度热流量dqdx r k T T o f s =--2π() (4)负号表示散热，T s 为环境土壤温度， k 为传热系数。
将(4)式带入(3)式，整理得dT dx T T A Bf s f=-+ (5) 式中A C r k Pm o =2π ,B dP dx gC v C dv dx J Pm Pm =--μθsin 。
式(5)为一阶线性微分方程, 对应的边界条件为:在出油管口处，当 x =0时 T f =T f0解方程(5)可得温度分布:T T T k k C x k T f f s Pm s=---+-()exp()012212πππΦΦ (6)式中Φ=--C dP dx g vdvdx Pm J μθsin 式(6)中除传热系数k 和压降dP/dx 外，其余参数一般均为常数。
管道压力降计算--气液两相
注意事项 1 适用于气液混合物中 气相在6% ～98%（体 积）范围内
2 在工程设 计中，一般 要求两相流 的流型为分 散流或环状 流；若选用 的管路经计 算后为柱状 流，应在压 力降允许的 情况下
气-液两 相流（非 闪蒸型）
uc PucM /(RT )
m
2
1 3
(
1
2)
M:
16 Kg/mol
ρ1:
2.842859335 Kg/m³
ρ2:
0.94977346 Kg/m³
ρm:
1.580802085 Kg/m³
基本与不可 压缩流体计 算方法一致
m
2
1 3
(
1
2)
（设定值）
PM / RT
控制误差 <1%
由Re准数 判断流型 输入： u: d:
μ:
w: Vf: ρ:
通过计算 P1：
选定P1:
432.8196316 Kpa 433 Kpa
等温流动声 速uc 输入：
T: M： 输出： uc：
298 K 16 Kg/mol
393.412951 m/s
声速下的临 界流量 输入：
uc: d: 输出：
Vuc:
393.41 m/s 0.307 m
104783.8384 m³/h
0 m/s 0 m/s 930 Kg/m³
0 Kpa
总压力降Δ P
ΔP:
259.2714412 Kpa
单相流 （可压缩 流体）
注意事项
1 压力较 低，压力降 较小的气体 管道，按等 温流动一般 计算式或不 可压缩流体 流动公式计 算，计算时 密度用平均 密度，对高 压
两相流计算方法l
蓝色流体网
旋风分离器可以除去大多数液体，设计中的主要问题不在于分离器能 否将液体分离出来，而是如何合理地处理已经被离心力甩向边壁的液 体相。
大部分气-液旋流分离器不采用锥体结构，而采用圆筒型结构。 分离器的容量应保证旋涡的尾部不与下部的液面接触。因此隔离板 （或消涡板或稳涡板）的作用是提供一个接触面，旋涡的末端在此面 上像陀螺一样不停地旋转，该板的目的不是破坏和干扰旋涡，而是避 免旋涡与液面接触，可以大大降低湍流和明显削弱旋涡强度。
vm A = vg Ag + vl AL
气-固两相流
vm A = vg Ag + vs As
液-固两相流
蓝色流体网
蓝色流体网
vm A = vs As + vl Al
工程中分相流速采用表观流速进行折算，即以分相流量除以管道
总截面的比值表示该相的分相速度。物理意义是当管道内的流体全是
防涡器一般安装在液体出口管的正上方。因为进入分离器的气液 两相混合物的旋转角动量使得储液池中的液体随之旋转。如果允许液 面上旋涡的存在，则来流中一部分气体就会在底流口排出，会对下游 泵产生气蚀或其他问题。这种旋涡还会产生射流阻塞效应，限制低部 出口管的液体流量。
防涡器通常为十字形金属平板或环形平板组成，位于排出管上方 一倍出口管直径处。
积之比。
气-液两相流：容积含气率
α = Vg = Vg （空隙率）
V Vg + Vl
气-固两相流：容积含固率
α = Vs = Vs V Vs + Vg
液-固两相流：容积含固率
α = Vs = Vs V Vs + Vl
D. 局部含率
如果把两相流体所占总体积 V 取得足够小，则在该微小体积上