📚 课程目录第4章相对挥发度的计算路径

4. 相对挥发度的计算路径

相对挥发度 α 是衡量两组分分离难易程度的核心参数，定义为气液平衡条件下两组分挥发度之比。工程上最常用的形式是 α_ij = (y_i/x_i)/(y_j/x_j)。它直接决定了达到指定分离要求所需的理论板数和回流比，因此“α 从哪里来、怎么算准”是精馏设计的第一道关口。实际工作中，我们面对的数据条件差别很大：有的体系有完整的活度系数模型参数，有的只有纯组分饱和蒸气压，有的则只有几张实验测定的 T-x-y 数据表。不同条件下应选择不同的计算路径，本节按三条主线展开。

路径一：理想体系——Raoult 定律直接计算

当体系满足“分子结构相似、非极性或弱极性、液相为理想溶液”的条件时（如苯—甲苯、正己烷—正庚烷、邻二甲苯—对二甲苯等），液相活度系数 γ_i ≈ 1，气相可近似为理想气体。此时组分 i 的挥发度可写为 K_i = P_i^sat/P，相对挥发度简化为两组分饱和蒸气压之比：

α_ij = P_i^sat / P_j^sat

饱和蒸气压用 Antoine 方程计算：
ln(P^sat) = A - B / (T + C)     P^sat 单位 mmHg，T 单位 ℃

示例（苯—甲苯，常压 P = 760 mmHg）：
苯   A=6.90565, B=1211.033, C=220.790
甲苯 A=6.95464, B=1344.800, C=219.480

取 T = 95 ℃：
P_苯^sat  = exp(6.90565 - 1211.033/(95+220.790)) ≈ 1176 mmHg
P_甲苯^sat = exp(6.95464 - 1344.800/(95+219.480)) ≈ 476 mmHg
α = 1176 / 476 ≈ 2.47

这条路径的关键认识是：α 随温度变化，但变化缓慢。上例中温度从 80 ℃ 升到 110 ℃，α 大致在 2.5～2.4 之间浮动，工程上常取塔顶、塔底温度的几何平均值作为全塔代表值。Raoult 路径的优点是只需纯组分 Antoine 常数，计算量极小；缺点是适用范围窄，一旦体系极性增强或存在缔合、氢键，误差会迅速放大，此时必须转入第二条路径。

路径二：非理想体系——活度系数模型修正

当液相偏离理想显著时，挥发度需写成 v_i = γ_iP_i^sat，相对挥发度变为：

α_ij = (γ_i · P_i^sat) / (γ_j · P_j^sat)

γ 由活度系数模型给出，常见选择有：Margules 方程（二元、参数少、适合简单非理想体系）、Van Laar 方程（适合极性差异较大的二元系）、Wilson 方程（能描述强非理想但无法处理液液分层）、NRTL 与 UNIQUAC（可处理部分互溶体系，工程上最常用）。以二元 Wilson 模型为例，给定无限稀释活度系数 γ^∞ 即可回归出模型参数 Λ₁₂、Λ₂₁，再逐点计算 γ₁、γ₂ 并代入上式。

这里有一个必须警惕的现象：共沸点。当 γ₁P₁^sat = γ₂P₂^sat 时，α = 1，常规精馏无法继续提纯。例如乙醇—水在常压下约 78.2 ℃ 形成共沸，此时 α 趋近 1，必须改用变压精馏、萃取精馏或加盐精馏等手段。因此用活度系数模型算出 α 后，应检查 α 是否在组成范围内穿越 1，这是判断分离可行性的关键一步。

路径三：实验数据回归——没有模型参数时的兜底方案

当体系缺乏可靠的模型参数，或模型外推风险较大时，最稳妥的做法是直接使用实验测定的 T-x-y 或 x-y 数据。对每一组平衡数据点，由 y_i、x_i 直接算出该组成下的 α，再拟合为组成的函数。常用经验形式为：

ln α = A + B·x_1 + C·x_1²   （或简化为 ln α = A + B·x_1）

示例（假设的某二元体系实验数据，仅作演示）：
x1       y1       α = (y1/x1)/((1-y1)/(1-x1))
0.10    0.28      3.50
0.30    0.55      2.85
0.50    0.72      2.57
0.70    0.85      2.43
0.90    0.95      2.11

对 ln α 与 x1 作线性回归，得 A ≈ 1.36，B ≈ -0.55，
即 α ≈ exp(1.36 - 0.55·x1)，可用于全塔逐板计算。

回归时要注意三点：一是数据点应覆盖实际操作的组成范围，避免远距离外推；二是热力学一致性检验（如 Herington 面积检验）不可省略，否则回归出的 α 可能系统性偏差；三是若数据呈现 α 随组成剧烈变化，说明体系强非理想，应优先回到路径二用 NRTL 等模型拟合，而不是强行用线性式。

三条路径的选择逻辑

把上述内容串起来，实际工作中的判断顺序可以概括为下图：先看体系是否近似理想，是则走 Raoult 快速估算；否则检查是否有可信的活度系数模型参数，有则用模型修正并检查共沸；两者都不满足时，才动用实验数据回归。三条路径并非互斥，工程上常用 Raoult 做初算、模型做精算、实验数据做校核，形成交叉验证。

最后强调一个易错点：α 不是常数。Raoult 路径下它随温度缓变，活度系数路径下它随组成显著变化，实验回归路径下它本身就是组成的函数。因此无论走哪条路径，都应明确“这个 α 对应什么温度和组成”，并在逐板计算中按实际条件更新，而不是全塔套用一个固定值。把这一点做扎实，后续的分离难度判断和塔板数估算才有可靠基础。