Investigation on Relative Heat Losses and Gains of Heating and Cooling Networks

Environmental and Climate Technologies 2021, vol. 25, no. 1, pp. 479–490 https://doi.org/10.2478/rtuect-2021-0035 https://content.sciendo.com 479 ©2021 Violeta Madan, Ingo Weidlich. This is an open access article licensed under the Creative Commons Attribution License (http://creativecommons.org/ licenses/by/4.0). Investigation on Relative Heat Losses and Gains of Heating and Cooling Networks Violeta MADAN1*, Ingo WEIDLICH2 1, 2 HafenCity University, Henning-Voscherau-Platz 1, Hamburg, 20457, Germany Abstract – The integration of district heating (DH) and cooling (DC) in the sustainable energy system of the future requires a significant reduction in operating temperatures. Supply temperatures below 70 °C are required for new 4th Generation DH. Main benefits are the use of low exergy heat sources and the reduction of heat losses. The reduction of heat losses is achieved by reducing the driving temperature difference between the medium pipe and the ground. The decrease of the return temperature level is limited by the consumer behaviour and the ground temperature level. As a consequence, the reduction of the supply temperature is accompanied by a reduction of the maximum transmittable heat flow. For energy efficiency and economic reasons, the relative heat losses are therefore an important design value for DH networks. The study proposes an approach to estimate the relative heat losses by using steadystate heat loss models and analyses the values for different DH generations. In particular, due to the rising of the near-surface soil temperature, the relative cold losses are also studied. Keywords – Cold losses; district cooling; district heating; heat losses Nomenclature A Cross-sectional area m² cp Heat capacity J/(kgK) D Half of the distance between the pipe axes m di, do Inner and outer insulation diameter m h Heat loss factor – T Temperature K Q  Heat flow rate W q Length related heat flow W/m ν Flow velocity m/s λg, λi Heat conductivity of the ground and insulation W/(mK) H Burial depth m R Thermal Resistance mK/W 1. INTRODUCTION District heating (DH) and district cooling (DC) can mainly contribute to a sustainable energy system based on renewable energy production. The temporal and local uncoupling of *Corresponding author. E-mail address: [email protected] Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 480 the demand and supply are advantageous for the integration of renewable energy sources and the sector coupling approach [1]. Other benefits are the significant improvement of the energy efficiency for high population density areas and more flexible exchange of the energy source. The development of DH extends over four generations and includes the optimization of materials, heat sources, pipe design and operating conditions. Decreasing of the operating temperature is one of the main objectives for newer DH generations, as a lower temperature level favours the integration of low temperature heat sources and the reduction of heat losses. Lund et al. define a supply temperature level below 100 °C for the 3rd and 30–70 °C for the 4th Generation DH [1]. The heat losses and gains are the main influencing parameters for design and efficiency of a DH or DC system. It is a major interest to minimize the heat losses (and gains in case of DC) for energy savings and economic reasons. Additional approaches to decrease heat losses are the improvement of the pipe insulation and the application of twin or triple pipes [2]. However, lower operation temperatures promise the higher impact compared to material changes in the pipe systems. For the determination of steady–state heat losses, analytic solutions [3] and explicit solutions for specific geometries are available [4], [5]. The explicit solution of Wallentén is implemented in DIN EN 13941–1, which represents the state of the art for the heat loss calculation in district heating [6]. Wallentén proposes an explicit solution for one or two insulated pipes and twin pipes in the ground by applying the multipole method [5]. Kvisgaard and Hadvig introduce a corrected depth for the pipe axis to consider the heat transfer coefficient at the ground surface [4]. This assumption is also adapted in DIN EN 13941–1 [6]. The model of Kvisgaard and Hadvig has been evaluated by Oppelt et al. for DH and DC alication. In [7] a maximum deviation of 2.0–2.7 % for DH and 3.8–6.2 % for DC was reported for the Kvisgaard/Hadvig model compared to simulations with the finite volume method in the case of two insulated single pipes. Bøhm reviewed different methods for the calculation of steady–state heat losses and observed the undisturbed temperature of the ground as a critical value for the existing steady– state models. In his study, he did experiments and FEM simulations to find the location of an undisturbed ground temperature that can be applied in the steady-state theory to achieve acceptable accuracy. He found the location of this temperature approximately at the top of the casing for preinsulated pipes [8]. However, other researchers showed a measurable increase of the temperature near to DH pipes in the field [9], [10]. Another factor is the change in the near–surface temperature field due to climate change, urbanization and thermal exposure of the soil. Henning and Limberg found a ground temperature rise of more than 4 K in urban regions compared to rural regions in Berlin [11]. Regarding low temperature DH, two effects have to be considered: The decrease of the supply (and return) temperature on the one hand and the decrease of temperature gradient between the return and supply pipe on the other hand. The reduction of the supply temperature level goes along with a significant reduction of the absolute heat losses of the system. The decrease of temperature gradient between the supply and return pipe reduces the transferred heat energy to the consumer. This loss of performance could be compensated by the trend of low energy buildings [1]. Strategies to supply existing buildings with temperatures below 70 °C are described in [12]. In this context the relative heat losses may be the main design value for new DH generations. The relative heat losses are usually calculated depending on the annual heat consumption [13] or by a detailed modelling of the specific network [14], [15]. If no measurements of the consumption data are available, the estimation of the annual demand for space heating and domestic hot water is prone to errors [13], [14]. The study proposes an approach to estimate Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 481 relative heat losses and gains by applying steady–state heat loss models. The calculation method can easily be adapted by the network operator to get a first impression of the network efficiency and capacity required for the heat production. In the following study, the relative heat losses are calculated with the proposed method and the results are compared for different generations DH networks and DC. Heating and cooling networks in Hamburg and surroundings have been selected to create realistic boundary conditions for the calculations. 2. METHODOLOGY In the following, the database, the heat loss models and the model parameters applied for the calculations will be described. 2.1. Database Hamburg’s DH network and several other local heating networks are covering about a quarter of the heating demand of the city. The thermal register of Hamburg provides an overview of the heating networks and operators [16]. For the investigation of the heat losses, different heat networks are selected that meet the temperature requirements for the second, third and fourth DH generations [1]. Hamburg’s DH network provides a supply temperature of 133 °C for the standard outer temperature of −12 °C for Hamburg according to the national annex of DIN EN 12831 [17]. Due to this high temperature level, Hamburg’s DH network is classified as 2nd Generation DH. To represent the 4th Generation DH (TS < 70 °C), the supply and return temperatures of a new residential area ‘Waisengärten’ in Schwerin, a city in northern Germany, have been selected for the study. The new residential buildings in the area of ‘Waisengärten’ is operated by Stadtwerke Schwerin GmbH, the public utility provider of Schwerin [18]. The energy concept for the new buildings combines geothermal heat supply with home stations to regulate the domestic hot water temperature [19]. The temperatures of the district and local heating networks selected for this study are taken from the official technical connection conditions of the network operators [18], [20]–[22] and are summarized in Table 1. The return temperatures in the official technical connection conditions are maximum return temperatures allowed by the operators. For the following investigation, a winter and summer case are considered to investigate the relative heat losses and gains of the systems. TABLE 1. SELECTED HEATING AND COOLING NETWORKS System operator Supply/return temperature winter case (−12 °C) Supply and return temperature summer case (20 °C) District Heating Generation according to Lund et al. [1] Wärme Hamburg GmbH 133/60 °C 90/20 °C 2nd Innogy SE 100/55 °C 69/20 °C 3rd HanseWerk Natur GmbH 80/45 °C 70/45 °C 3rd Stadtwerke Schwerin GmbH 56/35 °C 56/35 °C 4th Wärme Hamburg GmbH (Cooling network) 8/15 °C 6/12 °C – A higher accuracy for determination of heat losses and gains can be obtained by knowing the ground temperature profile [23]. Therefore, historical weather data of Hamburg listed in DIN 4710:2003–01 [24] has been applied in the steady state heat loss models presented in section 2.2. The ground temperatures in January and July are representative for the winter Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 482 and summer case. For the sake of a classification the resulting temperature profiles for Hamburg were drawn together with generalized curves according to Dahlem in Fig. 1. It is important to point out that the weather data summarized in DIN 4710:2003–01 [24] is the average temperature for the years 1961–1990. Current weather data from the German Weather Service (Deutscher Wetterdienst) [25] suggests that the ground temperatures slightly increased in the last decades and the ambient temperature can rise up to 35 °C during the summer season. Fig. 1. Ground temperature profiles of the city of Hamburg compared with curves according to Dahlem without ground water [23], [24]. 2.2. Model for Steady-State Heat Loss District heating networks usually consist of two insulated symmetrical buried pipes with a circular cross–section as shown in Fig. 2. For the calculation of the steady–state heat losses, pure heat conduction is usually assumed, while convection, moisture transport and phase changes are neglected. Other assumptions are homogeneous material properties and no thermal resistances of the medium pipe and casing [10]. Radial heat conduction through the insulation to the ground and between the pipes is the process described by the most steady– state heat loss models for DH application. In the following, the term ‘heat loss’ refers to the length-related heat loss in W/m if not otherwise specified in the text. Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 483 Fig. 2. Geometry of two insulated pipes in the ground. Wallentén describes the heat loss of each DH pipe as a superposition of symmetrical and anti–symmetrical heat losses referring to the multipole method [3]. Thus, the heat losses of the supply and return pipe can be calculated according to Eq. (1) and Eq. (2). ,, s sym s a s qq q= + (1) ,, r sym r a r qq q= − (2) The symmetrical and anti–symmetrical heat losses (3) and (4) can be calculated by applying the temperatures Tsym and Ta defined in Eq. (5) and Eq. (6). ( )2 sym sym g g sym q TT h= − ⋅ πλ ⋅ (3) 2 a a ga qT h= ⋅ πλ ⋅ (4) 2 sr sym TT T+ = (5) 2 sr a TT T− = (6) The first–order multipole formula to calculate the reciprocal of the heat loss factor for the symmetrical problem is defined in Eq. (7). The formula for the anti–symmetrical heat loss factor is given in [5]. The dimensionless thermal resistance parameter β is defined in Eq. (8). 22 2 00 0 222 1 2 0 4416( ) 4 ln ln 1 1 14 sym o dd d D H DH HH hdD d D −     ++       ⋅+       = +β+ + −     +β    + −β  (7) ln go ii d d λ β=  λ (8) Wallentén quantifies the relative errors of the calculated heat losses as less than 0.5 %, when the first order formulae are used [5]. For the following analysis, the Wallentén and Tg, λg Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 484 Kvisgaard/Hadvig model are applied to calculate the system heat losses. The two models differ in the calculation of the thermal soil resistance and in the value for the ground temperature. The corrected buried depth implemented by Kvisgaard and Hadvig to consider the heat transfer coefficient at the soil surface is defined in Eq. (9). For a detailed description of the Kvisgaard/Hadvig model refer to [4]. * 0.0685 g HH= + ⋅λ (9) The ground temperature applied for the heat loss calculations is a key difference between the Wallentén and Kvisgaard/Hadvig model. Wallentén proposes the temperature on the ground surface, while Kvisgaard and Hadvig apply the temperature at the ground level of the pipe axis at installation depth. To investigate the influence of the ground temperature on calculated heat losses, the results obtained by using the original models are presented among the results in section 3. Bøhm investigated the location of an undisturbed ground temperature for steady–state heat loss calculation to achieve accurate results. The undisturbed ground temperature position for preinsulated pipes found by Bøhm corresponds to the ground temperature at the top of the casing and is applied in this study. For this purpose, the function equation for the city of Hamburg derived using the weather data in DIN 4710:2003–01 is applied [24]. To calculate the relative heat losses, the sum of the length related heat losses (1) and (2) is divided by the maximum amount of heat transported in the network according to the basic equation of thermodynamics (10). The relevant cross–sectional area A is the cross–sectional area of the water pipe. Due to the stationary condition, only the full load case is considered. According to the AGFW (German District Heating Association), the steady–state condition is only applicable for network lengths below 10 km [26]. () p sr Q Avc T T=ρ⋅ ⋅ ⋅ ⋅ −  (10) The resulting parameter proposed in this study to estimate relative heat losses is a percentage heat loss per meter of pipe length according to Eq. (11). The corresponding unit is %/m.  100% sr qq q Q + = ⋅  (11) 2.3. Model Parameters The calculation parameters are selected according to typical operating conditions of DH networks. Key parameters are summarized in Table 2. The steel pipe dimension of DN 200 can be applied for district heating and cooling; the dimensions of the pipe and insulation comply with the standards DIN EN 253 and DIN EN 17415–1 [27], [28]. Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 485 TABLE 2. CALCULATION PARAMETERS Parameter Symbol Value Heat conductivity of the ground λg 1.2 W/(mK) Heat conductivity of the insulation λi 0.027 W/(mK) Cross-sectional area A 0.0346 m² Heat capacity of water cp 4190 J/(kgK) Inner insulation diameter di 0.219 m Outer insulation diameter do 0.315 m Distance between the pipe axes 2D 0.565 m Buried depth H 0.9575 m Flow velocity ν 1 m/s 3. RESULTS For the networks listed in Table 1 the absolute heat losses and gains calculated with the Wallentén and Kvisgaard/Hadvig model are visualized in Fig. 3 and Fig. 4, each for winter and summer case. The results presented in Fig. 2 and Fig. 3 are obtained by applying the ground temperatures proposed by the original authors to visualize the influence of the selected ground temperature on the heat losses. For the Wallentén model, a ground surface temperature of 0.8 °C (winter) and 18.5 °C (summer) and for the Kvisgaard/Hadvig model 4.6 °C (winter) and 14.9 °C (summer) are applied [24]. Fig. 3. Calculated heat losses – winter case. 0 10 20 30 40 50 60 70 80 0 10 20 30 40 50 60 70 80 133/60 100/55 80/45 56/35 8/15 Deltat |supply-return|, K Heat losses and cold gains, W/m Network supply/return temperature, °C Winter Wallenten model Kvisgaard/Hadvig model deltaT supply-return Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 486 Fig. 4. Calculated heat and cold losses – summer case. Fig. 5. Calculated relative heat losses – winter case. The resulting heat losses obtained with the Wallentén and Kvisgaard/Hadvig model show a constant deviation of about 3 W/m. The relative deviation between the calculated heat losses increases with decreasing temperature gradient between the water pipe and ground temperature. The relative deviation has a range of 4–12 % for the DH and reaches the maximum value of 38 % for DC. The relative deviation between the two models is less than 0.5 % if the same value is applied for the ground temperature. -20 0 20 40 60 80 -10 -5 0 5 10 15 20 25 30 35 40 90/20 69/20 70/45 56/35 6/12 Deltat |supply-return|, K Heat-and cold losses, W/m Network supply/return temperature, °C Summer Wallenten model Kvisgaard/Hadvig model deltaT supply-return 0 10 20 30 40 50 60 70 80 0 0.2 0.4 0.6 0.8 1 1.2 133/60 100/55 80/45 56/35 8/15 Deltat |supply-return|, K Relative heat losses and cold gains, %/km Network supply/return temperature, °C Winter rel. heat losses rel. cold gains delta T supply-return a Environmental and Climate Technologies ____________________________________________________________________________ 2021 / 25 487 The Wallentén model with the undisturbed ground temperature proposed by Bøhm is applied to determine the relative heat losses. The ground level of the top of the pipe casing corresponds to a depth of 0.8 m. The undisturbed ground temperature results in 4 °C for the winter case and 15.5 °C for the summer case. The resulting relative heat losses per kilometre pipe length are presented in Fig. 5 and Fig. 6. Fig. 6. Calculated relative heatand cold losses – summer case. 4. DISCUSSION In the previous section the results obtained by the Wallentén and Kvisgaard/Hadvig model for steady–state heat losses of different heating and one cooling network have been presented in Fig. 3 and Fig. 4. A constant deviation of about 3 W/m occurs between the results due to different ground temperatures implemented in the original models. As the heat losses are constantly decreasing with decreasing temperature gradient between the medium pipe and ground temperature, the relative deviation of the heat losses rises up to 38 % for the DC case. These results of the different steady–state models for low temperature DH and DC show the significance of an accurate selection of the undisturbed ground temperature for these applications. In this context, an optimized ground temperature position for the 4th Generation DH and for DC should be further investigated. The overall trend of decreasing heat losses for newer DH generations due to the decrease of the driving temperature gradient between the pipe and the ground is confirmed in the results for the winter and summer case. The outliners for the 90/20 °C and 69/20 °C networks in the summer can be explained by the lower return temperature of 20 °C compared to the other investigated heating networks. In this work, the percentage heat loss per meter of pipe length is proposed as a parameter to estimate the relative heat losses and gains of a system. In Fig. 5 and Fig. 6 the relative heat losses are presented for a winter and summer case. In the winter case, the relative heat losses of the investigated heating networks increase, while the temperature level is decreasing. As the absolute heat losses decrease as shown in Fig. 3, the influence of the slope of the transported heat flux 𝑄𝑄󰇗 causes the rising trend. For 0 20 40 60 80 100 0 0.2 0.4 0.6 0.8 1 90/20 69/20 70/45 56/35 6/12 DeltaT |supply-return|, K Relative heat-and cold losses, %/km Network supply/return temperature, °C Summer rel. heat losses rel. cold losses deltaT supply-return