Full Illinois Model

Snakefile

@@ -28,7 +28,8 @@ rule retrieve_load:
 rule retrieve_existing_generators:
     output: 
         generators = "data/existing_generators.csv",
-        aggregated = "data/aggregated_generators.csv"
+        aggregated = "data/aggregated_generators.csv",
+        build_year = "data/build_year_generators.csv"
     script: "scripts/retrieve_generators.py"
 
 rule retrieve_renewable_profiles:
@@ -39,6 +40,11 @@ rule retrieve_renewable_profiles:
         solar = "data/time_series/solar.csv"
     script: "scripts/retrieve_renewables.py"
 
+rule retrieve_co2_emissions:
+    output: 
+        emissions = "data/technology_emissions.csv"
+    script: "scripts/retrieve_emissions.py"
+
 rule build_topology:
     input: 
         supply_regions = "data/spatial_data/supply_regions.shp"
@@ -63,17 +69,27 @@ rule add_electricity:
         costs = "data/technology_costs.csv",
         wind_profile = "data/time_series/wind.csv",
         solar_profile = "data/time_series/solar.csv",
-        base_network = "data/networks/base_network.nc"
+        base_network = "data/networks/base_network.nc",
+        emissions = "data/technology_emissions.csv",
+        build_years = "data/build_year_generators.csv"
     output:
         elec_network = "data/networks/electricity_network.nc",
         final_costs = "data/final_costs.csv"
     script: "scripts/add_electricity.py"
 
-
 rule solve_network:
     input:
         elec_network = "data/networks/electricity_network.nc"
     output:
-        solved_network = "results/networks/illinois_solved.nc",
-        dispatch_figure = "results/figures/illinois_dispatch.png"
-    script: "scripts/solve_network.py"
+        solved_network = "results/networks/illinois_solved.nc"
+    script: "scripts/solve_network.py"
+
+rule plot_results:
+    input:
+        solved_network = "results/networks/illinois_solved.nc"
+    output: 
+        dispatch_figure = "results/figures/illinois_dispatch.png",
+        capacity_figure = "results/figures/illinois_capacity.png",
+        emissions_figure = "results/figures/illinois_emissions.png",
+        active_units_figure = "results/figures/illinois_active_units.png"
+    script: "scripts/plot_results.py"

config.yml

@@ -1,27 +1,53 @@
 state_abbr: 'IL'
 
 geo_res: 'rto'  # accepts: 'rto' or 'county'
-time_res: 4  # hours 
+time_res: 1  # hours 
+load_filter: 60e3  # MW, load above this level will be removed as outliers.
+total_demand: 136e6  # Annual MWh in the first year
+load_growth: 0.00  # % annual growth
+load_share:
+  MISO-Z4: 0.33
+  ComEd: 0.67
 
-model_year: 2025
+model_years: [2025, 2030, 2035, 2040, 2045]
+plot_year: 2045
 
-solar_year: 2019  # NSRDB only accepts years 1998-2020
-wind_year: 2011  # WTK only goes from 2009-2013
+co2_limits:
+  2020: 31.25
+  2025: 25  # MT
+  2030: 18.75
+  2035: 12.5
+  2040: 6.25
+  2045: 0.0
+
+solar_years: [2016,2017,2018,2019,2020]  # NSRDB only accepts years 1998-2020
+wind_years: [2009,2010,2011,2012,2013]  # WTK only goes from 2009-2013
 load_year: 2019
 
 rto_subba: ['0004','CE']  # MISO-Z4, ComEd
 region_names: ['MISO-Z4','ComEd']
 
 discount_rate: 0.07
 
+growth_rates:  # MW/year
+  Nuclear: 1e3
+  Biomass: 1e3
+  Natural Gas: 10e3
+  Coal: 10e3
+  Petroleum: 10e3
+  Batteries: 5e3
+  Solar: 5e3
+  Wind: 5e3
+
 lifetime:
   Nuclear: 80
   Biomass: 40
   Natural Gas: 40
-  Coal: 40
+  Coal: 50
   Solar: 20
   Wind: 20
   Batteries: 15
+  Petroleum: 40
 
 atb_params:
     atb_year: 2022  # the ATB publication year
@@ -56,11 +82,10 @@ atb_params:
 lines:  # this is an assumption
   v_nom: 400.
   s_nom: 0.
-  x: 1  # reactance
-  r: 1  # resistance
+  x: 0.001  # reactance
+  r: 0.001  # resistance
   s_nom_extendable: True
 
-
 turbine_params:
   cut_in: 3.0  # m/s
   cut_out: 25.0  # m/s
@@ -69,35 +94,25 @@ turbine_params:
   rated_power: 2.75  # MW
   air_density: 1.225  # kg/m^3
 
-# extendable_techs:  [
-#       'CCAvgCF', 
-#       'CTAvgCF', 
-#       'Biopower', 
-#       'Land-Based Wind',
-#       'Utility PV', 
-#       'IGCCAvgCF', 
-#       'LWR',
-#       'NuclearSMR', 
-#       '10Hr Battery Storage', 
-#       '8Hr Battery Storage',
-#       '6Hr Battery Storage', 
-#       '4Hr Battery Storage', 
-#       '2Hr Battery Storage'
-#       ]
+nuclear_params:  # for existing reactors, 2023 NEI costs in context
+  capacity_factor: 0.972
+  capital_cost: 6.88 # $/MWh
+  fixed_om: 18.68 # $/MWh
+  fuel: 5.37  # $/MWh
 
 extendable_techs:  [
       # 'CCAvgCF', 
       # 'CTAvgCF', 
-      # 'Biopower', 
-      # 'Land-Based Wind',
-      # 'Utility PV', 
       # 'IGCCAvgCF', 
+      # 'Biopower', 
+      'Land-Based Wind',
+      'Utility PV', 
       # 'LWR',
       # 'NuclearSMR', 
       # '10Hr Battery Storage', 
       # '8Hr Battery Storage',
       # '6Hr Battery Storage', 
-      # '4Hr Battery Storage', 
+      '4Hr Battery Storage', 
       # '2Hr Battery Storage'
       ]
 

functions/nrel_data_api.py

@@ -21,7 +21,7 @@
               'year':2019,
               'leap_day':'false',
               'selector':'POINT',
-              'utc':'false',
+              'utc':'true',
               'interval':'60',
               'attr_list':['ghi']}