4. 模型参数获取与回归基础

NRTL与UNIQUAC的方程形式是固定的，真正决定模型能否用于工程计算的是二元交互参数。对NRTL，每个二元对需要三个参数：τ₁₂、τ₂₁与α₁₂；对UNIQUAC，需要两个能量参数Δu₁₂、Δu₂₁（常以τ₁₂=exp(−Δu₁₂/RT)形式出现），而r、q纯组分结构参数由分子体积与表面积确定，不参与回归。参数回归的本质，是用实验测得的相平衡信息反推这些参数，使模型在关心的组成与温度范围内复现真实相行为。因此，数据源的选择、目标函数的构造与约束的施加，共同决定了参数的质量。

三类数据在回归中的角色分工

无限稀释活度系数（γ^∞）是最“干净”的数据。它对应溶质在溶剂中浓度趋于零的极限，此时模型方程大幅简化：NRTL中γ^∞只与τ₁₂、τ₂₁、α₁₂有关，UNIQUAC中只与τ₁₂、τ₂₁及r、q有关，不涉及组成变量。γ^∞通常由气液色谱、惰性气体气提、沸点计等实验获得，精度较高，特别适合定端点——把参数在x→0和x→1两端的值锚定住。若只用中段VLE数据回归，端点附近常出现外推失真，而端点恰恰是萃取精馏溶剂筛选、痕量组分汽提等场景最关心的区域。

二元VLE数据（T、P、x、y）是回归的主力数据。它覆盖全组成范围，能同时约束参数在两端与中段的形状。对NRTL，α₁₂常被固定为0.2～0.3（非随机性因子），仅回归两个τ；也可将α作为第三参数一并回归，但需警惕参数相关性导致的病态。对UNIQUAC，两个能量参数由VLE数据回归即可，结构参数r、q固定。VLE数据回归时要注意热力学一致性：用Herington面积检验或点检验判断数据是否可靠，不一致的数据会直接把偏差“传染”给参数。

液液平衡（LLE）数据用于部分互溶体系。当二元系存在液液分层时，VLE数据无法覆盖两相区，必须用LLE的互溶度数据（x₁^α、x₁^β）约束参数。LLE回归的判据是两相中各组分的活度相等：x_i^αγ_i^α=x_i^βγ_i^β。对NRTL，LLE对α₁₂较敏感，常需将α适当调低（如0.2）以产生分层；对UNIQUAC，LLE数据是回归能量参数的关键约束，否则模型可能预测出完全互溶的错误结论。

目标函数怎么构造

回归的核心是定义一个能反映“模型与实验差多少”的目标函数。最常用的是活度系数残差形式：

F = Σᵢ wᵢ · [ ln γᵢ^cal(T, x; θ) − ln γᵢ^exp ]²

其中θ为待回归参数向量，wᵢ为权重。用lnγ而非γ本身，是因为活度系数量级跨度大，取对数后残差分布更均匀，数值优化更稳定。对VLE数据，γ^exp由实验的T、P、x、y经修正后算出（需考虑气相非理想性与Poynting校正）；对γ^∞数据，直接取极限值；对LLE数据，则用两相活度相等构造残差：

F_LLE = Σ [ (x₁ᵅγ₁ᵅ − x₁ᵝγ₁ᵝ)² + (x₂ᵅγ₂ᵅ − x₂ᵝγ₂ᵝ)² ]

实际回归常把多类数据合并成一个总目标函数，用权重平衡不同数据源的数量与精度差异。例如γ^∞数据精度高但点数少，可给较大权重；VLE数据点多但含系统误差，权重适当降低。权重不是随意定的，应使各类数据的残差贡献大致相当，否则点数多的一类会“淹没”另一类。

约束与可辨识性

NRTL的α₁₂与τ参数存在较强相关性，若三者同时自由回归，容易得到多组数学上等价但物理意义不同的解。工程上通常固定α₁₂=0.2或0.3，只回归两个τ，既减少参数相关性，又保留模型对非理想性的描述能力。UNIQUAC的两个能量参数同样需要设定合理边界，例如限制τ在0.01～100之间，防止优化过程跑向无物理意义的极端值。

另一个关键约束来自数据本身的可辨识性。若只有中段VLE数据，端点参数可能无法唯一确定，此时应补充γ^∞数据；若体系存在液液分层而只用VLE数据回归，模型可能给出错误的互溶判断，必须引入LLE数据。回归完成后，不能只看残差平方和的大小，还要检查：端点γ^∞是否复现、LLE分层是否出现、多温度下参数是否稳定、留出数据点的预测偏差是否可接受。只有这些校验都通过，参数才具备工程可用性。