Skip to content

Couple TES temperatures to district heating supply temperatures#1612

Merged
fneum merged 65 commits intomasterfrom
e_max_pu_ptes
Mar 31, 2025
Merged

Couple TES temperatures to district heating supply temperatures#1612
fneum merged 65 commits intomasterfrom
e_max_pu_ptes

Conversation

@cpschau
Copy link
Copy Markdown
Contributor

@cpschau cpschau commented Mar 25, 2025 •
edited
Loading

Changes proposed in this Pull Request

This PR implicitly couples the top and bottom layer of PTES with the forward and return flow temperature of the local district heating network. Assuming a direct integration of the storage system, the bottom temperature equals the return flow temperature. The top temperature equals the forward flow temperature in the domain up to 90°C. At forward flow temperatures above 90°C the top temperature of the storage is kept at this operational limit of 90°C.

The variations of the network temperature lead to a dynamic storage capacity, that is described by the e_max_pu parameter/attribute of the Store component. According to Equation (1) in Sorknaes (2018), the dynamic storage capacity linearly depends on the deltaT between top and bottom temperature. The DEA catalogue indicates a deltaT of 55 K (top temperature of 90 and bottom temperature of 35 °C) as reference for the model. These reference values can be varied in the config settings sector:district_heating:ptes:max_top_temperature and sector:district_heating:ptes:min_bottom_temperature . As the spread between forward and return flow is smaller for a majority of systems, especially in summer, the effective storage capacity is smaller than the optimal value for the decision variable e_nom_opt.

2050, all countries (39 nodes), 8H resolution
image

The dynamic variation of the pit storage capacity can be turned off by setting the config parameter sector:district_heating:ptes:dynamic_capacity to false. The comparison of the SOC between a storage with a static and a dynamic storage capacity from an aggregated European perspective shows, that in the latter case the maximum SOC does never come close to the nominal value of the energetic capacity. It remains below 50%.

2050, all countries (39 nodes), 8H resolution
image

The changes in the storage operation are also characterized by a higher amount of cycles. Overall, the reduction in average PTES capacity makes it less attractive, leading to a reduction in PTES investments and an increase in total system costs.

System Costs [bn€/yr] PTES Capacity [TWh] PTES Charger Capacity [GW] PTES Cycles [1/yr]
Static 668.998 136.687 911.246 4.0134
Dynamic 674.144 97.2345 648.23 5.93258
Difference [%] 0.769166 -28.8633 -28.8633 47.8193

Checklist

  • I tested my contribution locally and it works as intended.
  • Code and workflow changes are sufficiently documented.
  • Changed dependencies are added to envs/environment.yaml.
  • Changes in configuration options are added in config/config.default.yaml.
  • Changes in configuration options are documented in doc/configtables/*.csv.
  • Sources of newly added data are documented in doc/data_sources.rst.
  • A release note doc/release_notes.rst is added.

TomKae00 and others added 30 commits February 12, 2025 16:52
@TomKae00
introduce etpr for TES using stores
55f4d0e
@TomKae00
Master into branch
6f684d3
@TomKae00
add constraints
2f67385
@TomKae00
copleted implementation of constraints for TES stores
8dedc43
@TomKae00
Merge branch 'master' into ltes_stores
2a44a18
@TomKae00
updated release notes
ca2d24c
@TomKae00
adjusted mock snakemake default
3576274
@TomKae00
library change
e40fcf6
@TomKae00
merge master into branch
a16db94
@TomKae00
update release notes
b46ec98
@pre-commit-ci
[pre-commit.ci] auto fixes from pre-commit.com hooks
ee3b5f9 
for more information, see https://pre-commit.ci
@TomKae00
Delete .DS_Store
84ebe30
@TomKae00
change retrieval of etpr values in etpr_constraint
df3fc81
@TomKae00
Merge branch 'ltes_stores' of https://github.com/TomKae00/pypsa-eur i…
9662bf2 
…nto ltes_stores # Bitte geben Sie eine Commit-Beschreibung ein, um zu erklären, warum dieser # Merge erforderlich ist, insbesondere wenn es einen aktualisierten # Upstream-Branch mit einem Thema-Branch zusammenführt. # # Zeilen,

die mit '#' beginnen, werden ignoriert, # und eine leere Beschreibung bricht den Commit ab.
@TomKae00
changed TES constraints check based on the configs
599867d
@TomKae00
constraint checks
69e0809
@TomKae00
enhance error handling for the tes constraints
68740c6
@TomKae00
Merge branch 'master' into ltes_stores
cf9c4e4
@TomKae00
small changes
b538ee7
@pre-commit-ci
[pre-commit.ci] auto fixes from pre-commit.com hooks
ef05a0e 
for more information, see https://pre-commit.ci
@TomKae00
changed mock snakemake back to default settings
ec0f962
@TomKae00
Merge branch 'ltes_stores' of https://github.com/TomKae00/pypsa-eur i…
7933ec2 
…nto ltes_stores
@amos-schledorn
Merge branch 'master' into ltes_stores
8ab2e84
@amos-schledorn
Merge branch 'master' into ltes_stores
bcc9f67
@amos-schledorn
feat: add marginal cost for water tank chargers
c45b19b
@pre-commit-ci
[pre-commit.ci] auto fixes from pre-commit.com hooks
b231dd3 
for more information, see https://pre-commit.ci
@amos-schledorn
docs: update add_retrofit_gas_boiler_constraint docs
53e778a
@amos-schledorn
docs: update return type annotation for get_heat_system_type method
3052766
@amos-schledorn
feat: add heat vents for decentral and rural heating
d76df97
@amos-schledorn
Merge remote-tracking branch 'origin/price-water-tank-chargers' into …
20b734c 
…price-water-tank-chargers
@cpschau cpschau requested a review from amos-schledorn March 25, 2025 14:05
@amos-schledorn
Copy link
Copy Markdown
Contributor

amos-schledorn commented Mar 25, 2025

Great feature and very useful plots! Striking by how much the storage capacity drops in summer.
Not sure I understand the time series plots. What does e_max_py (weighted) on the secondary y-axis mean exactly?

Code-wise, I would suggest to

  • move this to a separate rule that either returns flat ones or the dynamic e_max_pu
  • make 90/35C config settings or technology-data entries
  • rename the first e_max_pu in the if-clause to delta_t or similar (or just remove it and use e_max_pu = (t_top - t_bottom) / (t_top_nom - t_bottom_nom)
  • add the Sorknaes paper and other relevant sources as references

Could you also add release notes and update the docs?

@cpschau
Copy link
Copy Markdown
Contributor Author

cpschau commented Mar 26, 2025

Great feature and very useful plots! Striking by how much the storage capacity drops in summer. Not sure I understand the time series plots. What does e_max_py (weighted) on the secondary y-axis mean exactly?

Code-wise, I would suggest to

* move this to a separate rule that either returns flat ones or the dynamic `e_max_pu`

* make 90/35C config settings or `technology-data` entries

* rename the first `e_max_pu` in the if-clause to `delta_t` or similar (or just remove it and use `e_max_pu = (t_top - t_bottom) / (t_top_nom  - t_bottom_nom)`

* add the Sorknaes paper and other relevant sources as references

Could you also add release notes and update the docs?

I addressed all of your points @amos-schledorn. In the plot you asked about, the e_max_pu timeseries of the pit components are weighted by their optimal capacity values to aggregate them to a value between 0 and 1.

@cpschau cpschau requested a review from fneum March 26, 2025 14:07
amos-schledorn
amos-schledorn requested changes Mar 27, 2025
View reviewed changes
Copy link
Copy Markdown
Contributor

@amos-schledorn amos-schledorn left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks really good! I just have a few minor comments. Crazy by how much p_max_pu is reduced - we should look into validating DH supply temperatures.

remove unused function arguments and snakemake inputs, add filepath t…
28d8ece 
…o prevent call of snakemake object within function, revert mock_snakemake arguments
@cpschau cpschau requested a review from amos-schledorn March 27, 2025 14:37
fneum
fneum reviewed Mar 27, 2025
View reviewed changes
Copy link
Copy Markdown
Member

@fneum fneum left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks great! Just a few style comments.

amos-schledorn
amos-schledorn approved these changes Mar 28, 2025
View reviewed changes
Copy link
Copy Markdown
Contributor

@amos-schledorn amos-schledorn left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for the changes! Can be merged once Fabian's comments are addressed, imo.

@fneum fneum merged commit 6e25d96 into master Mar 31, 2025
18 of 19 checks passed
@fneum fneum deleted the e_max_pu_ptes branch March 31, 2025 11:42
@cpschau cpschau mentioned this pull request Jul 24, 2025
Closed
5 tasks
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Reviewers

@fneum fneum fneum approved these changes

@amos-schledorn amos-schledorn amos-schledorn approved these changes

Assignees

No one assigned

Labels

None yet

Projects

None yet

Milestone

No milestone

Development

Successfully merging this pull request may close these issues.

4 participants

@cpschau @amos-schledorn @fneum @TomKae00