Erhalt der Verteilung der Ergebnisse von LCIA

Question

Ich möchte die Verteilung für die Auswirkungen von Ecoinvent-Prozessen, unter Verwendung einer der bestehenden Folgenabschätzung Methoden, Monte-Carlo-Simulationen tun. Gibt es ein Beispiel Notebook oder Anweisungen dazu?

Answer 1

Hier ist der einfachste Weg, es zu tun (für eine zufällige Aktivität und Verfahren):

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

Die Monte-Carlo-Klasse ergibt sich aus der LCA-Klasse und instanziiert wird wie ein LCA-Objekt.
my_MC = MonteCarloLCA({ecoinvent.random():1}, methods.random())

Sagen Sie 1000 Proben erhalten möchten:
iterations = 1000

Sie können ein leeres numpy Array erstellen, die Ergebnisse zu sammeln:
scores = np.zeros([1, iterations])

Dann berechnen Sie Noten mit next auf Objekt:

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

In diesem Beispiel wird dies ein numpliges Array mit 1000 Elementen sein. Sie können dieses Array dann mit den statistischen Modulen analysieren, mit denen Sie sich auskennen.

Es gibt mehrere andere Monte-Carlo-basierte Klassen, die andere Funktionalitäten erlauben. Schauen Sie sich die source code an, vielleicht finden Sie etwas Nützliches.