# Assessment of Energy Potential and Benefit under the Ground Source Heat Pump Air Conditioning System in Anqing Area of Eastern China

Authors: Qing Zhang; Yunfeng Li; Jian Hua; Xiaonan Niu; Lili Hou; Zongfang Chen; Weiya Ge
Journal: Geofluids
 (2021)

---

## Abstract

Now ground source heat pump is a more efficient way to develop and utilize shallow geothermal energy because it is clean and environmentally friendly and has a relatively low energy cost. In order to optimize the planning layout and geographical space development in the eastern new town of Anqing city, which can realize the transformation and upgrading of real space and urban sustainable development, the exploration for shallow geothermal energy will be carried out in this area, so as to find out the comprehensive thermophysical parameters of the shallow rock-soil body and the heat transfer capacity of the vertical heat exchanger, etc. This paper takes the CBD in the eastern new town of Anqing as an example to provide the basis for the feasibility construction of the ground source heat pump project in the study area and evaluate the economic and environmental benefits of the expected project. According to the simulation test data of 5 working conditions of 4 geothermal exploration holes in the study area, we can clearly know that the energy cost per square meter of the ground source heat pump is 11.8 yuan for a building of one hundred thousand square meters in which the heat removal power is expected to be 9481 kW in summer and 3070 kW in winter. And the annual emission of carbon dioxide, sulfur dioxide, nitrogen oxides, suspended dust, and other air pollutants to the atmosphere can be reduced by 1442.5 t, and the solid waste ash and slag can be reduced by 59.7 t. The annual environmental treatment cost will be saved by 166000 yuan.


---

## Body

## 1. Introduction

Geothermal is a clean and renewable resource, which can be used as an effective substitute for fossil fuels [1]. Among the various technologies and methods that can be used to obtain geothermal energy, ground source heat pump (GSHP) system is considered an efficient and main means [2].The concept of “ground source heat pump” was first put forward by Swiss experts in 1912, and this technology was first put forward in Britain and the United States. Nordic countries mainly focus on winter heating, while the United States focuses on combined winter and summer heating. In recent years, the application of ground source heat pump technology for heating in China has been favored, and the application of ground source heat pump system in 2016 has been nearly six times larger than that in 2007 [3, 4]. Ground source heat pump system is considered an efficient geothermal energy acquisition technology because of its high efficiency and low operating cost, but its high initial cost limits the development of this technology. Many scholars have done a lot of research work on the ground source heat pump system in consideration of the performance [5, 6] and economy [7]. In terms of the performance of the ground source heat pump, Nouri et al. [8] suggested the vacuum tube collector and VGHP should be combined to meet the heating, cooling, and hot water demand of a residential building in Tabriz. Balbay and Esen [9, 10] carried out a preliminary feasibility study on the application of VGHP in snow melting of pavement and bridge deck for the first time. The results show that the snow on the pavement and bridge deck is effectively melted. For the economy of ground source heat pump, Esen et al. [5] reported the comparison between ground source heat pump system and ACHP system from the perspective of technology and economy. Lu et al. [11] used some economic indicators to compare the financial attractiveness of GSHP and ASHP in residential heating and cooling in Melbourne, and the results showed that the ground source heat pump system provided more benefits than the ASHP system for systems with more than 40 years of operation. Ground source heat pump system can extract geothermal energy for building heating and cooling. Compared with the traditional air-water heat pump (AWHP) and air-air heat pump (AAHP) system, the thermal performance of ground source heat pump system is usually higher, because the ground temperature is higher than the ambient air temperature in the heating season and lower than the ambient air temperature in the cooling season [12]. The ground source heat pump (GSHP) is a central heating and air conditioning system, which is composed of a water ground source heat pump unit, geothermal energy exchange system, and building system, with the rock-soil body, stratum soil, groundwater, or surface water as a low-temperature heat source. Therefore, it is very important to accurately grasp the thermal conductivity and other parameters of the heat source in situ for the evaluation and construction of a long-term effective operation geothermal heat pump system.Shallow geothermal energy investigation and ground source heat pump suitability zoning research have been carried out in many cities in China [13, 14]. However, due to the differences in geotechnical conditions, groundwater conditions, and climatic conditions in different regions, the resource characteristics of different regions have their unique characteristics [15, 16]. Moreover, the research content mostly focuses on the distribution characteristics of shallow geothermal resources, and the research on environmental benefit evaluation is less [17, 18]. In this paper, the self-developed equipment called FTPT11 for field thermal response test was used in the four geothermal wells, carrying out the simulation test of five working conditions, including the initial mean ground temperature test (case 1), high-power constant heat flux test (case 2), low-power constant heat flux test (case 3), simulated constant temperature test in summer (case 4), and simulated constant temperature test in winter (case 5). In this way, the development potential of shallow geothermal energy will be analyzed and calculated. What is more, the economic and environmental benefits of the proposed ground source heat pump system may be evaluated to provide a reference for later engineering construction.


## 2. Summary

### 2.1. Project Summary

The study area is located in the central business district (CBD) of the eastern new town of Anqing City, Anhui Province, Southeast China, with a total area of 1812 mu, as shown in Figure1. Anqing is located in the middle and low latitudes with four distinct seasons. It is hot in summer and cold in winter. In order to realize the sustainable development of the new urban area, it is planned to make large-scale use of “ground source heat pump air conditioning” technology to reduce the proportion of fossil fuels used in the new urban area.Figure 1
Geographical location and general situation of the study area.At present, there are four geothermal exploration boreholes in the study area, namely, 192k11, 192k12, 192k13, and 192k14. A large number of field tests and laboratory tests were carried out to obtain the formation parameters and geothermal distribution characteristics in the study area, so as to study the potential of vertical buried pipe heat exchange system.


### 2.2. Stratigraphic System

The strata in Anqing city are divided into two stratigraphic areas by the Tan Lu fault. The Dabie Mountains stratigraphic area is in the northwest while the Anqing stratigraphic area is in the southeast. The former consists of the metamorphic rocks of the Dabie Mountains group, Upper Archean, the metamorphic rocks of the Foziling group, Lower Proterozoic, and the middle basic volcanic rocks of the Upper Jurassic of the Mesozoic. The total thickness of the latter strata is more than 12000 m, but the outcrop is not good. It is sporadically distributed in Silurian, Devonian, Carboniferous, Permian, Triassic, Jurassic, Cretaceous, and tertiary. The Quaternary strata are widely distributed; the lithologic characteristics are shown in Figure2.Figure 2
General situation of strata in the study area.The rock of the Xuannan group, Upper Cretaceous (K2x), is purple-red; it covers medium-thick conglomerate, fine conglomerate, gravel bearing coarse sandstone, and argillaceous siltstone, which are interbedded. The gravel grain size is 1-50 mm, with general sorting and roundness. The gravels have a bedding structure and are arranged along the dip direction.In this study, cores of K11 and K14 geothermal holes were obtained for laboratory tests and analysis (Figure3). The results show that the average thermal conductivity of the hole K11 is 1.104~1.997 W/(m·K), and the average specific heat capacity is 0.95~1.697 kJ/(kg·K), while the average thermal conductivity of the hole K14 is 1.107~1.755 W/(m·K), and the average specific heat capacity is 1.074~1.589 kJ/(kg·K). On the whole, except for a few points, the specific heat capacity of the rock-soil body in the range of drilling depth has little change with depth, and the thermal conductivity has no obvious change with depth (Figure 4). Furthermore, due to the fact that holes K11 and K12 are near while holes K13 and K14 are near, the lithology of the stratum is consistent. Therefore, the average thermal conductivity of the whole hole section calculated by thickness weighting in the study area is about 1.575 W/(m·K), which is suitable for geothermal exploitation.Figure 3
Core drilling of K11 and K14 (a) and testing process (b).Figure 4
Variation curve of specific heat capacity (a) and thermal conductivity (b) of the rock-soil body around K11 and K14 with depth.


## 2.1. Project Summary

The study area is located in the central business district (CBD) of the eastern new town of Anqing City, Anhui Province, Southeast China, with a total area of 1812 mu, as shown in Figure1. Anqing is located in the middle and low latitudes with four distinct seasons. It is hot in summer and cold in winter. In order to realize the sustainable development of the new urban area, it is planned to make large-scale use of “ground source heat pump air conditioning” technology to reduce the proportion of fossil fuels used in the new urban area.Figure 1
Geographical location and general situation of the study area.At present, there are four geothermal exploration boreholes in the study area, namely, 192k11, 192k12, 192k13, and 192k14. A large number of field tests and laboratory tests were carried out to obtain the formation parameters and geothermal distribution characteristics in the study area, so as to study the potential of vertical buried pipe heat exchange system.


## 2.2. Stratigraphic System

The strata in Anqing city are divided into two stratigraphic areas by the Tan Lu fault. The Dabie Mountains stratigraphic area is in the northwest while the Anqing stratigraphic area is in the southeast. The former consists of the metamorphic rocks of the Dabie Mountains group, Upper Archean, the metamorphic rocks of the Foziling group, Lower Proterozoic, and the middle basic volcanic rocks of the Upper Jurassic of the Mesozoic. The total thickness of the latter strata is more than 12000 m, but the outcrop is not good. It is sporadically distributed in Silurian, Devonian, Carboniferous, Permian, Triassic, Jurassic, Cretaceous, and tertiary. The Quaternary strata are widely distributed; the lithologic characteristics are shown in Figure2.Figure 2
General situation of strata in the study area.The rock of the Xuannan group, Upper Cretaceous (K2x), is purple-red; it covers medium-thick conglomerate, fine conglomerate, gravel bearing coarse sandstone, and argillaceous siltstone, which are interbedded. The gravel grain size is 1-50 mm, with general sorting and roundness. The gravels have a bedding structure and are arranged along the dip direction.In this study, cores of K11 and K14 geothermal holes were obtained for laboratory tests and analysis (Figure3). The results show that the average thermal conductivity of the hole K11 is 1.104~1.997 W/(m·K), and the average specific heat capacity is 0.95~1.697 kJ/(kg·K), while the average thermal conductivity of the hole K14 is 1.107~1.755 W/(m·K), and the average specific heat capacity is 1.074~1.589 kJ/(kg·K). On the whole, except for a few points, the specific heat capacity of the rock-soil body in the range of drilling depth has little change with depth, and the thermal conductivity has no obvious change with depth (Figure 4). Furthermore, due to the fact that holes K11 and K12 are near while holes K13 and K14 are near, the lithology of the stratum is consistent. Therefore, the average thermal conductivity of the whole hole section calculated by thickness weighting in the study area is about 1.575 W/(m·K), which is suitable for geothermal exploitation.Figure 3
Core drilling of K11 and K14 (a) and testing process (b).Figure 4
Variation curve of specific heat capacity (a) and thermal conductivity (b) of the rock-soil body around K11 and K14 with depth.


## 3. Thermal Response Test

### 3.1. Installation of Vertical Heat Exchanger

Four groups of vertical heat exchangers were arranged in this survey, all of which were installed after secondary reaming of the coring hole (Figure5). The depth of coring in holes 192K11 and 192K12 is 60 m, and double-U-shaped PE pipe heat exchangers of 54 m and 60 m are installed, respectively. Besides, the depth of coring in holes 192K13 and 192K14 is 100 m, and 100 m double-U-shaped PE pipe heat exchangers are installed. The PE pipe should be left about 1 m long on the ground considering the need of connecting the thermal response test instrument.Figure 5
Installation of vertical heat exchanger.


### 3.2. Thermal Response Test

#### 3.2.1. Introduction to the Test

The main performance parameters of the equipment (Figure6) which is self-developed named FTPT11 for field thermal response test are shown in Table 1. The test equipment is mainly composed of a heat pump system, electric heater, water-replenishing tank, circulating pump, circulating pipeline, temperature and flow detection elements, etc. The flow meter is installed at the water supply pipe of the vertical buried pipe, and the temperature detection points are arranged at the water supply pipe, return pipe, and circulating pipe, which is used to automatically record the flow and temperature of circulating water in the pipe.Figure 6
Main machine (a) and refrigeration auxiliary machine (b).Table 1
Performance parameters of test equipment.
Equipment modelMaximum heating capacityMaximum cooling capacityMaximum powerCirculating water flowRangeFTPT118 kW7 kW12 kW1.8 m3/h10 mFour groups of thermal response tests have been completed in this exploration, and each group includes the initial mean ground temperature test (case 1), high-power constant heat flux test (case 2), low-power constant heat flux test (case 3), simulated constant temperature test in summer (case 4), and simulated constant temperature test in winter (case 5).The purpose of the low-power constant heat flux test and high-power constant heat flux test is to obtain the comprehensive thermal conductivity of the rock-soil body in the depth range of the vertical heat exchanger. The purpose of the simulated constant temperature test in summer and simulated constant temperature test in winter is to obtain the heat transfer power per meter length of the vertically buried tube heat exchanger under specific heat transfer conditions. Combined with the initial average ground temperature test results, the relationship between heat transfer power per unit length and medium temperature of vertical heat exchanger is further fitted, and the effective heat transfer coefficient of vertically buried tube heat exchanger is obtained.


#### 3.2.2. Results and Analysis

(1) The Initial Average Temperature of the Rock-Soil Body (Case 1). Figure 7 shows the test curve of the average initial temperature of the rock-soil body in the four geothermal holes. After a certain period of system circulation, the final stable average water temperature (initial average temperature of the rock-soil body) at the outlet of four test holes is 21.9°C, 20.8°C, 20.5°C, and 20.5°C, respectively.Figure 7
Initial average temperature of the rock-soil body in 4 geothermal test holes.(2) High- and Low-Power Constant Heat Flux Test (Case 2 and Case 3). The test duration, heating power, and flow rate of each test hole are 48 h, 1.910 kW, and 1.117 m3/h, respectively, in the low-power constant heat flux test. According to the test data, the curve of the temperature at the inlet and outlet of the buried pipe with time is drawn as shown in Figure 8. Each temperature measurement curve can be divided into two obvious stages: rapid temperature rise (stage I) and stable temperature (stage II). Stage I is in the first ten hours of the test, the temperature at the inlet and outlet of the buried pipe changes rapidly, and the temperature in the borehole (backfill material) rises rapidly. Stage II shows that the temperature change in the borehole decreases with the heating time, which can be regarded as a steady-state heat conduction process.Figure 8
Results of low-power constant heat flux test in four geothermal holes.The linear source model proposed by Hellström [16] is selected to calculate the thermal conductivity. The mean temperature of the circulating medium at time τ can be expressed as follows:
(1)Tf=Tff+ql⋅Rb+14πλs⋅Eidb2ρscs16λsτ,where Tf is the mean temperature of circulating medium (°C), Tff is the undisturbed ground temperature (°C), q1 is the heat flux per length (W), ρs is the mean density of the ground surrounding the BHE (kg/m3), τ is the time (s), cs is the mean specific heat capacity of the ground surrounding the BHE (J/(kg K)), db is the borehole diameter (m), and λs is the thermal conductivity of the ground surrounding the borehole (W/(m K)).When the heating time is sufficiently long, the following equation is obtained:(2)Eidb2ρscs16λsτ≈ln16λsτdb2ρscs−γ.Linear regression is presented between the borehole wall temperature and logarithmic time. Equation (4) can be simplified to
(3)Tf=mlnt+n,where m is the slope of the linear regression between temperature and logarithmic time.The equations determine the heat exchange between the circulating medium and the surrounding ground. Thus, the ground thermal conductivityλs can be determined. According to the test data of the temperature stable section of each test hole, the curve of constant heat flow heating Tf versus lnt of each test hole is drawn as Figure 9.Figure 9
The variation curve of low-power constant heat flow heatingTf with lnt in four geothermal test holes.The test duration, heating power, and flow rate of each test hole are 48 h, 3.725 kW, and 1.046 m3/h, respectively, in the high-power constant heat flux test. According to the test data, draw the curve of the temperature at the inlet and outlet of the buried pipe with time, as shown in Figure 10. It can be seen from the curve that the shape of the curve is basically the same as that of the low-power constant heat flow test. In the first ten hours, the temperature changes at the inlet and outlet of the buried pipe are relatively fast, and the temperature rises rapidly in the borehole (backfill material). With continuous heating, the range of temperature change is reduced in the hole decreases, and the heat transfer process in the hole is similar to the steady-state heat conduction process.Figure 10
Results of high-power constant heat flux test in four geothermal holes.According to the test data of the temperature stable section of each test hole, the curve of constant heat flow heatingTf with lnt is drawn in Figure 11, and the average thermal conductivity of the rock-soil body is calculated.Figure 11
The variation curve of high-power constant heat flow heatingTf with lnt in four geothermal test holes.In conclusion, the average thermal conductivity of the rock-soil body at the four exploration holes (Table2) is calculated based on the results of the high-power and low-power constant heat flow test, so as to predict the effect of geothermal exploitation in the study area more accurately. The results show that the average thermal conductivity of K11 and K12 is close, while the value of K13 and K14 is close. This is because the holes K11 and K12 are near while holes K13 and K14 are near, and the difference of lithology between adjacent boreholes is small.Table 2
Average thermal conductivity of rock-soil body under the high- and low-power constant heat flux test.
Test type192K11192K12192K13192K14Low-power constant heat1.93 W/(m·°C)1.80 W/(m·°C)2.04 W/(m·°C)2.10 W/(m·°C)High-power constant heat1.90 W/(m·°C)1.80 W/(m·°C)2.08 W/(m·°C)2.08 W/(m·°C)(3) Simulated Constant Temperature Test in Summer and Winter (Case 4 and Case 5). Because the ground source heat pump project is used for cooling in summer and heating in winter in the study area, it is necessary to test the operation effect of the heat pump system under reservoir conditions in the early stage of the construction. The variation of the temperature at the inlet and outlet of the buried pipe with time in the simulated winter and summer conditions of the four geothermal wells is shown in Figure 12. Among them, the simulation test time, the water supply temperature, and the flow rate are 24 h, 38.94°C, and 1.28 m3/h, respectively, in summer, while those are 50 h, 7.04°C, and 1.22 m3/h in winter. In the later stage of the test, the temperature and heat transfer at the inlet and outlet of the buried pipe are basically stable, which shows that the simulation test time meets the requirements. According to equation (2), the data within 12~24 h after the temperature curve reaches stable are taken to calculate the heat injection and extraction per unit of specific hole location.Figure 12
The variation of temperature with time of four geothermal holes under Case 4 and Case 5.According to the results of case 1, case 4, and case 5, the relationship curve between heat transferq with the average temperature Tf is determined, as shown in Figure 13(a). The single-hole heat transfer under different working conditions is further deduced based on the curve, which is shown in Figure 13(b). Among them, the temperature conditions of supply and return water under simulated winter and summer conditions are shown in Table 3.Figure 13
Relationship between the average temperature of supply water and heat transfer per unit hole depth of four test holes (a) and the relationship between the average temperature of supply water and single-hole heat transfer of four test holes (b). Note: “-” means taking heat from underground.Table 3
The temperature conditions of supply and return water under simulated winter and summer conditions.
ConditionsThe temperature of supply water (°C)The temperature of return water (°C)The average temperature of supply and return water (°C)Summer353032.5332830.5302527.5Winter5107.5586.5(1) Test Results. Based on the test data of the above holes, the average thermal conductivity, the heat injection, and extraction of each hole are calculated, as shown in Table 4. Therefore, the heat injection and extraction per unit of double-U heat exchanger in the shallow hole (54 m~60 m) and deep hole (100 m) are 37.9 (W/m), 48.1 (W/m), 43.7 (W/m), and 51.3 (W/m), respectively. The heat removal capacity and heat recovery capacity of the double-U heat exchanger in the deep hole have increased 15.3% and 6.7% higher, respectively, comparing the double-U heat exchanger in a shallow hole.Table 4
Test results of each test hole.
NameDepth (m)Backfill methodInitial temperature (°C)Low-power thermal conductivity (W/m·°C)High-power thermal conductivity (W/m·°C)Average power thermal conductivity (W/m·°C)The reference value of heat discharge (W/m)The reference value of calorific value (W/m)192K1154Bottom-up grouting21.91.931.901.9237.6-51.1192K1260Bottom-up grouting20.81.801.801.8038.1-45.0192K13100Bottom-up grouting20.52.042.082.0642.3-48.4192K14100Bottom-up grouting20.52.102.082.0945.0-54.2


## 3.1. Installation of Vertical Heat Exchanger

Four groups of vertical heat exchangers were arranged in this survey, all of which were installed after secondary reaming of the coring hole (Figure5). The depth of coring in holes 192K11 and 192K12 is 60 m, and double-U-shaped PE pipe heat exchangers of 54 m and 60 m are installed, respectively. Besides, the depth of coring in holes 192K13 and 192K14 is 100 m, and 100 m double-U-shaped PE pipe heat exchangers are installed. The PE pipe should be left about 1 m long on the ground considering the need of connecting the thermal response test instrument.Figure 5
Installation of vertical heat exchanger.


## 3.2. Thermal Response Test

### 3.2.1. Introduction to the Test

The main performance parameters of the equipment (Figure6) which is self-developed named FTPT11 for field thermal response test are shown in Table 1. The test equipment is mainly composed of a heat pump system, electric heater, water-replenishing tank, circulating pump, circulating pipeline, temperature and flow detection elements, etc. The flow meter is installed at the water supply pipe of the vertical buried pipe, and the temperature detection points are arranged at the water supply pipe, return pipe, and circulating pipe, which is used to automatically record the flow and temperature of circulating water in the pipe.Figure 6
Main machine (a) and refrigeration auxiliary machine (b).Table 1
Performance parameters of test equipment.
Equipment modelMaximum heating capacityMaximum cooling capacityMaximum powerCirculating water flowRangeFTPT118 kW7 kW12 kW1.8 m3/h10 mFour groups of thermal response tests have been completed in this exploration, and each group includes the initial mean ground temperature test (case 1), high-power constant heat flux test (case 2), low-power constant heat flux test (case 3), simulated constant temperature test in summer (case 4), and simulated constant temperature test in winter (case 5).The purpose of the low-power constant heat flux test and high-power constant heat flux test is to obtain the comprehensive thermal conductivity of the rock-soil body in the depth range of the vertical heat exchanger. The purpose of the simulated constant temperature test in summer and simulated constant temperature test in winter is to obtain the heat transfer power per meter length of the vertically buried tube heat exchanger under specific heat transfer conditions. Combined with the initial average ground temperature test results, the relationship between heat transfer power per unit length and medium temperature of vertical heat exchanger is further fitted, and the effective heat transfer coefficient of vertically buried tube heat exchanger is obtained.


### 3.2.2. Results and Analysis

(1) The Initial Average Temperature of the Rock-Soil Body (Case 1). Figure 7 shows the test curve of the average initial temperature of the rock-soil body in the four geothermal holes. After a certain period of system circulation, the final stable average water temperature (initial average temperature of the rock-soil body) at the outlet of four test holes is 21.9°C, 20.8°C, 20.5°C, and 20.5°C, respectively.Figure 7
Initial average temperature of the rock-soil body in 4 geothermal test holes.(2) High- and Low-Power Constant Heat Flux Test (Case 2 and Case 3). The test duration, heating power, and flow rate of each test hole are 48 h, 1.910 kW, and 1.117 m3/h, respectively, in the low-power constant heat flux test. According to the test data, the curve of the temperature at the inlet and outlet of the buried pipe with time is drawn as shown in Figure 8. Each temperature measurement curve can be divided into two obvious stages: rapid temperature rise (stage I) and stable temperature (stage II). Stage I is in the first ten hours of the test, the temperature at the inlet and outlet of the buried pipe changes rapidly, and the temperature in the borehole (backfill material) rises rapidly. Stage II shows that the temperature change in the borehole decreases with the heating time, which can be regarded as a steady-state heat conduction process.Figure 8
Results of low-power constant heat flux test in four geothermal holes.The linear source model proposed by Hellström [16] is selected to calculate the thermal conductivity. The mean temperature of the circulating medium at time τ can be expressed as follows:
(1)Tf=Tff+ql⋅Rb+14πλs⋅Eidb2ρscs16λsτ,where Tf is the mean temperature of circulating medium (°C), Tff is the undisturbed ground temperature (°C), q1 is the heat flux per length (W), ρs is the mean density of the ground surrounding the BHE (kg/m3), τ is the time (s), cs is the mean specific heat capacity of the ground surrounding the BHE (J/(kg K)), db is the borehole diameter (m), and λs is the thermal conductivity of the ground surrounding the borehole (W/(m K)).When the heating time is sufficiently long, the following equation is obtained:(2)Eidb2ρscs16λsτ≈ln16λsτdb2ρscs−γ.Linear regression is presented between the borehole wall temperature and logarithmic time. Equation (4) can be simplified to
(3)Tf=mlnt+n,where m is the slope of the linear regression between temperature and logarithmic time.The equations determine the heat exchange between the circulating medium and the surrounding ground. Thus, the ground thermal conductivityλs can be determined. According to the test data of the temperature stable section of each test hole, the curve of constant heat flow heating Tf versus lnt of each test hole is drawn as Figure 9.Figure 9
The variation curve of low-power constant heat flow heatingTf with lnt in four geothermal test holes.The test duration, heating power, and flow rate of each test hole are 48 h, 3.725 kW, and 1.046 m3/h, respectively, in the high-power constant heat flux test. According to the test data, draw the curve of the temperature at the inlet and outlet of the buried pipe with time, as shown in Figure 10. It can be seen from the curve that the shape of the curve is basically the same as that of the low-power constant heat flow test. In the first ten hours, the temperature changes at the inlet and outlet of the buried pipe are relatively fast, and the temperature rises rapidly in the borehole (backfill material). With continuous heating, the range of temperature change is reduced in the hole decreases, and the heat transfer process in the hole is similar to the steady-state heat conduction process.Figure 10
Results of high-power constant heat flux test in four geothermal holes.According to the test data of the temperature stable section of each test hole, the curve of constant heat flow heatingTf with lnt is drawn in Figure 11, and the average thermal conductivity of the rock-soil body is calculated.Figure 11
The variation curve of high-power constant heat flow heatingTf with lnt in four geothermal test holes.In conclusion, the average thermal conductivity of the rock-soil body at the four exploration holes (Table2) is calculated based on the results of the high-power and low-power constant heat flow test, so as to predict the effect of geothermal exploitation in the study area more accurately. The results show that the average thermal conductivity of K11 and K12 is close, while the value of K13 and K14 is close. This is because the holes K11 and K12 are near while holes K13 and K14 are near, and the difference of lithology between adjacent boreholes is small.Table 2
Average thermal conductivity of rock-soil body under the high- and low-power constant heat flux test.
Test type192K11192K12192K13192K14Low-power constant heat1.93 W/(m·°C)1.80 W/(m·°C)2.04 W/(m·°C)2.10 W/(m·°C)High-power constant heat1.90 W/(m·°C)1.80 W/(m·°C)2.08 W/(m·°C)2.08 W/(m·°C)(3) Simulated Constant Temperature Test in Summer and Winter (Case 4 and Case 5). Because the ground source heat pump project is used for cooling in summer and heating in winter in the study area, it is necessary to test the operation effect of the heat pump system under reservoir conditions in the early stage of the construction. The variation of the temperature at the inlet and outlet of the buried pipe with time in the simulated winter and summer conditions of the four geothermal wells is shown in Figure 12. Among them, the simulation test time, the water supply temperature, and the flow rate are 24 h, 38.94°C, and 1.28 m3/h, respectively, in summer, while those are 50 h, 7.04°C, and 1.22 m3/h in winter. In the later stage of the test, the temperature and heat transfer at the inlet and outlet of the buried pipe are basically stable, which shows that the simulation test time meets the requirements. According to equation (2), the data within 12~24 h after the temperature curve reaches stable are taken to calculate the heat injection and extraction per unit of specific hole location.Figure 12
The variation of temperature with time of four geothermal holes under Case 4 and Case 5.According to the results of case 1, case 4, and case 5, the relationship curve between heat transferq with the average temperature Tf is determined, as shown in Figure 13(a). The single-hole heat transfer under different working conditions is further deduced based on the curve, which is shown in Figure 13(b). Among them, the temperature conditions of supply and return water under simulated winter and summer conditions are shown in Table 3.Figure 13
Relationship between the average temperature of supply water and heat transfer per unit hole depth of four test holes (a) and the relationship between the average temperature of supply water and single-hole heat transfer of four test holes (b). Note: “-” means taking heat from underground.Table 3
The temperature conditions of supply and return water under simulated winter and summer conditions.
ConditionsThe temperature of supply water (°C)The temperature of return water (°C)The average temperature of supply and return water (°C)Summer353032.5332830.5302527.5Winter5107.5586.5(1) Test Results. Based on the test data of the above holes, the average thermal conductivity, the heat injection, and extraction of each hole are calculated, as shown in Table 4. Therefore, the heat injection and extraction per unit of double-U heat exchanger in the shallow hole (54 m~60 m) and deep hole (100 m) are 37.9 (W/m), 48.1 (W/m), 43.7 (W/m), and 51.3 (W/m), respectively. The heat removal capacity and heat recovery capacity of the double-U heat exchanger in the deep hole have increased 15.3% and 6.7% higher, respectively, comparing the double-U heat exchanger in a shallow hole.Table 4
Test results of each test hole.
NameDepth (m)Backfill methodInitial temperature (°C)Low-power thermal conductivity (W/m·°C)High-power thermal conductivity (W/m·°C)Average power thermal conductivity (W/m·°C)The reference value of heat discharge (W/m)The reference value of calorific value (W/m)192K1154Bottom-up grouting21.91.931.901.9237.6-51.1192K1260Bottom-up grouting20.81.801.801.8038.1-45.0192K13100Bottom-up grouting20.52.042.082.0642.3-48.4192K14100Bottom-up grouting20.52.102.082.0945.0-54.2


## 3.2.1. Introduction to the Test

The main performance parameters of the equipment (Figure6) which is self-developed named FTPT11 for field thermal response test are shown in Table 1. The test equipment is mainly composed of a heat pump system, electric heater, water-replenishing tank, circulating pump, circulating pipeline, temperature and flow detection elements, etc. The flow meter is installed at the water supply pipe of the vertical buried pipe, and the temperature detection points are arranged at the water supply pipe, return pipe, and circulating pipe, which is used to automatically record the flow and temperature of circulating water in the pipe.Figure 6
Main machine (a) and refrigeration auxiliary machine (b).Table 1
Performance parameters of test equipment.
Equipment modelMaximum heating capacityMaximum cooling capacityMaximum powerCirculating water flowRangeFTPT118 kW7 kW12 kW1.8 m3/h10 mFour groups of thermal response tests have been completed in this exploration, and each group includes the initial mean ground temperature test (case 1), high-power constant heat flux test (case 2), low-power constant heat flux test (case 3), simulated constant temperature test in summer (case 4), and simulated constant temperature test in winter (case 5).The purpose of the low-power constant heat flux test and high-power constant heat flux test is to obtain the comprehensive thermal conductivity of the rock-soil body in the depth range of the vertical heat exchanger. The purpose of the simulated constant temperature test in summer and simulated constant temperature test in winter is to obtain the heat transfer power per meter length of the vertically buried tube heat exchanger under specific heat transfer conditions. Combined with the initial average ground temperature test results, the relationship between heat transfer power per unit length and medium temperature of vertical heat exchanger is further fitted, and the effective heat transfer coefficient of vertically buried tube heat exchanger is obtained.


## 3.2.2. Results and Analysis

(1) The Initial Average Temperature of the Rock-Soil Body (Case 1). Figure 7 shows the test curve of the average initial temperature of the rock-soil body in the four geothermal holes. After a certain period of system circulation, the final stable average water temperature (initial average temperature of the rock-soil body) at the outlet of four test holes is 21.9°C, 20.8°C, 20.5°C, and 20.5°C, respectively.Figure 7
Initial average temperature of the rock-soil body in 4 geothermal test holes.(2) High- and Low-Power Constant Heat Flux Test (Case 2 and Case 3). The test duration, heating power, and flow rate of each test hole are 48 h, 1.910 kW, and 1.117 m3/h, respectively, in the low-power constant heat flux test. According to the test data, the curve of the temperature at the inlet and outlet of the buried pipe with time is drawn as shown in Figure 8. Each temperature measurement curve can be divided into two obvious stages: rapid temperature rise (stage I) and stable temperature (stage II). Stage I is in the first ten hours of the test, the temperature at the inlet and outlet of the buried pipe changes rapidly, and the temperature in the borehole (backfill material) rises rapidly. Stage II shows that the temperature change in the borehole decreases with the heating time, which can be regarded as a steady-state heat conduction process.Figure 8
Results of low-power constant heat flux test in four geothermal holes.The linear source model proposed by Hellström [16] is selected to calculate the thermal conductivity. The mean temperature of the circulating medium at time τ can be expressed as follows:
(1)Tf=Tff+ql⋅Rb+14πλs⋅Eidb2ρscs16λsτ,where Tf is the mean temperature of circulating medium (°C), Tff is the undisturbed ground temperature (°C), q1 is the heat flux per length (W), ρs is the mean density of the ground surrounding the BHE (kg/m3), τ is the time (s), cs is the mean specific heat capacity of the ground surrounding the BHE (J/(kg K)), db is the borehole diameter (m), and λs is the thermal conductivity of the ground surrounding the borehole (W/(m K)).When the heating time is sufficiently long, the following equation is obtained:(2)Eidb2ρscs16λsτ≈ln16λsτdb2ρscs−γ.Linear regression is presented between the borehole wall temperature and logarithmic time. Equation (4) can be simplified to
(3)Tf=mlnt+n,where m is the slope of the linear regression between temperature and logarithmic time.The equations determine the heat exchange between the circulating medium and the surrounding ground. Thus, the ground thermal conductivityλs can be determined. According to the test data of the temperature stable section of each test hole, the curve of constant heat flow heating Tf versus lnt of each test hole is drawn as Figure 9.Figure 9
The variation curve of low-power constant heat flow heatingTf with lnt in four geothermal test holes.The test duration, heating power, and flow rate of each test hole are 48 h, 3.725 kW, and 1.046 m3/h, respectively, in the high-power constant heat flux test. According to the test data, draw the curve of the temperature at the inlet and outlet of the buried pipe with time, as shown in Figure 10. It can be seen from the curve that the shape of the curve is basically the same as that of the low-power constant heat flow test. In the first ten hours, the temperature changes at the inlet and outlet of the buried pipe are relatively fast, and the temperature rises rapidly in the borehole (backfill material). With continuous heating, the range of temperature change is reduced in the hole decreases, and the heat transfer process in the hole is similar to the steady-state heat conduction process.Figure 10
Results of high-power constant heat flux test in four geothermal holes.According to the test data of the temperature stable section of each test hole, the curve of constant heat flow heatingTf with lnt is drawn in Figure 11, and the average thermal conductivity of the rock-soil body is calculated.Figure 11
The variation curve of high-power constant heat flow heatingTf with lnt in four geothermal test holes.In conclusion, the average thermal conductivity of the rock-soil body at the four exploration holes (Table2) is calculated based on the results of the high-power and low-power constant heat flow test, so as to predict the effect of geothermal exploitation in the study area more accurately. The results show that the average thermal conductivity of K11 and K12 is close, while the value of K13 and K14 is close. This is because the holes K11 and K12 are near while holes K13 and K14 are near, and the difference of lithology between adjacent boreholes is small.Table 2
Average thermal conductivity of rock-soil body under the high- and low-power constant heat flux test.
Test type192K11192K12192K13192K14Low-power constant heat1.93 W/(m·°C)1.80 W/(m·°C)2.04 W/(m·°C)2.10 W/(m·°C)High-power constant heat1.90 W/(m·°C)1.80 W/(m·°C)2.08 W/(m·°C)2.08 W/(m·°C)(3) Simulated Constant Temperature Test in Summer and Winter (Case 4 and Case 5). Because the ground source heat pump project is used for cooling in summer and heating in winter in the study area, it is necessary to test the operation effect of the heat pump system under reservoir conditions in the early stage of the construction. The variation of the temperature at the inlet and outlet of the buried pipe with time in the simulated winter and summer conditions of the four geothermal wells is shown in Figure 12. Among them, the simulation test time, the water supply temperature, and the flow rate are 24 h, 38.94°C, and 1.28 m3/h, respectively, in summer, while those are 50 h, 7.04°C, and 1.22 m3/h in winter. In the later stage of the test, the temperature and heat transfer at the inlet and outlet of the buried pipe are basically stable, which shows that the simulation test time meets the requirements. According to equation (2), the data within 12~24 h after the temperature curve reaches stable are taken to calculate the heat injection and extraction per unit of specific hole location.Figure 12
The variation of temperature with time of four geothermal holes under Case 4 and Case 5.According to the results of case 1, case 4, and case 5, the relationship curve between heat transferq with the average temperature Tf is determined, as shown in Figure 13(a). The single-hole heat transfer under different working conditions is further deduced based on the curve, which is shown in Figure 13(b). Among them, the temperature conditions of supply and return water under simulated winter and summer conditions are shown in Table 3.Figure 13
Relationship between the average temperature of supply water and heat transfer per unit hole depth of four test holes (a) and the relationship between the average temperature of supply water and single-hole heat transfer of four test holes (b). Note: “-” means taking heat from underground.Table 3
The temperature conditions of supply and return water under simulated winter and summer conditions.
ConditionsThe temperature of supply water (°C)The temperature of return water (°C)The average temperature of supply and return water (°C)Summer353032.5332830.5302527.5Winter5107.5586.5(1) Test Results. Based on the test data of the above holes, the average thermal conductivity, the heat injection, and extraction of each hole are calculated, as shown in Table 4. Therefore, the heat injection and extraction per unit of double-U heat exchanger in the shallow hole (54 m~60 m) and deep hole (100 m) are 37.9 (W/m), 48.1 (W/m), 43.7 (W/m), and 51.3 (W/m), respectively. The heat removal capacity and heat recovery capacity of the double-U heat exchanger in the deep hole have increased 15.3% and 6.7% higher, respectively, comparing the double-U heat exchanger in a shallow hole.Table 4
Test results of each test hole.
NameDepth (m)Backfill methodInitial temperature (°C)Low-power thermal conductivity (W/m·°C)High-power thermal conductivity (W/m·°C)Average power thermal conductivity (W/m·°C)The reference value of heat discharge (W/m)The reference value of calorific value (W/m)192K1154Bottom-up grouting21.91.931.901.9237.6-51.1192K1260Bottom-up grouting20.81.801.801.8038.1-45.0192K13100Bottom-up grouting20.52.042.082.0642.3-48.4192K14100Bottom-up grouting20.52.102.082.0945.0-54.2


## 4. Development Potential of Shallow Geothermal Energy in the Study Area

### 4.1. Limitation of Development Conditions of Shallow Geothermal Energy

The vertical ground source heat pump system can provide cooling in summer and heating in winter. The underground heat transfer power in summer and winter, respectively, refers to the heat transfer power emitted to the soil in summer and absorbed from soil in winter. Taking a building of one hundred thousand square meters as an example, assuming that the cooling load in summer is 80 W/m and the heating load in winter is 40 W/m when the ground source heat pump project adopts a double-U-type heat exchanger for100 meters, the cooling and heating load is 8000 kW and 4000 kW, respectively. The underground heat transfer powerQ1′ and Q2′ in summer and winter can be calculated by the following formula, respectively:
(4)Q1′=Q1×1+1COP1,(5)Q2′=Q2×1−1COP2,where Q1′ is the heat transfer power discharged to soil in summer (kW); Q1 is the total cooling load at the end of air conditioning in summer(kW); Q2′ is the heat transfer power absorbed from soil in winter (kW); Q2 is the total heat load at the end of air conditioning in winter(kW); COP1 is the refrigeration coefficient of water source heat pump unit under design condition; the value is 5.4; and COP2 is the heating coefficient of water source heat pump unit under design condition; the value is 4.3. According to the results, by calculating equations (4) and (5), it can be seen that the summer heat removal power should be 9481 kW, and the winter heat extraction power is 3070 kW, so as to meet the design target.


### 4.2. Layout of Geothermal Well

According to the results of various hot holes presented in Table4, it can be seen that the heat transfer capacity of a double-U heat exchanger in the shallow hole (54 m~60 m) is lower than that of the deep hole (100 m). Considering the need to reduce the floor area of a vertically buried pipe heat exchanger by improving the heat exchange capacity or the heat exchange capacity of the single hole, it is recommended to use a 100 m double-U-type heat exchanger in the later ground source heat pump project. Therefore, in order to meet the demand of cooling and heating if the cooling and heating load is calculated as 8000 kW and 4000 kW, respectively, of a building with an area of one hundred thousand square meters, the number of heat source wells to be arranged is shown in Table 5. To meet the heating and cooling requirements of the target area, about 2170 similar geothermal wells are required to operate in summer, while at least 600 geothermal wells are required to operate at the same time in winter.Table 5
Number of heat source wells required for a building with an area of one hundred thousand square meters.
Cooling load in summer (kW)Heat load in winter (kW)Heat removal power (kW)Heating extraction power (kW)Heat removal in single hole (kW)Heat extraction in single hole (kW)Number of wells required in summerNumber of wells required in winter80004000948130704.375.132170600The heating system runs for 90 days in winter and 18 hours a day. In summer, the refrigeration system operates for 120 days and 12 hours a day, and the load factor is 0.75 in winter and summer.The accumulated heat discharge in summer isQ1=9481×12×120×0.75×3.6/1000=36862 GJ.The accumulated heat discharge in winter isQ2=3070×18×90×0.75×3.6/1000=13428 GJ.


## 4.1. Limitation of Development Conditions of Shallow Geothermal Energy

The vertical ground source heat pump system can provide cooling in summer and heating in winter. The underground heat transfer power in summer and winter, respectively, refers to the heat transfer power emitted to the soil in summer and absorbed from soil in winter. Taking a building of one hundred thousand square meters as an example, assuming that the cooling load in summer is 80 W/m and the heating load in winter is 40 W/m when the ground source heat pump project adopts a double-U-type heat exchanger for100 meters, the cooling and heating load is 8000 kW and 4000 kW, respectively. The underground heat transfer powerQ1′ and Q2′ in summer and winter can be calculated by the following formula, respectively:
(4)Q1′=Q1×1+1COP1,(5)Q2′=Q2×1−1COP2,where Q1′ is the heat transfer power discharged to soil in summer (kW); Q1 is the total cooling load at the end of air conditioning in summer(kW); Q2′ is the heat transfer power absorbed from soil in winter (kW); Q2 is the total heat load at the end of air conditioning in winter(kW); COP1 is the refrigeration coefficient of water source heat pump unit under design condition; the value is 5.4; and COP2 is the heating coefficient of water source heat pump unit under design condition; the value is 4.3. According to the results, by calculating equations (4) and (5), it can be seen that the summer heat removal power should be 9481 kW, and the winter heat extraction power is 3070 kW, so as to meet the design target.


## 4.2. Layout of Geothermal Well

According to the results of various hot holes presented in Table4, it can be seen that the heat transfer capacity of a double-U heat exchanger in the shallow hole (54 m~60 m) is lower than that of the deep hole (100 m). Considering the need to reduce the floor area of a vertically buried pipe heat exchanger by improving the heat exchange capacity or the heat exchange capacity of the single hole, it is recommended to use a 100 m double-U-type heat exchanger in the later ground source heat pump project. Therefore, in order to meet the demand of cooling and heating if the cooling and heating load is calculated as 8000 kW and 4000 kW, respectively, of a building with an area of one hundred thousand square meters, the number of heat source wells to be arranged is shown in Table 5. To meet the heating and cooling requirements of the target area, about 2170 similar geothermal wells are required to operate in summer, while at least 600 geothermal wells are required to operate at the same time in winter.Table 5
Number of heat source wells required for a building with an area of one hundred thousand square meters.
Cooling load in summer (kW)Heat load in winter (kW)Heat removal power (kW)Heating extraction power (kW)Heat removal in single hole (kW)Heat extraction in single hole (kW)Number of wells required in summerNumber of wells required in winter80004000948130704.375.132170600The heating system runs for 90 days in winter and 18 hours a day. In summer, the refrigeration system operates for 120 days and 12 hours a day, and the load factor is 0.75 in winter and summer.The accumulated heat discharge in summer isQ1=9481×12×120×0.75×3.6/1000=36862 GJ.The accumulated heat discharge in winter isQ2=3070×18×90×0.75×3.6/1000=13428 GJ.


## 5. Economic and Environmental Benefits

### 5.1. Economic Evaluation

The operating economy of the ground source heat pump can be indirectly reflected by the energy consumption cost. Table6 shows the comparison of the costs of pure electric heating; boiler heating by coal fire, oil fire, and gas fire; traditional central air conditioning (dual-use cooling and heating); and ground source heat pump heating. From the perspective of energy conversion, 1 kWh electric energy can only produce 1 kWh of heat energy, and the thermal efficiency of a gas-fired boiler is only 70%~90%, which means it needs 4000~5143 kJ to produce the same heat. The traditional central air conditioning system in which COP is about 2.2~2.3 only needs 1565~1636 kJ energy to produce the same heat while 857~1091 kJ heat is needed to produce the same heat by the ground source heat pump air conditioning system in which COP is about 3.3~4.2.Table 6
Energy cost of different heating modes.
TypePriceConsumptionCost (yuan)Energy cost per square meter (yuan/m2)RemarksElectricity0.8 (yuan/kWh)4860000 (kWh)388800038.9Thermal efficiency is 1Coal600 (yuan/kg)1116172 (kg)6697036.7Thermal efficiency is 0.75Natural gas3 (yuan/m3)607500 (Nm3)182250018.2Thermal efficiency is 0.8Central air-conditioning (CAC)0.8 (yuan/kWh)2113043 (kWh)169043016.9Comprehensive energy efficiency is 2.3Ground source heat pump (GSHP)0.8 (yuan/kWh)1472727 (kWh)117818011.8Comprehensive energy efficiency is 3.3Remarks: electricity price is 0.8 yuan/kWh. Natural gas price is 3.0 yuan/m3, and its calorific value is 36.0 MJ/m3. The calorific value of coal is 20.9 MJ/kg. The heating season is 90 days, and the average load rate is 0.75.Suppose the heating area of the building is one hundred thousand square meters, and the total heat load in winter is 4000 kW. Taking the average load rate of the whole heating season as 0.75, the total heat demand in the heating season is 17496 (GJ). On this basis, the energy costs of different heating modes are calculated as shown in Table6.Comparing energy cost per square meter which is shown in Figure14, the cost of using coal-fired boiler is the lowest, followed by ground source heat pump, gas-fired boiler, central air conditioning, and electric heating. However, the coal-fired boiler will make the environmental pollution more serious by producing air pollutants and solid waste pollutants. The energy cost per square meter of the ground source heat pump is the lowest of the other three methods, and it is very good to protect the environment. At the same time, the ground source heat pump is a cold and heat source type in line with the national energy development direction, which can not only heat but also cool. Compared with other cold and heat source types, it should be vigorously promoted because of its obvious energy conservation.Figure 14
Heating cost of various energy sources.


### 5.2. Environmental Benefit

The environmental evaluation is carried out. Table7 shows the reference coefficient of converting various energy into standard coal, and the standard coal quantity (M) saved by using shallow geothermal energy to replace traditional coal for heating in winter is evaluated.Table 7
The reference coefficient of converting various energy into standard coal.
NameThe reference coefficientMean calorific valueRaw coal0.7143 kgce/kg20908 kJ/kgThermal power (equivalent value)/0.03412 kgce/MJBased on the data of saving standard coal (M) calculated by formula (6), the social and environmental benefits brought from coal saving under the same calorific value are evaluated. The calculation of emission reduction and treatment cost saving of air pollutants and solid wastes from coal combustion is shown in Table 8.
(6)M=∑Wt4.1868/7.Table 8
Pollutant emission reduction and cost-saving.
ItemCO2SO2NOxSuspended dustCoal ash residueReductionUnitt/at/at/at/at/aFormula(1)=2.386 M2=1.7%M3=0.6%M4=0.8%M5=0.1 MCostUnitTen thousand yuanTen thousand yuanTen thousand yuanTen thousand yuan/Formula0.1 yuan/kg1.1 yuan/kg2.4 yuan/kg0.8 yuan/kgFreight costsSuppose the heating area of the building is one hundred thousand square meters and the heating load is 4000 kW, the evaluation results of environmental benefits brought by the development and utilization of shallow geothermal energy for heating in winter are shown in Tables9 and 10.Table 9
Evaluation results of saving standard coal.
ProjectHeating load in winter (kW)Heat requirement for heating in winter (GJ)Saving standard coal (t)Heating400017496597Table 10
Evaluation results of emission reduction and cost-saving.
ItemCO2SO2NOxSuspended dustCoal ash residueReduction (t/a)142410.13.64.859.7Cost (ten thousand yuan/a)14.21.110.860.38Freight costsIt can be seen that the development and utilization of shallow geothermal energy resources for winter heating instead of conventional energy can reduce 59.7 tons of solid waste ash and the emission of 1442.5 tons of air pollutants such as carbon dioxide, sulfur dioxide, nitrogen oxides, and suspended dust every year when the heating area of the building is one hundred thousand square meters. Therefore, 166000 yuan can be saved on the environmental treatment (Figure15).Figure 15
Reducing emissions of pollutants.


## 5.1. Economic Evaluation

The operating economy of the ground source heat pump can be indirectly reflected by the energy consumption cost. Table6 shows the comparison of the costs of pure electric heating; boiler heating by coal fire, oil fire, and gas fire; traditional central air conditioning (dual-use cooling and heating); and ground source heat pump heating. From the perspective of energy conversion, 1 kWh electric energy can only produce 1 kWh of heat energy, and the thermal efficiency of a gas-fired boiler is only 70%~90%, which means it needs 4000~5143 kJ to produce the same heat. The traditional central air conditioning system in which COP is about 2.2~2.3 only needs 1565~1636 kJ energy to produce the same heat while 857~1091 kJ heat is needed to produce the same heat by the ground source heat pump air conditioning system in which COP is about 3.3~4.2.Table 6
Energy cost of different heating modes.
TypePriceConsumptionCost (yuan)Energy cost per square meter (yuan/m2)RemarksElectricity0.8 (yuan/kWh)4860000 (kWh)388800038.9Thermal efficiency is 1Coal600 (yuan/kg)1116172 (kg)6697036.7Thermal efficiency is 0.75Natural gas3 (yuan/m3)607500 (Nm3)182250018.2Thermal efficiency is 0.8Central air-conditioning (CAC)0.8 (yuan/kWh)2113043 (kWh)169043016.9Comprehensive energy efficiency is 2.3Ground source heat pump (GSHP)0.8 (yuan/kWh)1472727 (kWh)117818011.8Comprehensive energy efficiency is 3.3Remarks: electricity price is 0.8 yuan/kWh. Natural gas price is 3.0 yuan/m3, and its calorific value is 36.0 MJ/m3. The calorific value of coal is 20.9 MJ/kg. The heating season is 90 days, and the average load rate is 0.75.Suppose the heating area of the building is one hundred thousand square meters, and the total heat load in winter is 4000 kW. Taking the average load rate of the whole heating season as 0.75, the total heat demand in the heating season is 17496 (GJ). On this basis, the energy costs of different heating modes are calculated as shown in Table6.Comparing energy cost per square meter which is shown in Figure14, the cost of using coal-fired boiler is the lowest, followed by ground source heat pump, gas-fired boiler, central air conditioning, and electric heating. However, the coal-fired boiler will make the environmental pollution more serious by producing air pollutants and solid waste pollutants. The energy cost per square meter of the ground source heat pump is the lowest of the other three methods, and it is very good to protect the environment. At the same time, the ground source heat pump is a cold and heat source type in line with the national energy development direction, which can not only heat but also cool. Compared with other cold and heat source types, it should be vigorously promoted because of its obvious energy conservation.Figure 14
Heating cost of various energy sources.


## 5.2. Environmental Benefit

The environmental evaluation is carried out. Table7 shows the reference coefficient of converting various energy into standard coal, and the standard coal quantity (M) saved by using shallow geothermal energy to replace traditional coal for heating in winter is evaluated.Table 7
The reference coefficient of converting various energy into standard coal.
NameThe reference coefficientMean calorific valueRaw coal0.7143 kgce/kg20908 kJ/kgThermal power (equivalent value)/0.03412 kgce/MJBased on the data of saving standard coal (M) calculated by formula (6), the social and environmental benefits brought from coal saving under the same calorific value are evaluated. The calculation of emission reduction and treatment cost saving of air pollutants and solid wastes from coal combustion is shown in Table 8.
(6)M=∑Wt4.1868/7.Table 8
Pollutant emission reduction and cost-saving.
ItemCO2SO2NOxSuspended dustCoal ash residueReductionUnitt/at/at/at/at/aFormula(1)=2.386 M2=1.7%M3=0.6%M4=0.8%M5=0.1 MCostUnitTen thousand yuanTen thousand yuanTen thousand yuanTen thousand yuan/Formula0.1 yuan/kg1.1 yuan/kg2.4 yuan/kg0.8 yuan/kgFreight costsSuppose the heating area of the building is one hundred thousand square meters and the heating load is 4000 kW, the evaluation results of environmental benefits brought by the development and utilization of shallow geothermal energy for heating in winter are shown in Tables9 and 10.Table 9
Evaluation results of saving standard coal.
ProjectHeating load in winter (kW)Heat requirement for heating in winter (GJ)Saving standard coal (t)Heating400017496597Table 10
Evaluation results of emission reduction and cost-saving.
ItemCO2SO2NOxSuspended dustCoal ash residueReduction (t/a)142410.13.64.859.7Cost (ten thousand yuan/a)14.21.110.860.38Freight costsIt can be seen that the development and utilization of shallow geothermal energy resources for winter heating instead of conventional energy can reduce 59.7 tons of solid waste ash and the emission of 1442.5 tons of air pollutants such as carbon dioxide, sulfur dioxide, nitrogen oxides, and suspended dust every year when the heating area of the building is one hundred thousand square meters. Therefore, 166000 yuan can be saved on the environmental treatment (Figure15).Figure 15
Reducing emissions of pollutants.


## 6. Conclusions

Based on the analysis and calculation of the above test results, the following conclusions can be obtained, which can provide the basis for the later engineering construction:(1)
The initial average temperature of the shallow strata from 54 m to 60 m in the study area is 20.9°C~21.9°C. The average comprehensive thermal conductivity (average results of high- and low-power tests) of the strata is 1.80~1.92 W/(m. °C). The initial average temperature of the layer shallower than 100 m is about 20.5°C, and the average comprehensive thermal conductivity (average results of high- and low-power tests) of the rock-soil body in this layer is 2.06~2.09w/(m·°C)(2)
The heat removal capacity and heat recovery capacity of the double-U heat exchanger in the deep hole have increased 15.3% and 6.7% higher, respectively, comparing the double-U heat exchanger in the shallow hole(3)
The energy cost per square meter of the ground source heat pump is 11.8 yuan for a building of one hundred thousand square meters which the heat removal power is expected to be 9481 kW in summer and 3070 kW in winter. And the annual emission of carbon dioxide, sulfur dioxide, nitrogen oxides, suspended dust, and other air pollutants to the atmosphere can be reduced by 1442.5 t, and the solid waste ash and slag can be reduced by 59.7 t. The annual environmental treatment cost will be saved by 166000 yuan



---

*Source: 3271884-2021-11-09.xml*