Gen x retrofit mit

.gitignore

@@ -31,4 +31,5 @@ test/Logs
 **/TDR_results/
 
 # Test files
-test/TDR/TDR_results_test
+test/TDR/TDR_results_test
+!test/TDR/TDR_results_true

CHANGELOG.md

@@ -32,7 +32,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 - New settings parameter, VirtualChargeDischargeCost to test script and VREStor example case. The PR 608 attempts to 
   introduce this parameter as cost of virtual charging and discharging to avoid unusual results (#608). 
 - New settings parameter, StorageVirtualDischarge, to turn storage virtual charging and discharging off if desired by the user (#638).
-
+- Add module to retrofit existing resources with new technologies (#600).
 
 ### Fixed
 - Set MUST_RUN=1 for RealSystemExample/small_hydro plants (#517).

docs/make.jl

@@ -54,6 +54,7 @@ pages = OrderedDict(
                 "Long Duration Hydro" => "Model_Reference/Resources/hydro_inter_period_linkage.md"
             ],
             "Must Run" => "Model_Reference/Resources/must_run.md",
+            "Retrofit" => "Model_Reference/Resources/retrofit.md",
             "Storage" => [
                 "Storage" => "Model_Reference/Resources/storage.md",
                 "Investment Charge" => "Model_Reference/Resources/investment_charge.md",

docs/src/Model_Reference/Resources/retrofit.md

@@ -1,5 +1,6 @@
-## Retrofit
+# Retrofit
 
-```@docs
-GenX.retrofit
+```@autodocs
+Modules = [GenX]
+Pages = ["retrofits.jl"]
 ```

docs/src/User_Guide/model_input.md

@@ -156,6 +156,15 @@ Each file contains cost and performance parameters for various generators and ot
 |cluster | Number of the cluster when representing multiple clusters of a given technology in a given region.  |
 |**Required if electrolyzer is included in the model**|
 |Qualified_Hydrogen_Supply| {0,1}, Indicates that generator or storage resources is eligible to supply electrolyzers in the same zone (used for hourly clean supply constraint)|
+|**Required for retrofitting**|
+|Can\_Retrofit | {0, 1}, Flag for resource (storage, generation) eligibility for retrofit.|
+||Can\_Retrofit = 1: eligible for retrofit. |
+||Can\_Retrofit = 0: not eligible for retrofit.|
+|Retrofit | {0, 1}, Flag for resource retrofit technologies (i.e., retrofit options, e.g. CCS retrofit for coal plants).|
+||Retrofit = 1: is a retrofit technology. |
+||Retrofit = 0: is not a retrofit technology.|
+|Retrofit\_Id | Unique identifier to group retrofittable source technologies with retrofit options inside the same zone.|
+|Retrofit\_Efficiency | [0,1], Efficiency of the retrofit technology.|
 
 ##### Table 5b: Settings-specific columns in all resource .csv file
 ---
@@ -523,6 +532,7 @@ In addition to the files described above, the `resources` folder can contain add
 |Min\_Retired\_Cap\_MW  |Minimum required discharge capacity retirements in the current model period. This field can be used to enforce lifetime retirements of existing capacity. Note that for co-located VRE-STOR resources, this value pertains to the grid connection. |
 |Min\_Retired\_Energy\_Cap\_MW  |Minimum required energy capacity retirements in the current model period. This field can be used to enforce lifetime retirements of existing energy capacity. Note that for co-located VRE-STOR resources, this value pertains to the storage component.|
 |Min\_Retired\_Charge\_Cap\_MW  |Minimum required energy capacity retirements in the current model period. This field can be used to enforce lifetime retirements of existing charge capacity. |
+|Contribute\_Min\_Retirement | {0, 1}, Flag to indicate whether the (retrofitting) resource can contribute to the minimum retirement requirement.|
 |**co-located VRE-STOR resources only**|
 |Capital\_Recovery\_Period_DC  |Capital recovery period (in years) used for determining overnight capital costs from annualized investment costs for the inverter component. |
 |Capital\_Recovery\_Period_Solar  |Capital recovery period (in years) used for determining overnight capital costs from annualized investment costs for the solar PV component. |

docs/src/User_Guide/multi_stage_input.md

@@ -27,6 +27,7 @@ Instead of one set of input files, there is one directory of input files that ne
 | Lifetime                     | The operational lifespan in years of this technology after which it must be retired.                                                                       |
 | Capital\_Recovery\_Period      | The technology-specific period in years over which initial capital costs must be recovered. Note that for the co-located VRE-STOR module, this value represents the grid connection component.                                                               |
 | WACC                         | The technology-specific weighted average cost of capital. Note that for the co-located VRE-STOR module, this value represents the grid connection component.                                                                                                 |
+|Contribute\_Min\_Retirement | {0, 1}, Flag to indicate whether the (retrofitting) resource can contribute to the minimum retirement requirement.|
 
 |                              |                                                                   **co-located VRE-STOR resources only**                                                                  |
 |:------------------------------|:----------------------------------------------------------------------------------------------------------------------------------------------------------|
@@ -59,6 +60,53 @@ Instead of one set of input files, there is one directory of input files that ne
 | WACC                      | The line-specific weighted average cost of capital.                                                                                                  |
 
 
+!!! note "Allowing retrofitted capacity to not contribute to minimum retirement requirements (`myopic=0` only)"
+    Special considerations must be taken into account when utilizing the retrofit module alongside multi-stage planning, particularly when using a non zero value for the `Min_Retired_Cap_MW` column in the `Resource_multistage_data.csv` file.
+    1. When assigning a non-zero value to the `Min_Retired_Cap_MW` column in the `Resource_multistage_data.csv` file, the user can specify whether the model should consider the retrofitted capacity to contribute to the minimum retirement requirement. This is done by setting the `Contribute_Min_Retirement` column to 1 for the retrofit options in the same retrofit cluster (i.e., same `Retrofit_Id`).
+    2. By default, the model assumes that retrofitted capacity contributes to fulfilling minimum retirement requirements.
+    3. Should users wish to exclude retrofitted capacity from contributing to minimum retirement requirements, they must set the `Contribute_Min_Retirement` column to 0 for **all** retrofit options within the same retrofit cluster (i.e., sharing the same `Retrofit_Id`).
+    4. It's important to note that this additional functionality is not currently supported when `myopic=1`. In this case, the retrofit options are only allowed to contribute to the minimum retirement requirement.
+
+    Example 1: Retrofitted capacity is allowed to contribute to the minimum retirement requirement (i.e., retrofit options in the same cluster (`Retrofit_Id = 1`) all have `Contribute_Min_Retirement = 1`):
+    ```
+    Thermal.csv
+
+    Resource          │ Zone  | Retrofit | Can_Retrofit | Retrofit_Id | Retrofit_Efficiency
+    String            │ Int64 | Int64    | Int64        | Int64       | Float64
+    ─────────────────-┼───────┼─────────────────────────┼────────────-┼────────────────────
+    coal_1            │ 1     │ 0        │ 1            │ 1           │ 0
+    20_NH3_retrofit_1 │ 1     │ 1        │ 0            │ 1           │ 0.85
+    20_NH3_retrofit_2 │ 1     │ 1        │ 0            │ 1           │ 0.85
+    ```
+
+    ```
+    Resource_multistage_data.csv
+
+    Resource          │ Min_Retired_Cap_MW | Contribute_Min_Retirement
+    String            │ Float64            | Float64 
+    ─────────────────-┼────────────────────┼──────────────────────────
+    coal_1            │ 4500               │ 0
+    20_NH3_retrofit_1 │ 0                  │ 1                         <---------
+    20_NH3_retrofit_2 │ 0                  │ 1                         <---------
+    ```
+
+    Example 2: Retrofitted capacity is not allowed to contribute to the minimum retirement requirement (i.e., none of the retrofit options in the same cluster (`Retrofit_Id = 1`) contribute to the minimum retirement requirement (`myopic=0`)):
+
+    Thermal.csv: same as Example 1.
+
+    ```
+    Resource_multistage_data.csv
+
+    Resource          │ Min_Retired_Cap_MW | Contribute_Min_Retirement
+    String            │ Float64            | Float64
+    ─────────────────-┼────────────────────┼──────────────────────────
+    coal_1            │ 4500               │ 0
+    20_NH3_retrofit_1 │ 0                  │ 0                         <---------
+    20_NH3_retrofit_2 │ 0                  │ 0                         <---------
+    ```
+
+    And the case where some retrofit options contribute to the minimum retirement requirement and some do not is not currently supported and will be addressed in a future release.
+
 !!! warning "Warning"
     If `New_Build` and `Can_Retire` are both set to 0, the model will not transfer built capacity from one stage to the next, but will instead set capacity to the value of existing capacity from the input files for each stage. Therefore, the user must ensure that the capacity is correctly set in the input files for each stage. Not following this guideline may result in incorrect or unexpected results, particularly when setting a a non-zero value for the `Min_Retired_Cap_MW` parameter. 
 

example_systems/3_three_zone_w_co2_capture/Run.jl

@@ -1,4 +1,3 @@
 using GenX
-using Gurobi
 
-run_genx_case!(dirname(@__FILE__), Gurobi.Optimizer)
+run_genx_case!(dirname(@__FILE__))

example_systems/8_three_zone_w_retrofit/README.md

@@ -0,0 +1,32 @@
+# Three Zones with Retrofit
+
+This is a one-year example with hourly resolution which contains zones representing Massachusetts, Connecticut, and Maine. The twenty-two represented resources include natural gas, solar PV, wind, lithium-ion battery, and coal power plants. 
+This examples shows the usage of the retrofit module of GenX, and the model will be allowed to retire as well as retrofit the existing coal power plants and replacing the coal with blue ammonia with 85% efficiency. 
+To run the model, first navigate to the example directory:
+
+- Using a Julia REPL:
+
+```bash
+$ julia
+julia> cd("example_systems/8_three_zone_w_retrofit/")
+```
+
+- Using a terminal or command prompt:
+```bash
+$ cd example_systems/8_three_zone_w_retrofit/
+``` 
+
+Next, ensure that your settings in `settings/genx_settings.yml` are correct. The linear clustering unit commitment method (settings["UCommit"] = 2) is used. The default settings use the solver `HiGHS`, time domain reduced input data (`TimeDomainReduction: 1`) and minimum capacity requirement policy (`MinCapReq: 1`) as specified in the `policies/Minimum_capacity_requirement.csv` file. Other optional policies include a capacity reserve margin, an energy share requirement (such as renewable portfolio standard (RPS) or clean electricity standard (CES) policies), a CO2 emissions cap, and a maximum capacity requirement policy (see the documentation for more details). A mass-based carbon cap of 50 gCO<sub>2</sub> per kWh is specified in the `policies/CO2_cap.csv` input file.
+
+Once the settings are confirmed, run the model with the `Run.jl` script in the example directory:
+
+- Using a Julia REPL (recommended)
+```julia
+julia> include("Run.jl")
+```
+- Using a terminal or command prompt:
+```bash
+$ julia Run.jl
+```
+
+Once the model has completed, results will write to the `results` directory.

example_systems/8_three_zone_w_retrofit/Run.jl

@@ -0,0 +1,3 @@
+using GenX
+
+run_genx_case!(dirname(@__FILE__))

example_systems/8_three_zone_w_retrofit/policies/CO2_cap.csv

@@ -0,0 +1,4 @@
+,Network_zones,CO_2_Cap_Zone_1,CO_2_Cap_Zone_2,CO_2_Cap_Zone_3,CO_2_Max_tons_MWh_1,CO_2_Max_tons_MWh_2,CO_2_Max_tons_MWh_3,CO_2_Max_Mtons_1,CO_2_Max_Mtons_2,CO_2_Max_Mtons_3
+MA,z1,1,0,0,0.05,0,0,0.018,0,0
+CT,z2,0,1,0,0,0.05,0,0,0.025,0
+ME,z3,0,0,1,0,0,0.05,0,0,0.025

example_systems/8_three_zone_w_retrofit/policies/Capacity_reserve_margin.csv

@@ -0,0 +1,4 @@
+,Network_zones,CapRes_1
+MA,z1,0.156
+CT,z2,0.156
+ME,z3,0.156

example_systems/8_three_zone_w_retrofit/policies/Energy_share_requirement.csv

@@ -0,0 +1,4 @@
+,Network_zones,ESR_1,ESR_2
+MA,z1,0.259,0.348
+CT,z2,0.44,0.44
+ME,z3,0.776,0.776

example_systems/8_three_zone_w_retrofit/policies/Minimum_capacity_requirement.csv

@@ -0,0 +1,4 @@
+MinCapReqConstraint,ConstraintDescription,Min_MW
+1,MA_PV,5000
+2,CT_Wind,10000
+3,All_Batteries,6000

example_systems/8_three_zone_w_retrofit/resources/Storage.csv

@@ -0,0 +1,4 @@
+Resource,Zone,Model,LDS,New_Build,Can_Retire,Existing_Cap_MW,Existing_Cap_MWh,Existing_Charge_Cap_MW,Max_Cap_MW,Max_Cap_MWh,Max_Charge_Cap_MW,Min_Cap_MW,Min_Cap_MWh,Min_Charge_Cap_MW,Inv_Cost_per_MWyr,Inv_Cost_per_MWhyr,Inv_Cost_Charge_per_MWyr,Fixed_OM_Cost_per_MWyr,Fixed_OM_Cost_per_MWhyr,Fixed_OM_Cost_Charge_per_MWyr,Var_OM_Cost_per_MWh,Var_OM_Cost_per_MWh_In,Heat_Rate_MMBTU_per_MWh,Fuel,Cap_Size,Start_Cost_per_MW,Start_Fuel_MMBTU_per_MW,Ramp_Up_Percentage,Ramp_Dn_Percentage,Hydro_Energy_to_Power_Ratio,Min_Power,Self_Disch,Eff_Up,Eff_Down,Min_Duration,Max_Duration,Reg_Max,Rsv_Max,Reg_Cost,Rsv_Cost,MGA,Resource_Type,region,cluster
+MA_battery,1,1,0,1,0,0,0,0,-1,-1,-1,0,0,0,19584,22494,0,4895,5622,0,0.15,0.15,0.0,None,0,0,0,1.0,1.0,0,0.0,0,0.92,0.92,1,10,0.0,0.0,0,0,0,battery_mid,MA,0
+CT_battery,2,1,0,1,0,0,0,0,-1,-1,-1,0,0,0,19584,22494,0,4895,5622,0,0.15,0.15,0.0,None,0,0,0,1.0,1.0,0,0.0,0,0.92,0.92,1,10,0.0,0.0,0,0,0,battery_mid,CT,0
+ME_battery,3,1,0,1,0,0,0,0,-1,-1,-1,0,0,0,19584,22494,0,4895,5622,0,0.15,0.15,0.0,None,0,0,0,1.0,1.0,0,0.0,0,0.92,0.92,1,10,0.0,0.0,0,0,0,battery_mid,ME,0

example_systems/8_three_zone_w_retrofit/resources/Thermal.csv

@@ -0,0 +1,16 @@
+Resource,Zone,Model,Retrofit_Id,New_Build,Can_Retire,Retrofit,Can_Retrofit,Retrofit_efficiency,Existing_Cap_MW,Max_Cap_MW,Min_Cap_MW,Inv_Cost_per_MWyr,Fixed_OM_Cost_per_MWyr,Var_OM_Cost_per_MWh,Heat_Rate_MMBTU_per_MWh,Fuel,Cap_Size,Start_Cost_per_MW,Start_Fuel_MMBTU_per_MW,Up_Time,Down_Time,Ramp_Up_Percentage,Ramp_Dn_Percentage,Min_Power,Reg_Max,Rsv_Max,Reg_Cost,Rsv_Cost,MGA,Resource_Type,region,cluster
+MA_natural_gas_combined_cycle,1,1,0,1,0,0,0,0,0,-1,0,65400,10287,3.55,7.43,MA_NG,250,91,2,6,6,0.64,0.64,0.468,0.25,0.5,0,0,1,natural_gas_fired_combined_cycle,MA,1
+CT_natural_gas_combined_cycle,2,1,0,1,0,0,0,0,0,-1,0,65400,9698,3.57,7.12,CT_NG,250,91,2,6,6,0.64,0.64,0.338,0.133332722,0.266665444,0,0,1,natural_gas_fired_combined_cycle,CT,1
+ME_natural_gas_combined_cycle,3,1,0,1,0,0,0,0,0,-1,0,65400,16291,4.5,12.62,ME_NG,250,91,2,6,6,0.64,0.64,0.474,0.033333333,0.066666667,0,0,1,natural_gas_fired_combined_cycle,ME,1
+MA_coal_super_1,1,1,1,0,1,0,1,0,6000,0,0,112000,10287,0,7.43,coal,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,Coal_super,MA,1
+MA_coal_super_2,1,1,2,0,1,0,1,0,6000,0,0,112000,10287,0,7.43,coal,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,Coal_super,MA,1
+CT_coal_super_1,1,1,1,0,1,0,1,0,6000,0,0,112000,10287,0,7.43,coal,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,Coal_super,CT,1
+CT_coal_super_2,1,1,2,0,1,0,1,0,6000,0,0,112000,10287,0,7.43,coal,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,Coal_super,CT,1
+ME_coal_super_1,1,1,1,0,1,0,1,0,6000,0,0,112000,10287,0,7.43,coal,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,Coal_super,ME,1
+ME_coal_super_2,1,1,2,0,1,0,1,0,6000,0,0,112000,10287,0,7.43,coal,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,Coal_super,ME,1
+MA_20A_retrofit_1,1,1,1,1,0,1,0,0.85,0,-1,0,2342,10287,0,5,blue_ammonia,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,ammonia,MA,1
+MA_20A_retrofit_2,1,1,2,1,0,1,0,0.85,0,-1,0,2342,10287,0,5,blue_ammonia,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,ammonia,MA,1
+CT_20A_retrofit_1,1,1,1,1,0,1,0,0.85,0,-1,0,2342,10287,0,5,blue_ammonia,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,ammonia,CT,1
+CT_20A_retrofit_2,1,1,2,1,0,1,0,0.85,0,-1,0,2342,10287,0,5,blue_ammonia,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,ammonia,CT,1
+ME_20A_retrofit_1,1,1,1,1,0,1,0,0.85,0,-1,0,2342,10287,0,5,blue_ammonia,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,ammonia,ME,1
+ME_20A_retrofit_2,1,1,2,1,0,1,0,0.85,0,-1,0,2342,10287,0,5,blue_ammonia,300,17,2,4,4,0.64,0.64,0.55,0.25,0.25,0,0,1,ammonia,ME,1

example_systems/8_three_zone_w_retrofit/resources/Vre.csv

@@ -0,0 +1,5 @@
+Resource,Zone,LDS,Num_VRE_Bins,New_Build,Can_Retire,Existing_Cap_MW,Max_Cap_MW,Min_Cap_MW,Inv_Cost_per_MWyr,Fixed_OM_Cost_per_MWyr,Var_OM_Cost_per_MWh,Heat_Rate_MMBTU_per_MWh,Fuel,Cap_Size,Start_Cost_per_MW,Start_Fuel_MMBTU_per_MW,Ramp_Up_Percentage,Ramp_Dn_Percentage,Hydro_Energy_to_Power_Ratio,Min_Power,Reg_Max,Rsv_Max,Reg_Cost,Rsv_Cost,MGA,Resource_Type,region,cluster
+MA_solar_pv,1,0,1,1,0,0,-1,0,85300,18760,0.0,9.13,None,0,0,0,1.0,1.0,0,0.0,0.0,0.0,0,0,1,solar_photovoltaic,MA,1
+CT_onshore_wind,2,0,1,1,0,0,-1,0,97200,43205,0.1,9.12,None,0,0,0,1.0,1.0,0,0.0,0.0,0.0,0,0,1,onshore_wind_turbine,CT,1
+CT_solar_pv,2,0,1,1,0,0,-1,0,85300,18760,0.0,9.16,None,0,0,0,1.0,1.0,0,0.0,0.0,0.0,0,0,1,solar_photovoltaic,CT,1
+ME_onshore_wind,3,0,1,1,0,0,-1,0,97200,43205,0.1,9.12,None,0,0,0,1.0,1.0,0,0.0,0.0,0.0,0,0,1,onshore_wind_turbine,ME,1

example_systems/8_three_zone_w_retrofit/resources/policy_assignments/Resource_minimum_capacity_requirement.csv

@@ -0,0 +1,11 @@
+Resource,Min_Cap_1,Min_Cap_2,Min_Cap_3
+MA_natural_gas_combined_cycle,0,0,0
+MA_solar_pv,1,0,0
+CT_natural_gas_combined_cycle,0,0,0
+CT_onshore_wind,0,1,0
+CT_solar_pv,0,0,0
+ME_natural_gas_combined_cycle,0,0,0
+ME_onshore_wind,0,0,0
+MA_battery,0,0,1
+CT_battery,0,0,1
+ME_battery,0,0,1

example_systems/8_three_zone_w_retrofit/settings/clp_settings.yml

@@ -0,0 +1,14 @@
+# Clp Solver parameters https://github.com/jump-dev/Clp.jl
+# Common solver settings
+Feasib_Tol: 1e-5              # Primal/Dual feasibility tolerance
+TimeLimit: -1.0               # Terminate after this many seconds have passed. A negative value means no time limit
+Pre_Solve: 0                  # Set to 1 to disable presolve
+Method: 5                     # Solution method: dual simplex (0), primal simplex (1), sprint (2), barrier with crossover (3), barrier without crossover (4), automatic (5)
+
+#Clp-specific solver settings
+DualObjectiveLimit: 1e308     # When using dual simplex (where the objective is monotonically changing), terminate when the objective exceeds this limit
+MaximumIterations: 2147483647 # Terminate after performing this number of simplex iterations
+LogLevel: 1                   # Set to 1, 2, 3, or 4 for increasing output. Set to 0 to disable output
+InfeasibleReturn: 0           # Set to 1 to return as soon as the problem is found to be infeasible (by default, an infeasibility proof is computed as well)
+Scaling: 3                    # 0 -off, 1 equilibrium, 2 geometric, 3 auto, 4 dynamic(later)
+Perturbation: 100             # switch on perturbation (50), automatic (100), don't try perturbing (102)

example_systems/8_three_zone_w_retrofit/settings/cplex_settings.yml

@@ -0,0 +1,10 @@
+# CPLEX Solver Parameters
+Feasib_Tol: 1.0e-05 # Constraint (primal) feasibility tolerances.
+Optimal_Tol: 1e-5 # Dual feasibility tolerances.
+Pre_Solve: 1 # Controls presolve level.
+TimeLimit: 110000 # Limits total time solver.
+MIPGap: 1e-3 # Relative (p.u. of optimal) mixed integer optimality tolerance for MIP problems (ignored otherwise).
+Method: 2 # Algorithm used to solve continuous models (including MIP root relaxation).
+BarConvTol: 1.0e-08 # Barrier convergence tolerance (determines when barrier terminates).
+NumericFocus: 0 # Numerical precision emphasis.
+SolutionType: 2 # Solution type for LP or QP.