3. UNIQUAC模型方程拆解

UNIQUAC（Universal Quasi-Chemical）模型把液体混合物的过量吉布斯自由能拆成两个物理来源截然不同的部分：组合项（combinatorial）描述分子大小和形状差异引起的熵效应，剩余项（residual）描述分子间相互作用能差异引起的焓效应。理解这一拆分，是判断UNIQUAC在什么工况下比NRTL更合适的关键。

一、组合项：只由纯组分参数r、q决定

组合项来自Flory-Huggins晶格理论，其表达式为：

g^C / RT = Σ xi·ln(Φi/xi) + (z/2)·Σ qi·xi·ln(θi/Φi)

Φi = ri·xi / Σ rj·xj        （体积分数）
θi = qi·xi / Σ qj·xj        （表面积分数）
z = 10                       （配位数，通常固定）

这里体积参数r反映分子在晶格中占据的格点数，近似正比于分子范德华体积；表面积分数参数q反映分子暴露在外的接触面积。r与q的比值r/q刻画分子形状：r/q接近1为球形小分子（如水、甲醇），r/q显著大于1为链状或大分子（如正己烷、长链醇）。组合项只依赖纯组分r、q和组成，不含任何二元交互参数，因此对同分异构体、链长差异大的体系，组合项本身就能贡献可观的非理想性。

二、剩余项：由二元交互参数τij承载

剩余项采用准化学近似：

g^R / RT = -Σ qi·xi·ln( Σ θj·τji )

τij = exp( -(uij - ujj) / RT ) = exp( -aij / T )

τij是二元交互参数，物理上表示i-j接触对与j-j接触对的能量差。注意下标方向：τij描述i周围出现j的相互作用，τij≠τji，两个方向需分别拟合。工程实现中常把τij写成温度函数τij=exp(-aij/T)，aij为可拟合常数；若温度范围窄，也可用固定τij。剩余项通过θi加权，因此q不仅影响组合项，也通过表面积分数进入剩余项——这是UNIQUAC中q的双重作用，也是它与NRTL结构上的重要差异。

三、τij的拟合逻辑

τij不能由纯组分数据得到，必须用二元相平衡数据回归。典型流程是：先用UNIFAC基团贡献法估算r、q（或查DECHEMA表），再固定r、q，用二元VLE的T、P、x、y数据以最大似然或最小二乘拟合aij、aji。由于组合项已固定，拟合实际只调整剩余项的两个参数，自由度低、外推较稳。对液液平衡（LLE）体系，τij对互溶度极其敏感，需用LLE数据单独拟合，不能直接套用VLE参数。

四、与NRTL的数学结构对应

两者都采用“局部组成”思想，剩余项形式高度相似：NRTL的τij=exp(-αij·τij^NRTL)中的指数项与UNIQUAC的τij=exp(-aij/T)在数学上同构，都是Boltzmann型温度依赖。差别在于：NRTL用摩尔分数加权的局部组成，参数为(gij-gjj)与αij三个；UNIQUAC用表面积分数加权，参数为aij、aji两个，且组合项由r、q显式给出。因此UNIQUAC对分子大小差异大的体系（如聚合物-溶剂、长链-短链）通常优于NRTL；而NRTL对强缔合、部分互溶体系通过αij可灵活调节，拟合能力更强。两者在二元VLE上精度接近，多元外推时UNIQUAC因参数少、组合项物理约束强，往往更稳健。

五、工程判断要点

选型时可按以下逻辑判断：若体系以尺寸/形状差异为主导（聚合物溶液、脂肪酸-烷烃、长链醇-水），优先UNIQUAC，因为组合项显式捕捉了这部分熵贡献，剩余项只需拟合较小的相互作用修正；若体系以强缔合或部分互溶为主导（醇-水、酮-水、胺-烃），NRTL的αij提供了额外的非随机性调节自由度，拟合更灵活。两者在常规二元VLE上精度相当，差异主要体现在多元外推和参数可获取性上：UNIQUAC参数少、物理约束强，外推通常更稳；NRTL参数多、拟合残差更小，但外推风险也更高。实际工程中应优先使用经过实验数据验证的参数集，而非盲目追求模型形式。