diff --git a/core/declarations.gms b/core/declarations.gms index 99fa4ac30..556f80571 100644 --- a/core/declarations.gms +++ b/core/declarations.gms @@ -529,6 +529,8 @@ q_PE_histCap(ttot,all_regi,all_enty,all_enty) "require minimum of fossil q_PE_histCap_NGCC_2020_up(ttot,all_regi,all_enty,all_enty) "require maximum of gas power capacity in 2020 based on 2015 historical data and growth rate assumptions" q_shbiofe_up(ttot,all_regi,all_enty,emi_sectors,all_emiMkt) "constrain share of biomass in hydrocarbons of sectoral final energy in historical years based on historical data (upper bound)" q_shbiofe_lo(ttot,all_regi,all_enty,emi_sectors,all_emiMkt) "constrain share of biomass in hydrocarbons of sectoral final energy in historical years based on historical data (lower bound)" +q_shbiofe_indst_up(ttot,all_regi,all_enty,emi_sectors) "constrain share of biomass in solids of industry final energy in historical years based on historical data (upper bound)" +q_shbiofe_indst_lo(ttot,all_regi,all_enty,emi_sectors) "constrain share of biomass in solids of industry final energy in historical years based on historical data (lower bound)" *** energy service layer equations (only relevant for transport) q_transFe2Es(ttot,all_regi,all_enty,all_esty,all_teEs) "conversion from final energy to energy services" diff --git a/core/equations.gms b/core/equations.gms index 9e541d22b..8d3e80033 100644 --- a/core/equations.gms +++ b/core/equations.gms @@ -1140,7 +1140,7 @@ q_limitCapFeH2BI(t,regi,sector)$(SAMEAS(sector,"build") OR SAMEAS(sector,"indst" ***--------------------------------------------------------------------------- *' Enforce historical data biomass share per carrier in sector final energy for transport and buildings (+- 2%) -*' Exempt industry solids as these are covered by the equation q37_limitBioSolidsIndst +*' Exempt industry solids as they are covered below ***--------------------------------------------------------------------------- q_shbiofe_up(t,regi,entyFe,sector,emiMkt)$(pm_secBioShare(t,regi,entyFe,sector) and sector2emiMkt(sector,emiMkt) and NOT (sameas(sector,"indst") and sameas(entyFe,"fesos"))).. @@ -1159,6 +1159,35 @@ q_shbiofe_lo(t,regi,entyFe,sector,emiMkt)$(pm_secBioShare(t,regi,entyFe,sector) sum((entySeBio,te)$se2fe(entySeBio,entyFe,te), vm_demFeSector_afterTax(t,regi,entySeBio,entyFe,sector,emiMkt)) ; +***--------------------------------------------------------------------------- +*' Enforce historical data biomass share per carrier in sector final energy for industry (+- 2%) +*' Applies to the sum of emiMkt +***--------------------------------------------------------------------------- + +q_shbiofe_indst_up(t,regi,entyFe,sector)$(pm_secBioShare(t,regi,entyFe,sector) and (sameas(sector,"indst") and sameas(entyFe,"fesos"))).. + (pm_secBioShare(t,regi,entyFe,sector) + 0.02) + * + sum(emiMkt$sector2emiMkt(sector,emiMkt), + sum((entySe,te)$se2fe(entySe,entyFe,te), + vm_demFeSector_afterTax(t,regi,entySe,entyFe,sector,emiMkt))) + =g= + sum(emiMkt$sector2emiMkt(sector,emiMkt), + sum((entySeBio,te)$se2fe(entySeBio,entyFe,te), + vm_demFeSector_afterTax(t,regi,entySeBio,entyFe,sector,emiMkt))) +; + +q_shbiofe_indst_lo(t,regi,entyFe,sector)$(pm_secBioShare(t,regi,entyFe,sector) and (sameas(sector,"indst") and sameas(entyFe,"fesos"))).. + (pm_secBioShare(t,regi,entyFe,sector) - 0.02) + * + sum(emiMkt$sector2emiMkt(sector,emiMkt), + sum((entySe,te)$se2fe(entySe,entyFe,te), + vm_demFeSector_afterTax(t,regi,entySe,entyFe,sector,emiMkt))) + =l= + sum(emiMkt$sector2emiMkt(sector,emiMkt), + sum((entySeBio,te)$se2fe(entySeBio,entyFe,te), + vm_demFeSector_afterTax(t,regi,entySeBio,entyFe,sector,emiMkt))) +; + ***--------------------------------------------------------------------------- *' Penalty for secondary energy share deviation in sectors ***---------------------------------------------------------------------------