بهره وری از حلقه پمپ حرارتی زمین گرمایی بسته: تجزیه و تحلیل حساسیت

Geothermal heat pumps are becoming more and more popular as the price of fossil fuels is increasing and a strong reduction of anthropogenic CO2 emissions is needed. The energy performances of these plants are closely related to the thermal and hydrogeological properties of the soil, but a proper design and installation also plays a crucial role. A set of flow and heat transport simulations has been run to evaluate the impact of different parameters on the operation of a GSHP. It is demonstrated that the BHE length is the most influential factor, that the heat carrier fluid also plays a fundamental role, and that further improvements can be obtained by using pipe spacers and highly conductive grouts. On the other hand, if the physical properties of the soil are not surveyed properly, they represent a strong factor of uncertainty when modelling the operation of these plants. The thermal conductivity of the soil has a prevailing importance and should be determined with in-situ tests (TRT), rather than assigning values from literature. When groundwater flow is present, the advection should also be considered, due to its positive effect on the performances of BHEs; by contrast, as little is currently known about thermal dispersion, relying on this transport mechanism can lead to an excessively optimistic design.

Ground Source Heat Pumps (GSHP) are space heating and cooling plants which exploit the soil as a thermal source or sink, through the circulation of a heat carrier fluid in a closed pipe loop. Different pipe arrangements are available, among which the most common is the Borehole Heat Exchanger, a vertical pipe loop reaching depths of 50–200 m (Fig. 1). Below a depth of a few meters from the ground surface, the seasonal variation of the air temperature disappears due to the large thermal inertia of the soil. Therefore, if compared to the air, the soil is a warmer source for heating during winter and a cooler sink for cooling during summer, and higher system efficiencies can therefore be achieved compared to Air Source Heat Pumps. Full-size image (43 K) Fig. 1. Scheme of a Ground Source Heat Pump (GSHP): the Borehole Heat Exchanger (BHE) exchanges heat between the surrounding soil and the heat pump. A thermal storage tank reduces the frequency of start-up and stop of the heat pump. Radiant panels and fan coils are the most diffused heating terminals for GSHPs. If present, groundwater flow enhances the heat transport around the BHE, permitting to achieve better energy performances. Figure options GSHPs are rapidly spreading in Europe, China and USA, and have a great potential for energy, cost and CO2 emission saving [1]. About 100,000 low-enthalpy geothermal plants are installed every year in Europe, mainly for new dwellings in Sweden, Germany and France [2] and [3]. According to Saner et al. [4], the use of GSHP in place of methane furnaces allows the CO2 emissions to be reduced by up to 84%, depending on the sources used for the production of electricity. From the economic point of view, the geothermal heat pumps lead to a considerable reduction of the maintenance costs and, although their installation is more expensive than the other heating and cooling plants, the payback periods proved to be reasonable, i.e. less than 10 years [5], [6] and [7]. Since the thermal exploitation of the soil induces a gradual temperature drift, an accurate heat transport modelling of soil and aquifer systems is essential for a correct design of GSHPs. Indeed, the efficiency of the heat pump is strongly influenced by the temperature of the heat carrier fluid, which in turns depends on the temperature of the surrounding soil. To estimate the thermal impact of BHEs and the working temperatures of the heat carrier fluid, different methods have been developed, which can be divided into analytical, semi-analytical and numerical. The Kelvin infinite line source [8] and the infinite cylindrical source [9] are the simplest analytical methods for estimating the thermal disturbance induced by a BHE, since they rely on the assumption of a purely conductive and radial heat transport. Their main limitation is that of not accounting for the vertical thermal gradient and fluxes [10] and for the heterogeneity of the heat exchange over the length. Moreover, the advective and dispersive heat transport occurring in aquifer systems is also neglected. Nevertheless, these analytical solutions are still widely used for the interpretation of Thermal Response Tests [11], since they last for a short time (48÷72 h) and therefore the vertical heat transport can be neglected. The subsurface flow and the seasonal changes of groundwater levels can significantly alter the results of a TRT, as pointed out by Bozdaǧ et al. [12]. To overcome this problem, Wagner et al. [13] recently developed a method for the interpretation of TRTs in the presence of strong groundwater flow. The semi-analytical method proposed by Eskilson [14] takes into account the finite length of the exchanger and different BHE field layouts, but the advection and the dispersion are neglected. This method is applied by two of the most popular BHE design software programmes, Earth Energy Design [15] and GLHEPRO [16]. Analytical models which take into account the beneficial effects of groundwater flow [17], of the finite length of the BHE [18], and both them together [19] have been developed in the last few years, and they could be used in the future for the dimensioning of BHE fields. Recently, numerical modelling has often been applied to the design of BHE fields. The finite-difference modelling software MODFLOW can be used coupled with the solute transport package MT3D (or MT3DMS) and by applying the analogy between heat and solute transport [20] and [21], or with the specific heat transport package SEAWAT [22]. On the other hand, the finite element software FEFLOW includes a special package for the simulation of BHEs [23] and [24] which is particularly suitable for non conventional BHE field layouts and for taking into account the thermal advection and dispersion in aquifer systems. The heat transport simulation of GSHPs permits the assessment of their performances, which are influenced by the properties of the exchanger and the thermo-hydrogeological parameters of the soil. According to Chiasson et al. [25], groundwater flow significantly enhances the performances of BHEs, and the Peclet number is a good indicator for whether advective transport needs to be taken into account or neglected. Wang et al. [26] have developed a method to estimate the velocity of groundwater movement measuring the temperature profiles in a BHE. Lee [27] has investigated the effect of vertical heterogeneities of the soil thermal conductivity, concluding that the adoption of depth-averaged thermal parameters is appropriate. Chung and Choi [28] have found that an increase of the fluid flow rate reduces the heat transfer rate per unit length. Delaleux et al. [29] have studied the increase of the thermal conductivity of grouts with the addition of graphite flakes, concluding that a noticeable heat transfer improvement is achieved by BHEs. Jun et al. [30] have evaluated the influence of running time, pipe spacing, grout conductivity, borehole depth, fluid flow rate, inlet fluid temperature and soil type on the heat transfer length and on the thermal resistance of borehole and soil. Michopoulos and Kyriakis [31] have found a non-linear relation between the BHE length and the heat pump consumption, which can be used for optimization processes in the dimensioning of large plants. The aforementioned studies deal with single or few parameters, but a thorough comparative analysis of all these factors together is still missing, and constitutes the objective of this work. The functioning of a single BHE was simulated for 30 years, using a benchmark cyclic thermal load and changing the operational parameters of the scenario. The resulting fluid temperatures at the end of the BHE were processed and used to estimate the COP of the heat pump and its annual energy consumption under different conditions. On the basis of the results it is possible to draw some practical conclusions on the margins of improvement of BHEs and on the proper choice of soil parameters for the simulations.

In this work, the most important parameters which influence the performances of Ground Source Heat Pumps have been thoroughly analysed, running long-term simulations and estimating the energy consumption of the heat pump for each setting. Most of these factors have been already analysed in other works, but none of them considered all the parameters together, using the same modelling framework and considering the effect on the lifetime of a GSHP. The analysis of the BHE design parameters (length, pipe spacing, fluid, grout) permits to understand which are the margins of improvement, while the physical parameters of the soil (thermal conductivity and dispersivity, groundwater flow) have been studied in order to understand their effect on the uncertainty in the project phase. The results of the simulations prove that the length of the BHE is the most important parameter in the design of a GSHP. Indeed, increasing the borehole depth results in a reduction of the thermal disturbance in the subsoil and a higher efficiency of the heat pump is achieved, but also a larger investment is needed for the installation. An optimum length should be found, which minimizes the total cost over the plant lifetime, considering also the trend of increase of the unit cost of electricity. While the drilled depth has an appraisable impact on the initial investment, there are also other important factors to be considered for the optimization of BHEs, like the pipe arrangement, the grout and the heat carrier fluid. A large pipe spacing and a highly conductive grout, reducing the heat losses in the heat exchange with the soil, achieves an appraisable reduction of the energy costs for the heat pump with a negligible expense, compared to the borehole drilling. For the circulation pump, a trade-off can be found for the choice of the correct flow rate for the heat carrier fluid, allowing the minimization of both the energy losses due to the thermal resistance and the friction losses due to the circulation of the fluid. The antifreeze and its concentration heavily influence the energy performance of GSHPs, in particular the borehole resistance and the power consumed by the auxiliary plants. The saline solutions, with a smaller viscosity compared to ethanol and glycols, permit to reduce all these energy losses, although special components are needed to avoid corrosion problems. Optimizing the design and the installation of BHEs is useless without a thorough characterization of the subsoil, which has a large influence on the performances of these systems. When no groundwater flow occurs, the thermal conductivity is the most important parameter for the dimensioning of BHEs. The technical literature provides wide ranges of the thermal conductivity for each lithology, which can vary due to porosity, saturation and other factors; in-situ Thermal Response Tests are therefore strongly advised for large plants to avoid under or over dimensioning. The advection enhances the performances of GSHP, and the groundwater flow should be taken into account using conservative values of hydraulic conductivity and gradient, unless they are known by field tests. On the other hand, it is risky to consider also the beneficial effect of heat dispersion, because the thermal dispersivity is still scarcely known in real-scale BHEs. In situ tests to estimate these parameters would be highly desirable to simulate the behaviour of BHE fields with a better precision.