Brightway2与Activity Browser的LCA结果差异排查求助

问题背景 #

• 运行环境：搭载M1芯片的MacBook Pro
• 使用工具：Brightway2、Activity Browser
• 基础数据库：Ecoinvent-3.8(apos)自定义参数化数据库

核心差异现象 #

• 计算结果：Brightway2输出值为82439873.99596766，Activity Browser为65659359.41451239，后者与Simapro结果更匹配
• 差异范围：仅在**非长期化的CML v4.8 2016 abiotic depletion potential (ADP)**方法下差异显著，GWP100等其他LCIA方法结果基本一致
• 已确认前提：二者使用完全相同的数据库和LCIA方法参数
• 警告信息：Brightway2运行时触发UmfpackWarning: (almost) singular matrix! (estimated cond. number: 2.52e+13)

用到的Brightway2代码 #

lca=bw.LCA({bw.get_activity(('ei38', 'fleet, vehicle cycle')):1},('CML v4.8 2016 no LT',
  'material resources: metals/minerals no LT',
  'abiotic depletion potential (ADP): elements (ultimate reserves) no LT')
)
lca.lci()

lca.lcia()
lca.score

可能的原因排查方向 #

1. 矩阵奇异性的直接影响：
   Umfpack警告说明LCI矩阵接近奇异，数值求解稳定性极差。ADP方法对资源消耗系数敏感度高，微小的计算误差会被极端放大，直接导致结果偏差。Activity Browser可能采用了不同的线性求解器或数值稳定策略，规避了该问题。

2. LCIA方法的实现细节差异：
   即使参数一致，二者对CML v4.8 no LT方法的底层处理可能存在分歧：

   • 非长期化因子的应用逻辑不同
   • 金属/矿物资源的归类、权重计算方式有细微差别
   • 对低贡献但系数极端的小众流处理规则不一致

3. 自定义数据库的参数同步问题：
   虽然确认数据库一致，但参数化过程中，Brightway2可能在导入或计算时对部分参数（如默认值、单位转换）的处理与Activity Browser不同，尤其针对Ecoinvent apos系统模型中的特定活动。

4. 求解器配置差异：
   Brightway2默认使用Umfpack求解器，Activity Browser可能切换了更稳定的求解器（如MUMPS）。可尝试在Brightway2中手动更换求解器验证：

   from brightway2.solvers import choose_solver
   choose_solver("mumps")  # 需提前安装MUMPS相关依赖
   

内容的提问来源于stack exchange，提问作者Lgs