企业官网建设流程全解析

gPROMS二次开发教程（17）：复杂反应网络与自定义反应器单元实战——网络组织、热点飞温与多重稳态

版本声明块

• 工具/软件：gPROMS 桌面建模环境gPROMS ModelBuilder；检索期官方发布锚点gPROMS Process 2022.1.0，适用版本以官方发布说明为准
• 语言/环境：gPROMS 语言（equation-oriented modelling）+ Python 3 + NumPy（离线复算，不调用任何 gPROMS 组件）
• 本文目标：把"论文里的反应动力学"落成"能算、能估参数、能优化、能解释"的复杂反应网络与可复用反应器单元，并说清热点、飞温与多重稳态这三件事为什么必须靠仿真而不是靠经验判断

一句话结论：复杂反应网络的建模秩序是四件套——组分清单 → 反应清单 → 速率式 → 能量耦合，其中"能量耦合"是决定热点与飞温的最后一环；把反应器封装成带PORT的自定义单元后，它才能被流程复用；而多重稳态的存在意味着"同一个反应器、同一组参数，可能停在两个完全不同的工况"，因此开车路径不是操作习惯问题，而是安全性问题——结论必须落在"哪条路径安全、哪条危险"，而不只是给出一条曲线。

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

• Q1：复杂反应网络该怎么"组织"才不乱？四件套分别负责什么，为什么顺序不能颠倒？
• Q2：从"文献动力学"到"可运行模型"的翻译链有哪些环节？每一步最容易引入什么错误？
• Q3：自定义反应器单元（绝热 / 非等温 / 固定床 / 气液）怎么封装？端口设计有哪些专用考虑？
• Q4：热点、飞温、多重稳态是怎么出现的？为什么它们不能靠"经验值"判断？
• Q5：怎么把复杂反应网络与参数估计、动态优化串成一条闭合工作流？工程结论该怎么写才可解释？

一、机制解析

1.1 反应网络的组织四件套：顺序不可颠倒

一个复杂反应网络，本质上是"多个反应在同一个控制体内同时发生，彼此通过组分与温度耦合"。要把它写得不乱，必须按固定顺序落四件套：

① 组分清单 ── 先定义"有哪些物种"，每个物种一个状态量（浓度/分率） │ 纪律：组分名与化学式一一对应，不重不漏 ▼ ② 反应清单 ── 再定义"发生哪些反应"，逐条写出化学计量 │ 纪律：每条反应只写一次；正逆反应要么合并、要么显式并列 ▼ ③ 速率式 ── 每条反应一条速率式（Arrhenius × 浓度依赖） │ 纪律：速率式的"净生成速率"方向必须与化学计量一致 ▼ ④ 能量耦合 ── 把反应热接进能量衡算 纪律：放热项与组分裂解必须来自同一套速率式

为什么顺序不能颠倒：组分没定清就写速率式，会出现"某个组分被生成却没人消耗"的隐性漏洞；反应清单没定清就写能量衡算，会出现"放热算少了/算重了"的静默错误。而这两类漏洞都不会报语法错——它们只会让结果"看起来正常但不对"。

铁律 8（建模先算自由度）：四件套写完后必须立即做自由度核对：待求量个数 = 方程条数。反应网络最容易出的错，就是"多写了一个速率常数或净生成速率，却没给它配方程"。

1.2 从文献动力学到可运行模型：翻译链

铁律 10（数值与参数不臆造）：指前因子、活化能、反应热一律取自可引用文献或官方示例并标注来源。凭记忆填的动力学参数，会让整条反应网络的可信度归零。

1.3 自定义反应器单元的封装：端口设计的专用考虑

把反应器封装成可复用单元（第 09 篇的封装规范）时，除了通用的"输入清单 / 输出清单"，反应器单元还有两个专用考虑：

1. 必须暴露"热"通道：反应器与外界的热交换（夹套、冷却盘管、预热器）是操作与优化最常动的量。若把热负荷藏在内部不暴露，优化器与控制器就"动不到"它（与第 15 篇"耦合点必须可见"同源）。
2. 必须暴露"压降/机理性诊断量"：固定床的压降、热点温度这类量，不暴露就无法在后处理阶段做安全校核。

反应器单元端口设计（反应器专用补充） ┌──────────────────────────────────────────────┐ │ feed IN ：F_in, C_Ain, T_in, P_in │ 入口（含压力） │ prod OUT：F_in, C_A, C_B, T, P_out │ 出口（含出口压力） │ duty OUT：Q（热负荷，供外部读写） │ 热通道（优化/控制的抓手） │ diag OUT：T_max, P_drop │ 诊断量（安全校核用） └──────────────────────────────────────────────┘ 纪律：端口成员要"够用"但不能"泛滥"——下游用不到的量不进端口

1.4 热点、飞温与多重稳态：机理比经验可靠

这三件事都源于能量衡算的"双线"结构：

放热线（产热） G(T)：随温度上升而"指数变陡"（Arrhenius 特征） 移热线（移热） R(T)：随温度上升而"近似线性上升"（传热温差特征） T ▲ │ G(T) ← 指数曲线 │ ／ │ R(T) ——— 直线 │ ／ ＼ │ ／ ＼ │ ／ ＼ └──┴──┬───┬──┴──────────▶ T 稳态1 稳态2 交点 = 稳态 · 一个交点：唯一稳态（安全） · 两个交点：一个低温稳定解 + 一个高温解 ⇒ 多重稳态 · 三个交点：典型的"低稳定 - 中不稳定 - 高稳定" · 飞温：当"产热增量始终大于移热增量"时，温度单调冲高（失控） · 热点：分布参数反应器中温度沿轴向的局部最高点（局部飞温）

关键认知：这两条线的位置由参数决定，而参数有不确定性（第 12 篇），因此"是否有危险稳态"是一个需要扫描的问题，不是一个可以口算的问题。工程师常说的"这个反应器会飞温"，其严格表述是"能量衡算存在高温稳定解，且开车路径可能进入它"。

铁律 4（商业闭源软件的诚实边界）：gPROMS 对多重稳态的数值求解与分支跟踪策略属闭源实现，本文只给机理与判定方法，具体求解设置以官方《gPROMS ModelBuilder Guide》为准。

1.5 闭合工作流：反应网络 → 参数估计 → 动态优化

① 反应网络（本篇） 写物理：组分/反应/速率/能量 │ ▼ ② 参数估计（第 11 篇） 用实验数据反演 k0、Ea、dH │ 官方能力名：Estimating model parameters from experimental data ▼ ③ 动态仿真（第 08 篇） 看动态行为：热点、飞温、稳态切换 │ ▼ ④ 动态优化（第 15 篇） 找安全最优操作：温控曲线 / 开车路径 │ 官方能力名：Dynamic flowsheeting and optimization ▼ ⑤ 参数不确定性（第 12 篇）GSA：稳态是否随参数漂移而改变？ │ 官方能力名：Global System Analysis ▼ ⑥ 工程结论（本篇） "哪条路径安全、裕度多少、风险在哪"

闭环的意义：①→②→③→④→⑤ 跑通一次，你就有了"从动力学到安全操作"的完整证据链；而⑥是这条链的输出——一份能被评审的结论，而不是一堆曲线。

1.6 工程结论的可解释性

一份合格的复杂反应器仿真结论，必须回答三问：“最危险的点在哪（位置+时刻）？”“它离安全边界还剩多少裕度？”“换一条操作路径能不能拉大裕度？”只给"最高温度是 623 K"是不够的——没有位置、没有时刻、没有裕度的结论，不可用于决策。

二、完整代码/模型与逐行剖析

代码 2-1：复杂反应网络模型（绝热，连串-平行反应）

% =========================================================================== % 模型名 : reaction_network_adiabatic 模型类型 : LUMPED（集中参数 / 稳态） % 物理背景 : 绝热反应器中的复杂反应网络（连串-平行） % 反应 R1（主）: A -> B 速率 r1 = k1*C_A 目标产物 B % 反应 R2（副）: A -> C 速率 r2 = k2*C_A 平行副反应，耗 A % 反应 R3（副）: B -> C 速率 r3 = k3*C_B 连串副反应，耗 B % 组织四件套 : 组分清单 / 反应清单 / 速率式 / 能量耦合 % 说明 : "结构示意"。MODELS / PARAMETER / VARIABLE 为已确证关键字；数组与下标的 % 声明语法、方程段的段名以官方《gPROMS ModelBuilder Guide》为准（铁律 4）。 % 本例以"逐组分独立变量"写法呈现，便于逐条核对；组分较多时应改用数组 % 与批量方程（第 06 篇），语法以官方文档为准。 % 数值 : 教学示意值，正式使用须替换为可引用文献值（铁律 10） % =========================================================================== MODELS PARAMETER R = 8.314 % 通用气体常数 [J/(mol*K)]，文献值 k01 = 4.0e6 % R1 指前因子 [1/s]（示意） k02 = 2.0e6 % R2 指前因子 [1/s]（示意） k03 = 1.0e6 % R3 指前因子 [1/s]（示意） Ea1 = 5.0e4 % R1 活化能 [J/mol]（示意） Ea2 = 6.0e4 % R2 活化能 [J/mol]（示意） Ea3 = 7.0e4 % R3 活化能 [J/mol]（示意） dH1 = -8.0e7 % R1 反应热 [J/kmol]，放热取负（示意） dH2 = -1.0e8 % R2 反应热 [J/kmol]（示意） dH3 = -5.0e7 % R3 反应热 [J/kmol]（示意） rho_cp = 2.0e6 % 体积热容 [J/(m^3*K)] tau = 600.0 % 停留时间 [s] C_Ain = 1.0 % A 进料浓度 [kmol/m^3] T_in = 600.0 % 进料温度 [K] VARIABLE C_A % A 浓度（状态量） C_B % B 浓度（状态量，目标产物） C_C % C 浓度（状态量，副产物） T % 反应器温度（状态量） r1 % R1 净速率 [kmol/(m^3*s)] r2 % R2 净速率 r3 % R3 净速率 k1 % R1 速率常数 [1/s] k2 % R2 速率常数 k3 % R3 速率常数 %% ===== ① 组分清单 ===== %% 物种：A、B、C（各一个状态量 C_A / C_B / C_C，不重不漏） %% ===== ② 反应清单 ===== %% R1: A -> B ； R2: A -> C ； R3: B -> C（每条只写一次） %% ===== ③ 速率式 ===== k1 - k01*EXP(-Ea1/(R*T)) = 0 % E1 Arrhenius：R1 速率常数 k2 - k02*EXP(-Ea2/(R*T)) = 0 % E2 Arrhenius：R2 速率常数 k3 - k03*EXP(-Ea3/(R*T)) = 0 % E3 Arrhenius：R3 速率常数 r1 - k1*C_A = 0 % E4 R1 速率（对 A 一级） r2 - k2*C_A = 0 % E5 R2 速率（对 A 一级） r3 - k3*C_B = 0 % E6 R3 速率（对 B 一级） %% ===== 组分衡算：净积累 = 生成 - 消耗（稳态用停留时间形式）===== (C_A - C_Ain)/tau + (r1 + r2) = 0 % E7 A：被 R1、R2 消耗 (C_B - 0.0)/tau - (r1 - r3) = 0 % E8 B：由 R1 生成、被 R3 消耗（B 无进料） (C_C - 0.0)/tau - (r2 + r3) = 0 % E9 C：由 R2、R3 生成（C 无进料） %% ===== ④ 能量耦合（绝热：反应放热全部转为显热升温）===== rho_cp*(T - T_in)/tau - (r1*(-dH1) + r2*(-dH2) + r3*(-dH3)) = 0 %% E10 量纲核对：左右均为 [J/(m^3*s)] %% rho_cp*(T-T_in)/tau : [J/(m^3*K)]*[K]/[s] = [J/(m^3*s)] %% r_i*(-dH_i) : [kmol/(m^3*s)]*[J/kmol] = [J/(m^3*s)] ✓ %% ===== 自由度核对（铁律 8）===== %% 待求量 10 个：C_A,C_B,C_C,T,r1,r2,r3,k1,k2,k3 %% 方程 10 条：E1~E6（6）+ E7~E9（3）+ E10（1） ⇒ 自由度 = 0 ✓

逐行剖析：

• 四件套在代码里是"分区分块"出现的：①组分清单是注释块（对应三个状态量），②反应清单是注释块（对应三条化学计量），③速率式是E1~E6，④能量耦合是E10。这种"分层注释 + 对应的方程块"写法，让"某条反应有没有漏掉放热"变成一眼可查的事——这正是长模型可维护性的落地形态（第 05 篇）。
• E8中 B 的进料浓度写成0.0而不是省略：显式写出"B 无进料"这条假设，比"心算掉"要安全得多。假设显式化，将来换成"B 有循环回料"时只需改这一个数（与第 15 篇"简化假设也要写成方程"同源）。
• E7~E9的净生成速率方向逐条可核：A 被R1+R2消耗、B 由R1生成且被R3消耗、C 由R2+R3生成。方向写反是最隐蔽的错误之一——它不报错，只会给出"转化率很高但收率很低"这种"看起来可能是对的"结果。
• E10的量纲核对写在注释里：三项逐一写出单位并确认一致。这是本篇最重要的肌肉动作——因为能量耦合一旦量纲错位，"热点温度"这个结论就整体失效。
• 自由度核对写在模型末尾：10 个待求量对 10 条方程。把账算平写进注释，比事后在求解器报错里找漏项快得多（铁律 8）。

代码 2-2：自定义固定床反应器单元（分布参数 + 端口化封装，结构示意）

% =========================================================================== % 模型名 : fixed_bed_reactor 模型类型 : DISTRIBUTED（分布参数） % 用途 : 沿床层轴向变化的固定床反应器，封装为可被流程复用的单元（第 09 篇） % 说明 : "结构示意"。UNIT / PORT / PARAMETER / VARIABLE / DISTRIBUTED 为已确证 % 关键字；轴向坐标声明、空间导数记号与离散方式以官方文档为准（铁律 4）。 % 数值 : 教学示意值；压降关联式系数须取自可引用文献（铁律 10） % =========================================================================== MODELS PARAMETER L = 2.0 % 床层长度 [m] D = 0.05 % 床层内径 [m] rho_cp = 2.0e6 % 气相体积热容 [J/(m^3*K)] u = 0.5 % 表观气速 [m/s] dH1 = -8.0e7 % 主反应反应热 [J/kmol]（示意） dH2 = -1.0e8 % 副反应反应热 [J/kmol]（示意） k01 = 4.0e6 % 主反应指前因子 [1/s]（示意） k02 = 2.0e6 % 副反应指前因子 [1/s]（示意） Ea1 = 5.0e4 % 主反应活化能 [J/mol]（示意） Ea2 = 6.0e4 % 副反应活化能 [J/mol]（示意） R = 8.314 % 通用气体常数 [J/(mol*K)] VARIABLE Z % 轴向坐标（记号以官方文档为准） C_A % A 的轴向分布浓度 C_B % B 的轴向分布浓度 T % 轴向分布温度 P % 轴向分布压力 r1 % 主反应速率（轴向分布） r2 % 副反应速率（轴向分布） k1 % 主反应速率常数（轴向分布） k2 % 副反应速率常数（轴向分布） T_max % 床层最高温度（诊断量，供安全校核） P_drop % 床层压降（诊断量，供流程/能耗核算） %% ===== 对外接口：端口（PORT）——反应器单元的专用补充（见 1.3）===== PORT feed IN % 入口端口 F_in, C_Ain, T_in, P_in % 含压力：分布参数模型的边界条件之一 PORT prod OUT % 出口端口 F_in, C_A, C_B, T, P_out% F_in 亦在此 => 表达"流量连续" PORT duty OUT % 热通道端口：供外部读写热负荷（优化/控制抓手） Q % 夹套/盘管换热量 PORT diag OUT % 诊断量端口：供后处理做安全校核 T_max, P_drop % 热点温度、压降 %% ===== 单元内部物理：分布参数衡算（沿 Z）===== k1 - k01*EXP(-Ea1/(R*T)) = 0 % 本构：主反应 Arrhenius k2 - k02*EXP(-Ea2/(R*T)) = 0 % 本构：副反应 Arrhenius r1 - k1*C_A = 0 % 本构：主反应速率 r2 - k2*C_A = 0 % 本构：副反应速率 %% 轴向物料衡算（稳态、以"对流项 + 净生成项"形式书写；记号以官方为准） u*( 轴向一阶导数(C_A) ) + (r1 + r2) = 0 % A 沿程被消耗 u*( 轴向一阶导数(C_B) ) - r1 = 0 % B 沿程由主反应生成 %% 轴向能量衡算（对流项 = 反应放热项 + 外部换热项） rho_cp*u*( 轴向一阶导数(T) ) - (r1*(-dH1) + r2*(-dH2)) + Q/L = 0 %% 轴向压降（经验关联式）；系数须取自可引用文献（铁律 10） 轴向一阶导数(P) + ( 经验阻力系数 )*u*u / D = 0 %% 诊断量定义 T_max - ( 沿 Z 的 T 的最大值 ) = 0 % 热点温度（供安全校核） P_drop - ( P_in - P_out ) = 0 % 压降（供能耗核算） %% 说明：空间导数的具体记号与离散方式以官方文档为准；此处仅表达"物理身份"

逐行剖析：

• 端口设计体现了 1.3 的两条专用考虑：duty端口暴露热负荷Q，diag端口暴露T_max与P_drop。前者让优化器能"拧"热量，后者让报告能"报"热点与压降——两者缺一，单元就只能"算自己"，无法支撑决策。
• 空间导数项写成"轴向一阶导数(量)"的中文占位：这是本系列的诚实姿态——具体记号（如是否用PARTIAL、坐标声明方式）属闭源细节，一律以官方文档为准（铁律 4）。宁可留一个明确标注的占位，也不臆造一个看起来对的语法。
• 能量衡算里"外部换热项Q/L"显式出现：把热通道接进物理，"联用/优化/控制能动的量"与"方程里的量"就对齐了（与第 16 篇"Q 必须暴露"是同一条纪律）。
• T_max与P_drop被定义为"诊断量"而非内部中间量：定义成可寻址的名字，报告的"最高温度"才有唯一口径——否则每个人从曲线上肉眼读值，结论会不一致。
• 压降关联式标注"系数须取自可引用文献"：直接落地铁律 10。固定床压降系数是典型的"凭记忆填数就错"的地方。

代码 2-3：热点与多重稳态分析脚本（纯离线，可完整运行）

# -*- coding: utf-8 -*-""" 复杂反应网络的热点与多重稳态分析（纯离线，不调用任何 gPROMS 组件）。 两件事： (A) 多重稳态判定：在绝热稳态能量衡算的"双线结构"上找交点； f(T) = 放热升温(T) - 实际升温(T)，f 的符号变化次数 = 稳态个数。 (B) 热点搜索：沿床层轴向用显式步进做一次"温度-转化"演化，找局部最高点。 纪律： · 参数与模型（代码 2-1 / 2-2）保持同一套，便于与 gPROMS 结果对照（铁律 9）； · 本脚本是"机理判读器"，用于在提交 gPROMS 求解前先建立数量级直觉。 """importnumpyasnp# --- 与代码 2-1 一致的参数（复核必须回到原模型，铁律 9） ----------------------R_GAS=8.314K01,K02=4.0e6,2.0e6# 指前因子 [1/s]（示意值，铁律 10）EA1,EA2=5.0e4,6.0e4# 活化能 [J/mol]（示意值）DH1,DH2=-8.0e7,-1.0e8# 反应热 [J/kmol]（放热取负）RHO_CP=2.0e6# 体积热容 [J/(m^3*K)]TAU=600.0# 停留时间 [s]C_AIN=1.0# 进料浓度 [kmol/m^3]T_IN=600.0# 进料温度 [K]defsteady_state_at(T):"""给定温度 T，解出该温度下"绝热稳态"所需的实际升温与放热升温。 只保留两条反应（R1、R2 耗 A），用解析式以便快速扫描。"""k1=K01*np.exp(-EA1/(R_GAS*T))k2=K02*np.exp(-EA2/(R_GAS*T))# 稳态浓度（CSTR 型解析解，A 被两条反应消耗）C_A=C_AIN/(1.0+TAU*(k1+k2))r1,r2=k1*C_A,k2*C_A# 放热升温：由反应放热折算成温升dT_gen=TAU*(r1*(-DH1)+r2*(-DH2))/RHO_CPreturndT_gen,C_Adeffind_steady_states(T_lo=550.0,T_hi=900.0,n=4001):"""扫描 f(T) = 放热升温(T) - (T - T_in)，用符号变化定位所有稳态。 返回交点温度列表：长度 1 = 唯一稳态；>=2 = 多重稳态。"""Ts=np.linspace(T_lo,T_hi,n)f=np.array([steady_state_at(T)[0]-(T-T_IN)forTinTs])roots=[]foriinrange(len(Ts)-1):iff[i]==0.0:roots.append(Ts[i])eliff[i]*f[i+1]<0.0:# 符号变化 => 有交点# 线性内插给出交点温度（足够用于判读）t=Ts[i]-f[i]*(Ts[i+1]-Ts[i])/(f[i+1]-f[i])roots.append(round(float(t),3))returnrootsdefhot_spot_profile(T0=T_IN,n=200,tau_total=TAU):"""沿"停留进度"步进，给出温度-转化轨迹，返回最高温与其位置。 这是固定床轴向热点的一维类比（代码 2-2 的简化版），用于数量级判读。"""taus=np.linspace(0.0,tau_total,n)# 进度坐标（类比轴向位置）T=T0 dtau=tau_total/(n-1)T_hist=[T]for_inrange(n-1):k1=K01*np.exp(-EA1/(R_GAS*T))k2=K02*np.exp(-EA2/(R_GAS*T))C_A=C_AIN/(1.0+T*(k1+k2))# 局部近似浓度dT=(k1*C_A*(-DH1)+k2*C_A*(-DH2))/RHO_CP*dtau T=T+dT# 放热升温（绝热）T_hist.append(T)T_hist=np.array(T_hist)i_max=int(np.argmax(T_hist))returntaus,T_hist,float(T_hist[i_max]),float(taus[i_max])if__name__=="__main__":# ---- (A) 多重稳态判定 ----roots=find_steady_states()print("=== (A) 稳态交点扫描 ===")fori,Tinenumerate(roots,start=1):print(f" 稳态{i}: T ={T:.2f}K")iflen(roots)>=2:print(f" 判读：存在{len(roots)}个稳态 ⇒ 多重稳态。"f"最低温解与最高温解都可能被实际路径'落到'。")eliflen(roots)==1:print(" 判读：唯一稳态，未发现危险的高温分支（仅限本参数下）。")else:print(" 判读：扫描区间内无交点，检查参数或扫描范围。")# ---- (B) 热点搜索 ----taus,T_hist,T_max,tau_at_max=hot_spot_profile()print("\n=== (B) 热点搜索（一维类比）===")print(f" 最高温度 T_max ={T_max:.2f}K")print(f" 出现在进度 tau ={tau_at_max:.1f}s（占总停留 "f"{100.0*tau_at_max/TAU:.1f}%）")print(f" 终点温度 T_end ={T_hist[-1]:.2f}K")print(" 判读：热点在'前段升温迅速、后段趋于平缓'时成立；""若 T_max 逼近材料上限，需在代码 2-2 中显式加路径约束。")

逐行剖析：

• find_steady_states用"符号变化 + 线性内插"定位交点：这就是 1.4 节"双线结构"的数值实现。它的价值在于把"是否存在危险稳态"从口号变成一次扫描——交点个数 ≥2 就是多重稳态的硬证据。这类判读应在提交 gPROMS 求解之前先做，用来建立数量级直觉。
• steady_state_at用 CSTR 型解析式求局部浓度：这是简化假设，其作用是"快速扫描"，不是"替代 gPROMS"。在注释里写明它是简化，避免读者误以为它等价于代码 2-1 的完整模型（这也是规避"两套模型结果不一致"的正规做法，第 14 篇）。
• hot_spot_profile用"停留进度"作为轴向坐标的类比：把分布参数问题降成一维步进，用于回答"热点大概在哪、有多高"。真实结论仍须由代码 2-2 的分布参数模型给出——脚本只是"提交求解前的侦察"。
• 输出同时给出"最高温度"与"出现位置/占比"：直接落地 1.6 节的可解释性要求。没有位置与裕度的最高温度，不是一条可用结论。

三、常见报错与排查

报错 3-1：转化率很高，但目标产物收率很低。
现象：结果"看起来正常"却不合理。根因：净生成速率方向写反（如把 B 的消耗写成了生成），或漏掉了连串副反应 R3。解法：回到代码 2-1 的②反应清单，逐条核对每个组分的"生成项 / 消耗项"是否齐全、方向是否正确（铁律 8）。

报错 3-2：热点温度算得异常高（或异常低）。
现象：温度明显偏离物理直觉。根因：能量耦合量纲错位，或某个反应的放热项漏乘化学计量数。解法：逐项做量纲走查（代码 2-1 的E10注释格式）；核对每条反应的放热是否都进入了能量衡算（铁律 10）。

报错 3-3：换个初值结果跳到另一个稳态。
现象：同一模型两套结果。根因：多重稳态（1.4 节）。解法：用代码 2-3 的扫描确认真实稳态个数；把"低稳定解 / 高稳定解"分别记录并给出各自的开车路径，不要只报一个解。

报错 3-4：求解器报"欠定"或存在自由变量。
现象：模型不收敛。根因：多写了一个速率常数或净生成速率却没配方程（四件套顺序颠倒的典型后果）。解法：重列待求量清单与方程清单，核对条数一致。

报错 3-5：分布参数模型报空间导数相关错误。
现象：与轴向坐标或导数记号有关的报错。根因：空间坐标/导数记号写法与版本不符（属闭源细节）。解法：以官方《gPROMS ModelBuilder Guide》核对坐标声明与离散设置；不得凭推断补全语法（铁律 4）。参考官方能力名Implementing a distributed model in gPROMS的对应资源。

铁律 9（改完必验）：反应网络的任何参数调整（尤其活化能与反应热）都必须复跑固定算例，并与回归基线比对。反应网络对活化能极其敏感——Ea 差一点，热点差几十开。

四、动手练习

• 练习 1（四件套落地）：把一篇文献里的三反应网络翻译成代码 2-1 格式的完整模型。判定：组分清单、反应清单、速率式、能量耦合四块齐全；自由度核对通过；所有参数标注来源。
• 练习 2（方向核对）：把代码 2-1 中E8的符号故意改错，重跑。判定：观察到"转化率高、收率低"的现象；能说明净生成速率方向错误为何不报错却给错结果；随后改回并复跑。
• 练习 3（多重稳态）：运行代码 2-3 的 (A) 段，分别用两组不同的活化能。判定：给出交点个数并判读是否存在多重稳态；指出"哪一组参数下存在危险的高温分支"。
• 练习 4（热点定位）：运行代码 2-3 的 (B) 段。判定：给出最高温度、出现位置与占比；在代码 2-2 中补一条T ≤ T_limit的路径约束（第 15 篇），说明它如何改变最优操作。
• 练习 5（封装检查）：为你的反应器写一份代码 2-2 格式的端口清单。判定：入口/出口/热通道/诊断量四类端口齐全；每个端口成员都能说明"谁给它赋值、谁读它"；下游用不到的量未进端口。

五、小结与下一篇预告

本篇把"文献动力学"变成"可运行、可复用、可解释"的资产：反应网络的组织秩序是四件套（组分清单 → 反应清单 → 速率式 → 能量耦合），顺序颠倒会埋下"不报错却算错"的隐性漏洞；翻译链上最容易出错的是单位、方向、漏项三处；反应器封装在通用端口之外还要暴露热通道与诊断量；热点、飞温、多重稳态都源于能量衡算的**“双线结构”——它们是需要扫描的问题，不是可以口算的问题；最后把反应网络 → 参数估计 → 动态仿真 → 动态优化 → 参数不确定性 → 工程结论**串成一条闭合工作流，而结论必须回答"危险点在哪、裕度多少、换路径能否拉大裕度"。

第 18 篇《动态过程与批次/开停车场景实战》把"有情节的过程"讲透：Batch flowsheeting与连续动态流程的本质差别、TASK的PARAMETER/VARIABLE/SCHEDULE三段如何分工、事件与状态切换怎么表达（以官方文档为准）、开停车过渡过程为什么必须专门仿真，以及批次模型如何用"参数化配方"实现可复用化。

本篇认知问题回显（FAQ）

Q1：复杂反应网络该怎么组织才不乱？四件套为什么顺序不能颠倒？
A：按固定顺序落四件套：先组分清单定义有哪些物种、每个物种一个状态量；再反应清单逐条写出化学计量；再对每条反应写一条速率式；最后把反应热接进能量衡算。顺序不能颠倒的原因是，组分没定清就写速率式会出现某组分被生成却无人消耗的隐性漏洞，反应清单没定清就写能量衡算会出现放热算少或算重的静默错误，这两类漏洞都不报语法错，只会给出看起来正常但错误的结果。

Q2：从文献动力学到可运行模型的翻译链有哪些环节，哪一步最容易出错？
A：五个环节：抄写文献的反应式与速率式，对齐参数并统一单位，补全各组分净生成速率，把速率式耦合进能量衡算的放热项，最后做自由度核对。最容易出错的是单位与方向与漏项三处：论文常用 mol/L/s 而模型用 kmol/m3/s，活化能存在 J/mol 与 cal/mol 混用，反应热符号放热须取负，净生成速率方向容易写反，放热项容易漏掉某个反应或漏乘化学计量数。

Q3：自定义反应器单元怎么封装？端口设计有哪些专用考虑？
A：在通用封装规范之外，反应器单元有两条专用考虑。必须暴露热通道，即夹套或盘管或预热器的热负荷，因为它是操作与优化最常动的量，藏在内部会让优化器与控制器动不到它。必须暴露压降与热点温度这类机理诊断量，否则无法在后处理阶段做安全校核。典型端口集合为入口端口含压力、出口端口含出口压力与流量、热通道端口暴露热负荷、诊断量端口暴露最高温度与压降。

Q4：热点、飞温与多重稳态怎么出现？为什么不能靠经验判断？
A：三者都源于能量衡算的双线结构：产热随温度呈指数上升，移热随温度近似线性上升。交点即稳态，一个交点表示唯一稳态，两个或三个交点表示多重稳态，其中高温解可能被实际路径落到；当产热增量始终大于移热增量时温度单调冲高即飞温；分布参数反应器中温度沿轴向的局部最高点即热点。不能靠经验判断的原因是这两条线的位置由参数决定，而参数存在不确定性，因此是否有危险稳态属于需要扫描的问题，一个严谨表述是能量衡算存在高温稳定解且开车路径可能进入它。

Q5：怎么把反应网络与参数估计、动态优化串成闭合工作流？结论怎么写才可解释？
A：工作流为六步：先写反应网络的组分与反应与速率与能量耦合，再用实验数据反演指前因子与活化能与反应热，再做动态仿真观察热点飞温与稳态切换，再做动态优化寻找安全最优操作或开车路径，再做全局系统分析判断稳态是否随参数漂移而改变，最后输出工程结论。一份合格结论必须回答三问：最危险的点在哪即位置与时刻，它离安全边界还剩多少裕度，换一条操作路径能否拉大裕度。只给一个最高温度数值而没有位置、时刻与裕度的结论不可用于决策。

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