Table of Contents

In a glance

This module generates a shapefile of potential district heating areas based on an assessment of the heat distribution costs. Inputs to the module are heat demand and gross floor area density maps, costs of network expansion, development of heat demand and connection rates, depreciation time, interest rate and a threshold for the accepted heat distribution costs. Furthermore, it calculates the costs of transmission lines between identified district heating areas.

To Top


This calculation module uses the European heat density map (EHDM) and a European gross floor area map (EGFAM), both of which were developed in course of the Hotmaps project, to propose a GIS-based method for determining potential DH areas with specific focus on district heating (DH) grid costs. The DH areas are determined via performing sensitivity analyses on the EHDM under consideration of predefined upper bound of the average distribution costs. The approach additionally allows for the estimation of length and diameter of transmission lines and their associated costs. The outputs are GIS layers that illustrate areas that are economically viable for construction of DH as well as the cost-minimal transmission lines connecting these regions to each other. The calculation module can be used to study the impact of parameters like grid costs ceiling and market share on potential and on expansion and extension of the DH systems.

To Top

Inputs and outputs

The input parameters and layers, as well as output layers and parameters, are as follows.

Input layers and parameters are:

  • First year of investment
  • Last year of investment
  • Depreciation time in years
  • Expected accumulated energy saving: The ratio of energy that is being saved in the last year of investment due to e.g. retrofitting the buildings to the heat demand in the first year of the investment
  • DH market share at the beginning of the investment period
  • Expected DH market share at the end of the investment period
  • Interest rate
  • DH grid cost ceiling in EUR/MWh: In potential DH areas, the distribution grid cost may not exceed the distribution grid cost ceiling.
  • Construction cost constant in EUR/m as well as Construction cost coefficient in EUR/m2 for:
    • Inner city
    • Outer city
    • Park
  • Full load hours
  • Heat density map and gross floor area density map (by default is provided by the toolbox)
    • in raster format (*.tif)
    • with 1 hectare resolution
    • demand densities in MWh/ha and gross floor area densities in m2/ha

Output layers and parameters are:

  • Total demand in the selected region in the first year of investment in MWh
  • Total demand in the selected region in the last year of investment in MWh
  • Maximum potential of DH system through the investment period in MWh
  • Energetic specific DH grid costs in EUR/MWh
  • Energetic specific DH distribution grid costs in EUR/MWh
  • Energetic specific DH transmission grid costs in EUR/MWh
  • Specific DH distribution grid costs per meter in EUR/m
  • Specific DH transmission grid costs per meter in EUR/m
  • Total grid costs - annuity in EUR
  • Total distribution grid costs - annuity in EUR
  • Total transmission grid costs - annuity in EUR
  • Total distribution grid trench length in km
  • Total transmission grid trench length in km
  • Total number of coherent areas
  • Number of economic coherent areas
  • Heat demand density in the last year of the investment period (considering energy savings) in raster format
  • DH areas (both economic and non-economic ones) in shapefile format
  • Transmission lines and their capacities in shapefile format

To Top


Here, a brief explanation of the methodology is provided. For a more complete explanation of the methodology and formulations, please to the paper published about this calculation module [1].

The aim of the calculation module is to find regions in which DH system can be built without exceeding a user-defined average specific cost ceiling in EUR/MWh. This is done under the following assumptions:

  • The economic DH area with highest heat demand is considered as the only available heat source. It produces the heat for itself and all other economic coherent areas.
  • between two DH areas, heat can flow in one direction,
  • the annual DH demand is considered to remain constant after the last year of investment period
  • market share or energy-saving has the same percentages within cells of a DH area and also within different DH areas.
  • The model creates only one connected DH system. It is not possible to have two networks.
  • The input parameter grid cost ceiling is multiplied by ~95% to get distribution grid cost ceiling. This value is used for the determination of potential DH areas.

The determination of economic DH areas is done in three steps. For more details refer to the provided test runs.

STEP 1: Calculation of distribution grid costs based on heat demand and plot ratio using EHDM and EGFAM

STEP 2: Determination of potential DH areas

STEP 3: Determination of economic DH areas and transmission line capacities and configuration required to connect these areas to each other.

Solver Options

This calculation module uses Gurobi solver for solving the optimization problem. In order to guarantee a stable functionality of the calculation module, we have introduced several options for solving the optimization problem. These options are as follows:

  • The gap between the lower and upper objective bound is set to 0.01 (MIPGap = 1e-2).
  • The relative difference between the primal and dual objective value was set to 0.0001 (BarConvTol = 1e-4).
  • The solver focus is set to 1 to find the feasible solutions. Here, the focus is neither optimality nor objective bound (MIPFocus = 1).
  • We have limited the amount of the used RAM to 500 MB in order not to enter to critical situations in case of concurrent runs by users (NodefileStart = 0.5).

To Top

GitHub repository of this calculation module

Here you get the bleeding-edge development for this calculation module.

To Top

Sample run

Here, the calculation module is run for the case study of Vienna, Austria. First, use the "Go To Place" bar to navigate to Vienna and select the city. Click on the "Layers" button to open the "Layers" window and then click on the "CALCULATION MODULE" tab. In the list of calculation modules, select "CM - District heating potential: economic assessment".

To Top

Test Run: default input values for the case study of Vienna

The provided default values in the toolbox are basically suitable for Vienna, i.e. it may not suit for other regions and should be adapted depending on your case study. The calculation is done for the period from 2018 to 2030 (2018 is the year 0 and 2030 is the year 12 and the investment period will be 12 years). The expected accumulated energy saving ratio shows the reduction of heat demand compared to the beginning of the investment period (year 2018). The DH market share refers to the market share within the DH areas. Its value at the beginning of the investment period (year 2018), shows the actual market share (usually known). The expected market share at the investment period is what you expect to reach. This value comes from road maps, scenarios, policies and etc. For the default case, we consider the interest rate of 5 per cent. The DH grid cost ceiling is multiplied by ~95% to get a cost ceiling for the distribution grid. Using this value, the potential DH areas are obtained. Within the potential areas, the average distribution grid cost may not exceed the distribution grid cost ceiling. The value of full load hours is used to estimate the peak load and find a suitable dimension for the transmission grid.

The construction cost constant and the construction cost coefficient originate from reference [2]. The obtained regions are very sensitive to these values. Therefore, as a general comment, we suggest to calculate with these values first and only if you think these values lead to an over- or underestimation of your results, then modify them.

By default, the heat density map and the gross floor area density map that are provided by the toolbox are used for the calculation. You can use your own uploaded layers for running the calculation. Here, we use default layers.

Now, press the run button and wait until the calculation is finished.

IMPORTANT NOTE: Please note that this calculation module may take several minutes to find the final solution. If your calculation takes very long (more than 10 minutes), select a smaller region for your calculation. Also, using arbitrary values can lead to a long calculation time. Therefore, make sure that your provided values are suitable for the selected region.

The following figure shows the obtained results for the given input parameters in Vienna. The most important indicators are demonstrated in the RESULTS window. Additionally, you can get some indicators by pressing on each single potential areas on the map

Figure 1

The output layers will appear in the LAYERS window under the Calculation module section.

To Top


[1]. Fallahnejad M, Hartner M, Kranzl L, Fritz S. Impact of distribution and transmission investment costs of district heating systems on district heating potential. Energy Procedia 2018;149:141–50. doi:10.1016/j.egypro.2018.08.178.

[2]. Persson U, Werner S. Heat distribution and the future competitiveness of district heating. Appl Energy 2011;88:568–76. doi:10.1016/j.apenergy.2010.09.020.

To Top

How to cite

Mostafa Fallahnejad, in Hotmaps-Wiki, CM-District-heating-grid-costs (April 2019)

To Top

Authors and reviewers

This page is written by Mostafa Fallahnejad*.

  • [ ] This page was reviewed by Lukas Kranzl*.

* Energy Economics Group - TU Wien Institute of Energy Systems and Electrical Drives Gusshausstrasse 27-29/370 1040 Wien

To Top


Copyright © 2016-2019: Mostafa Fallahnejad

Creative Commons Attribution 4.0 International License This work is licensed under a Creative Commons CC BY 4.0 International License.

SPDX-License-Identifier: CC-BY-4.0


To Top


We would like to convey our deepest appreciation to the Horizon 2020 Hotmaps Project (Grant Agreement number 723677), which provided the funding to carry out the present investigation.

To Top

View in another language:


* machine translated