能源系统优化中沼气厂的负荷依赖原料混合效率与消化器稳定性建模：基于ADM1的产率与VFA/TIC约束

《Energy Conversion and Management-X》：Biogas plants in energy system optimization: modelling loading-dependent feedstock mix efficiency and digester stability

　　

沼气厂在能源系统优化中：建模负荷依赖的原料混合效率与消化器稳定性。用于沼气能源系统优化的厌氧消化模型1号（ADM1，Anaerobic Digestion Model No. 1）导出产率与挥发性脂肪酸/总无机碳（VFA/TIC）约束。甲烷产率约束将有机负荷率（OLR）与原料组成相联系。VFA/TIC约束产生用于稳定性的内生最小粪污份额。与ADM1相比，线性模型将甲烷产率高估至多56 %。凸公式在99 %的运行中与ADM1甲烷产量的偏差保持在10 %以内。

研究背景与问题：厌氧消化通过水解、酸化、产乙酸和产甲烷等耦合微生物步骤将有机物转化为甲烷（CH₄）和二氧化碳（CO₂），各步骤速率和微生物群不同，挥发性脂肪酸（VFA）若不能被及时消耗，会消耗碱度缓冲并抑制产甲烷。消化器状态不仅由原料决定，还受水力停留时间（HRT）、有机负荷率（OLR）、pH及VFA/总无机碳（TIC）等运行量影响。ADM1及其系统降阶模型（如ADM1-R2）能较准确再现厌氧过程，但保持非线性、非凸数学性质。能源系统与技术经济优化模型依赖凸性以处理大决策空间，现有沼气建模常采用时间不变沼气潜力、批次实验产率、相关模型或Wei?bach潜力估算，将年度物料投入作为决策变量并以固定产率计算甲烷，或将底物混合与沼气调度分离；生物质供应链和选址模型也常按原料类型固定转化效率，或以容量成本项表示沼气厂。过程稳定性则常用OLR和HRT线性不等式保障，或以碳氮比和估算pH约束组成，但这些量不描述VFA积累与缓冲平衡，同一混合原料可在某一负荷下稳定、在另一负荷下酸化。实验证据显示OLR升高会显著降低比沼气产率，因此缺失的是可被凸优化模型承载的运行点依赖表示：甲烷产率需在运行范围内连续近似，稳定性需作为可行与不可行运行的边界。为此，Mathias Heiker、Matthias Kraume和Sandra Rosenberger开展研究，推导可纳入二次约束规划（QCP）的效率约束和VFA/TIC稳定性约束，二者均为凸约束并可形成二阶锥规划（SOCP），使调度问题可全局最优求解；系数对ADM1-R2稳态运行点网格校准，并在最小调度优化模型中与线性规划（LP）及ADM1-R2耦合非线性规划（NLP）比较。其意义在于把厌氧消化的运行点依赖以可扩展凸约束形式嵌入能源系统与投资模型。

主要关键技术与方法（不超过250字）：研究人员采用ADM1-R2降阶模型描述单级连续搅拌釜反应器（CSTR，38 °C），构建四维实验设计（DoE）稳态网格，覆盖OLR、玉米青贮可再生份额、牛粪和固体粪污份额，共13,629个唯一运行点；以最小二乘和坐标下降校准甲烷效率约束，以逻辑回归校准VFA/TIC稳定性约束，校准走廊为1≤OLR≤6 kg oDM/m³/d且VFA/TIC≤0.8；将约束作为QCP/SOCP嵌入热电联产（CHP）调度模型，与LP及ADM1-R2耦合NLP比较，并用ADM1-R2对LP和QCP所选日投料进行事后复算。未涉及人群样本队列。

研究结果：
ADM1网格表征：通过ADM1-R2稳态网格积分发现，纯底物甲烷产率效率随OLR升高近似线性下降，斜率按池分组而非按单一底物分组；混合原料效率接近纯底物加性预测，平均残差+0.007、均方根偏差0.0134，存在小幅协同；VFA/TIC受OLR、液体污泥/固体粪污分配和可再生份额驱动，液体污泥最敏感，固体粪污提供碱度，可再生底物残余酸少。玉米单消化在7.5 kg oDM/m³/d以上崩溃，模型可靠性限于6 kg oDM/m³/d校准上限。

校准参数与约束验证：通过最小二乘和坐标下降校准效率约束，留出验证的均方根偏差为0.0165—0.0278，平均偏差为负，纯底物过估最大+0.017，负荷上端最大+0.032；稳定性约束逻辑回归对校准走廊97.6 %的网格点分类正确，124个不稳定点中1个被允许，6,359个稳定点中157个被拒绝；纯池OLR限为可再生21.2、固体粪污5.7、液体污泥1.6 kg oDM/m³/d。

示范：对CHP调度的效应：在固定上网电价（FT）和EPEX现货（EPEX spot）下，LP几乎单选固体粪污，QCP与NLP随价格混合并缩放投料且达玉米上限；LP高估甲烷54—63 %（FT）和59—60 %（spot），QCP偏差为?1至?2 %；LP报告利润高于ADM1参考99—120千欧/年（FT）和69—98千欧/年（spot），QCP与参考差7—10千欧/年（FT）和小于2千欧/年（spot）。跨情景中LP在669/671个运行无例外高估，平均+55.6 %，无一在±10 %内；QCP在64 %运行中±5 %内、99 %运行中±10 %内，中位+4 %；三池可用时QCP偏差?1.8至?2.6 %，受限两池时+3.8至+7.3 %。

运行时间：LP中位22.6 s，QCP中位25.4 s（约增加12 %），NLP中位1,499 s（约66倍），QCP为22—40 s，NLP为144—31,380 s；凸公式与初值无关并达全局最优，NLP解时强烈依赖初值。

讨论部分总结：研究人员将贡献定位为更易集成的凸约束，而非仅更快数值求解；约束仅在校准范围内、仅对校准变量提供信息，不含中间产物或范围外运行信息，但补足了线性模型省略的运行点依赖。系数源于ADM1-R2模拟并继承其参数偏差，玉米单消化崩溃阈值7.75 kg oDM/m³/d与实验约3.7 kg oDM/m³/d不一致，故校准限于6。对碳水化合物一阶水解常数k_CH、乙酸营养型产甲烷菌最大生长速率μ_m,ac及其半饱和常数K_S,ac的敏感性分析显示，稳定性约束负荷系数约0.038—0.060 (kg oDM/m³/d)^-1，分类正确率97.0—99.3 %，效率截距约变1 %，低μ_m,ac下污泥主导进料无法校准。运行柔性增加时，日变化因数周停留时间被平均，稳态校准仍可沿用；底物选择转向残余物和废弃物时稳定性边界更关键，线性模型的产率误判会进入投资决策。局限包括五底物三池、单温度、单ADM1-R2参数集和单稳定性阈值；其他底物需新网格；稳定性约束中仅液体污泥限在校准走廊内被识别，固体粪污和可再生限为外推；阈值0.6—1.0变化保持池排序但绝对限变化；液体污泥效率截距保守使受限运行偏差达3.8—4.3 %。

研究结论部分翻译：研究人员推导了两个凸约束，将甲烷产率和过程稳定性对厌氧消化运行点的依赖纳入优化模型。效率约束通过底物特异截距和负荷依赖下降限制甲烷产量；稳定性约束通过VFA/TIC稳定性比，将允许负荷率限制为池组成的函数。二者均可表示为二阶锥，使调度优化问题可用内点法全局最优求解。ADM1-R2降阶网格支持效率约束的结构假设：纯底物效率随负荷升高近线性下降，斜率按池分组，混合物遵循潜在加性预测，平均残差+0.007，与其线性化误差同量级；因此可分离形式并非混合物行为的近似，其残差已被量化并偏向低估甲烷产率。在调度示范中，传统线性公式对每种底物组合平均高估甲烷产量55.6 %，所提公式中位偏差4 %；线性公式报告利润高于非线性参考99—120/69—98千欧/年，却实现最高实际产气成本，凸公式利润与其自身底物选择结果一致。凸公式增加的过程细节约使求解时间为线性模型的112 %，而ADM1耦合非线性模型约为66倍，因此适合嵌入调度和投资模型。后续方向包括以实测运行数据校准系数、为差异显著的底物类别建立新池，以及因约束对任意非负系数均保持二阶锥可表示而将其用于选址和原料物流模型。