推进湿氧化的广义动力学模型及其在水热液化水相中的应用

《Journal of Water Process Engineering》：Advancing the generalized kinetic model for wet oxidation and its application to hydrothermal liquefaction aqueous phases

　　

研究人员提出了一种精简的动力学模型，用于描述水热液化（HTL）水相（APs）的湿式氧化（WO）行为。该模型基于实验数据进行参数化，并在Aspen Plus^?中实现。研究人员预测了两种HTL水相的关键水质指标（如化学需氧量（COD）和总有机碳（TOC））以及热释放。该模型成功捕捉了不同组成水相中TOC和COD的变化趋势。该模型支持集成HTL-WO精炼概念的系统级评估与优化。

**论文解读文章**

**研究背景与问题**
水热液化（HTL）是一种将湿生物质转化为生物原油的技术，该生物原油可进一步升级为运输燃料（包括航空燃料）。然而，连续HTL工艺通常需要处理大量水，因为进料浆料的干物质（DM）含量受限于10–20 wt%。由此产生的水相（APs）含有大量可溶性有机物，既需要处理后才能排放，也为回收有价值的副产品提供了机会。目前讨论的三种主要水处理技术包括催化水热气化（cHTG）、厌氧消化（AD）和湿式氧化（WO）。cHTG虽能有效转化有机污染物，但操作温度高（350–700°C）、压力大（20–35 MPa），且催化剂易失活，难以在宽原料范围内可靠运行。AD是一种成熟且节能的技术，但HTL-AP中高浓度的抑制性化合物（如酚类、铵）需要复杂的预处理策略，尚无法可靠处理多种原料。相比之下，WO通常在中等温度（150–370°C）和高压（50–200 bar）下使用氧气、空气或过氧化氢作为氧化剂，能有效降低毒性，且可通过调节工艺条件和氧化剂供应来控制转化程度，例如设计耦合HTL-WO反应器方案中的热释放。尽管WO的单一反应机理已有广泛研究，但实际应用中常涉及复杂混合物（如废水或HTL-AP），其多样组成使得单独建模每种底物的消耗不切实际。因此，需要一种通用动力学模型（GKM）来描述整体氧化动力学。Li等人开发的GKM基于简化反应方案，考虑了中间产物的生成与分解，但尚未针对HTL-AP进行验证，也未利用乙酸浓度数据进行模型验证。

**研究内容与结论**
本研究在Li等人GKM的基础上进行了多项关键改进。首先，将底物从废水扩展到HTL-AP，包括用于模型开发的污水污泥来源的AP，以及两种来自牛粪与小麦秸秆（50:50混合）的AP，用于验证模型的通用性和预测能力。其次，将模型预测的物质属性从总有机碳（TOC）和化学需氧量（COD）扩展到乙酸（AA）和总氮（TN）浓度。第三，预测了WO放热过程的热释放量。最后，将改进后的GKM在Aspen Plus^?中实现，便于未来集成到耦合HTL-WO绿色精炼框架中，支持过程优化及生命周期评估（LCA）和技术经济评估（TEA）。研究结果表明，该模型能有效描述TOC、COD及热释放的变化趋势，对乙酸和铵浓度的预测精度较低但总体趋势合理。通过参数拟合并用两个独立样本验证，模型在不同组成的水相中表现出良好的预测能力，支持系统级热集成分析。该论文发表在《Journal of Water Process Engineering》。

**关键技术方法（不超过250字）**
研究人员采用的主要关键技术方法包括：（1）连续湿式氧化实验，在实验室规模连续活塞流反应器中进行，HTL-AP样品来自污水污泥（模型开发）以及牛粪/小麦秸秆混合物（验证）；（2）基于Peng-Robinson（PR）状态方程（EoS）描述高压气液相平衡（VLE），并回归了水-氧和水-乙酸二元交互参数（k_ij）；（3）建立通用动力学模型（GKM），将有机质（OM）和乙酸（AA）作为主要反应物，采用一级反应动力学（氧气为过量供应，零级反应）；（4）在Aspen Plus^?中实现过程模型，模拟包括预热、反应、冷却和分离的完整流程；（5）参数回归使用Britt-Luecke算法和最小二乘法优化，模型验证通过三个独立HTL-AP样本（HTL-AP PF、HTL-AP P1、HTL-AP P2）进行，并使用均方根误差（NRMSE）、决定系数（R²）和平均绝对误差（MAE）进行定量评估。

**研究结果**
**3.1 二元交互参数k_ij**
通过回归水-氧和水-乙酸系统的二元交互参数，改进了PR状态方程对典型WO条件下非理想行为的描述。水-氧系统的k_ij为0.533（无温度依赖性），水-乙酸系统的k_ij = -0.248 + 2.3×10^-4·T/K，R²为0.93。引入这些参数后，模型更准确地预测了乙酸在高温下的气液相平衡，避免了乙酸被错误分配至气相。

**3.2 动力学参数拟合**
通过优化确定了三个模型反应的动力学参数：反应1（OM完全氧化）的活化能E_a,1 = 66.0 kJ mol^-1，反应2（OM部分氧化为AA）的E_a,2 = 75.0 kJ mol^-1，反应3（AA氧化）的E_a,3 = 168.0 kJ mol^-1。这些值与文献报道一致，尤其AA的活化能与Foussard等人报道的167.7 kJ mol^-1高度吻合。参数稳健性分析显示，OM转化反应的参数识别性较低，而AA氧化反应的活化能约束较窄。模型对TOC的预测（平均R² = 0.87，平均NRMSE = 12%）优于对AA的预测（平均R² = 0.41，平均NRMSE = 29%）。

**3.3 模型验证**
利用参数拟合后未调整的模型，对HTL-AP PF样本的COD和热释放进行验证，并预测了HTL-AP P1和HTL-AP P2的TOC、AA、COD、铵浓度和热释放。结果表明，模型对COD和TOC的预测与实验值高度一致（平均NRMSE约17%，平均R² = 0.74）。热释放预测的偏差主要源于归一化热释放值的差异（模拟值-394 kJ mol^-1 O₂ vs 文献值-435 kJ mol^-1 O₂）。对于预测样本，HTL-AP P1的AA浓度预测较差（R² = -1.03），但第二样本HTL-AP P2的整体预测性能更好（平均NRMSE = 15%，平均R² = 0.78）。铵浓度在中等停留时间出现系统高估，提示存在难氧化的氮中间体。

**3.4 与现有WO动力学模型的比较**
直接定量比较困难，因为大多数WO动力学研究未报告标准化性能指标。Schuck等人针对HTL工艺水的模型（含13种代表性化合物）在高温下表现良好（绝对误差<8%），但本模型结构更精简，且评估了更宽的操作条件范围。与Kindsigo等人对木质素溶液的模型（R² = 0.93）相比，本模型处理的HTL-AP基质更复杂。

**3.5 模型局限性与未来改进**
当前模型在铵浓度和AA浓度预测中存在局限：铵的高估表明有机氮向铵的转化路径过于直接，可能忽略了难氧化的含氮中间体；AA的预测低估表明有机质不能始终用单一团块（lumped）表示，其组成差异导致AA生成与降解行为不同。未来改进方向包括引入额外的难氧化氮中间体，以及将有机质分为多个子团块以区分反应活性。同时，需要更多不同原料和操作条件下的HTL-AP样本来扩展验证基础。

**结论部分翻译**
本工作提出了一种简单而有效的动力学模型，用于描述HTL水相的湿式氧化行为，该模型采用最小参数化，并在组成不同的样品上进行了验证。视觉趋势分析与定量性能指标（R²和NRMSE）的结合为评估模型保真度和指导未来改进提供了稳健基础。模型以显著的一致性捕捉了关键水质指标（COD、TOC）以及热释放的降解趋势，表明其在不同进料组成和废水来源下具有预测能力和稳健性。然而，更具体的参数（如乙酸和铵浓度变化曲线）揭示了在捕捉中间体形成和氧化动力学方面的局限性。本研究的主要结论如下：
- 该模型仅需简单的输入参数（TOC、AA、NH₄⁺、TN和COD）以及六个动力学参数（用于三个反应）即可描述WO化学，使其高效且实用。
- 模型能够高精度预测新HTL-AP的TOC、COD和热释放（平均R² ~ 0.86，平均NRMSE ~ 15%），但对AA和铵浓度的预测精度较低。
- 提高模型对AA和铵浓度预测的灵活性需要更深入地理解底层化学机理、额外的输入参数以及显著增加的模型复杂度。
- 铵形成的实验数据表明存在难氧化的氮物种，这与含氮芳香族化合物（如吡啶）的高活化能观察结果一致，但当前模型无法捕捉。
- 两个预测样本之间AA生成和降解的显著差异可能源于HTL-AP的组成差异，这些差异改变了反应环境和主导路径，导致生成了更高浓度的优先形成乙酸作为难氧化组分的物种。
- 来自相同HTL原料的APs在反应行为上可能表现出显著差异，即使其理化性质看似相似。

在HTL与WO集成到绿色精炼概念的更广泛背景下，该模型有助于过程级理解，支持系统级评估（如热集成、生命周期评估和技术经济分析）。重要的是，模型高精度地捕捉了WO过程的热释放，使其可直接应用于能量平衡计算和热集成研究。其简洁性使其特别适用于需要计算效率的大型系统模型。未来工作将集中于更好地理解乙酸和铵的行为以扩展模型保真度，以及集成到精炼概念中。