brightway - LCIA からの結果の分布の取得

Question

既存の影響評価方法の 1 つを使用して、モンテカルロ シミュレーションを行い、ecoinvent プロセスの影響の分布を希望します。これを行うためのノートブックまたは手順の例はありますか?

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 配列になります。次に、使い慣れた統計モジュールを使用して、この配列を分析できます。

他にもいくつかのモンテカルロ ベースのクラスがあり、他の機能を使用できます。ソースコードを見てください。役に立つものが見つかるかもしれません。