ReEDS-2.0/d2_unfix_op.gms at main · NatLabRockies/ReEDS-2.0 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
* load
LOAD.lo(r,h,t_unfix) = 0 ;
LOAD.up(r,h,t_unfix) = +inf ;
FLEX.lo(flex_type,r,h,t_unfix)$Sw_EFS_flex  = 0 ;
FLEX.up(flex_type,r,h,t_unfix)$Sw_EFS_flex  = +inf ;
DROPPED.lo(r,h,t_unfix)$[(yeart(t_unfix)<Sw_StartMarkets) or (Sw_PCM=1)] = 0 ;
DROPPED.up(r,h,t_unfix)$[(yeart(t_unfix)<Sw_StartMarkets) or (Sw_PCM=1)] = +inf ;
EXCESS.lo(r,h,t_unfix)$[(yeart(t_unfix)<Sw_StartMarkets) or (Sw_PCM=1)] = 0 ;
EXCESS.up(r,h,t_unfix)$[(yeart(t_unfix)<Sw_StartMarkets) or (Sw_PCM=1)] = +inf ;

* generation and storage
GEN.lo(i,v,r,h,t_unfix)$valgen(i,v,r,t_unfix) = 0 ;
GEN.up(i,v,r,h,t_unfix)$valgen(i,v,r,t_unfix) = +inf ;
GEN_PLANT.lo(i,v,r,h,t_unfix)$[storage_hybrid(i)$(not csp(i))$valgen(i,v,r,t_unfix)$Sw_HybridPlant] = 0 ;
GEN_PLANT.up(i,v,r,h,t_unfix)$[storage_hybrid(i)$(not csp(i))$valgen(i,v,r,t_unfix)$Sw_HybridPlant] = +inf ;
GEN_STORAGE.lo(i,v,r,h,t_unfix)$[storage_hybrid(i)$(not csp(i))$valgen(i,v,r,t_unfix)$Sw_HybridPlant] = 0 ;
GEN_STORAGE.up(i,v,r,h,t_unfix)$[storage_hybrid(i)$(not csp(i))$valgen(i,v,r,t_unfix)$Sw_HybridPlant] = +inf ;
CURT.lo(r,h,t_unfix)$Sw_CurtMarket = 0 ;
CURT.up(r,h,t_unfix)$Sw_CurtMarket = +inf ;
MINGEN.lo(r,szn,t_unfix)$Sw_Mingen = 0 ;
MINGEN.up(r,szn,t_unfix)$Sw_Mingen = +inf ;
STORAGE_IN.lo(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$(storage_standalone(i) or hyd_add_pump(i))] = 0 ;
STORAGE_IN.up(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$(storage_standalone(i) or hyd_add_pump(i))] = +inf ;
STORAGE_IN_PLANT.lo(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$storage_hybrid(i)$(not csp(i))$Sw_HybridPlant] = 0 ;
STORAGE_IN_PLANT.up(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$storage_hybrid(i)$(not csp(i))$Sw_HybridPlant] = +inf ;
STORAGE_IN_GRID.lo(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$storage_hybrid(i)$(not csp(i))$Sw_HybridPlant] = 0 ;
STORAGE_IN_GRID.up(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$storage_hybrid(i)$(not csp(i))$Sw_HybridPlant] = +inf ;
STORAGE_LEVEL.lo(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$storage(i)] = 0 ;
STORAGE_LEVEL.up(i,v,r,h,t_unfix)$[valgen(i,v,r,t_unfix)$storage(i)] = +inf ;
STORAGE_INTERDAY_LEVEL.lo(i,v,r,allszn,t_unfix)$[valgen(i,v,r,t_unfix)$storage_interday(i)$numpartitions(allszn)] = 0 ;
STORAGE_INTERDAY_LEVEL.up(i,v,r,allszn,t_unfix)$[valgen(i,v,r,t_unfix)$storage_interday(i)$numpartitions(allszn)] = +inf ;
AVAIL_SITE.lo(x,h,t_unfix)$[Sw_SpurScen$xfeas(x)] = 0 ;
AVAIL_SITE.up(x,h,t_unfix)$[Sw_SpurScen$xfeas(x)] = +inf ;
RAMPUP.lo(i,r,h,hh,t_unfix)$[Sw_StartCost$startcost(i)$numhours_nexth(h,hh)$valgen_irt(i,r,t_unfix)] = 0 ;
RAMPUP.up(i,r,h,hh,t_unfix)$[Sw_StartCost$startcost(i)$numhours_nexth(h,hh)$valgen_irt(i,r,t_unfix)] = +inf ;

* flexible CCS
CCSFLEX_POW.lo(i,v,r,h,t_unfix)$[ccsflex(i)$valgen(i,v,r,t_unfix)$(Sw_CCSFLEX_BYP OR Sw_CCSFLEX_STO OR Sw_CCSFLEX_DAC)] = 0 ;
CCSFLEX_POW.up(i,v,r,h,t_unfix)$[ccsflex(i)$valgen(i,v,r,t_unfix)$(Sw_CCSFLEX_BYP OR Sw_CCSFLEX_STO OR Sw_CCSFLEX_DAC)] = +inf ;
CCSFLEX_POWREQ.lo(i,v,r,h,t_unfix)$[ccsflex_sto(i)$valgen(i,v,r,t_unfix)$Sw_CCSFLEX_STO] = 0 ;
CCSFLEX_POWREQ.up(i,v,r,h,t_unfix)$[ccsflex_sto(i)$valgen(i,v,r,t_unfix)$Sw_CCSFLEX_STO] = +inf ;
CCSFLEX_STO_STORAGE_LEVEL.lo(i,v,r,h,t_unfix)$[ccsflex_sto(i)$valgen(i,v,r,t_unfix)$Sw_CCSFLEX_STO] = 0 ;
CCSFLEX_STO_STORAGE_LEVEL.up(i,v,r,h,t_unfix)$[ccsflex_sto(i)$valgen(i,v,r,t_unfix)$Sw_CCSFLEX_STO] = +inf ;

* trade
FLOW.lo(r,rr,h,t_unfix,trtype)$routes(r,rr,trtype,t_unfix) = 0 ;
FLOW.up(r,rr,h,t_unfix,trtype)$routes(r,rr,trtype,t_unfix) = +inf ;
OPRES_FLOW.lo(ortype,r,rr,h,t_unfix)$[Sw_OpRes$opres_model(ortype)$opres_routes(r,rr,t_unfix)$opres_h(h)] = 0 ;
OPRES_FLOW.up(ortype,r,rr,h,t_unfix)$[Sw_OpRes$opres_model(ortype)$opres_routes(r,rr,t_unfix)$opres_h(h)] = +inf ;
PRMTRADE.lo(r,rr,trtype,ccseason,t_unfix)$[routes(r,rr,trtype,t_unfix)$routes_prm(r,rr)] = 0 ;
PRMTRADE.up(r,rr,trtype,ccseason,t_unfix)$[routes(r,rr,trtype,t_unfix)$routes_prm(r,rr)] = +inf ;

* operating reserve
OPRES.lo(ortype,i,v,r,h,t_unfix)$[Sw_OpRes$valgen(i,v,r,t_unfix)$reserve_frac(i,ortype)$opres_h(h)] = 0 ;
OPRES.up(ortype,i,v,r,h,t_unfix)$[Sw_OpRes$valgen(i,v,r,t_unfix)$reserve_frac(i,ortype)$opres_h(h)] = +inf ;

* variable fuel amounts
GASUSED.lo(cendiv,gb,h,t_unfix)$[(Sw_GasCurve=0)$h_rep(h)] = 0 ;
GASUSED.up(cendiv,gb,h,t_unfix)$[(Sw_GasCurve=0)$h_rep(h)] = +inf ;
VGASBINQ_NATIONAL.lo(fuelbin,t_unfix)$[Sw_GasCurve=1] = 0 ;
VGASBINQ_NATIONAL.up(fuelbin,t_unfix)$[Sw_GasCurve=1] = +inf ;
VGASBINQ_REGIONAL.lo(fuelbin,cendiv,t_unfix)$[Sw_GasCurve=1] = 0 ;
VGASBINQ_REGIONAL.up(fuelbin,cendiv,t_unfix)$[Sw_GasCurve=1] = +inf ;
BIOUSED.lo(bioclass,r,t_unfix)$[sum{(i,v)$(bio(i) or cofire(i)), valgen(i,v,r,t_unfix) }] = 0 ;
BIOUSED.up(bioclass,r,t_unfix)$[sum{(i,v)$(bio(i) or cofire(i)), valgen(i,v,r,t_unfix) }] = +inf ;

* RECS
RECS.lo(RPSCat,i,st,ast,t_unfix)$[stfeas(st)$RecMap(i,RPSCat,st,ast,t_unfix)$(stfeas(ast) or sameas(ast,"voluntary"))$Sw_StateRPS] = 0 ;
RECS.up(RPSCat,i,st,ast,t_unfix)$[stfeas(st)$RecMap(i,RPSCat,st,ast,t_unfix)$(stfeas(ast) or sameas(ast,"voluntary"))$Sw_StateRPS] = +inf ;
ACP_PURCHASES.lo(RPSCat,st,t_unfix)$[(stfeas(st) or sameas(st,"voluntary"))$Sw_StateRPS] = 0 ;
ACP_PURCHASES.up(RPSCat,st,t_unfix)$[(stfeas(st) or sameas(st,"voluntary"))$Sw_StateRPS] = +inf ;

* transmission
CONVERSION.lo(r,h,intype,outtype,t_unfix)$Sw_VSC = 0 ;
CONVERSION.up(r,h,intype,outtype,t_unfix)$Sw_VSC = +inf ;
CONVERSION_PRM.lo(r,ccseason,intype,outtype,t_unfix)$Sw_VSC = 0 ;
CONVERSION_PRM.up(r,ccseason,intype,outtype,t_unfix)$Sw_VSC = +inf ;

* water
WAT.lo(i,v,w,r,h,t_unfix)$[i_water(i)$valgen(i,v,r,t_unfix)$Sw_WaterMain] = 0 ;
WAT.up(i,v,w,r,h,t_unfix)$[i_water(i)$valgen(i,v,r,t_unfix)$Sw_WaterMain] = +inf ;

* H2 and DAC production
PRODUCE.lo(p,i,v,r,h,t_unfix)$[consume(i)$i_p(i,p)$valcap(i,v,r,t_unfix)$h_rep(h)$Sw_Prod] = 0 ;
PRODUCE.up(p,i,v,r,h,t_unfix)$[consume(i)$i_p(i,p)$valcap(i,v,r,t_unfix)$h_rep(h)$Sw_Prod] = +inf ;
H2_FLOW.lo(r,rr,h,t_unfix)$[h2_routes(r,rr)$(Sw_H2=2)] = 0 ;
H2_FLOW.up(r,rr,h,t_unfix)$[h2_routes(r,rr)$(Sw_H2=2)] = +inf ;
H2_STOR_IN.lo(h2_stor,r,h,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)] = 0 ;
H2_STOR_IN.up(h2_stor,r,h,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)] = +inf ;
H2_STOR_OUT.lo(h2_stor,r,h,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)] = 0 ;
H2_STOR_OUT.up(h2_stor,r,h,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)] = +inf ;
H2_STOR_LEVEL.lo(h2_stor,r,actualszn,h,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)$(Sw_H2_StorTimestep=2)] = 0 ;
H2_STOR_LEVEL.up(h2_stor,r,actualszn,h,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)$(Sw_H2_StorTimestep=2)] = +inf ;
H2_STOR_LEVEL_SZN.lo(h2_stor,r,actualszn,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)$(Sw_H2_StorTimestep=1)] = 0 ;
H2_STOR_LEVEL_SZN.up(h2_stor,r,actualszn,t_unfix)$[(h2_stor_r(h2_stor,r))$(Sw_H2=2)$(Sw_H2_StorTimestep=1)] = +inf ;
CREDIT_H2PTC.lo(i,v,r,h,t_unfix)$[valgen_h2ptc(i,v,r,t_unfix)$h_rep(h)$Sw_H2_PTC$h2_ptc_years(t_unfix)$(yeart(t_unfix)>=h2_demand_start)] = 0 ;
CREDIT_H2PTC.up(i,v,r,h,t_unfix)$[valgen_h2ptc(i,v,r,t_unfix)$h_rep(h)$Sw_H2_PTC$h2_ptc_years(t_unfix)$(yeart(t_unfix)>=h2_demand_start)] = +inf ;

* CO2 capture and storage
CO2_CAPTURED.lo(r,h,t_unfix)$Sw_CO2_Detail = 0 ;
CO2_CAPTURED.up(r,h,t_unfix)$Sw_CO2_Detail = +inf ;
CO2_STORED.lo(r,cs,h,t_unfix)$[Sw_CO2_Detail$r_cs(r,cs)] = 0 ;
CO2_STORED.up(r,cs,h,t_unfix)$[Sw_CO2_Detail$r_cs(r,cs)] = +inf ;
CO2_FLOW.lo(r,rr,h,t_unfix)$[Sw_CO2_Detail$co2_routes(r,rr)] = 0 ;
CO2_FLOW.up(r,rr,h,t_unfix)$[Sw_CO2_Detail$co2_routes(r,rr)] = +inf ;

* Positive/negative variables
EMIT.lo(etype,e,r,t_unfix)$emit_modeled(e,r,t_unfix) = -inf ;
EMIT.up(etype,e,r,t_unfix)$emit_modeled(e,r,t_unfix) = +inf ;
STORAGE_INTERDAY_DISPATCH.lo(i,v,r,h,t_unfix)$[storage_interday(i)$valgen(i,v,r,t_unfix)] = -inf ;
STORAGE_INTERDAY_DISPATCH.up(i,v,r,h,t_unfix)$[storage_interday(i)$valgen(i,v,r,t_unfix)] = +inf ;
STORAGE_INTERDAY_LEVEL_MAX_DAY.lo(i,v,r,allszn,t_unfix)$[storage_interday(i)$valgen(i,v,r,t_unfix)] = -inf ;
STORAGE_INTERDAY_LEVEL_MAX_DAY.up(i,v,r,allszn,t_unfix)$[storage_interday(i)$valgen(i,v,r,t_unfix)] = +inf ;
STORAGE_INTERDAY_LEVEL_MIN_DAY.lo(i,v,r,allszn,t_unfix)$[storage_interday(i)$valgen(i,v,r,t_unfix)] = -inf ;
STORAGE_INTERDAY_LEVEL_MIN_DAY.up(i,v,r,allszn,t_unfix)$[storage_interday(i)$valgen(i,v,r,t_unfix)] = +inf ;