Dev/feedstocks review

config/default.cfg

@@ -27,7 +27,7 @@ cfg$regionmapping <- "config/regionmappingH12.csv"
 ### Additional (optional) region mapping, so that those validation data can be loaded that contain the corresponding additional regions.
 cfg$extramappings_historic <- ""
 #### Current input data revision (<mainrevision>.<subrevision>) ####
-cfg$inputRevision <- "6.501"
+cfg$inputRevision <- "6.54APT-2023-06-11"
 
 #### Current CES parameter and GDX revision (commit hash) ####
 cfg$CESandGDXversion <- "7c6249fa1d775d0e1a4ce1602fc321b47c2a33be"

core/datainput.gms

@@ -1461,6 +1461,40 @@ pm_emifac(ttot,regi,"seliqfos","fehos","tdfoshos","co2") = p_ef_dem(regi,"fehos"
 pm_emifac(ttot,regi,"seliqfos","fepet","tdfospet","co2") = p_ef_dem(regi,"fepet") / (sm_c_2_co2*1000*sm_EJ_2_TWa); !! GtC/TWa
 pm_emifac(ttot,regi,"seliqfos","fedie","tdfosdie","co2") = p_ef_dem(regi,"fedie") / (sm_c_2_co2*1000*sm_EJ_2_TWa); !! GtC/TWa
 
+***------ Read in emission factors for process emissions in chemicals sector---
+*' calculated using IEA data on feedstocks flows and UNFCCC data on chem sector process emissions
+***TO DO: remove the "*" by modif the input file (fegas instead of gases, etc.)
+*** "non-energy emission factors [GtC per ZJ]"
+parameter f_nechem_emissionFactors(ttot,all_regi,*)  
+/
+$ondelim
+$include "./core/input/f_nechem_emissionFactors.cs4r"
+$offdelim
+/
+;
+
+pm_emifacNonEnergy(ttot,regi,'segafos','fegas','indst','co2') = f_nechem_emissionFactors(ttot,regi,"gases") / s_zj_2_twa;
+pm_emifacNonEnergy(ttot,regi,'seliqfos','fehos','indst','co2') = f_nechem_emissionFactors(ttot,regi,"liquids") / s_zj_2_twa;
+pm_emifacNonEnergy(ttot,regi,'sesofos','fesos','indst','co2') = f_nechem_emissionFactors(ttot,regi,"solids") / s_zj_2_twa;
+
+pm_emifacNonEnergy(ttot,regi,'segabio','fegas','indst','co2') = f_nechem_emissionFactors(ttot,regi,"gases") / s_zj_2_twa;
+pm_emifacNonEnergy(ttot,regi,'seliqbio','fehos','indst','co2') = f_nechem_emissionFactors(ttot,regi,"liquids") / s_zj_2_twa;
+pm_emifacNonEnergy(ttot,regi,'sesobio','fesos','indst','co2') = f_nechem_emissionFactors(ttot,regi,"solids") / s_zj_2_twa;
+
+pm_emifacNonEnergy(ttot,regi,'segasyn','fegas','indst','co2') = f_nechem_emissionFactors(ttot,regi,"gases") / s_zj_2_twa;
+pm_emifacNonEnergy(ttot,regi,'seliqsyn','fehos','indst','co2') = f_nechem_emissionFactors(ttot,regi,"liquids") / s_zj_2_twa;
+
+***------ Read in projections for incineration rates of plastic waste---
+*** "incineration rates [fraction]"
+parameter f_incinerationShares(ttot,all_regi)         "incineration rate of plastic waste"
+/
+$ondelim
+$include "./core/input/f_incinerationShares.cs4r"
+$offdelim
+/
+;
+pm_incinerationRate(ttot,all_regi)=f_incinerationShares(ttot,all_regi);
+
 *** some balances are not matching by small amounts;
 *** the differences are cancelled out here!!!
 pm_cesdata(ttot,regi,in,"offset_quantity")$(ttot.val ge 2005)       = 0;

core/declarations.gms

@@ -122,11 +122,16 @@ p_extRegiccsinjecrateRegi(ext_regi)                         "Regional CCS inject
 pm_dataeta(tall,all_regi,all_te)                            "regional eta data"
 p_emi_quan_conv_ar4(all_enty)                               "conversion factor for various gases to GtCeq"
 pm_emifac(tall,all_regi,all_enty,all_enty,all_te,all_enty)  "emission factor by technology for all types of emissions in emiTe"
+pm_emifacNonEnergy(ttot,all_regi,all_enty,all_enty,emi_sectors,all_enty)                "emission factor for non-energy fedstocks. For now only for Chemicals Industry [GtC per TWa]"
+pm_incinerationRate(ttot,all_regi)                          "share of plastic waste that gets incinerated [fraction]"
 pm_omeg (all_regi,opTimeYr,all_te)                          "technical depreciation parameter, gives the share of a capacity that is still usable after tlt. [none/share, value between 0 and 1]"
 p_aux_lifetime(all_regi,all_te)                             "auxiliary parameter for calculating life times, calculated externally in excel sheet"
 pm_pedem_res(ttot,all_regi,all_te)                          "Demand for pebiolc residues, needed for enhancement of residue potential [TWa]"
 p_ef_dem(all_regi,all_enty)                                 "Demand side emission factor of final energy carriers [MtCO2/EJ]"
 
+*** share of plastic waste that is incinerated REDUNDANT; DELETE ME
+pm_incinerationShares(ttot,all_regi)                        "Share of plastic waste that gets incinerated [fraction]"
+
 pm_secBioShare(ttot,all_regi,all_enty,emi_sectors)           "share of biomass per carrier for each sector"
 
 p_avCapFac2015(all_regi,all_te)                             "average capacity factor of non-bio renewables in 2015 in REMIND"
@@ -389,6 +394,7 @@ vm_prodPe(ttot,all_regi,all_enty)                    "pe production. [TWa, Urani
 vm_demSe(ttot,all_regi,all_enty,all_enty,all_te)     "se demand. [TWa]"
 vm_prodSe(tall,all_regi,all_enty,all_enty,all_te)    "se production. [TWa]"
 vm_prodFe(ttot,all_regi,all_enty,all_enty,all_te)    "fe production. [TWa]"
+vm_demFENonEnergySector(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt) "energy flows of non-energy feedstocks [TWa]"
 vm_demFeSector(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt)          "fe demand per sector and emission market. Taxes should be applied to this variable or variables closer to the supply side whenever possible so the marginal prices include the tax effects. [TWa]"
 vm_demFeSector_afterTax(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt) "fe demand per sector and emission market after tax. Demand sectors should use this variable in their fe balance equations so demand side marginals include taxes effects. [TWa]"
 v_costFu(ttot,all_regi)                              "fuel costs"
@@ -421,6 +427,12 @@ v_shGasLiq_fe(ttot,all_regi,emi_sectors)             "share of gases and liquids
 
 vm_emiCdrAll(ttot,all_regi)                          "all CDR emissions"
 
+vm_FeedstocksCarbon(ttot,all_regi,all_enty,all_enty,all_emiMkt)             "carbon contained in chemical feedstocks [GtC]"
+vm_plasticsCarbon(ttot,all_regi,all_enty,all_enty,all_emiMkt)               "carbon contained in plastics [GtC]"
+vm_feedstockEmiUnknownFate(ttot,all_regi,all_enty,all_enty,all_emiMkt)      "carbon contained in plastics with unknown fate (assumed to go back into the atmosphere) [GtC]"
+vm_incinerationEmi(ttot,all_regi,all_enty,all_enty,all_emiMkt)              "Emissions from incineration of plastic waste [GtC]"
+vm_nonIncineratedPlastics(ttot,all_regi,all_enty,all_enty,all_emiMkt)       "carbon contained in plastics that are not incinerated [GtC]"
+
 v_changeProdStartyearAdj(ttot,all_regi,all_te)       "Absolute effect size of changing output with respect to the reference run for each te"
 vm_changeProdStartyearCost(ttot,all_regi,all_te)     "Costs for changing output with respect to the reference run for each te"
 

core/equations.gms

@@ -511,26 +511,33 @@ q_emiTe(t,regi,emiTe(enty))..
 *' from secondary to final energy transformation (some air pollutants), or
 *' transformations within the chain of CCS steps (Leakage).
 ***-----------------------------------------------------------------------------
-
-q_emiTeDetailMkt(t,regi,enty,enty2,te,enty3,emiMkt)$(emi2te(enty,enty2,te,enty3) OR (pe2se(enty,enty2,te) AND sameas(enty3,"cco2")) ) ..
+q_emiTeDetailMkt(t,regi,enty,enty2,te,enty3,emiMkt)$(
+                           emi2te(enty,enty2,te,enty3) 
+                        OR (pe2se(enty,enty2,te) AND sameas(enty3,"cco2")) ) ..
   vm_emiTeDetailMkt(t,regi,enty,enty2,te,enty3,emiMkt)
   =e=
     sum(emi2te(enty,enty2,te,enty3),
-      (
-	    sum(pe2se(enty,enty2,te),
-		  pm_emifac(t,regi,enty,enty2,te,enty3)
-		  * vm_demPE(t,regi,enty,enty2,te)
-		  )
-	    + sum((ccs2Leak(enty,enty2,te,enty3),teCCS2rlf(te,rlf)),
-		    pm_emifac(t,regi,enty,enty2,te,enty3)
-		    * vm_co2CCS(t,regi,enty,enty2,te,rlf)
-		  )
-	  )$(sameas(emiMkt,"ETS"))
-	  + sum(se2fe(enty,enty2,te),
+      ( sum(pe2se(enty,enty2,te),
           pm_emifac(t,regi,enty,enty2,te,enty3)
-		  * sum(sector$(entyFe2Sector(enty2,sector) AND sector2emiMkt(sector,emiMkt)), vm_demFeSector(t,regi,enty,enty2,sector,emiMkt))
-		)
-	)
+        * vm_demPE(t,regi,enty,enty2,te)
+      )
+    + sum((ccs2Leak(enty,enty2,te,enty3),teCCS2rlf(te,rlf)),
+        pm_emifac(t,regi,enty,enty2,te,enty3)
+      * vm_co2CCS(t,regi,enty,enty2,te,rlf)
+      )
+    )$( sameas(emiMkt,"ETS") )
+  + sum(se2fe(enty,enty2,te),
+      pm_emifac(t,regi,enty,enty2,te,enty3)
+    * sum(sector$(    entyFe2Sector(enty2,sector) 
+                  AND sector2emiMkt(sector,emiMkt) ), 
+        vm_demFeSector(t,regi,enty,enty2,sector,emiMkt)
+        !! substract FE used for non-energy purposes (as feedstocks) so it does
+        !! not create energy-related emissions
+      - sum(entyFe2sector2emiMkt_NonEn(enty2,sector,emiMkt),
+          vm_demFENonEnergySector(t,regi,enty,enty2,sector,emiMkt))
+        )
+      )
+    )
 ;
 
 ***--------------------------------------------------
@@ -558,44 +565,76 @@ q_emiEnFuelEx(t,regi,emiTe(enty))..
 ***--------------------------------------------------
 *' Total energy-emissions per emission market, region and timestep  
 ***--------------------------------------------------
-q_emiTeMkt(t,regi,emiTe(enty),emiMkt)..
+q_emiTeMkt(t,regi,emiTe(enty),emiMkt) ..
   vm_emiTeMkt(t,regi,enty,emiMkt)
   =e=
-***   emissions from fuel combustion
+    !! emissions from fuel combustion
     sum(emi2te(enty2,enty3,te,enty),     
       vm_emiTeDetailMkt(t,regi,enty2,enty3,te,enty,emiMkt)
     )
-***   energy emissions fuel extraction
-	+ v_emiEnFuelEx(t,regi,enty)$(sameas(emiMkt,"ETS"))
-***   Industry CCS emissions
-	- ( sum(emiMac2mac(emiInd37_fuel,enty2),
-		  vm_emiIndCCS(t,regi,emiInd37_fuel)
-		)$( sameas(enty,"co2") )
-	)$(sameas(emiMkt,"ETS"))
-***   LP, Valve from cco2 capture step, to mangage if capture capacity and CCU/CCS capacity don't have the same lifetime
-  + ( v_co2capturevalve(t,regi)$( sameas(enty,"co2") ) )$(sameas(emiMkt,"ETS"))
-***  JS CO2 from short-term CCU (short term CCU co2 is emitted again in a time period shorter than 5 years)
+    !! energy emissions fuel extraction
+  + v_emiEnFuelEx(t,regi,enty)$(sameas(emiMkt,"ETS"))
+    !! Industry CCS emissions
+  - ( sum(emiMac2mac(emiInd37_fuel,enty2),
+        vm_emiIndCCS(t,regi,emiInd37_fuel)
+      )$( sameas(enty,"co2") )
+    )$(sameas(emiMkt,"ETS") )
+    !! substract carbon from biogenic or synthetic origin contained in
+    !! plastics that don't get incinerated ("plastic removals")
+  - sum(entyFe2sector2emiMkt_NonEn(entyFe,"indst",emiMkt),
+      sum(se2fe(entySe,entyFe,te)$( entySeBio(entySe) OR entySeSyn(entySe) ),
+        vm_nonIncineratedPlastics(t,regi,entySe,entyFe,emiMkt)
+      )
+    )$( sameas(enty,"co2") )
+    !! add emissions from plastics incineration. CHECK FOR DOUBLE-COUNTING RISK
+  + sum(entyFe2sector2emiMkt_NonEn(entyFe,"indst",emiMkt),
+      sum(se2fe(entySe,entyFe,te),
+        vm_incinerationEmi(t,regi,entySe,entyFe,emiMkt)
+      )
+    )$( sameas(enty,"co2") )
+    !! add emissions from chemical feedstock with unknown fate
+  + sum(entyFe2sector2emiMkt_NonEn(entyFe,"indst",emiMkt),
+      sum(se2fe(entySe,entyFe,te),
+        vm_feedstockEmiUnknownFate(t,regi,entySe,entyFe,emiMkt)
+      )
+    )$( sameas(enty,"co2") )
+    !! Valve from cco2 capture step, to mangage if capture capacity and CCU/CCS
+    !! capacity don't have the same lifetime
+  + v_co2capturevalve(t,regi)$( sameas(enty,"co2") AND sameas(emiMkt,"ETS") )
+    !! CO2 from short-term CCU (short term CCU co2 is emitted again in a time
+    !! period shorter than 5 years)
   + sum(teCCU2rlf(te2,rlf),
-		vm_co2CCUshort(t,regi,"cco2","ccuco2short",te2,rlf)$( sameas(enty,"co2") ) 
-	)$(sameas(emiMkt,"ETS"))
+      vm_co2CCUshort(t,regi,"cco2","ccuco2short",te2,rlf)$( sameas(enty,"co2") ) 
+    )$(sameas(emiMkt,"ETS"))
 ;
 
 ***--------------------------------------------------
 *' Total emissions
 ***--------------------------------------------------
-q_emiAllMkt(t,regi,emi,emiMkt)..
+q_emiAllMkt(t,regi,emi,emiMkt) ..
   vm_emiAllMkt(t,regi,emi,emiMkt)
-	=e=
-	vm_emiTeMkt(t,regi,emi,emiMkt)
-*** Non-energy sector emissions. Note: These are emissions from all MAC curves. 
-*** So, this includes fugitive emissions, which are sometimes also subsumed under the term energy emissions. 
-	+	sum(emiMacSector2emiMac(emiMacSector,emiMac(emi))$macSector2emiMkt(emiMacSector,emiMkt),
-   	vm_emiMacSector(t,regi,emiMacSector)
-  )
-*** CDR from CDR module
-	+ vm_emiCdr(t,regi,emi)$(sameas(emi,"co2") AND sameas(emiMkt,"ETS")) 
-*** Exogenous emissions
-  +	pm_emiExog(t,regi,emi)$(sameas(emiMkt,"other"))
+  =e=
+    vm_emiTeMkt(t,regi,emi,emiMkt)
+    !! Non-energy sector emissions. Note: These are emissions from all MAC
+    !! curves.  So, this includes fugitive emissions, which are sometimes also
+    !! subsumed under the term energy emissions. 
+  + sum((emiMacSector2emiMac(emiMacSector,emiMac(emi)),
+         macSector2emiMkt(emiMacSector,emiMkt)),
+      vm_emiMacSector(t,regi,emiMacSector)
+    )
+    !! CDR from CDR module
+  + vm_emiCdr(t,regi,emi)$( sameas(emi,"co2") AND sameas(emiMkt,"ETS") ) 
+    !! Exogenous emissions
+  + pm_emiExog(t,regi,emi)$( sameas(emiMkt,"other") )
+    !! non energy emi from chem sector (process emissions from feedstocks):
+  + sum((entyFe2sector2emiMkt_NonEn(entyFe,sector,emiMkt), 
+         se2fe(entySe,entyFe,te)), 
+      vm_demFENonEnergySector(t,regi,entySe,entyFe,sector,emiMkt)
+    * pm_emifacNonEnergy(t,regi,entySe,entyFe,sector,emi)
+    )
+    !!emissions from plastics incineration
+
+    !!emissions from chemical feedstock with unknown fate (assumed to go into the atmosphere)
 ;
 
 

core/input/files

@@ -21,6 +21,7 @@ p_adj_deltacapoffset.cs4r
 p_boundCapCCSindicator.cs4r
 p_boundEmi.cs4r
 p_emiFossilFuelExtr.cs4r
+f_nechem_emissionFactors.cs4r
 p_emineg_econometric.cs3r
 pm_histCap.cs3r
 p_histCapFac.cs4r

core/input/generisdata_emi.prn

@@ -1,9 +1,5 @@
 *** SOF ./core/input/generisdata_emi.prn
-*TA* This file includes the emission factors of secondary energy production technologies given in GtC per ZJ (ie 10^9tC per 10^9TJ) 
-*RP* the emission factors for oil were reduced from 20 to 18.4 to represent 
-***  that a share of oil is not burned but rather converted into long-lived
-***  goods (bitumen->tar for streets, durable plastics - most of plastics is 
-***  quickly burned, but some remains unburnt for decades) - "carbonstored".
+*TA* This file includes the emission factors of secondary energy production technologies given in GtC per ZJ (ie 10^9tC per 10^9TJ)
 ***  NB: reduction calculated based on values in manual to IEA CO2 emissions statistics (2007).
 *JeS* emission factors for n2o are in Tg N per TWa (native remind units)
 *DK* biomass N2O emissions are taken care via p_efFossilFuelExtr in q_macBase in core/equations.gms
@@ -35,15 +31,15 @@ pebiolc.seel.bioigccc         -20.0      20.0
 pebiolc.seel.bioigcc
 pebiolc.segabio.biogasc        -9.5       9.5
 pebiolc.seliqbio.bioftcrec    -12.0      12.0 
-pebiolc.seh2.bioh2c           -22.5      22.5   
-pebiolc.seel.biochp                             
-pebiolc.seh2.bioh2                              
-pebiolc.sehe.biohp                              
-pebiolc.sesobio.biotr                           
-pebiolc.sesobio.biotrmod                        
-pebiolc.segabio.biogas                          
-peoil.seliqfos.refliq          18.4               
-peoil.seel.dot                 18.4             0.01892
+pebiolc.seh2.bioh2c           -22.5      22.5   0.12614
+pebiolc.seel.biochp                             0.12614
+pebiolc.seh2.bioh2                              0.12614
+pebiolc.sehe.biohp                              0.12614
+pebiolc.sesobio.biotr                           0.12614
+pebiolc.sesobio.biotrmod                        0.12614
+pebiolc.segabio.biogas                          0.12614
+peoil.seliqfos.refliq            20               
+peoil.seel.dot                   20             0.01892
 seh2.segasyn.h22ch4                                              0.395
 seh2.seliqsyn.MeOH                                               0.53
 

core/sets.gms

@@ -1737,9 +1737,8 @@ entySeBio(all_enty)       "biomass secondary energy types"
 
 entySeSyn(all_enty)       "synfuel secondary energy types"
 /
-	seliqbio     "secondary energy liquids from biomass"
-	sesobio      "secondary energy solids from biomass"
-	segabio      "secondary energy gas from biomass"
+	seliqsyn     "secondary energy synthetic liquids from H2"
+	segasyn      "secondary energy synthetic gas from H2"
 /
 
 entySeFos(all_enty) "secondary energy types from fossil primary energy"
@@ -2099,6 +2098,14 @@ macSector2emiMkt(all_enty,all_emiMkt)  "mapping mac sectors to emission markets"
         co2luc.other
         co2cement_process.ETS
 /
+
+entyFE2sector2emiMkt_NonEn(all_enty,emi_sectors,all_emiMkt)                     "combinations of FE type, sector and emissions markets existing for FE non-energy use"
+/
+  fesos . indst . ETS
+  fehos . indst . ETS
+  fegas . indst . ETS
+/
+
 ccsCo2(all_enty)    "only cco2 (???)"
 /
         cco2
@@ -2489,8 +2496,14 @@ se2fe(all_enty,all_enty,all_te)   "map secondary energy to end-use energy using
         seliqsyn.fedie.tdsyndie
         seh2.feh2t.tdh2t
 /
+
+
 sefe(all_enty,all_enty) "map secondary energy to final energy"
+//
+
 fete(all_enty,all_te) "map final energy to technologies"
+//
+
 fe2ue(all_enty,all_enty,all_te)    "map FE carriers to ES via appliances"
 //
 

modules/02_welfare/utilitarian/declarations.gms

@@ -23,6 +23,8 @@ $endif.inconv
 $IFTHEN.INCONV_bioSwitch "%cm_INCONV_PENALTY_FESwitch%" == "on"
 p02_inconvPen_Switch_Track(ttot,all_regi)                       "Parameter to track magnitude of inconvenience penalty for bio/synfuel share switching [share of consumption]"
 $ENDIF.INCONV_bioSwitch
+
+p02_inconvPen_NonEn_Switch_Track(ttot,all_regi)                  "Parameter to track magnitude of inconvenience penalty for bio/synfuel non-energy fuel switching [share of consumption]"
 ;
 
 ***-------------------------------------------------------------------------------
@@ -44,6 +46,9 @@ $IFTHEN.INCONV_bioSwitch "%cm_INCONV_PENALTY_FESwitch%" == "on"
 v02_NegInconvPenFeBioSwitch(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt) "Negative inconvenience penalty in the welfare function for bio/synfuel shares switch between sectors and emissions markets"
 v02_PosInconvPenFeBioSwitch(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt) "Positive inconvenience penalty in the welfare function for bio/synfuel shares switch between sectors and emissions markets"
 $ENDIF.INCONV_bioSwitch
+
+v02_NegInconvPenNonEnSwitch(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt) "Negative inconvenience penalty in the welfare function for fuel switching between biomass/synfuel in non-energy use"
+v02_PosInconvPenNonEnSwitch(ttot,all_regi,all_enty,all_enty,emi_sectors,all_emiMkt) "Positive inconvenience penalty in the welfare function for fuel switching between biomass/synfuel in non-energy use"
 ;
 
 ***-------------------------------------------------------------------------------
@@ -62,6 +67,7 @@ $IFTHEN.INCONV_bioSwitch "%cm_INCONV_PENALTY_FESwitch%" == "on"
 q02_inconvPenFeBioSwitch(ttot,all_regi,all_enty,all_enty,all_te,emi_sectors,all_emiMkt)  "Calculate the inconvenience penalty to avoid switching biomass and synfuel shares in hydrocarbons in buildings, transport and industry and emissions markets if costs are relatively close"
 $ENDIF.INCONV_bioSwitch
 
+q02_inconvPenNonEnSwitch(ttot,all_regi,all_enty,all_enty,all_te,emi_sectors,all_emiMkt) "Calculate the inconvenience penalty to avoid switching biomass and synfuel shares in FE non-energy use"    
 ;
 
 *** EOF ./modules/02_welfare/utilitarian/declarations.gms