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 很小时误差增大,需谨慎使用。

Gilliland关联:回流比与理论板数的相互制约 X = (R − R_min)/(R + 1) Y = (N − N_min)/(N + 1) 0 0.25 0.5 0.75 1.0 0 0.25 0.5 0.75 1.0 R→R_min:N→∞(夹点) R≫R_min:N→N_min(全回流) 工作点:X=0.5, Y≈0.37 读图逻辑 ① 由Underwood求 R_min ② 由Fenske求 N_min ③ 选定R,算X,查Y ④ 反算 N,校核经济性

三、Fenske-Underwood-Gilliland简捷法计算流程

简捷法(FUG法)是Aspen Plus中DSTWU模块的算法基础,其计算顺序固定为“Fenske→Underwood→Gilliland”,具体步骤如下:

  1. 确定关键组分与分离要求:指定轻关键组分在塔顶的回收率、重关键组分在塔底的回收率(或塔顶/塔底组成)。
  2. Fenske求 Nmin:用塔顶、塔底关键组分组成与几何平均相对挥发度计算最小理论板数,并可用Fenske方程估算非关键组分的分配。
  3. Underwood求 Rmin:由进料组成、q值与相对挥发度解 θ,再求 Rmin。
  4. 选定操作回流比:按 R = (1.1~2.0)Rmin 取初值,计算 X。
  5. Gilliland求 N:由X查图或Eduljee式求Y,反算理论板数 N,并确定进料板位置(常用Kirkbride经验式)。
简捷法得到的是理论板数与回流比的初步估计,未考虑板效率、非理想热力学与严格水力学。工程上通常以FUG结果作为严格模型(RadFrac)初值,再逐步校正。

四、示例:苯—甲苯常压精馏的简捷估算

设进料流量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_minXY理论板数 N相对再沸负荷
A1.20.100.72约331.00
B1.50.210.62约251.25
C2.00.340.50约191.67

表中“相对再沸负荷”以方案A为基准,按回流比近似线性估算,仅用于说明趋势。实际选型需在设备投资(塔板数、塔径)与操作费用(再沸器蒸汽、冷凝器冷却水)之间做总费用最小化权衡,通常最优点落在 R/Rmin=1.1~1.5 区间。

五、在Aspen Plus中的对应操作

在Aspen Plus中新建DSTWU模块,输入进料条件、压力、关键组分及回收率,并指定“Reflux ratio”为1.1~2.0倍最小回流比,软件即按FUG法输出最小理论板数、最小回流比、实际理论板数与进料板位置。若需更精确结果,可将DSTWU结果作为RadFrac的初值,再启用严格计算并考察板效率与温度分布。理解FUG三式的物理含义,才能在软件给出结果时判断其合理性——例如当Rmin计算值异常偏大时,往往提示相对挥发度取值或关键组分设定有误。

注意:Gilliland关联与Eduljee式均为经验拟合,对非理想性强的物系(如共沸、强极性)误差较大;此时应以严格模型或实验数据为准,简捷法仅作初值参考。