Endogenize steel and cement production and add existing industry sites

config/config.default.yaml

@@ -973,6 +973,14 @@ sector:
       methanol: 121
       gas: 122
       oil: 125
+  endogenous_sectors:
+  - steel
+  - cement
+  hbi_relocation: false
+  steel_bof:
+    pledge: true
+    pledge_delay: 0
+    default_phase_out: 2035
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#industry
 industry:

config/config.pr1719.yaml

@@ -0,0 +1,24 @@
+%YAML 1.1
+---
+# yaml-language-server: $schema=./schema.json
+
+run:
+  name: pr1719
+
+foresight: myopic
+
+scenario:
+  clusters:
+    - 39
+  planning_horizons:
+    - 2030
+    - 2040
+
+sector:
+  methanol:
+    regional_methanol_demand: true
+  ammonia: regional
+
+clustering:
+  temporal:
+    resolution_sector: 365H

config/plotting.default.yaml

@@ -712,3 +712,19 @@ plotting:
     import NH3: '#e2ed74'
     import oil: '#93eda2'
     import methanol: '#87d0e6'
+    steel: '#22333B'
+    steel emission: '#5E503F'
+    steel emission CC: '#A9927D'
+    hbi: '#CCDBDC'
+    DRI: '#8E7C93'
+    EAF: '#586357'
+    H2 DRI: '#618ab0'
+    gas DRI: '#baaf68'
+    BOF: '#0A0908'
+    grey methanol: '#7C6C77'
+    cement: '#c98349'
+    cement emission: '#119DA4'
+    cement emission CC: '#0C7489'
+    cement kiln: '#D90368'
+    cement finishing: '#820263'
+    clinker: '#F5FDC6'

config/schema.json

@@ -4758,6 +4758,49 @@
               "type": "object"
             }
           }
+        },
+        "endogenous_sectors": {
+          "description": "Sectors modelled endogenously.",
+          "items": {
+            "enum": [
+              "steel",
+              "cement"
+            ],
+            "type": "string"
+          },
+          "type": "array"
+        },
+        "hbi_relocation": {
+          "default": false,
+          "description": "Allow the relocation of iron ore reduction step.",
+          "type": "boolean"
+        },
+        "steel_bof": {
+          "description": "Configuration for `sector.steel_bof` settings.",
+          "properties": {
+            "pledge": {
+              "default": true,
+              "description": "Decommissions coal blast furnaces at their pledged phase out year.",
+              "type": "boolean"
+            },
+            "pledge_delay": {
+              "default": 0,
+              "description": "Delays the decommission year by number of years.",
+              "type": "integer"
+            },
+            "default_phase_out": {
+              "anyOf": [
+                {
+                  "type": "boolean"
+                },
+                {
+                  "type": "integer"
+                }
+              ],
+              "default": 2035,
+              "description": "Decommissions the coal blast furnaces that do not have a pledged phase out year."
+            }
+          }
         }
       }
     },
@@ -7433,6 +7476,33 @@
         }
       }
     },
+    "_SteelBOFConfig": {
+      "description": "Configuration for `sector.steel_bof` settings.",
+      "properties": {
+        "pledge": {
+          "default": true,
+          "description": "Decommissions coal blast furnaces at their pledged phase out year.",
+          "type": "boolean"
+        },
+        "pledge_delay": {
+          "default": 0,
+          "description": "Delays the decommission year by number of years.",
+          "type": "integer"
+        },
+        "default_phase_out": {
+          "anyOf": [
+            {
+              "type": "boolean"
+            },
+            {
+              "type": "integer"
+            }
+          ],
+          "default": 2035,
+          "description": "Decommissions the coal blast furnaces that do not have a pledged phase out year."
+        }
+      }
+    },
     "_TechnologyMappingConfig": {
       "description": "Configuration for `electricity.estimate_renewable_capacities.technology_mapping` settings.",
       "properties": {
@@ -10868,6 +10938,49 @@
               "type": "object"
             }
           }
+        },
+        "endogenous_sectors": {
+          "description": "Sectors modelled endogenously.",
+          "items": {
+            "enum": [
+              "steel",
+              "cement"
+            ],
+            "type": "string"
+          },
+          "type": "array"
+        },
+        "hbi_relocation": {
+          "default": false,
+          "description": "Allow the relocation of iron ore reduction step.",
+          "type": "boolean"
+        },
+        "steel_bof": {
+          "description": "Configuration for `sector.steel_bof` settings.",
+          "properties": {
+            "pledge": {
+              "default": true,
+              "description": "Decommissions coal blast furnaces at their pledged phase out year.",
+              "type": "boolean"
+            },
+            "pledge_delay": {
+              "default": 0,
+              "description": "Delays the decommission year by number of years.",
+              "type": "integer"
+            },
+            "default_phase_out": {
+              "anyOf": [
+                {
+                  "type": "boolean"
+                },
+                {
+                  "type": "integer"
+                }
+              ],
+              "default": 2035,
+              "description": "Decommissions the coal blast furnaces that do not have a pledged phase out year."
+            }
+          }
         }
       }
     },

rules/build_sector.smk

@@ -1444,6 +1444,27 @@ def input_heat_source_power(w):
     }
 
 
+rule build_industry_plants:
+    params:
+        countries=config_provider("countries"),
+    input:
+        regions_onshore=resources("regions_onshore_base_s_{clusters}.geojson"),
+        ammonia="data/ammonia_plants.csv",
+        isi_database="data/1-s2.0-S0196890424010586-mmc2.xlsx",
+        gem_gcct="data/gem/Global-Cement-and-Concrete-Tracker_July-2025.xlsx",
+    output:
+        industry_plants=resources("industry_plants_{clusters}.csv"),
+    threads: 1
+    resources:
+        mem_mb=2000,
+    log:
+        logs("build_industry_plants_{clusters}.log"),
+    benchmark:
+        benchmarks("build_industry_plants_{clusters}")
+    script:
+        "../scripts/build_industry_plants.py"
+
+
 rule prepare_sector_network:
     params:
         time_resolution=config_provider("clustering", "temporal", "resolution_sector"),
@@ -1603,6 +1624,9 @@ rule prepare_sector_network:
             if config_provider("sector", "district_heating", "ates", "enable")(w)
             else []
         ),
+        industry_sector_ratios=resources(
+            "industry_sector_ratios_{planning_horizons}.csv"
+        ),
     output:
         resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}.nc"

rules/solve_myopic.smk

@@ -12,6 +12,8 @@ rule add_existing_baseyear:
         costs=config_provider("costs"),
         heat_pump_sources=config_provider("sector", "heat_pump_sources"),
         energy_totals_year=config_provider("energy", "energy_totals_year"),
+        MWh_NH3_per_tNH3=config_provider("industry", "MWh_NH3_per_tNH3"),
+        MWh_MeOH_per_tMeOH=config_provider("industry", "MWh_MeOH_per_tMeOH"),
     input:
         network=resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}.nc"
@@ -25,6 +27,7 @@ rule add_existing_baseyear:
             "existing_heating_distribution_base_s_{clusters}_{planning_horizons}.csv"
         ),
         heating_efficiencies=resources("heating_efficiencies.csv"),
+        industry_plants=resources("industry_plants_{clusters}.csv"),
     output:
         resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}_brownfield.nc"

rules/solve_perfect.smk

@@ -10,6 +10,8 @@ rule add_existing_baseyear:
         costs=config_provider("costs"),
         heat_pump_sources=config_provider("sector", "heat_pump_sources"),
         energy_totals_year=config_provider("energy", "energy_totals_year"),
+        countries=config_provider("countries"),
+        MWh_NH3_per_tNH3=config_provider("industry", "MWh_NH3_per_tNH3"),
     input:
         network=resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}.nc"
@@ -27,6 +29,9 @@ rule add_existing_baseyear:
         ),
         existing_heating="data/existing_infrastructure/existing_heating_raw.csv",
         heating_efficiencies=resources("heating_efficiencies.csv"),
+        regions_onshore=resources("regions_onshore_base_s_{clusters}.geojson"),
+        ammonia="data/ammonia_plants.csv",
+        isi_database="data/1-s2.0-S0196890424010586-mmc2.xlsx",
     output:
         resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}_brownfield.nc"