|
1
|
|
|
"""The central module containing all code dealing with power plant data. |
|
2
|
|
|
""" |
|
3
|
|
|
from pathlib import Path |
|
4
|
|
|
|
|
5
|
|
|
from geoalchemy2 import Geometry |
|
6
|
|
|
from sqlalchemy import BigInteger, Column, Float, Integer, Sequence, String |
|
7
|
|
|
from sqlalchemy.dialects.postgresql import JSONB |
|
8
|
|
|
from sqlalchemy.ext.declarative import declarative_base |
|
9
|
|
|
from sqlalchemy.orm import sessionmaker |
|
10
|
|
|
import geopandas as gpd |
|
11
|
|
|
import pandas as pd |
|
12
|
|
|
|
|
13
|
|
|
from egon.data import config, db |
|
14
|
|
|
from egon.data.datasets import Dataset |
|
15
|
|
|
from egon.data.datasets.power_plants import assign_voltage_level |
|
16
|
|
|
from egon.data.datasets.storages.home_batteries import ( |
|
17
|
|
|
allocate_home_batteries_to_buildings, |
|
18
|
|
|
) |
|
19
|
|
|
from egon.data.datasets.storages.pumped_hydro import ( |
|
20
|
|
|
apply_voltage_level_thresholds, |
|
21
|
|
|
get_location, |
|
22
|
|
|
match_storage_units, |
|
23
|
|
|
select_mastr_pumped_hydro, |
|
24
|
|
|
select_nep_pumped_hydro, |
|
25
|
|
|
) |
|
26
|
|
|
|
|
27
|
|
|
Base = declarative_base() |
|
28
|
|
|
|
|
29
|
|
|
|
|
30
|
|
View Code Duplication |
class EgonStorages(Base): |
|
|
|
|
|
|
31
|
|
|
__tablename__ = "egon_storages" |
|
32
|
|
|
__table_args__ = {"schema": "supply"} |
|
33
|
|
|
id = Column(BigInteger, Sequence("storage_seq"), primary_key=True) |
|
34
|
|
|
sources = Column(JSONB) |
|
35
|
|
|
source_id = Column(JSONB) |
|
36
|
|
|
carrier = Column(String) |
|
37
|
|
|
el_capacity = Column(Float) |
|
38
|
|
|
bus_id = Column(Integer) |
|
39
|
|
|
voltage_level = Column(Integer) |
|
40
|
|
|
scenario = Column(String) |
|
41
|
|
|
geom = Column(Geometry("POINT", 4326)) |
|
42
|
|
|
|
|
43
|
|
|
|
|
44
|
|
|
class Storages(Dataset): |
|
45
|
|
|
def __init__(self, dependencies): |
|
46
|
|
|
super().__init__( |
|
47
|
|
|
name="Storages", |
|
48
|
|
|
version="0.0.4", |
|
49
|
|
|
dependencies=dependencies, |
|
50
|
|
|
tasks=( |
|
51
|
|
|
create_tables, |
|
52
|
|
|
allocate_pumped_hydro_eGon2035, |
|
53
|
|
|
allocate_pumped_hydro_eGon100RE, |
|
54
|
|
|
allocate_pv_home_batteries_to_grids, |
|
55
|
|
|
allocate_home_batteries_to_buildings, |
|
56
|
|
|
), |
|
57
|
|
|
) |
|
58
|
|
|
|
|
59
|
|
|
|
|
60
|
|
|
def create_tables(): |
|
61
|
|
|
"""Create tables for power plant data |
|
62
|
|
|
Returns |
|
63
|
|
|
------- |
|
64
|
|
|
None. |
|
65
|
|
|
""" |
|
66
|
|
|
|
|
67
|
|
|
cfg = config.datasets()["storages"] |
|
68
|
|
|
db.execute_sql(f"CREATE SCHEMA IF NOT EXISTS {cfg['target']['schema']};") |
|
69
|
|
|
engine = db.engine() |
|
70
|
|
|
db.execute_sql( |
|
71
|
|
|
f"""DROP TABLE IF EXISTS |
|
72
|
|
|
{cfg['target']['schema']}.{cfg['target']['table']}""" |
|
73
|
|
|
) |
|
74
|
|
|
|
|
75
|
|
|
db.execute_sql("""DROP SEQUENCE IF EXISTS pp_seq""") |
|
76
|
|
|
EgonStorages.__table__.create(bind=engine, checkfirst=True) |
|
77
|
|
|
|
|
78
|
|
|
|
|
79
|
|
|
def allocate_pumped_hydro_eGon2035(export=True): |
|
80
|
|
|
"""Allocates pumped_hydro plants for eGon2035 scenario and either exports |
|
81
|
|
|
results to data base or returns as a dataframe |
|
82
|
|
|
|
|
83
|
|
|
Parameters |
|
84
|
|
|
---------- |
|
85
|
|
|
export : bool |
|
86
|
|
|
Choose if allocated pumped hydro plants should be exported to the data |
|
87
|
|
|
base. The default is True. |
|
88
|
|
|
If export=False a data frame will be returned |
|
89
|
|
|
|
|
90
|
|
|
Returns |
|
91
|
|
|
------- |
|
92
|
|
|
power_plants : pandas.DataFrame |
|
93
|
|
|
List of pumped hydro plants in 'eGon2035' scenario |
|
94
|
|
|
""" |
|
95
|
|
|
|
|
96
|
|
|
carrier = "pumped_hydro" |
|
97
|
|
|
|
|
98
|
|
|
cfg = config.datasets()["power_plants"] |
|
99
|
|
|
|
|
100
|
|
|
nep = select_nep_pumped_hydro() |
|
101
|
|
|
mastr = select_mastr_pumped_hydro() |
|
102
|
|
|
|
|
103
|
|
|
# Assign voltage level to MaStR |
|
104
|
|
|
mastr["voltage_level"] = assign_voltage_level( |
|
105
|
|
|
mastr.rename({"el_capacity": "Nettonennleistung"}, axis=1), cfg |
|
106
|
|
|
) |
|
107
|
|
|
|
|
108
|
|
|
# Initalize DataFrame for matching power plants |
|
109
|
|
|
matched = gpd.GeoDataFrame( |
|
110
|
|
|
columns=[ |
|
111
|
|
|
"carrier", |
|
112
|
|
|
"el_capacity", |
|
113
|
|
|
"scenario", |
|
114
|
|
|
"geometry", |
|
115
|
|
|
"MaStRNummer", |
|
116
|
|
|
"source", |
|
117
|
|
|
"voltage_level", |
|
118
|
|
|
] |
|
119
|
|
|
) |
|
120
|
|
|
|
|
121
|
|
|
# Match pumped_hydro units from NEP list |
|
122
|
|
|
# using PLZ and capacity |
|
123
|
|
|
matched, mastr, nep = match_storage_units( |
|
124
|
|
|
nep, mastr, matched, buffer_capacity=0.1, consider_carrier=False |
|
125
|
|
|
) |
|
126
|
|
|
|
|
127
|
|
|
# Match plants from NEP list using plz, |
|
128
|
|
|
# neglecting the capacity |
|
129
|
|
|
matched, mastr, nep = match_storage_units( |
|
130
|
|
|
nep, |
|
131
|
|
|
mastr, |
|
132
|
|
|
matched, |
|
133
|
|
|
consider_location="plz", |
|
134
|
|
|
consider_carrier=False, |
|
135
|
|
|
consider_capacity=False, |
|
136
|
|
|
) |
|
137
|
|
|
|
|
138
|
|
|
# Match plants from NEP list using city, |
|
139
|
|
|
# neglecting the capacity |
|
140
|
|
|
matched, mastr, nep = match_storage_units( |
|
141
|
|
|
nep, |
|
142
|
|
|
mastr, |
|
143
|
|
|
matched, |
|
144
|
|
|
consider_location="city", |
|
145
|
|
|
consider_carrier=False, |
|
146
|
|
|
consider_capacity=False, |
|
147
|
|
|
) |
|
148
|
|
|
|
|
149
|
|
|
# Match remaining plants from NEP using the federal state |
|
150
|
|
|
matched, mastr, nep = match_storage_units( |
|
151
|
|
|
nep, |
|
152
|
|
|
mastr, |
|
153
|
|
|
matched, |
|
154
|
|
|
buffer_capacity=0.1, |
|
155
|
|
|
consider_location="federal_state", |
|
156
|
|
|
consider_carrier=False, |
|
157
|
|
|
) |
|
158
|
|
|
|
|
159
|
|
|
# Match remaining plants from NEP using the federal state |
|
160
|
|
|
matched, mastr, nep = match_storage_units( |
|
161
|
|
|
nep, |
|
162
|
|
|
mastr, |
|
163
|
|
|
matched, |
|
164
|
|
|
buffer_capacity=0.7, |
|
165
|
|
|
consider_location="federal_state", |
|
166
|
|
|
consider_carrier=False, |
|
167
|
|
|
) |
|
168
|
|
|
|
|
169
|
|
|
print(f"{matched.el_capacity.sum()} MW of {carrier} matched") |
|
170
|
|
|
print(f"{nep.c2035_capacity.sum()} MW of {carrier} not matched") |
|
171
|
|
|
|
|
172
|
|
|
if nep.c2035_capacity.sum() > 0: |
|
173
|
|
|
|
|
174
|
|
|
# Get location using geolocator and city information |
|
175
|
|
|
located, unmatched = get_location(nep) |
|
176
|
|
|
|
|
177
|
|
|
# Bring both dataframes together |
|
178
|
|
|
matched = matched.append( |
|
179
|
|
|
located[ |
|
180
|
|
|
[ |
|
181
|
|
|
"carrier", |
|
182
|
|
|
"el_capacity", |
|
183
|
|
|
"scenario", |
|
184
|
|
|
"geometry", |
|
185
|
|
|
"source", |
|
186
|
|
|
"MaStRNummer", |
|
187
|
|
|
] |
|
188
|
|
|
], |
|
189
|
|
|
ignore_index=True, |
|
190
|
|
|
) |
|
191
|
|
|
|
|
192
|
|
|
# Set CRS |
|
193
|
|
|
matched.crs = "EPSG:4326" |
|
194
|
|
|
|
|
195
|
|
|
# Assign voltage level |
|
196
|
|
|
matched = apply_voltage_level_thresholds(matched) |
|
197
|
|
|
|
|
198
|
|
|
# Assign bus_id |
|
199
|
|
|
# Load grid district polygons |
|
200
|
|
|
mv_grid_districts = db.select_geodataframe( |
|
201
|
|
|
f""" |
|
202
|
|
|
SELECT * FROM {cfg['sources']['egon_mv_grid_district']} |
|
203
|
|
|
""", |
|
204
|
|
|
epsg=4326, |
|
205
|
|
|
) |
|
206
|
|
|
|
|
207
|
|
|
ehv_grid_districts = db.select_geodataframe( |
|
208
|
|
|
f""" |
|
209
|
|
|
SELECT * FROM {cfg['sources']['ehv_voronoi']} |
|
210
|
|
|
""", |
|
211
|
|
|
epsg=4326, |
|
212
|
|
|
) |
|
213
|
|
|
|
|
214
|
|
|
# Perform spatial joins for plants in ehv and hv level seperately |
|
215
|
|
|
power_plants_hv = gpd.sjoin( |
|
216
|
|
|
matched[matched.voltage_level >= 3], |
|
217
|
|
|
mv_grid_districts[["bus_id", "geom"]], |
|
218
|
|
|
how="left", |
|
219
|
|
|
).drop(columns=["index_right"]) |
|
220
|
|
|
power_plants_ehv = gpd.sjoin( |
|
221
|
|
|
matched[matched.voltage_level < 3], |
|
222
|
|
|
ehv_grid_districts[["bus_id", "geom"]], |
|
223
|
|
|
how="left", |
|
224
|
|
|
).drop(columns=["index_right"]) |
|
225
|
|
|
|
|
226
|
|
|
# Combine both dataframes |
|
227
|
|
|
power_plants = pd.concat([power_plants_hv, power_plants_ehv]) |
|
228
|
|
|
|
|
229
|
|
|
# Delete existing units in the target table |
|
230
|
|
|
db.execute_sql( |
|
231
|
|
|
f""" DELETE FROM {cfg ['target']['schema']}.{cfg ['target']['table']} |
|
232
|
|
|
WHERE carrier IN ('pumped_hydro') |
|
233
|
|
|
AND scenario='eGon2035';""" |
|
234
|
|
|
) |
|
235
|
|
|
|
|
236
|
|
|
# If export = True export pumped_hydro plants to data base |
|
237
|
|
|
|
|
238
|
|
|
if export: |
|
239
|
|
|
# Insert into target table |
|
240
|
|
|
session = sessionmaker(bind=db.engine())() |
|
241
|
|
|
for i, row in power_plants.iterrows(): |
|
242
|
|
|
entry = EgonStorages( |
|
243
|
|
|
sources={"el_capacity": row.source}, |
|
244
|
|
|
source_id={"MastrNummer": row.MaStRNummer}, |
|
245
|
|
|
carrier=row.carrier, |
|
246
|
|
|
el_capacity=row.el_capacity, |
|
247
|
|
|
voltage_level=row.voltage_level, |
|
248
|
|
|
bus_id=row.bus_id, |
|
249
|
|
|
scenario=row.scenario, |
|
250
|
|
|
geom=f"SRID=4326;POINT({row.geometry.x} {row.geometry.y})", |
|
251
|
|
|
) |
|
252
|
|
|
session.add(entry) |
|
253
|
|
|
session.commit() |
|
254
|
|
|
|
|
255
|
|
|
else: |
|
256
|
|
|
return power_plants |
|
257
|
|
|
|
|
258
|
|
|
|
|
259
|
|
|
def allocate_pumped_hydro_eGon100RE(): |
|
260
|
|
|
"""Allocates pumped_hydro plants for eGon100RE scenario based on a |
|
261
|
|
|
prox-to-now method applied on allocated pumped-hydro plants in the eGon2035 |
|
262
|
|
|
scenario. |
|
263
|
|
|
|
|
264
|
|
|
Parameters |
|
265
|
|
|
---------- |
|
266
|
|
|
None |
|
267
|
|
|
|
|
268
|
|
|
Returns |
|
269
|
|
|
------- |
|
270
|
|
|
None |
|
271
|
|
|
""" |
|
272
|
|
|
|
|
273
|
|
|
carrier = "pumped_hydro" |
|
274
|
|
|
cfg = config.datasets()["power_plants"] |
|
275
|
|
|
boundary = config.settings()["egon-data"]["--dataset-boundary"] |
|
276
|
|
|
|
|
277
|
|
|
# Select installed capacity for pumped_hydro in eGon100RE scenario from |
|
278
|
|
|
# scenario capacities table |
|
279
|
|
|
capacity = db.select_dataframe( |
|
280
|
|
|
f""" |
|
281
|
|
|
SELECT capacity |
|
282
|
|
|
FROM {cfg['sources']['capacities']} |
|
283
|
|
|
WHERE carrier = '{carrier}' |
|
284
|
|
|
AND scenario_name = 'eGon100RE'; |
|
285
|
|
|
""" |
|
286
|
|
|
) |
|
287
|
|
|
|
|
288
|
|
|
if boundary == "Schleswig-Holstein": |
|
289
|
|
|
# Break capacity of pumped hydron plants down SH share in eGon2035 |
|
290
|
|
|
capacity_phes = capacity.iat[0, 0] * 0.0176 |
|
291
|
|
|
|
|
292
|
|
|
elif boundary == "Everything": |
|
293
|
|
|
# Select national capacity for pumped hydro |
|
294
|
|
|
capacity_phes = capacity.iat[0, 0] |
|
295
|
|
|
|
|
296
|
|
|
else: |
|
297
|
|
|
raise ValueError(f"'{boundary}' is not a valid dataset boundary.") |
|
298
|
|
|
|
|
299
|
|
|
# Get allocation of pumped_hydro plants in eGon2035 scenario as the |
|
300
|
|
|
# reference for the distribution in eGon100RE scenario |
|
301
|
|
|
allocation = allocate_pumped_hydro_eGon2035(export=False) |
|
302
|
|
|
|
|
303
|
|
|
scaling_factor = capacity_phes / allocation.el_capacity.sum() |
|
304
|
|
|
|
|
305
|
|
|
power_plants = allocation.copy() |
|
306
|
|
|
power_plants["scenario"] = "eGon100RE" |
|
307
|
|
|
power_plants["el_capacity"] = allocation.el_capacity * scaling_factor |
|
308
|
|
|
|
|
309
|
|
|
# Insert into target table |
|
310
|
|
|
session = sessionmaker(bind=db.engine())() |
|
311
|
|
|
for i, row in power_plants.iterrows(): |
|
312
|
|
|
entry = EgonStorages( |
|
313
|
|
|
sources={"el_capacity": row.source}, |
|
314
|
|
|
source_id={"MastrNummer": row.MaStRNummer}, |
|
315
|
|
|
carrier=row.carrier, |
|
316
|
|
|
el_capacity=row.el_capacity, |
|
317
|
|
|
voltage_level=row.voltage_level, |
|
318
|
|
|
bus_id=row.bus_id, |
|
319
|
|
|
scenario=row.scenario, |
|
320
|
|
|
geom=f"SRID=4326;POINT({row.geometry.x} {row.geometry.y})", |
|
321
|
|
|
) |
|
322
|
|
|
session.add(entry) |
|
323
|
|
|
session.commit() |
|
324
|
|
|
|
|
325
|
|
|
|
|
326
|
|
|
def home_batteries_per_scenario(scenario): |
|
327
|
|
|
"""Allocates home batteries which define a lower boundary for extendable |
|
328
|
|
|
battery storage units. The overall installed capacity is taken from NEP |
|
329
|
|
|
for eGon2035 scenario. The spatial distribution of installed battery |
|
330
|
|
|
capacities is based on the installed pv rooftop capacity. |
|
331
|
|
|
|
|
332
|
|
|
Parameters |
|
333
|
|
|
---------- |
|
334
|
|
|
None |
|
335
|
|
|
|
|
336
|
|
|
Returns |
|
337
|
|
|
------- |
|
338
|
|
|
None |
|
339
|
|
|
""" |
|
340
|
|
|
|
|
341
|
|
|
cfg = config.datasets()["storages"] |
|
342
|
|
|
dataset = config.settings()["egon-data"]["--dataset-boundary"] |
|
343
|
|
|
|
|
344
|
|
|
if scenario == "eGon2035": |
|
345
|
|
|
|
|
346
|
|
|
target_file = ( |
|
347
|
|
|
Path(".") |
|
348
|
|
|
/ "data_bundle_egon_data" |
|
349
|
|
|
/ "nep2035_version2021" |
|
350
|
|
|
/ cfg["sources"]["nep_capacities"] |
|
351
|
|
|
) |
|
352
|
|
|
|
|
353
|
|
|
capacities_nep = pd.read_excel( |
|
354
|
|
|
target_file, |
|
355
|
|
|
sheet_name="1.Entwurf_NEP2035_V2021", |
|
356
|
|
|
index_col="Unnamed: 0", |
|
357
|
|
|
) |
|
358
|
|
|
|
|
359
|
|
|
# Select target value in MW |
|
360
|
|
|
target = capacities_nep.Summe["PV-Batteriespeicher"]*1000 |
|
361
|
|
|
|
|
362
|
|
|
else: |
|
363
|
|
|
target = db.select_dataframe( |
|
364
|
|
|
f""" |
|
365
|
|
|
SELECT capacity |
|
366
|
|
|
FROM {cfg['sources']['capacities']} |
|
367
|
|
|
WHERE scenario_name = '{scenario}' |
|
368
|
|
|
AND carrier = 'battery'; |
|
369
|
|
|
""" |
|
370
|
|
|
).capacity[0] |
|
371
|
|
|
|
|
372
|
|
|
pv_rooftop = db.select_dataframe( |
|
373
|
|
|
f""" |
|
374
|
|
|
SELECT bus, p_nom, generator_id |
|
375
|
|
|
FROM {cfg['sources']['generators']} |
|
376
|
|
|
WHERE scn_name = '{scenario}' |
|
377
|
|
|
AND carrier = 'solar_rooftop' |
|
378
|
|
|
AND bus IN |
|
379
|
|
|
(SELECT bus_id FROM {cfg['sources']['bus']} |
|
380
|
|
|
WHERE scn_name = '{scenario}' AND country = 'DE' ); |
|
381
|
|
|
""" |
|
382
|
|
|
) |
|
383
|
|
|
|
|
384
|
|
|
if dataset == "Schleswig-Holstein": |
|
385
|
|
|
target = target / 16 |
|
386
|
|
|
|
|
387
|
|
|
battery = pv_rooftop |
|
388
|
|
|
battery["p_nom_min"] = target * battery["p_nom"] / battery["p_nom"].sum() |
|
389
|
|
|
battery = battery.drop(columns=["p_nom"]) |
|
390
|
|
|
|
|
391
|
|
|
battery["carrier"] = "home_battery" |
|
392
|
|
|
battery["scenario"] = scenario |
|
393
|
|
|
|
|
394
|
|
|
if scenario == "eGon2035": |
|
395
|
|
|
source = "NEP" |
|
396
|
|
|
|
|
397
|
|
|
else: |
|
398
|
|
|
source = "p-e-s" |
|
399
|
|
|
|
|
400
|
|
|
battery[ |
|
401
|
|
|
"source" |
|
402
|
|
|
] = f"{source} capacity allocated based in installed PV rooftop capacity" |
|
403
|
|
|
|
|
404
|
|
|
# Insert into target table |
|
405
|
|
|
session = sessionmaker(bind=db.engine())() |
|
406
|
|
|
for i, row in battery.iterrows(): |
|
407
|
|
|
entry = EgonStorages( |
|
408
|
|
|
sources={"el_capacity": row.source}, |
|
409
|
|
|
source_id={"generator_id": row.generator_id}, |
|
410
|
|
|
carrier=row.carrier, |
|
411
|
|
|
el_capacity=row.p_nom_min, |
|
412
|
|
|
bus_id=row.bus, |
|
413
|
|
|
scenario=row.scenario, |
|
414
|
|
|
) |
|
415
|
|
|
session.add(entry) |
|
416
|
|
|
session.commit() |
|
417
|
|
|
|
|
418
|
|
|
|
|
419
|
|
|
def allocate_pv_home_batteries_to_grids(): |
|
420
|
|
|
|
|
421
|
|
|
home_batteries_per_scenario("eGon2035") |
|
422
|
|
|
home_batteries_per_scenario("eGon100RE") |
|
423
|
|
|
|