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

industry: subtract today's ind elec demand, add back new demand

Browse Source 
Since today's industrial electricity demand is distributed by
population and GDP, subtract this from the regular electricity demand
(which already has space/water heating subtracted).

Now regular electricity demand is only non-heating electricity demand
in residential and tertiary sectors.

Add back new industry electricity demand at the correct locations, as
determined using the hotmaps database.
This commit is contained in:
Tom Brown 2020-10-12 14:56:41 +02:00
parent 80cbe98630
commit e8b923e076
1 changed files with 10 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
scripts/prepare_sector_network.py
View File

@ -666,7 +666,7 @@ def insert_electricity_distribution_grid(network):
capital_cost=costs.at['electricity distribution grid','fixed']*snakemake.config["sector"]['electricity_distribution_grid_cost_factor']) capital_cost=costs.at['electricity distribution grid','fixed']*snakemake.config["sector"]['electricity_distribution_grid_cost_factor'])
#this catches regular electricity load and "industry new electricity" #this catches regular electricity load and "industry electricity"
loads = network.loads.index[network.loads.carrier.str.contains("electricity")] loads = network.loads.index[network.loads.carrier.str.contains("electricity")]
network.loads.loc[loads,"bus"] += " low voltage" network.loads.loc[loads,"bus"] += " low voltage"
@ -1635,12 +1635,18 @@ def add_industry(network):
carrier="low-temperature heat for industry", carrier="low-temperature heat for industry",
p_set=industrial_demand.loc[nodes,"low-temperature heat"]/8760.) p_set=industrial_demand.loc[nodes,"low-temperature heat"]/8760.)
#remove today's industrial electricity demand by scaling down total electricity demand
for ct in n.buses.country.unique():
loads = n.loads.index[(n.loads.index.str[:2] == ct) & (n.loads.carrier == "electricity")]
factor = 1 - industrial_demand.loc[loads,"current electricity"].sum()/n.loads_t.p_set[loads].sum().sum()
n.loads_t.p_set[loads] *= factor
network.madd("Load", network.madd("Load",
nodes, nodes,
suffix=" industry new electricity", suffix=" industry electricity",
bus=nodes, bus=nodes,
carrier="industry new electricity", carrier="industry electricity",
p_set = (industrial_demand.loc[nodes,"electricity"]-industrial_demand.loc[nodes,"current electricity"])/8760.) p_set=industrial_demand.loc[nodes,"electricity"]/8760.)
network.madd("Bus", network.madd("Bus",
["process emissions"], ["process emissions"],