在brightway2中创建非常简单的LCIA方法

时间:2017-01-04 14:28:25

标签: python brightway

我对bw2很新,我试图在不使用数据库的情况下制作简单的LCA,而只是使用手动定义的库存和LCIA方法。我使用了“生命周期评估的计算结构”一书的第11章中的示例中的值。

我能够创建库存并运行LCI计算:

t_db = Database("testdb")

t_db.write({
    ("testdb", "Electricity production"):{
        'name':'Electricity production',
        'unit': 'kWh', 
        'exchanges': [{
                'input': ('testdb', 'Fuel production'),
                'amount': 2,
                'unit': 'kg',
                'type': 'technosphere'
            },{
                'input': ('testdb', 'Carbon dioxide'),
                'amount': 1,
                'unit': 'kg',
                'type': 'biosphere'
            },{
                'input': ('testdb', 'Sulphur dioxide'),
                'amount': 0.1,
                'unit': 'kg',
                'type': 'biosphere'
            },{
                'input': ('testdb', 'Electricity production'), #important to write the same process name in output
                'amount': 10,
                'unit': 'kWh',
                'type': 'production'
            }]
        },
    ('testdb', 'Fuel production'):{
        'name': 'Fuel production',
        'unit': 'kg',
        'exchanges':[{
                'input': ('testdb', 'Carbon dioxide'),
                'amount': 10,
                'unit': 'kg',
                'type': 'biosphere'
            },{
                'input': ('testdb', 'Sulphur dioxide'),
                'amount': 2,
                'unit': 'kg',
                'type': 'biosphere'
            },{
                'input': ('testdb', 'Crude oil'),
                'amount': -50,
                'unit': 'kg',
                'type': 'biosphere'
            },{
                'input': ('testdb', 'Fuel production'),
                'amount': 100,
                'unit': 'kg',
                'type': 'production'
            }]
    },
    ('testdb', 'Carbon dioxide'):{'name': 'Carbon dioxide', 'unit':'kg', 'type': 'biosphere'},
    ('testdb', 'Sulphur dioxide'):{'name': 'Sulphur dioxide', 'unit':'kg', 'type': 'biosphere'},
    ('testdb', 'Crude oil'):{'name': 'Crude oil', 'unit':'kg', 'type': 'biosphere'}

    })

functional_unit = {t_db.get("Electricity production") : 1000}
lca = LCA(functional_unit) 
lca.lci()
print(lca.inventory)

然而,当我创建虚构的LCIA方法时,问题就开始了(为简单起见,所有CF都设置为1)。这是我使用的代码,但很明显它不起作用。关键问题似乎是未能将交易所从库存链接到LCIA方法。

myLCIAdata = [[('biosphere', 'Carbon dioxide'), 1.0], 
         [('biosphere', 'Sulphur dioxide'), 1.0],
         [('biosphere', 'Crude oil'), 1.0]]

method_key = ('simplemethod', 'imaginaryendpoint', 'imaginarymidpoint')
Method(method_key).validate(myLCIAdata) #returns "TRUE"
Method(method_key).register() 
Method(method_key).write(myLCIAdata)
Method(method_key).load() #check everything works
lca = LCA(functional_unit, method_key) #run LCA calculations again with method
lca.characterized_inventory

结果是<3x2 sparse matrix of type '<class 'numpy.float64'>' with 0 stored elements in Compressed Sparse Row format> 所以一个空矩阵。我知道我犯了什么错误吗?我是否应该像现有数据库一样为每个交易所获取唯一标识符?我已经检查了本网站上的bw2教程,文档和以前的问题,但找不到答案。提前谢谢。

1 个答案:

答案 0 :(得分:0)

你非常接近。在定义CF时,请执行以下操作:

myLCIAdata = [[('biosphere', 'Carbon dioxide'), 1.0], 
              [('biosphere', 'Sulphur dioxide'), 1.0],
              [('biosphere', 'Crude oil'), 1.0]]

取第一行,这意味着在数据库biosphere中,会有一个代码为Carbon dioxide的流。但是您没有biosphere数据库,您将所有内容放入名为testdb的数据库中:

    ('testdb', 'Carbon dioxide'):{'name': 'Carbon dioxide', 'unit':'kg', 'type': 'biosphere'},

因此,要么在特征化因子列表中使用正确的数据库名称,要么创建一个名为biosphere的单独生物圈数据库。

请注意,而不是:

method_key = ('simplemethod', 'imaginaryendpoint', 'imaginarymidpoint')
Method(method_key).validate(myLCIAdata) #returns "TRUE"
Method(method_key).register() 
Method(method_key).write(myLCIAdata)

你应该这样做:

method_key = ('simplemethod', 'imaginaryendpoint', 'imaginarymidpoint')
my_method = Method(method_key)
my_method.validate(myLCIAdata)
my_method.register() 
my_method.write(myLCIAdata)

(你的代码没有被破坏,但可能更优雅。)