Merge pull request #127 from REAM-lab/database/transmission_fix

Script to fix derating factor for transmission lines
switch-model · Jan 29, 2023 · b73e37f · b73e37f
2 parents 2921d7a + 0f435df
commit b73e37f
Show file tree

Hide file tree

Showing 3 changed files with 106 additions and 12 deletions.
diff --git a/database/2022-10-12_Transmission_scenario_fix.sql b/database/2022-10-12_Transmission_scenario_fix.sql
@@ -0,0 +1,78 @@
+-- October 12, 2022
+-- Author: pesap
+
+-- Reason: We detected a type in the transmission derating factor for all the transmission
+-- lines for the WECC model. This is the description in the switch module of what the
+-- transmission derate factor does: trans_derating_factor[tx in TRANSMISSION_LINES] is an
+-- overall derating factor for each transmission line that can reflect forced outage
+-- rates, stability or contingency limitations. This parameter is optional and defaults to
+-- 1. This parameter should be in the range of 0 to 1, being 0 a value that disables the
+-- line completely.
+-- In our scenarios we have been using 0.59 for no apparent reason and this limited the
+-- amount of electricity that flow throught the transmission line. The new value that we
+-- will use for the scenarios is 0.95 or a 5% derate.
+
+alter table 
+  switch.transmission_lines 
+add 
+  column transmission_scenario_id int;
+
+alter table 
+  switch.scenario 
+add 
+  column transmission_scenario_id int;
+
+UPDATE switch.transmission_lines SET transmission_scenario_id = 1;
+UPDATE switch.scenario SET transmission_scenario_id = 1;
+
+alter table switch.transmission_lines drop constraint transmission_lines_pkey;
+
+insert into switch.transmission_lines(
+  transmission_line_id, start_load_zone_id,
+  end_load_zone_id, trans_length_km,
+  trans_efficiency, existing_trans_cap_mw,
+  new_build_allowed, derating_factor,
+  terrain_multiplier, transmission_cost_econ_multiplier,
+  transmission_scenario_id
+)
+select
+  transmission_line_id,
+  start_load_zone_id,
+  end_load_zone_id,
+  trans_length_km,
+  trans_efficiency,
+  existing_trans_cap_mw,
+  new_build_allowed,
+  0.95 as derating_factor,
+  terrain_multiplier,
+  transmission_cost_econ_multiplier,
+  2 as transmission_scenario_id
+from
+  switch.transmission_lines;
+
+insert into switch.scenario
+values (179,
+		'[LDES and SWITCH_WEAP]',
+		'Baseline scenario with 0 carbon cap for WECC and CA by 2045, 15% PRM, updated 2 days per month sampling and tranmission_line scenario',
+		17, -- study_timeframe_id
+		16, -- time_sample_id
+		115, -- demand_scenario_id
+		4, -- fuel_simple_price_scenario, without Bio_Solid costs, because they are provided by supply curve
+		23, -- generation_plant_scenario_id
+		25, -- generation_plant_cost_scenario_id
+		21, -- generation_plant_existing_and_planned_scenario_id
+		23, -- hydro_simple_scenario_id
+		92, -- carbon_cap_scenario_id
+		2, -- supply_curve_scenario_id
+		1, -- regional_fuel_market_scenario_id
+		NULL, -- rps_scenario_id
+		NULL, --enable_dr
+		NULL, --enable_ev
+		2, --transmission_base_capital_cost_scenario_id
+		NULL, --ca_policies_scenario_id
+		1, --enable_planning_reserves
+		2, --generation_plant_technologies_scenario_id
+		3, --variable_o_m_cost_scenario_id
+		NULL,--wind_to_solar_ratio
+    2 --tranmission_scenario_id
+		);
diff --git a/switch_model/wecc/get_inputs/get_inputs.py b/switch_model/wecc/get_inputs/get_inputs.py
@@ -315,16 +315,32 @@ def query_db(config, skip_cf):
             "trans_terrain_multiplier",
             "trans_new_build_allowed",
         ],
-        """
-         SELECT start_load_zone_id || '-' || end_load_zone_id, t1.name, t2.name,
-             trans_length_km, trans_efficiency, existing_trans_cap_mw, transmission_line_id,
-            derating_factor, terrain_multiplier * transmission_cost_econ_multiplier as terrain_multiplier,
-            new_build_allowed
-         FROM transmission_lines
-             join load_zone as t1 on(t1.load_zone_id=start_load_zone_id)
-             join load_zone as t2 on(t2.load_zone_id=end_load_zone_id)
-         WHERE start_load_zone_id < end_load_zone_id
-         ORDER BY 2,3;
+        f"""
+        SELECT
+          start_load_zone_id || '-' || end_load_zone_id,
+          t1.name,
+          t2.name,
+          trans_length_km,
+          trans_efficiency,
+          existing_trans_cap_mw,
+          transmission_line_id,
+          derating_factor,
+          terrain_multiplier * transmission_cost_econ_multiplier as terrain_multiplier,
+          new_build_allowed
+        FROM
+          switch.transmission_lines
+          join switch.load_zone as t1 on(
+            t1.load_zone_id = start_load_zone_id
+          )
+          join switch.load_zone as t2 on(
+            t2.load_zone_id = end_load_zone_id
+          )
+        WHERE
+          start_load_zone_id < end_load_zone_id and
+          transmission_scenario_id = {params.transmission_scenario_id}
+        ORDER BY
+          2,
+          3;
          """,
     )
 
@@ -337,10 +353,9 @@ def query_db(config, skip_cf):
             "trans_lifetime_yrs",
             "trans_fixed_om_fraction",
         ],
-        # See Issue #80 for reasoning behind the 85 year lifetime.
         f"""
         SELECT trans_capital_cost_per_mw_km,
-            85 as trans_lifetime_yrs,
+            25 as trans_lifetime_yrs,
             0.03 as trans_fixed_om_fraction
         FROM transmission_base_capital_cost
         WHERE transmission_base_capital_cost_scenario_id = {params.transmission_base_capital_cost_scenario_id}

diff --git a/switch_model/wecc/get_inputs/scenario.py b/switch_model/wecc/get_inputs/scenario.py
@@ -25,6 +25,7 @@ def __init__(self, scenario_id=None):
         self.generation_plant_technologies_scenario_id = None
         self.variable_o_m_cost_scenario_id = None
         self.wind_to_solar_ratio = None
+        self.transmission_scenario_id = None
 
 
 def load_scenario_from_config(config, db_cursor) -> ScenarioParams: