bjoern/pypsa-eur
1
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases Wiki Activity

solve_network: fix net vres share

Browse Source
This commit is contained in:
Fabian 2022-11-09 16:50:05 +01:00
parent 1e80edd265
commit 992efcaa91
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
scripts/solve_network.py
View File

@ -212,20 +212,24 @@ def add_SAFE_constraints(n, config):
def add_minRenew_constraints(n, config, o): def add_minRenew_constraints(n, config, o):
''' '''
Adds the constraint to have a minimum share of renewable energy production. Adds the constraint to have a minimum share of renewable energy production.
As renewable carriers the renewables from the configs listed under renewable are taken. As renewable carriers the renewables from the configs listed under renewable are taken.
To use this constraint simply add the wildcard RE{share} in the opts wildcard like RE0.8 for a 80% renewable share To use this constraint simply add the wildcard RE{share} in the opts wildcard like RE0.8 for a 80% renewable share
''' '''
import operator import operator
renewables=list(config["electricity"]["renewable_aim"].keys()) renewables=list(config["electricity"]["renewable_aim"].keys())
if len(o)>2: if len(o)>2:
share=float(o[2:]) share=float(o[2:])
renewables_i = n.generators[n.generators.carrier.str.contains("|".join(renewables))].index renewables_b = n.generators.carrier.str.contains("|".join(renewables))
renewables_i = n.generators[renewables_b].index
conventionals_i = n.generators[~renewables_b].index
weightings = n.snapshot_weightings.generators weightings = n.snapshot_weightings.generators
coeff = pd.DataFrame({c: weightings for c in renewables_i}) coeff = pd.DataFrame({c: weightings for c in renewables_i})
vars = get_var(n, "Generator", "p")[renewables_i] vres = get_var(n, "Generator", "p")[renewables_i]
lhs = linexpr((coeff, vars)).sum().sum() conv = get_var(n, "Generator", "p")[conventionals_i]
rhs = share * weightings @ n.loads_t.p_set.sum(1) lhs = linexpr(((1 - share) * coeff, vres)).sum().sum()
lhs += linexpr((- share * coeff, conv)).sum().sum()
rhs = 0
define_constraints(n, lhs, '>=', rhs, 'Carrier', 'min_generation_renewables') define_constraints(n, lhs, '>=', rhs, 'Carrier', 'min_generation_renewables')
else: else:
for tech in renewables: for tech in renewables: