废弃咖啡渣协同厌氧消化制备可再生天然气的技术经济分析与生命周期评价：欧洲与北美生物炼制厂比较

《ACS Sustainable Chemistry & Engineering》：Techno-Economic and Life Cycle Assessment of Renewable Natural Gas Production from Spent Coffee Grounds Co-Digestion: A European and North American Biorefinery Comparison

　　

废弃咖啡渣（spent coffee grounds，SCGs）是一类在全球范围内大量产生的有机残余物，其挥发性固体（volatile solids，VS）超过总固体（total solids，TS）的90%，生物甲烷潜力为250–350 mL CH₄/g VS。本研究将厌氧消化模型1号（Anaerobic Digestion Model No. 1，ADM1）动力学整合到BioSTEAM过程模拟框架中，对两种总进料量均为5000 kg/h的协同消化情景下由SCGs生产可再生天然气（renewable natural gas，RNG）进行耦合技术经济分析（techno-economic analysis，TEA）与生命周期评价（life cycle assessment，LCA）：（i）马德里情景，即SCGs与来自西班牙Mercamadrid的果蔬废弃物（fruit and vegetable waste，FVW）协同消化；（ii）爱荷华情景，即SCGs与美国爱荷华州混合草原生物质协同消化。尽管这两个情景针对特定地点进行了参数化设定，但也可更广泛地视为SCG分别与食物废弃物生物质和草本生物质协同消化的可迁移示例。沼气通过单乙醇胺（monoethanolamine，MEA）吸收法提纯为RNG，并回收水热碳化（hydrothermal carbonization，HTC）生物炭作为联产品。马德里情景实现了186.3 kg/h的RNG产率、$41.8/GJ的最低售价（minimum selling price，MSP）、8.0%的内部收益率（internal rate of return，IRR）以及2.3 M$的净现值（net present value，NPV）。爱荷华情景的RNG产率为180.6 kg/h，MSP为$55.9/GJ，IRR为3.3%，NPV为–7.6 M$；蒙特卡罗分析（Monte Carlo analysis，n = 2000）表明，D_VS是两个情景中主导RNG产率不确定性的因素，而SCG价格对MSP影响最强。在TRACI 2.1框架下并采用系统扩展法的LCA结果显示，马德里与爱荷华情景的全球变暖潜势（global warming potential，GWP）分别为每千克RNG对应–10.96和–12.31 kg CO₂-eq，10个影响类别中分别有9个和8个呈现净负值。结果表明，基于SCG的协同消化生物炼制厂能够带来显著环境效益，但其经济性具有情景依赖性，且并不一定与最优气候效益相一致。

该文发表于《ACS Sustainable Chemistry 》。论文围绕废弃咖啡渣（spent coffee grounds，SCGs）资源化利用展开，聚焦于通过协同厌氧消化（anaerobic codigestion）与后续沼气提纯构建可再生天然气（renewable natural gas，RNG）生物炼制体系，并同步评价其经济可行性与环境可持续性。研究背景在于，全球咖啡消费持续增长，使SCGs成为数量庞大但尚未充分利用的有机固体废弃物。SCGs挥发性固体（volatile solids，VS）含量高，具备较高生物甲烷潜力，理论上适宜进入厌氧消化（anaerobic digestion，AD）体系；但其木质纤维素结构、脂质与酚类抑制成分又使单独消化面临降解受限、甲烷产率不稳定与运行风险升高等问题。现有研究多集中于AD动力学、单一工艺经济性或孤立环境评价，对于将机理模型、技术经济分析（techno-economic analysis，TEA）和生命周期评价（life cycle assessment，LCA）耦合起来，进而在不同地理与原料背景下开展可比性分析的文献仍较少。因此，开展本研究的必要性在于：一方面，构建可迁移的SCG协同消化评价框架；另一方面，澄清不同共底物类型对RNG产量、经济性及环境收益的综合影响。

研究人员据此构建了两个具有代表性的生物炼制情景。其一为马德里情景，以SCGs与Mercamadrid果蔬废弃物（fruit and vegetable waste，FVW）按质量比50:50协同消化；其二为爱荷华情景，初始筛选为草原生物质与SCGs 50:50，后经灵敏度分析修正为40:60，以提高SCG占比并改善经济表现。两种情景均设定总进料量5000 kg/h，采用中温连续搅拌釜式反应器（continuous stirred-tank reactor，CSTR）在36 °C、30 d水力停留时间（hydraulic retention time，HRT）下运行，产出经MEA吸收提纯后的RNG，同时将消化残渣进行水热碳化（hydrothermal carbonization，HTC）以制备生物炭（biochar），并回收液态CO₂。研究最终得出：马德里情景在经济上更具可行性，而爱荷华情景在气候变化缓解方面表现更优；然而最佳环境绩效并不自动对应最佳经济绩效，这体现了RNG生物炼制系统中显著的TEA-LCA权衡关系。

在方法层面，研究人员主要采用以下关键技术。首先，利用ADM1动力学在经验证的Python环境中模拟厌氧消化过程，并与BioSTEAM生物炼制仿真平台集成，实现全流程物料与能量衡算及TEA计算。其次，采用TRACI 2.1方法开展LCA，并使用系统扩展法处理生物炭与回收CO₂等联产品替代效应。再次，基于拉丁超立方抽样的蒙特卡罗分析（n = 2000）以及确定性敏感性分析识别RNG产率、MSP与GWP的关键驱动因素。样本与供给数据来源于西班牙马德里自治区、Mercamadrid以及美国爱荷华州的官方统计和行业资料，用于估算SCG、FVW和草原生物质的区域可获得性。

研究结果部分可按原文小标题概括如下。

Preliminary SCG Composition Sensitivity Analysis
研究人员首先对SCG质量分数进行初步敏感性分析，以确定更合理的进料配比。结果显示，在爱荷华情景中，随着SCG占比从10%提升至90%，MSP由$109.8/GJ显著下降至$36.5/GJ，NPV由–$43.7 M升至$5.9 M，表现出明显的单调经济改善趋势。基于这一结果，研究将爱荷华基准情景由原始50:50修订为40:60（草原生物质:SCG）。马德里情景中，提高SCG比例同样可降低MSP并提高NPV，但原始50:50方案已低于$45/GJ市场参考线且NPV为正，因此继续保留50:50作为基准。该部分结果说明，共底物比例是影响系统经济性的关键设计参数。

ADM1 Simulation Validation: Specific Methane Yield
在正式报告TEA与LCA结果前，研究人员以比甲烷产率（specific methane yield，SMY）验证模型可信度。马德里情景模拟得到249 mL CH₄/g VS，虽低于部分文献中SCG与食物废弃物共消化的批量BMP结果，但仍处于相关咖啡残渣共消化文献范围内。爱荷华50:50筛选情景模拟SMY为146 mL CH₄/g VS，略低于已有草原生物质相关研究报道值。研究据此认为，该ADM1-BioSTEAM实现能够较好反映SCG协同消化体系的量级特征，同时指出原料脱脂状态、批式与连续式运行差异会影响直接对比。该部分为后续系统评价提供了模型基础。

Mass and Energy Balances
物料与能量衡算显示，马德里情景产沼气695.0 kg/h，甲烷摩尔分数50.1%，RNG产率186.3 kg/h，生物炭产率797.6 kg/h；爱荷华情景产沼气659.1 kg/h，甲烷摩尔分数50.8%，RNG产率180.6 kg/h，生物炭产率1165.9 kg/h。爱荷华情景生物炭产量更高，主要源于草原生物质消化残渣总固体含量和木质素含量较高。两种方案的电力与热需求均相对有限，说明情景差异更多来自消化残渣分配和产物结构，而非提纯能耗的大幅差异。

Techno-Economic Analysis
技术经济分析表明，马德里情景的MSP为$41.8/GJ，IRR为8.0%，NPV为2.3 M$，简单回收期（simple payback period，SPP）为12.0年，显示其在基准条件下具有经济吸引力。其优势核心在于FVW被设定为零成本原料，甚至在现实情境中还有可能带来处置费减免或入厂费收益。相比之下，爱荷华情景MSP为$55.9/GJ，IRR仅3.3%，NPV为–7.6 M$，SPP为16.8年，贴现回收期超过20年寿命期，说明其在当前假设下经济性较弱。草原生物质采购成本为其主要负担之一。研究还指出，生物炭是两个情景的重要收入来源，特别是爱荷华情景因生物炭产量更高而具备更大的联产品收入潜力。

Monte Carlo Uncertainty Analysis
通过蒙特卡罗不确定性传播，研究人员发现两个情景的经济稳健性差异明显。马德里情景MSP均值为$46.4/GJ，略高于基准值，NPV均值为–$1.0 M，表明虽具潜力但仍受不确定性影响。爱荷华情景MSP均值为$60.5/GJ，NPV均值为–$10.8 M，整体更为脆弱。标准化回归系数（standardized regression coefficient，SRC）显示，D_VS是RNG产率最主要的不确定性来源；而SCG价格是MSP方差的最强驱动因子，在爱荷华情景中草原生物质价格也具有重要影响。该结果说明，提升原料品质表征与降低采购价格波动，是提高系统稳健性的两条不同路径。

Sensitivity Analysis
确定性敏感性分析进一步表明，HRT与温度变化对RNG产量的影响仅为0.1%–0.6%，远小于原料组成和价格因素。这意味着，在既定中温连续消化条件下，运行参数微调对总体绩效改善有限，而原料配比和供应链经济性才是决定性因素。

Life Cycle Assessment
生命周期评价方面，采用TRACI 2.1与系统扩展法后，两种情景均表现出显著净负全球变暖潜势。马德里情景GWP为–10.96 kg CO₂-eq/kg RNG，爱荷华情景为–12.31 kg CO₂-eq/kg RNG。负排放贡献主要来自SCG避免填埋的原料信用以及生物炭替代传统木炭的系统扩展信用；正环境负担则主要包括甲烷泄漏、电力消耗以及爱荷华情景中特有的草原生物质种植过程。全TRACI 2.1结果显示，马德里情景在10个中点影响类别中有9个为净负值，仅臭氧消耗潜势（ozone depletion potential，ODP）略为正值；爱荷华情景则有8个类别为净负值，ODP与富营养化（eutrophication，EUTR）为正，其中EUTR负担主要与草原种植及其上游供应链有关。分配方法敏感性分析显示，不论采用系统扩展、质量分配还是经济分配，两情景GWP均保持净负值，从而增强了主要环境结论的稳健性。

Uncertainty Analysis: Standardized Regression Coefficients
环境不确定性分析表明，D_VS不仅主导RNG产率，也主导GWP变化。由于基准GWP为负值，D_VS增加虽然提升RNG产量，却会使GWP向较少负值方向移动，即气候优势略有减弱。与此相对，MSP主要受SCG价格及当地共底物价格控制。研究由此指出，经济驱动因素与环境驱动因素并不一致，RNG系统的不确定性控制需针对具体输出指标实施差异化管理。

Uncertainty Analysis: TEA-LCA
TEA-LCA耦合分析所展示的帕累托前沿（Pareto front）进一步证明，较低MSP与更负GWP之间存在下降型权衡关系。也就是说，若追求更强气候收益，通常需接受更高的成本；若追求最低成本，则往往会牺牲部分环境优势。研究认为，决策上更适合关注“膝点”区域，即在成本温和增加时可换取较明显环境收益提升的操作区间。

Allocation-Method Comparison
分配方法比较结果显示，相较于系统扩展法，质量分配与经济分配会使RNG承担更多系统环境负担，因此GWP绝对值减小、表现得不那么负。但无论何种方法，两种情景的GWP均为负值，说明其气候效益结论具有方法学稳健性。

Cross-Scenario Comparison and Contextualization
跨情景比较是本文的核心贡献之一。研究表明，马德里情景经济性优于爱荷华情景，原因主要在于零成本FVW共底物与更宽松的SCG供给边界；爱荷华情景虽在GWP上更有优势，但其SCG回收要求接近实际可获得性的上限，同时还伴随富营养化负担。作者特别强调，这两个案例不仅是地理上具体的西欧与北美案例，也可被视为“SCG+食物废弃物”和“SCG+草本生物质”两种可转移的协同消化原型，因此具有更广泛的方法学推广意义。

讨论部分集中体现出三点。其一，SCG协同消化确实能够在循环生物经济（circular bioeconomy）框架下兼顾有机废弃物管理、低碳燃料生产与联产品增值，但经济性高度依赖于共底物成本与区域供应链结构。其二，较优的环境绩效并不必然转化为较优的投资回报，说明未来工艺优化不能仅围绕甲烷产率展开，还需综合考虑原料组织、联产品市场和政策支持。其三，爱荷华情景暴露出草本生物质路线在供应连续性、回收效率和营养盐环境负担上的现实约束，而马德里情景则证明城市食物废弃物流更适合作为SCG协同消化的支撑底物。

研究结论部分可译为：本研究在BioSTEAM过程模拟框架内，构建并验证了两种基于SCG协同消化的RNG生物炼制路线的耦合TEA与LCA模型。由于FVW原料成本为零，马德里情景在基准条件下保持经济吸引力（MSP = $41.8/GJ；NPV = $2.3 M）；而爱荷华情景未达到7%的经济目标（MSP = $55.9/GJ；NPV = –$7.6 M），仍需更低的原料成本、入厂费收入或政策支持。虽然研究以马德里与爱荷华为背景，但更应将其理解为SCG分别与食物废弃物生物质和草本生物质协同消化的可迁移示例，而非仅限于特定场址设计。两种情景均实现了显著净负的全球变暖潜势，分别为–10.96和–12.31 kg CO₂-eq/kg RNG；在TRACI 2.1的10项影响类别中，分别有9项和8项呈现净负值，且三种分配方法下GWP结论一致。D_VS始终是两种情景中RNG产率不确定性的主导因素，而SCG价格主导MSP波动。原料供应链方面，马德里情景具备较可信的区域供给基础，而爱荷华情景仅在SCG高效回收的严格假设下才表现可行。HTC生物炭在两个系统中均是重要联产品，尤其在爱荷华情景中更值得进一步实验表征。总体而言，研究支持将SCG协同消化视为欧洲与北美背景下兼顾有机废弃物治理、可再生燃料生产和气候变化缓解的可行循环生物经济路径。