气液两相流流型预测与段塞流治理:流型图、Python模型与工程实践
2026/9/18 15:23:27 网站建设 项目流程

简介:这是一份面向两相流研究学习者的流型专题PPT课件,聚焦教材第3章内容,系统介绍气液两相流中流型的基本概念、分类方法、识别技术及工程意义。课件从Baker等人上世纪50年代的研究切入,梳理了由实验观察、流型图绘制到建立数学物理模型的方法演进,并具体对比了基于两相外形分布与基于数学模型的两种分类体系,涵盖泡状流、柱塞流、段塞流、环状流、分层流等常见流型,同时分析了高速摄影、射线衰减、电容法、压降脉动分析及电导法等识别手段的适用性与局限,帮助读者理解流型转变的复杂边界与过渡状态。压缩包内含1个PPTX文件,整体大小5.1MB,页面以要点插图形式编排,便于课堂教学展示与课后自学梳理。目前已有134人学习下载,适合能源动力、化学工程、石油储运等专业的学生、教师及工程技术人员作为两相流基础知识的配套参考资料。

1. 流型不是教科书名词:气液两相流的形态直接决定管道振幅和仪表读数

一条水平输油管道的某个坡度段持续震动,压差变送器读数每隔十几秒规律性跳变,超声波流量计偶尔直接丢信号。换了仪表、加了减震垫都没有用,后来在透明管段上用高速摄影才看清:液体不是均匀裹着气体,而是聚成一个个液塞,周期性地撞向弯头。

这个场景就是气液两相流中的流型(flow pattern)在起作用。流型指气液两相在管道内流动时的空间分布结构,它决定计算压降该选哪种模型、流量计能不能测准、管壁会不会被冲蚀,也是多相流仿真里每一项边界条件的起点。无论是工艺设计、仪表选型,还是正在准备流型课件的人,先把分类和判别方法吃透,后面所有结论才立得住。下面的内容就从流型分类讲到可运行的预测代码,再落到工程验证上。

2. 水平管与垂直管的流型分类:从泡状流到环状流,边界条件是工程判据

2.1 六种基本流型的形态与工程风险

按空间构型划分流型是最常用、也最容易被误解的做法。水平管与垂直管的流型差异很大,同一组气液流量在水平管里可能是层状流,换成垂直管就变成弹状流,所以先明确管子的空间方位再谈分类。

水平管中常见六种流型:泡状流、塞状流、段塞流、层状流、波状流和环状流,气量继续增大还会进入雾状流。泡状流的气泡直径远小于管径,常见于高压下气相体积占比低、或液相流量很高的工况;段塞流则是液体周期性堵塞整个管道截面,后方跟着一个拉长的气相区,这是工程上最需要警惕的流型。垂直管更常见泡状流、弹状流、搅动流和环状流,弹状流的特征是泰勒气泡(Taylor bubble)接近管径、把液相分成一段一段。

流型常见工况形态特征主要工程风险
泡状流高液量、高压、低气量小气泡分散在液相中风险较低,压降可用均相模型近似
段塞流/弹状流中等气液量、水平或微倾管液塞堵满管截面,周期性冲击机械振动、仪表读数跳变、分离器液位波动
层状流/波状流低气液量、水平或近水平气液分层,界面平直或有波浪底部腐蚀集中,液位界面不稳定
环状流高气量、垂直或水平液膜沿壁面环形流动,高速气相在中心液膜夹带、弯头冲蚀、传热不稳定
雾状流超高气量液滴悬浮在连续气相中湿气测量偏差大,高速液滴冲蚀

表里的风险不是绝对的。层状流在大直径管道中也可能过渡到段塞流,完全取决于气液相表观速度落在哪个区间。这里先建立一个基础概念:表观速度就是把两相分别按管道全截面计算的速度,符号是 jl 和 jg。两相各自的实际流速是未知的,因为截面含气率还没算出来,但体积流量和管截面积是设计值,可以先算作流型图坐标,这个顺序不要搞反。

2.2 流型图是把经验变成坐标:Mandhane、Baker 与 Weisman

上世纪七十年代到八十年代,研究者把大量实验数据画成了流型图,最常用的是 Mandhane 图。横坐标是气相表观速度 Vsg,纵坐标是液相表观速度 Vsl,图里用边界线划分出泡状、段塞、层状、波状、环状区域。它直观,绘图时算出两点坐标再落点就行。缺点是数据点来自常压空气-水系统,液相黏度、气液密度比偏离较多时,边界会明显移动。

所以工程上还有另外两种选择。Baker 图改用气液相质量流速做坐标,并加入了液相表面张力、气液密度相关的物性修正系数,对油品伴生气等石油化工介质更友好;Weisman 则把密度比、黏度比、表面张力整理成无量纲组合后再做坐标,适用范围比 Mandhane 宽。垂直管常用基于动量通量 G²/ρ 的 Hewitt-Roberts 图,判断逻辑类似,但坐标系换成动量通量后,高低压工况的迁移性更好。

提示:流型图本质上都是经验图。用 Mandhane 图做初步判断没有问题,但正式的管道设计不能只靠它,要结合物性修正或商用瞬态软件做段塞频率分析。如果手头只有一张课件里的流型图,先问一句它是哪个系统的数据。

鉴别流型图适用范围有一个快速方法:看横纵坐标是不是体积通量或质量通量。凡是直接用 Vsg/Vsl 的图,大多来自空气-水实验;标了物性修正系数的图,才考虑了不同介质的迁移。课件里常见把两种图混用,这是最容易出错的地方。

2.3 流型转变的物理边界:界面失稳、漂移通量与截面含气率

流型不是随机出现的,它由气液界面上的受力平衡决定。水平管里,层状界面上气相速度增加,波动振幅会增大,当气相对于波峰的动压超过重力恢复力,界面失稳进入波状流。波状流的波峰再增高,波峰上方气速增大、静压降低,伯努利效应会把波峰“吸”起来,最终堵塞管道成为段塞流,这是 Kelvin-Helmholtz 失稳在管道流里的经典表现。

垂直管则是另一种机制。小气泡在液相里以相对速度 Vgj 上升,气量增大时气泡浓度上升,合并成接近管径的泰勒气泡,形成弹状流;气量进一步增大,泰勒气泡之间的液相段被搅碎,进入搅动流;气相变成连续相、液相变成液膜,就过渡到环状流。无论水平还是垂直,核心都是气液两相之间的相对速度。

相对速度的工程表达常用漂移通量模型:

vg = C0 * j + Vgj

其中 j 是气液混合物总体积通量,等于 jl 加 jg。C0 是分布系数,圆管泡状流常见取值约 1.2;Vgj 是气相在停滞液体中的漂移速度,垂直管中单个气泡约 0.25 m/s,泰勒气泡可按 0.35*sqrt(gD) 估算,水平管中浮力方向与流动方向垂直,漂移速度的概念弱一些,更多由界面剪切主导。

用漂移通量模型可以反推截面含气率 α = jg / (C0*j + Vgj),它是两相流压降模型的前提,也是后面代码里交叉验证的关键参数。这里先记住结论:α 接近 1 时气相已经连续,泡状流、段塞流这类“液相连续”的判断基本不成立。

3. 用 Python 实现流型预测:表观速度计算与漂移通量模型

流型图看着简单,真要在项目里反复算,最省力的方式是把物性、流量和几何参数写成一段可复用代码,输入工况、输出流型区间和截面含气率。下面给出一个最小实现,跑通后再按自己的介质改参数。

3.1 最小代码:把质量流量转成表观速度并落在流型区间

import math # ---- 工况输入参数 ---- D = 0.05 # 管道内径 m(不是公称直径 DN) rho_l = 998.0 # 液相密度 kg/m3 rho_g = 1.2 # 气相密度 kg/m3(工况压力下的密度) m_l = 2.0 # 液相质量流量 kg/s m_g = 0.02 # 气相质量流量 kg/s # ---- 计算表观速度 ---- A = math.pi * D**2 / 4 j_l = (m_l / rho_l) / A j_g = (m_g / rho_g) / A j_m = j_l + j_g # ---- 简化水平管流型判据(空气-水,常压附近)---- if j_g > 3.0: flow = "环状流 annular" elif j_l > 0.5 and j_g < 1.0: flow = "泡状流 bubbly" elif j_g < 0.3 and j_l < 0.01: flow = "层状流 stratified" elif j_g < 0.3 and j_l < 0.1: flow = "波状流 wavy" else: flow = "段塞流 slug" print(f"液相表观速度 jl = {j_l:.3f} m/s") print(f"气相表观速度 jg = {j_g:.3f} m/s") print(f"混合物表观速度 j = {j_m:.3f} m/s") print(f"预测流型:{flow}")

这段代码做了三件事。第一,由质量流量除以密度得到体积流量,再除以管道截面积,得到表观速度,这是流型图的横纵坐标。第二,用几个边界区间简化 Mandhane 图,得到初步流型。第三,把 j_m 一并打印,供下一步漂移通量计算使用。

简化判据的写法上,我把高气速的环状流放最前面,因为气相表观速度超过 3 m/s 后优先怀疑环状流,再根据液量大小核实液膜是否夹带;段塞判据放最后,覆盖中间的过渡区间。这套边界对空气-水常压系统能给出合理结果,换成油品或高压天然气就必须换边界,生产环境建议直接查原始流型图数据或加物性修正。

3.2 漂移通量模型:用截面含气率交叉验证流型判断

import math D = 0.05 g = 9.81 # 垂直管段塞流中泰勒气泡的漂移速度估算 Vgj = 0.35 * math.sqrt(g * D) # 约 0.245 m/s C0 = 1.2 alpha = j_g / (C0 * j_m + Vgj) alpha = min(alpha, 0.95) # 限制物理上限 print(f"漂移通量模型截面含气率 alpha = {alpha:.3f}") if alpha > 0.7: print("提示:含气率偏高,气相趋于连续,核实是否应为环状流/雾状流")

这里直接使用上一段算出的 j_g 和 j_m。C0 取 1.2 是圆管湍流泡状流的常用分布系数;Vgj 用 0.35 乘以根号下 gD 估算垂直管泰勒气泡的漂移速度,水平管计算时可以把 Vgj 调到 0.1 m/s 量级。alpha 算出来是截面含气率,即气体占截面的平均份额。如果 alpha 超过 0.7 而流型判据给出的却是泡状流或段塞流,两个结论互相矛盾,优先怀疑判据或物性输入有误。

提示:漂移通量模型假设气泡均匀分布在液相中并整体向上漂移,更适合垂直管。水平管需要引入更复杂的 Taitel-Dukler 判据,代码量会大一个量级,工程上先用模型交叉验证再决定要不要走全瞬态仿真。

3.3 输入参数与单位一致性检查清单

参数符号单位示例值常见错误来源
管道内径Dm0.050用 DN 公称直径代替内径,小管径误差可达 5%-10%
液相密度rho_lkg/m3998没按工况温度查密度
气相密度rho_gkg/m31.2把标方密度直接当工况密度,高压下差 10 倍以上
液相质量流量m_lkg/s2.0从体积流量换算时没乘操作密度
气相质量流量m_gkg/s0.02标方流量没换算到工况
液相表观速度j_lm/s1.02用实际流速代替表观速度
气相表观速度j_gm/s8.49同上

单位一致性是流型预测里最常翻车的点。化工装置里气相流量经常给标方每小时,必须先换算成工况压力、温度下的实际体积流量,再除以管截面积。不进行工况换算直接用标方,算出的表观速度会大到一两个数量级,流型图直接落在环状流或雾状流区,后面的设计动作全部失真。拿到物性表先确认密度是不是在 P、T 工况点上,这是跑代码之前就该做好的事。

4. 流型结果用在管道设计、仪表选型与段塞流缓解上

流型判断本身不是终点。实际项目中拿到一个“段塞流”结论后,马上要回答两个问题:现场的仪表还能不能用,管道结构扛不扛得住。下面按仪表、管线分段、缓解措施三个层次讲。

4.1 段塞流对差压式与涡街流量计的典型影响

差压式流量计的流量模型建立在稳态流动假设上,差压与流量之间是平方根关系。段塞流进入测量管段时,差压信号随液塞和气泡交替出现大幅波动,瞬时差压波动很容易超过满量程的 50%。简单加阻尼滤波会增大响应迟滞,液塞周期变化慢时读数仍然跟随波动,体积流量积分误差不可忽略。

涡街流量计的问题在低流速间隙。段塞与段塞之间是低速液膜区,流速低于旋涡脱落的最低雷诺数条件时,涡街信号直接消失,输出掉零,形成周期性的锯齿波。我遇到这类工况,第一反应不是调滤波参数,而是建议把流量计移到垂直管段,或者在仪表上游装整流器。

4.2 管线分段计算流程与流型判据

管路不可能整条用一个流型。倾角、管径、阀门和弯头都会改变局部气液表观速度。落地的做法是先把路由拆开再逐段判断:

  1. 按照倾角变化点、管径变化点和控制阀位置,把管线切成若干段。
  2. 每一段取入口压力估算该段的实际气液密度,特别注意气体密度随压降变化。
  3. 计算该段的气液相表观速度 jl、jg。
  4. 按水平或垂直选对应流型图或代码判据落点。
  5. 把段塞流、弹状流出现的段标记为高风险段,优先做管架、仪表测点复核。

分段不必过细,实际工程中 5 到 10 米一段已经足够。重点是不要只看管线的平均压力,高压段和低压段的气体体积流量可能差一倍,流型也就完全变了。管线压力从 6 MPa 降到 2 MPa 的过程里,气相体积膨胀约三倍,同样的液相流量下,流型完全可能从段塞流迁移到环状流。

判出流型设计动作
段塞流/弹状流设段塞捕集器、加固管架、仪表避开低点、增加防振支架
环状流/雾状流弯头背面加厚或加防冲板、控制气速上限、湿气流量计做修正
层状流/波状流低点增加排凝、底部包覆耐腐蚀材料、液位计选型核对
泡状流常规设计即可,注意气泡析出后流型迁移

4.3 缓解段塞流的结构设计与仪表参数调整

结构上的根本缓解办法是段塞捕集器。常见做法是来料先进入一个大直径缓冲段或一组指状管,把液塞的液相暂时存下来,让气相在顶部自然分离,再以平稳流量送往后续处理。段塞流对下游气液分离器液位的影响也很大,液塞进入分离器会让液位骤升触发高液位联锁,分离器设计时要把段塞体积计入持液容积。

仪表层面能做的有限。差压变送器阻尼时间一般按段塞周期的 3 到 5 倍设置,但这会让真实瞬态信号被平滑掉,只能作为临时手段。段塞频率低于 0.1 Hz 时,30 秒阻尼会产生明显相位滞后,控制回路容易不稳定。更可靠的做法是重新评估测量点位置:段塞流最平稳的位置通常出现在垂直上升段和分离器入口附近,水平低点处液塞聚集最明显,不是装差压表或涡街表的首选位置。

5. 验证流型预测的三种手段与三个常见误用

5.1 压差信号的频谱特征验证

现场没有高速摄影时,压差信号本身就能验证段塞流。段塞流有准周期特征,液塞通过测量段时差压升高,气泡段通过时差压回落,在频谱上表现为集中的主峰。我一般用 100 Hz 采样率采集差压或管道振动信号,取 10 分钟数据做 FFT:

import numpy as np # dp: 差压变送器同步采集的时间序列(numpy 数组) # fs: 采样率,单位 Hz dp = dp - np.mean(dp) # 去掉直流分量 spec = np.abs(np.fft.rfft(dp)) ** 2 # 功率谱 freq = np.fft.rfftfreq(len(dp), 1.0 / fs) peak_f = freq[np.argmax(spec[1:]) + 1] # 去掉直流后找主峰 print(f"差压信号主频:{peak_f:.3f} Hz")

如果主频落在 0.01 到 0.5 Hz 之间且频谱有清晰的单峰,说明周期波动确实存在,流型判断里的段塞流得到现场旁证。配合时域波形看,段塞流的差压峰值往往是快速上升、缓慢下降,和环状流的高频随机波动很容易区分。

5.2 透明管段、高速摄影与可视化检查

实验室和中试装置上的验证更直接:换一段透明管段,用 500 fps 以上的高速摄影拍 10 秒。段塞流的标志是液塞几乎堵满管截面、气液界面前后清晰、液塞后跟一个长气泡;层状流的标志是界面平直,波浪偶尔出现但波峰不会跨越管顶。

5.3 三个最容易翻车的误用

第一,用均相模型包打天下。均相模型把气液当成均匀混合物,在高气速的环状流或雾状流下误差不大,但在层状流、段塞流下会显著低估真实压降。判断流型不是走形式,它决定后续选哪种压降关联式。

第二,忽略倾角和管径迁移。水平管流型图用于 5° 以上的倾斜管时,层状流到段塞流的边界会明显移动;垂直管用水平图判断更是错得离谱。任何迁移都要先确认管子和流型图实验条件一致。

第三,物性取错。把常压空气密度用到 1 MPa 工况,jg 直接差出 10 倍,流型图落点完全变了。最后留一个交叉验证技巧:把漂移通量模型算出的截面含气率 α 与流型判据同时摆出来,α 大于 0.7 时液相不再是连续相,任何泡状流或段塞流的结论都站不住,先回头查气体密度和表观速度的单位。

本文还有配套的精品资源,点击获取

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询