1a_definitions.gms 15.6 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
$ontext
This file is part of Backbone.

Backbone is free software: you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

Backbone is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public License
along with Backbone.  If not, see <http://www.gnu.org/licenses/>.
$offtext

18
19
* --- Model parameters, features and switches ---------------------------------
Sets  // Model related selections
20
    mType "model types in the Backbone"
21
        /invest, storage, schedule, realtime, building, schedule_dispatch/
22
    mSetting "setting categories for models"
23
        /t_start, t_jump, t_horizon, t_forecastStart, t_forecastLength, t_forecastJump, t_end, t_reserveLength, samples, forecasts, readForecastsInTheLoop, intervalEnd, intervalLength, IntervalInHours, t_aggregate/
24
25
26
    counter "general counter set"
        /c000*c999/

27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
    solveInfoAttributes
        /
          modelStat
          solveStat
          totalTime
          iterations
          nodes
          numEqu
          numDVar
          numVar
          numNZ
          sumInfes
          objEst
          objVal
        /

43
    // Efficiency approximation related sets
44
45
46
    op "Operating points in the efficiency curves, also functions as index for data points"
        /op00*op12/ // IMPORTANT! Has to equal the same param_unit!
    eff "Effiency for the corresponding operating point ('op') in the efficiency curves, also used for data indexing"
47
48
49
        /eff00*eff12/ // IMPORTANT! Has to equal the same param_unit!
    lambda "Lambda approximation indeces"
        /lambda01*lambda12/ // IMPORTANT! Has to equal effLambda!
Juha Kiviluoma's avatar
Juha Kiviluoma committed
50
    effSelector "Select equations and lambdas/slope for efficiency calculations"
51
        / directOff, directOn, lambda01*lambda12 /
Juha Kiviluoma's avatar
Juha Kiviluoma committed
52
53
54
55
56
57
58
59
    effDirect(effSelector) "Using direct input to output equation"
        / directOff, directOn /
    effDirectOff(effSelector) "Using direct input to output equation without online variable, i.e. constant efficiency"
        / directOff /
    effDirectOn(effSelector) "Using direct input to output equation with online variable"
        / directOn /
    effLambda(effSelector) "Lambdas in use for part-load efficiency representation"
        / lambda01*lambda12 /
60
61
    effOnline(effSelector) "Efficiency selectors that use online variables"
        / directOn, lambda01*lambda12 / // IMPORTANT! Online variables are generated based on this, so keep it up to date!
Juha Kiviluoma's avatar
Juha Kiviluoma committed
62
63
    effLevel "Pre-defined levels for efficiency representation that can start from t_solve + x"
        / level1*level9 /
64
65
66
67
68
69

    // Variable direction related sets
    up_down "Direction set used by some variables, e.g. reserve provisions and generation ramps"
        / up, down /
    inc_dec "Increase or decrease in dummy, or slack variables"
       / increase, decrease /
70
71
    min_max "Minimum and maximum"
        / min, max /
72
73
;

74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
Sets //Reserve type sets
    restype "Reserve types"
        / primary "Automatic frequency containment reserves"
          secondary "Fast frequency restoration reserves"
          tertiary "Replacement reserves"
        /
    restypeDirection(restype, up_down) "Different combinations of reserve types and directions"
        / primary.up
          primary.down
          secondary.up
          secondary.down
          tertiary.up
          tertiary.down
        /
;
Juha Kiviluoma's avatar
Juha Kiviluoma committed
89

90
Sets //Startup related sets
91
92
93
94
95
96
97
    starttype "Startup types"
        / hot "Hot start"
          warm "Warm start"
          cold "Cold start"
        /
    starttypeConstrained(starttype) "Startup types with constrained maximum non-opearational time"
        / hot, warm /
98
99
100
101
    cost_consumption "Startup cost or startup fuel consumption"
        / cost, consumption /
    unit_capacity "Unit or capacity based parameter"
        / unit, capacity /
102
103
;

104
105
* Numeric parameters
Parameter
106
    settings(mSetting)
107
    mSettings(mType, mSetting)
Juha Kiviluoma's avatar
Juha Kiviluoma committed
108
    mSettingsEff(mtype, effLevel)
109
    mInterval(mType, mSetting, counter)
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
    t_skip_counter
;


Parameter params(*) /
$if exist 'params.inc' $include 'params.inc'
/;


* Model features
Set feature "Set of optional model features" /
    findStorageStart "Solve for optimal storage start levels"
    storageValue     "Use storage value instead of fixed control"
    storageEnd       "Expected storage end levels greater than starting levels"
    addOn            "Use StoSSch as a storage add-on to a larger model"
    extraRes         "Use extra tertiary reserves for error in elec. load during time step"
    rampSched        "Use power based scheduling"
/;

Set active(feature) "Set membership tells active model features" /
$if exist 'features.inc' $include 'features.inc'
/;

* --- Parse command line options and store values -----------------------------

* Features
$if set findStorageStart active('findStorageStart') = %findStorageStart%;
$if set storageValue active('storageValue') = %storageValue%;
$if set storageEnd active('storageEnd') = %storageEnd%;
$if set addOn active('addOn') = %addOn%;
$if set extraRes active('extraRes') = %extraRes%;
$if set rampSched active('rampSched') = %rampSched%;


Juha Kiviluoma's avatar
Juha Kiviluoma committed
144
145
146
* --- Set definitions for parameters -----------------------------------------------------
Sets

147
param_gn  "Possible parameters for grid, node" /
148
    selfDischargeLoss "Self discharge rate of the node [MW/v_state]"
149
    energyStoredPerUnitOfState "A possible unit conversion if v_state uses something else than MWh"
150
151
152
153
154
    boundStart    "A flag to bound the first t in the run using reference constant or time series"
    boundStartAndEnd "A flag that both start and end are bound using reference constant or time series"
    boundEnd      "A flag to bound last t in each solve based on the reference constant or time series"
    boundAll      "A flag to bound the state to the reference in all time steps"
    boundStartToEnd  "Force the last states to equal the first state"
155
156
    boundCyclic   "A flag to impose cyclic bounds for the first and the last states within a sample"
    boundCyclicBetweenSamples   "A flag to impose cyclic bounds for the last and first states between samples"
157
    forecastLength "Length of forecasts in use for the node (hours). After this, the node will use the central forecast."
158
    capacityMargin "Capacity margin used in invest mode (MW)"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
159
160
/

161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
param_gnBoundaryTypes "Types of boundaries that can be set for a node with a state variable" /
    upwardLimit   "Absolute maximum state of the node (unit depends on energyCapacity)"
    downwardLimit "Absolute minimum energy in the node (unit depends on energyCapacity)"
    upwardSlack01*upwardSlack20 "A threshold after which a specific cost co-efficient is applied (unit depends on energyCapacity)"
    downwardSlack01*downwardSlack20 "A threshold after which a specific cost co-efficient is applied (unit depends on energyCapacity)"
    reference     "Reference value for a state that can be used to bound a state (unit depends on energyCapacity)"
    maxSpill      "Maximum spill rate from the node (MWh/h)"
    minSpill      "Minimum spill rate from the node (MWh/h)"
/

param_gnBoundaryProperties "Properties that can be set for the different boundaries" /
    useTimeSeries "A flag to use time series to set state bounds and limits"
    useConstant   "A flag to use constant to set state bounds and limits"
    deltaFromReference "The constant or the time series indicate how much the boundary deviates from reference (instead of being an absolute number)"
    constant      "A constant value for the boundary or the reference"
    slackCost     "The cost of exceeding the slack boundary"
    multiplier    "A multiplier to change the value (either constant or time series), default 1"
/

slack(param_gnBoundaryTypes) "Categories for slack variables"
       / upwardSlack01*upwardSlack20, downwardSlack01*downwardSlack20 /
upwardSlack(slack) "Set of upward slacks"
       / upwardSlack01*upwardSlack20 /
downwardSlack(slack) "Set of downward slacks"
       / downwardSlack01*downwardSlack20 /
stateLimits(param_gnBoundaryTypes) "set of upward and downward state limits"
       / upwardLimit, downwardLimit /
spillLimits(param_gnBoundaryTypes) "set of upward and downward state limits"
       / maxSpill, minSpill /
useConstantOrTimeSeries(param_gnBoundaryProperties) "useTimeSeries and useConstant property together"
       / useTimeSeries, useConstant /

Juha Kiviluoma's avatar
Juha Kiviluoma committed
193
194
param_gnn "Set of possible data parameters for grid, node, node (nodal interconnections)" /
    transferCap "Transfer capacity limits"
195
    transferCapBidirectional "Total bidirectional transfer capacity limit"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
196
197
198
    transferLoss "Transfer losses"
    diffCoeff   "Coefficients for energy diffusion between nodes"
    boundStateOffset "Offset parameter for relatively bound node states"
199
    boundStateMaxDiff "Maximum difference of node state pairs"
200
201
202
    transferCapInvLimit "Capacity limit for investments (MW)"
    investMIP   "Choice of making integer investment instead of continous investment (MW versus number of links)"
    unitSize    "Size of one link for integer investments (MW)"
203
    invCost "Investment cost (€/MW)"
204
    annuity "Investment annuity"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
205
206
207
208
209
210
211
/

param_gnu "Set of possible data parameters for grid, node, unit" /
    maxGen      "Maximum output capacity (MW)"
    maxCons     "Maximum loading capacity (MW)"
    cB          "Share of output"
    cV          "Reduction in primary output when increasing secondary output, e.g. reduction of electricity generation due to heat generation in extraction CHP (MWh_e/MWh_h)"
212
213
    maxRampUp   "Speed to ramp up (p.u. / min)"
    maxRampDown "Speed to ramp down (p.u. / min)"
214
    rampUpCost  "Wear and tear cost of ramping up (€/MW)"
215
    rampDownCost "Wear and tear cost of ramping down (€/MW)"
216
217
    maxGenCap   "Maximum output capacity investment (MW)"
    maxConsCap  "Maximum loading capacity investment (MW)"
218
219
    minGenCap   "Minimum output capacity investment (MW)"
    minConsCap  "Minimum loading capacity investment (MW)"
220
    upperLimitCapacityRatio  "Ratio of the upper limit of the node state and the unit capacity investment"
221
222
    unitSizeGen "Output capacity of one unit for integer investments (MW)"
    unitSizeCons  "Loading capacity of one unit for integer investments (MW)"
223
224
    unitSizeTot "Sum of output and loading capacity of one unit (MW)"
    unitSizeGenNet "Output minus loading capacity of one unit (MW)"
225
226
    invCosts    "Investment costs (€/MW)"
    annuity     "Investment annuity"
227
228
    inertia  "Inertia of the unit (s)"
    unitSizeMVA "Generator MVA rating (MVA)"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
229
/
230

Topi Rasku's avatar
Topi Rasku committed
231
param_unit "Set of possible data parameters for units" /
232
    unitCount   "Number of subunits if aggregated"
233
    outputCapacityTotal "Output capacity of the unit, calculated by summing all the outputs together by default, unless defined in data"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
234
    availability "Availability of given energy conversion technology (p.u.)"
235
    omCosts     "Variable operation and maintenance costs (€/MWh)"
236
237
238
239
240
241
242
243
244
245
246
247
    startCostCold "Variable start-up costs for cold starts excluding fuel costs (€/subunit) NOT IN USE"
    startCostWarm "Variable start-up costs for warm starts excluding fuel costs (€/subunit) NOT IN USE"
    startCost "Variable start-up costs for starts excluding fuel costs (€/subunit)"
    startFuelConsCold "Consumption of start-up fuel per subunit started up (MWh_fuel/subunit) NOT IN USE"
    startFuelConsWarm "Consumption of start-up fuel per subunit started up (MWh_fuel/subunit) NOT IN USE"
    startFuelCons "Consumption of start-up fuel per subunit started up (MWh_fuel/subunit)"
    startCostCold_MW "Variable start-up costs for cold starts excluding fuel costs (€/MW) NOT IN USE"
    startCostWarm_MW "Variable start-up costs for warm starts excluding fuel costs (€/MW) NOT IN USE"
    startCost_MW "Variable start-up costs for starts excluding fuel costs (€/MW)"
    startFuelConsCold_MW "Consumption of start-up fuel per capacity started up (MWh_fuel/MW) NOT IN USE"
    startFuelConsWarm_MW "Consumption of start-up fuel per capacity started up (MWh_fuel/MW) NOT IN USE"
    startFuelCons_MW "Consumption of start-up fuel per capacity started up (MWh_fuel/MW)"
248
249
    startCold   "Offline hours after which the start-up will be a cold start (h) NOT IN USE"
    startWarm   "Offline hours after which the start-up will be a warm start (h) NOT IN USE"
250
251
    minOperationTime "Minimum operation time (h), prevents shutdown after startup until the defined amount of time has passed"
    minShutDownTime "Minimum shut down time (h), prevents starting up again after the defined amount of time has passed"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
252
253
254
255
256
    SO2         "SO2 emissions (tonne per MWh_fuel)"
    NOx         "NOx emissions (tonne per MWh_fuel)"
    CH4         "CH4 emissions (tonne per MWh_fuel)"
    resTimelim  "How long should a storage be able to provide reserve (h)"
    eff00 * eff12    "Efficiency of the unit to convert input to output/intermediate product"
257
258
    opFirstCross "The operating point where the real efficiency curve and approximated efficiency curve cross"
    op00 * op12     "Right border of the efficiency point"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
259
    level1 * level9 "Level of simplification in the part-load efficiency representation"
260
    useTimeseries "Uses time series form input for unit parameters whenever possible"
261
    section     "Possibility to define a no load fuel use for units with zero minimum output"
262
263
264
    investMIP   "Choice of making integer investment instead of continous investment (number of units versus MW)"
    maxUnitCount "Maximum number of units when making integer investments"
    minUnitCount "Minimum number of units when making integer investments"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
265
/
266

Juha Kiviluoma's avatar
Juha Kiviluoma committed
267
268
269
270
271
272
273
274
275
276
277
278
param_fuel "Parameters for fuels" /
    emissionIntensity "Intensity of emission from fuel (kg/MWh_fuel)"
    main        "Main fuel of the unit - unless input fuels defined as grids"
    startup     "Startup fuel of the unit, if exists. Can be the same as main fuel - consumption using startupFuelCons"
/

param_unitFuel "Parameters for fuel limits in units" /
    maxFuelCons "Maximum absolute fuel consumption in a unit"
    maxFuelFraction "Maximum share of a fuel in the consumption mix"
/

param_policy "Set of possible data parameters for grid, node, regulation" /
279
    emissionTax "Emission tax (EUR/tonne)"
280
281
    update_frequency "Frequency of updating reserve contributions"
    gate_closure "Number of timesteps ahead of dispatch that reserves are fixed"
282
    use_time_series "Flag for using time series data. !!! REDUNDANT WITH useTimeseries, PENDING REMOVAL !!!"
283
    reserveContribution "Reliability parameter of reserve provisions"
284
    emissionCap "Emission limit (tonne)"
285
286
287
    instantaneousShareMax "Maximum instantaneous share of generation and import from a particular group of units and transfer links"
    energyShareMax "Maximum energy share of generation from a particular group of units"
    energyShareMin "Minimum energy share of generation from a particular group of units"
288
    kineticEnergyMin "Minimum system kinetic energy (MWs)"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
289
290
/

291
292
293
294
param_gnugnu "Set of possible data parameters for grid, node, unit, grid, node, unit" /
    capacityRatio "Fixed ratio of the capacity of the first unit to the second"
/

295
296
297
298
299
300
301
302
303
304
305
306
307
308
solve_info "Containers for solve information" /
    modelStat "Status of the model after solve"
    solveStat "Status of the solve"
    totalTime "Total solve time"
    iterations "Number of iteration"
    nodes "Number of nodes in the solve"
    numEqu "Number of equations in the problem"
    numDVar "Number of D variables in the problem"
    numVar "Number of variables in the problem"
    numNZ "Number of non-zeros in the problem"
    sumInfes "Sum of infeasibilities"
    objEst "Estimate for the best possible objective value"
    objVal "Objectiv value"
/
309
; // End parameter set declarations