Merge pull request #732 from yerbol-akhmetov/fix_380kV

fix hardcoded 380kv

Author: pz-max

config.default.yaml

@@ -118,6 +118,7 @@ load_options:
   scale: 1  # scales all load time-series, i.e. 2 = doubles load
 
 electricity:
+  base_voltage: 380.
   voltages: [220., 300., 380.]
   co2limit: 7.75e+7  # European default, 0.05 * 3.1e9*0.5, needs to be adjusted for Africa
   co2base: 1.487e+9  # European default, adjustment to Africa necessary

config.tutorial.yaml

@@ -132,6 +132,7 @@ load_options:
   scale: 1  # scales all load time-series, i.e. 2 = doubles load
 
 electricity:
+  base_voltage: 380.
   voltages: [220., 300., 380.]
   co2limit: 1.487e+9
   co2base: 1.487e+9

doc/configtables/electricity.csv

@@ -1,4 +1,5 @@
 ,Unit,Values,Description
+base_voltage, kV, float, "Base voltage to which all lines are simplified/aggregated. Simplification preserves transmission capacities."
 voltages, kV, "Any subset of {220., 300., 380.}", "Voltage levels considered."
 co2limit,:math:`t_{CO_2-eq}/a`, float, "Cap on system total annual carbon dioxide equivalent emissions."
 co2base,:math:`t_{CO_2-eq}/a`, float, "Reference value of system total annual carbon dioxide equivalent emissions. Used only if relative emission reduction target is specified in ``{opts}`` wildcard."

doc/release_notes.rst

@@ -111,6 +111,8 @@ PyPSA-Earth 0.2.0
 
 * Update instructions on how to write documentation. `PR #720 <https://github.com/pypsa-meets-earth/pypsa-earth/pull/720>`__
 
+* Fix hard-coded simplification of lines to 380kV. `PR #732 <https://github.com/pypsa-meets-earth/pypsa-earth/pull/732>`__
+
 PyPSA-Earth 0.1.0
 =================
 

scripts/simplify_network.py

@@ -110,21 +110,20 @@
 logger = logging.getLogger(__name__)
 
 
-def simplify_network_to_380(n, linetype):
+def simplify_network_to_base_voltage(n, linetype, base_voltage):
     """
-    Fix all lines to a voltage level of 380 kV and remove all transformers.
+    Fix all lines to a voltage level of base voltage level and remove all transformers.
     The function preserves the transmission capacity for each line while updating
     its voltage level, line type and number of parallel bundles (num_parallel).
     Transformers are removed and connected components are moved from their
     starting bus to their ending bus. The corresponding starting buses are
     removed as well.
     """
-    logger.info("Mapping all network lines onto a single 380kV layer")
+    logger.info(f"Mapping all network lines onto a single {int(base_voltage)}kV layer")
 
-    linetype_380 = linetype
-    n.lines["type"] = linetype_380
-    n.lines["v_nom"] = 380
-    n.lines["i_nom"] = n.line_types.i_nom[linetype_380]
+    n.lines["type"] = linetype
+    n.lines["v_nom"] = base_voltage
+    n.lines["i_nom"] = n.line_types.i_nom[linetype]
     # Note: s_nom is set in base_network
     n.lines["num_parallel"] = n.lines.eval("s_nom / (sqrt(3) * v_nom * i_nom)")
 
@@ -715,7 +714,8 @@ def merge_isolated_nodes(n, threshold, aggregation_strategies=dict()):
 
     n = pypsa.Network(snakemake.input.network)
 
-    linetype = snakemake.config["lines"]["types"][380.0]
+    base_voltage = snakemake.config["electricity"]["base_voltage"]
+    linetype = snakemake.config["lines"]["types"][base_voltage]
 
     aggregation_strategies = snakemake.config["cluster_options"].get(
         "aggregation_strategies", {}
@@ -725,7 +725,7 @@ def merge_isolated_nodes(n, threshold, aggregation_strategies=dict()):
         p: {k: getattr(pd.Series, v) for k, v in aggregation_strategies[p].items()}
         for p in aggregation_strategies.keys()
     }
-    n, trafo_map = simplify_network_to_380(n, linetype)
+    n, trafo_map = simplify_network_to_base_voltage(n, linetype, base_voltage)
 
     Nyears = n.snapshot_weightings.objective.sum() / 8760