Constrain every efuel trade link individually (#42)

* more specific import bans instead of export bans

* remove duplicate code

* .png -> .pdf

* [GW] in title of Stromnetzausbau

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

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

* temporarily fix pypsa version, because of inconsistency with DC projects

---------

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>

Author: Micha

Snakefile

@@ -583,9 +583,9 @@ rule plot_ariadne_variables:
         transmission_investment_csv=RESULTS + "ariadne/transmission_investment.csv",
         trassenlaenge_csv=RESULTS + "ariadne/trassenlaenge.csv",
         Kernnetz_Investment_plot=RESULTS + "ariadne/Kernnetz_Investment_plot.png",
-        elec_trade=RESULTS + "ariadne/elec-trade-DE.png",
-        h2_trade=RESULTS + "ariadne/h2-trade-DE.png",
-        trade_balance=RESULTS + "ariadne/trade-balance-DE.png",
+        elec_trade=RESULTS + "ariadne/elec-trade-DE.pdf",
+        h2_trade=RESULTS + "ariadne/h2-trade-DE.pdf",
+        trade_balance=RESULTS + "ariadne/trade-balance-DE.pdf",
     log:
         RESULTS + "logs/plot_ariadne_variables.log",
     script:
@@ -687,11 +687,11 @@ rule plot_ariadne_report:
         ),
     output:
         elec_price_duration_curve=RESULTS
-        + "ariadne/report/elec_price_duration_curve.png",
-        elec_price_duration_hist=RESULTS + "ariadne/report/elec_price_duration_hist.png",
-        backup_capacity=RESULTS + "ariadne/report/backup_capacity.png",
-        backup_generation=RESULTS + "ariadne/report/backup_generation.png",
-        elec_prices_spatial_de=RESULTS + "ariadne/report/elec_prices_spatial_de.png",
+        + "ariadne/report/elec_price_duration_curve.pdf",
+        elec_price_duration_hist=RESULTS + "ariadne/report/elec_price_duration_hist.pdf",
+        backup_capacity=RESULTS + "ariadne/report/backup_capacity.pdf",
+        backup_generation=RESULTS + "ariadne/report/backup_generation.pdf",
+        elec_prices_spatial_de=RESULTS + "ariadne/report/elec_prices_spatial_de.pdf",
         results=directory(RESULTS + "ariadne/report"),
         elec_transmission=directory(RESULTS + "ariadne/report/elec_transmission"),
         h2_transmission=directory(RESULTS + "ariadne/report/h2_transmission"),
@@ -710,6 +710,6 @@ rule plot_ariadne_report:
 rule ariadne_report_only:
     input:
         expand(
-            RESULTS + "ariadne/report/elec_price_duration_curve.png",
+            RESULTS + "ariadne/report/elec_price_duration_curve.pdf",
             run=config_provider("run", "name"),
         ),

config/config.yaml

@@ -432,7 +432,7 @@ solving:
         H2 pipeline retrofitted: 0.05
       fractional_last_unit_size: true
   constraints:
-    efuel_export_ban: true
+    efuel_export_ban: false
     limits_capacity_max:
       Generator:
         onwind:

envs/environment.yaml

@@ -17,7 +17,7 @@ dependencies:
 - jpype1
 
 # Inhouse packages
-- pypsa>=0.32
+- pypsa<=0.31
 - atlite>=0.2.9
 - linopy
 - powerplantmatching>=0.5.15,<0.6

scripts/pypsa-de/additional_functionality.py

@@ -640,33 +640,42 @@ def add_h2_derivate_limit(n, investment_year, limits_volume_max):
             ]
         ].index
 
-        incoming_p = (
-            n.model["Link-p"].loc[:, incoming] * n.snapshot_weightings.generators
-        ).sum()
-        outgoing_p = (
-            n.model["Link-p"].loc[:, outgoing] * n.snapshot_weightings.generators
-        ).sum()
+        carrier_idx_dict = {
+            "renewable_oil": 0,
+            "methanol": 1,
+            "renewable_gas": 2,
+            "H2_derivate": [0, 1, 2],
+        }
+        for carrier, idx in carrier_idx_dict.items():
+            cname = f"{carrier}_import_limit-{ct}"
+
+            incoming_p = (
+                n.model["Link-p"].loc[:, incoming[idx]]
+                * n.snapshot_weightings.generators
+            ).sum()
+            outgoing_p = (
+                n.model["Link-p"].loc[:, outgoing[idx]]
+                * n.snapshot_weightings.generators
+            ).sum()
 
-        lhs = incoming_p - outgoing_p
+            lhs = incoming_p - outgoing_p
 
-        cname = f"H2_derivate_import_limit-{ct}"
+            n.model.add_constraints(lhs <= limit, name=f"GlobalConstraint-{cname}")
 
-        n.model.add_constraints(lhs <= limit, name=f"GlobalConstraint-{cname}")
+            if cname in n.global_constraints.index:
+                logger.warning(
+                    f"Global constraint {cname} already exists. Dropping and adding it again."
+                )
+                n.global_constraints.drop(cname, inplace=True)
 
-        if cname in n.global_constraints.index:
-            logger.warning(
-                f"Global constraint {cname} already exists. Dropping and adding it again."
+            n.add(
+                "GlobalConstraint",
+                cname,
+                constant=limit,
+                sense="<=",
+                type="",
+                carrier_attribute="",
             )
-            n.global_constraints.drop(cname, inplace=True)
-
-        n.add(
-            "GlobalConstraint",
-            cname,
-            constant=limit,
-            sense="<=",
-            type="",
-            carrier_attribute="",
-        )
 
     # Export bans on efuels are implemented in modify_prenetwork by restricting p_max_pu of the DE -> EU links