Include interest area in AC clustering

etrago/analyze/plot.py

@@ -51,6 +51,7 @@
     import tilemapbase
 
     from etrago.execute import import_gen_from_links
+    from etrago.tools.utilities import find_buses_area
 
 __copyright__ = (
     "Flensburg University of Applied Sciences, "
@@ -2476,6 +2477,8 @@ def plot_grid(
         * 'q_flow_max': maximal reactive flows
         * 'dlr': energy above nominal capacity
         * 'grey': plot all lines and DC links grey colored
+        * 'interest_area': plot all AC buses inside the interest area larger
+          than buses outside the interest area
 
     bus_sizes : float, optional
         Size of buses. The default is 0.001.
@@ -2909,6 +2912,22 @@ def plot_grid(
         bus_sizes[bus_sizes != "AC"] = 0
         bus_sizes[bus_sizes == "AC"] = 1 * bus_scaling
         bus_scaling = bus_sizes
+    elif bus_colors == "interest_area":
+        bus_scaling = bus_sizes
+        # only plot AC buses
+        bus_sizes = pd.Series(
+            data=network.buses.carrier, index=network.buses.index
+        )
+        bus_sizes[bus_sizes != "AC"] = 0
+        bus_sizes[bus_sizes == "AC"] = 1 * bus_scaling
+        # only plot buses inside interest area
+        buses_interest_area = find_buses_area(self, "AC")
+        buses_outside = [
+            _ for _ in bus_sizes.index if _ not in buses_interest_area
+        ]
+        bus_sizes.loc[buses_outside] = bus_sizes.loc[buses_outside] * 0.3
+        bus_scaling = bus_sizes
+        bus_colors = coloring()["AC"]
     else:
         logger.warning("bus_color {} undefined".format(bus_colors))
 

etrago/appl.py

@@ -85,8 +85,7 @@
     "model_formulation": "kirchhoff",  # angles or kirchhoff
     "scn_name": "eGon2035",  # scenario: eGon2035, eGon100RE or status2019
     # Scenario variations:
-    "scn_extension": None,  # None or array of extension scenarios
-    "scn_decommissioning": None,  # None or decommissioning scenario
+    "scn_extension": ["nep2021_c2035"],  # None or array of extension scenarios
     # Export options:
     "lpfile": False,  # save pyomo's lp file: False or /path/to/lpfile.lp
     "csv_export": "results",  # save results as csv: False or /path/tofolder
@@ -113,15 +112,17 @@
     "extra_functionality": {},  # Choose function name or {}
     # Spatial Complexity:
     "delete_dispensable_ac_buses": True,  # bool. Find and delete expendable buses
+    "interest_area": False,  # False, path to shapefile or list of nuts names of the area that is excluded from the clustering. By default the buses inside remain the same, but the parameter "n_cluster_interest_area" inside "network clustering" defines if it should be clustered to a certain number of buses.
     "network_clustering_ehv": {
         "active": False,  # choose if clustering of HV buses to EHV buses is activated
         "busmap": False,  # False or path to stored busmap
     },
     "network_clustering": {
         "active": True,  # choose if clustering is activated
         "method": "kmedoids-dijkstra",  # choose clustering method: kmeans or kmedoids-dijkstra
-        "n_clusters_AC": 30,  # total number of resulting AC nodes (DE+foreign)
+        "n_clusters_AC": 30,  # total number of resulting AC nodes (DE+foreign-interest_area)
         "cluster_foreign_AC": False,  # take foreign AC buses into account, True or False
+        "n_cluster_interest_area": False, # False or number of buses.
         "method_gas": "kmedoids-dijkstra",  # choose clustering method: kmeans or kmedoids-dijkstra
         "n_clusters_gas": 14,  # total number of resulting CH4 nodes (DE+foreign)
         "cluster_foreign_gas": False,  # take foreign CH4 buses into account, True or False
@@ -263,40 +264,20 @@ def run_etrago(args, json_path):
     scn_name : str
          Choose your scenario. Currently, there are two different
          scenarios: "eGon2035", "eGon100RE". Default: "eGon2035".
-    scn_extension : None or str
-        This option does currently not work!
+    scn_extension : None or list of str
 
         Choose extension-scenarios which will be added to the existing
-        network container. Data of the extension scenarios are located in
-        extension-tables (e.g. model_draft.ego_grid_pf_hv_extension_bus)
-        with the prefix 'extension\_'.
-        There are three overlay networks:
-
-        * 'nep2035_confirmed' includes all planed new lines confirmed by the
-          Bundesnetzagentur
-        * 'nep2035_b2' includes all new lines planned by the
-          Netzentwicklungsplan 2025 in scenario 2035 B2
-        * 'BE_NO_NEP 2035' includes planned lines to Belgium and Norway and
-          adds BE and NO as electrical neighbours
-
+        network container. In case new lines replace existing ones, these are
+        dropped from the network. Data of the extension scenarios is located in
+        extension-tables (e.g. grid.egon_etrago_extension_line)
+        There are two overlay networks:
+
+        * 'nep2021_confirmed' includes all planed new lines confirmed by the
+          Bundesnetzagentur included in the NEP version 2021
+        * 'nep2021_c2035' includes all new lines planned by the
+          Netzentwicklungsplan 2021 in scenario 2035 C
         Default: None.
-    scn_decommissioning : NoneType or str
-        This option does currently not work!
-
-        Choose an extra scenario which includes lines you want to decommission
-        from the existing network. Data of the decommissioning scenarios are
-        located in extension-tables
-        (e.g. model_draft.ego_grid_pf_hv_extension_bus) with the prefix
-        'decommissioning\_'.
-        Currently, there are two decommissioning_scenarios which are linked to
-        extension-scenarios:
-
-        * 'nep2035_confirmed' includes all lines that will be replaced in
-          confirmed projects
-        * 'nep2035_b2' includes all lines that will be replaced in
-          NEP-scenario 2035 B2
 
-        Default: None.
     lpfile : bool or str
         State if and where you want to save pyomo's lp file. Options:
         False or '/path/tofile.lp'. Default: False.
@@ -404,6 +385,16 @@ def run_etrago(args, json_path):
         connected lines are merged. This reduces the spatial complexity without
         losing any accuracy.
         Default: True.
+
+    interest_area: False, list, string
+        Area of especial interest that will be not clustered, except when
+        n_cluster_interest_area is provided. It is by default set to false.
+        When an interest_area is provided, the given value for n_clusters_AC
+        will mean the total of AC buses outside the area.The area can be
+        provided in two ways: list of nuts names e.G.
+        ["Cuxhaven", "Bremerhaven", "Bremen"] or a string with a path to a
+        shape file (.shp).
+
     network_clustering_ehv : dict
         Choose if you want to apply an extra high voltage clustering to the
         electrical network.
@@ -435,12 +426,14 @@ def run_etrago(args, json_path):
             * "kmeans": considers geographical locations of buses
             * "kmedoids-dijkstra":  considers electrical distances between
             buses
-
             Default: "kmedoids-dijkstra".
-        * "n_clusters_AC" : int
+        * "n_clusters_AC" : int, False
             Defines total number of resulting AC nodes including DE and foreign
             nodes if `cluster_foreign_AC` is set to True, otherwise only DE
-            nodes.
+            nodes. When using the interest_area parameter, n_clusters_AC could
+            be set to False, which means that only the buses inside the 
+            provided area are clustered.
+
             Default: 30.
         * "cluster_foreign_AC" : bool
             If set to False, the AC buses outside Germany are not clustered
@@ -449,13 +442,18 @@ def run_etrago(args, json_path):
             as well and included in number of clusters specified through
             ``'n_clusters_AC'``.
             Default: False.
+
+        * "n_cluster_interest_area": False, int
+            Number of buses to cluster all the electrical buses in the area
+            of interest. Method provided in the arg "method" is used. If
+            it is set to False, the area of interest is not clustered.
+            Default: False.
         * "method_gas" : str
             Method used for gas clustering. You can choose between two
             clustering methods:
             * "kmeans": considers geographical locations of buses
             * "kmedoids-dijkstra":  considers 'electrical' distances between
             buses
-
             Default: "kmedoids-dijkstra".
         * "n_clusters_gas" : int
             Defines total number of resulting CH4 nodes including DE and

etrago/args.json

@@ -28,7 +28,6 @@
     "model_formulation": "kirchhoff",
     "scn_name": "eGon2035",
     "scn_extension": null,
-    "scn_decommissioning": null,
     "lpfile": false,
     "csv_export": "results",
     "extendable": {
@@ -62,6 +61,7 @@
     "generator_noise": 789456,
     "extra_functionality": {},
     "delete_dispensable_ac_buses": true,
+    "interest_area": false,
     "network_clustering_ehv": {
         "active": false,
         "busmap": false
@@ -71,6 +71,7 @@
         "method": "kmedoids-dijkstra",
         "n_clusters_AC": 30,
         "cluster_foreign_AC": false,
+        "n_cluster_interest_area": false,
         "method_gas": "kmedoids-dijkstra",
         "n_clusters_gas": 17,
         "cluster_foreign_gas": false,