Merge pull request PyPSA#35 from PyPSA/landfill

allow landfill of plastics, and option for waste-to-energy

Author: Micha

config/config.yaml

@@ -4,7 +4,7 @@
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#run
 run:
-  prefix: 20240524mastrchps
+  prefix: 20240527wasteCHP
   name:
   # - CurrentPolicies
   - KN2045_Bal_v4
@@ -259,6 +259,16 @@ industry:
     2040: 0.12
     2045: 0.16
     2050: 0.20
+  HVC_environment_sequestration_fraction:
+    2020: 0.1
+    2025: 0.1
+    2030: 0.12
+    2035: 0.15
+    2040: 0.18
+    2045: 0.20
+    2050: 0.20
+  waste_to_energy: true
+  waste_to_energy_cc: true
 
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#solving

workflow/scripts/export_ariadne_variables.py

@@ -215,6 +215,8 @@ def _get_capacities(n, region, cap_func, cap_string="Capacity|"):
     var[cap_string + "Electricity|Biomass"] = \
         var[cap_string + "Electricity|Biomass|Solids"]
 
+    var[cap_string + "Electricity|Non-Renewable Waste"] = \
+        capacities_electricity.filter(like="waste CHP").sum()
 
     var[cap_string + "Electricity|Coal|Hard Coal"] = \
         capacities_electricity.filter(like="coal").sum()                                             
@@ -289,9 +291,6 @@ def _get_capacities(n, region, cap_func, cap_string="Capacity|"):
     var[cap_string + "Electricity|Hydrogen"] = \
         var[cap_string + "Electricity|Hydrogen|FC"]
 
-    # var[cap_string + "Electricity|Non-Renewable Waste"] = 
-    # ! Not implemented
-
     var[cap_string + "Electricity|Nuclear"] = \
         capacities_electricity.get("nuclear", 0)
 
@@ -409,6 +408,7 @@ def _get_capacities(n, region, cap_func, cap_string="Capacity|"):
             cap_string + "Electricity|Hydro",
             cap_string + "Electricity|Hydrogen",
             cap_string + "Electricity|Nuclear",
+            cap_string + "Electricity|Non-Renewable Waste",
             ]].sum()
 
     # Test if we forgot something
@@ -460,6 +460,9 @@ def _get_capacities(n, region, cap_func, cap_string="Capacity|"):
         var[cap_string + "Heat|Biomass"],
         capacities_central_heat.filter(like="biomass").sum()
     )
+
+    var[cap_string + "Heat|Non-Renewable Waste"] = \
+        capacities_central_heat.filter(like="waste CHP").sum()
 
     var[cap_string + "Heat|Resistive heater"] = \
         capacities_central_heat.filter(like="resistive heater").sum()
@@ -507,7 +510,8 @@ def _get_capacities(n, region, cap_func, cap_string="Capacity|"):
         var[cap_string + "Heat|Gas"] +
         var[cap_string + "Heat|Processes"] +
         #var[cap_string + "Heat|Hydrogen"] +
-        var[cap_string + "Heat|Heat pump"]
+        var[cap_string + "Heat|Heat pump"] +
+        var[cap_string + "Heat|Non-Renewable Waste"]
     )
 
     # This check requires further changes to n.statistics
@@ -783,6 +787,20 @@ def get_primary_energy(n, region):
     # )
     # ! There are CC sub-categories that could be used
 
+
+    waste_CHP_E_usage, waste_CHP_H_usage = get_CHP_E_and_H_usage(
+        n, "non-sequestered HVC", region)
+
+    var["Primary Energy|Waste|Electricity"] = \
+        waste_CHP_E_usage
+    var["Primary Energy|Waste|Heat"] = \
+        waste_CHP_H_usage
+    var["Primary Energy|Waste"] = (
+        var["Primary Energy|Waste|Electricity"]
+        + var["Primary Energy|Waste|Heat"]
+    )
+
+
     coal_usage = n.statistics.withdrawal(
         bus_carrier=["lignite", "coal"], 
         **kwargs,
@@ -1013,6 +1031,9 @@ def get_secondary_energy(n, region):
         electricity_supply.get("H2 Fuel Cell", 0)
     # ! Add H2 Turbines if they get implemented
 
+    var["Secondary Energy|Electricity|Waste"] = \
+        electricity_supply.filter(like="waste CHP").sum()
+
     var["Secondary Energy|Electricity|Curtailment"] = \
         n.statistics.curtailment(
             bus_carrier=["AC", "low voltage"], **kwargs
@@ -1065,6 +1086,7 @@ def get_secondary_energy(n, region):
         #+ var["Secondary Energy|Electricity|Transmission Losses"]
         #+ var["Secondary Energy|Electricity|Storage Losses"]
         + var["Secondary Energy|Electricity|Hydrogen"]
+        + var["Secondary Energy|Electricity|Waste"]
     )
 
     assert isclose(
@@ -1116,6 +1138,8 @@ def get_secondary_energy(n, region):
         var["Secondary Energy|Heat|Electricity|Heat Pumps"] 
         + var["Secondary Energy|Heat|Electricity|Resistive"] 
     )
+    var["Secondary Energy|Heat|Waste"] = \
+        heat_supply.filter(like="waste CHP").sum()
     var["Secondary Energy|Heat|Other"] = \
         heat_supply.reindex(
             [
@@ -1127,6 +1151,7 @@ def get_secondary_energy(n, region):
             ]
         ).sum()
     # TODO remember to specify in comments
+
     var["Secondary Energy|Heat"] = (
         var["Secondary Energy|Heat|Gas"]
         + var["Secondary Energy|Heat|Biomass"]
@@ -1135,6 +1160,7 @@ def get_secondary_energy(n, region):
         + var["Secondary Energy|Heat|Electricity"]
         + var["Secondary Energy|Heat|Other"]
         + var["Secondary Energy|Heat|Coal"]
+        + var["Secondary Energy|Heat|Waste"]
     )
     assert isclose(
         var["Secondary Energy|Heat"],
@@ -1563,6 +1589,16 @@ def get_final_energy(n, region, _industry_demand, _energy_totals):
     # var["Final Energy|Geothermal"] = \
     # ! Not implemented
 
+    waste_withdrawal = n.statistics.withdrawal(
+        bus_carrier=["non-sequestered HVC"], 
+        **kwargs,
+    ).filter(
+        like=region,
+    ).groupby("carrier").sum().multiply(MWh2PJ)
+
+    var["Final Energy|Waste"] = \
+        waste_withdrawal.filter(like="waste CHP").sum()
+
     var["Final Energy incl Non-Energy Use incl Bunkers"] = (
         var["Final Energy|Industry"]
         + var["Final Energy|Residential and Commercial"]
@@ -1600,6 +1636,11 @@ def get_emissions(n, region, _energy_totals):
         bus_carrier="co2",**kwargs
     ).filter(like=region).filter(like="CHP").multiply(t2Mt)
 
+    # exclude waste CHPs because they are accounted separately
+    CHP_emissions = CHP_emissions[
+        ~CHP_emissions.index.get_level_values(
+            "carrier").str.contains("waste")]
+
     CHP_E_to_H =  (
         n.links.loc[CHP_emissions.index.get_level_values("name")].efficiency 
         / n.links.loc[CHP_emissions.index.get_level_values("name")].efficiency2
@@ -1794,7 +1835,11 @@ def get_emissions(n, region, _energy_totals):
         ).sum()
 
     var["Emissions|CO2|Supply|Non-Renewable Waste"] = \
-        co2_emissions.get("naphtha for industry")
+        co2_emissions.reindex([
+            "HVC to air",
+            "waste CHP",
+            "waste CHP CC",
+        ]).sum()
 
     # var["Emissions|CO2|Energy|Supply|Liquids"] = \
     # Our only Liquid production is Fischer-Tropsch