1a_definitions.gms 17.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
20
21
22
23
24
* =============================================================================
* --- Model Related Set Definitions -------------------------------------------
* =============================================================================

Sets

* --- Model Related Selections ------------------------------------------------
25

26
27
28
29
    mType "model types in the Backbone" /
        building,
        invest,
        schedule
30
        /
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53

    mSetting "setting categories for models" /

        // General Time Structure
        t_start, // First time step for the start of simulation
        t_jump, // Number of time steps realized with each solve
        t_horizon, // Length of the simulation horizon (central forecast)
        t_end, // Last time step of the simulation
        intervalEnd, // End index of a time step interval
        intervalLength, // Number of time steps aggregated within interval
        IntervalInHours, // Length of one time step in hours

        // Samples and Forecasts
        samples, // Number of active samples
        forecasts, // Number of active forecasts
        t_forecastStart, // Time step for first reading the forecasts (not necessarily t_start)
        t_forecastLength, // Length of forecasts in time steps
        t_forecastJump, // Number of time steps between each update of the forecasts
        readForecastsInTheLoop, // Flag for updating forecasts in the solve loop

        // Features
        t_reserveLength, // Length of reserve provision horizon in time steps
        t_aggregate // Unit aggregation threshold time index
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
    // Solve info
    solveInfoAttributes "Information about model solves" /
        modelStat
        solveStat
        totalTime
        iterations
        nodes
        numEqu
        numDVar
        numVar
        numNZ
        sumInfes
        objEst
        objVal
        /

    // !!! REDUNDANT SETS PENDING REMOVAL !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    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"
    /

* --- Efficiency Approximation Related Sets -----------------------------------

    // Efficiency Levels and Categories
    effLevel "Pre-defined levels for efficiency representation that can start from t_solve + x"
        / level1*level9 /
Juha Kiviluoma's avatar
Juha Kiviluoma committed
93
    effSelector "Select equations and lambdas/slope for efficiency calculations"
94
        / lambda01*lambda12, directOff, directOnLP, directOnMIP / // NOTE! Lambdas required first!
Juha Kiviluoma's avatar
Juha Kiviluoma committed
95
    effDirect(effSelector) "Using direct input to output equation"
96
        / directOff, directOnLP, directOnMIP /
Juha Kiviluoma's avatar
Juha Kiviluoma committed
97
98
99
    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"
100
        / directOnLP, directOnMIP /
Juha Kiviluoma's avatar
Juha Kiviluoma committed
101
102
    effLambda(effSelector) "Lambdas in use for part-load efficiency representation"
        / lambda01*lambda12 /
103
    effOnline(effSelector) "Efficiency selectors that use online variables"
104
        / directOnLP, directOnMIP, lambda01*lambda12 / // IMPORTANT! Online variables are generated based on this, so keep it up to date!
105

106
107
108
109
110
111
112
113
114
115
116
    // Efficiency Approximation related Counters
    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"
        /eff00*eff12/ // IMPORTANT! Has to equal the same param_unit!
    lambda "Lambda approximation indeces"
        /lambda01*lambda12/ // IMPORTANT! Has to equal effLambda!

* --- General and Directional Sets --------------------------------------------

    // General Counter
117
    counter "General counter set"
118
        /c000*c999/
119
    cc(counter) "Temporary subset of counter used for calculations"
120
121

    // Directional Sets
122
123
124
    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"
125
        / increase, decrease /
126
127
    min_max "Minimum and maximum"
        / min, max /
128

129
130
131
132
133
134
135
* --- Model Feature Sets ------------------------------------------------------

    // Unit Startup Related Sets
    starttype "Startup types" /
        hot "Hot start"
        warm "Warm start"
        cold "Cold start"
136
137
138
        /
    starttypeConstrained(starttype) "Startup types with constrained maximum non-opearational time"
        / hot, warm /
139
140
141
    cost_consumption "Startup cost or startup fuel consumption"
        / cost, consumption /
    unit_capacity "Unit or capacity based parameter"
142
        / unit, capacity / // Is capacity needed anymore?
143

144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
    // Other Features
    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"
        /
; // END Sets

* =============================================================================
* --- Model Parameter Definitions ---------------------------------------------
* =============================================================================

* --- Numeric Model Parameters ------------------------------------------------

// General model parameter arrays
162
Parameter
163
164
165
166
167
*    settings(mSetting)
    mSettings(mType, mSetting) "Model settings array"
    mSettingsEff(mtype, effLevel) "Model efficiency approximation array"
    mInterval(mType, mSetting, counter) "Model interval array"
    t_skip_counter "Numerical counter for solve time steps"
168
169
;

170
// Include additional parameters if found
171
172
173
174
Parameter params(*) /
$if exist 'params.inc' $include 'params.inc'
/;

175
// Activate model features if found
176
Set active(mType, feature) "Set membership tells active model features" /
177
178
179
$if exist 'features.inc' $include 'features.inc'
/;

180
// Parse command line options and store values for features
181
182
183
184
185
186
187
$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%;

188
189
190
* =============================================================================
* --- Parameter Set Definitions -----------------------------------------------
* =============================================================================
191

Juha Kiviluoma's avatar
Juha Kiviluoma committed
192
193
Sets

194
195
* --- Parameter Data Related Sets ---------------------------------------------

196
param_gn  "Possible parameters for grid, node" /
197
    selfDischargeLoss "Self discharge rate of the node [MW/v_state]"
198
    energyStoredPerUnitOfState "A possible unit conversion if v_state uses something else than MWh"
199
200
201
202
203
    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"
204
205
    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"
206
    forecastLength "Length of forecasts in use for the node (hours). After this, the node will use the central forecast."
207
    capacityMargin "Capacity margin used in invest mode (MW)"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
208
209
/

210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
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"
/

Juha Kiviluoma's avatar
Juha Kiviluoma committed
229
230
param_gnn "Set of possible data parameters for grid, node, node (nodal interconnections)" /
    transferCap "Transfer capacity limits"
231
    transferCapBidirectional "Total bidirectional transfer capacity limit"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
232
233
    transferLoss "Transfer losses"
    diffCoeff   "Coefficients for energy diffusion between nodes"
234
    boundStateMaxDiff "Maximum difference of node state pairs"
235
236
237
    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)"
238
    invCost "Investment cost (€/MW)"
239
    annuity "Investment annuity"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
240
241
242
243
244
245
246
/

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)"
247
248
    maxRampUp   "Speed to ramp up (p.u. / min)"
    maxRampDown "Speed to ramp down (p.u. / min)"
249
    rampUpCost  "Wear and tear cost of ramping up (€/MW)"
250
    rampDownCost "Wear and tear cost of ramping down (€/MW)"
251
252
    maxGenCap   "Maximum output capacity investment (MW)"
    maxConsCap  "Maximum loading capacity investment (MW)"
253
254
    minGenCap   "Minimum output capacity investment (MW)"
    minConsCap  "Minimum loading capacity investment (MW)"
255
    upperLimitCapacityRatio  "Ratio of the upper limit of the node state and the unit capacity investment"
256
257
    unitSizeGen "Output capacity of one unit for integer investments (MW)"
    unitSizeCons  "Loading capacity of one unit for integer investments (MW)"
258
259
    unitSizeTot "Sum of output and loading capacity of one unit (MW)"
    unitSizeGenNet "Output minus loading capacity of one unit (MW)"
260
    invCosts    "Investment costs (€/MW)"
261
262
263
    annuity     "Investment annuity factor"
    fomCosts    "Fixed operation and maintenance costs (€/MW/a)"
    inertia     "Inertia of the unit (s)"
264
    unitSizeMVA "Generator MVA rating (MVA)"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
265
/
266

Topi Rasku's avatar
Topi Rasku committed
267
param_unit "Set of possible data parameters for units" /
268
    unitCount   "Number of subunits if aggregated"
269
    outputCapacityTotal "Output capacity of the unit, calculated by summing all the outputs together by default, unless defined in data"
270
    unitOutputCapacityTotal "Output capacity of the unit, calculated by summing all the unit output sizes together by default"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
271
    availability "Availability of given energy conversion technology (p.u.)"
272
    omCosts     "Variable operation and maintenance costs (€/MWh)"
273
    startCostCold "Variable start-up costs for cold starts excluding fuel costs (€/MW)"
274
275
    startCostWarm "Variable start-up costs for warm starts excluding fuel costs (€/MW)"
    startCostHot "Variable start-up costs for hot starts excluding fuel costs (€/MW)"
276
    startFuelConsCold "Consumption of start-up fuel per cold subunit started up (MWh_fuel/MW)"
277
278
    startFuelConsWarm "Consumption of start-up fuel per warm subunit started up (MWh_fuel/MW)"
    startFuelConsHot "Consumption of start-up fuel per hot subunit started up (MWh_fuel/MW)"
279
280
    startColdAfterXhours   "Offline hours after which the start-up will be a cold start (h)"
    startWarmAfterXhours   "Offline hours after which the start-up will be a warm start (h)"
281
    rampSpeedToMinLoad "Ramping speed from start-up to minimum load (p.u. / min.)"
282
283
    minOperationHours "Minimum operation time (h), prevents shutdown after startup until the defined amount of time has passed"
    minShutdownHours "Minimum shut down time (h), prevents starting up again after the defined amount of time has passed"
Juha Kiviluoma's avatar
Juha Kiviluoma committed
284
285
286
287
288
    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"
289
290
    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
291
    level1 * level9 "Level of simplification in the part-load efficiency representation"
292
    useTimeseries "Uses time series form input for unit parameters whenever possible"
293
    section     "Possibility to define a no load fuel use for units with zero minimum output"
294
295
296
    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
297
/
298

Juha Kiviluoma's avatar
Juha Kiviluoma committed
299
300
301
302
303
304
305
306
307
308
309
310
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" /
311
    emissionTax "Emission tax (EUR/tonne)"
312
313
    update_frequency "Frequency of updating reserve contributions"
    gate_closure "Number of timesteps ahead of dispatch that reserves are fixed"
314
    use_time_series "Flag for using time series data. !!! REDUNDANT WITH useTimeseries, PENDING REMOVAL !!!"
315
    reserveContribution "Reliability parameter of reserve provisions"
316
    emissionCap "Emission limit (tonne)"
317
318
319
    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"
320
    kineticEnergyMin "Minimum system kinetic energy (MWs)"
321
322
    constrainedCapMultiplier "Multiplier a(i) for unit investments in equation Sum(i, a(i)*v_invest(i)) <= b"
    constrainedCapTotalMax "Total maximum b for unit investments in equation Sum(i, a(i)*v_invest(i)) <= b"
323
324
    constrainedOnlineMultiplier "Multiplier a(i) for online units in equation Sum(i, a(i)*v_online(i)) <= b"
    constrainedOnlineTotalMax "Total maximum b for online units in equation Sum(i, a(i)*v_online(i)) <= b"
325
326
/

327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
* --- Counters and Directional Sets -------------------------------------------

// Slack categories
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 /

// Flags for boundaries
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 /
; // END parameter set declarations