From ac0d52e1c3542fd6d4fa0ebd639905bac0578a12 Mon Sep 17 00:00:00 2001
From: Tom Brown <tom@nworbmot.org>
Date: Fri, 1 Mar 2024 18:08:28 +0100
Subject: [PATCH] Allow some plastics sequestered in landfill or built
 environment

I.e. assume some fraction of plastics don't end up being burned or
decaying into air.
---
 config/config.default.yaml        |  1 +
 scripts/prepare_sector_network.py | 13 +++++++++----
 2 files changed, 10 insertions(+), 4 deletions(-)

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,
     )