Feature/#814 modeling on interest area

etrago/appl.py

@@ -48,16 +48,17 @@
     from etrago import Etrago
 
 args = {
+    "nuts_3_map" : "germany-de-nuts-3-regions.geojson",
     # Setup and Configuration:
-    "db": "egon-data",  # database session
+    "db": "egon-data-wam02",  # database session # "egon-data-wam02"
     "gridversion": None,  # None for model_draft or Version number
     "method": {  # Choose method and settings for optimization
         "type": "lopf",  # type of optimization, 'lopf' or 'sclopf'
         "n_iter": 4,  # abort criterion of iterative optimization, 'n_iter' or 'threshold'
         "formulation": "linopy",
         "market_optimization":
             {
-                "active": True,
+                "active": False,
                 "market_zones": "status_quo", # only used if type='market_grid'
                 "rolling_horizon": {# Define parameter of market optimization
                     "planning_horizon": 168, # number of snapshots in each optimization
@@ -72,19 +73,19 @@
         "q_allocation": "p_nom",  # allocate reactive power via 'p_nom' or 'p'
     },
     "start_snapshot": 1,
-    "end_snapshot": 168,
+    "end_snapshot": 8760,
     "solver": "gurobi",  # glpk, cplex or gurobi
     "solver_options": {
-        "BarConvTol": 1.0e-5,
-        "FeasibilityTol": 1.0e-5,
+        "BarConvTol": 1.0e-9,
+        "FeasibilityTol": 1.0e-9,
         "method": 2,
-        "crossover": 0,
+        "crossover": 1,
         "logFile": "solver_etrago.log",
         "threads": 4,
         "BarHomogeneous": 1,
     },
     "model_formulation": "kirchhoff",  # angles or kirchhoff
-    "scn_name": "eGon100RE",  # scenario: eGon2035, eGon100RE or status2019
+    "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
@@ -111,16 +112,19 @@
         },
     },
     "generator_noise": 789456,  # apply generator noise, False or seed number
-    "extra_functionality": {},  # Choose function name or {}
+    "extra_functionality": {
+       "add_chp_ratio_constraint": {}
+    },  # Choose function name or {}
     # Spatial Complexity:
+    "delete_dispensable_ac_buses": True,  # bool. Find and delete expendable 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
         "CPU_cores": 4,  # number of cores used during clustering, "max" for all cores available.
     },
     "network_clustering": {
         "method": {
-            "focus_region": None,  # None, shape-file or list with string for Kreise
+            "focus_region": ["Ingolstadt"],  # None, shape-file or list with string for Kreise
             "algorithm": "kmedoids-dijkstra",  # choose clustering method: kmeans or kmedoids-dijkstra
             "remove_stubs": False,  # remove stubs before kmeans clustering
             "use_reduced_coordinates": False,  # if True, do not average cluster coordinates in kmeans
@@ -133,7 +137,7 @@
         },
         "electricity_grid": {
             "active": True,  # choose if clustering is activated
-            "cluster_within_focus": False,  # False for no clustering within focus region
+            "cluster_within_focus": 1,  # False for no clustering within focus region
             "n_clusters": 30,  # total number of resulting AC nodes
             "cluster_foreign": False,  # take foreign AC buses into account, True or False
             "k_elec_busmap": False,  # False or path/to/busmap.csv
@@ -164,7 +168,7 @@
         "n_clusters": 5,  # number of periods - only relevant for 'typical_periods'
         "n_segments": 5,  # number of segments - only relevant for segmentation
     },
-    "skip_snapshots": 5,  # False or number of snapshots to skip
+    "skip_snapshots": 3,  # False or number of snapshots to skip
     "temporal_disaggregation": {
         "active": False,  # choose if temporally full complex dispatch optimization should be conducted
         "no_slices": 8,  # number of subproblems optimization is divided into
@@ -662,6 +666,8 @@ def run_etrago(args, json_path):
     """
     etrago = Etrago(args, json_path=json_path)
 
+    print(datetime.datetime.now())
+
     # import network from database
     etrago.build_network_from_db()
 
@@ -681,6 +687,53 @@ def run_etrago(args, json_path):
     # skip snapshots
     etrago.skip_snapshots()
 
+    # set interest components to extendable and add additional components
+
+    etrago.print_capital_costs()
+
+    etrago.adjust_capital_costs()
+
+    etrago.print_capital_costs()
+
+    etrago.add_extendable_solar_generators_to_interest_area()
+
+    etrago.reset_gas_CHP_capacities()
+
+    etrago.replace_gas_links_with_extendable()
+
+    etrago.add_biogas_CHP_extendable()
+
+    etrago.add_biomass_CHP_extendable()
+
+    etrago.add_extendable_heat_pumps_to_interest_area()
+
+    etrago.add_waste_CHP_ingolstadt()
+
+    #etrago.add_biomass_boiler_extendable()
+
+    etrago.set_battery_parameter_interest_area()
+    #etrago.set_battery_and_heat_store_parameters_interest_area()
+
+    etrago.set_cyclic_constraints()
+
+    buses_ing = etrago.find_interest_buses()
+
+    gens_ing = etrago.network.generators[etrago.network.generators.bus.isin(buses_ing.index)]
+    gcolumns = ["bus", "carrier", "p_nom", "marginal_cost", "capital_cost", "p_nom_extendable"]
+
+    links_ing = etrago.find_links_connected_to_interest_buses()
+    lcolumns = ["bus0", "bus1", "carrier", "p_nom", "p_nom_extendable", "capital_cost", "marginal_cost"]
+
+    with pd.option_context("display.max_columns", None):
+        print(buses_ing)
+        print(links_ing[lcolumns])
+        print(gens_ing[gcolumns])
+
+    etrago.network.export_to_netcdf("base_network_Scenario_1f.nc")
+
+    #import pdb
+    #pdb.set_trace()
+
     # start linear optimal powerflow calculations
     etrago.optimize()
 

etrago/network.py

@@ -113,6 +113,21 @@
     set_random_noise,
     set_trafo_costs,
     update_busmap,
+    find_interest_buses,
+    find_links_connected_to_interest_buses,
+    add_extendable_solar_generators_to_interest_area,
+    add_extendable_heat_pumps_to_interest_area,
+    set_battery_parameter_interest_area,
+    set_battery_and_heat_store_parameters_interest_area,
+    add_waste_CHP_ingolstadt,
+    reset_gas_CHP_capacities,
+    add_biogas_CHP_extendable,
+    add_biomass_CHP_extendable,
+    add_biomass_boiler_extendable,
+    adjust_capital_costs,
+    print_capital_costs,
+    replace_gas_links_with_extendable,
+    set_cyclic_constraints
 )
 
 logger = logging.getLogger(__name__)
@@ -381,6 +396,36 @@ def __init__(
 
     levelize_abroad_inland_parameters = levelize_abroad_inland_parameters
 
+    find_interest_buses = find_interest_buses
+
+    find_links_connected_to_interest_buses = find_links_connected_to_interest_buses
+
+    add_extendable_solar_generators_to_interest_area = add_extendable_solar_generators_to_interest_area
+
+    add_extendable_heat_pumps_to_interest_area = add_extendable_heat_pumps_to_interest_area
+
+    set_battery_parameter_interest_area = set_battery_parameter_interest_area
+
+    set_battery_and_heat_store_parameters_interest_area = set_battery_and_heat_store_parameters_interest_area
+
+    add_waste_CHP_ingolstadt = add_waste_CHP_ingolstadt
+
+    reset_gas_CHP_capacities = reset_gas_CHP_capacities
+
+    add_biogas_CHP_extendable = add_biogas_CHP_extendable
+
+    add_biomass_CHP_extendable = add_biomass_CHP_extendable
+
+    add_biomass_boiler_extendable = add_biomass_boiler_extendable
+
+    adjust_capital_costs = adjust_capital_costs
+
+    print_capital_costs = print_capital_costs
+
+    replace_gas_links_with_extendable = replace_gas_links_with_extendable
+
+    set_cyclic_constraints = set_cyclic_constraints
+
     def dc_lines(self):
         return self.filter_links_by_carrier("DC", like=False)