Merge pull request #781 from openego/features/adjust_PtH2_model

Implement method für adjusting PtH2-Model with multiple_link

etrago/appl.py

@@ -696,7 +696,7 @@ def run_etrago(args, json_path):
 
     # spatial clustering
     etrago.spatial_clustering()
-    etrago.spatial_clustering_gas()
+    etrago.spatial_clustering_gas()    
 
     # snapshot clustering
     etrago.snapshot_clustering()

etrago/network.py

@@ -110,6 +110,7 @@
     set_random_noise,
     set_trafo_costs,
     update_busmap,
+    adjust_PtH2_model,
 )
 
 logger = logging.getLogger(__name__)
@@ -372,6 +373,8 @@ def __init__(
 
     sclopf = iterate_sclopf
 
+    adjust_PtH2_model = adjust_PtH2_model
+
     def dc_lines(self):
         return self.filter_links_by_carrier("DC", like=False)
 

etrago/tools/utilities.py

@@ -3352,3 +3352,127 @@ def export_to_shapefile(pypsa_network, shape_files_path=None, srid=4326):
             v.loc[:, save_cols[k]].to_file(shp_filename)
 
     return components_dict
+
+
+def adjust_PtH2_model(self):
+    """
+    Adjust the modelling of electrolyzer with waste-heat and
+    O2-utilisation. The method creates a multiple-link-model
+    out of the single links from the dataset (power_to_H2-,
+    PtH2_waste_heat-, PtH2_O2-links) for each location where
+    coupling-product usage is possible. The resulting model
+    consists of a multiple link, additional buses and stores for
+    waste_heat and O2-utilisation, and the connection link to the
+    final Heat- and O2-Bus.
+
+    Parameters
+    ----------
+    etrago : :class:`Etrago
+        Overall container of Etrago
+
+    Returns
+    -------
+    None.
+
+    """
+    PtH2_links = self.network.links[
+        self.network.links.carrier == "power_to_H2"
+    ].index
+    PtH2_AC_buses = self.network.buses.loc[
+        self.network.links.loc[PtH2_links, "bus0"]
+    ].index.unique()
+
+    from itertools import count
+
+    store_ids = self.network.stores.index
+    numeric_parts = [int(idx.split()[0]) for idx in store_ids]
+    next_bus_id = max(self.network.buses.index.astype(int)) + 1
+    next_store_id = max(numeric_parts) + 1
+    bus_id_counter = count(start=next_bus_id)
+    store_id_counter = count(start=next_store_id)
+
+    for bus in PtH2_AC_buses:
+        power_to_h2_links = self.network.links[
+            (self.network.links.carrier == "power_to_H2")
+            & (self.network.links.bus0 == bus)
+        ].index
+        power_to_heat_links = self.network.links[
+            (self.network.links.carrier == "PtH2_waste_heat")
+            & (self.network.links.bus0 == bus)
+        ].index
+        power_to_o2_links = self.network.links[
+            (self.network.links.carrier == "PtH2_O2")
+            & (self.network.links.bus0 == bus)
+        ].index
+
+        link_data = []
+        if len(power_to_heat_links) > 0:
+            link_data.append(
+                {
+                    "links": power_to_heat_links,
+                    "efficiency": 0.2,
+                    "carrier": "waste_heat",
+                }
+            )
+        if len(power_to_o2_links) > 0:
+            link_data.append(
+                {
+                    "links": power_to_o2_links,
+                    "efficiency": 0.015,
+                    "carrier": "O2",
+                }
+            )
+
+        for idx, data in enumerate(link_data, start=2):
+            links = data["links"]
+            efficiency = data["efficiency"]
+            carrier = data["carrier"]
+            new_bus_id = next(bus_id_counter)
+            new_store_id = next(store_id_counter)
+            store_name = f"{new_store_id} PtH2_{carrier}"
+
+            # Add new heat/o2-bus and -store
+            self.network.add(
+                "Bus", name=new_bus_id, carrier=f"PtH2_extra_bus_{carrier}"
+            )
+            self.network.add(
+                "Store",
+                name=store_name,
+                bus=new_bus_id,
+                carrier=f"PtH2_{carrier}",
+                e_nom_extendable=True,
+                e_nom=100000,
+                marginal_cost=0,
+                capital_cost=0,
+                e_cyclic=False,
+                standing_loss=1,
+            )
+
+            self.network.stores.loc[store_name, "scn_name"] = self.args[
+                "scn_name"
+            ]
+            self.network.buses.loc[str(new_bus_id), "scn_name"] = self.args[
+                "scn_name"
+            ]
+            self.network.buses.loc[str(new_bus_id), "country"] = "DE"
+
+            # set coordinates for new bus for the clustering method
+            bus0_index = self.network.links.loc[power_to_h2_links, "bus0"]
+            self.network.buses.loc[str(new_bus_id), "x"] = (
+                self.network.buses.loc[bus0_index, "x"].values[0]
+            )
+            self.network.buses.loc[str(new_bus_id), "y"] = (
+                self.network.buses.loc[bus0_index, "y"].values[0]
+            )
+
+            # Connect multiple link with new heat/o2-bus
+            self.network.links.loc[power_to_h2_links, f"bus{idx}"] = str(
+                new_bus_id
+            )
+            self.network.links.loc[power_to_h2_links, f"efficiency{idx}"] = (
+                efficiency
+            )
+
+            # Adjust originals waste_heat- and oxygen-links with new bus0
+            for link in links:
+                self.network.links.loc[link, "bus0"] = str(new_bus_id)