3. 物性方法选择与验证

在 Aspen Plus 中，物性方法决定了逸度系数、活度系数、焓与熵的计算路径，也直接决定闪蒸结果是否可信。选错方法，后续的收敛、灵敏度分析和设计都会被系统性偏差污染。本节围绕三个问题展开：如何判别体系理想性、如何在状态方程法与活度系数法之间取舍、以及如何查询与回归二元交互参数并用案例验证。

一、先判别体系：理想与非理想的三个判据

判别不是凭经验拍板，而是看三个量。第一是操作压力与临界压力之比：当 P < 0.2·Pc 且温度远离临界区时，气相近似理想气体，液相非理想性成为主要矛盾。第二是极性：若组分含羟基、羧基、胺基或强氢键，液相活度系数 γ 显著偏离 1，必须用活度系数模型。第三是共沸与液液分相：存在共沸点或互溶度间隙时，任何单一状态方程都难以准确描述，需活度系数法或带 GE 混合规则的 EOS。

一个实用的量化判据是无限稀释活度系数 γ∞。若所有组分的 γ∞ 落在 0.9～1.1 之间，体系接近理想，可用 Raoult 定律或理想模型；若 γ∞ 大于 3 或小于 0.3，则强非理想，必须选 NRTL、UNIQUAC 或 UNIFAC。例如乙醇—水在 1 atm 下 γ∞(乙醇) 约 3～4，属典型强非理想体系，用理想模型计算的气液平衡会有明显偏差。

二、两类方法的适用边界

状态方程法（PR、SRK）以逸度系数为基础，同时描述气液两相，天然适合高压、含轻气体（H₂、N₂、CO₂、CH₄）的体系，如加氢反应器、气液分离器、天然气脱水。PR 对液相密度和临界区优于 SRK，是炼化与气体处理的首选。其短板在于强极性、强缔合体系：常规 PR 难以准确给出乙醇—水的共沸组成，除非引入 Wong-Sandler 或 Huron-Vidal 混合规则并回归 kij。

活度系数法（NRTL、UNIQUAC）以液相活度系数描述非理想性，气相用理想气体或 EOS 修正，适合低压至中压（一般 < 10 bar）的极性、共沸、液液萃取体系，如精馏、吸收、萃取。NRTL 对强非理想和部分互溶体系拟合能力强，参数少、易回归；UNIQUAC 的组合项使其对分子大小差异大的体系更稳健，且可外推到多元。两者都需要二元交互参数，缺参数时可用 UNIFAC 基团贡献法估算作为初值。

三、二元交互参数的查询与回归

Aspen Plus 内置 APV 数据库（VLE-IG、VLE-RK、LLE-ASPEN 等），在 Properties → Parameters → Binary Interaction 中按组分对检索。以 PR 为例，需要的是 kij；以 NRTL 为例，需要 Aij、Aji、αij（αij 常取 0.2～0.3，强非理想可取 0.3）。若数据库无记录，先用 UNIFAC 估算，再以实验数据回归。

回归流程：在 Data Regression 中建立数据集，输入 T、P、x、y（或 γ∞），选择目标函数（如最大似然），勾选待回归参数，运行后查看残差与参数置信区间。示例：乙醇(1)—水(2) 在 1 atm 下，NRTL 回归得到 A12 ≈ 0.28、A21 ≈ 1.42（示例量级，实际以回归结果为准），α12 = 0.3，此时共沸点预测约 78.2 ℃、x₁ ≈ 0.89，与实验值接近。若残差呈系统趋势，说明模型形式不足，应改用 UNIQUAC 或引入温度相关项。

四、验证案例：乙醇—水闪蒸

设进料 100 kmol/h，乙醇摩尔分数 0.4，水 0.6，温度 85 ℃，压力 1 atm，绝热闪蒸。分别用 NRTL 与 PR 计算，观察气相组成与温度。NRTL 给出气相乙醇约 0.62、闪蒸温度约 82 ℃；PR 因未准确描述氢键，气相乙醇偏低约 0.05～0.08，温度偏差 1～2 K。再用 Txy 图叠加实验点，NRTL 曲线与实验点吻合，PR 在富乙醇区明显偏离。结论：该体系应选 NRTL，若必须用 PR，则需回归 kij 并检查共沸预测。

最后给出可执行的检查清单：确认体系极性与压力区间；查 APV 参数并记录来源；缺参数时用 UNIFAC 估算并标注“估算值”；回归后检查残差随机性；用 Txy 与实验对比；在报告中写明方法、参数版本与验证偏差。做到这五步，物性方法的选择才是有据可依的工程决策，而非默认选项。