c3151902f6
* address geopandas and pandas deprecations/warnings * compatibility with geopandas v1 * [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci --------- Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
91 lines
2.8 KiB
Python
91 lines
2.8 KiB
Python
# -*- coding: utf-8 -*-
|
|
# SPDX-FileCopyrightText: : 2020-2024 The PyPSA-Eur Authors
|
|
#
|
|
# SPDX-License-Identifier: MIT
|
|
"""
|
|
Build hourly heat demand time series from daily heat demand.
|
|
|
|
Water and space heating demand profiles are generated using intraday profiles from BDEW. Different profiles are used for the residential and services sectors as well as weekdays and weekend.
|
|
|
|
The daily heat demand is multiplied by the intraday profile to obtain the hourly heat demand time series. The rule is executed in ``build_sector.smk``.
|
|
|
|
|
|
Relevant Settings
|
|
-----------------
|
|
|
|
.. code:: yaml
|
|
|
|
snapshots:
|
|
drop_leap_day:
|
|
|
|
Inputs
|
|
------
|
|
|
|
- ``data/heat_load_profile_BDEW.csv``: Intraday heat profile for water and space heating demand for the residential and services sectors for weekends and weekdays.
|
|
- ``resources/daily_heat_demand_<scope>_elec_s<simpl>_<clusters>.nc``: Daily heat demand per cluster.
|
|
|
|
Outputs
|
|
-------
|
|
|
|
- ``resources/hourly_heat_demand_<scope>_elec_s<simpl>_<clusters>.nc``:
|
|
"""
|
|
|
|
from itertools import product
|
|
|
|
import pandas as pd
|
|
import xarray as xr
|
|
from _helpers import generate_periodic_profiles, get_snapshots, set_scenario_config
|
|
|
|
if __name__ == "__main__":
|
|
if "snakemake" not in globals():
|
|
from _helpers import mock_snakemake
|
|
|
|
snakemake = mock_snakemake(
|
|
"build_hourly_heat_demand",
|
|
scope="total",
|
|
simpl="",
|
|
clusters=5,
|
|
)
|
|
set_scenario_config(snakemake)
|
|
|
|
snapshots = get_snapshots(
|
|
snakemake.params.snapshots, snakemake.params.drop_leap_day
|
|
)
|
|
|
|
daily_space_heat_demand = (
|
|
xr.open_dataarray(snakemake.input.heat_demand)
|
|
.to_pandas()
|
|
.reindex(index=snapshots, method="ffill")
|
|
)
|
|
|
|
intraday_profiles = pd.read_csv(snakemake.input.heat_profile, index_col=0)
|
|
|
|
sectors = ["residential", "services"]
|
|
uses = ["water", "space"]
|
|
|
|
heat_demand = {}
|
|
for sector, use in product(sectors, uses):
|
|
weekday = list(intraday_profiles[f"{sector} {use} weekday"])
|
|
weekend = list(intraday_profiles[f"{sector} {use} weekend"])
|
|
weekly_profile = weekday * 5 + weekend * 2
|
|
intraday_year_profile = generate_periodic_profiles(
|
|
daily_space_heat_demand.index.tz_localize("UTC"),
|
|
nodes=daily_space_heat_demand.columns,
|
|
weekly_profile=weekly_profile,
|
|
)
|
|
|
|
if use == "space":
|
|
heat_demand[f"{sector} {use}"] = (
|
|
daily_space_heat_demand * intraday_year_profile
|
|
)
|
|
else:
|
|
heat_demand[f"{sector} {use}"] = intraday_year_profile
|
|
|
|
heat_demand = pd.concat(heat_demand, axis=1, names=["sector use", "node"])
|
|
|
|
heat_demand.index.name = "snapshots"
|
|
|
|
ds = heat_demand.stack(future_stack=True).to_xarray()
|
|
|
|
ds.to_netcdf(snakemake.output.heat_demand)
|