exam 6 sol

Page 1

EXAM N3 Sets plant /'coal','gas','solar',’dac’/ t /2020*2060/; Sets t0(t); #crea un subset del set t t0(t) = yes$(ord(t) eq 1); #riempie il subset: t0 è il 2020 Parameters c_fx(plant) 'unit cost of installing a plant [$/kW]' #costo fisso c_op(plant) 'unit cost of operating a plant [$/kWh]' #costo operativo cf(plant) 'capacity factor of plant [hours equivalent/(hours*yr)]' lifetime(plant) 'lifetime of the plant' delta(plant) 'capital depreciation factor for the plant'; Scalar prtp 'pure rate of time preference' /0.03/; Scalar tstep /1/; #ogni t vale 1 anno ora *assume all relevant data, BUT the new ones you need to introduce, are loaded here $batinclude data.gms Parameters dmd(t) 'average (yearly) energy demand' rr(t) 'discount factor' ; Parameter emi_st(plant); emi_st(‘coal’) = 1; emi_st(‘gas’) = 0.5; emi_st(‘solar’) = 0; Parameter cprice(t); cprice(t)$(t.val ge 2025) = 100*(1+0.02)*(t.val-2025); Variables I(plant,t) investment in plant [kW*yr-1] K(plant,t) capacity of plant [kW] E(plant,t) energy generated by plant [kWh*yr-1] C(plant,t) cost of plant [$*yr-1] OBJ objective variable EMI(plant,t) K_DAC(plant,t); Positive variables E,K,C,I; K.l(plant,t) = 1; K.up(plant,t) = 100; I.up(plant,t) = 10; K.fx("gas",t0) = 0.5; K.fx("coal",t0) = 2; K.fx("solar",t0) = 0; Equations eq_stock eq_en eq_dmd eq_cost eq_obj; eq_emi; eq_rem; eq_stock(plant,t+1).. K(plant,t+1) =E= K(plant,t)*(1-delta(plant))**tstep + I(plant,t)*tstep; #solita equazione tipo del capitale, è quello che mi interessa eq_en(plant,t).. E(plant,t) =L= K(plant,t)*cf(plant)*8760; #vincolo eq_en(‘dac’,t).. E(‘dac’,t) =G= en_sp(‘dac’,t)*E(‘dac’,t); eq_cost(plant,t).. C(plant,t) =E= c_fx(plant)*I(plant,t)*lifetime(plant)/tstep + E(plant,t)*c_op(plant) + EMI(plant,t)*cprice(t) ; # costo = fisso + di operazione…