#--------------------------------------------------------------------------------------------# #THE SCENARIO TREE #--------------------------------------------------------------------------------------------# param years:=5; #Number of years between scenario tree levels param inv_levels:=3; #Number of levels at which investments can be made param blocks:=3; #Number of scenario tree building blocks (Low High CCS) param rate:=0.093; #Discount rate. param lifetime:=30; #Calculation horizon in years = economic lifetime of investments set SCEN:= BAU S1 S2 S3; #Set of all scenarios (or paths through the scenario tree) #Set of all nodes in a scenario set SCEN_NODES["BAU"]:=0 1 3 6 10 14 18; set SCEN_NODES["S1"]:= 0 1 3 7 11 15 19; set SCEN_NODES["S2"]:= 0 1 4 8 12 16 20; set SCEN_NODES["S3"]:= 0 2 5 9 13 17 21; #Defines the building block of each node param: block_no:= 0 1 1 1 2 2 3 1 4 2 5 2 6 1 7 3 8 3 9 3 10 1 11 3 12 3 13 3 14 1 15 3 16 3 17 3 18 1 19 3 20 3 21 3; #The probability of each scenario to occur param scen_prob:= BAU 0.25 S1 0.25 S2 0.25 S3 0.25; #--------------------------------------------------------------------------------------------# #THE MODEL OF THE MILL #--------------------------------------------------------------------------------------------# #Definitions of sets #Technologies using a heat surplus at the mill to enable energy exports #(Back-pressure turbine, condensing turbine, lignin extraction, district heating 100gC, district heating LP steam, district heating 60gC) set USE:=BP CT LIG DH100 DHLP DH60 ; #Steam pressure levels set PRESSURE:=HP MP LP; #Steam saving measures set MEASURE:=Flash HWWS WoodYard Press BlowOut BlowDown ConVap7el ConVap7lig PIvap7el PIvap7lig ConVap RBU Piping; #Technology options causing pressure reductions set REDUC:=BP Valves; #Energy "products" for export set PRODUCTS:=Elec Lignin Heat; #Export options for which a finer time resolution is needed, e.g. all options affected by a varying district heating demand. set VARIED:=CT DH100 DHLP DH60; #--------------------------------------------------------------------------------------------# # Definitions of parameters #...........Economic parameters - Prices, costs, etc..........# #Export product of each "USE" option param product:= BP Elec CT Elec LIG Lignin DH100 Heat DHLP Heat DH60 Heat; #Price of product param price: #[M€/GWh] 1 2 3:= Elec 0.073 0.074 0.067 Lignin 0.023 0.031 0.031 Heat 0.025 0.034 0.034 ; #Operating/running cost param operating_cost: #[M€/GWh] 1 2 3 := BP 0 0 0 CT 0 0 0 LIG 0.0065 0.0065 0.0065 DH100 0 0 0 DHLP 0 0 0 DH60 0.020 0.022 0.025 ; #Investment cost for fixed size/cost measures param: measure_cost := # [M€] Flash 3.5 HWWS 0.6 WoodYard 0 Press 6 BlowOut 0 BlowDown 0 ConVap7el 9.3 ConVap7lig 9.7 PIvap7el 11.7 PIvap7lig 10.3 ConVap 4.4 RBU 29.8 Piping 6.6; #Parameters for calculations of linearized investment cost data for energy export options #Base cost and number of linearization intervals param: base_cost ints:= #[M€] [] BP 0.274 10 CT 0.746 8 LIG 0.882 10 DH100 0.400 1 DHLP 0.570 1 DH60 0.1089 2; #Linearization interval sizes param int_size: #[MW] 1 2 3 4 5 6 7 8 9 10:= BP 0.1 0.1 0.3 0.9 2.3 4.1 6.5 9.6 13.4 12.7 CT 0.2 0.3 0.8 1.4 2.3 3.4 4.8 5.8 . . LIG 1 1 3 5 8 12 17 23 30 25 DH100 16.3 . . . . . . . . . DHLP 125 . . . . . . . . . DH60 1 99 . . . . . . . .; #Parameters a, b and c used to calculate the investment cost according to # cost[M€]= aa * (size[MW] + cc)^bb param: aa bb cc:= BP 1.090 0.6 0 CT 1.960 0.6 0 LIG 0.8817 0.6 0 DH100 0.400 0 0 DHLP 0.030 1 19 DH60 0.1089 1 0; #........Steam balance parameters........# #Steam savings at different steam pressure levels param steam: LP MP HP := # [tonnes/h][tonnes/h][tonnes/h] Flash 3.125 16.125 0 HWWS 12.5 0 0 WoodYard 2.625 0 0 Press 11.75 0 0 BlowOut 13.875 0 0 BlowDown -78.375 -39.5 102.5 ConVap7el 59.5 -19.25 9.25 ConVap7lig 59.5 -19.25 9.25 PIvap7el 63.875 -14.25 9.25 PIvap7lig 63.875 -14.25 9.25 ConVap 0 0 0 RBU 0 0 0 Piping 0 0 0; #States the conversion relation from steam flow to energy export param conv default 0: LP MP HP:= BP . 0.04064 0.0976 #[(MW elec)/(tonnes/h)] CT 0.1243 . . #[(MW elec)/(tonnes/h)] LIG . . 0.8110 #[(MW lignin)/(tonnes/h)] DH100 . . . #[(MW heat)/(tonnes/h)] DHLP 0.647 . . #[(MW heat)/(tonnes/h)] DH60 . . .; #[(MW heat)/(tonnes/h)] #Factor used to include quench water in the steam balances. Quench water is mixed with the steam to saturate it when the steam pressure #is reduced in e.g. a turbine or a pressure-drop valve. param quench: MP LP:= BP 1.040 1.010 #[] Valves 1.174 1.027; #[] #Parameter stating the maximum output which can be achieved without any steam use param max_free default 0:= #[MW] BP 1.3 DH100 16.25 DH60 83.6; #Parameter stating the maximum decrease of steam flow to existing equipment, i.e., the current steam flow. param max_loss (tr) default 0: BP :=# [tonnes/h] HP 201 MP 157 LP 0 ; #The current capacity of existing equipment param existing_size default 0:= BP 24.7; # [MW] #Parameter to convert MW to GWh/year param MW_to_GWh:=7.848; #[GWh/year/MW] #.....Yearly variations parameters........# #Number of periods per year for parameter/variables affected by yearly variations param per:=20; #Length of period [h]. param hours:= 1 720 2 24 3 648 4 24 5 720 6 24 7 720 8 744 9 696 10 24 11 720 12 24 13 720 14 24 15 720 16 744 17 696 18 24 19 720 20 24; #Max output for options affected by yearly variations (typically district heating demand) [MW]. param demand:= 1 20.8971428571429 2 41.04 3 20.8971428571429 4 41.04 5 20.8971428571429 6 41.04 7 10.29 8 3.315 9 3.315 10 2.94 11 1.49 12 1.14 13 3.315 14 2.94 15 10.29 16 10.29 17 20.8971428571429 18 41.04 19 20.8971428571429 20 41.04; #.............CO2 parameters.............# #Reduction of CO2 emissions associated with energy export param co2: #[ktonnes/GWh] 1 2 3:= BP 0.723 0.723 0.136 CT 0.723 0.723 0.136 LIG 0.329 0.329 0.329 DH100 0.387 0.387 0.387 DHLP 0.387 0.387 0.387 DH60 0.129 0.129 0.339; #Requirement on expected "discounted" reduction (relative year 0) during the calculation horizon param co2_limit:=0; # [ktonnes] #CO2 "discount rate" param co2_discount:=0.093; #--------------------------------------------------------------------------------------------# #Sets used in non-general constraints set PIVAPS:=PIvap7el PIvap7lig; set CONVAPS:= ConVap7el ConVap7lig ; set LIGVAPS:= ConVap7lig PIvap7lig; set EVAPS:=PIvap7el PIvap7lig ConVap7el ConVap7lig ConVap; set DH:=DH100 DHLP DH60; set ANTI_PIVAP:=Flash DH100 DH60;