Skip to content

平均吸收因子法

定义/概述

平均吸收因子法是工程上常用的多组分气体吸收简捷计算方法。通过定义吸收因子,建立吸收率、吸收因子与理论板数的关系。

吸收因子

吸收因子 AiA_i 定义为操作线斜率与平衡线斜率之比:

Ai=LKiV(9.1)A_i = \frac{L}{K_i V} \tag{9.1}

其中 LL 为液相流量,VV 为气相流量,KiK_i 为相平衡常数(由 Antoine 方程 → Psat/PP^{\text{sat}}/P 计算)。

物理意义:Ai>1A_i > 1 表示操作线斜率大于平衡线斜率,组分倾向于被液相吸收;Ai<1A_i < 1 表示组分倾向于留在气相。关键组分的吸收因子通常在 1~2 之间。

Horton-Franklin 方程(通用多板方程)

Kremser 方程的推导起点是 Horton-Franklin 方程,它描述了多板逆流吸收塔中任意组分在任意板的物料衡算关系。

对于 NN 板吸收塔,组分 ii 在塔顶(板 0)和塔底(板 NN)之间的总物料衡算可表达为:

vN+1v1vN+1=A1A2AN+A2A3AN++ANA1A2AN+A2A3AN++AN+1(9.2)\frac{v_{N+1} - v_1}{v_{N+1}} = \frac{A_1 A_2 \cdots A_N + A_2 A_3 \cdots A_N + \cdots + A_N}{A_1 A_2 \cdots A_N + A_2 A_3 \cdots A_N + \cdots + A_N + 1} \tag{9.2}

其中 A1,A2,,ANA_1, A_2, \ldots, A_N 为各板的吸收因子。此方程是精确形式,但计算繁琐——需要知道每块板上的 L/VL/VKiK_i 才能求各板的 AA 值。

Kremser 方程(平均吸收因子简化)

当各板吸收因子取平均值 AA 时(即假设全塔 L/VL/VKK 恒定),Horton-Franklin 方程简化为 Kremser 方程:

vN+1v1vN+1v1=AN+1AAN+11(9.3)\frac{v_{N+1} - v_1}{v_{N+1} - v_1^*} = \frac{A^{N+1} - A}{A^{N+1} - 1} \tag{9.3}

其中:

  • vN+1v_{N+1} — 进塔气中组分 ii 的流量,kmol/h
  • v1v_1 — 出塔气中组分 ii 的流量,kmol/h
  • v1v_1^* — 与入口吸收剂平衡时,组分 ii 在气相中的理论流量,kmol/h。由 v1=KiL0xi0v_1^* = K_i \cdot L_0 x_{i0} 计算
  • L0L_0 — 塔顶进入的吸收剂(贫液)总流量,kmol/h
  • xi0x_{i0} — 吸收剂中组分 ii 的摩尔分数(新鲜吸收剂通常 xi00x_{i0} \approx 0

当吸收剂不含溶质时 v1=0v_1^* = 0,相对吸收率等于吸收率 ϕ\phi

理论板数计算

已知关键组分的吸收因子 AkeyA_{key} 和吸收率 ϕkey\phi_{key},由 Kremser 方程反解:

N=log(Akeyϕkey1ϕkey)logAkey1(9.4)N = \frac{\log\left(\frac{A_{key} - \phi_{key}}{1 - \phi_{key}}\right)}{\log A_{key}} - 1 \tag{9.4}

最小液气比

在无穷多塔板条件下(NN \to \infty),吸收因子等于吸收率:Akey,min=ϕkeyA_{key,\min} = \phi_{key}

(LV)min=Kkeyϕkey(9.5)\left(\frac{L}{V}\right)_{\min} = K_{key} \cdot \phi_{key} \tag{9.5}

实际液气比取最小液气比的 1.2~2.0 倍:

LV=(1.22.0)(LV)min(9.6)\frac{L}{V} = (1.2 \sim 2.0) \left(\frac{L}{V}\right)_{\min} \tag{9.6}

非关键组分的吸收率

各组分液气比相同,非关键组分吸收因子由关键组分换算:

Ai=AkeyKkeyKi(9.7)A_i = A_{key} \frac{K_{key}}{K_i} \tag{9.7}

吸收率由 Kremser 方程计算:

ϕi=AiN+1AiAiN+11(9.8)\phi_i = \frac{A_i^{N+1} - A_i}{A_i^{N+1} - 1} \tag{9.8}

平均有效吸收因子法(Edmister)

前述 Kremser 方程假设全塔 AA 恒定,但实际塔中塔顶和塔底温度不同,KiK_i 值差异显著。Edmister 提出用塔顶和塔底条件的吸收因子计算一个"有效"平均值。

塔顶和塔底的吸收因子

A1=(LKV)top,AN=(LKV)bottom(9.9)A_1 = \left(\frac{L}{K V}\right)_{\text{top}},\quad A_N = \left(\frac{L}{K V}\right)_{\text{bottom}} \tag{9.9}

其中塔顶温度低、KK 值小 → A1A_1 较大;塔底温度高、KK 值大 → ANA_N 较小。

有效吸收因子推导

Edmister 通过分析全塔物料衡算,发现有效吸收因子可以由塔顶和塔底的 AA 值组合得到。推导分两步:

第一步:用塔底条件定义辅助量 AeA_e'——该值可理解为"以塔底吸收因子为基准、经塔顶修正后的加权平均":

Ae=AN(A1+1)AN+1(9.10)A_e' = \frac{A_N (A_1 + 1)}{A_N + 1} \tag{9.10}

AeA_e' 本身并非最终使用的有效因子,它是构建最终公式的中间变量。

第二步:将 AeA_e'A1A_1 作几何平均并加经验修正项 0.25,得到最终的 Edmister 有效吸收因子:

Ae=A1(AN+1)+0.250.5(9.11)\boxed{A_e = \sqrt{A_1 (A_N + 1) + 0.25} - 0.5} \tag{9.11}

A1ANA_1 \gg A_N(塔顶吸收远强于塔底,多数吸收塔的实际情况)时,式 (9.10) 可作为 AeA_e 的简化近似直接使用:

AeAe=ANA1+1AN+1(9.12)A_e \approx A_e' = A_N \cdot \frac{A_1 + 1}{A_N + 1} \tag{9.12}

使用方式:将 AeA_e 替代 Kremser 方程 (9.3) 中的恒定 AA,即 ϕ=AeN+1AeAeN+11\phi = \dfrac{A_e^{N+1} - A_e}{A_e^{N+1} - 1}。当塔顶塔底温差显著时,Edmister 修正大幅提高计算精度。

关键参数

符号含义单位
AA吸收因子 A=L/(KV)A = L/(KV)无量纲
ϕ\phi吸收率无量纲
NN理论板数无量纲
L/VL/V液气比无量纲
AeA_eEdmister 有效吸收因子无量纲

前置知识

  • 相平衡常数 Ki=Pisat/PK_i = P_i^{\text{sat}} / P,由 Antoine 方程 在操作温度下计算 → 4.1.1 相平衡常数
  • 吸收塔与精馏塔同属气液传质设备,但吸收无冷凝器和再沸器 → 第 7 章 多组分精馏

关联条目


← 返回目录

讨论与反馈

对本文内容有疑问或建议?请在此留言。评论使用 GitHub Discussions 存储。