4. 苯-甲苯精馏塔RadFrac严格模型搭建

苯-甲苯精馏塔的严格模型对象与自由度

本例用 RadFrac 搭建苯-甲苯二元精馏塔的严格平衡级模型。与简捷法不同,RadFrac 对每一块理论板同时求解 MESH 方程(物料平衡 M、相平衡 E、摩尔分率加和 S、焓平衡 H),因此必须把塔的“几何与操作”信息补全:理论板数、进料板位置、进料状态、塔顶压力与全塔压降、回流比或馏出率、以及两个产品规格之一。苯-甲苯为理想近似体系,活度系数接近 1,用 Raoult 定律加理想气相即可,物性方法选 IDEAL 或 RK-SOAVE 均可,本例用 IDEAL 便于对照手算。

RadFrac 的自由度关系是搭建时的核心约束:给定进料(流量、组成、温度/压力、气相分率)后,还需指定 理论板数 N、进料板位置 NF、塔顶压力 Ptop、全塔压降 ΔP、以及两个操作/规格变量(如回流比 R 与馏出率 D/F,或 R 与塔顶苯纯度)。少给一个规格,模型欠定;多给一个,模型过定并报自由度错误。

RadFrac 塔结构(N=12,进料第 7 板) 全凝器(第1板) 再沸器(第12板) 2 3 4 5 6 7 8 9 进料 F 第7板 馏出 D 回流 L 釜液 B 上升蒸气 V RadFrac 输入规格(本例取值) 理论板数 N = 12(含全凝器与再沸器) 进料板 N_F = 7(自上往下计) 塔顶压力 P_top = 101.3 kPa,全塔压降 ΔP = 0 进料:F = 100 kmol/h,z_F = 0.50 苯,饱和液体 操作规格:回流比 R = 2.0,馏出率 D/F = 0.45 物性方法:IDEAL(苯-甲苯近似理想) 自由度校核 已指定:N、N_F、P_top、ΔP、F、z_F、进料热状态、R、D/F 待求:各板 T、P、x_i、y_i、L、V 及 D、B 组成 规格数 = 2(R 与 D/F)→ 满足 RadFrac 定解要求 若改为指定 R 与塔顶苯纯度,则 D/F 变为结果 若同时给 R、D/F、塔顶纯度 → 过定,报自由度错误 若只给 R 不给第二规格 → 欠定,无法收敛 结果查看:Profiles → T / P / x / y 沿板分布 塔顶 T≈80.1℃(纯苯沸点),塔底 T≈110.6℃(甲苯沸点)

搭建步骤与关键设置

  1. 新建流程,从 Columns 选 RadFrac 放置,连接进料流 FEED、塔顶 DIST、塔底 BOT,并连接冷凝器回流与再沸器返回线(RadFrac 内部自动生成,无需外接换热器)。
  2. 在 Setup → Specifications 中设 Number of stages = 12,Condenser = Total,Reboiler = Kettle,Valid phases = Vapor-Liquid。
  3. 在 Streams 页设进料板 Feed stage = 7,Above-Stage Pressure = 101.3 kPa,Stage Pressure Drop = 0。
  4. 在 Setup → Specifications 的 Design Specs 区设两个操作变量:Reflux ratio = 2.0、Distillate rate = 45 kmol/h(对应 D/F=0.45)。
  5. 进料流设 F = 100 kmol/h、P = 101.3 kPa、Vapor fraction = 0、苯摩尔分率 0.50,物性方法选 IDEAL。
  6. 点击 Run,收敛后在 Blocks → RadFrac → Profiles 查看 T、P、x、y 沿板分布。

示例结果与判读

以下为按上述规格运行后应出现的预期趋势(示例值,用于对照判读,非实测数据):塔顶温度接近纯苯常压沸点 80.1 ℃,塔底接近甲苯沸点 110.6 ℃;温度沿板自上而下单调升高,进料板附近出现斜率变化(进料冷液或饱和液使该板温度略降后再升)。塔顶馏出液中苯摩尔分率约 0.95 以上,塔底釜液中苯摩尔分率约 0.03 以下,具体数值取决于 N、R、N_F 的组合。

观察项位置预期现象判读要点
温度 T第 1 板 → 第 12 板80.1 ℃ → 110.6 ℃ 单调升高若出现温度倒挂,检查压降设置或进料板位置
苯液相分率 x塔顶 / 塔底≈0.95 / ≈0.03与 D/F、R 是否匹配;偏离过大说明规格不合理
苯气相分率 y各板y > x(苯为轻组分)验证相平衡方向正确
回流比 R冷凝器等于设定值 2.0若被 Design Spec 调整,说明规格冲突
收敛状态Control PanelConverged,迭代次数 < 20不收敛先查自由度与规格
若 RadFrac 报“Degrees of freedom”错误,优先检查是否同时指定了 R、D/F 和产品纯度三个规格;若报“Not converged”,可先把 R 调大(如 3.0)或把 N 增至 15 再试,收敛后再回调到目标规格。
本例用 IDEAL 物性方法仅适用于苯-甲苯这类近似理想体系;换成醇-水等强非理想体系时,必须改用 NRTL/UNIQUAC 等活度系数模型,否则塔顶组成与温度分布会明显失真。