企业官网建设流程全解析

gPROMS二次开发教程（07）：物性与热力学——把真物性接进你的模型

版本声明块

• 工具/软件：gPROMS 桌面建模环境gPROMS ModelBuilder（gPROMS ModelBuilder为桌面建模工作台）；检索期官方发布锚点gPROMS Process 2022.1.0，适用版本以官方发布说明为准
• 语言/环境：gPROMS 语言（equation-oriented modelling）；实体MODELS，段PARAMETER/VARIABLE
• 本文目标：把"物性"从"模型里的几个常数"还原成"方程组的一半"，学会判断物性接入是否可靠、物性参数是否可引用

一句话结论：在 equation-oriented 模型里，物性不是一个装饰性的查询表，而是被联立进方程组的一批约束——一个闪蒸模型（F*z(i) = L*x(i) + V*y(i)配y(i) = K(i)*x(i)）里，K(i)随P、T和组成变化，物性一旦失真，物料衡算会照样"守恒"，结论却整体偏移；也因此，"热力学一致性"比"单点物性精确"更优先。

〇、本篇要解决的认知问题

• Q1：物性在模型里的"位置"到底在哪？它为什么不是常数，而是方程组的一部分？
• Q2：模型与物性之间的抽象层（物性包 / 属性调用）是如何工作的？为什么这层抽象能让你"换物性而不改模型"？
• Q3：相平衡与热力学一致性，为什么会直接影响反应与分离模型的结论？不一致有哪些可观察的症状？
• Q4：物性参数与文献数据要遵循什么来源与引用纪律？为什么"查到一个数就填进去"是高风险动作？
• Q5：物性错误会导致哪些典型求解失败？到第 11 篇做参数估计时，物性参数该怎么处理？

一、机制解析

1.1 物性在模型里的"位置"：它是方程的一部分，不是常数

初学者的心智模型通常是：“物性 = 一串常数，填进PARAMETER段即可”。这在教学例子里勉强成立，在真实模型里会立刻失真：

【错误心智模型】 【正确心智模型】 模型 + 物性常数 = 可解模型 模型方程 ┐ 物性错了 ⇒ 数值偏一点 物性约束 ├→ 联立求解 ⇒ 唯一解 求解器 ┘ 物性错了 ⇒ 非线性行为变了 ⇒ 解的位置与稳定性都可能变

关键推论：物性失真不会让模型"报错"，只会让它的非线性结构悄悄改变——唯一稳态可能变成多重稳态，收敛的迭代可能开始振荡。这类失败在现象上与"数值设置不对"极像，因此排查时必须先排除物性。

经验法则：拿到"莫名不收敛"的模型，先冻结物性再看：把强非线性物性临时换成工作点常数，若模型立刻稳定收敛，问题就在物性侧。

1.2 抽象层：模型只认"属性调用"，不认物性细节

gPROMS 平台在模型与物性之间隔了一层抽象：你的方程调用"属性"，物性系统负责按当前的 T、P、组成算出这个属性的值。这层抽象带来三个直接好处：

┌──────────────────────────────────────────────┐ │ 你的模型（方程） │ │ 守恒方程 / 本构方程 / 相平衡关系 │ └───────────────┬──────────────────────────────┘ │ 属性调用（向量：T, P, 组成 → 属性值） ┌───────────────▼──────────────────────────────┐ │ 物性 / 热力学 抽象层 │ │ 物性方法选择 · 参数集 · 混合规则 │ └───────────────┬──────────────────────────────┘ ┌───────────────▼──────────────────────────────┐ │ 物性系统的具体实现（平台提供 / 外部接入） │ └──────────────────────────────────────────────┘ 收益 1：换物性方法 ≠ 改模型。模型里只是"调用属性"，换了底层物性方法， 方程一字不改；这让"同一模型跑多种热力学假设"成为低成本实验。 收益 2：参数集集中管理。物性参数（如二元交互参数）集中在一处维护， 不必散落在各个模型里，避免"同一个物性在不同模型里取值不同"。 收益 3：可追溯。属性调用点就是模型对物性的全部依赖，做证据审查时 能一条条列出来核对（铁律 2：证据分级）。

铁律 4（商业闭源软件的诚实边界）：属性的具体名称、物性方法的选择方式、参数集的声明与覆盖语法、以及外部物性与高级热力学包的接入方式，全部属于 gPROMS 语言与平台细节，一律以随附官方《gPROMS ModelBuilder Guide》与平台文档为准。本文只讲"抽象层的职责划分"与"方程如何消费物性结果"，不点任何未在本系列事实底账中登记的物性产品或方法名。

1.3 相平衡与热力学一致性：为什么它决定结论

相平衡是分离与反应模型的核心非线性来源。以气液平衡为例，模型里有三类约束同时存在：

这三类约束的耦合方式决定了一个重要事实：

物料衡算永远会"守恒"，因此物性错误不会被守恒方程拦住——它只会把解推向一个"守恒但物理上错"的位置。这就是"物性错误最难查"的根本原因。

而"热力学一致性"指的是：不同物性属性之间必须相互自洽（例如相平衡关系与焓、熵的推导链条不能互相矛盾）。一致性被破坏后，典型症状有三种：

1. 能算出负浓度或摩尔分数越界：归一化被强行满足，但某个x(i)落到 0 以下；
2. 能量衡算与相平衡互相打架：按相平衡算出的组成，代入焓模型后能量衡算残差始终不降；
3. 同一个体系在不同初值下收敛到不同解：非线性结构被人为扭曲，出现伪多重稳态。

最佳实践：“物性一致"优先于"单点物性精确”。一个"精确的"纯组分常数若与体系其余物性推导链不匹配，往往比一套自洽但稍粗的物性更危险。工程做法是：整套方法与参数集一起选定、一起记录、一起做回归。

1.4 接入路线：平台内与平台外

物性来源可以分两条路线，两条路线的风险点完全不同：

铁律 7（外部连接必声明前提）：凡涉及外部连接（第 14 篇将系统展开gO:Python/gO:MATLAB/AdaptorFO/ gPROMS Foreign Process Interface 等机制），必须列出前置条件并声明参数细节以官方文档为准。物性接入是外部连接的一个高频场景，连接方式的具体支持范围与配置细节以官方文档为准，本文不做点名式推荐。

1.5 物性参数的来源与引用纪律

物性参数是模型里最容易被"凭印象填"的部分，因此本系列把它列为铁律级问题：

铁律 10（数值与参数不臆造）：动力学常数、物性参数、经验系数一律取自可引用文献或官方示例并标注来源，不得凭记忆填数。

执行层面，要求每个物性参数都能回答三个问题：（1）值是多少？（2）从哪来？（3）在什么范围内有效？第三问最常被忽略——很多参数只在有限温区内回归得到，把模型跑到温区之外，“不收敛"其实是"外推失效”。

1.6 物性引发的典型求解失败

1.7 与第 11 篇参数估计的分工

参数估计时最大的陷阱是把物性参数和动力学参数混在一起估：两者在数据上的可辨识性往往很差（"速率常数调大一点 + 相平衡常数调小一点"可能给出几乎相同的曲线），估出来的参数组即使残差很小，也可能毫无物理意义。本系列的处置策略是：

第一步：物性参数 → 用独立的物性数据（纯组分性质、二元相平衡数据）单独确定， 不放进动力学回归的决策变量集合。 第二步：动力学参数 → 在物性参数固定的前提下，用反应/过程数据回归（第 11 篇）。 第三步：一致性检查 → 把两步的结果放回同一模型，核对是否仍能复现全部数据。

经验法则：一次只估一类参数。参数估计的自由度是稀缺资源（铁律 8 的另一种表现）：把它花在物性上，动力学就失去了可辨识性。物性侧通常有独立数据可用，不该占用过程数据的解释力。

二、完整模型与逐行剖析

代码 2-1：气液平衡闪蒸模型（结构示意）

% --------------------------------------------------------------------------- % 说明：本块为"结构示意"。MODELS / PARAMETER / VARIABLE 为已确证关键字； % 属性调用（物性系统 → 模型）与相平衡属性的确切写法、物性方法的选择语法 % 以官方《gPROMS ModelBuilder Guide》为准，此处仅示意"方程如何消费物性"。 % 数值为教学示例值，正式使用须替换为可引用文献或官方示例数据（铁律 10）。 % 物理背景：单级等温闪蒸。进料 F（组成 z）在给定 T、P 下分成液相 L 与气相 V。 % --------------------------------------------------------------------------- MODELS PARAMETER N = 2 % 组分数：示意用二元体系 F = 1.0 % 进料摩尔流量（示意单位 mol/s） z(i) % 进料摩尔分数（向量参数：Σz = 1） T = 350.0 % 闪蒸温度（示例值） P = 1.0e5 % 闪蒸压力（示例 Pa） VARIABLE x(i) % 液相摩尔分数 y(i) % 气相摩尔分数 K(i) % 相平衡常数：由物性系统按 T、P、组成给出（物性侧量） L % 液相摩尔流量 V % 气相摩尔流量 %% ===== 物性侧：相平衡常数 K(i) ===== %% K(i) 由物性系统按当前 T、P 与组成计算（其内部走逸度系数或活度系数路径）。 %% 属性调用与物性方法选择的语法以官方《gPROMS ModelBuilder Guide》为准。 %% 本处只声明"K(i) 是物性侧的量"，并把它当作模型的输入属性消费。 %% —— 这一行是"科研复现"的关键：模型对物性的依赖必须显式写出来。 %% ===== 物料衡算（每组分一条）===== F*z(i) = L*x(i) + V*y(i) % E1 组分 i 的摩尔衡算（进料 = 液相 + 气相） %% ===== 相平衡关系 ===== y(i) = K(i)*x(i) % E2 相平衡：气相 = 相平衡常数 × 液相 %% ===== 归一化 ===== SUM(x(i)) = 1.0 % E3 液相摩尔分数归一（SUM 为求和算子示意写法） SUM(y(i)) = 1.0 % E4 气相摩尔分数归一 %% ===== 总物料衡算 ===== F = L + V % E5 总摩尔衡算

逐行剖析：

• z(i)写在PARAMETER段：进料组成是已知量（建模前确定）。注意它是向量参数，其维度必须与x(i)、y(i)一致（第 06 篇的维度一致性纪律）。
• K(i)写在VARIABLE段，并在注释中明确它是物性侧的量。这不是形式主义：读模型的人必须一眼看出"这个模型依赖哪些物性属性"；把依赖藏起来，模型就变成不可审查的黑箱。
• E1是逐组分物料衡算，E5是总物料衡算。两者不是独立的：把 E1 对全部组分求和即得 E5（在归一化与总流量定义下）。写全两条可以让物理意义更清楚，但自由度核算时必须意识到这一点。
• E3/E4两条归一化：这两条之中只有一条是独立约束（下面自由度核对里会证明）。这是闪蒸模型的经典结构性陷阱。
• 全式采用残差式（第 05 篇的规范）：F*z(i) - L*x(i) - V*y(i) = 0的写法把"进料减两股出料"原样保留，便于逐项核对是否漏项。

自由度核对（铁律 8）——这是本篇最值得带走的一段：

归一化冗余证明（对 N = 2，可逐项推广）： 把 E1 对 i 求和： F*Σz(i) = L*Σx(i) + V*Σy(i) 代入 Σz(i) = 1（进料归一）、Σx(i) = 1（E3）： F = L + V*Σy(i) 再代入 E5（F = L + V）： Σy(i) = 1 ← 即 E4 结论：E4 可由 E1(全部) + E3 + E5 导出，因此 E3 与 E4 中只有一个带来新信息。

这意味着什么：若把 E3、E4 都当作独立约束硬写进去，模型看起来"方程比未知数多一条"，实际自由度却是正确的——因为求解器面对的是线性相关的方程组，不会报错，但会给出奇怪的诊断信息（如"方程冗余"或残差不降）。初学者常在此处把一条归一化删掉，结果模型真正变成欠定。正确处置是：保留两条（提高可读性与鲁棒性），但在模型头里写明"两条归一化中仅一条独立"。

铁律 8（建模先算自由度）：自由度核算必须识别线性相关/可导出的方程。方程数不等于独立约束数；闪蒸归一化、总衡算与分衡算之间的依赖关系，是最常见的"数不对但解得了"的场合。

代码 2-2：物性消费清单——把"模型对物性的依赖"显式登记（结构示意）

% =========================================================================== % 物性消费清单（模型对物性的全部依赖，逐条登记） % 用途：换物性方法时按此清单核对；做证据审查时按此清单追溯；写论文时按此 % 清单在"物性假设"一节逐条声明。 % % 序 | 模型中的量 | 消费方式 | 依赖变量 | 参数来源 | 适用范围 % ---+------------------+-----------------+-----------------+----------------+----------- % 1 | K(i) 相平衡常数 | 属性调用 | T, P, x(i), y(i)| 参数集 A | 273~450 K % 2 | 液相密度 rho_L | 属性调用 | T, P, x(i) | 参数集 A | 常压附近 % 3 | 气相黏度 mu_V | 属性调用 | T, P, y(i) | 参数集 B | 300~600 K % 4 | 反应焓 dH_rxn | 属性调用（生成焓差）| T | 参数集 C | 标准态附近 % % --------------------------------------------------------------------------- % 换物性方法时的核对动作（三步） % 第 1 步：清单第 4 列的"依赖变量"是否与新方法一致？（如新方法还依赖压力） % 第 2 步：清单第 5 列的"参数来源"是否齐全？缺失参数必须先补齐，再试算。 % 第 3 步：清单第 6 列的"适用范围"是否覆盖全部工况？超范围必须改方法或缩工况。 % % --------------------------------------------------------------------------- % 回归基线（铁律 9：改完必验） % 物性方法/参数集一经变更，必须用固定算例复跑并对比下列量： % 基线量：液相组成 x(i)、气相组成 y(i)、两相流量 L 与 V、目标转化率/回收率 % 对比阈值：由你的工程精度要求决定，但"变更物性后结果一致"必须留下记录 % ===========================================================================

逐行剖析：这张清单把"物性依赖"从"散落在方程里的若干调用"变成一张可审查的表。它的三个用途分别对应三类真实场景：换物性方法时的核对表、证据审查时的追溯链（铁律 2：官方源优于第三方）、论文写作时的"物性假设"小节素材。最后一段的回归基线要求，则是把铁律 9 落到物性这个最容易被"悄悄改掉"的位置上——物性变更往往被当成"无关紧要的调整"，而它恰恰是能让所有历史结论失效的那类变更。

代码 2-3：含相平衡的分馏/吸收类模型的物性接入骨架（结构示意）

% --------------------------------------------------------------------------- % 说明：本块为"结构示意"，示范"反应 + 相平衡 + 能量衡算"三者耦合时的物性消费 % 骨架。属性调用、物性方法选择与参数集声明的语法以官方《gPROMS % ModelBuilder Guide》为准。数值为教学示例值（铁律 10）。 % 物理背景：单个平衡级（理论板），进料与两股出料，含一股热负荷 Q。 % --------------------------------------------------------------------------- MODELS PARAMETER N = 2 % 组分数 z(i) % 进料摩尔分数（Σz = 1） F = 1.0 % 进料摩尔流量 P = 1.0e5 % 级压力 Q = 0.0 % 级热负荷（示意 W；0 表示绝热级） VARIABLE x(i) % 液相摩尔分数（出液相） y(i) % 气相摩尔分数（出气相） K(i) % 相平衡常数（物性侧） h_L % 液相摩尔焓（物性侧） h_V % 气相摩尔焓（物性侧） h_F % 进料摩尔焓（物性侧） T % 级温度（由相平衡与能量衡算共同决定） L % 液相流量 V % 气相流量 %% ===== 1) 物性侧量：全部由物性系统按 (T, P, 组成) 给出 ===== %% K(i)、h_L、h_V、h_F 的调用语法以官方《gPROMS ModelBuilder Guide》为准； %% 它们是本模型对物性的全部依赖，已在"物性消费清单"中登记。 %% ===== 2) 相平衡（把物性结果消费进方程）===== y(i) = K(i)*x(i) % E1 相平衡：连接气相与液相组成 %% ===== 3) 逐组分物料衡算 ===== F*z(i) = L*x(i) + V*y(i) % E2 组分 i：进料 = 出液相 + 出气相 %% ===== 4) 能量衡算（焓必须与相平衡来自一致的热力学链条）===== F*h_F + Q = L*h_L + V*h_V % E3 级能量衡算（残差式） %% ===== 5) 归一化 ===== SUM(x(i)) = 1.0 % E4 液相归一 SUM(y(i)) = 1.0 % E5 气相归一（与 E4 中仅一条独立，见模型头） %% ===== 6) 总物料衡算 ===== F = L + V % E6 总摩尔衡算

逐行剖析：

• E3是本块的教学重点：四个焓值全部来自物性侧。若相平衡方法（产生K(i)）与焓模型来自不同的推导链或参数集，E3与E1、E2就会互相矛盾——残差在数学上可以下降，物理上得到的温度与组成却是"焓或相平衡其中一方说了假话"的结果。这就是 §1.3"一致性优先于单点精确"的具体后果。
• T在本模型中不是参数而是变量：闪蒸温度由能量衡算与相平衡共同决定。这解释了为什么物性参数必须带"适用范围"——温度是被解出来的，你无法事先保证它落在参数回归的温区内。
• Q = 0.0表示绝热级；把Q改为非零值即转成"带热负荷的级"，这正是精馏塔逐板计算的单元级原型。结构实验只需改一行（第 06 篇"参数化"思路的延续）。
• 注释里再次出现"以官方《gPROMS ModelBuilder Guide》为准"——这不是啰嗦，而是本系列的双层证据制在每一处涉及闭源语言细节的位置都必须显式落地。

三、常见报错与排查

报错 3-1：相平衡不收敛 / 残差长时间抖动。
现象：迭代过程残差震荡或长时间不下降，尤其出现在高温高压或接近临界区。根因：物性强非线性（相平衡关系在工作点附近陡峭），或初值远离工作点。解法：先做"冻结物性二分法"（把所有强非线性物性临时替换为工作点常数）判断问题在物性侧还是求解器侧；其次给T、x(i)、y(i)设接近预期的初值；最后按第 10 篇的初始化思路分步求解（官方能力名Initialization procedures in gPROMS）。

报错 3-2：摩尔分数越界或出现负值。
现象：求解"成功"，但x(i)或y(i)出现负值或大于 1。根因：相平衡关系与实际工况不匹配（物性方法选错、参数超出回归范围），归一化被强行满足而某一相组分被"挤"成负值。解法：核对物性方法的适用体系与温压范围；检查参数集是否缺项；在物性消费清单里补登"适用范围"并据此缩小工况。

报错 3-3：能量衡算残差始终不降。
现象：物料衡算与相平衡看起来都满足，唯独能量衡算残差居高不下。根因：焓与相平衡来自不一致的热力学链条（例如焓用一套参数、相平衡用另一套）。解法：把焓与相平衡统一到同一物性方法与同一参数集；若必须混用，则退回"等温/绝热简化"版本先跑通，再逐步放开。

报错 3-4：报"方程冗余"或"结构性异常"，但模型其实可解。
现象：诊断信息提示方程数多于未知数，或提示某些方程线性相关。根因：闪蒸/平衡级模型中归一化方程与总物料衡算、逐组分衡算之间存在可导出关系（本篇代码 2-1 已给证明）。解法：在模型头里写明"两条归一化中仅一条独立"；不要盲目删掉归一化——删掉会让模型真正欠定，正确做法是保留并把冗余关系记录在案。

报错 3-5：物性参数缺失或超出适用范围（外推失效）。
现象：报错提示某物性参数不可用，或模型在某个工况下突然失败。根因：该组分/二元对的参数未在参数集中提供，或解出的温度超出参数回归温区。解法：按物性消费清单逐条核对参数完整性；对超范围工况，改为适用范围更宽的方法，或在结论中显式声明物性假设的边界（并禁止用超范围结果做工程判断）。

报错 3-6：求解时间异常增长。
现象：加了物性之后求解时间成倍增加。根因：物性计算被高频调用，且属性本身非线性强。解法：分析阶段先用常物性做趋势判断，正式计算再打开全物性；把"物性简化假设"与"正式计算假设"分别记录（铁律 9：两者结果必须对比留痕）。

四、动手练习

• 练习 1（自由度与冗余）：为代码 2-1 做完整自由度核对，要求写出"归一化冗余"的三行证明（对 N = 2）。判定：证明步骤完整；能说明为什么不能简单删掉 E4；模型求解成功且Σx、Σy均等于 1（数值误差范围内）。
• 练习 2（物性冻结二分法）：把代码 2-1 中的K(i)临时替换为两个常数（取一个合理工作点值），重新求解并与全物性版本对比x(i)、y(i)、L、V。判定：记录两组结果的差异幅度；能据此判断"该体系对相平衡的敏感程度"，并把差异写入工程记录（铁律 9）。
• 练习 3（物性消费清单落地）：为你的课题里的模型，按代码 2-2 格式写一份物性消费清单。判定：至少 5 行；每行都有"依赖变量、参数来源、适用范围"三列且非空；换物性方法的"三步核对动作"能逐条执行。
• 练习 4（范围越界观察）：把代码 2-1 的T从示例值逐步提高到清单中的适用范围之外，观察求解行为变化。判定：能记录到至少一种异常（不收敛、越界、时间暴涨）；能说明"外推失效"与"物性方法选错"在现象上的区别。
• 练习 5（思考题，无标准答案）：为什么"热力学一致性"比"单点物性精确"更优先？请用"残差可以下降但结论错误"这一现象给出两条理由。

五、小结与下一篇预告

本篇把物性从"模型里的常数"还原为"方程组里的一半"：物性是被联立进方程的一批约束，它的失真不会破坏物料衡算的守恒性，只会把解推到"守恒但错误"的位置；模型与物性之间的抽象层（属性调用）让你能"换物性而几乎不改模型"，前提是把模型对物性的全部依赖显式登记；相平衡与焓必须来自一致的热力学链条，一致性优先于单点精确；物性参数必须回答"值多少、从哪来、在什么范围内有效"三问（铁律 10）；闪蒸类模型里"两条归一化仅一条独立"是必须写进模型头的结构性事实（铁律 8）。最后，参数估计时一次只估一类参数——物性参数用独立物性数据先定，动力学参数留给过程数据（第 11 篇）。

第 08 篇进入动态世界：第一个动态模型——从LUMPED到DISTRIBUTED。我们讲清动态与稳态的本质差别（时间自变量、初值、微分方程）、集中参数（LUMPED）与分布参数（DISTRIBUTED）模型的分野（对应官方能力名Implementing a lumped model in gPROMS与Implementing a distributed model in gPROMS）、空间离散与网格/阶数选择对精度与耗时的影响，以及动态仿真对初值的敏感性——那是从"会算稳态"到"会算过程"的分水岭。

本篇认知问题回显（FAQ）

Q1：物性在模型里的位置到底在哪？为什么它不是常数？
A：只有纯组分基础物性（分子量、临界参数、生成焓）近似常数，可放进 PARAMETER 段；密度、黏度、热容随温度乃至组成与压力变化，相平衡常数、逸度系数、活度系数则是 P、T 与全组成的强非线性函数，反应焓也随温度变化。因此物性不是查询表，而是被联立进方程组的约束。它失真的后果不是报错，而是非线性结构悄悄改变，例如唯一稳态变多重稳态、收敛变振荡，所以排查时要先排除物性。

Q2：模型与物性之间的抽象层是如何工作的？为什么能换物性而不改模型？
A：抽象层让模型只做属性调用，即把当前的 T、P、组成作为输入、取回属性值，由物性系统负责具体计算。这带来三个好处：换物性方法时模型一字不改，同一模型跑多种热力学假设成为低成本实验；物性参数集集中维护，避免同一物性在不同模型里取值不同；属性调用点即模型对物性的全部依赖，便于逐条追溯审查。属性的具体名称、物性方法选择、参数集覆盖语法与外部物性接入方式均以官方《gPROMS ModelBuilder Guide》为准。

Q3：相平衡与热力学一致性为什么会直接影响结论？不一致有哪些症状？
A：因为物料衡算无论物性对错都永远成立，物性错误不会被守恒方程拦住，只会把解推向守恒但物理上错的位置。热力学一致性要求相平衡关系与焓、熵等属性相互自洽。不一致的三种典型症状是：摩尔分数越界甚至为负，归一化被强行满足而某组分被挤成负值；能量衡算残差始终不降，因为焓与相平衡来自不同推导链；同一体系在不同初值下收敛到不同解，出现伪多重稳态。

Q4：物性参数与文献数据要遵循什么纪律？
A：按铁律 10，物性参数必须取自可引用文献或官方示例并标注来源，不得凭记忆填数。执行上每个参数都要回答三问：值是多少、从哪来、在什么范围内有效。第三问最常被忽略，因为动态模型里的温度往往是被解出来的变量，无法事先保证落在参数回归温区内，把模型跑到温区之外时所谓不收敛其实是外推失效。此外所有物性依赖应按清单登记，并在物性方法或参数集变更后按铁律 9 用固定算例复跑对比基线。

Q5：物性错误会造成哪些典型失败？参数估计时物性参数怎么处理？
A：典型失败包括相平衡不收敛或残差抖动、摩尔分数越界、能量衡算残差不降、诊断提示方程冗余而模型实际可解、物性参数缺失或超出适用范围、求解时间异常增长。参数估计的处置是三步：先用独立的物性数据（纯组分性质、二元相平衡数据）单独确定物性参数，不放进动力学回归的决策变量集合；再在物性固定的前提下用过程数据回归动力学参数；最后把两步结果放回同一模型，核对能否复现全部数据。原则是一次只估一类参数。

本文涉及的语言细节以所购版本随附的官方 gPROMS 文档为准。