diff --git a/core/declarations.gms b/core/declarations.gms index 6ffb720daab8ed02a7af8b89eade13bfe02405c5..29febcbcdbf63f8a615dcfc86e7c63a6fe6a92b1 100644 --- a/core/declarations.gms +++ b/core/declarations.gms @@ -432,6 +432,8 @@ s_macChange "maximum yearly increase o s_tgn_2_pgc "conversion factor 100-yr GWP from TgN to PgCeq" s_tgch4_2_pgc "conversion factor 100-yr GWP from TgCH4 to PgCeq" +sm_MtCH4_2_TWa "Energy content of methane. MtCH4 --> TWa: 1 MtCH4 = 1.23 * 10^6 toe * 42 GJ/toe * 10^-9 EJ/GJ * 1 TWa/31.536 EJ = 0.001638 TWa (BP statistical review)" /0.001638/ + sm_D2015_2_D2005 "Convert $2015 to $2005 by dividing by 1.2: 1/1.2 = 0.8333" /0.8333/ sm_DptCO2_2_TDpGtC "Conversion multiplier to go from $/tCO2 to T$/GtC: 44/12/1000" /0.00366667/ diff --git a/core/equations.gms b/core/equations.gms index 4f8c6539f00d30663d71080b7f325840b62f7652..987f1e53d02c1f672bf4f1ff060ec1f5a0777051 100644 --- a/core/equations.gms +++ b/core/equations.gms @@ -120,7 +120,7 @@ q_balSe(t,regi,enty2)$( entySE(enty2) AND (NOT (sameas(enty2,"seel"))) ).. ) *** add (reused gas from waste landfills) to segas to not account for CO2 *** emissions - it comes from biomass - + ( 0.001638 + + ( sm_MtCH4_2_TWa * ( vm_macBase(t,regi,"ch4wstl") - vm_emiMacSector(t,regi,"ch4wstl") ) diff --git a/core/presolve.gms b/core/presolve.gms index 4629f60437462b4a5f7260a4d963197da3571245..d55cb4ace552d5412a9e3a63b1b3fe03ed09918e 100644 --- a/core/presolve.gms +++ b/core/presolve.gms @@ -269,9 +269,8 @@ pm_macStep(ttot,regi,enty)$(MacSector(enty)) = min(801, ceil(p_priceCO2forMAC(ttot,regi,enty) / sm_dmac) + 1); *** If the gas price increase since 2005 is higher than the CO2 price, it will drive CH4 emission abatement. -*** Conversion: pm_priceCO2 [$/tCeq]; T$/TWa = 1e6 M$/TWa * 0.001638 TWa/MtCH4 * 1 MtCH4/s_gwpCH4 MtCO2eq * (44/12) MtCO2eq/MtCeq -*** conversion factor MtCH4 --> TWa: 1 MtCH4 = 1.23 * 10^6 toe * 42 GJ/toe * 10^-9 EJ/GJ * 1 TWa/31.536 EJ = 0.001638 TWa (BP statistical review) -p_priceGas(ttot,regi)=q_balPe.m(ttot,regi,"pegas")/(qm_budget.m(ttot,regi)+sm_eps) * 1000000 * 0.001638 * (1/s_gwpCH4) * 44/12; +*** Conversion: pm_priceCO2 [$/tCeq]; T$/TWa = 1e6 M$/TWa * sm_MtCH4_2_TWa * 1 MtCH4/s_gwpCH4 MtCO2eq * (44/12) MtCO2eq/MtCeq +p_priceGas(ttot,regi)=q_balPe.m(ttot,regi,"pegas")/(qm_budget.m(ttot,regi)+sm_eps) * 1000000 * sm_MtCH4_2_TWa * (1/s_gwpCH4) * 44/12; pm_macStep(ttot,regi,"ch4gas") = min(801, ceil(max(pm_priceCO2(ttot,regi) * (25/s_gwpCH4), max(0,(p_priceGas(ttot,regi)-p_priceGas("2005",regi))) ) / sm_dmac) + 1); @@ -403,9 +402,8 @@ pm_macCost(ttot,regi,emiMacSector(enty)) *JeS* add 50% of abated CH4 from coal MACs to PEprod. CH4 from oil production is usually flared and not re-used. CH4 from gas production is mostly avoided losses from leakages. *** These losses are not accounted for, so neither are the avoided losses. -*** conversion factor MtCH4 --> TWa: 1 MtCH4 = 1.23 * 10^6 toe * 42 GJ/toe * 10^-9 EJ/GJ * 1 TWa/31.536 EJ = 0.001638 (BP statistical review) p_macPE(ttot,regi,enty) = 0.0; -p_macPE(ttot,regi,"pegas")$(ttot.val gt 2005) = 0.001638 * 0.5 * (vm_macBase.l(ttot,regi,"ch4coal")-vm_emiMacSector.l(ttot,regi,"ch4coal")); +p_macPE(ttot,regi,"pegas")$(ttot.val gt 2005) = sm_MtCH4_2_TWa * 0.5 * (vm_macBase.l(ttot,regi,"ch4coal")-vm_emiMacSector.l(ttot,regi,"ch4coal")); ***------------ adjust adjustment costs for advanced vehicles according to CO2 price in the previous time step ---------------------- diff --git a/modules/37_industry/fixed_shares/postsolve.gms b/modules/37_industry/fixed_shares/postsolve.gms index 314d19e086469f72336934b500d1e755bafb1edf..2d7a5f21726ff7857130c70f7d92324252a78da7 100644 --- a/modules/37_industry/fixed_shares/postsolve.gms +++ b/modules/37_industry/fixed_shares/postsolve.gms @@ -30,7 +30,7 @@ loop (enty$( sameas(enty,"co2") OR sameas(enty,"cco2") ), vm_prodSE.l(ttot,regi,entyPE2,entySE,te2) ) !! SE CH4 from waste bypass - + ( 0.001638 + + ( sm_MtCH4_2_TWa * ( vm_macBase.l(ttot,regi,"ch4wstl") - vm_emiMacSector.l(ttot,regi,"ch4wstl") ) diff --git a/modules/37_industry/off/not_used.txt b/modules/37_industry/off/not_used.txt index ffee9f1da0e9c7e5a69bdc7712291b62064c610c..916bc26b0fe672f0ca750efc6d5d0ea10399f630 100644 --- a/modules/37_industry/off/not_used.txt +++ b/modules/37_industry/off/not_used.txt @@ -24,3 +24,4 @@ cm_startyear,input,questionnaire fm_dataemiglob,input,questionnaire pm_delta_kap,input,questionnaire pm_macCostSwitch,input,questionnaire +sm_MtCH4_2_TWa,input,questionnaire diff --git a/modules/37_industry/subsectors/not_used.txt b/modules/37_industry/subsectors/not_used.txt index b1320a459ff89139655770c650f6b77c56eb18a2..653bb331285282f1b8a53320015b578cae78ceaf 100644 --- a/modules/37_industry/subsectors/not_used.txt +++ b/modules/37_industry/subsectors/not_used.txt @@ -28,3 +28,4 @@ cm_CCS_chemicals,input,questionnaire cm_CCS_steel,input,questionnaire fm_dataemiglob,input,questionnaire pm_macCostSwitch,input,questionnaire +sm_MtCH4_2_TWa,input,questionnaire