2. 热力学基础与物性方法选择

在 Aspen Plus 中做能量衡算,真正决定热流数值是否可信的,往往不是流程图画得对不对,而是两件事:焓的基准与计算路径,以及物性方法是否匹配物系。同一股物流,换一个状态方程,焓值可能相差几个百分点;对换热负荷动辄数百千瓦的单元,这个偏差足以让设备选型失准。本节围绕“焓、熵、热容的基本关系”和“如何为物系挑选状态方程或活度系数模型”两条主线展开,并给出可操作的判断依据。

一、焓、熵、热容:能量衡算的三块基石

焓是状态函数,工程上关心的是焓差而非绝对值。Aspen Plus 对每个组分都定义了参考态(理想气体在 298.15 K 或 0 K 等,取决于物性数据库与选项),因此软件输出的摩尔焓是“相对参考态的绝对值”,做能量衡算时必须保证所有物流、所有单元使用同一参考态,否则相减得到的 ΔH 没有意义。这一点在多股进料、循环回路和撕裂流收敛时尤其容易出问题。

热容是焓对温度的导数:Cp = (∂H/∂T)p。理想气体热容通常用多项式关联(如 Cp = a + bT + cT² + dT³),而真实气体还需加上剩余焓(departure enthalpy),即压力对焓的贡献。液相则通过活度系数模型或状态方程计算超额焓。因此,一个完整的焓计算链条是:

H(T,p) = Hig(T) + [H(T,p) − Hig(T)]剩余

其中理想气体部分由热容积分得到,剩余部分由状态方程或活度系数模型给出。熵同理,且熵直接进入吉布斯自由能,影响相平衡判断。理解这条链条,就能明白为什么“物性方法选错,热流一定错”——错误的方法首先扭曲了剩余焓和相平衡,再传导到换热负荷。

二、物性方法选择的判断逻辑

Aspen Plus 的物性方法大致分三类:状态方程法(如 PR、RK-SOAVE、SRK)、活度系数法(如 NRTL、UNIQUAC、UNIFAC、WILSON)和专用模型(如电解质 ELECNRTL、聚合物 POLYNRTL、石油 API 等)。选择的核心判断是:物系是否非理想、是否含极性组分、是否处于高压或近临界区、是否有强缔合或电解质。

工程上常用的经验规则是:

  • 烃类、轻气体、超临界或高压体系:优先状态方程,PR 或 RK-SOAVE 适用面广,对气液相平衡和剩余焓都有较好描述。
  • 极性组分(醇、酮、水、酸)在低压液相:活度系数模型更合适,NRTL 对强非理想混合物的液液平衡和汽液平衡精度高;UNIFAC 适合缺少二元交互参数时的估算。
  • 含水、含电解质或强缔合体系:ELECNRTL 或带 Hayden-O'Connell 处理的模型,否则气相缔合会被严重低估。
  • 宽沸程、跨临界流程:单一模型往往不够,需在 Aspen 中按温度/压力区间切换或使用 PR 配合适当的 alpha 函数。

需要强调的是,二元交互参数是活度系数模型的灵魂。NRTL 若缺少准确的 Aij、Aji、αij,计算结果可能比理想溶液还差。Aspen 内置的 VLE-IG、VLE-RK 等数据库可自动检索,但工程上仍应核对数据来源的温度范围和组分匹配性。

三、从物性方法到热流:一个可复现的示例

假设一股甲醇-水混合液,流量 100 kmol/h,甲醇摩尔分数 0.4,温度 60 °C,压力 1 atm,需预热至泡点后部分汽化。这是一个典型的极性-缔合体系,低压液相非理想性强。若误用 PR 状态方程,甲醇-水的交叉缔合被忽略,预测的泡点和汽化潜热会偏低;改用 NRTL(配合理想气体热容和 Watson 关联的汽化潜热),泡点约 78–80 °C(示例估算值,实际以软件计算为准),换热负荷明显高于 PR 结果。

操作上,在 Aspen 中应做三件事:第一,在 Methods 中选定 NRTL,并检查 Binary Parameters 是否检索到甲醇-水交互参数;第二,在 Report Options 中勾选摩尔焓、摩尔熵输出,确认参考态一致;第三,用 Sensitivity 分析对比两种方法下的换热负荷差异,差异超过 5% 即说明方法选择对结论有实质影响,必须在报告中说明依据。

判断物性方法是否“够用”的实用做法:对关键单元做一次方法敏感性扫描(如 PR vs NRTL),若换热负荷、相分率或产品纯度变化在工程允许误差内,可认为方法稳健;若差异显著,则需回到二元交互参数和数据来源核查。

四、常见误区与排障要点

第一,参考态混用:部分物流来自不同物性方法或不同数据库,焓基准不一致,导致能量衡算“凭空多出热量”。排障时检查所有物流的 enthalpy basis 是否统一。第二,忽略压力对液相焓的影响:低压下可忽略,但高压液相进料若用理想模型,剩余焓缺失会造成明显偏差。第三,活度系数模型用于气相:NRTL 只描述液相非理想性,气相仍需状态方程或理想气体假设,Aspen 默认组合需确认。第四,超临界区误用活度模型:此时液相概念消失,必须切换状态方程。

下面这张图把“物系特征 → 方法选择 → 焓计算链条 → 热流结果”的判断关系集中呈现,可作为选型时的快速对照。

物系特征 → 物性方法 → 焓计算链条 → 热流结果 物系特征判断 烃类 / 轻气体 高压、近临界 极性组分液相 醇 / 酮 / 水 电解质 / 强缔合 酸 / 胺 / 盐 宽沸程 / 跨临界 需分区切换 二元交互参数 是否齐全、温度范围 物性方法选择 PR / RK-SOAVE 状态方程法 NRTL / UNIQUAC 活度系数法 ELECNRTL 电解质专用 UNIFAC 估算 缺参数时 参数校验 VLE 数据回归 焓计算链条与热流 H_ig(T) = ∫Cp_ig dT 理想气体热容积分 H_剩余 = H(T,p) − H_ig(T) 由 EOS 或活度模型给出 H = H_ig + H_剩余 统一参考态后相减得 ΔH Q = ΔH × 流量 换热负荷 / 汽化潜热 方法敏感性校验 PR vs NRTL 对比换热负荷 差异 > 5% 需说明依据 偏差大则回查参数

把这张图当作选型检查表使用:先定位物系落在哪一类,再确认所选方法能否给出剩余焓和相平衡,最后用敏感性对比验证热流结论。对甲醇-水这类体系,NRTL 是起点而非终点,二元参数的来源和适用温度区间才是决定热流精度的最后一环。掌握这条从“物系特征”到“热流数值”的传导链,后续的换热器建模、塔器再沸器负荷核算才有可靠基础。