Experimental and Theoretical Investigations of a Ground Source Heat Pump System for Water and Space Heating Applications in Kazakhstan

The ground source heat pump heating system is considered as one of the best solutions for the transition towards green heating under the continental climate conditions like Kazakhstan. In this paper, experimental and theoretical investigations were carried out to develop a ground source heat pump-based heating system under the weather conditions in Kazakhstan and to evaluate its thermodynamic performance. The water-to-water heat pump heating system, integrated with a ground source heat exchanger and used refrigerant R134a, was designed to provide hot water to meet the requirements for space heating. The predicted values of the coefficient of performance and the experimental results were found to be in good agreement within 6.2%. The thermodynamic performance of the system was also assessed using various environment-friendly refrigerants, such as R152a, R450A, R513A, R1234yf and R1234ze, as potential replacements for R134a. Although R152a is found to be a good alternative for R134a in terms of coefficient of performance and total equivalent warming impact, its flammability hinders its application. The heating system using refrigerants R450A, R513A, R1234yf and R1234ze shows 2-3% lower coefficient of performance than that of R134a. The highest exergy destruction is found to be attributed to the compressor, followed by the expansion valve, evaporator, and condenser. Considering their low flammability and low environmental impact, R450A, R513A, R1234yf and R1234ze are identified as valuable replacements for R134a.


Introduction
Heat pump systems are identified as energy efficient heating equipment used for water heating and space heating in continental climatic regions. The performance of heat pumps is improved by integrating with renewable energy sources like solar, geothermal, and solar-geothermal hybrid modes. Due to poor solar availability during winter in continental climates, the concept of integrating solar energy with heat pump systems is not feasible. Presently, the ground source heat pumps (GSHPs) are considered as the most efficient space heating system in continental climatic regions in Kazakhstan [1][2][3][4]. Low-temperature underground thermal energy source at [30][31][32][33][34][35][36][37][38][39][40] °C is available at a depth of 200 to 400 m [5,6]. To maintain thermal comfort of the room during winter, thermal energy is extracted from the ground; whereas, during summer, the excess thermal energy from the room is pumped and stored into the ground. This is especially valuable for Kazakhstan, whose majority of territory belongs to the continental climate, which is characterized by the presence of contrasted seasons from cold winters to hot summers [7]. Ambient air temperature significantly affects not only the building heat demand but also the shallow geothermal energy characteristics. For example, in severe winter conditions, since the building heat demand is high, the operation of a GSHP causes a quick depletion of the soil thermal energy around borehole heat exchanger (BHE) [8,9]. This leads to a gradual decrease in soil temperature. Such problems may arise for the central, eastern, and northern regions of Kazakhstan when severe winter occurs. In order to restore the thermal energy that has been depleted from the soil during the summer heating season, the excess solar thermal and waste heat from space cooling is pumped into the soil through ground heat exchangers [8][9][10]. Depending on the solar irradiance and the room cooling demand, this thermal energy may not be sufficient to restore the thermal energy depleted during the heating season. For the southern and western (Caspian) regions of Kazakhstan, where milder winter conditions are observed, ground thermal depletion may not be significant. In this case, when surplus solar heat and waste heat from the heat pump cooling mode in summer are pumped into the ground, not only the depleted thermal energy is restored but also excess thermal energy may be stored. Seasonal underground thermal energy storage is implemented using the borehole thermal energy storage (BTES) and aquifer thermal energy storage (ATES) methods [6,11]. In the case of BTES, underground closed-loop type BHEs are used, usually U-shaped vertical heat exchangers. In contrast to BTES system, ATES is an open-loop type system, where groundwater from shallow bore holes is circulated between two or more wells. During the heating season, the groundwater is circulated from heat source well(s) to heat sink well(s) using a pump; in summer, it is circulated in the opposite direction.
The two types of ground heat exchangers integrated with heat pump system are: (a) horizontal, and (b) vertical types. The horizontal ground heat exchangers (GHEs) are buried in trenches with a depth of 1-2 m, where the distance between the trenches is 2-3 m and the trench length is 24-48 m [12]. Due to low cost and easy installation, such systems are popular in warm regions for building heating and cooling. The horizontal type GHEs are recommended in the southern and western parts of Kazakhstan. In these locations, the higher atmospheric temperature and more solar irradiance availability positively influence the utilization of horizontal GHEs. A horizontal GHE is a closed-loop type system, where a heat transfer fluid (HTF) circulates inside to exchange heat only with the ambient soil. However, the horizontal heat exchanger is not suitable for the extreme winter conditions with low solar irradiance [13,14].
The most common type of GHE is a vertical type that can be used on a limited surface area without disturbing the landscape. The BHEs and geothermal energy piles (GEPs) [1] can harvest natural ground temperature, which is used as a heat source for heating in winter and a heat sink for cooling in summer. The GEPs serve as load-bearing deep foundations for buildings as well as GHEs for GSHPs. The most commonly used BHEs are Ushaped pipes made of high-density polyethylene (HDPE), which are installed in a well (of various diameters) backfilled with cement/bentonite grout. These vertical boreholes are buried at a depth of about 200 m, which corresponds to a shallow geothermal energy. Depending on the soil or rock type, the initial investment for drilling wells is always high. To reduce costs while maintaining the thermal potential of the soil, it is recommended to drill several wells down to 100 m. A spacing of 5-8 m is required to minimize thermal interference between several BHEs [15,16]. The BHE (GEP) is a closed-loop type GHE, where HTF circulates heat between heat pump evaporator and the ground. The BHEs are integrated with a water-to-water heat pump for efficient use of the shallow geothermal energy. Such U-shaped BHE-based GSHPs are recommended to be used in almost the entire territory of Kazakhstan, especially where severe winters are observed. The influence of the atmospheric air temperature on vertical GHEs is found to be insignificant [17][18][19][20].
Another type of GHEs with a closed-loop structure are shallow ground basket-type heat exchangers, which are separate pipe loops arranged in a spiral or in a special twisted shape, buried at a shallow depth of 1-6 m [21,22]. These type of shallow subsurface GHEs, also called energy baskets, can be immersed in a water body such as a dam or a lake away from the building [21][22][23][24][25]. Like horizontal GHE based GSHPs, these energy baskets significantly reduce the capital cost and are suitable for mild winter conditions or for buildings near lakes/dams. Such geothermal heat exchangers are not popular for the continental climatic conditions of Kazakhstan.
Closed-loop systems have been discussed above. In contrast to these systems, where only heat transfer occurs, both heat and mass transfer occurs in GHEs with an open-loop structure. An open-water loop system involves underground water that is being abstracted from shallow geothermal aquifer layers, and then returned after heat is extracted or added. Therefore, such a system consists of at least two wells where underground water from the aquifer is taken and reinjected from one well to another. One well is used for heat storage and the other for cold storage, according to the season. The wells are located at several tens of meters from each other, excluding the mutual influence of warm and cold "bells", and are connected on the ground to a water-to-water type heat pump by a pipeline [26]. In most systems, the winter temperature of the reinjected water is 6-9 °C, while the summer temperature is 15-25 °C [26]. One of the important advantages of such a system is the large volumes of water abstracted, which results in a large heat exchange capacity and can be used in large scale applications [27][28][29][30].
Among these types of GSHPs, heat pumps with BHEs and wells from aquifers are the most mature technologies for utilizing shallow geothermal energy. The BHEs do not require special preliminary survey of geological or hydrogeological structure on the presence of aquifers and are applicable to all climates. The seasonal variation of the ground temperature due to seasonal change in the ambient air temperature only occurs at the depth of 15 m below ground surface for cold regions such as Kazakhstan, Russia, Northern China, Mongolia, Northern USA and Canada. Therefore, the use of such vertical BHE (below 50-150 m depth) as part of the GSHP will be efficient and widespread. In this paper, an experimental and theoretical study of a GSHP is carried out using U-shape BHEs under the weather conditions of the city of Almaty, in the Republic of Kazakhstan.
The vapor compression water-to-water heat pump is used to enhance the low-grade heat extracted through GHEs from shallow geothermal resources. The vapor-compression heat pump cycle can be described by the reverse Carnot cycle, an ideal circular process consisting of two adiabatic and two isothermal processes [31,32]. The main working fluid is a refrigerant, where the heat exchange process takes place with phase change. In an indirect expansion system, a water-to-water or brine-to-water heat pump is a main element [33,34]. GSHPs with aquifers as heat sources use a water-to-water heat pump. BHEbased GSHPs use both water-to-water and brine-to-water heat pumps, depending on the source and sink temperature ranges. Recent works in this area are reported in [35][36][37].
Many research initiatives have been devoted to heat pump systems with vertical and horizontal ground heat exchangers under various weather climates [2,3,7,8,19]. However, to the best of our knowledge, thermodynamic performances for the specific weather conditions of Kazakhstan have not yet been reported. Moreover, halogen-based refrigerants like R134a, R407C and R410A have been widely used in heat pump systems integrating with ground heat exchangers. These refrigerants have a high GWP, and many researchers are looking for low GWP alternatives [35][36][37]. Hence, the possibility of using ground source heat pumps using environment friendly refrigerant mixtures for hot water and space heating applications in Kazakhstan is considered in this research. A detailed energy and exergy analysis is carried out to assess both energy conversion and energy losses in the system.

Apparatus
The experimental setup installed in Almaty is depicted in Figure 1. Two vertical Ushaped BHEs are used to extract the low-grade ground thermal energy using a single Upipe and double U-pipe of 50 m depth. A water-to-water heat pump with a rated heating capacity of 5 kW is developed and integrated with BHEs for water heating and space heating. A hot water storage tank of 300 L is used to store the thermal energy. To install the two BHEs, two wells of depth 50 m have been drilled. High-density polyethylene (HDPE) piping material of 32 mm diameter is used for the U-shaped heat exchangers. After immersing the heat exchangers into the well, the voids are filled with grout, which is a solution of water (70%), cement (24%) and bentonite (6%). After installing the GHEs, pressure drops and air in the hydraulic loop are checked. A thermal response test (TRT) is carried out to determine the soil thermal conductivity and the BHE thermal resistance. A water-to-water heat pump operates for two modes (heating and cooling) using a reversing valve. The performance of the water-to-water heat pump, integrated with ground heat exchanger operating in heating mode, is tested for water heating and space heating. Figure 2 shows the schematic of the water-to-water heat pump integrated with GHE. The system has three loops: GHE loop, heat pump refrigerant loop and hot water storage loop. The GHE circuit consists of a U-shaped BHE, a circulation pump, an expansion tank, and an evaporator of the heat pump. The heat pump refrigerant loop consists of a compressor, a condenser, an expansion device, and an evaporator. The detailed specifications of the components of the water-to-water heat pump are given in Tables 1-3. Accessories such as filter-drier, liquid receiver and sight glass are installed in the liquid line. The reversing valve is used to divert the refrigerant flow according to the requirement. The hot water storage tank loop consists of a condenser, a circulation pump, a 300-L storage tank and an expansion tank. The temperature and pressure of refrigerants in the water-to-water heat pump loop are measured with thermocouples and pressure gauges (see Figure 2).    The present work reports the experimental data of a GSHP heating system and the thermodynamic analysis based on the 1st and 2nd laws of thermodynamics [31,32]. Figure  3 shows the vapor-compression heat pump cycle on the pressure-enthalpy diagram. The HTF leaving the BHE enters the evaporator at point 1, and the heat is transferred from the HTF to the refrigerant. The heat transfer of the HTF in both the BHE and the evaporator is single-phase convective heat transfer. The HTF is heated in the condenser from point 7 to 8 and then transfers heat to water in the thermal energy storage tank. The heat transfer of the HTF in both the condenser and the tank is also single-phase convective heat transfer. The processes 3-4, 4-5, 5-6, and 6-3 represent the compression, condensation, expansion, and evaporation processes, respectively. The processes 5′-5 and 3′-3 represent the sub-cooling of the refrigerant in the condenser and the superheating in the evaporator, respectively. Process 3-4s corresponds to the ideal compression process [31,32], while process 3-4 represents the actual compression process in the compressor. States 3′, 4′ and 5′ are the saturation states. The thermodynamic analysis was carried out using the Engineering Equation Solver Pro V10 (EES) software [38].

Energy, Exergy and Environment Impact Analysis
This section provides a detailed presentation of the set of equations for evaluating the energy and exergy efficiency of the system, its components, and the system and vapor compression cycle . The description of thermodynamic processes with the corresponding points has been presented in the previous section and in Figure 3. The equations presented in this section are derived based on the assumptions that the compression process is isentropic, the evaporation and condensation processes are isobaric, and the expansion process is isenthalpic [2][3][4]. Energy and exergy analyses are based on the 1st and 2nd laws of thermodynamics [31,32].

Evaporator
The heat source for GSHP is the ground, which is extracted using BHE. The temperature of the HTF at point 1 ( ) was assumed to be equal to the source temperature ( ). The approach ( ) and superheating ( ) temperatures in the evaporator area were taken from experimental observations. It follows that the temperature in the compressor suction line ( ) is given by Equation (1) [31,32]: (1) In this case, the HTF temperature at point 2 is determined by Equation (2)  ( The amount of heat received by the evaporator from the low-grade ground heat source is:

Compressor
According to the definition of isentropic efficiency of the compressor, the refrigerant enthalpy after compression (state 4) can be calculated as follows: In this case, the isentropic efficiency is determined by Equation (5) The power and electricity consumption of the compressor can be determined as: , , , where is combination of mechanical, transmission and electrical motor efficiencies [31,32], and is the ratio of ideal isentropic work to real electricity consumption. In this paper, coefficient (in %) is taken as a function of and according to the compressor's manufacturer (Copeland ZR28K3E-PFJ) Equation (9):

Condenser
The heat sink for GSHP is the water storage tank. The temperature of the HTF at point 7 ( ) was taken equal to the sink temperature ( ). The approach ( ) and subcooling temperatures in the condenser area were given by experimental observations. Thus, condenser outlet temperature can be determined by Equation (10) [31,32]: (10) In this case, the temperature of the HTF at point 8 is determined by Equation (11) [2,3]: , . ( The amount of useful heat from the condenser is determined by Equation (12):

Expansion Device
Since the expansion process is assumed to take place at a constant enthalpy, the enthalpy at the exit of the expansion valve is equal to that at the inlet, i.e., ℎ ℎ .

System Parameters
The refrigerant mass flow rate is obtained from the condenser heating capacity by Equation (13): The overall electrical power input for running the system is the sum of the power consumptions of the compressor and the two circulating pumps: The volumetric refrigerating capacity at the outlet of the evaporator is calculated as: The of the vapor compression heat pump cycle (process 3-4-5-6) and the of the GSHP system are given by Equations (16) and (17), respectively: .
The of the reverse Carnot cycle is expressed in terms of the source and sink temperatures as: .

Exergy Destruction
Exergy analysis provides a method for estimating the maximum work extracted from a substance relative to its environment as a reference state (the "dead state") [2,3]. Exergy analysis is widely used as a method for the thermodynamic optimization of thermal energy systems, which helps determine the inefficient components of GSHP [31,32].
Exergy destruction of vapor compression heat pump cycle components is determined by Equation (19), which is derived from the steady-state form of the exergy rate balance: .
The first and second terms of the right-hand-side of Equation (19) are the flow exergy at the inlet and outlet of a specific system component; the third term is related to the work done by the component, and the fourth term corresponds to the heat exchange between environment and the component.
Specific flow exergy and flow exergy at the th point are determined based on enthalpy and entropy of this point and the dead state: ℎ ℎ .
For each component of the heat pump, the exergy destruction is given by Equation (22) to Equation (25) for the evaporator, compressor, condenser, and expansion valve, respectively: , , , , .
The overall cycle exergy destruction is determined by Equation (26), as the sum of th components exergy destruction: The relative exergy destruction of the th component is determined by Equation (27):

Exergy Efficiency
The evaporator and condenser are brazed plate heat exchangers, in which heat is transferred with counter-current flows (see Table 3). The exergy efficiency of the evaporator and condenser can be determined by Equations (28) and (29), respectively: . .
The vapor compression process of the refrigerant in the compressor is assumed to be adiabatic. Therefore, its exergy efficiency can be determined by Equation (30): The expansion process of the refrigerant is assumed to be isenthalpic. The exergy efficiency of the expansion valve can be determined by the ratio of the exergy output to the exergy input [40]: .
The exergy efficiencies of the heat pump cycle and GSHP system can be calculated by Equations (32) and (33) [2,3]: .
The second-law efficiency of the cycle is defined as the ratio of the cycle to the of the reverse Carnot cycle [31,32]: .

Total Equivalent Warming Impact
The total equivalent warming impact ( ) includes direct emission due to the system refrigerant leakage losses and indirect emission due to the heat pump electricity consumption.
is assessed by Equation (35) [41]: where is the Global Warming Potential of the refrigerant (in CO2 eq.), is the charge of the refrigerant (in kg), is the leakage rate per year (in %/year , is system operating life (in years); is the recovery/recycling factor (in %), is the energy consumption per year (in kWh/year), is the indirect emission factor (in kg/kWh).

Results and Discussions
This section presents experimental observations, validation of theoretical models, and thermodynamic performance comparison of water-to-water heat pump using environment friendly refrigerants. In addition, the total equivalent warming impact of the system working with R134a, and its possible alternatives is discussed.

Experimental Results
The measurements were carried out using the measuring devices shown in Figure 2. As seen from Figure 2, temperature and pressure sensors are installed at the nodal points of the GSHP. In addition, in the hydraulic loops between BHE and the evaporator, as well as between the condenser and the water storage tank, liquid flow meters are installed. The observations are depicted in Figures 4-8 for the water heating in the storage tank.      As shown in Figure 6, the source temperature changes from 12.45 °C to 13.83 °C. The ground type of the local terrain is mainly composed of gravel and till/loam. For the specified type of ground, the operation of the constructed GSHP gives the specified regression dependence between two temperatures, which is then used for validation (see Table 4). A similar relationship for measuring the compressor electricity consumption , is shown in Figure 7.  Furthermore, to compare the measured and calculated results, the following parameters such as superheating , sub-cooling , evaporator approach temperature , and condenser approach temperature are used and encountered in Equations (1) and (10). Based on the superheating and sub-cooling results, the operation of the vapor compression heat pump system, the amount of refrigerant mass charge in the system, and the thermostatic expansion valve, adjustments are ascertained. Figure 8 Table 4).

Thermodynamic Model Validation
This subsection contains the comparison of experimental and theoretical results on temperatures in the heat pump cycle, evaporator, and condenser heat fluxes, of the cycle and system. For the thermodynamic calculations based on Equations (1)-(18), the parameters indicated in Table 4 are taken. The theoretically predicted thermodynamic performance results are compared with experimental results for validation. Figure 9 illustrates the temperature at typical nodal points as the function of the sink temperature . The experimental and calculated values of the refrigerant temperatures at the compressor suction line , as well as at the condenser outlet , are in very good agreement. The temperature difference at the evaporator inlet is 3.16-4.01 °C with an average value 3.66 °C. The difference at the compressor discharge line is 0.06-6.12 °C. This difference between the experiment and the model is because the model does not consider the presence of a four-way reverse valve.     (17) vary in the range of 2.28-2.98. The minimum, maximum and average differences between the experimental and theoretical values are 4.9%, 7.5% and 6.2%, respectively. Therefore, the thermodynamic model developed in the present work is validated and can be used to further analyze the thermal performance of the GSHP system. Figure 11. Comparisons of predicted and measured of the cycle and system with the temperature of the heat sink.

Thermodynamic and Environmental Calculations
Using the theoretical model, the thermodynamic performance of the heat pump is simulated for the environment friendly refrigerant options. The principle behind selecting alternative refrigerants is affordability, similar thermodynamic and thermo-physical properties, low values, and non-or mild flammability. Table 5 shows the thermodynamic properties of selected refrigerants in this work. The thermodynamic analysis is performed using the model of Equations (1)-(35) and the operating conditions given in Table 4. However, instead of the given input power consumption, a value of = 5 kW is fixed as an input condition. Figure 12 shows the relationship between required of the system and the specific latent heat of vaporization of refrigerants. The highest value of the specific latent heat of vaporization corresponds to refrigerant R152a with the lowest value of while the lowest value of the specific latent heat refers to refrigerant R1234yf with the highest value of . As seen from Figure 12, the refrigerant mass flow rate required with the increase in the specific latent heat of vaporization of refrigerants.  Figure 13 shows the comparison of the required refrigerant mass flow rate and volumetric refrigeration capacity ( ) for a given heat capacity for refrigerant R134a and environment-friendly alternative refrigerants. For the same heating capacity as refrigerant R134a, the low refrigerant mass flow rate ( ) with high results in a smaller size compressor. In this regard, refrigerant R152a is seen to be the best option. For R152a, the refrigerant mass flow rate ( ) required is estimated to be 17.5 g/s and 1949 kJ/m 3 . The mass flow rates required for R450A and R1234ze(E) are close to that of R134a, with the values of 31.4 g/s and 32.0 g/s, respectively. However, refrigerant R450A and R1234ze(E) have lower values of than that of R134a, which requires a large volume compressor. In addition, R513A and R1234yf have closer values of to that of R134a at higher mass flow rates.  Figures 14 and 15 show the comparison of the exergy destructions and relative exergy destructions of the system components and cycle for different refrigerants, respectively. Among the selected refrigerants, the minimum system exergy destruction is observed for R152a with a value of 965 W and the maximum exergy destruction of 1119 W is observed for R1234yf. Among components, the maximum exergy destruction is observed in the compressor, followed by the expansion valve, evaporator, and condenser.  As is seen from Figure 15, the expansion valve can apparently affect the growth of system exergy destruction for these refrigerants. The share of the expansion device RE increases from 16% to 23%, while the relative exergy destructions of other components vary slightly by 1-3%. Referring to Table 4, since and HTF flow rate are fixed for the evaporator, for these refrigerants, there is no significant difference in . Although and for the condenser are fixed, is different. This is due to the different discharge temperature depending on the thermodynamic properties of the refrigerant. Figure 16 shows the effect of discharge temperature on the condenser exergy efficiency . The highest value of corresponds to refrigerant R1234yf with the lowest value of whereas the lowest value of refers to refrigerant R152a with the highest value of . It is found that a linear relationship between and . Thus, this dependence shows that with an increase in the discharge temperature, the condenser efficiency decreases in this operation mode. For these refrigerants, the condenser exergy efficiency varies between 94.9% and 99.3%. The exergy efficiencies of other components are given in Table 6. The exergy efficiencies of the evaporator, compressor and expansion valve are in the range of 50.4-53.6%, 58.3-58.8% and 79.9-85.7%, respectively. The second law , cycle , and overall system efficiencies are calculated using Equations (32)- (34). Figure 17 shows the evaluating results of these efficiencies. These efficiencies for all the refrigerants appear a similar trend in the ranges of 32.3-35.3%, 29.4-32.1%, and 27.8-30.2%, respectively. As is seen from Table 6, the maximal values of the efficiencies correspond to refrigerant R152a, and the minimal values to refrigerant R1234yf. Among the selected refrigerants, only R152a is superior to R134a while R1234ze(E) has the closest values to those of R134a. The values of and vary similarly.  Figure 17. Second law, cycle, and system efficiencies for different refrigerants. Figure 18 shows the total equivalent warming impact ( ) of a water-to-water heat pump using R134a and environment-friendly alternative refrigerants. The includes both direct and indirect emissions depending on the of the refrigerant. R134a appears the maximum direct emission due to its high than other refrigerants, followed by R513A, R450A, and R152a having the direct impacts of 1.14, 0.84, and 0.12 tons of CO2 equivalent, respectively. The lowest value of direct emissions is observed for HFO refrigerants like R1234yf and R1234ze(E) with a value of 0.01. The lowest value of indirect emissions is observed for R152a with 31.64-ton eq. CO2 while the highest value of 34.32ton eq. CO2 is associated with R1234yf due to its high compressor power consumption. The carbon-dioxide emission factor ( ) is assumed to be 0.8 kg of CO2 per kWh for coalbased electric power generation system [41].  for different refrigerants. R134a is used as the base refrigerant and others refrigerants are considered as the low alternatives. HFCs such as R134a will be phased out soon according to the Paris protocol 2016 (Kigali agreement). The second quadrant is the best choice for both cases in Figure 19 (a) and (b). In this respect, a normalized lower than one, and a greater than one is preferable. Similarly, a normalized greater than one, and a lower than one is preferable.
(a) (b) Figure 19. Normalized cycle parameter dependences for different refrigerants: (a) Coefficient of performance vs. total equivalent warming impact; (b) Cycle efficiency vs. exergy destruction.
As shown in Figure 19, R152a is a good choice with higher and than R134a; moreover, its values of and are the lowest. The next refrigerant with lower is R1234ze(E), while the remaining parameters are close to those of R134a. The lowest efficiency parameters are observed for R1234yf with a lower than that of R134a, but with the highest exergy destruction. Refrigerant R450A and R513A display similar properties, with efficiencies lower than that of R134a despite their low . Refrigerant R152a is seen to be a good alternative to replace R134a. However, due to its flammability, the mass charge of R152a is limited [40]. The HFO refrigerants such as R1234yf and R1234ze(E) are considered to be the key substitutes for the HFC refrigerants.
Even though these HFOs have low global warming potential, the thermo-physical properties are found to be lower than those of R134a, which results in poor thermodynamic performance. Moreover, fast decomposition of HFO refrigerants in the atmosphere results in formation of trifluoroacetate (TFA) [42]. For example, the level of TFA for R1234yf is almost five times higher than that for R134a [42]. The release of TFA in the atmosphere leads to the pollution of aquatic systems [2,3]. Therefore, HFC/HFO blends such as R450A and R513A are being considered as alternative refrigerants. Refrigerant R450A and R513A are non-flammable with a 60% lower than that of R134a. For these reasons, R450A and R513A can be recommended as good alternatives.

Conclusions
A prototype of a ground source heat pump heating system has been developed and installed in the Koksai Mosque, in the Almaty region of Kazakhstan. Experiments have been performed for refrigerant R134a. The energy and exergy efficiencies of the ground source heat pump heating system operating under the weather conditions of Kazakhstan have been analyzed. The thermodynamic model was validated by comparison with experimental data. Based on the verified model, calculation of energy and exergy efficiencies for environmentally friendly refrigerants were carried out as alternatives to R134a. The environmental impact of the heating system using different refrigerants has been estimated using the parameter. The following major conclusions can be drawn from this study:  The experimental values of vary in the range of 2.24-2.97, while vary in the range of 2.11-2.76.  The predicted and experimental results are found to be in good agreement within 6.2%.  The thermodynamic model describes the energy and exergy efficiencies of the system with sufficient accuracy.  As input data for model calculations, we have used measurements from the borehole heat exchanger, that takes into account the local ground conditions.  The maximal exergy destruction is attributed to the compressor, followed by the expansion valve, evaporator, and condenser.  The exergy efficiency of the condenser depends on the compressor discharge temperature for different refrigerants, and it is maximal in the case of R1234yf.  The R152a is found to be a good alternative to R134a in terms of , , and with estimated values of 3.09, 2.91, and 31.76-ton eq. CO2, respectively.  Refrigerant R450A and R513A are found to be safer alternatives to R134a with slightly (of order 3%) lower , , and .  Refrigerant R1234yf and R1234ze are identified as interim replacements to R134a due to the formation of trifluoroacetic acid.