3. 塔板数与回流比的理论关系
在精馏塔设计中,塔板数与回流比是一对相互制约的核心参数:回流比越大,达到相同分离要求所需的塔板数越少,但再沸器负荷和操作费用上升;回流比越小,塔板数增加,设备投资上升。Gilliland关联图、Underwood方程与Fenske-Underwood-Gilliland简捷法正是定量描述这一制约关系的经典工具。本节从三个层次展开:先建立最小回流比与最小塔板数两个极限,再用Gilliland关联把两者与操作回流比、理论板数联系起来,最后给出可执行的简捷法计算流程。
一、两个极限:最小回流比与最小塔板数
最小回流比 Rmin对应的是“无穷多塔板”的操作极限。当回流比降低到某一数值时,精馏段与提馏段的操作线在平衡曲线上的某一点(夹点)相交,此时达到规定分离所需的塔板数为无穷大。Rmin是操作回流比的下限,实际回流比通常取 R = (1.1~2.0)Rmin。
对于相对挥发度 α 近似恒定的二元或窄沸程物系,Underwood方程给出最小回流比的解析解:
Σ [ αᵢ·x_F,i / (αᵢ - θ) ] = 1 - q
R_min + 1 = Σ [ αᵢ·x_D,i / (αᵢ - θ) ]
其中 q 为进料热状态参数(饱和液体进料 q=1,饱和蒸汽进料 q=0),θ 为根,位于轻、重关键组分相对挥发度之间。求解时先由第一式解出 θ,再代入第二式得到 Rmin。多组分精馏中,Underwood方程按关键组分处理,非关键组分按分配估算。
最小塔板数 Nmin对应全回流(R→∞)极限,此时操作线与对角线重合,分离能力最大。Fenske方程给出:
N_min = lg[ (x_D,LK / x_D,HK) · (x_B,HK / x_B,LK) ] / lg α_LK,HK
式中 LK、HK 分别代表轻、重关键组分,α 取塔顶与塔底相对挥发度的几何平均值。Nmin 是理论板数的下限,实际理论板数 N 必大于 Nmin。
二、Gilliland关联图:把两个极限连起来
Gilliland关联图以两个无量纲坐标为轴:横坐标 X = (R - Rmin)/(R + 1),纵坐标 Y = (N - Nmin)/(N + 1)。图中一条经验曲线把 X 与 Y 对应起来,其物理含义是:当操作回流比接近最小回流比时(X→0),所需理论板数趋于无穷(Y→1);当回流比远大于最小回流比时(X→1),理论板数趋近最小塔板数(Y→0)。
工程上常用Eduljee拟合式替代查图,便于编程计算:
Y = 1 - exp[ (1 + 54.4X)/(11 + 117.2X) · (X - 1)/√X ]
由 Y 反算理论板数:N = (Nmin + Y)/(1 - Y)。该式在 0.01 ≤ X ≤ 1 范围内精度较好,X 很小时误差增大,需谨慎使用。
三、Fenske-Underwood-Gilliland简捷法计算流程
简捷法(FUG法)是Aspen Plus中DSTWU模块的算法基础,其计算顺序固定为“Fenske→Underwood→Gilliland”,具体步骤如下:
- 确定关键组分与分离要求:指定轻关键组分在塔顶的回收率、重关键组分在塔底的回收率(或塔顶/塔底组成)。
- Fenske求 Nmin:用塔顶、塔底关键组分组成与几何平均相对挥发度计算最小理论板数,并可用Fenske方程估算非关键组分的分配。
- Underwood求 Rmin:由进料组成、q值与相对挥发度解 θ,再求 Rmin。
- 选定操作回流比:按 R = (1.1~2.0)Rmin 取初值,计算 X。
- Gilliland求 N:由X查图或Eduljee式求Y,反算理论板数 N,并确定进料板位置(常用Kirkbride经验式)。
四、示例:苯—甲苯常压精馏的简捷估算
设进料流量100 kmol/h,组成 xF=0.45(苯摩尔分数),饱和液体进料(q=1),塔顶 xD=0.98,塔底 xB=0.02,平均相对挥发度 α=2.45。以下数值为教学示例假设,非实测结果。
Fenske:Nmin = lg[(0.98/0.02)·(0.98/0.02)] / lg2.45 = lg(49×49)/0.389 ≈ 2×1.690/0.389 ≈ 8.7,取9块理论板。
Underwood:q=1,解 Σ[αxF/(α−θ)] = 0,得 θ≈1.62;代入 Rmin+1 = Σ[αxD/(α−θ)],得 Rmin≈1.05。
取 R=1.5Rmin≈1.58,则 X=(1.58−1.05)/(1.58+1)=0.205。由Eduljee式得 Y≈0.62,反算 N=(8.7+0.62)/(1−0.62)≈24.5,取25块理论板。若改取 R=2.0Rmin≈2.10,X≈0.34,Y≈0.50,N≈18.4,取19块。可见回流比从1.5Rmin提高到2.0Rmin,理论板数由25降至19,但再沸器负荷相应增加约33%。
| 方案 | R/R_min | X | Y | 理论板数 N | 相对再沸负荷 |
|---|---|---|---|---|---|
| A | 1.2 | 0.10 | 0.72 | 约33 | 1.00 |
| B | 1.5 | 0.21 | 0.62 | 约25 | 1.25 |
| C | 2.0 | 0.34 | 0.50 | 约19 | 1.67 |
表中“相对再沸负荷”以方案A为基准,按回流比近似线性估算,仅用于说明趋势。实际选型需在设备投资(塔板数、塔径)与操作费用(再沸器蒸汽、冷凝器冷却水)之间做总费用最小化权衡,通常最优点落在 R/Rmin=1.1~1.5 区间。
五、在Aspen Plus中的对应操作
在Aspen Plus中新建DSTWU模块,输入进料条件、压力、关键组分及回收率,并指定“Reflux ratio”为1.1~2.0倍最小回流比,软件即按FUG法输出最小理论板数、最小回流比、实际理论板数与进料板位置。若需更精确结果,可将DSTWU结果作为RadFrac的初值,再启用严格计算并考察板效率与温度分布。理解FUG三式的物理含义,才能在软件给出结果时判断其合理性——例如当Rmin计算值异常偏大时,往往提示相对挥发度取值或关键组分设定有误。