# Study on Engineering Design and Cost Performance of Borehole Heat Exchangers

Authors: Lei Yu; Jiemin Cheng
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The ground source heat pump (GSHP) is an environmentally friendly and energy-efficient technology, and the investigations on the simulation models of borehole heat exchanger (BHE) have made progress. However, the researches on engineering design and economic performance based on the simulation models of BHE are a few. This paper firstly describes the interior and exterior heat transfer according to the existing achievements, and the corresponding characteristics of models are explained. Afterwards, the detailed design procedures of BHE size are provided and the cost performance is highlighted. Accordingly, an actual engineering project which employs GSHP system as the air conditioning system is researched; the necessary expenses on BHE, heat pump units, and indoor end equipment are all taken into account. The initial investment, operating cost, and other factors that exert influences on the cost performance are demonstrated. In addition, comparisons between different systems are conducted to show the advantages of GSHP. By means of thermal transfer analysis and the design of engineering project, the practicability, energy saving, and environmental protection of GSHP are brought into play, which means the promotion of this technology is the general trend of development of air condition system.


---

## Body

## 1. Introduction

The ground source heat pump (GSHP) system is the most advanced green air-conditioning technology in the world, because GSHP is favorable to achieve environmental protection and energy-saving. The underground low-grade energy is utilized with the help of heat pump unit and therefore the heating, cooling, and domestic hot water of buildings are all achieved [1]. More specifically, the thermal transfer is conducted from the circulating liquid of the heat exchange tube to the underground medium, which means the renewable energy existing in underground space can be made use of. Underground mediums locating at the depth of more than ten meters keep nearly constant temperature all the year round because they are free of the thermal disturbance that comes from the outside environment. The depth of underground mediums can achieve one hundred meters while GSHP is employed. Underground temperatures are lower and higher than outdoor air’s temperatures in summer and in winter, respectively [2, 3]. Therefore, the technical barriers such as frosting that appears at air source heat pump can be overcome and then the system efficiency is improved significantly. For example, the 4 kW heat quantity is acquired if only 1 kW electric energy is consumed; accordingly, the energy efficiency ratio is larger than 4. In addition, GSHP that employs borehole heat exchanger (BHE) as the underground heat transfer elements does not extract groundwater and thus the system is a new energy-saving technology keeping sustainable development. The whole system usually consists of several parts including BHEs, the heat pump unit, and end equipment, and sometimes domestic hot water is contained either.Also, the advantages of GSHP should be followed with interest and they are summarized as several aspects. Firstly, high-efficiency and energy-saving are undeniable because running efficiency of GSHP is 30~50% larger than that of traditional air-conditioning system. Secondly, there is no boiler and the corresponding room and therefore the problems of both air pollution and greenhouse gases emission are prevented. Thirdly, an equipment can serve several purposes. Next, the objectives of beautifying building and household metering are both realized when the GSHP is adopted in buildings. Last but not least, the characteristics of the safety and reliability should be emphasized as the life spans of BHE and heat pump unit, respectively, attain 50 and 20 years. Figure1 shows the schematic diagram of BHEs that are responsible for the heat or cooling load of the building [4].Figure 1
The schematic diagram of underground heat exchanger employed for cooling or heating.It can be observed that the underground heat transfer components of one engineering project are composed of a number of BHEs. No matter whether to study the theoretical models or to analyze the reasonable design of GSHP system, it is related to the BHEs.The simulation models describing the heat transfer inside and outside borehole have gained more and more attention, and they have been put to use in theoretical calculation or engineering design. Apart from that, the design procedures of obtaining the BHE’s reasonable size are becoming increasingly applicable. It is a pity that the analyses on engineering projects based on both thermal transmission models and design procedures are a few, and this is a noticeable research gap. The novelty and originality of the paper are to investigate design procedure and cost performance of engineering projects according to the heat transfer simulation models. Afterwards, the detailed design plane for an actual project is conducted to certify the necessity of applying GSHP system; the comparison between GSHP and other systems is made to show the significant economic performance of GSHP system.


## 2. Simulation Models of Heat Transfer

### 2.1. Heat Transfer Inside Borehole GHE

Borehole is drilled and then the U-tubes are installed into it [5]. The vertical BHE is produced if the backfilling materials are poured, and these materials have the effects of sealing and improving heat transfer performance. The heat transfer domain from U-tube to underground medium is divided into two parts, that is, the areas inside and outside borehole [6]. The steady-state mode can be put into use to deal with the heat exchange inside GHE because the geometric size and thermal capacity are relatively small compared with those outside borehole. The equivalent diameter is employed to simplify the model and it is usually applied to explore the heat transfer between U-tube and the borehole wall [7]; the principle is that all branch pipes are simplified as a pipe with larger diameter; therefore, the two-dimensional conduction can be regarded as a one-dimension radial heat exchange.The heat conductivities of backfilling material, pipe material, and underground medium are, respectively,kb,kp, and ks. The radius of borehole, external radius of pipe, and inner radius of pipe are, respectively, rb,r0, and ri. The layout drawing of U-tube inside borehole and the corresponding equivalent diameter are both shown in Figure 2.Figure 2
The diameter and the corresponding equivalent diameter of U-tube.The equivalent diameter of U-tubes is approximatelynd0 in case the number of branch tubes is n; thus, the thermal-convection resistance from the circulating fluid to the inner wall of tube is expressed in (1) while the convection coefficient is h:(1)Rf=12πrih.Then, the thermal resistance between inner wall and external wall of tube is concluded in(2)Rt=12πkpln⁡nr0nr0-r0-ri.Afterwards, the heat is further conveyed to the borehole wall when the backfilling material acts as the thermal conductor, and the corresponding resistance is shown in(3)Rb=12πkbln⁡rbnr0.The total resistance inside borehole should be summarized according to the analyses above, and the expression of energy equation is given in (4) while the temperatures of circulating liquid and the borehole wall are, respectively, tf and tb. There is no doubt that the parameters such as driving force of heat transfer, thermal resistance and heat flux are shown clearly. Three components constitute the total resistance inside BHE, and the original source of heat energy is the circulating liquid. ql is the heat transfer rate of BHE, that is, the heat transfer quantity per meter BHE:(4)tf-tb=Rf+Rt+Rbql.


### 2.2. Heat Transfer Outside Borehole GHE

#### 2.2.1. Line Heat Source Model

Obviously, the heat continues to be transmitted to far locations of underground medium. The heat transfer process outside BHE plays a significant role in carrying energy, because the borehole’s diameter is so small that the heat capacity inside BHE is little and the thermal resistance outside borehole is much larger than that inside borehole. The investigation on heat transfer outside borehole is becoming increasingly important and therefore the reasonable mathematical simulation models should be taken into account. The borehole’s depth is far larger than its diameter and thus BHE can be regarded as a line heat source [8], which means the radius of BHE is ignored when the heat transfer outside borehole is studied. The diagram of line heat source is described in Figure 3.Figure 3
The line heat source adopted for BHE.Accordingly, a most widely used one-dimensional model for this purpose is Kelvin’s line source model, in which the borehole is replaced by a line heat source with its radial dimension neglected. However, the actual depth of any BHE is finite and then the finite line heat source is proposed; the Kelvin’s model lays a firm foundation for the research of the finite line heat source [9]. In order to investigate the conduction around the borehole, Green’s function theory is employed and the corresponding expression is given in (5) showing the temperature response of any point (x,y,z) except heat source in the underground medium at time τ, and the point heat source locating at (x′,y′,z′) emits heat at time τ′ [10]:(5)Gx,y,z,τ;x′,y′,z′,τ′=18πaτ-τ′3/2exp⁡-x-x′2+y-y′2+z-z′24aτ-τ′.After that, the temperature response induced by the finite line heat source is obtained in (6), where erfc is the complementary error function and t0 is the initial temperature of the underground medium [11]:(6)t-t0=ql4πk∫0Herfcx-x′2+y-y′2+z-z′2/2aτx-x′2+y-y′2+z-z′2-erfcx-x′2+y-y′2+z-z′2/2aτx-x′2+y-y′2+z+z′2dz′.To facilitate the expression, the nondimensional parameters are employed and the corresponding forms areΘ=k(t-t0)/ql, X=x/rb, Y=y/rb, Z=z/rb, X′=x′/rb, Y′=y′/rb, Z′=z′/rb, H=h/rb, Fo=aτ/rb2. Equation (6) is translated into(7)Θ=14π∫0HerfcX-X′2+Y-Y′2+Z-Z′2/2FoX-X′2+Y-Y′2+Z-Z′2-erfcX-X′2+Y-Y′2+Z+Z′2/2FoX-X′2+Y-Y′2+Z+Z′2dZ′.The temperature responses with the time are shown in Figure4 according to the calculation results of (7).Figure 4
The temperature responses of finite line heat source with the time.It is found that all the temperature responses will be stable at last. Furthermore, the larger the length, the fiercer the response.


#### 2.2.2. Cylindrical Heat Source Model

The heat capacity inside borehole can be ignored as the diameter is extremely small compared to its depth, but the existence of inner heat transfer is undisputed. Accordingly, the line heat source model is not an ideal choice for simulating the heat transfer of BHE. The solid cylindrical heat sources have been proposed in recent years; in this model, it is supposed that the cylinder is no longer a cavity but is filled with the medium identical to that out of the cylinder [12–14]. The whole infinite domain is composed of a homogeneous medium and therefore the heat capacity of borehole is taken into consideration. The diagram of cylindrical heat source is shown in Figure 5.Figure 5
The cylindrical heat source adopted for BHE.The Green function is still the basic method to obtain the temperature response caused by the heat source, and the temperature response of any point in the underground medium except heat source is shown in(8)t-t0=ql8πk∫0τ1τ-τ′I0x2+y2rb2aτ-τ′exp⁡-x2+y2+rb24aτ-τ′·erfcz-H2aτ-τ′+erfcz+H2aτ-τ′-2erfcz2aτ-τ′dτ′,where I0 is the zero-order modified Bessel function and the borehole is supposed to be arranged at z-axis. Again, dimensionless method is used to simplify the expression and the corresponding formula is given in(9)Θ=18π∫0FodFo′Fo-Fo′I0X2+Y22Fo-Fo′exp⁡-X2+Y2+14Fo-Fo′·erfcZ-H2Fo-Fo′+erfcZ+H2Fo-Fo′-2erfcZ2Fo-Fo′.The temperature responses with the time are shown in Figure6, and the regular patterns and trends are similar with those of line heat source; that is, all responses must attain the steady states and the response value increase with the borehole length if the time is constant.Figure 6
The temperature responses of cylindrical heat source with the time.As far as the differences between line and cylindrical sources are concerned, the temperature responses changing with time are described in Figure7 when the radius and depth of borehole are constant.Figure 7
The comparison of temperature responses between line and cylindrical heat sources.It is noticed that the temperature response of cylindrical heat source is always larger than that of line heat source. The dimensionless temperature response isΘ=k(t-t0)/ql, and the heat transfer rate of line heat source becomes larger while the actual temperature rise is stable. Therefore, the calculational length of BHEs by means of line heat source is shorter than that depending on cylindrical heat source when the design of BHEs is conducted.


## 2.1. Heat Transfer Inside Borehole GHE

Borehole is drilled and then the U-tubes are installed into it [5]. The vertical BHE is produced if the backfilling materials are poured, and these materials have the effects of sealing and improving heat transfer performance. The heat transfer domain from U-tube to underground medium is divided into two parts, that is, the areas inside and outside borehole [6]. The steady-state mode can be put into use to deal with the heat exchange inside GHE because the geometric size and thermal capacity are relatively small compared with those outside borehole. The equivalent diameter is employed to simplify the model and it is usually applied to explore the heat transfer between U-tube and the borehole wall [7]; the principle is that all branch pipes are simplified as a pipe with larger diameter; therefore, the two-dimensional conduction can be regarded as a one-dimension radial heat exchange.The heat conductivities of backfilling material, pipe material, and underground medium are, respectively,kb,kp, and ks. The radius of borehole, external radius of pipe, and inner radius of pipe are, respectively, rb,r0, and ri. The layout drawing of U-tube inside borehole and the corresponding equivalent diameter are both shown in Figure 2.Figure 2
The diameter and the corresponding equivalent diameter of U-tube.The equivalent diameter of U-tubes is approximatelynd0 in case the number of branch tubes is n; thus, the thermal-convection resistance from the circulating fluid to the inner wall of tube is expressed in (1) while the convection coefficient is h:(1)Rf=12πrih.Then, the thermal resistance between inner wall and external wall of tube is concluded in(2)Rt=12πkpln⁡nr0nr0-r0-ri.Afterwards, the heat is further conveyed to the borehole wall when the backfilling material acts as the thermal conductor, and the corresponding resistance is shown in(3)Rb=12πkbln⁡rbnr0.The total resistance inside borehole should be summarized according to the analyses above, and the expression of energy equation is given in (4) while the temperatures of circulating liquid and the borehole wall are, respectively, tf and tb. There is no doubt that the parameters such as driving force of heat transfer, thermal resistance and heat flux are shown clearly. Three components constitute the total resistance inside BHE, and the original source of heat energy is the circulating liquid. ql is the heat transfer rate of BHE, that is, the heat transfer quantity per meter BHE:(4)tf-tb=Rf+Rt+Rbql.


## 2.2. Heat Transfer Outside Borehole GHE

### 2.2.1. Line Heat Source Model

Obviously, the heat continues to be transmitted to far locations of underground medium. The heat transfer process outside BHE plays a significant role in carrying energy, because the borehole’s diameter is so small that the heat capacity inside BHE is little and the thermal resistance outside borehole is much larger than that inside borehole. The investigation on heat transfer outside borehole is becoming increasingly important and therefore the reasonable mathematical simulation models should be taken into account. The borehole’s depth is far larger than its diameter and thus BHE can be regarded as a line heat source [8], which means the radius of BHE is ignored when the heat transfer outside borehole is studied. The diagram of line heat source is described in Figure 3.Figure 3
The line heat source adopted for BHE.Accordingly, a most widely used one-dimensional model for this purpose is Kelvin’s line source model, in which the borehole is replaced by a line heat source with its radial dimension neglected. However, the actual depth of any BHE is finite and then the finite line heat source is proposed; the Kelvin’s model lays a firm foundation for the research of the finite line heat source [9]. In order to investigate the conduction around the borehole, Green’s function theory is employed and the corresponding expression is given in (5) showing the temperature response of any point (x,y,z) except heat source in the underground medium at time τ, and the point heat source locating at (x′,y′,z′) emits heat at time τ′ [10]:(5)Gx,y,z,τ;x′,y′,z′,τ′=18πaτ-τ′3/2exp⁡-x-x′2+y-y′2+z-z′24aτ-τ′.After that, the temperature response induced by the finite line heat source is obtained in (6), where erfc is the complementary error function and t0 is the initial temperature of the underground medium [11]:(6)t-t0=ql4πk∫0Herfcx-x′2+y-y′2+z-z′2/2aτx-x′2+y-y′2+z-z′2-erfcx-x′2+y-y′2+z-z′2/2aτx-x′2+y-y′2+z+z′2dz′.To facilitate the expression, the nondimensional parameters are employed and the corresponding forms areΘ=k(t-t0)/ql, X=x/rb, Y=y/rb, Z=z/rb, X′=x′/rb, Y′=y′/rb, Z′=z′/rb, H=h/rb, Fo=aτ/rb2. Equation (6) is translated into(7)Θ=14π∫0HerfcX-X′2+Y-Y′2+Z-Z′2/2FoX-X′2+Y-Y′2+Z-Z′2-erfcX-X′2+Y-Y′2+Z+Z′2/2FoX-X′2+Y-Y′2+Z+Z′2dZ′.The temperature responses with the time are shown in Figure4 according to the calculation results of (7).Figure 4
The temperature responses of finite line heat source with the time.It is found that all the temperature responses will be stable at last. Furthermore, the larger the length, the fiercer the response.


### 2.2.2. Cylindrical Heat Source Model

The heat capacity inside borehole can be ignored as the diameter is extremely small compared to its depth, but the existence of inner heat transfer is undisputed. Accordingly, the line heat source model is not an ideal choice for simulating the heat transfer of BHE. The solid cylindrical heat sources have been proposed in recent years; in this model, it is supposed that the cylinder is no longer a cavity but is filled with the medium identical to that out of the cylinder [12–14]. The whole infinite domain is composed of a homogeneous medium and therefore the heat capacity of borehole is taken into consideration. The diagram of cylindrical heat source is shown in Figure 5.Figure 5
The cylindrical heat source adopted for BHE.The Green function is still the basic method to obtain the temperature response caused by the heat source, and the temperature response of any point in the underground medium except heat source is shown in(8)t-t0=ql8πk∫0τ1τ-τ′I0x2+y2rb2aτ-τ′exp⁡-x2+y2+rb24aτ-τ′·erfcz-H2aτ-τ′+erfcz+H2aτ-τ′-2erfcz2aτ-τ′dτ′,where I0 is the zero-order modified Bessel function and the borehole is supposed to be arranged at z-axis. Again, dimensionless method is used to simplify the expression and the corresponding formula is given in(9)Θ=18π∫0FodFo′Fo-Fo′I0X2+Y22Fo-Fo′exp⁡-X2+Y2+14Fo-Fo′·erfcZ-H2Fo-Fo′+erfcZ+H2Fo-Fo′-2erfcZ2Fo-Fo′.The temperature responses with the time are shown in Figure6, and the regular patterns and trends are similar with those of line heat source; that is, all responses must attain the steady states and the response value increase with the borehole length if the time is constant.Figure 6
The temperature responses of cylindrical heat source with the time.As far as the differences between line and cylindrical sources are concerned, the temperature responses changing with time are described in Figure7 when the radius and depth of borehole are constant.Figure 7
The comparison of temperature responses between line and cylindrical heat sources.It is noticed that the temperature response of cylindrical heat source is always larger than that of line heat source. The dimensionless temperature response isΘ=k(t-t0)/ql, and the heat transfer rate of line heat source becomes larger while the actual temperature rise is stable. Therefore, the calculational length of BHEs by means of line heat source is shorter than that depending on cylindrical heat source when the design of BHEs is conducted.


## 2.2.1. Line Heat Source Model

Obviously, the heat continues to be transmitted to far locations of underground medium. The heat transfer process outside BHE plays a significant role in carrying energy, because the borehole’s diameter is so small that the heat capacity inside BHE is little and the thermal resistance outside borehole is much larger than that inside borehole. The investigation on heat transfer outside borehole is becoming increasingly important and therefore the reasonable mathematical simulation models should be taken into account. The borehole’s depth is far larger than its diameter and thus BHE can be regarded as a line heat source [8], which means the radius of BHE is ignored when the heat transfer outside borehole is studied. The diagram of line heat source is described in Figure 3.Figure 3
The line heat source adopted for BHE.Accordingly, a most widely used one-dimensional model for this purpose is Kelvin’s line source model, in which the borehole is replaced by a line heat source with its radial dimension neglected. However, the actual depth of any BHE is finite and then the finite line heat source is proposed; the Kelvin’s model lays a firm foundation for the research of the finite line heat source [9]. In order to investigate the conduction around the borehole, Green’s function theory is employed and the corresponding expression is given in (5) showing the temperature response of any point (x,y,z) except heat source in the underground medium at time τ, and the point heat source locating at (x′,y′,z′) emits heat at time τ′ [10]:(5)Gx,y,z,τ;x′,y′,z′,τ′=18πaτ-τ′3/2exp⁡-x-x′2+y-y′2+z-z′24aτ-τ′.After that, the temperature response induced by the finite line heat source is obtained in (6), where erfc is the complementary error function and t0 is the initial temperature of the underground medium [11]:(6)t-t0=ql4πk∫0Herfcx-x′2+y-y′2+z-z′2/2aτx-x′2+y-y′2+z-z′2-erfcx-x′2+y-y′2+z-z′2/2aτx-x′2+y-y′2+z+z′2dz′.To facilitate the expression, the nondimensional parameters are employed and the corresponding forms areΘ=k(t-t0)/ql, X=x/rb, Y=y/rb, Z=z/rb, X′=x′/rb, Y′=y′/rb, Z′=z′/rb, H=h/rb, Fo=aτ/rb2. Equation (6) is translated into(7)Θ=14π∫0HerfcX-X′2+Y-Y′2+Z-Z′2/2FoX-X′2+Y-Y′2+Z-Z′2-erfcX-X′2+Y-Y′2+Z+Z′2/2FoX-X′2+Y-Y′2+Z+Z′2dZ′.The temperature responses with the time are shown in Figure4 according to the calculation results of (7).Figure 4
The temperature responses of finite line heat source with the time.It is found that all the temperature responses will be stable at last. Furthermore, the larger the length, the fiercer the response.


## 2.2.2. Cylindrical Heat Source Model

The heat capacity inside borehole can be ignored as the diameter is extremely small compared to its depth, but the existence of inner heat transfer is undisputed. Accordingly, the line heat source model is not an ideal choice for simulating the heat transfer of BHE. The solid cylindrical heat sources have been proposed in recent years; in this model, it is supposed that the cylinder is no longer a cavity but is filled with the medium identical to that out of the cylinder [12–14]. The whole infinite domain is composed of a homogeneous medium and therefore the heat capacity of borehole is taken into consideration. The diagram of cylindrical heat source is shown in Figure 5.Figure 5
The cylindrical heat source adopted for BHE.The Green function is still the basic method to obtain the temperature response caused by the heat source, and the temperature response of any point in the underground medium except heat source is shown in(8)t-t0=ql8πk∫0τ1τ-τ′I0x2+y2rb2aτ-τ′exp⁡-x2+y2+rb24aτ-τ′·erfcz-H2aτ-τ′+erfcz+H2aτ-τ′-2erfcz2aτ-τ′dτ′,where I0 is the zero-order modified Bessel function and the borehole is supposed to be arranged at z-axis. Again, dimensionless method is used to simplify the expression and the corresponding formula is given in(9)Θ=18π∫0FodFo′Fo-Fo′I0X2+Y22Fo-Fo′exp⁡-X2+Y2+14Fo-Fo′·erfcZ-H2Fo-Fo′+erfcZ+H2Fo-Fo′-2erfcZ2Fo-Fo′.The temperature responses with the time are shown in Figure6, and the regular patterns and trends are similar with those of line heat source; that is, all responses must attain the steady states and the response value increase with the borehole length if the time is constant.Figure 6
The temperature responses of cylindrical heat source with the time.As far as the differences between line and cylindrical sources are concerned, the temperature responses changing with time are described in Figure7 when the radius and depth of borehole are constant.Figure 7
The comparison of temperature responses between line and cylindrical heat sources.It is noticed that the temperature response of cylindrical heat source is always larger than that of line heat source. The dimensionless temperature response isΘ=k(t-t0)/ql, and the heat transfer rate of line heat source becomes larger while the actual temperature rise is stable. Therefore, the calculational length of BHEs by means of line heat source is shorter than that depending on cylindrical heat source when the design of BHEs is conducted.


## 3. The Validation of Heat Transfer Experiments

Line and cylindrical heat source models are usually adopted to simulate the heat transfer of BHE, and their rationalities can be validated by means of experiments. A hot water tank with electronic heater is employed and the power of tank is constant in the process of experiment. Figure8 shows the diagram of experiment and the BHE with U-tube is used. The circulating fluid enters the tank and is heated by electronic heater and then flows through the U-tube of borehole to exchange heat with the underground medium. The water pump drives the circulating liquid to flow circularly. Thermal sensors are set at the surface of borehole to record the temperature response with the time. The initial temperature and thermophysical properties of underground medium are known before conducting the experiment. In addition, the inlet and outlet temperatures of circulating liquid are recorded to calculate the heat transfer rate ql. There are three sensors, one is set at the depth’s middle location of the outer surface of borehole, the other two are employed to record the inlet and outlet temperatures of circulating liquid, and the accuracy of thermal sensors is ±0.1%. The diameter of BHE is 130 mm and the flow rate of fluid is around 0.3liters/s.Figure 8
Diagram of experiment of temperature response.Θ
=
k
(
t
-
t
0
)
/
q
l, and the dimensionless temperature response can be obtained based on the actual temperature recorded by the thermal sensor with the time. Then, the theoretical dimensional temperature responses of line and cylindrical heat sources can be obtained according to (7) and (9), respectively. Accordingly, the comparisons between experimental data and simulation results are shown in Figure 9.Figure 9
The comparisons between theoretical results and experimental data.
(a)
(b)The linear graph stands for the calculation results by means of simulation models and the scattergram means the experimental data. It is found that the deviation between theoretical result and experimental value is smaller if the cylindrical heat source model is employed; this is because the inner heat capacity of BHE is considered. As a result, cylindrical heat source model is more reasonable for the study of the engineering design or cost performance.


## 4. The Engineering Design for Borehole GHEs

Section2 presents the necessary theoretical knowledge of simulating heat transfer of BHE, and these contents illustrate the fundamental principles of heat transfer of BHE. There is no doubt that engineering design can put this technology into practice, and the significant design procedures and the corresponding design parameters should be emphasized to optimize the whole system. There are a number of considerations that need to be studied and they are listed in the following.


### 4.1. The Test of Thermal Property of Underground Medium

The thermal physical properties are very important because they can be used to calculate the size of BHE, and the field measurement combined with parameters estimation is employed to obtain the actual thermal physical parameters [15]. The on-site measuring devices connected with U-tube are shown in Figure 10.Figure 10
The equipment recording the parameters of circulating liquid.The object function reporting the sum of variance is established and the expression is given in(10)f=∑i=1nTcal,i-Texp,i,where Tcal,i and Texp,i are, respectively, the circulating liquid’s average temperatures of mathematical model and experimental recorded data [16]. Because the parameters of thermal physical parameters of underground medium are contained in the simulation model, the value of thermal physical parameters can be obtained by means of reverse calculation method while the value of f achieves the minimum.


### 4.2. The Cooling and Heat Load Assumed by UHEs

The heat of the building is abstracted and then is discharged into underground in summer, thereby the maximal heat release appears at the time of the maximal cooling load of buildings. The expression of the heat release is the maximal cooling load × (1 + 1/COP), where COP is the coefficient of performance of the heat pump unit.In winter, the ultimate value of heat absorption occurs while the heat load of buildings reaches the maximum, and the corresponding formula is the maximal heating load × (1 − 1/COP), where COP means the heating performance coefficient. Comparing the heat release with the heat absorption, the final total length of BHE is determined according to the large one.


### 4.3. The Calculation of Pressure Drop

The purpose of calculating pressure drop is to choose the water pump so that the liquid can achieve circulation flow. As far as the BHEs are concerned, to obtain the linear head loss and partial resistance is the basis of determining the circulating pump.W,A, and d are, respectively, the rate of flow, area of cross section, and inner diameter of tube; therefore, the velocity V is obtained according to(11)V=W3600×A=4W3600×π×d2.Because the turbulence is prerequisite while the liquid flows during the running process, the corresponding Reynolds number, that is,Re, should be larger than 2300 to ensure this state. The expression of Re is shown in(12)Re=ρVdμ,where ρ and μ are, respectively, the density and coefficient of liquid’s dynamic viscosity and the calculation of Re is significant before studying the head loss. Next, (13) offers the linear heat loss and it involves several parameters:(13)PL=0.158×ρ0.75×μ0.25×d1.25×V1.75×L,where L denotes the length of tube of liquid.Partial loss usually occurs when fluid flows through the elbows, equipment, and so on, and (14) shows the detailed expression(14)Pp=0.158×ρ0.75×μ0.25×d1.25×V1.75×E,where E means the equivalent length of the elbows, equipment, and so on.Afterwards, the total pressure drop is equal to the sum of linear head loss and partial loss [17, 18].


### 4.4. The Joining of BHEs

It is unquestionable that there are a number of BHEs for engineering project, which means BHEs are distributed row after row or column after column; the connection of BHEs is a technology that importance should be attached to. The subcatchment devices are commonly employed to connect a group of BHEs and there are many groups. Then, all distributors are linked to the room of heat pump. The liquids of a group of BHEs are connected by collecting breeching and then enters the room of heat pump. Later, the liquid returns and then is distributed by the separator to flow through tubes of boreholes in parallel. The material of U-tube is high-density polyethylene and this material has favorable conductivity, corrosion resistance, and long service life. The diagram describing the connection of a group of BHEs is shown in Figure11.Figure 11
The connection of a group of BHEs.As a rule, the external diameter of every U-tube is usually 25 mm or 32 mm, and the HDPE tubes which have the external diameter of 50 mm or 63 mm are often used as the main tube of a group of BHEs.


### 4.5. The Backfilling Material inside Borehole

Backfilling material lies in the area between U-tube and borehole wall, and it enhances the heat transfer and stops the surface water from groundwater permeation. The reasonable backfilling material guarantees the heat transfer performance of BHEs. The heat resistance inside borehole is inevitably increased in case the backfilling material has undesirable conduction performance; therefore, the total lengths of BHEs are lengthened and both investment and operating costs are added. Researchers have paid more and more attention to the improvement of backfilling material with the development of GSHP technology.The main characteristics of backfilling material consist of conduction coefficient, uniformity, stability, and antihot pressing [19–21]. Cement, sand, and bentonite are usually employed to compose the material and the corresponding conduction performance depends on the ratios of these elements. The overseas research indicated that the heat conductivity can arrive at 2.0 W/(m·K), and in recent years the civil experts or scholars found that the value can be further improved and even achieved 2.1 W/(m·K) [22]. The influence that backfilling material exerts on the length of BHEs is obvious even if the value of heat conductivity is increased a little.


## 4.1. The Test of Thermal Property of Underground Medium

The thermal physical properties are very important because they can be used to calculate the size of BHE, and the field measurement combined with parameters estimation is employed to obtain the actual thermal physical parameters [15]. The on-site measuring devices connected with U-tube are shown in Figure 10.Figure 10
The equipment recording the parameters of circulating liquid.The object function reporting the sum of variance is established and the expression is given in(10)f=∑i=1nTcal,i-Texp,i,where Tcal,i and Texp,i are, respectively, the circulating liquid’s average temperatures of mathematical model and experimental recorded data [16]. Because the parameters of thermal physical parameters of underground medium are contained in the simulation model, the value of thermal physical parameters can be obtained by means of reverse calculation method while the value of f achieves the minimum.


## 4.2. The Cooling and Heat Load Assumed by UHEs

The heat of the building is abstracted and then is discharged into underground in summer, thereby the maximal heat release appears at the time of the maximal cooling load of buildings. The expression of the heat release is the maximal cooling load × (1 + 1/COP), where COP is the coefficient of performance of the heat pump unit.In winter, the ultimate value of heat absorption occurs while the heat load of buildings reaches the maximum, and the corresponding formula is the maximal heating load × (1 − 1/COP), where COP means the heating performance coefficient. Comparing the heat release with the heat absorption, the final total length of BHE is determined according to the large one.


## 4.3. The Calculation of Pressure Drop

The purpose of calculating pressure drop is to choose the water pump so that the liquid can achieve circulation flow. As far as the BHEs are concerned, to obtain the linear head loss and partial resistance is the basis of determining the circulating pump.W,A, and d are, respectively, the rate of flow, area of cross section, and inner diameter of tube; therefore, the velocity V is obtained according to(11)V=W3600×A=4W3600×π×d2.Because the turbulence is prerequisite while the liquid flows during the running process, the corresponding Reynolds number, that is,Re, should be larger than 2300 to ensure this state. The expression of Re is shown in(12)Re=ρVdμ,where ρ and μ are, respectively, the density and coefficient of liquid’s dynamic viscosity and the calculation of Re is significant before studying the head loss. Next, (13) offers the linear heat loss and it involves several parameters:(13)PL=0.158×ρ0.75×μ0.25×d1.25×V1.75×L,where L denotes the length of tube of liquid.Partial loss usually occurs when fluid flows through the elbows, equipment, and so on, and (14) shows the detailed expression(14)Pp=0.158×ρ0.75×μ0.25×d1.25×V1.75×E,where E means the equivalent length of the elbows, equipment, and so on.Afterwards, the total pressure drop is equal to the sum of linear head loss and partial loss [17, 18].


## 4.4. The Joining of BHEs

It is unquestionable that there are a number of BHEs for engineering project, which means BHEs are distributed row after row or column after column; the connection of BHEs is a technology that importance should be attached to. The subcatchment devices are commonly employed to connect a group of BHEs and there are many groups. Then, all distributors are linked to the room of heat pump. The liquids of a group of BHEs are connected by collecting breeching and then enters the room of heat pump. Later, the liquid returns and then is distributed by the separator to flow through tubes of boreholes in parallel. The material of U-tube is high-density polyethylene and this material has favorable conductivity, corrosion resistance, and long service life. The diagram describing the connection of a group of BHEs is shown in Figure11.Figure 11
The connection of a group of BHEs.As a rule, the external diameter of every U-tube is usually 25 mm or 32 mm, and the HDPE tubes which have the external diameter of 50 mm or 63 mm are often used as the main tube of a group of BHEs.


## 4.5. The Backfilling Material inside Borehole

Backfilling material lies in the area between U-tube and borehole wall, and it enhances the heat transfer and stops the surface water from groundwater permeation. The reasonable backfilling material guarantees the heat transfer performance of BHEs. The heat resistance inside borehole is inevitably increased in case the backfilling material has undesirable conduction performance; therefore, the total lengths of BHEs are lengthened and both investment and operating costs are added. Researchers have paid more and more attention to the improvement of backfilling material with the development of GSHP technology.The main characteristics of backfilling material consist of conduction coefficient, uniformity, stability, and antihot pressing [19–21]. Cement, sand, and bentonite are usually employed to compose the material and the corresponding conduction performance depends on the ratios of these elements. The overseas research indicated that the heat conductivity can arrive at 2.0 W/(m·K), and in recent years the civil experts or scholars found that the value can be further improved and even achieved 2.1 W/(m·K) [22]. The influence that backfilling material exerts on the length of BHEs is obvious even if the value of heat conductivity is increased a little.


## 5. The Analysis on the Actual Projects of GSHP System

Although Sections2 and 3, respectively, present the theoretical simulation models and engineering design of BHEs, the practicability, economical efficiency, and energy-saving of GSHP technology should be validated by actual engineering projects.


### 5.1. The Project Information

This building is a library and it belongs to public construction, the floor and air-conditioning areas are, respectively, 26800 m2 and 21300 m2, and the height is nearly 30 m. Table 1 lists the information of the areas and air-conditioning loads.Table 1
The areas and corresponding loads.
Building nature
Building area(m2)
Air-conditioning area(m2)
Cooling load(kW)
Heating load(kW)
Library
26800
21300
2378
2308


### 5.2. The First Design Plan

Both cooling and heat loads are assumed by BHEs so that the unity of cold and heat sources is achieved. There are two heat pump units that are responsible for cooling and heating and the corresponding parameters are given in Table2.Table 2
Parameters of the heat pump unit of GSHP system.
Operating mode
Cooling/heatingcapacity(kW)
Temperature of U-tube circulating liquid(°C)
Temperature of cooling water(°C)
Pressure drop of evaporator(kPa)
Pressure drop of condenser(kPa)
Electricity consumption(kW)
Cooling
1130
12/7
15/23
98
13.7
193
Heating
1473
50/45
15/10
98
13.7
336The total length of BHEs is 36640 m according to the analysis and calculation; there are 458 boreholes with the diameter of 140 mm and every one has the depth of 80 m, and the row and column intervals are both 4 m.


### 5.3. The Second Design Plan

The direct-fired lithium-bromide refrigerating machine (DFLBRM) is employed for cooling and heating. The cooling tower is set in the open air to act as the heat discharge equipment for the cooling condition, and the gas is fired to become the driving heat source so that the heating objective is achieved. The characteristics of this plan can be summarized as several points: firstly, it can save electricity but not save energy; secondly, the cold and heat sources are united; thirdly, the requirements of the automatic control and safety are very high.According to the cooling and heat loads, two DFLBRMs are selected and the corresponding parameters are listed in Table3.Table 3
Parameters of the direct-fired lithium-bromide refrigerating machine.
Cooling capacity (kW)
Heating capacity (kW)
Cold water
Cooling water
The maximum amount of fuel
Distribution capacity
Solution capacity
Running weight
Rate of flow(m3/h)
Pressure drop(MPa)
Diameter of joint pipe(mm)
Rate of flow(m3/h)
Pressure drop(MPa)
Diameter of joint pipe(mm)
Cooling 104 kcal
Heating104 kcal
kW
ton
kg
1454
1121
250
0.07
200
304
0.11
250
94.6
104
9.1
4.1
20.4


### 5.4. Comparisons of Costs between the Planes

Tables4 and 5, respectively, describe the investments of the first and the second design planes.Table 4
Initial cost of GSHP system.
Equipment
Number
Unit price
Sum price
(¥104)
(¥104)
Heat pump units
2
119.365
238.73
Cold water pump
2
1.55
3.1
Cooling water pump
2
1.55
3.1
Water separator(outdoor circulation)
1
0.98
0.98
Water collector(outdoor circulation)
1
0.98
0.98
Barrier diaphragm-type pressurized water tank
1
0.69
0.69
Installation of UHEs
492.85
Total
740.43Table 5
Initial cost of the direct-fired lithium-bromide refrigerating machine system.
Equipment
Number
Unit price
Sum price
(¥104)
(¥104)
DFLBRM
2
135
270
Cooling tower
2
11.05
22.1
Cold water pump
2
1.55
3.1
Cooling water pump
2
1.55
3.1
Gas
29779.9
0.0025
74.4498
Gas pipe laying  work
100
Total
472.75In this project, the thermal conductivity of backfilling material is 1.5 W/(m·K) and the corresponding installation cost of BHEs is ¥492.85 × 104. The thermal conductivity can exert influence on the cost, and the corresponding cost decreases with the increase of thermal conductivity, the detailed information is shown in Figure 12.Figure 12
The installation cost of GHEs changes with the thermal conductivity of backfilling material.It should be admitted that the initial cost of GSHP is obviously higher than that of DFLBRMs, but the running cost of these two different modes can embody the advantage of GSHP technology. And the initial cost of GSHP can be further reduced while the thermal conductivity of backfilling material is increased, and the installation cost of BHEs is obtained according to the minimal value of conductivity. The detailed information of operating costs is shown in Tables6 and 7.Table 6
Operating cost of GSHP system.
Equipment
Number
Electricityconsumption
Runningtime
Energyconsumption
Unit price
Operating cost per day
Operating cost per month
Operating months
Operating cost per  year
kW
h
kWh
¥/kWh
¥104
¥104
¥104
Heat pump units
Cooling2
193
12
4632
0.6
0.27792
8.3376
3
25.0128
Heating2
336
12
8064
0.6
0.48384
14.5152
4
58.0608
Cold water pump
2
37
12
888
0.6
0.05328
1.5984
7
11.1888
Cooling    water pump
2
37
12
888
0.6
0.05328
1.5984
7
11.1888
Total
105.4512Table 7
Operating cost of DFLBRMs system.
Number
Electricity (gas)consumption
Running time
Energy consumption
Unit price
Operating cost per day
Operating cost per month
Operating months
Operating cost per year
kW (m3/h)
h
kW
¥
¥104
¥104
¥104
DFLBRMs
Cooling2
9.1
12
218.4
0.6
0.013104
0.39312
3
1.17936
Heating2
9.1
12
218.4
0.6
0.013104
0.39312
4
1.57248
Cooling tower
2
25
12
600
0.6
0.036
1.08
3
3.24
Cold water pump
2
37
12
888
0.6
0.05328
1.5984
3
4.7952
Cooling water pump
2
37
12
888
0.6
0.05328
1.5984
7
11.1888
Gas cooling
267
12
2.2
0.70488
21.1464
3
63.4392
Heating
288
12
2.2
0.76032
22.8096
4
91.2384
Total
176.6534The significant data of tables make clear that the initial cost of GSHP system is expensive; however, the favorable benefits in terms of energy-saving and environmental protection are certainly achieved [23]. In addition, the most attention-getting point is that the running cost is obviously lower and the payback period is less than 4 years compared with the DFLBRMs. What is more, the installation cost of BHEs can be reduced while the thermal conductivity of backfilling material increases, and this can further improve the economic performance of the whole GSHP system.


## 5.1. The Project Information

This building is a library and it belongs to public construction, the floor and air-conditioning areas are, respectively, 26800 m2 and 21300 m2, and the height is nearly 30 m. Table 1 lists the information of the areas and air-conditioning loads.Table 1
The areas and corresponding loads.
Building nature
Building area(m2)
Air-conditioning area(m2)
Cooling load(kW)
Heating load(kW)
Library
26800
21300
2378
2308


## 5.2. The First Design Plan

Both cooling and heat loads are assumed by BHEs so that the unity of cold and heat sources is achieved. There are two heat pump units that are responsible for cooling and heating and the corresponding parameters are given in Table2.Table 2
Parameters of the heat pump unit of GSHP system.
Operating mode
Cooling/heatingcapacity(kW)
Temperature of U-tube circulating liquid(°C)
Temperature of cooling water(°C)
Pressure drop of evaporator(kPa)
Pressure drop of condenser(kPa)
Electricity consumption(kW)
Cooling
1130
12/7
15/23
98
13.7
193
Heating
1473
50/45
15/10
98
13.7
336The total length of BHEs is 36640 m according to the analysis and calculation; there are 458 boreholes with the diameter of 140 mm and every one has the depth of 80 m, and the row and column intervals are both 4 m.


## 5.3. The Second Design Plan

The direct-fired lithium-bromide refrigerating machine (DFLBRM) is employed for cooling and heating. The cooling tower is set in the open air to act as the heat discharge equipment for the cooling condition, and the gas is fired to become the driving heat source so that the heating objective is achieved. The characteristics of this plan can be summarized as several points: firstly, it can save electricity but not save energy; secondly, the cold and heat sources are united; thirdly, the requirements of the automatic control and safety are very high.According to the cooling and heat loads, two DFLBRMs are selected and the corresponding parameters are listed in Table3.Table 3
Parameters of the direct-fired lithium-bromide refrigerating machine.
Cooling capacity (kW)
Heating capacity (kW)
Cold water
Cooling water
The maximum amount of fuel
Distribution capacity
Solution capacity
Running weight
Rate of flow(m3/h)
Pressure drop(MPa)
Diameter of joint pipe(mm)
Rate of flow(m3/h)
Pressure drop(MPa)
Diameter of joint pipe(mm)
Cooling 104 kcal
Heating104 kcal
kW
ton
kg
1454
1121
250
0.07
200
304
0.11
250
94.6
104
9.1
4.1
20.4


## 5.4. Comparisons of Costs between the Planes

Tables4 and 5, respectively, describe the investments of the first and the second design planes.Table 4
Initial cost of GSHP system.
Equipment
Number
Unit price
Sum price
(¥104)
(¥104)
Heat pump units
2
119.365
238.73
Cold water pump
2
1.55
3.1
Cooling water pump
2
1.55
3.1
Water separator(outdoor circulation)
1
0.98
0.98
Water collector(outdoor circulation)
1
0.98
0.98
Barrier diaphragm-type pressurized water tank
1
0.69
0.69
Installation of UHEs
492.85
Total
740.43Table 5
Initial cost of the direct-fired lithium-bromide refrigerating machine system.
Equipment
Number
Unit price
Sum price
(¥104)
(¥104)
DFLBRM
2
135
270
Cooling tower
2
11.05
22.1
Cold water pump
2
1.55
3.1
Cooling water pump
2
1.55
3.1
Gas
29779.9
0.0025
74.4498
Gas pipe laying  work
100
Total
472.75In this project, the thermal conductivity of backfilling material is 1.5 W/(m·K) and the corresponding installation cost of BHEs is ¥492.85 × 104. The thermal conductivity can exert influence on the cost, and the corresponding cost decreases with the increase of thermal conductivity, the detailed information is shown in Figure 12.Figure 12
The installation cost of GHEs changes with the thermal conductivity of backfilling material.It should be admitted that the initial cost of GSHP is obviously higher than that of DFLBRMs, but the running cost of these two different modes can embody the advantage of GSHP technology. And the initial cost of GSHP can be further reduced while the thermal conductivity of backfilling material is increased, and the installation cost of BHEs is obtained according to the minimal value of conductivity. The detailed information of operating costs is shown in Tables6 and 7.Table 6
Operating cost of GSHP system.
Equipment
Number
Electricityconsumption
Runningtime
Energyconsumption
Unit price
Operating cost per day
Operating cost per month
Operating months
Operating cost per  year
kW
h
kWh
¥/kWh
¥104
¥104
¥104
Heat pump units
Cooling2
193
12
4632
0.6
0.27792
8.3376
3
25.0128
Heating2
336
12
8064
0.6
0.48384
14.5152
4
58.0608
Cold water pump
2
37
12
888
0.6
0.05328
1.5984
7
11.1888
Cooling    water pump
2
37
12
888
0.6
0.05328
1.5984
7
11.1888
Total
105.4512Table 7
Operating cost of DFLBRMs system.
Number
Electricity (gas)consumption
Running time
Energy consumption
Unit price
Operating cost per day
Operating cost per month
Operating months
Operating cost per year
kW (m3/h)
h
kW
¥
¥104
¥104
¥104
DFLBRMs
Cooling2
9.1
12
218.4
0.6
0.013104
0.39312
3
1.17936
Heating2
9.1
12
218.4
0.6
0.013104
0.39312
4
1.57248
Cooling tower
2
25
12
600
0.6
0.036
1.08
3
3.24
Cold water pump
2
37
12
888
0.6
0.05328
1.5984
3
4.7952
Cooling water pump
2
37
12
888
0.6
0.05328
1.5984
7
11.1888
Gas cooling
267
12
2.2
0.70488
21.1464
3
63.4392
Heating
288
12
2.2
0.76032
22.8096
4
91.2384
Total
176.6534The significant data of tables make clear that the initial cost of GSHP system is expensive; however, the favorable benefits in terms of energy-saving and environmental protection are certainly achieved [23]. In addition, the most attention-getting point is that the running cost is obviously lower and the payback period is less than 4 years compared with the DFLBRMs. What is more, the installation cost of BHEs can be reduced while the thermal conductivity of backfilling material increases, and this can further improve the economic performance of the whole GSHP system.


## 6. Conclusions

The vertical BHEs are typical underground components that are in charge of heat release and heat abstraction; the paper demonstrates the fundamental principles of heat transfer inside and outside borehole. The research on theoretical knowledge lays a firm foundation for the engineering design of GSHP projects. The construction of underground installation is complicated and therefore a series of technological process are explained in detail, including steps and parameters. The discussions of theoretical heat transfer and actual engineering knowledge enable researchers to understand how to make full use of GSHP. In addition, an actual engineering project that employs GSHP system is analyzed; as far as GSHP and DFLBRMs are concerned, not only the investments but also the operating costs are compared with each other. By means of comparisons, it should be highlighted that the initial expense of GSHP is higher but the running cost is lower, which means the added part of initial cost can be saved by way of saving running cost. The initial cost of GSHP can be reduced further while the thermal conductivity of backfilling material is increased, and this can let the cost performance of GSHP be improved further. According to the detailed illustration, GSHP will make greater progress and gain more and more attention in the future.



---

*Source: 870380-2015-10-18.xml*