Skip to content

例 9.1 — 多组分吸收塔设计

题目

操作压力 1.06 MPa,进料气体流量 100 kmol/h,温度 32°C,组成:

组分CH₄C₂H₆C₃H₈n-C₄H₁₀n-C₅H₁₂
摩尔分数0.700.120.080.060.04
分类NKNKNKKeyNK

吸收剂为不挥发油,温度 32°C。假设全塔平均温度 37°C,操作压力 1.06 MPa。 要求 n-C₄H₁₀ 的回收率为 99.5%

计算:

  1. 最小液气比 (L/V)min(L/V)_{\min}
  2. 实际液气比 L/V=2.0(L/V)minL/V = 2.0 (L/V)_{\min} 时的理论板数
  3. 各组分的吸收率和塔顶尾气的流量与组成
  4. 塔顶应加入的吸收剂量

已知条件整理

参数
总压 PP1.06 MPa
温度 TT37°C(310.15 K)
进料 FF100 kmol/h
关键组分n-C₄H₁₀
要求吸收率 ϕkey\phi_{key}0.995
L/VL/V 倍数2.0 × 最小

求解过程

步骤 1:计算各组分的相平衡常数

关键Ki=Pisat/PK_i = P_i^{\text{sat}} / P,由 Antoine 方程在 37°C 下求得饱和蒸气压。

组分PisatP_i^{\text{sat}} (37°C)Ki=Pisat/1.06MPaK_i = P_i^{\text{sat}} / 1.06\text{MPa}
CH₄71.53
C₂H₆4.59
C₃H₈1.19
n-C₄H₁₀ (Key)0.33
n-C₅H₁₂0.099

步骤 2:最小液气比

在无穷多塔板条件下(N=N = \infty),吸收因子等于吸收率:

Akey,min=ϕkey=0.995A_{key,\min} = \phi_{key} = 0.995

(LV)min=KkeyAkey,min=0.33×0.995=0.328\left(\frac{L}{V}\right)_{\min} = K_{key} \cdot A_{key,\min} = 0.33 \times 0.995 = \mathbf{0.328}

步骤 3:实际液气比和理论板数

实际液气比:

LV=2.0×0.328=0.656\frac{L}{V} = 2.0 \times 0.328 = \mathbf{0.656}

关键组分的实际吸收因子:

Akey=LKkeyV=0.6560.33=1.99A_{key} = \frac{L}{K_{key} V} = \frac{0.656}{0.33} = \mathbf{1.99}

由 Kremser 方程反算理论板数:

N=log(Akeyϕkey1ϕkey)logAkey1N = \frac{\log\left(\dfrac{A_{key} - \phi_{key}}{1 - \phi_{key}}\right)}{\log A_{key}} - 1

N=log(1.990.99510.995)log1.991=log(199)0.2991=2.2990.2991=7.691=6.69N = \frac{\log\left(\dfrac{1.99 - 0.995}{1 - 0.995}\right)}{\log 1.99} - 1 = \frac{\log(199)}{0.299} - 1 = \frac{2.299}{0.299} - 1 = 7.69 - 1 = 6.69

圆整:N=7\boxed{N = 7} 块理论板。

步骤 4:非关键组分的吸收率和尾气组成

非关键组分的吸收因子由关键组分的 AkeyA_{key} 换算(液气比相同):

Ai=AkeyKkeyKi=1.99×0.33KiA_i = A_{key} \frac{K_{key}}{K_i} = 1.99 \times \frac{0.33}{K_i}

CH₄

A1=1.99×0.3371.53=0.0092A_1 = 1.99 \times \frac{0.33}{71.53} = 0.0092

ϕ1=A1N+1A1A1N+11=0.009280.00920.0092810.0092\phi_1 = \frac{A_1^{N+1} - A_1}{A_1^{N+1} - 1} = \frac{0.0092^8 - 0.0092}{0.0092^8 - 1} \approx 0.0092

被吸收量:70×0.0092=0.6470 \times 0.0092 = 0.64 kmol/h

C₂H₆

A2=1.99×0.334.59=0.143A_2 = 1.99 \times \frac{0.33}{4.59} = 0.143

ϕ2=0.14380.1430.143810.143\phi_2 = \frac{0.143^8 - 0.143}{0.143^8 - 1} \approx 0.143

被吸收量:12×0.143=1.7212 \times 0.143 = 1.72 kmol/h

C₃H₈

A3=1.99×0.331.19=0.552A_3 = 1.99 \times \frac{0.33}{1.19} = 0.552

ϕ3=0.55280.5520.55281=0.548\phi_3 = \frac{0.552^8 - 0.552}{0.552^8 - 1} = 0.548

被吸收量:8×0.548=4.388 \times 0.548 = 4.38 kmol/h

n-C₄H₁₀ (Key)

A4=1.99A_4 = 1.99

被吸收量:6×0.995=5.976 \times 0.995 = 5.97 kmol/h

n-C₅H₁₂

A5=1.99×0.330.099=6.63A_5 = 1.99 \times \frac{0.33}{0.099} = 6.63

ϕ5=6.6386.636.63811.0\phi_5 = \frac{6.63^8 - 6.63}{6.63^8 - 1} \approx 1.0

被吸收量:4×1.0=4.004 \times 1.0 = 4.00 kmol/h


步骤 5:物料衡算 — 求吸收剂流量

关键认知:吸收塔中气液两相的流量沿塔高是变化的——气体上升过程不断损失溶质(流量递减),液体下降过程不断吸收溶质(流量递增)。因此需要用平均流量来关联液气比。

记塔顶和塔底的液体流量分别为 L0L_0(贫液)和 LL(富液),气体流量分别为 VV(塔顶尾气)和 FF(塔底进料气)。

总物料衡算(气体减少量 = 液体增加量 = 总吸收量):

FV=LL0=Δ=16.71 kmol/hF - V = L - L_0 = \Delta = 16.71\ \text{kmol/h}

平均流量

Vavg=F+V2=100+83.292=91.64 kmol/hV_{avg} = \frac{F + V}{2} = \frac{100 + 83.29}{2} = \mathbf{91.64}\ \text{kmol/h}

由液气比 L/V=0.656L/V = 0.656 得平均液体流量:

Lavg=0.656×Vavg=0.656×91.64=60.12 kmol/hL_{avg} = 0.656 \times V_{avg} = 0.656 \times 91.64 = \mathbf{60.12}\ \text{kmol/h}

由平均液体流量反推塔顶加入量

Lavg=L0+L2=L0+(L0+Δ)2=L0+Δ2L_{avg} = \frac{L_0 + L}{2} = \frac{L_0 + (L_0 + \Delta)}{2} = L_0 + \frac{\Delta}{2}

L0=LavgΔ2=60.1216.712=51.76 kmol/h\therefore\quad L_0 = L_{avg} - \frac{\Delta}{2} = 60.12 - \frac{16.71}{2} = \mathbf{51.76}\ \text{kmol/h}

塔底富液流量:L=L0+Δ=51.76+16.71=68.47 kmol/hL = L_0 + \Delta = 51.76 + 16.71 = 68.47\ \text{kmol/h}


结果汇总

参数参数
最小液气比 (L/V)min(L/V)_{\min}0.328理论板数 NN7
实际液气比 L/VL/V0.656塔顶尾气 VV83.29 kmol/h
关键组分吸收率99.5%吸收剂流量 L0L_051.76 kmol/h
组分KiK_iAiA_iϕi\phi_i尾气 (kmol/h)吸收量 (kmol/h)
CH₄71.530.0090.00969.360.64
C₂H₆4.590.1430.14310.281.72
C₃H₈1.190.5520.5483.624.38
n-C₄H₁₀0.331.990.9950.035.97
n-C₅H₁₂0.0996.63~1.0~04.00

方法总结:吸收塔简捷计算流程

  1. 由 Antoine 方程计算各组分 PisatP_i^{\text{sat}}Ki=Pisat/PK_i = P_i^{\text{sat}} / P
  2. 关键组分 Akey,min=ϕkeyA_{key,\min} = \phi_{key}(L/V)min=Kkeyϕkey(L/V)_{\min} = K_{key} \cdot \phi_{key}
  3. L/V=(1.22.0)(L/V)minL/V = (1.2 \sim 2.0)(L/V)_{\min}Akey=L/(KkeyV)A_{key} = L/(K_{key}V)
  4. Kremser 方程 → NN
  5. Ai=AkeyKkey/KiA_i = A_{key} \cdot K_{key}/K_iϕi\phi_i → 物料衡算

关联条目


← 返回目录

讨论与反馈

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