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

Reduce retrofit potential in myopic optimization (#1258)

Browse Source 
* reduce potential of gas network that can be retrofitted to H2 when gas grid is not resolved in myopic optimization

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* adding comment

---------

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
This commit is contained in:
Toni Seibold 2024-09-25 13:23:50 +02:00 committed by GitHub
parent c5ebb66f72
commit 79cea53146
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 29 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
scripts/add_brownfield.py
View File

@ -87,30 +87,48 @@ def add_brownfield(n, n_p, year):
n.import_series_from_dataframe(c.pnl[tattr], c.name, tattr) n.import_series_from_dataframe(c.pnl[tattr], c.name, tattr)
# deal with gas network # deal with gas network
pipe_carrier = ["gas pipeline"]
if snakemake.params.H2_retrofit: if snakemake.params.H2_retrofit:
# subtract the already retrofitted from today's gas grid capacity # subtract the already retrofitted from the maximum capacity
h2_retrofitted_fixed_i = n.links[ h2_retrofitted_fixed_i = n.links[
(n.links.carrier == "H2 pipeline retrofitted") (n.links.carrier == "H2 pipeline retrofitted")
& (n.links.build_year != year) & (n.links.build_year != year)
].index ].index
gas_pipes_i = n.links[n.links.carrier.isin(pipe_carrier)].index h2_retrofitted = n.links[
CH4_per_H2 = 1 / snakemake.params.H2_retrofit_capacity_per_CH4 (n.links.carrier == "H2 pipeline retrofitted")
fr = "H2 pipeline retrofitted" & (n.links.build_year == year)
to = "gas pipeline" ].index
# today's pipe capacity
pipe_capacity = n.links.loc[gas_pipes_i, "p_nom"] # pipe capacity always set in prepare_sector_network to todays gas grid capacity * H2_per_CH4
# and is therefore constant up to this point
pipe_capacity = n.links.loc[h2_retrofitted, "p_nom_max"]
# already retrofitted capacity from gas -> H2 # already retrofitted capacity from gas -> H2
already_retrofitted = ( already_retrofitted = (
n.links.loc[h2_retrofitted_fixed_i, "p_nom"] n.links.loc[h2_retrofitted_fixed_i, "p_nom"]
.rename(lambda x: x.split("-2")[0].replace(fr, to) + f"-{year}") .rename(lambda x: x.split("-2")[0] + f"-{year}")
.groupby(level=0) .groupby(level=0)
.sum() .sum()
) )
remaining_capacity = pipe_capacity - CH4_per_H2 * already_retrofitted.reindex( remaining_capacity = pipe_capacity - already_retrofitted.reindex(
index=pipe_capacity.index index=pipe_capacity.index
).fillna(0) ).fillna(0)
n.links.loc[gas_pipes_i, "p_nom"] = remaining_capacity n.links.loc[h2_retrofitted, "p_nom_max"] = remaining_capacity
# reduce gas network capacity
gas_pipes_i = n.links[n.links.carrier == "gas pipeline"].index
if not gas_pipes_i.empty:
# subtract the already retrofitted from today's gas grid capacity
pipe_capacity = n.links.loc[gas_pipes_i, "p_nom"]
fr = "H2 pipeline retrofitted"
to = "gas pipeline"
CH4_per_H2 = 1 / snakemake.params.H2_retrofit_capacity_per_CH4
already_retrofitted.index = already_retrofitted.index.str.replace(fr, to)
remaining_capacity = (
pipe_capacity
- CH4_per_H2
* already_retrofitted.reindex(index=pipe_capacity.index).fillna(0)
)
n.links.loc[gas_pipes_i, "p_nom"] = remaining_capacity
n.links.loc[gas_pipes_i, "p_nom_max"] = remaining_capacity
def disable_grid_expansion_if_limit_hit(n): def disable_grid_expansion_if_limit_hit(n):