Multiple emission constraints

Phil · 31 May 2021 15:47

Hi there again,

I am still simulating the German electricity and heat market in a simplified form and I’m struggling with a problem again. Due to the EU-ETS-System the power-plants with electricity outputs got emission constraints. In the next years the power-plants with a thermal output are regulated with a CO2 taxes. And later on with their own emission limit.
Additional there is a methanation transformer with flows to both markets in the system, which reduces the CO2-Emissions.

How can I implement two different emission constraints? (The CO2-taxes are not a problem)

currently i’m using:

constraints.emission_limit(model, limit = Emissionslimit) #unit: [t]

but i don’t know how to split the emissions in this case.

I tried to understand the documentation but I still don’t know exactly how pyomo is working and what is necessary to implement these two Constraints with the methanation in the background.

github.com

oemof/oemof-examples/blob/master/oemof_examples/oemof.solph/v0.3.x/flexible_modelling/add_constraints.py

# -*- coding: utf-8 -*-
"""
General description
-------------------
This script shows how to add an individual constraint to the oemof solph
OperationalModel.
The constraint we add forces a flow to be greater or equal a certain share
of all inflows of its target bus. Moreover we will set an emission constraint.
If you want to model an emission constraint and do not want to get into
the details of creating constraints, there is another example 

https://github.com/oemof/oemof-examples/blob/master/oemof_examples/
oemof.solph/v0.3.x/emission_constraint/emission_constraint.py

that shows how to use the built-in oemof function.

Installation requirements
-------------------------
This example requires the latest version of oemof. Install by:

This file has been truncated. show original

I hope you can help me with my problem.

Greetings,
Philipp

PyMoF · 4 June 2021 08:15

Hi Phil,
maybe you can implement just another keyword for the emissions of the flows of the second system.
keyword of flow of system 1:
inflow_args.update({‘emission_factor_1’: x})

keyword of flow of system 2:
inflow_args.update({‘emission_factor_2’: y})

parts of the constraint:
generic_integral_limit(om,
keyword=‘emission_factor_1’,
flows=flows,
limit=limit_1)
generic_integral_limit(om,
keyword=‘emission_factor_2’,
flows=flows,
limit=limit_2)

And then use the documented “generic_integral_limit” twice.

Kind regards,
Mathias