Biogas

config/config.yaml

@@ -4,7 +4,7 @@
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#run
 run:
-  prefix: 20240829mergemaster
+  prefix: 20240904_biogas_baseload
   name:
   # - CurrentPolicies
   - KN2045_Bal_v4
@@ -62,7 +62,7 @@ scenario:
   ll:
   - vopt
   clusters:
-  - 22 #current options: 22, 44
+  - 23 #current options: 22, 44
   opts:
   - ''
   sector_opts:
@@ -552,7 +552,7 @@ co2_price_add_on_fossils:
   2020: 25
   2025: 60
 
-must_run_biomass:
+must_run_biogas:
   enable: false
-  p_min_pu: 0.7
+  p_min_pu: 0.6
   regions: ['DE']

workflow/Snakefile

@@ -237,7 +237,7 @@ rule modify_prenetwork:
             "sector", "H2_retrofit_capacity_per_CH4"
         ),
         transmission_costs=config_provider("costs", "transmission"),
-        biomass_must_run=config_provider("must_run_biomass"),
+        biogas_must_run=config_provider("must_run_biogas"),
         clustering=config_provider("clustering", "temporal", "resolution_sector"),
         H2_plants=config_provider("electricity", "H2_plants_DE"),
         land_transport_electric_share=config_provider(

workflow/scripts/export_ariadne_variables.py

@@ -481,31 +481,26 @@ def _get_capacities(n, region, cap_func, cap_string="Capacity|", costs=None):
     #     capacities_electricity = \
     #         capacities_electricity.multiply(MW2GW)
 
-    capacities_biomass = capacities_electricity.filter(like="biomass")
+    var[cap_string + "Electricity|Biomass|w/ CCS"] = capacities_electricity.get(
+        "urban central solid biomass CHP CC ", 0
+    )
 
-    var[cap_string + "Electricity|Biomass|w/ CCS"] = capacities_biomass.filter(
-        like="CC"
+    var[cap_string + "Electricity|Biomass|w/o CCS"] = capacities_electricity.reindex(
+        ["urban central solid biomass CHP", "solid biomass", "biogas"]
     ).sum()
 
-    var[cap_string + "Electricity|Biomass|w/o CCS"] = capacities_biomass.reindex(
-        ["urban central solid biomass CHP", "solid biomass"]
+    var[cap_string + "Electricity|Biomass|Solids"] = capacities_electricity.filter(
+        like="solid biomass"
     ).sum()
 
-    var[cap_string + "Electricity|Biomass|Solids"] = var[
-        [
-            cap_string + "Electricity|Biomass|w/ CCS",
-            cap_string + "Electricity|Biomass|w/o CCS",
-        ]
-    ].sum()
-
-    # Ariadne does no checks, so we implement our own?
-    assert isclose(
-        var[cap_string + "Electricity|Biomass|Solids"], capacities_biomass.sum()
+    var[cap_string + "Electricity|Biomass|Gases and Liquids"] = (
+        capacities_electricity.get("biogas", 0)
     )
 
-    var[cap_string + "Electricity|Biomass"] = var[
-        cap_string + "Electricity|Biomass|Solids"
-    ]
+    var[cap_string + "Electricity|Biomass"] = (
+        var[cap_string + "Electricity|Biomass|Solids"]
+        + var[cap_string + "Electricity|Biomass|Gases and Liquids"]
+    )
 
     var[cap_string + "Electricity|Non-Renewable Waste"] = capacities_electricity.filter(
         like="waste CHP"
@@ -1208,8 +1203,8 @@ def get_primary_energy(n, region):
         ~biomass_usage.index.str.contains("CC")
     ].sum()
 
-    var["Primary Energy|Biomass|Electricity"] = biomass_CHP_E_usage + biomass_usage.get(
-        "solid biomass", 0
+    var["Primary Energy|Biomass|Electricity"] = (
+        biomass_CHP_E_usage + biomass_usage.reindex(["solid biomass", "biogas"]).sum()
     )
     var["Primary Energy|Biomass|Heat"] = biomass_CHP_H_usage + biomass_usage.get(
         "urban central solid biomass boiler", 0
@@ -1355,12 +1350,18 @@ def get_secondary_energy(n, region, _industry_demand):
         + var["Secondary Energy|Electricity|Coal"]
     )
 
-    var["Secondary Energy|Electricity|Biomass|w/o CCS"] = electricity_supply.get(
-        "urban central solid biomass CHP", 0
-    ) + electricity_supply.get("solid biomass", 0)
+    var["Secondary Energy|Electricity|Biomass|w/o CCS"] = electricity_supply.reindex(
+        ["urban central solid biomass CHP", "solid biomass", "biogas"]
+    ).sum()
     var["Secondary Energy|Electricity|Biomass|w/ CCS"] = electricity_supply.get(
         "urban central solid biomass CHP CC", 0
     )
+    var["Secondary Energy|Electricity|Biomass|Solid"] = electricity_supply.filter(
+        like="solid biomass"
+    ).sum()
+    var["Secondary Energy|Electricity|Biomass|Gaseous and Liquid"] = (
+        electricity_supply.get("biogas")
+    )
     var["Secondary Energy|Electricity|Biomass"] = (
         var["Secondary Energy|Electricity|Biomass|w/o CCS"]
         + var["Secondary Energy|Electricity|Biomass|w/ CCS"]

workflow/scripts/modify_prenetwork.py

@@ -584,16 +584,15 @@ def transmission_costs_from_modified_cost_data(
     n.links.loc[dc_b, "overnight_cost"] = overnight_cost
 
 
-def must_run_biomass(n, p_min_pu, regions):
+def must_run_biogas(n, p_min_pu, regions):
     """
-    Set p_min_pu for biomass generators to the specified value.
+    Set p_min_pu for biogas generators to the specified value.
     """
     logger.info(
-        f"Must-run condition enabled: Setting p_min_pu = {p_min_pu} for biomass generators."
+        f"Must-run condition enabled: Setting p_min_pu = {p_min_pu} for biogas generators."
     )
     links_i = n.links[
-        (n.links.carrier == "solid biomass")
-        & (n.links.bus0.str.startswith(tuple(regions)))
+        (n.links.carrier == "biogas") & (n.links.bus0.str.startswith(tuple(regions)))
     ].index
     n.links.loc[links_i, "p_min_pu"] = p_min_pu
 
@@ -857,11 +856,11 @@ def force_retrofit(n, params):
         snakemake.params.length_factor,
     )
 
-    if snakemake.params.biomass_must_run["enable"]:
-        must_run_biomass(
+    if snakemake.params.biogas_must_run["enable"]:
+        must_run_biogas(
             n,
-            snakemake.params.biomass_must_run["p_min_pu"],
-            snakemake.params.biomass_must_run["regions"],
+            snakemake.params.biogas_must_run["p_min_pu"],
+            snakemake.params.biogas_must_run["regions"],
         )
 
     if snakemake.params.H2_plants["enable"]: