2. 共沸体系热力学基础

要判断一个共沸体系能否用普通精馏分开、该选变压精馏还是萃取精馏，第一步不是查工艺手册，而是回到热力学：这套体系在给定压力下到底有没有共沸点？共沸组成是多少？改变压力后共沸点会往哪个方向移动？这些问题的答案，全部藏在汽液平衡（VLE）关系和剩余曲线里。本节把这三件事串起来讲清楚：VLE 与剩余曲线提供"地图"，共沸点判定与组成计算给出"坐标"，活度系数模型的选择决定这张地图画得准不准。

一、汽液平衡关系：共沸判定的出发点

对二元体系，汽液平衡的基本关系是修正后的拉乌尔定律：

y_i · P = γ_i · x_i · P_i^sat(T) · Φ_i
i = 1, 2；  y_1 + y_2 = 1，x_1 + x_2 = 1

其中 x_i、y_i 为液相与汽相摩尔分数，P 为总压，P_i^sat 为纯组分饱和蒸气压（由 Antoine 方程给出），γ_i 为液相活度系数，Φ_i 为汽相校正因子（中低压下常近似为 1）。共沸点的定义是汽液两相组成相同，即 y_i = x_i。代入上式可得共沸判据：

γ_1 · P_1^sat = γ_2 · P_2^sat = P_az

这个式子说明：共沸是否出现，取决于"活度系数之比"与"饱和蒸气压之比"是否恰好抵消。定义相对挥发度 α_12 = (y_1/x_1)/(y_2/x_2) = γ_1 P_1^sat / (γ_2 P_2^sat)。当 α_12 在某个组成处等于 1，该组成就是共沸点。若 α_12 在全组成范围恒大于 1 或恒小于 1，则无共沸，普通精馏可分离；若 α_12 在两端分别大于 1 和小于 1（中间穿过 1），则存在共沸点。这是判断"要不要走特殊精馏"的第一道门槛。

二、剩余曲线：把 VLE 变成可读的"分离地图"

剩余曲线（residue curve）是间歇简单蒸馏中釜液组成随时间变化的轨迹，其微分方程为：

dx_i/dξ = x_i - y_i(x)，  ξ 为无量纲蒸馏时间

把不同初始组成对应的剩余曲线画在同一张三角图上，就得到剩余曲线图。它有三个关键读法：

• 节点与鞍点：剩余曲线的起点（不稳定节点）通常是低沸点纯组分或低沸共沸物，终点（稳定节点）是高沸点纯组分或高沸共沸物，鞍点则是介于两者之间的共沸点。
• 蒸馏边界：不同稳定节点吸引的曲线之间形成分界线，即蒸馏边界。边界两侧的物料不能用普通精馏跨越，这正是共沸体系"分不开"的几何本质。
• 可行性判断：进料点与目标产物若落在同一蒸馏区域内，普通精馏理论上可行；若分处两个区域，就必须借助变压、萃取、加盐等方式移动边界或改变共沸点。

剩余曲线图还有一个实用价值：它把"共沸点在哪、有几个、是最高还是最低共沸"一次性呈现出来。对二元最低共沸体系，剩余曲线从纯组分出发汇聚到共沸点；对最高共沸体系则相反，曲线从共沸点发散向两端纯组分。看到图形形态，就能快速判断体系类型。

三、共沸点判定与共沸组成计算

工程上判定共沸点通常走"先粗判、再精算"两步。粗判可用沸点比较：若两组分沸点相差很大（例如超过 30 K）且无强缔合，出现共沸的概率较低；若沸点接近或存在氢键、极性差异，则需认真核算。精算则解共沸判据方程。以乙醇-水为例（示例参数，非实测）：常压下乙醇沸点约 351.4 K，水约 373.15 K，二者形成最低共沸物，共沸温度约 351.3 K，共沸组成乙醇摩尔分数约 0.89。计算时给定 P = 101.325 kPa，用 Antoine 方程算 P_i^sat，用 NRTL 或 UNIQUAC 算 γ_i，令 y_i = x_i 联立求解，得到 T_az 与 x_az。

压力是移动共沸点最直接的杠杆。多数最低共沸体系的共沸组成随压力变化而漂移，例如乙醇-水在 101.3 kPa 下共沸乙醇约 0.89（摩尔分数），降至约 10 kPa 时共沸组成会明显向乙醇含量降低的方向移动。这个"漂移量"正是变压精馏（PSD）的设计依据：两个压力下的共沸组成差越大，PSD 越容易实现。因此计算共沸组成时，必须同时算至少两个压力点，并记录 Δx_az，作为方案比选的关键参数。

四、活度系数模型选择：预测准不准的分水岭

共沸预测的误差，八成来自 γ_i 模型选错。常用模型各有适用边界：

选择原则可以概括为：有可靠 VLE 实验数据时，优先用 NRTL 或 UNIQUAC 回归二元交互参数，并用另一压力或另一温度的数据做验证；无数据时用 UNIFAC 估算，但必须把结果标注为"估算值"，并在方案设计前补做至少一组 VLE 实验校核。特别提醒：乙醇-水这类强缔合体系，若误用 Margules 且参数取自稀溶液数据，预测的共沸组成可能偏差 0.05 以上，足以让后续 PSD 的压力选择失效。

五、把三者串成一条可执行的校核链

实际做方案选择时，建议按下面的顺序推进，而不是一上来就打开流程模拟软件：先用 Antoine 参数和选定的活度系数模型计算 α_12 随组成的变化，确认是否穿过 1；若有共沸，解共沸判据得到 T_az 与 x_az；再换一个压力重算，记录 Δx_az；然后用剩余曲线图确认进料与目标产物是否处于同一蒸馏区域；最后用一组独立的 VLE 实验数据（不同温度或不同压力）验证模型外推是否可靠。这条链走完，才能有把握地说"这个体系适合变压精馏"或"必须上萃取精馏"。