厌氧共消化中拮抗行为建模：一种基于相互作用的方法应用于树篱修剪残渣与果蔬垃圾

《Bioresource Technology Reports》：Modeling antagonistic behavior in anaerobic co-digestion: An interaction-based approach applied to hedge pruning residues and fruit and vegetable waste

　　

厌氧共消化（Anaerobic co-digestion, AcoD）可通过混合具有互补特性的底物来提高甲烷产量。然而，使用传统动力学模型（将每种底物混合物独立处理）往往难以解释协同或拮抗行为的相互作用。该研究评估了中温条件下树篱修剪残渣（Hedge pruning residues, HPR）（包括 Duranta erecta、Teucrium fruticans 和 Myrtus communis）与果蔬垃圾（Fruit and vegetable waste, FVW）的生化甲烷潜力（Biochemical methane potential, BMP）。果蔬垃圾表现出高生物降解性（453 NmL CH₄ gVS^-1），生物降解指数（Biodegradability index, BDI）达97.7%，而树篱修剪残渣因其木质纤维素顽固性导致甲烷产量较低。Myrtus 对厌氧消化（Anaerobic digestion, AD）表现出抑制作用。相反，Duranta 和 Teucrium 混合物表现出拮抗行为，甲烷产量和共消化性能指数（Co-digestion performance index, CPI）随树篱修剪残渣比例增加而降低。实验数据采用传统动力学模型和基于两分数一级动力学（Two-fraction first-order, TFFO）推导的现象学相互作用模型进行分析。传统模型对甲烷产量曲线拟合极佳，但需对每种混合物独立进行参数估计。而基于相互作用的公式通过单一组底物特异性参数和相互作用项同时描述了所有单底物和共消化实验。基于相互作用模型获得的拟合优度指标与传统公式相当（全局 R² 高达0.998）。独立于参数估计的共消化实验被保留用于预测评估。基于相互作用的模型和两分数一级动力学模型均产生了同样精确的预测。结果表明，拮抗共消化行为可通过可转移的底物级参数化表示，同时保持与传统模型相当的预测性能。

论文解读文章

厌氧消化（Anaerobic digestion, AD）作为一种可持续处理有机废弃物并产生可再生能源（沼气）的技术，在循环经济框架中日益重要。然而，单一底物消化常因营养失衡或抑制物积累而效率受限。厌氧共消化（Anaerobic co-digestion, AcoD）通过混合互补底物可提升甲烷产量，但传统动力学模型将每种混合物视为独立系统，难以解析底物间协同或拮抗相互作用的机制。树篱修剪残渣（Hedge pruning residues, HPR）是城市绿色废弃物的重要组分，其木质纤维素结构及潜在生物活性次级代谢物（如多酚、萜类）可能抑制微生物活动；而果蔬垃圾（Fruit and vegetable waste, FVW）则具有高生物降解性。目前针对两者共消化的系统研究较少，且缺乏能够显式捕捉拮抗效应的简化建模工具。为此，研究人员开展本研究，旨在通过生化甲烷潜力（Biochemical methane potential, BMP）实验评估HPR与FVW的共消化行为，并提出一种基于相互作用的动力学模型，以统一参数化描述底物间的非加性效应。该研究发表于《Bioresource Technology Reports》。

在关键技术方法方面，研究人员从某批发市场收集了FVW（每日约10吨有机废弃物），从某城市园林维护点采集了三种HPR（Duranta erecta、Teucrium fruticans、Myrtus communis），所用接种污泥来自某市政污水处理厂的中温消化器（35±1°C，水力停留时间约20天）。BMP实验按照VDI 4630标准在500 mL玻璃瓶中进行，接种物与底物比（以挥发性固体计）为4:1，实验温度35°C，持续27天。甲烷产量通过排水法监测，气体成分由气相色谱分析。理论甲烷潜力采用修正的Buswell-Mueller方程基于元素组成计算。性能指标包括生物降解指数（Biodegradability index, BDI）和共消化性能指数（Co-digestion performance index, CPI）。动力学建模采用五种传统模型（一级、两分数一级、Monod型、修正Gompertz、叠加模型）以及基于两分数一级动力学衍生的相互作用模型。参数估计通过非线性回归（RMSE最小化）并采用Akaike信息准则（Akaike Information Criterion, AIC）比较模型性能，预测评估使用独立保留的60:40混合比例实验。

研究结果如下：

3.1 接种物与底物表征
通过总固体（Total solids, TS）、挥发性固体（Volatile solids, VS）及元素分析得出，HPR的有机含量（VS/TS>0.88）和C/N比（15.4-27.1）均处于厌氧消化适宜范围内，但FVW的C/N比为20.5，更接近最优值（20-30）。除Duranta外，各底物C/N比接近最优范围，表明共消化中可能的协同效应并非源于营养平衡改善。

3.2 实验BMP
FVW单底物消化显示高甲烷潜力（453 NmL CH₄ gVS^-1），生物降解指数达97.7%，无滞后期。三种HPR中，Duranta和Teucrium的甲烷产量较低（分别为167和131 NmL CH₄ gVS^-1），BDI分别为35.6%和28.9%，归因于木质纤维素顽固性。Myrtus则表现出抑制效应，实测BMP为负值（-125 NmL gVS^-1），表明其完全抑制了产甲烷过程。对于Duranta和Teucrium与FVW的共消化，所有混合比例下CPI均小于1（0.77-0.94），证实拮抗行为，且拮抗程度在低HPR比例时更显著。

3.3 BMP建模
传统动力学模型对甲烷累积曲线拟合极佳（R²>0.99），但需独立估计每种混合物的参数。其中，固定理论BMP的两分数一级模型（TFFO*）在AIC评价中表现最优，但其动力学常数（k_f、k_s）随混合比变化缺乏系统性趋势（Spearman等级相关系数不显著），限制了可推广性。相比，提出的相互作用模型通过一组底物特异性参数（9个）和三个相互作用参数（k_int、n₁、n₂）同时描述了5条实验曲线，全局R²达0.996-0.998，AIC在Duranta系统中略高于TFFO*（484 vs 406），但Teucrium系统中则更低（385 vs 395）。参数分析显示，FVW的固有动力学常数在两种共消化系统中高度一致，表明相互作用项有效地捕获了拮抗效应。相互作用指数揭示不同拮抗模式：Duranta系统中n₁=0.76、n₂=0.63，表明拮抗依赖于两种底物浓度；Teucrium系统中n₂=0，表明拮抗主要受FVW调控，不随HPR浓度变化。预测验证（60:40混合比）表明，相互作用模型与TFFO*具有几乎相同的预测精度（R² 0.960-0.997，RMSE 3-13 NmL gVS^-1）。

讨论部分强调，相互作用模型在保持预测能力的同时，显著减少了参数数量（从15个降至9个），且参数可跨混合物转移，克服了传统模型需逐混合物重调参数的局限性。该框架为解释共消化系统中的非加性行为提供了一种简洁且物理一致的方法，有助于未来BMP研究中底物相互作用的理解和预测。

研究结论翻译如下：
果蔬垃圾与树篱修剪残渣的厌氧共消化主要表现为拮抗行为，反映为甲烷产量降低和共消化性能指数（CPI）低于1。在所研究的三种树篱物种中，Myrtus communis表现出强抑制效应，导致负的生化甲烷潜力（BMP）值，无法实现有意义的共消化。传统动力学模型成功复现了累积甲烷产量曲线，其中固定理论BMP的两分数一级（TFFO*）模型具有优异的描述性能。然而，独立拟合的动力学参数在不同混合物组成间一致性有限，使底物相互作用的解释变得困难。相比之下，所提出的相互作用模型通过单一组底物特异性动力学参数和三个相互作用参数成功描述了所有单底物和共消化实验，同时保持了相当的精度。底物级动力学常数与独立回归结果一致，表明相互作用项捕获了拮抗效应，而未显著改变各底物的固有降解特征。拟合的相互作用指数揭示了Duranta和Teucrium系统中不同的拮抗模式，体现了现象学框架在不预设底物间函数关系的情况下再现不同拮抗行为的能力。使用独立共消化实验的预测评估进一步表明，相互作用模型与TFFO*公式实现了同样准确的预测。总体而言，相互作用模型在预测精度上并未超越最佳传统公式，但通过大幅减少参数集实现了相当的性能。TFFO*需15个独立拟合参数描述每种共消化系统的5条实验曲线，而相互作用模型采用统一的框架，仅用9个底物特异性和相互作用参数。模型复杂度的降低未以牺牲拟合优度或预测能力为代价。因此，所提出的公式为描述厌氧共消化系统中的非加性行为提供了一种物理一致且简约的方法，并可能促进未来BMP研究中底物相互作用的解释和预测。