diff --git a/config/config.default.yaml b/config/config.default.yaml index 42fe5e2a..d4a8d21d 100644 --- a/config/config.default.yaml +++ b/config/config.default.yaml @@ -667,6 +667,7 @@ industry: 2040: 0.12 2045: 0.16 2050: 0.20 + HVC_environment_sequestration_fraction: 0. sector_ratios_fraction_future: 2020: 0.0 2025: 0.1 diff --git a/scripts/prepare_sector_network.py b/scripts/prepare_sector_network.py index 432d9a82..c308f786 100755 --- a/scripts/prepare_sector_network.py +++ b/scripts/prepare_sector_network.py @@ -2865,7 +2865,7 @@ def add_industry(n, costs): if demand_factor != 1: logger.warning(f"Changing HVC demand by {demand_factor*100-100:+.2f}%.") - p_set_plastics = ( + p_set_naphtha = ( demand_factor * industrial_demand.loc[nodes, "naphtha"].rename( lambda x: x + " naphtha for industry" @@ -2874,7 +2874,7 @@ def add_industry(n, costs): ) if not options["regional_oil_demand"]: - p_set_plastics = p_set_plastics.sum() + p_set_naphtha = p_set_naphtha.sum() n.madd( "Bus", @@ -2889,7 +2889,7 @@ def add_industry(n, costs): spatial.oil.naphtha, bus=spatial.oil.naphtha, carrier="naphtha for industry", - p_set=p_set_plastics, + p_set=p_set_naphtha, ) # some CO2 from naphtha are process emissions from steam cracker @@ -2900,6 +2900,11 @@ def add_industry(n, costs): ) emitted_co2_per_naphtha = costs.at["oil", "CO2 intensity"] - process_co2_per_naphtha + non_sequestered = 1 - get( + snakemake.config["industry"]["HVC_environment_sequestration_fraction"], investment_year + ) + + n.madd( "Link", spatial.oil.naphtha, @@ -2909,7 +2914,7 @@ def add_industry(n, costs): bus3=spatial.co2.process_emissions, carrier="naphtha for industry", p_nom_extendable=True, - efficiency2=emitted_co2_per_naphtha, + efficiency2=emitted_co2_per_naphtha*non_sequestered, efficiency3=process_co2_per_naphtha, )