bjoern/pypsa-eur
1
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases Wiki Activity

add regionalised oil load for process emissions from naphtha as feedstock

Browse Source
This commit is contained in:
chrstphtrs 2023-10-26 11:17:57 +02:00
parent e2b2eafbc1
commit 2ad9ca8f7b
1 changed files with 27 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
scripts/prepare_sector_network.py
View File

@ -2782,9 +2782,25 @@ def add_industry(n, costs):
# convert process emissions from feedstock from MtCO2 to energy demand # convert process emissions from feedstock from MtCO2 to energy demand
# need to aggregate potentials if oil not nodally resolved # need to aggregate potentials if oil not nodally resolved
if options["co2_budget_national"]: if options["co2_budget_national"]:
p_set = demand_factor * (industrial_demand.loc[nodes, "naphtha"] - industrial_demand.loc[nodes, "process emission from feedstock"] / costs.at["oil", "CO2 intensity"]) / nhours p_set_plastics = demand_factor * (industrial_demand.loc[nodes, "naphtha"] - industrial_demand.loc[nodes, "process emission from feedstock"] / costs.at["oil", "CO2 intensity"]) / nhours
else: else:
p_set = demand_factor * (industrial_demand.loc[nodes, "naphtha"] - industrial_demand.loc[nodes, "process emission from feedstock"] / costs.at["oil", "CO2 intensity"]).sum() / nhours p_set_plastics = demand_factor * (industrial_demand.loc[nodes, "naphtha"] - industrial_demand.loc[nodes, "process emission from feedstock"] / costs.at["oil", "CO2 intensity"]).sum() / nhours
if options["co2_budget_national"]:
p_set_process_emissions = (
demand_factor
* (industrial_demand.loc[nodes, "process emission from feedstock"]
/ costs.at["oil", "CO2 intensity"])
/ nhours
)
else:
p_set_process_emissions = (
demand_factor
* (industrial_demand.loc[nodes, "process emission from feedstock"]
/ costs.at["oil", "CO2 intensity"]
).sum()
/ nhours
)
n.madd( n.madd(
"Bus", "Bus",
@ -2799,7 +2815,15 @@ def add_industry(n, costs):
spatial.oil.naphtha, spatial.oil.naphtha,
bus=spatial.oil.naphtha, bus=spatial.oil.naphtha,
carrier="naphtha for industry", carrier="naphtha for industry",
p_set=p_set, p_set=p_set_plastics,
)
n.madd(
"Load",
["naphtha for industry into process emissions from feedstock"],
bus=spatial.oil.nodes,
carrier="naphtha for industry",
p_set=p_set_process_emissions,
) )
n.madd( n.madd(