工业蒸汽甲烷重整装置催化剂失活的混合现象学-经验建模研究

《Processes》：Hybrid Phenomenological-Empirical Modeling of Catalyst Deactivation in an Industrial Steam Methane Reformer

　　

工业蒸汽甲烷重整（steam methane reforming, SMR）仍是大规模制氢的主导路线，然而长期催化剂失活与管壁温度演变很少被同时纳入工业反应器模型中。本研究提出了一种针对工业顶燃式蒸汽甲烷重整炉的混合现象学-经验模型（hybrid phenomenological–empirical model），该模型将一维拟均相反应器模型（one-dimensional pseudo-homogeneous reactor model）、具有残余活性的催化剂失活模型（catalyst deactivation model with residual activity）以及基于常规测量过程变量的经验管壁温度模型（empirical wall-temperature model）相结合。模型参数基于超过1000天的工业运行数据，采用混合全局-局部优化策略（hybrid global–local optimization strategy）同时估计，随后使用残差自助法（residual bootstrap）评估参数不确定性。所提出的框架准确复现了甲烷转化率和重整炉出口温度，平均绝对百分比误差（mean absolute percentage error, MAPE）低于5%，出口温度预测无显著系统性偏差，甲烷转化率在整个运行周期内仅有?2.00%的适度系统性低估。估计的失活参数与文献报道的Ni基蒸汽重整催化剂烧结相关老化机制一致。尽管模型复杂度有所降低，该模型成功捕捉了近三年工业运行中催化剂活性与炉膛热行为耦合演变的特征。所提出的方法为工业蒸汽甲烷重整炉的长期过程监测、催化剂健康评估、运行周期分析和预测性维护提供了一种实用且计算高效的分析框架。

论文解读：工业蒸汽甲烷重整装置催化剂失活与炉膛热行为的耦合建模研究

一、研究背景与问题提出

氢能作为化学工业的重要原料和清洁能源载体，其大规模生产主要依赖天然气蒸汽甲烷重整（steam methane reforming, SMR）路线。目前全球约80%的氢气产量来自化石资源，其中天然气占比接近60%，SMR技术凭借其技术成熟度和经济竞争力（生产成本约1–2.5 USD kg?1）成为当前全球氢气供应的主导途径。在工业SMR装置中，脱硫后的天然气与蒸汽混合后进入重整炉，在镍基催化剂作用下于高温（750–900 °C）和中压（20–30 bar）条件下发生强吸热的甲烷重整反应，生成富含氢气和一氧化碳的合成气，随后经水煤气变换（water–gas shift, WGS）和变压吸附（pressure swing adsorption, PSA）等工序获得高纯氢气。

然而，工业重整炉的运行涉及复杂的多物理场耦合：化学反应动力学、炉膛辐射与对流传热、颗粒内扩散以及沿催化剂填充管的轴向和径向温度梯度等强相互作用现象交织在一起。更为关键的是，在长达数年的工业运行周期中，催化剂活性因烧结、积碳和中毒等机制而持续衰减，这些失活效应随时间改变有效反应速率和沿管热分布，使SMR构成一个高度非线性和时变的反应系统。

尽管蒸汽重整过程的数学建模在过去几十年取得了显著进展，但现有反应器模型大多假设催化剂活性恒定，未显式表征工业长期运行中催化性能的渐进损失。少数考虑催化剂失活的研究也存在明显局限：Taji等人开发的一维非均相重整炉模型采用简单一阶失活规律，但炉膛热行为和管壁温度未纳入参数估计过程；Jokar等人提出的动态现象学模型采用幂律失活表达式，但其研究对象为膜重整炉，管壁温度作为已知边界条件而非估计变量，且其失活模型假设催化剂最终完全失活，未考虑工业催化剂常见的残余活性。此外，多维和多尺度反应器模型虽具有较高的现象学严谨性，但计算需求巨大，难以用于常规工业应用。

因此，现有研究在利用常规工业运行数据同时估计催化剂活性和炉膛热行为方面存在明显空白。催化剂活性通常被假设为常数或采用简化失活模型描述，管壁温度则直接取自工厂测量值，催化剂老化与炉膛热响应在长期工业运行中的耦合演变尚未得到充分研究。

二、研究内容与核心贡献

针对上述空白，本研究提出了一种混合现象学-经验建模框架（hybrid phenomenological–empirical modeling framework），将催化剂活性和管壁温度从超过1000天的工业运行数据中同时估计。该框架结合一维拟均相重整炉模型、具有残余活性的广义幂律失活模型（generalized power law expression, GPLE）和经验管壁温度关联式，在参数估计全过程中保持热-动力学耦合。与以往研究不同，该框架无需实测管壁温度作为模型输入，仅依赖常规监测的工厂变量，为长期性能评估、催化剂状态监测和工业运行周期分析提供了计算高效且保持工程预测能力的解决方案。

三、关键技术与方法

研究人员采用的核心方法包括：一维拟均相准稳态反应器模型，基于Xu和Froment提出的Langmuir–Hinshelwood动力学表达式描述蒸汽重整和WGS反应；具有残余活性的广义幂律催化剂失活模型，失活速率常数采用以参考温度重参数化的Arrhenius形式；经验管壁温度模型，通过对燃料气流量和空燃比进行一阶指数平滑构造累积变量，以表征重整炉的高热惯性特征；加权最小二乘目标函数，采用粒子群优化（particle swarm optimization, PSO）与内点法（interior-point algorithm）相结合的混合全局-局部优化策略进行参数估计；残差自助法（residual bootstrap）用于参数不确定性评估。工业数据来自某顶燃式工业蒸汽甲烷重整装置，涵盖超过1000天的运行记录。

四、研究结果

4.1 新鲜催化剂条件下模型验证

研究人员在假设新鲜催化剂（活性a=1）且采用光学高温计实测管壁温度的条件下，独立评估了现象学模型的预测能力。结果表明，甲烷转化率的最大偏差为4.15%，出口压力偏差为0.96%，出口温度偏差为0.46%。这一验证确认了现象学模型在新鲜催化剂条件下能够捕捉重整炉的整体行为，为后续纳入催化剂失活和管壁温度子模型提供了可靠基础。

4.2 时间序列模型拟合

研究人员将完整模型拟合至整个工业运行周期的数据，结果显示模型成功复现了甲烷转化率和重整炉出口温度在1000余天运行中的时间演变，既包括运行波动也包括长期趋势。甲烷转化率的平均绝对百分比误差（MAPE）为4.79%，偏差为?2.00%；出口温度的MAPE为0.35%，偏差为1.25 °C。值得注意的是，该工业装置通过调整操作条件将甲烷转化率维持在平均值附近的窄区间内以补偿催化剂活性损失，这一行为被模型成功复现。出口各组分摩尔分数的预测与测量也呈现良好一致性，表明所采用的动力学方案在催化剂活性递减条件下仍能一致捕捉重整与WGS反应之间的耦合。

4.3 催化剂活性模型参数估计

研究人员从工业数据中同时估计了催化剂失活模型和管壁温度模型的参数。估计的失活活化能（E_d = 127 kJ mol?1）落在文献报道的Ni/α-Al₂O₃催化剂在富氢气氛中烧结失活的活化能范围（31–159 kJ mol?1）内，与工业重整炉报道值（93和94 kJ mol?1）以及Ni/α-Al₂O₃水热烧结活化能（126 kJ mol?1）和Ni–OH物种表面迁移活化能（137 kJ mol?1）高度吻合。失活速率常数在参考温度下为k_d0 = 2.26 × 10?? d?1，与文献报道值同数量级。失活级数m = 1.18接近1，在GPLE框架内这一数值与Ostwald熟化机制相关联。残余活性a_s = 0.42表明催化剂在长期运行后仍保留约42%的活性。预测的催化剂活性曲线显示，运行初期数月内活性从1快速下降至约0.5，随后以递减速率趋近残余活性水平，这一行为与新鲜镍催化剂初始小晶粒快速烧结、随后因晶粒增大而烧结速率降低的机理一致。

4.4 管壁温度模型参数估计

管壁温度模型估计的特征时间常数τ = 27天，证实了重整炉与炉膛耐火衬里、催化剂床层和循环气体相关的热惯性假设。基线管壁温度参数β₀ = 1051.80 °C代表炉膛平均热状态。标准化效应分析表明，平滑空燃比对管壁温度的敏感性最高（7.43 °C），其次是进料流量（4.07 °C）和累积燃料输入（3.45 °C）。模型预测的管壁温度在整个运行周期的大部分时间内落在光学高温计测量的标准差带内，表明模型能够捕捉平均管壁温度趋势。研究还揭示了两个不同的运行阶段：运行初期管壁温度逐渐上升至近似稳态水平，此后保持恒定并仅有小幅波动。

4.5 参数相关性分析

研究人员通过残差自助法获得了参数相关矩阵。Arrhenius重参数化显著降低了失活指前因子与活化能之间的相关性（相关系数仅0.03），表明以参考温度为中心的参数化策略有效改善了参数可辨识性。失活级数m与k_d0之间存在中等负相关（?0.53），这是非线性动力学模型中常见的参数补偿效应。残余活性a_s与管壁温度模型参数β₃之间存在较强交叉相关（0.70），反映了两个子模型之间固有的物理耦合——管壁温度强烈影响失活动力学，而催化剂活性损失又改变维持重整炉性能所需的热负荷。

五、讨论与结论

研究人员指出，所提出的混合框架成功实现了催化剂活性与炉膛热行为的同步估计，避免了热参数与动力学参数之间的补偿效应，从而获得了更符合物理实际的长期重整过程描述。管壁温度模型的建立使得从常规监测变量连续估计重整炉热状态成为可能，补充了光学高温计测量在时间连续性上的不足。估计的催化剂活性可作为催化剂状态的宏观指标，服务于催化剂老化评估、运行周期分析、性能监测和维护规划。

该框架作为详细多维炉膛和反应器模型的计算高效替代方案，仍存在若干局限：催化剂活性以集总变量表示，未区分具体失活机制；管壁温度关联式具有经验性，依赖参数估计所覆盖的运行范围；失活活化能因工业运行温度范围较窄而仅具中等可辨识性；反应器模型采用一维拟均相准稳态假设，未解析径向梯度、颗粒内效应或快速瞬态动力学。

尽管如此，参数估计、不确定性分析和模型性能均表明，所采用的模型复杂度足以捕捉长期重整炉运行中催化剂老化和热漂移时间尺度上的主导相互作用。该研究弥合了详细反应器建模与长期工业催化剂管理之间的差距，为工业蒸汽甲烷重整炉的在线监测、催化剂更换规划、运行优化以及向实时过程控制和数字孪生应用扩展奠定了方法论基础。