Merge pull request #710 from openego/features/implement-sclopf-pypsa0262

Implement sclopf

etrago/appl.py

@@ -28,6 +28,7 @@
 import datetime
 import os
 import os.path
+import pandas as pd
 
 __copyright__ = (
     "Flensburg University of Applied Sciences, "
@@ -51,7 +52,7 @@
     "db": "egon-data",  # database session
     "gridversion": None,  # None for model_draft or Version number
     "method": {  # Choose method and settings for optimization
-        "type": "market_grid",  # type of optimization, 'lopf' or 'market_grid'
+        "type": "market_grid",  # type of optimization, 'lopf', 'sclopf' or 'market_grid'
         "n_iter": 1,  # abort criterion of iterative optimization, 'n_iter' or 'threshold'
         "pyomo": True,  # set if pyomo is used for model building
         "formulation": "pyomo",
@@ -187,8 +188,9 @@ def run_etrago(args, json_path):
         Choose method and settings for optimization.
         The provided dictionary can have the following entries:
 
-        * "lopf" : str
-            Type of optimization, currently only "lopf". Default: "lopf".
+        * "type" : str
+            Choose the type of optimization. Current options: "lopf", "sclopf"
+            or "market_grid". Default: "market_grid".
         * "n_iter" : int
             In case of extendable lines, several LOPFs have to be performed.
             You can either set "n_iter" and specify a fixed number of
@@ -736,6 +738,14 @@ def run_etrago(args, json_path):
     # skip snapshots
     etrago.skip_snapshots()
 
+    # Temporary drop DLR as it is currently not working with sclopf
+    if etrago.args["method"] != "lopf":
+        etrago.network.lines_t.s_max_pu = pd.DataFrame(
+            index=etrago.network.snapshots,
+            columns=etrago.network.lines.index,
+            data=1.0,
+        )
+
     # start linear optimal powerflow calculations
 
     etrago.network.storage_units.cyclic_state_of_charge = True

etrago/execute/__init__.py

@@ -418,6 +418,14 @@ def optimize(self):
 
         self.grid_optimization()
 
+    elif self.args["method"]["type"] == "sclopf":
+        self.sclopf(
+            post_lopf=False,
+            n_process=4,
+            delta=0.01,
+            n_overload=0,
+            div_ext_lines=False,
+        )
     else:
         print("Method not defined")
 
@@ -579,7 +587,8 @@ def import_gen_from_links(network, drop_small_capacities=True):
     df["control"] = "PV"
     df.reset_index(inplace=True)
 
-    df.index = df.bus + " " + df.carrier
+    if not df.empty:
+        df.index = df.bus + " " + df.carrier
 
     # Aggregate disptach time series for new generators
     gas_to_add["bus1_carrier"] = gas_to_add.bus + " " + gas_to_add.carrier

etrago/execute/grid_optimization.py

@@ -51,7 +51,13 @@ def grid_optimization(self):
     #         self.network.links.index.str.contains('ramp')
     #         ].index, inplace=True)
     logger.info("Start solving grid optimization model")
-    self.lopf()
+    self.sclopf(
+        post_lopf=False,
+        n_process=4,
+        delta=0.01,
+        n_overload=0,
+        div_ext_lines=False,
+    )
 
 
 def fix_chp_generation(self):