streamline regret_runs; use only one co2 price in ST model

Author: lindnemi

scripts/pypsa-de/additional_functionality.py

@@ -58,17 +58,6 @@ def add_capacity_limits(
 
                 cname = f"capacity_{sense}-{ct}-{c.name}-{carrier.replace(' ', '-')}"
 
-                if snakemake.params.get("regret_run"):
-                    logger.info(
-                        f"Skipping capacity limit adjustment for {c.name} {carrier} with planning horizons {investment_year}, because of regret run."
-                    )
-                    if cname in n.global_constraints.index:
-                        logger.warning(
-                            f"Global constraint {cname} already exists. Dropping it."
-                        )
-                        n.global_constraints.drop(cname, inplace=True)
-                    continue
-
                 if cname in n.global_constraints.index:
                     logger.warning(
                         f"Global constraint {cname} already exists. Dropping and adding it again."
@@ -842,44 +831,46 @@ def additional_functionality(n, snapshots, snakemake):
     investment_year = int(snakemake.wildcards.planning_horizons[-4:])
     constraints = snakemake.params.solving["constraints"]
 
-    add_capacity_limits(
-        n, investment_year, constraints["limits_capacity_min"], snakemake, "minimum"
-    )
-
-    add_capacity_limits(
-        n, investment_year, constraints["limits_capacity_max"], snakemake, "maximum"
-    )
+    if not snakemake.config.get("regret_run"):
+        add_capacity_limits(
+            n, investment_year, constraints["limits_capacity_min"], snakemake, "minimum"
+        )
 
-    add_power_limits(n, investment_year, constraints["limits_power_max"])
+        add_capacity_limits(
+            n, investment_year, constraints["limits_capacity_max"], snakemake, "maximum"
+        )
 
-    if snakemake.wildcards.clusters != "1":
-        h2_import_limits(n, investment_year, constraints["limits_volume_max"])
+        add_power_limits(n, investment_year, constraints["limits_power_max"])
 
-        electricity_import_limits(n, investment_year, constraints["limits_volume_max"])
+        if snakemake.wildcards.clusters != "1":
+            h2_import_limits(n, investment_year, constraints["limits_volume_max"])
 
-    if investment_year >= 2025:
-        h2_production_limits(
-            n,
-            investment_year,
-            constraints["limits_volume_min"],
-            constraints["limits_volume_max"],
-        )
+            electricity_import_limits(
+                n, investment_year, constraints["limits_volume_max"]
+            )
 
-    add_h2_derivate_limit(n, investment_year, constraints["limits_volume_max"])
+        if investment_year >= 2025:
+            h2_production_limits(
+                n,
+                investment_year,
+                constraints["limits_volume_min"],
+                constraints["limits_volume_max"],
+            )
+        add_h2_derivate_limit(n, investment_year, constraints["limits_volume_max"])
+
+        if isinstance(constraints["co2_budget_national"], dict):
+            add_national_co2_budgets(
+                n,
+                snakemake,
+                constraints["co2_budget_national"],
+                investment_year,
+            )
+        else:
+            logger.warning("No national CO2 budget specified!")
 
     # force_boiler_profiles_existing_per_load(n)
     force_boiler_profiles_existing_per_boiler(n)
 
-    if isinstance(constraints["co2_budget_national"], dict):
-        add_national_co2_budgets(
-            n,
-            snakemake,
-            constraints["co2_budget_national"],
-            investment_year,
-        )
-    else:
-        logger.warning("No national CO2 budget specified!")
-
     if isinstance(constraints.get("decentral_heat_pump_budgets"), dict):
         add_decentral_heat_pump_budgets(
             n,

scripts/pypsa-de/solve_regret.py

@@ -177,10 +177,56 @@ def fix_capacities(realization, decision, scope="DE", strict=False):
     )
     np.random.seed(solve_opts.get("seed", 123))
 
+    strict = False
+
     n = fix_capacities(
-        realization, decision, scope=snakemake.params.scope_to_fix, strict=False
+        realization, decision, scope=snakemake.params.scope_to_fix, strict=strict
+    )
+    if strict:
+        snakemake.params.solving["options"]["load_shedding"] = True
+
+    n.add("Carrier", "H2 vent", color="#dd2e23", nice_name="H2 vent")
+
+    n.add(
+        "Generator",
+        n.buses.query("carrier=='H2'").index,
+        " H2 vent",
+        bus=n.buses.query("carrier=='H2'").index,
+        carrier="H2 vent",
+        sign=-1e-3,  # Adjust sign to measure p and p_nom in kW instead of MW
+        marginal_cost=100,  # Eur/kWh
+        p_nom=1e9,  # kW
     )
 
+    snakemake.config["regret_run"] = True
+
+    if snakemake.config["regret_run"]:
+        constraint_names = [  # TODO assert really everything gets dropped
+            "H2_production_limit_upper-DE",
+            "H2_production_limit_lower-DE",
+            "capacity_minimum-DE-Generator-solar",
+            "capacity_maximum-DE-Generator-onwind",
+            "capacity_maximum-DE-Generator-offwind",
+            "capacity_maximum-DE-Generator-solar",
+            "capacity_maximum-DE-Store-co2-sequestered",
+            "capacity_maximum-DE-Link-HVC-to-air",
+            "H2_import_limit-DE",
+            "H2_export_ban-DE",
+            "Electricity_import_limit-DE",
+            "renewable_oil_import_limit-DE",
+            "methanol_import_limit-DE",
+            "renewable_gas_import_limit-DE",
+            "H2_derivate_import_limit-DE",
+            "H2_derivate_meoh_gas_import_limit-DE",
+            "H2_derivate_oil_meoh_import_limit-DE",
+            "H2_derivate_oil_meoh_gas_import_limit-DE",
+            "H2_derivate_oil_gas_import_limit-DE",
+        ]
+        logger.info("Regret run detected. Dropping the following constraints:")
+        logger.info(constraint_names)
+
+        n.global_constraints.drop(constraint_names, errors="ignore", inplace=True)
+
     if solve_opts["post_discretization"].get("enable") and not solve_opts.get(
         "skip_iterations"
     ):
@@ -221,9 +267,13 @@ def fix_capacities(realization, decision, scope="DE", strict=False):
             p_max_pu=0,
             p_nom_extendable=True,
             carrier="co2",
-            marginal_cost=realization.global_constraints.loc["CO2Limit", "mu"],
+            marginal_cost=(
+                realization.global_constraints.loc["CO2Limit", "mu"]
+                + realization.global_constraints.loc["co2_limit-DE", "mu"]
+            ),
         )
         n.global_constraints.drop("CO2Limit", inplace=True)
+        n.global_constraints.drop("co2_limit-DE", inplace=True)
 
     with memory_logger(
         filename=getattr(snakemake.log, "memory", None), interval=logging_frequency