在化工、环保、能源、制药等众多工业领域中,气相传质系数是衡量气体组分从气相向液相转移效率的核心参数。它不仅影响吸收塔、解吸塔、洗涤器等设备的设计与运行,更直接关联到产品收率、能耗与排放控制。本文以气相传质系数公式为主线,系统梳理其理论框架、经典模型、典型应用与常见误区,结合真实工业案例与工程经验,为技术人员提供一套可落地、可验证、可优化的实践指南。
立即探索原理体系理解气相传质过程,需跳出直觉,从分子运动与界面动力学切入
气相传质系数(通常记为 kG 或 ky)定义为:单位推动力(如气相分压差或摩尔分数差)下,单位面积上气体组分的传质速率。其物理本质是反映气体分子穿越气膜(或主流区)向相界面扩散能力的强弱。数值越大,说明传质阻力越小,效率越高。
年,Lewis与Whitman提出的双膜理论是理解气液传质的经典模型。该理论假设:在气液界面两侧存在两个静止的气膜与液膜,传质阻力集中于膜内,且界面处无传质阻力、已达平衡。虽然实际流动中不存在真正静止的膜,但该模型因其物理图像清晰、数学处理简便,至今仍是工程估算与设备设计的基石。
其中:
• NA:组分A的传质通量(mol·m⁻²·s⁻¹)
• kG:以分压为推动力的气相传质系数(mol·m⁻²·s⁻¹·Pa⁻¹)
• pA:气相中A的分压
• pA:与液相浓度平衡的气相分压(由亨利定律确定)
气体分子密度低、扩散系数大(约10⁻⁵ m²/s),看似传质快;但实际中,当气体溶解度极低(如O₂、N₂在水中的溶解度仅为10⁻³量级)或反应速率慢时,气膜阻力常占总阻力的70%以上。此时,即使强化搅拌或提高液相浓度,效果也有限——必须从气相传质系数本身入手,如增大气液接触面积、增强湍流程度、优化分布器设计。
在自然通风冷却塔中,水蒸气从液相进入气相,看似是“蒸发”,实则属于气相传质(水蒸气从界面扩散至主流空气)。由于水蒸气在空气中扩散系数大(约2.6×10⁻⁵ m²/s),且气相湍流强烈,实际气相传质系数可达0.02~0.05 m/s量级,远高于低溶解度气体(如CO₂在水中kG约10⁻⁵~10⁻⁴ m/s)。因此,只要湿球温度允许,蒸发速率主要受气相传质控制——这正是冷却塔高效运行的关键物理机制。
不同工况需匹配不同模型,选错公式 = 工程事故
Hatta于1927年提出该模型,用于处理快反应吸收过程——即气体一旦接触液相,瞬间即被反应消耗。此时,气相传质系数需结合反应动力学进行修正。
• Ha:哈塞尔巴尔赫数(无量纲)
• kG:气相传质系数
• CA0:气液界面处气体浓度
• DL:组分在液相中的扩散系数
• kr:二级反应速率常数
当Ha > 2时,反应速率远大于扩散速率,反应“前沿”紧贴界面,传质被气膜控制;当Ha < 0.02时,反应极慢,可视为物理吸收;当0.02 < Ha < 2时,需引入“反应增强因子”修正。
CO₂与NaOH反应属快反应(kr≈10⁴ L·mol⁻¹·s⁻¹),按Hatta数计算,当[NaOH] > 2 mol/L时,Ha > 5,反应前沿极薄。此时若继续提高碱浓度,反应速率已不再提升(受限于CO₂扩散速率),反而因溶液黏度增大、扩散系数下降,导致气相传质系数降低,整体吸收速率不增反减——这正是工业中为何常采用1~1.5 mol/L NaOH的原因。
经典双膜理论忽略界面湍动、表面张力梯度(Marangoni效应)及界面更新速率。实际中,常引入“表面更新模型”思想进行修正:将气相传质系数表示为:
• DG:气体在气相中的扩散系数
• f:表面更新频率(Hz)
• tc:接触时间(s)
该式表明:提高液体湍流程度(增大f)比单纯增加气体流速更有效——因为f与搅拌功率的0.4~0.6次方成正比,而气体流速仅线性影响kG。
规整填料(如Mellapak)通过规则通道形成薄液膜,大幅提高界面更新频率f。实验测得其气相传质系数比拉西环高30%~60%,且压降更低。这正是双膜理论修正模型的直接工程体现——通过结构设计优化表面更新动力学。
Danckwerts指出:液相表面分子不断被新分子替换,传质是动态过程。他提出表面更新模型,其中气相传质系数可表示为:
• S:表面更新频率(m/s)
• DG:气相扩散系数(m²/s)
该模型强调:传质速率与扩散系数的平方根成正比,且与表面更新频率的平方根成正比。这意味着,仅靠提高气体流速(线性提升S)效果有限,而通过填料结构诱导涡流、增强界面扰动(非线性提升S)更有效。
| 设备类型 | kG (m/s) | S (m/s) |
|---|---|---|
| 鼓泡塔(无搅拌) | 0.001~0.003 | 0.01~0.05 |
| 搅拌釜(P=1 kW/m³) | 0.005~0.015 | 0.1~0.5 |
| 规整填料塔 | 0.01~0.04 | 0.3~1.2 |
数据来源:Perry's Chemical Engineers' Handbook, 9th Ed., Table 13-3
Higbie假设液相表面短暂暴露于气相,接触时间为tc,则气相传质系数为:
该模型适用于喷雾塔、文丘里洗涤器等气体与液体短时接触的设备。例如,在喷雾塔中,液滴停留时间约0.1~1秒,计算得kG约0.002~0.006 m/s,与实测值吻合良好。
已知:气体流速30 m/s,喉管直径0.2 m,液气比0.005 L/m³,气体为含SO₂空气(DG=1.6×10⁻⁵ m²/s)。根据经验公式:
tc ≈ 0.1 s,代入得:
kG = √(1.6×10⁻⁵ / (π×0.1)) ≈ 0.0071 m/s
实测值约0.006~0.008 m/s,误差<15%——证明该模型在高湍流、短接触时间工况下具有工程实用性。
脱离案例的公式是纸上谈兵,脱离公式的案例是经验主义
过程描述:SO₂氧化为SO₃后,用98.3% H₂SO₄吸收SO₃生成发烟硫酸。反应为:
SO₃(g) + H₂O(l) → H₂SO₄(l) (瞬间完成)
问题:传统设计中,若按纯物理吸收计算气相传质系数,会高估吸收效率,导致尾气SO₃超标。
解决方案:引入Hatta数校正。实测[ H₂SO₄ ]=98.3%时,Ha≈8.2,需将kG乘以增强因子β≈1.8。据此优化填料层高度后,吸收率从99.2%提升至99.95%,年增效益超200万元/万吨酸。
过程描述:在高压(0.2~0.3 MPa)下,NH₃与CO₂在饱和盐水中反应生成NH₄HCO₃,再与NaCl反应析出NaHCO₃。其中NH₃吸收是限速步。
误区:工程师常误认为高压直接提高溶解度,无需优化气相传质系数。
真相:根据亨利定律,NH₃溶解度虽随压力上升,但气相传质系数kG与压力平方根成反比(因气体密度↑、扩散系数↓)。计算得:压力从0.1 MPa升至0.25 MPa时,kG下降约22%!因此,必须同步强化气相湍流(如改进分布器、增加挡板),才能抵消kG下降的负面影响。
背景:SCR脱硝反应为:
4NO + 4NH₃ + O₂ → 4N₂ + 6H₂O
低温(<150℃)时,反应速率慢,NH₃易逃逸。
优化路径:通过CFD模拟发现,低温下烟气黏度增大,导致气相传质阻力上升。实测kG下降约35%。对策:在催化剂前增设湍流促进器,使kG恢复至设计值,NH₃逃逸率从8 ppm降至3 ppm以下。
| 设备类型 | kG (m/s) | 压降 (Pa/m) | 适用场景 |
|---|---|---|---|
| 鼓泡塔 | 0.001~0.003 | 50~200 | 慢反应、高液气比 |
| 填料塔(散堆) | 0.005~0.012 | 80~300 | 中等反应、清洁体系 |
| 填料塔(规整) | 0.01~0.04 | 30~150 | 快反应、高效率要求 |
| 喷雾塔 | 0.002~0.008 | 100~500 | 高粉尘、急冷急热 |
注:工业现场多采用动态吸收法+经验公式校正,误差约±15%,已满足设计需求。
忽略任一因素,都可能让精心设计的系统“失灵”
温度升高 → 气体扩散系数DG↑(按Stokes-Einstein关系,DG ∝ T1.75/P)→ kG↑
但同时 → 溶解度↓(亨利常数↑)→ 推动力↓ → 传质通量可能下降!
例:25℃→50℃,T从298K→323K,T1.75增长约14%,故kG理论提升14%;但CO₂在水中的溶解度下降约30%,导致整体吸收率反而降低——这正是碳酸饮料在夏天“更易跑气”的物理本质。
增压可提高气体分压(↑推动力),但:
• 气体密度↑ → 扩散系数DG↓(DG ∝ 1/P)→ kG↓
• 液体黏度↑ → 传质阻力↑
因此,存在一个最优操作压力。以合成氨原料气净化(脱CO₂)为例:压力从2 MPa升至4 MPa,吸收率仅增3%,但压缩功耗翻倍,经济性恶化。
| 气体 | DG (10⁻⁵ m²/s) | 25℃水中的溶解度 (mol/m³·Pa) | 相对kG(同设备下) |
|---|---|---|---|
| He | 6.2 | 3.8×10⁻⁴ | 1.8 |
| H₂ | 4.5 | 7.8×10⁻⁴ | 1.2 |
| N₂ | 2.0 | 6.5×10⁻⁴ | 1.0 |
| CO₂ | 1.6 | 3.3×10⁻² | 0.8 |
| NH₃ | 1.4 | 0.55 | 0.3 |
注:NH₃虽扩散慢(DG小),但因强极性、高反应性,实际吸收极快,整体传质速率高。这再次说明:气相传质系数仅反映气相阻力,需结合液相传质系数(kL)和平衡关系综合判断。
空分装置中,分子筛吸附器需同时去除H₂O、CO₂、C₂H₂。若进气温度骤降(如暴雨天),气体中水蒸气过饱和,瞬间凝结成液滴。此时:
• 液滴堵塞通道 → 气体流速局部↑ → 湍流↑ → kG短期↑
• 但液滴覆盖吸附剂表面 → 有效吸附面积↓ → 实际传质速率↓↓
更严重的是,水在分子筛孔道内结冰(-40℃环境),导致永久性孔道堵塞。因此,设计时必须考虑:
1、进气温度控制(≥10℃)
2、设置预冷+气液分离器
3、选用疏水性改性分子筛(如Si/Al比>25的X型)
这些措施,本质上都是为了保障气相传质系数在长时间运行中稳定。
这些问题,90%的新人会问,80%的老师傅也答不全
这是最常见的“理论-实践偏差”。原因通常有:
建议:现场测试时,务必用多点取样法测量进出口浓度分布,结合CFD验证流场均匀性。
可以,但需注意:
这取决于搅拌类型与体系性质:
经验公式(搅拌釜):
kGa = 0.025·(Ng)0.7·(P/V)0.4·(μL/ρL)-0.2
其中P/V为单位体积功率(W/m³),μL为黏度(Pa·s)。
这是单位制与推动力选择导致的:
• kG:以分压为推动力,单位 mol·m⁻²·s⁻¹·Pa⁻¹
• ky:以摩尔分数为推动力,单位 mol·m⁻²·s⁻¹(注意:此时推动力是y−y,无量纲)
二者关系:
ky = kG·P
(P为总压,单位Pa)
换算示例:若kG=1.0×10⁻⁵ mol·m⁻²·s⁻¹·Pa⁻¹,P=101325 Pa,则ky=1.013×10⁻⁰ ≈ 1.0 mol·m⁻²·s⁻¹(数值巧合,但方法通用)。
建议:统一使用kG,避免混淆;或明确标注单位与推动力定义。