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

Merge branch 'master' of https://github.com/LukasFrankenQ/pypsa-eur

Browse Source
This commit is contained in:
LukasFrankenQ 2023-11-14 01:37:00 +00:00
commit c1f44ea86f
2 changed files with 11 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
scripts/prepare_sector_network.py
View File

@ -3473,12 +3473,13 @@ def add_enhanced_geothermal(
n.links.at[bus + " geothermal organic rankine cycle", "efficiency"] = efficiency_orc
if snakemake.params.sector["enhanced_geothermal_flexible"]:
# this StorageUnit represents flexible operation using the geothermal reservoir.
# Hence, it is counter-intuitive to install it at the surface bus,
# this is however the more lean and computationally efficient solution.
max_hours = snakemake.params.sector["enhanced_geothermal_reservoir_max_hours"]
max_hours = snakemake.params.sector[
"enhanced_geothermal_reservoir_max_hours"
]
boost = snakemake.params.sector["enhanced_geothermal_reservoir_max_boost"]
max_hours = max_hours * boost
@ -3488,7 +3489,7 @@ def add_enhanced_geothermal(
bus=f"{bus} geothermal heat surface",
carrier="geothermal heat",
p_nom_extendable=True,
p_min_pu=-1. - boost,
p_min_pu=-1.0 - boost,
max_hours=max_hours,
)

14
scripts/solve_network.py
View File

@ -779,15 +779,13 @@ def add_geothermal_chp_constraint(n):
def add_flexible_egs_constraint(n):
well_index = n.links.loc[n.links.carrier == 'geothermal heat'].index
storage_index = n.storage_units.loc[n.storage_units.carrier == 'geothermal heat '].index
well_index = n.links.loc[n.links.carrier == "geothermal heat"].index
storage_index = n.storage_units.loc[
n.storage_units.carrier == "geothermal heat "
].index
p_nom_rhs = (
n.model["Link-p_nom"].loc[well_index]
)
p_nom_lhs = (
n.model["StorageUnit-p_nom"].loc[storage_index]
)
p_nom_rhs = n.model["Link-p_nom"].loc[well_index]
p_nom_lhs = n.model["StorageUnit-p_nom"].loc[storage_index]
n.model.add_constraints(
p_nom_lhs <= p_nom_rhs,