brightway - 从 LCIA 获取结果分布

Question

我想使用现有的影响评估方法之一，进行蒙特卡罗模拟，了解生态发明过程的影响分布。是否有示例笔记本或执行此操作的说明？

Answer 1

这是最简单的方法（对于随机活动和方法）：

from brightway2 import *
import numpy as np
ecoinvent = Database("ecoinvent 3.2 cutoff")

MonteCarlo 类派生自 LCA 类，并像 LCA 对象一样被实例化。
my_MC = MonteCarloLCA({ecoinvent.random():1}, methods.random())

假设您想获得 1000 个样本：
iterations = 1000

您可以创建一个空的 numpy 数组来收集结果：
scores = np.zeros([1, iterations])

然后，您使用对象计算分数next：

for iteration in range(1000):
    next(my_MC)
    scores[0, iteration] = my_MC.score

在此示例中，这将是一个包含 1000 个元素的 numpy 数组。然后，您可以使用您熟悉的任何统计模块来分析此数组。

还有其他几个基于 Monte Carlo 的类，允许其他功能。看看源代码，你可能会发现一些有用的东西。