Skip to main content
  • Original article
  • Open access
  • Published:

Bioenergy plants’ potential for contributing to heat generation in Germany



The Paris Climate Agreement requires a rapid and efficient shift to renewable energies and a decarbonization of the energy system. Combined heat and power provision from biomass is one way to efficiently provide renewable heat. Despite this, many bioenergy plants in Germany are mainly used to generate electricity and the provision of externally usable heat still has untapped potential. In this study, we investigated gross quantities as well as the economically viable potential of Germany’s current bioenergy plant stock in supplying renewable heat.


We used a top-down GIS modeling approach to spatially and explicitly assess the heat demand of three different categories of heat sinks at sub-municipal level. These included residential, commercial, and industrial areas, as well as large individual heat consumers. We then calculated the plant-specific heat sales potential for two different district heating network options. In addition, we developed a method for assessing the economic efficiency of the previously identified technical heat sales volume for a set of 20 different clusters of bioenergy plants.


The results show that about 50% of the bioenergy plants have potential heat consumers in their immediate proximity. The overall technical heat sales potential for all three categories totals around 150 TWhth/a. However, this potential is not evenly distributed throughout Germany. Certain regions appear to be more favorable for investing in district heating networks powered by heat from biomass. The economically viable heat sales potential related to electrical energy generation ranges from − 0.128 to 0.160 €/kWhth.


We concluded that, under certain conditions such as location or supply and demand structure, German bioenergy plants have the potential to provide a significant share to renewable energies in the heating sector. In addition, the heat sales potential is highly relevant for plant operators as the importance of heat as a business segment is set to increase. Furthermore, bioenergy plants could contribute 2.1% (16.3 TWhth/a) to the total demand for space heating in Germany (765 TWhth/a) when considering certain technical and economic constraints.


Heat accounts for half of the final energy consumption in Germany, but currently only 13.9 % of heat is supplied by renewable energy sources (RES) [1]. Thus, the share of renewable energy sources in heat supply needs to increase if the greenhouse gas reduction targets of the Paris Climate Agreement are to be met. In order to achieve this goal, the European Union (EU) intends to analyze different ways to promote a transition towards using RES in supplying heat as part of its strategy on heating and cooling (COM 2016). One option is to decarbonize existing district heating networks (DHNs) and to supply more heat consumers through new, yet-to-be-built DHNs. At the moment, district heating provides 9% of the EU’s heat [2]. In 2012, fossil fuels were the primary energy sources for heating and cooling in these networks (75%). The most important renewable fuel for DHNs, with a share of 12% in 2015, is heat generated from biomass [3].

In 2012, Germany’s share in final energy consumption for heating and cooling supplied by district heating systems was 8-9% [2]. The German district heating systemFootnote 1 has a connected load of 51 GW and a transformation output of around 470 PJ. The overall fuel input amounts to 522.3 PJ and consists of coal (29%), lignite (11%), mineral oil (1%), natural gas (42%), waste (11%), and biomass (6%) [4]. However, in order to achieve the greenhouse gas reduction targets set by the German government, DHNs supplied by RES will need to provide 23% of the final energy consumption for building heat by 2050 [5]. One reason for this is that regions without a natural gas network and with predominantly decentralized oil heating systems also need to have renewable alternatives for their so far fossil-based heat supply [6]. Furthermore, with the foreseeable shutdown of coal-fired power plants in Germany by 2038 [7], the 50 TWhth/a of heat [8] currently produced by this sector needs to be replaced by low-carbon fuels. Fritz and Pehnt [8] recently stated that biomass must also partly contribute to this as well.

The largest share of heat from RES in Germany is generated by biomass (87.3%), mainly in small combustion plants [9]. Additionally, Germany has an installed capacity of 5.6 GWel [10] in biogas plants and 1.6 GWel [10] in plants using solid biomass fuels, which are currently not utilizing all of the available heat [11].

DHNs are a suitable way of increasing the use of RES in the heating sector, as Lund et al. point out [12]. This particularly applies to biomass. Bio-fueled combined heat and power units (CHPU) can achieve an overall efficiency of up to 90% [13]. The co-generated heat can be used to supply DHNs. In order to fully exploit this potential, however, suitable heat consumers (heat sinks) must be located near the plants; otherwise, the generated heat cannot be supplied properly due to transport-related energy losses. The plant-specific and, consequently, the total heat supply potential of the bioenergy plants (BEPs) in Germany thus crucially depends on plant location.

In current energy scenarios for Germany, there is little information available specifically on local district heating networks and the fuels being used there. Koch et al. [14] stated an overall contribution to the generation of heat of 16 TWh with natural gas being the dominant fuel. Biomass only plays a minor role in the supply of these networks but with increasing proportions until 2050 [14].

Against this backdrop, the question arises as to the extent to which the current BEP installations in Germany can contribute to the provision of heat to external users by means of district heating networks. In this study, we investigate the gross amount of heat demand and supply of the German BEP stock as well as the economically viable potential. For this purpose, we seek to answer the following research questions:

  • How many BEPs have nearby heat sinks suitable for a DHN?

  • What is the extent of the technical heat sales potential (tHSP) (in TWhth/a) of these BEPs?

  • How can the economic viability of tHSP be assessed?

  • How high is the economically viable heat sales potential (eHSP)?

To answer these questions, we modeled the heat demand for three different categories of heat sinks at a sub-municipal spatial resolution: the residential sector (1), the industrial sector as well as the sector trade, commerce and services (TCS) (2), and large individual heat consumers (LIHC) such as hospitals and schools (3). We then evaluated which BEPs are capable of delivering heat at all, as well as the plant-specific and overall heat sales potential per category.

Additionally, we developed a method to evaluate the economic viability of the previously assessed technical heat sales for a set of 20 different clusters of BEPs.

In the context of modeling future energy systems, Weinand et al. emphasize the need for data at a high spatial resolution [15]. While the general potential for DHNs in Germany has been analyzed before [16], an investigation into the specific potential of bio-fueled CHPUs has only been carried out in some states, such as Saarland [17]. A Germany-wide assessment has yet to be made.


Modeling scheme

In order to determine the technical heat sales potential (tHSP) for each plant individually, it is necessary to identify suitable heat sinks close to the respective plant. Consequently, the tHSP is defined as the heat demand of suitable heat sinks within a defined maximum distance to the respective BEP. While spatial information on BEPs is available through the bioenergy plant database of the German Biomass Research Centre (DBFZ) [18], the spatially explicit heat demand related to these plants is missing and needs to be identified. The modeling consisted of two separate steps. First, we modeled the heat demand at a spatial resolution below the municipal level. For the residential sector, this means the level of individual residential blocks and for the industrial sector the level of industrial and commercial areas as recorded in the ATKIS Basis Digital Landscape Model (DLM) [19]. The ATKIS Basis DLM offers a description of these topographic elements in vector data format for all of Germany on a scale of 1:25,000 [19]. For the LIHC, this means their exact geographical location. In a second step, we calculated the tHSP for each plant individually and concluded the overall tHSP.

Based on the tHSP, we developed an estimation approach to assess the economically viable heat sales potential (eHSP, see Fig. 5). The basic method for calculating tHSP as well as the input data used is illustrated in Fig. 1.

Fig. 1
figure 1

Modeling scheme to identify plant-specific heat supply potential. The models of the input data are further described in Table 1

Most of the data used in the model is publicly available. The content of the respective datasets, their origin and the latest update are listed in Table 1. Only the municipal level heat demand (later referred to as “IZES heat demand”) and the DBFZ BEP database contain internal data of the respective institutions and are therefore not publicly available.

Table 1 Data used for modeling heat demand

Modeling the heat demand at a high spatial resolution

Heat is required and consumed in three different sectors: residential, TCS, and industrial. As the German Environment Agency states, heat in the residential and TCS sectors is primarily required for space heating, while the industrial sector requires the largest share for process heat [1]. Consequently, the methods for determining the total heat demand per sector and for its spatial modeling vary.

For the residential, TCS and industrial sectors, we followed a top-down modeling approach as described by Fleiter et al. and applied by Baur et al. [3, 17]. The basic principle of this approach is to disaggregate the heat demand from the existing aggregated level into smaller spatial units. In our case, we used data on the heat demand per sector for each municipality, calculated by Baur et al. [20], and allocated it to the residential areas and the commercial and/or industrial areas of each community.

We simultaneously applied a bottom-up approach for LIHCs, such as hospitals, commercial greenhouses, schools, and public outdoor swimming pools, which are particularly suitable for DHNs fueled by heat generated from biomass [26]. This means that we calculated their yearly heat demand on the basis of an object-specific key such as kWhtha/bed (hospitals), kWhtha/student (schools), kWhtha/m2 (commercial greenhouses), and kWhtha/m2 water surface (public outdoor swimming pools).

Residential sector

The heat demand of residential units on a municipal level was assessed using building stock grouped by building type (single family dwelling, townhouse, etc.), year of construction, and municipality. This data was extracted from the 2011 census database [27]. We multiplied the number of buildings, divided into building class and municipality, by the corresponding heat demand in kWhth/a and respective residential unit, derived from the German building typology of IWU 2003 [28]. Additionally, we took an average annual domestic hot water demand of 500 kWhth per capita into account, based on the municipal population on December 31, 2011 [27]. The overall methodology of assessing the heat demand is discussed in more detail in Baur et al. [20].

Within the ATKIS Basis DLM, the utilization of a given area is defined by unique codes related to object type. Subsequently, we selected all objects that are assigned to the residential sector including the object types 2102 (dwelling), 2111 (residential area), and 2113 (area of mixed use). These objects were then labeled as residential area objects (RAO) and are hereinafter referred to as such. All other objects (e.g., parks, bodies of water, train tracks, and stations) were no longer considered in the model and were therefore removed from our analysis. We then calculated the residential area factor (RAF), which represents the share of each individual RAO in the total residential area of the municipality, with

$$ RAF=\frac{\ {A}_{\mathrm{RAO}}}{A_{\mathrm{total}\ \mathrm{residential}\ \mathrm{area}}\ } $$

where RAF is the residential area factor, RAO is the residential area object, and A is the area (m2).

Since the heat demand in the residential sector also depends on the number of inhabitants per area, it was necessary to determine the number of inhabitants for each RAO individually. To do so, we used population data from the 2011 census which is available in a nationwide grid with raster cell sizes of 100 m × 100 mFootnote 2 [25]. The geodata for the grid is provided by the German Federal Agency for Cartography and Geodesy (BKG) [22]. We intersected the RAOs with the population data and calculated the population figures per RAOFootnote 3. Subsequently, we determined the population factor (PF), which is the ratio of any given RAO to the municipality’s total population with

$$ PF=\frac{\ \mathrm{Inhabitants}\ \mathrm{per}\ \mathrm{RAO}}{\mathrm{Inhabitants}\ \mathrm{of}\ \mathrm{municipality}} $$

where PF is the population factor and RAO is the residential area object.

In the next step, we calculated the RAO-specific heat demand using the heat demand per municipality provided by IZES in combination with the previously introduced factors RAF and PF. For RAOs containing information on the number of inhabitants, the specific heat demand (Q) is

$$ {Q}_{\mathrm{RAO}}={Q}_{\mathrm{municipality}}\times PF $$

where Q is the heat demand (GWhth/a), PFis the population factor, and RAO is the residential area object.

For RAOs with missing information on the number of inhabitants it is

$$ {Q}_{RAO}={Q}_{municipality}\ast RAF $$

where Q is the heat demand (GWhth/a), PF is the residential area factor, and RAO is the residential area object.

Finally, the heat density (kWhth/m2a) was determined for each RAO by dividing the RAO-specific heat demand (GWhth/a) by its area.

$$ {HD}_{RAO}=\left(\frac{Q_{RAO}}{A_{RAO}}\right) $$

where HD is the heat density (kWhth/m2a), Q is the heat demand (kWhth/a), A is the area (m2), and RAO is the residential area object.

The result is a dataset for the whole of Germany with very detailed spatial information on the heat demand and heat density of residential areas.

TCS and industrial sector

The heat demand of the industrial sector on a district level (NUTS 3) was estimated by multiplying the number of employees in the manufacturing sectorFootnote 4 by the corresponding average heat demand per employee and industry derived from the national statistic on industrial energy demand [31, 32]. Due to the lack of corresponding data for the TCS sector, we determined the heat demand of the TCS sector by using the gross added values of the related economic sectorsFootnote 5 [33]. They were multiplied by the average heat demand of the corresponding sectors derived from the reference scenario, which is described by Kirchner et al. [34]. In order to also determine the heat demand of both the TCS and industrial sector on a municipal level (LAU 2), we multiplied the calculated heat demands by the ratio of the number of employees at a municipal and district level [35].

As with the residential sector, the DLM was the basis for modeling the heat demand for the TCS and industrial sectors. Here, we used the object types “Industrial and commercial area” (2112) and “Area with special functional characteristics” (2114). We recognized that the DLM alone is only partially suitable for spatially modeling the heat demand of industrial and TCS areas. There are numerous objects (areas) assigned in the classification 2112 that have not been developed, as well as objects which do not have actual heat demand, like landfills or substations. To exclude such areas from the assessment, we added the OpenStreetMap (OSM) buildings dataset to the model [24]. During the modeling process, we only considered the DLM objects which intersect with the OSM buildings dataset. We then created new area objects by calculating minimum boundary polygons around the buildings inside the objects classified as 2112 and 2114 and labeled them as “TCS area objects.” Analogous to the residential sector, we then assessed an area factor for those TCS areas and assigned them their object-specific heat demand (GWhth/a) and heat density (kWhth/m2a).

The quality of the DLM dataset varies strongly for objects within object type 2112. Only the federal states of Brandenburg, Berlin, Saarland and Hamburg provide detailed information about how an area is actually used, e.g., for industrial purposes or for commercial functions like trade or storage. We were only able to distinguish between the heat demand for the industrial sector and the TCS sector in these states. Since no distinction could be made in the other states, we used the combined heat demand of the industrial and the TCS sector to model the heat demand of TCS areas.

Large individual heat consumers

We identified four categories of LIHC which are particularly suited for DHNs fueled by heat generated from biomass: schools, hospitals, public outdoor swimming pools, and commercial greenhouses [26]. For the categories “schools” and “hospitals,” we used the POI-Bund Dataset which contains both spatial and content-related information on schools and hospitals such as the number of students per school or the number of beds per hospital [23].

For the category “public outdoor swimming pools,” we used the DLM which contains outdoor swimming pools (object type 2345). The values for heat demand per square meter of water surface vary in the literature between 280 kWhth/m2a and 700 kWhth/m2a [36, 37]. In order to define an appropriate value for our model, we evaluated municipal energy reports on ten public outdoor swimming pools and used the mean value of 435Footnote 6 kWhth/m2a.

The spatial information for the category “greenhouses” has been extracted from the OpenStreetMap buildings dataset. As this dataset does not contain any information on whether a greenhouse is used commercially, we have set the threshold of 100 m2 floor spaceFootnote 7 as an indicator for this.

The heat demand for each object of those four categories was determined on the basis of object-specific key figures and object-specific reference units. Those key figures and their respective reference units are shown in Table 2. This enabled us to model the heat demand for a total of 55,465 individual objects.

Table 2 Heat demand key figures of large individual heat consumers

Assessment of the BEP-specific heat sales potential

To determine the BEP-specific heat supply potentials, we looked at two different options on how BEPs can feed heat into a potential DHN. The spatial and technical information needed to do this was taken from the DBFZ bioenergy plant database. This database is continuously updated and contains detailed information on BEPs in Germany. It is based on information which the Federal Network Agency (Bundesnetzagentur, BNetzA) publishes on BEPs as part of the German Renewable Energy Act (Erneuerbare Energien Gesetz (EEG)) remuneration scheme and explained in more detail in Scheftelowitz et al. [18]. It is not possible to derive from the data whether a BEP already feeds into a DHN or not. For the following analyses, this aspect is not taken into consideration. It is also not possible to determine from the BEP database whether a BEP is actually a production unit or if it is a satellite CHPU which cannot be operated independently of its production unit. For this reason, all investigated BEPs are considered to be production units, an assumption which distorts the actual situation.

Distance-related heat supply scenarios

To determine the plant-specific tHSP, we looked at two different design options for DHNs: A and B. These are defined as follows: in Option A the heat produced at the BEP site is transported to the consumers via a directly connected local heating network, whereas in Option B biogas is transported via a raw biogas pipeline to a CHPU located directly on the site of the heat sink and the DHN commences at this point. Since the way heat is transported in Option A involves losses of up to 10%, the heat sink must be located in the immediate vicinity of the BEP [48]. The minimum overall heat density for a DHN in year-round operation is defined by the Quality Management (QM) system of the German/Swiss initiative “QM Holzheizwerke” as being 2 MWh/line meter and year [48]. In a comparable calculation for a region in Rhineland-Palatinate made in 2010, IZES and the Institut für Energie- und Umweltforschung (ifeu) [49] stated that there should be a heat sink within a radius of 1000 m around a BEP to reach this value. Daniel et al. [50] found a higher figure for the whole grid and concluded that losses must not exceed 20% to ensure sufficient profitability. Losses up to this value occur with total grid lengths of between 1000 m and 5000 m and depend on the capacity of the respective CHPU [50]. Taking this assumption into account, we defined for our model a maximum distance of 1500 m between the BEP and the heat sink. This option is applicable for all BEPs.

Option B is applicable for biogas plants only. Since biogas can be transported over longer distances with only minor losses, the distance between the heat source and heat sink can be greater for this option. Obviously, this distance depends again on the capacity of the installation. IZES and ifeu [49] reported that BEPs with an installed capacity of 500 kWel can reach a distance of up to 5000 m for raw biogas pipelines without losing their economic feasibility [49]. In our case, we followed the assumption of this maximum distance between BEP and heat sink (5000 m) even though the mean value of the installed capacity of biogas plants throughout Germany is only 435 kWel [51]. A visual illustration of the two options for district heating networks under consideration is provided in Fig. 2.

Fig. 2
figure 2

District network options for bioenergy plants

To identify the plant-specific heat sales potential, we then used the spatial information on BEPs from the DBFZ bioenergy plant database and conducted a spatial query into whether there are suitable heat sinks at the respective distances. For residential areas and TCS areas, we defined that the minimum heat densityFootnote 8 of the nearest heat sink has to be higher than 50 kWhth/m2 × aFootnote 9 in order to be suitable for the DHN design options presented above. For the LIHCs, we did not define such a precondition, but instead considered the total heat demand of the respective LIHC.

Categorizing the heat supply potential using line density modeling

So far in the spatial modeling approach, all heat sinks with a corresponding heat density of 50 kWhth/m2a in the vicinity of the respective BEP are initially selected as potential heat sinks. However, the gross heat demand (in GWhth/a) of the heat sink can be insufficient for achieving an economically viable DHN. Therefore, the spatial model required enhancement. To do this, we introduced the line density factor LDF. This factor represents the ratio of the heat demand of a heat sink in relation to the total potential network length and the funding conditions of the German Kreditanstalt für Wiederaufbau (KfW) for DHNs in GermanyFootnote 10 and was calculated as follows:

$$ LDF=\left(\frac{Q_{RAO}}{D\times NLF}\right)\times 500\ kWh/a $$

where LDF is the line density factor, Q is the heat demand (kWhth/a), D is the linear distance bioenergy plant-heat sink (m), NLF is the network length factor, and RAO is the residential area object.

We assumed that if the annual heat demand of a given heat sink was unable to generate the required heat turnover of 500 kWhth per line meter [53], it was highly unlikely that a DHN would be built. This much more simplified method is a suitable way to eliminate inadequate heat sinks from the model. Values for a line density factor below 1 indicate that the annual heat turnover of a potential DHN would be too low. All heat sinks with a line density factor < 1 were therefore excluded from the model, which subsequently only contained heat sinks suitable for DHNs.

In order to determine the line density factor, it is necessary to know the total length of the DHN. Within the GIS model, however, only linear distances from the heat source to the potential first consumer are known. The total network length must be estimated. For this purpose, we introduced the network length factor NLF. It describes the ratio of the total network length of a DHN to the distance between the heat source (the BEP) and the initial heat consumer and was calculated using Formula 7.

$$ NLF=\frac{L_{tot.}}{D} $$

where NLF is the network length factor, Ltot.. is the total network length (m), and D is the linear distance bioenergy plant to heat sink (m).

The NLF had to be assessed empirically. For this purpose, we analyzed 20 existing biomass-fueled DHNs for which the total network length is known [54,55,56,57,58,59,60,61,62,63,64,65,66,67]. We measured the distance between the heat-generating plant and the potential first consumer of these DHNs and calculated the NLF. Our data on the respective DHNs only includes cases where the distance between the plant and the first consumer is less than 500 m. For these systems, the median value of the NLF is 26.1.Footnote 11 For systems where there is a distance > 500 m between the heat source and the initial consumer, which occurs in our GIS data, the NLF was determined with the corresponding degression (see Fig. 3).

Fig. 3
figure 3

Network length factor for twenty biomass-fueled district heating networks

Figure 4 illustrates the correction process described above. If heat sink 1 is located within the search radius of 1500 m around a BEP and has a heat demand of 1.9 GWhth/a then, to develop the area, the heating network would have to be 6500 m long (assumption). However, to receive financial support, an annual heat turnover of at least 3.25 GWhth (6500 m × 500 kWhth/a = 3,250,000 kWhth/m×a) is required. Heat sink 1 does not have a sufficient heat demand to operate the DHN economically and can therefore not be considered suitable for a local heating network (line density factor < 1). Heat sinks 2 and 3, on the other hand, have a sufficient heat demand and can therefore be considered for a DHN (line density factor > 1).

Fig. 4
figure 4

Illustration of the model correction using the line density factor

Economic viability of heat sales potential

This section presents the methodology used for assessing the economic heat sales potential (eHSP), since a second objective of this study was to evaluate the economic viability of the specific heat supply potential of the current German BEP stock. Here, the same dataset for BEPs was used as shown in Table 1.

Even though this study aimed to gain initial insights into the market potential of heat as a by-product of electricity generation, a full techno-economic feasibility assessment was not within the scope of our study. In our approach, economic viability depends on the potential of heat sales (demand/supply) and the associated estimated costs for heat distribution infrastructure (pipelines, peak load boilers, and storage facilities) as well as the earnings from heat sales.

Evaluation scheme and main assumptions

In order to assess the eHSP, we developed an evaluation scheme that consisted of four main steps (Fig. 5). We first established individual BEP-heat sink relationships (Step 1) and then derived relevant indicators based on plant-specific data (Step 2). Before transferring all the information onto a bigger scale (Step 3), we grouped all the individual BEPs into clusters according to their EEG remuneration key. This enabled us to obtain the desired overview of the German BEP stock without compromising plant-specific data and provided the basis for the final economic assessment (Step 4).

Fig. 5
figure 5

Scheme to evaluate economic heat sales potential of BEPs in Germany, from prioritization (1) to derivation of plant-specific indicators (2), their transfer to a bigger cluster scale (3), and their economic evaluation (4)

In addition to the information on location and the related potential heat demand from spatial heat demand modeling (see above), we used energy-related BEP inventory data (Table 3) obtained from two databases: publicly available data published by the German Federal Network Agency (BNetzA) and datasets retrieved from annual plant operator surveys by the DBFZ (see [51]).

Table 3 Mean cluster-specific energy (supply) capacities based on BEP-specific values from German inventory data

In this approach, we assumed four factors to be determinant for the economic viability of heat sales from BEPs (Table 4). Whereas the distance between the heat sink and the BEP (1), the type of heat sink (2) and the heat supply-demand ratio (3) all govern the choice and size of the heat distribution infrastructure (pipelines, peak-load boilers, and storage facilities), the heat market price (4) ultimately determines income from heat sales. Other factors (e.g., local conditions underground which influence the excavation costs or regional differences between temperature profiles) may also have an influence on economic viability but their consideration would go far beyond the scope of this study.

Table 4 Key determinants for assessing the economic feasibility of BEP heat sales potential (eHSP)

Steps 1-3: assessment indicators

In the first step of this evaluation scheme (Fig. 5), we used the spatial and energy-related information (derived from the above presented assessment of tHSP) to establish single BEP-heat sink relationships. We assumed specific heat delivery costs would decrease (economies of scale) as potential heat demand increased. Therefore, we linked each BEP to a single heat sink within a radius of 1500 m (Option A of the DHN design options) that had the highest potential heat demand (see Fig. 4 for an illustration). Since we expected that heat distribution costs would significantly affect economic viability, we excluded the option of exploring a larger radius or possible connections between multiple producers and consumers We deliberately excluded biogas plants with potential heat sinks within a radius of 1500-5000 m and the construction of a raw biogas pipeline from this assessment (Option B of the previously described DHN design options). Also, our methodology is only suitable for a 1:1 ratio between BEP and heat sinks, so we were unable to incorporate cross-linked networks.

We operationalized the BEP-heat sink relationship by introducing the “Q-factor” (FQ) which is heat demand (Qsnk) divided by net heat production (Qsrc) (Formula 8). Furthermore, we decided to use this indicator as a main reference point for the economic assessment.

$$ \mathrm{FQ}=\frac{Q_{\mathrm{snk}}}{Q_{\mathrm{src}}} $$

where FQ is the heat source-sink-ratio, Qsnk is the heat demand of a heat sink (GWhth/a), and Qsrc is the net heat production of BEP (GWhth/a).

To refine this assessment, we considered seasonal fluctuation in heat demand depending on the characteristics of the type of heat sink (TyHS). Here, we distinguished between three TyHS for which we created synthetic heat load profiles according to a methodology that was presented by [68] (Table 5). For reasons of simplicity, we allocated the annual heat demand over the course of a year (8760 h) and decided not to create regionally explicit heat load profiles but to use only one climate reference dataset (weather station in Potsdam [69]).

Table 5 Construction of synthetic heat load profiles

In the second step of the evaluation scheme (Fig. 5), we extended the plant-specific data and derived other indicators that influence the economic viability of heat sales, i.e., the maximum cogeneration utilization rate (CURmax) and a heat pipe capacity factor (Ftk).

The maximum cogeneration utilization rate (CURmax) indicates how much of the BEP’s gross heat generation could be sold on the market, taking into account a seasonal fluctuation in heat demand. For reasons of simplicity, we do not consider flexible operation of BEPs but assume 24-h operations that result in a constant heat supply at full capacity (net heat production) throughout the year. When both power and heat have to be considered at a temporal resolution, the analysis becomes much more complex. We focused on varying heat demands and used heat load profiles, we had generated previously (Table 5) to calculate the CURmax for each TyHS and for the FQ at intervals of 0.1 until a CURmax of nearly 100% was reached (see Fig. 6).

Fig. 6
figure 6

Maximum cogeneration utilization rate (CURmax) in relation to FQ for three types of heat sinks (TyHS)

Furthermore, we assumed that the investment costs of heat pipelines depend not only on their length but also on their diameter (heat transmission capacity). When heat storage facilities (HST) and bio-fueled peak load boilers (PLB) were used to supplement cogeneration and maintain heat supply, we also needed to adjust the heat pipeline diameter accordingly. Therefore, we introduced a heat pipe capacity factor (Ftk), defined as the ratio of net heat supply to annual peak loads as observed in the generated load profiles (Table 5).

In the third step of the evaluation scheme (Fig. 5), we transferred the single-source sink-related data to the cluster level by aggregating the previously gathered information (FQ, CURmax, Ftk) and the cluster-related information (Table 3). To do this, we divided each cluster into five groups each of which represented 20% of the BEPs in the respective cluster (quintiles) and which referred to equally sized classes of FQ. The FQ values assigned to each previous BEP-heat sink pair (Formula 8) were used for this statistical division. Whereas quintile 1 represents the 20% of BEPs with the smallest FQ values, quintile 5 represents the 20% of BEPs with the largest FQ values. Using this approach, we then assigned each BEP-heat sink pair to a cluster quintile and calculated the arithmetic mean of the CURmax and Ftk.

Step 4: economic assessment

In the fourth and final step of the evaluation scheme (Fig. 5), we assessed the economic heat sales potential (eHSP) of German BEPs for all cluster quintiles based on the previously generated dataset. The remainder of this section will describe how we calculated the capital expenditures (CAPEX) for the heat distribution infrastructure and how we related these costs to expected earnings from heat sales.

The CAPEX in this assessment included investment costs such as annuities for heat pipelines, peak load boilers (PLB), and heat storage facilities (HST). For all cost calculations, we fixed the interest rate at 5% and assumed no price increase during a depreciation period of 15 years for PLBs and HSTs, and 25 years for heat pipelines. As a requirement for heat sales from BEPs, we first calculated the CAPEX for heat pipelines for all BEP-heat sink pairs. Then, we defined the PLBs and HSTs to those pairs that fulfilled a specific set of criteria (Table 6) and calculated the additional CAPEX from using them (including fuel costs for PLB). Finally, we calculated the mean CAPEX for each BEP cluster quintile and added the expected earnings from heat sales. Hence, we determined the economic viability of heat sales from German BEPs, including additional earnings from using PLBs and HSTs.

Table 6 Rules for deciding on the heat sales marketing concept based on the FQ

To calculate the CAPEX for heat pipelines (Cpipe), we derived specific pipeline costs and multiplied them by the length of the heat pipelines required to connect each BEP with a heat sink (Formula 9). We introduced a correction factor (CF) to account for the fact that a heat pipeline will need to deviate from a direct line in order to overcome infrastructural barriers like residential areas, roads, or rail lines. We multiplied this CF by the distance between each BEP and heat sink to derive the required pipeline length (Formula 9). However, we did not consider the local distribution and transmission infrastructure that would be required, e.g., for district heating in residential areas [70].

$$ {L}_{\mathrm{T}}=D\times CF $$

where LT is the length of a heat pipeline that is equal to the distance between a BEP-heat sink pair (Step 1, Fig. 5) multiplied by a correction factor (CF). CF is equal to the square root of 2 which results from applying the Pythagorean theorem. Streamlining deviations, we divided the alternative way (the sum of the two legs of a right-angled triangle) by the direct distance between the BEP and heat sink (hypotenuse).

$$ {C}_{\mathrm{pipe}}=\left({C}_{\mathrm{T}}+{C}_{\mathrm{c}}\right)\times {L}_{\mathrm{T}} $$

where Cpipe is the cost of the heat pipeline calculated for each BEP-heat sink pair (Step 1, Fig. 5). CT is the cost of the heat pipelines per meter (as a function of Pth and Ftk). Cc represents other construction-related costs per meter (underground work and planning), and LT is the length of heat pipeline (m).

Serving peak loads or storing heat in times of low demand might be economically reasonable if the expected additional earnings outweigh the additional CAPEX for HSTs and PLBs (including fuel costs for a PLB). Furthermore, the use of HSTs and PLBs would have an impact on the CAPEX for heat pipelines, affecting the required heat transmission capacity (diameter). We therefore established a set of rules for three different heat marketing concepts based on the level of FQ (Table 6). According to these rules, we decided whether to use peak load boilers (PLB) and/or heat storage (HST) for each BEP-heat sink pair and adjusted the required heat pipeline capacity accordingly.

Following a “full supply” marketing concept, we used HST and PLB and placed them close to the source (here the BEP). By adapting the pipeline capacity, all peak loads of heat could be potentially served via the heat pipeline (Ftk×Pth). The “basic supply” marketing concept differs as it prescribes only the use of an HST. Due to periods of peak demand, heat pipelines should be able to transport heat produced by the BEP and discharge the capacity of the storage. We therefore proposed that a heat pipeline’s capacity should be double the value of the BEP’s thermal capacity (Ftk = 2) while both the charge and discharge rates are equal to the BEP heat capacity. The “full feed-in” option assumes there is no need for a PLB or HST so that we based the transmission capacity of the pipelines solely to the BEP’s capacity (Ftk = 1).

In addition to considering the use of PLB and HST for the “full supply” and “basic supply” concepts, we optimized their size aiming to find the most cost-efficient combination of both. For this purpose, we used a cost function (Formula 12) that builds upon an optimization approach which suggests that the investment costs for HST and PLB depend on their volume and power capacity [71].

The optimization approach for an ideal combination of PLB and HST is described in [71]. Here, we only accounted for the different heat load profiles for each quintile in every cluster (based on the previous TyHS and CURmax) in order to find minimum overall costs.

$$ {C}_{\mathrm{opt}.}=f(V)+f(P)+\left({C}_{\mathrm{fuel}}-{C}_{\mathrm{h}@\mathrm{POS}}\right)\times {P}_{th}\to \mathit{\min} $$

where Copt. is the cost optimum, f(V) and f(P) are functions of the CAPEX for HSTs and PLBs [71] and Pth is a cluster’s power capacity (see Table 3), Cfuel represents fuel costs and Ch@POS is the heat price at the point of sale. For simplification purposes, we fixed the potential heat losses from heat storage at a value of 10%, the PLB efficiency (η) at 85%, the Ch@POS at 50 €/MWhth and the Cfuel at 30 €/MWhth.

We are aware of the fact that we created an inherent incentive to use PLB (Ch@POS > Cfuel) and real prices for heat would be influenced by many factors such as the locality, market characteristics, and the type of heat sink. For instance, factors that affect local demand for heat from BEPs include existing access to the gas grid, availability (and abundance) of solid biogenic fuel (from wood), and a fossil-fueled heat source that is already connected to existing heat pipelines. However, due to the scope of this study, we did not specify and adjust the calculations to these factors. Finally, we determined the economic heat sales potential (eHSP) from German BEPs for all the cluster quintiles by calculating the annual heat sales potential in terms of heat earnings (Esale) (Formula 13) minus the financial expenditure for heat distribution (PLB, HST, and heat pipelines) (Formula 14).

$$ {E}_{\mathrm{sale}}={C}_{\mathrm{h}@\mathrm{POS}}\times {W}_{\mathrm{th}}\times {\mathrm{CUR}}_{\mathrm{max}} $$

where Esale is the annual heat earnings (€), Ch@POS is the price of heat at the point of sale fixed at 50 €/MWhth, Wth is the net heat production (Table 3), and CURmax is the maximum cogeneration utilization rate for each BEP cluster quintile.

$$ \mathrm{eHSP}={E}_{\mathrm{sale}}-{C}_{\mathrm{opt}}-{C}_{\mathrm{pipe}} $$

where eHSP is the economically viable heat sales potential, Esale is the annual heat earnings (€), Copt. is the cost optimum, and Cpipe is the cost of heat pipeline from BEP to heat sink (€).

Since we considered mainly heat as a by-product of electricity generation, we divided the eHSP (Formula 14) by the current electricity generated by the German BEP stock (as kWh of electrical power per BEP cluster [51]). This enabled us to compare the option of heat sales to other CHP options.


The outcomes of this study can be divided into three categories. Firstly, the spatially explicit heat demands for different relevant sectors, secondly the heat sales potential of the BEPs as a whole and on a plant-by-plant basis, and finally the estimation of the valuable heat sales potential.

Spatially explicit data on heat demand

As a result of the previously described modeling process, spatially explicit data is now available on the heat demand of various types of heat sinks (residential, TCS and industrial, large individual heat consumers) at a sub-municipal level for the whole of Germany. The data includes the type of heat sink (residential area objects, TCS/industrial area objects, LIHC), their gross heat demand in GWhth/a and their heat density in kWhth/m2a. Table 7 lists the residential area objects and TCS/industrial area objects based on their heat density. Here, the data indicates a clear difference between the two types of heat sinks. TCS/industrial area objects show a generally higher heat density than residential area objects. This is particularly noticeable in the number of objects with a heat density higher than 50 kWhth/m2a (the previously defined minimum heat density for a DHN). The vast majority (91%) of the TCS/industrial area objects have a heat density above this value; 41% even exhibit values of more than 200 kWhth/m2a. For the residential area objects, on the other hand, only approximately 7% of all objects have a heat density of 50 kWhth/m2a or higher. Almost 75% exhibit heat density values of less than 25 kWhth/m2a.

Table 7 Number of residential area objects and TCS/industrial area objects classified by their heat density

A visual representation of this data for the categories residential sector and TCS/industrial sector is presented in Fig. 7. This section of a city depicts how the variations in urban fabric produce different values for heat density, making them either more or less suitable for district heating networks. The inner-city areas show higher heat densities (usually more densely built up, more inhabitants per area) whereas the outer areas exhibit lower values (often residential areas with a less dense development).

Fig. 7
figure 7

Example of the heat demand of residential area objects and TCS/industrial area objects modeled at a sub-municipal level

Technical heat sales potential

Total heat sales potential

Using the previously presented methods, we were able to identify those BEPs which have suitable heat sinks nearby and to quantify the tHSP per category and on the whole. Of the total 14,236 BEPs analyzed, 52% (7383) have tHSP. The category TCS/industrial is most dominant with 5875 BEPs located next to suitable heat sinks, which is about 80% of all plants with tHSP. Two thousand four hundred twenty-nine plants have tHSP in the residential sector and 2,096 plants could serve as heat suppliers for LIHC. The detailed figures are presented in Table 8.

Table 8 Number of BEP with option for district heating networks and heat sales potential per category

When considering the two different options for DHNs, Option A (DHN commencing directly at the BEP site) is suitable for most BEPs. For the categories TCS/industrial and LIHC, Option A is suitable for 85%, and 91%, of the plants respectively, while Option B is only suitable for 15% and 9% respectively. In contrast, for the residential sector, the ratio between Option A and Option B (DHN starting from the BEP, raw biogas pipeline to CHPU located at heat sink) is more evenly distributed. Forty-one percent of the BEPs with tHSP could theoretically apply Option B. This is conclusive as the previous chapter has shown that TCS/industrial areas are better suited for DHNs or better qualified as heat consumers, as heat demand and heat density are generally higher in these areas.

Looking at the heat sales potential per category, a similar picture unfolds. The overall tHSP is 153.61 TWhth/a of which more than two-thirds (69% or 105.94 TWhth/a) correspond to the category TCS/industrial. The residential sector accounts for 30% (46.15 TWhth/a), and 1% (1.52 TWhth/a) is attributable to LIHC.

The total tHSP of 153 TWhth/a corresponds to 17% of the final energy demand for heating (only space heating and hot water; reference year 2011, which is the point in time of the IZES heat demand figures).

For the residential and TCS/industrial sector categories, this results in a proportion of final energy consumption in the heating sector of 7% and 38% respectively.Footnote 12

However, this potential is not evenly distributed among all BEPs with a DHN option. Taking into account all BEPs for which we could determine sufficient tHSP, and using the heat density of the respective heat sinks as an indicator for the feasibility of a DHN, it is possible to identify statistical hot and cold spots. The hot spots represent plants whose heat sinks have heat densities that have a high statistical significance and are surrounded by plants whose heat sinks have equally high heat densities. For the cold spots, the opposite is true. Consequently, it is possible to identify spatial clusters of plants which could be preferable for a heating network. Figure 8 shows all BEPs with tHSP for each category as well as statistical hot and cold spots.

Fig. 8
figure 8

Statistical hot and cold spots of bioenergy plants with tHSP in relation to the heat density of their nearest heat sink

The identified hot spots point to conditions in the corresponding regions that tend to be advantageous for DHNs. In the residential sector, these include a high population density and corresponding housing typologies (i.e., multi-story buildings) [73]. For the category TCS/industrial, we observed hot and cold spots. The hot spots indicate that the BEPs are close to commercial or industrial areas with a high heat demand, whereas the cold spots imply that the heat demand in the areas under observation is low in comparison to surrounding areas.

Economic heat sales potential

As described in the methodology section, we analyzed the German BEP inventory data to estimate the economic heat sales potential (eHSP) in a step-by-step process (see Fig. 5). Since we limited the scope of this study to looking at the bigger picture, we did not conduct regionally explicit assessments. Therefore, we aggregated single plant data into a set of 20 BEP clusters and further refined the view by splitting each cluster into 20% shares (quintiles) based on a heat source-sink-ratio (FQ). This ratio represents the available heat generated by an operating BEP divided by the potential heat demand of a heat sink within a maximum distance of 1500 m.

The results for all clusters and their quintiles (Qn) are presented in Table 9 and show that for some of the quintiles the costs for heat distribution exceed the net earnings and vice versa. The results are quite heterogeneous with deficits reaching 12.8 eurocents per kWh of electricity produced by BEPs (Cluster 18, Q2) and potential earnings of up to 16 eurocents per kWh (Cluster 10, Q5). Heat sales could be economically viable for 78% of the assessed BEPs but would generate deficits in 13% of the cases. Nine percent of the BEP quintiles would not be considered for heat distribution since they only contain BEPs with no potential heat customer in their vicinity (1500 m). Furthermore, heat sales seem to be economically more viable for the upper quintiles in each cluster (highest FQ).

Table 9 Annual heat sales potential: net earnings minus annuity of distribution costs

One of the reasons why the eHSP in the upper quintiles is higher, it is due to the fact that these mostly contain BEP-heat sink-pairs with higher FQ values. Since the quintile distribution was based on the FQ, the proportion of salable heat from BEPs (waste heat) is therefore higher in the upper quintiles and affects the eHSP. Once the FQ exceeds a value of 5, we assumed the maximum cogeneration utilization rate (CURmax) to be 100% (Fig. 6). In Table 10, we can observe that the 4th and 5th quintiles show values from 92-100% for CURmax in all the clusters. This indicates that at least 40% of the BEPs within each cluster can possibly supply all, or nearly all, of their net heat production to nearby heat sinks.

Table 10 Results for the gross exhaust heat amount, the individual CURmax of the 20 BEP plant clusters

Table 9 and Table 10 also display the actual supply potential of the German BEP stock. The overall exhaust heat amount of all 20 clusters totals to a gross value of 37.7 TWhth/a (Table 9). Taking into consideration the different rations for CURmax, the externally usable heat amounts to 19.1 TWhth/a (Table 10) of which 16.3 TWhth/a could be identified economically viable heat sales potential (Table 9).


As the results indicate, the presented methods are a suitable approach for determining the BEP-specific tHSP and subsequently for estimating the eHSP. The required data could be generated by mostly publicy available data sources as described in the method section. The data indicates that there is currently unused potential especially with regard to tHSP but for eHSP too. Most of the BEPs display positive values concerning this matter, especially in the quintiles Q4 and Q5 of the 20 BEP-clusters.

The total tHSP of 153 TWhth/a corresponds to 17%Footnote 13 of the final energy demand for space heating [72]. Bioenergy plants in contrast, generate a gross heat amount of 37.7 TWhth/a, of which 19.1 TWhth/a could be used externally (net heat amount). Our results show, that the economically viable supply potential of 16.3 TWhth/a (see Table 9) almost covers the net heat generation.

According to most recent energy statistics, biogas plants already generate heat in an order of 15.5 TWhth/a [10]. We estimate that the economical heat sales potential for those clusters of plants would be about 11.2 TWh. Thus, assuming that co-generated heat is not solely used for space heating, the order of the calculated potential is closely related to the recent amount. The official energy statistics [74] do not contain disaggregated data for plants fueled by solid biomass, so a similar comparison is not possible here.

These numbers exceed the figures which have been found in energy scenarios for heat supply by local district heating networks in Germany [14]. However, it also displays the high potential of biomass fueled co-generation plants, which could also provide a noteworthy contribution to renewable heat in Germany.

To unlock this currently unused potential, further investigations on specific constraints are necessary, because the input dataset has some weaknesses. First, the heat demand was estimated using primary data representing the status quo in 2011. For modeling the heat demand, influencing factors such as regional climate differences, the refurbishment status of the buildings, or the impact of demographic change were not taken into account. In particular, the latter is likely to have an impact on settlement structures and, as energetic refurbishment progresses; the heat demand in the building sector will decrease. However, the stagnating energy demand for heating and hot water [74] shows that energetic refurbishment of buildings (existing stock and new constructions) is progressing slowly. The impact of this on the model results should be low. Second, data for the TCS/industrial sector does not include the sector-specific demand pattern (temperature level, daily and seasonal demand patterns). In order to assign the exact heat demand to each individual industrial or commercial area, it would also be necessary to know exactly what kind of companies (in terms of industry) are located in a specific area. Detailed feasibility studies, in which the heat demand of the consumers is precisely determined, are always necessary in order to make concrete investment decisions in heating networks. Furthermore, there are a number of hard and soft conditions that a municipality must fulfill in order to ensure that the prerequisites for the economic operation of a DHN are met. These include setting the right regulatory framework, e.g., securing areas in spatial planning or requiring connection to the DHN in certain areas [6].

Third, temperature level constraints, in terms of both the forward and backward flow of potential DHNs, were not considered during the modeling process. This is due to the fact that the DBFZ BEP plant database does not include technical specifications of the respective CHPUs or the temperature levels they are able to provide. Thus, it cannot be assumed that the heat demand of heat consumers requiring higher temperature levels (e.g., process heat for industrial processes) can be completely covered by the heat generated by the corresponding BEP. Nevertheless, the analyzed BEPs should at least be able to contribute to the generation of the required space heating or domestic hot water.

With these uncertainties in mind, the previously stated tHSP for BEP is over- rather than underestimated.

When discussing the results for estimating economic viability, some inadequateness of the methodology should first be clarified. It is important to mention that the described approach for calculating the heat sales potential for all BEP in Germany is based on spatial information, which shows some unavoidable blurring. For example, the location of biomethane CHPs and their corresponding heat sinks may produce results (Q1-Q3) which imply that a significant proportion of that cluster cannot ensure a CURmax of 100%. However, it is assumed that all biomethane-CHPs utilize all of their heat, since the remuneration scheme they use within the EEG demands complete cogeneration utilization. It is also important to keep in mind that the results presented here represent mean values calculated for the specific quintiles and therefore only describe the conditions for a group of BEPs. This does not necessarily mean that each individual plant in a cluster can be precisely characterized by the values for the quintiles. Also, the assumption of a direct link between the heat source and the heat sink underestimates real world conditions. Usually, there would be many barriers which lead to additional demands for a heat pipeline. Up to this point, we do not have a suitable approach for identifying empirical data for the extra distance needed. In order to interpret the tHSP results, it needs to be mentioned that we did not include competing heat providers in our analysis. This must be especially be taken into account for the TCS/industrial category where the calculated heat demand is likely to be met by other sources such as excess industrial heat.

The methodology was developed within the project “Bioenergy - potentials, long-term perspectives and strategies for power generation plants after 2020 (BE20plus)” which pursues the aim of evaluating the economic viability of BEPs after they drop out of the EEG’s first remuneration period after 20 years. Subsequently, we will use these results to establish which proportion of the recent BEP-portfolio can achieve profitability in a (potentially) second operational period. The findings show that heat sales are not a worthwhile option for some quintiles of the BEP clusters. This is partly due to the fact that there are no suitable heat sinks within the search radius (1500 m) and that the estimated costs for heat distribution exceed the assumed earnings from heat sales. This applies to cases where the heat demand of a heat sink is insufficient and the distance between BEP and heat sink is too far.

In general, even though decreasing electrical efficiency would increase eHSP, power plant operators would not compromise their electricity generation for greater (or exclusive) heat production, since the economic value of heatFootnote 14 is lower than that of electricity within the EEG remuneration scheme.


In our study, we analyzed the heat sales potential of BEPs within the EEG remuneration scheme in Germany. The results show that, based on our assumptions, 52% of the BEPs are suitable for distributing heat through a DHN. The overall technical heat sales potential for these BEPs is approximately 153 TWhth/a. About two-thirds of this is attributed to the TCS/industrial sector and about one-third to the residential sector. The share of large individual heat consumers in this potential is only about one percent. This total technical heat sales potential is equivalent to almost 12% of the final energy consumption in the German heating sector. However, this analyzed technical heat sales potential of 153 TWhth/a cannot be fully covered by the BEPs identified with the potential for a DHN option. Their gross supply potential for externally usable heat adds up to 19.1 TWhth/a and an economically viable supply potential of 16.3 TWhth/a. The figures on heat sales potential reveal two different aspects. On the one hand, it becomes clear that, in a best-case scenario, only half of the existing BEPs in Germany have the possibility of distributing heat by means of a DHN and generating revenues by doing so. This must be especially taken into account when looking at future possible business segments for BEPs after their EEG remuneration period ends, since those plants are not competitive by generating power only. On the other hand, we identified considerable potentials for DHNs for those BEPs which have heat sinks nearby. These networks can either be fully supplied with biomass or supplemented with other technologies suited for grid-bound heat supply (e.g., larger heat pumps). We conclude that biomass has further potential as a renewable fuel in supplying district heating networks and that the current BEP stock can contribute to increasing the share of RES in the heating sector. Decisive factors in this context are the BEP’s relationship to potential heat customers and their specific requirements for heat supply. In order to fully exploit this potential, the BEPs have to be embedded in local or regional heating strategies which take local conditions into account.

With regard to economic viability, it can be concluded that even if there are potential heat consumers within an acceptable distance to the BEP, heat sales are often not an attractive opportunity for all of them since the costs for distributing the heat can exceed the potential earnings. This might change in the medium term with the obligation to include RES into new heating systems, the introduction of CO2-prices and the ban of oil-fired heating systems from 2026 on [76]. These measures are included in a law passed by the German Federal Government at the end of 2019 [77] to achieve its Paris Climate Agreement goals.

The question of heat sales potential is crucial for future business concepts of BEPs, since heat is the most valuable by-product for most of those plants. We are looking forward to integrating the results of this work in a wider approach to cover all factors relevant for BEP profitability.

Availability of data and materials

The datasets generated and/or analyzed as part of the current study are not publicly available due to the privacy policy of the institutions involved, but are available from the corresponding author upon request.


  1. Without rural local district heating networks fueled by heat generated from biomass.

  2. To protect privacy, grid cells with one person are identified as grid cells without persons and grid cells with only two persons are identified as grid cells with three persons [29].

  3. 75% of the residential area objects are covered by grid cells with more than one inhabitant.

  4. The manufacturing sector includes the economic sectors classified by WZ08-10 to WZ08-33 according to the German Classification of Economic Activities, 2008 edition [30].

  5. Meant here are all economic sectors except the manufacturing sector.

  6. For the original values see [38,39,40,41,42,43,44,45]

  7. We assumed that the floor space of privately used greenhouses does not exceed 100 m2.

  8. Heat density is a function of the area (m2) of a heat sink and its annual heat demand in kWhth.

  9. 50 kWh/m2a can be regarded as the minimum heat density required for a local heating network [52].

  10. The German Kreditanstalt für Wiederaufbau (KfW) is Germany’s largest development bank and grants loans to construct DHNs. In return, the bank demands a minimum yearly heat turnover of 500 kWhth/a per line meter [53].

  11. The original values for the total network length, the linear distance from the BEP to the potential first heat customer and the resulting NLF can be found in Table 11 in the Appendix.

  12. The final energy consumption in the heating sector for space heating and hot water (2011) was 916.1 TWhth/a of which 631.6 TWhth/a were consumed by the residential sector and 281.4 TWhth/a by the TCS and industrial sectors [72].

  13. As previously mentioned, this is the figure for 2011. Latest energy statistics (2018) state an energy demand for space heating and hot water of 764.4 TWhth/a [74]. In this case tHSP corresponds to 20%.

  14. Biogas plant operators received 2.6 ct/kWhth on average [75].



Amtliches Topographisch-Kartographisches Informations System


Bioenergy plant


Combined heat and power


Combined heat and power unit

CURmax :

Maximum cogeneration utilization rate


District heating network


Digital landscape model


Economical heat sales potential


Heat demand-supply factor


Heat storage




Population factor


Heat price at the point of sale


Peak load boiler


Annual heat demand of a heat sink


Q-factor (ratio of heat demand and gross heat supply of a given source)


Residential area factor


Residential area object


Trade, commerce and services


Technical heat sales potential


Type of heat sink


  1. UBA [German Environment Agency] (2018) Energieverbrauch für fossile und erneuerbare Wärme.

  2. European Commission [EC] (2016) Communication from the Commission to the European Parliament, the Council, The European Economic and Social Committee and the Committee of the Regions. An EU Strategy on Heating and Cooling. COM(2016) 51 final on 16.2.2016. EC, Brussels

  3. Fleiter T, Elsland R, Rehfeld M, Steinbach J, Reiter U, Catenazzi G, jakob M, Rutten C, Harmsen R, Dittmann F, Rivière P, Stabat P (2017) Heat roadmap Europe. A low-carbon heating and cooling strategy. Profile of heating and cooling demand in 2015. Fraunhofer Institut für System- und Innovationsforschung (ISI), TEP Energy GmbH, University Utrecht, ARMINES, Karlsruhe

  4. AGFW [Der Effizienzverband für Wärme Kälte und KWK e.V] (2018) AGFW - Hauptbericht 2017. AGFW, Frankfurt am Main

  5. Gerhardt N, Sandau F, Becker S, Scholz A, Schumacher P, Schmidt D (2017) Wärmewende 2030. Schlüsseltechnologien zur Erreichung der mittel- und langfristigen Klimaschutzziele im Gebäudesektor. Studie im Auftrag von Agora Energiewende. Fraunhofer-Institut für Windenergie und Energiesystemtechnik (IWES), Fraunhofer-Institut für Bauphysik (IBP), Berlin

  6. Baur F, Noll F, Vogler C, Wern B, Mees M, Steinert M, Wiedemeyer N, Beyer S, Hill A, Gebauer M, Pasuchinger T, Kornmann A, Menke C, Barton M (2016) Wärmestudie Region Eifel und Trier. Studie im Auftrag des Ministeriums für Wirtschaft, Klimaschutz, Energie und Landesplanung des Landes Rheinland-Pfalz. IZES gGmbH, SWT Stadtwerke Trier AöR, Energieagentur Region Trier GmbH, Kompetenzzentrum Solar der Hochschule Trier, ECOSCOP Gesellschaft für Umweltberatung und -recherche mbH Trier, Solites -Steinbeis Forschungsinstitut für solare und zukunftsfähige thermische Energiesysteme (Zuständigkeit Wärmenetze), Saarbrücken

  7. BMWi [Federal Ministry for Economic Affairs and Energy] (2019) Kohleausstieg und Strukturwandel. BMWi. Accessed 22.09.2019

  8. Fritz S, Pehnt M (2019) Der Kohleausstieg und die Auswirkungen auf die betroffenen Wärmenetze. Kurzstudie im Auftrag des Bundesministeriums für Umwelt. Naturschutz und nukleare Sicherheit. ifeu, Heidelberg

    Google Scholar 

  9. BMWi [Federal Ministry for Economic Affairs and Energy] (2018) Erneuerbare Energien in Zahlen. Nationale und internationale Entwicklung im Jahr 2017. BMWi, Berlin

  10. Zeitreihen zur Entwicklung der erneuerbaren Energien in Deutschland unter Verwendung von Daten der Arbeitsgruppe Erneuerbare Energien-Statistik (AGEE-Stat). (2019);jsessionid=906F1084ED959B25AC4C4B27C0C43A70?__blob=publicationFile&v=22.

  11. Hoffstede U, Stelzer M, Hahn H, Beil M, Krautkremer B, Kasten J, Beyrich W, Holzhammer U (2018) Vorbereitung und Begleitung bei der Erstellung eines Erfahrungsberichts gemäß § 97 Erneuerbare-Energien-Gesetz. Teilvorhaben II a: Biomasse. Zwischenbericht. Erstellt im Auftrag des Bundesministeriums für Wirtschaft und Energie. Fraunhofer-Institut für Energiewirtschaft und Energiesystemtechnik (IEE), Kassel

  12. Lund R, Østergaard DS, Yang X, Mathiesen BV (2017) Comparison of low-temperature district heating concepts in a long-term energy system perspective. International Journal of Sustainable Energy Planning and Management 12:5-18.

  13. FNR [Fachagentur Nachwachsende Rohstoffe e.V.] (2019) Biogas: Faustzahlen. FNR. Accessed 19.08.2019

  14. Koch M, Hennenberg K, Hünecke K, Haller M, Hesse T (2018) Rolle der Bioenergie im Strom- und Wärmemarkt bis 2050 unter Einbeziehung des zukünftigen Gebäudebestandes. Öko-Institut e.V, Freiburg, Darmstadt

    Google Scholar 

  15. Weinand JM, McKenna R, Fichtner W (2019) Developing a municipality typology for modelling decentralised energy systems. Util Policy 57:75–96.

    Article  Google Scholar 

  16. Esch T, Taubenböck H, Geiß C, Schillings C, Nast M, Metz A, Heldens W, Keil M (2011) Potenzialanalyse zum Aufbau von Wärmenetzen unter Auswertung siedlungsstruktureller Merkmale. Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR), Deutsches Fernerkundungsdatenzentrum (DFD), Oberpfaffenhofen

  17. IZES gGmbH [Institut für ZukunftsEnergie und Stoffstromsysteme] (2017) Wärmekataster Saarland. Abschlussbericht im Auftrag des MWAEV des Saarlandes. IZES gGmbH, Saarbrücken

  18. Scheftelowitz M, Becker R, Thrän D (2018) Improved power provision from biomass: a retrospective on the impacts of German energy policy. Biomass and Bioenergy 111:1-12. doi:

    Article  Google Scholar 

  19. BKG [Federal Agency for Cartography and Geodesy] Digitale Landschaftsmodelle und Landbedeckungsmodelle von Deutschland. Digitales Basis-Landschaftsmodell (Basis-DLM). Accessed 21 Aug 2019

  20. Baur F, Noll F, Wern B, Vogler C, Weiler K, Arnold K, Carpantier R, Pietzner K, Samadi S, Wane S, Hiebel M, Dresen B, Nühlen J (2015) Nachhaltige Integration von Bioenergiesystemen im Kontext einer kommunalen Entscheidungsfindung (KomInteg). Abschlussbericht. Institut für ZukunftsEnergie und Stoffstromsysteme (IZES gGmbH), Wuppertal Institut für Klima, Umwelt und Energie GmbH, Fraunhofer UMSICHT, Saarbrücken, Wuppertal, Oberhausen

  21. Verwaltungsgebiete 1:250 000 mit Einwohnerzahlen (Ebenen), Stand 31.12. (VG250-EW 31.12.) (2018) BKG [Federal Agency for Cartography and Geodesy]. Accessed 10.06.2018

  22. Geographical Grid of Germany in Lambert-Projection (2015) BKG [Federal Agency for Cartography and Geodesy].

  23. Points of Interest (POI-Bund/POI-BPol) (2018) BKG [Federal Agency for Cartography and Geodesy]. Accessed 21 Aug 2019

  24. OpenStreetMap data for Germany - Buildings (2018) Geofrabrik GmbH and OpenStreetMap Contributors.

  25. Ergebnisse des Zensus 2011: Bevölkerung im 100 Meter-Gitter (2015) Destatis [Statistisches Bundesamt].

  26. Gaderer M, Lautenbach M, Fischer T, Ebertsch G (2007) Wärmenutzung bei kleinen landwirtschaftlichen Biogasanlagen. Bayerisches Landesamt für Umwelt (LfU), Bayerisches Zentrum für Angewandte Energieforschung e. V. (ZAE Bayern), Augsburg

  27. Ergebnisse des Zensus 2011: Bevölkerung nach Jahren und Geschlecht für Gemeinden. (2014) Destatis [Statistisches Bundesamt].

  28. IWU [Institut für Wohnen und Umwelt] (2003) Deutsche Gebäudetypologie: Systematik und Datensätze. IWU, Darmstadt

  29. Zensusdatenbank - Einwohnerzahl je Hektar. Ergebnisse des Zenus am 09. Mai 2011 pro ha (2015) Destatis [Statistisches Bundesamt].,,,.

  30. Destatis [Statistisches Bundesamt] (2008) Klasifikation der Wirtschaftszweige. Mit Erläuterungen. 2008. Destatis, Wiesbaden

  31. Destatis [Statistisches Bundesamt] (2008) Bruttoinlandsprodukt/Bruttowertschöpfung – Jahressumme – regionale Tiefe: Kreise und krfr. Städte. Regionaldatenbank Deutschland der Statistischen Ämter des Bundes und der Länder. Destatis, Wiesbaden

  32. Destatis [Statistisches Bundesamt] (2011) Sozialversicherungspflichtig Beschäftigte am Arbeitsort nach Geschlecht und Nationalität – Stichtag 30.06. – regionale Tiefe: Gemeinden, Samt-/Verbandsgemeinden. Regionaldatenbank Deutschland der Statistischen Ämter des Bundes und der Länder. Destatis, Wiesbaden

  33. Destatis [Statistisches Bundesamt] (2010) Bruttoinlandsprodukt/Bruttowertschöpfung – Jahressumme – regionale Tiefe: Kreise und krfr. Städte. Regionaldatenbank Deutschland der Statistischen Ämter des Bundes und der Länder. Destatis, Wiesbaden

  34. Almut Kirchner, Schlesinger M, Weinmann B, Hofer P, Rits V, Wünsch M, Koepp M, Kemper L, Zweers U, Straßburg S, Matthes F, Busche J, Graichen V, Zimmer W, Hermann H, Penninger G, Mohr L, Ziesing H-J (2009) Modell Deutschland. Klimaschutz bis 2050. Vom Ziel her denken. Endbericht. Im Auftrag des WWF Deutschland. Prognos AG, Öko-Institut e.V., Basel, Berlin

  35. Destatis [Statistisches Bundesamt] (2011) Verarbeitendes Gewerbe, Bergbau u. Gewinnung von Steinen und Erden: Betriebe u. Beschäftigte nach Wirtschaftsabteilungen – Stichtag 30.09. – regionale Tiefe: Kreise u. krfr. Städte. Regionaldatenbank Deutschland der Statistischen Ämter des Bundes und der Länder. Destatis, Wiesbaden

  36. VDI [Verband Deutscher Ingenieure] (2014) Richtlinie VDI 3807: Verbrauchskennwerte für Gebäude. Blatt 2: - Verbrauchskennwerte für Heizenergie, Strom und Wasser. vol 3807. VDI, Düsseldorf

  37. Schlesiger G (2001) Energie- und wassersparende Maßnahmen in Bädern. 2001 edn. BISp [Federal Institute for Sports Science], Bonn

  38. Merker H-J, Nißen S (2006) Energiebericht der Stadt Voerde für das Jahr 2006. Kommunalbetrieb Voerde - Gebäudemanagement, Voerde

    Google Scholar 

  39. Stadt Dinslaken (2014) Energiebericht. Stadt Dinslaken Dinslaken

    Google Scholar 

  40. Kerpern S (2012) Energiebericht 2012. Stadt Kerpern, Kerpen

    Google Scholar 

  41. Stadt Rheinfelden (Baden) (2016) Energiebericht 2016. Stadt Rheinfelden (Baden), Rheinfelden (Baden)

  42. Stadt Wiesloch (2005) Energiebericht. Kommunale Liegenschaften. Berichtsjahr 2005. Stadt Wiesloch, Wiesloch

  43. Stadtverwaltung Paffenhofen a. d. Ilm (2014) Energiebericht 2014. Stadtverwaltung Paffenhofen a. d. Ilm, Paffenhofen a. d. Ilm

  44. Oertli Rohleder Wärmetechnik GmbH (2010) In Sachen Technik. Planungsanleitung Gas-Schwimmbadheizer GX 01 150 – 800. Oertli Rohleder Wärmetechnik GmbH, Möglingen

  45. Innovations- und Transferinstitut Bingen GmbH (2002) Wärmeversorgung für das Freibad und die Sporthalle in Maikammer. Forstliche Versuchs- und Forschungsanstalt Bingen

  46. Blesl M, Kempe S, Ohl M, Fahl U, König A, Jenssen T, Eltrop L (2009) Wärmeatlas Baden-Württemberg - Erstellung eines Leitfadens und Umsetzung für Modellregionen. Universität Stuttgart - Institut für Energiewirtschaft und Rationelle Energieanwendung (IER), Stuttgart

  47. Energieagentur Weststeiermark (2003) Energiebedarf der Gartenbaubetriebe – Auswertung der Gartenbaubetriebe. Zusammenfassung. Basierend auf der Projektstudie: Grundlagen für den Einsatz von erneuerbaren Energieträgern wie Biomasse und Solarenergie zur Beheizung und Trocknung im Gartenbau (Endbericht 31. Oktober 2003). Energieagentur Weststeiermark, Stainz

  48. Arbeitsgemeinschaft QM Holzheizwerke (2004) Q-Leitfaden. Checkliste Standard, MS1 – MS5. 2004. Centrales Agrar-Rohstoff Marketing- und Energie-Netzwerk e.V. (C.A.R.M.E.N. e.V.), Straubing

  49. IZES gGmbH [Institut für ZukunftsEnergie und Stoffstromsysteme], ifeu [Institut für Energie- und Umweltforschung GmbH] (2010) Nachhaltige Biogasnutzung im Raum Trier. Abschlussbericht für die Stadtwerke Trier. IZES gGmbH, ifeu, Trier

  50. Daniel J, Scholwin F, Vogt R (2008) Materialband D. Biogasnutzung, im Rahmen des BMU Forschungsvorhabens „Optimierung für einen nachhaltigen Ausbau der Biogaserzeugung und –nutzung in Deutschland, FKZ: 0327544. Heidelberg

  51. Scheftelowitz M, Rensberg N, Denysenko V, Daniel-Gromke J, Stinner W, Hillebrand K, Naumann K, Peetz D, Henning C, Thrän D, Beil M, Kasten J, Vogel L (2015) Stromerzeugung aus Biomasse (Vorhaben IIa Biomasse). Zwischenbericht Mai 2015. DBFZ [Deutsches Biomasseforschungszentrum gemeinnützige GmbH], Fraunhofer Institut für Windenergie und Energiesystemtechnik (IWES), INL GmbH, Helmholtz-Zentrum für Umweltforschung - UFZ, Leipzig

  52. C.A.R.M.E.N. e.V. [Centrales Agrar-Rohstoff Marketing- und Energie-Netzwerk e.V.] Nahwärmenetze und Bioenergieanlagen. Ein Beitrag zur effizienten Wärmenutzung und zum Klimaschutz. C.A.R.M.E.N. e.V., Straubing

  53. KfW [Kreditanstalt für Wiederaufbau] (2018) Merkblatt: Erneuerbare Energien Premium. Finanzierung von Maßnahmen zur Nutzung Erneuerbarer Energien im Wärmemarkt. KfW, Frankfurt am Main

  54. Ökoprojekte Gronbach GmbH (2019) Nahwärme: Wärmenetz Uttenhofen.

  55. Gemeinde Buttenwiesen (2019) Das Gute liegt so nah: Nahwärme Buttenwiesen. Gemeinde Buttenwiesen.

  56. Fachzentrum für Energie & Landtechnik (2016) Nahwärmeversorgung Kastl. Gemeinde Kastl.

  57. Fachverband Biogas e.V. (2019) Bioenergiedorf Reiffenhausen – Gemeinsinn als Schlüssel zum Erfolg.

  58. Energievision Frankenwald e.V. (2019) Technische Daten zum Nahwärmenetz in Gössersdorf. Energievision Frankenwald e.V.

  59. EnergieAgentur.NRW GmbH (2012) Biogas von Stadtteil zu Stadtteil – neues Nahwärmenetz in Schmallenberg. EnergieAgentur.NRW GmbH. Accessed 13.11.2019

  60. DVL [Deutscher Verband für Landschaftspflege e.V.] (2019) Biogasanlage Ribbeck, Brandenburg. DVL [Deutscher Verband für Landschaftspflege e.V.]. Accessed 13 Nov 2019

  61. Dorfgemeinschaft Ebbinghof (2019) Nahwärmenetz. Dorfgemeinschaft Ebbinghof.

  62. Bioenergiedorf Oberrosphe eG (2019) Bioenergiedorf Oberrosphe eG. Bioenergiedorf Oberrosphe eG. Accessed 13.11.2019

  63. Berghamer O (2013) Nahwärmeversorgung in Mardorf.

  64. Krah A (2017) Bebauungsplan mit integriertem Grünordnungsplan: SO Biogasanlage Oberindling. Pocking

    Google Scholar 

  65. Kompetenzzentrum HessenRohstoffe (HeRo) e.V. (2009) Bioenergiedörfer in Hessen – mit Erneuerbaren Energien in die Zukunft. Kompetenzzentrum HessenRohstoffe (HeRo) e.V.,, Witzenhausen

  66. Fuchs G (2013) Praxisbeispiel - Satelitten BHKW: Abwärme für Wohnhäuser undfür einen Gasthof. Lengdorf

    Google Scholar 

  67. Agentur für Erneuerbare Energien e.V. (2018) Energie vom Land hält Warm. Agentur für Erneuerbare Energien e.V, Berlin

    Google Scholar 

  68. Ahlers E, Seidel M, Burg C, Huber R, Hügging T, Kettrup H, Makrlik N, Renken A, Sachse U, Straub F, Sztochaj T, Geiger B, Kleeberger H, Pätzold K (2007) Abwicklung von Standardlastprofilen. BGW [Bundesverband der deutschen Gas- und Wasserwirtschaft], Bonn

  69. Aktualisierte und erweitere Testreferenzjahre (TRY) von Deutschland für mittlere und extreme Witterungsverhältnisse (2019) BBSR [Bundesinstitut für Bau,- Stadt,- und Raumforschung im Bundesamt für Bauwesen und Raumordnung]. Accessed 13 Sep 2019

  70. Fallahnejad M, Hartner M, Kranzl L (2018) Fritz S (2018) Impact of distribution and transmission investment costs of district heating systems on district heating potential. Energy Procedia 149:141–150.

    Article  Google Scholar 

  71. Dotzauer M, Lauer M, Schneider J Servant of two masters? Trade-offs for agricultural biogas plants between flexible power generation and heat cogeneration. In: 24th European Biomass Conference and Exhibition, Amsterdam, 2016. doi:10.5071/24thEUBCE2016-5AO.9.1

  72. AGEB [AG Energiebilanzen e.V.] (2013) Anwendungsbilanzen für die Endenergiesektoren in Deutschland in den Jahren 2010 und 2011, Studie beauftragt vom Bundesministerium für Wirtschaft und Technologie Projektnummer: 23/11. AGEB, Berlin

  73. Wolff D, Jagnow K (2011) Untersuchung von Nah- und Fernwärmenetzen. Überlegungen zu Einsatzgrenzen und zur Gestaltung einer zukünftigen. Fern- und Nahwärmeversorgung, Wolfenbüttel, Braunschweig

    Google Scholar 

  74. Zahlen und Fakten. Energiedaten (2019) BMWi [Federal Ministry for Economic Affairs and Energy].

  75. Herbes C, Halbherr V, Braun L (2018) Preise für die Abgabe von Wärme von Biogasanlagen an Dritte. Agrarbetrieb 2018(1):12–15

    Google Scholar 

  76. The German Federal Government (2019) Klimaschutzprogramm 2030 der Bundesregierung zur Umsetzung des Klimaschutzplans 2050. The German Federal Government,

    Google Scholar 

  77. The German Federal Government (2019) CO2-Ausstoß verbindlich senken. The German Federal Government.

Download references


Joint research on this paper has been undertaken by the institutions UFZ, DBFZ, and IZES in the context of the project “Bioenergy - potentials, long-term perspectives and strategies for power generation plants after 2020 (BE20plus)”.


This research received funding from the Agency for Renewable Resources e.V. (FNR) and the Federal Ministry for Food and Agriculture (BMEL). The funding bodies had no influence over the design of the study, the interpretation of the data, or the writing of the manuscript.

Author information

Authors and Affiliations



Several authors drafted the different parts of the initial version of this manuscript. All authors contributed to the revising and finalizing of the article. All authors read and approved the final manuscript.

Corresponding author

Correspondence to M. Steubing.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.



Table 11 Total network length, distance of BEP to first potential customer and resulting network length factor of 20 existing DHNs fueled by heat generated from biomass in Germany

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Steubing, M., Dotzauer, M., Zakaluk, T. et al. Bioenergy plants’ potential for contributing to heat generation in Germany. Energ Sustain Soc 10, 14 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: