Study of Heat Transfer Characteristics and Economic Analysis of a Closed Deep Coaxial Geothermal Heat Exchanger Retrofitted from an Abandoned Oil Well

: It is economical to transform abandoned oil/geothermal wells into closed deep geothermal heat exchangers with coaxial tubes. A numerical model of a coaxial geothermal heat exchanger (CGHE) with varying borehole diameters is established according to an abandoned well in Northern China. The finite difference method is adopted to solve the temperature distribution, and the accuracy of the model is validated with experimental data. Based on the existing structure of the abandoned well with different depths, the feasibility of its conversion into a deep CGHE is discussed, and this study uses the orthogonal experimental method to analyze the influence of four main factors and their significance level on the average heat extraction rate, with the heat extraction rate up to 422.18 kW in the optimal combination. This study also integrates with actual project considerations and conducts an economic analysis to determine the most appropriate circulation fluid flow rate. The results highlight the key factors on the heat transfer performance of the CGHE, with the inlet water temperature to the CGHE being the most significant, followed by the configuration of the CGHE retrofitted from abandoned. From the economic perspective, given that the CGHE in this study is retrofitted from the abandoned oil Wells, the drilling cost can be reduced by up to CNY 1800 thousand, and the flow rate design of 35 m 3 /h is the optimal choice, ensuring a cost-effective system operation while meeting the operational requirements of the deep CGHE.


Introduction
It is estimated that coal remains the primary energy source for 92% of space heating systems in China [1], precipitating a substantial contribution to greenhouse gas emissions.Consequently, it is necessary to utilize renewable energy for building heating to reduce coal consumption and alleviate air pollution.Geothermal energy has shown significant potential in space heating because of its environmental friendliness and reduced energy consumption [2,3].The installed capacity of geothermal heat pumps for space heating and cooling exceeds 26,450 MWt in China by the end of 2020 [4].
Deep geothermal space heating systems have recently aroused growing interest in academic and engineering fields because of their enhanced thermal capacities and decreased installation space requirements in contrast to shallow geothermal systems [5,6].The most crucial component of a deep geothermal heating system is the deep coaxial geothermal heat exchanger (CGHE), also known as coaxial borehole heat exchanger or wellbore heat exchanger.Nevertheless, drilling a borehole incurs the most significant financial investment of any geothermal system, accounting for up to 70% of the total investment costs on average, and therefore makes these systems economically unprofitable [7].This prompts the desire to explore alternative solutions for efficient and economical acquisition of geothermal energy.
In recent years, retrofitting an abandoned oil well into a CGHE has emerged as a promising method to enhance geothermal energy applications.For instance, Morita et al. [8,9] reported the pioneering CGHE experiment in an abandoned well on the island of Hawaii.Following this, Kohl et al. [10,11] reported the performance of a CGHE at Weissbad and Weggis, Switzerland.Abandoned oil and gas wells (AOGWs) with high bottom-hole temperatures contain abundant geothermal energy.This energy can be harnessed by retrofitting these wells into novel geothermal systems for various uses without the need for costly drilling [12].Furthermore, Cheng et al. [13,14] successfully employed the organic Rankine cycle to exploit the geothermal energy of abandoned oil wells.They also presented a novel method to enhance geothermal utilization efficiency through the development of thermal reservoirs.
Previous studies have demonstrated the feasibility of harvesting geothermal energy from existing abandoned oil and gas wells.A wide range of wells across different countries have been covered, such as horizontal wells [15], enhanced geothermal systems (EGS) and low-temperature deep borehole heat exchangers (DBHE) in California [16], U-tube heat exchanger located in Southern Iran [17], and double-pipe geothermal and Jachowka K-2 [18] heat exchanger in the Persian Gulf [19].Some studies also proposed mathematical models to analyze key factors affecting heat outputs from CGHEs.A software package (ANSYS 2019) was developed based on the finite difference method for thermal analysis of a CGHE [20].An unsteady-state heat transfer model [21] was also established for a deep CGHE and validated with experimental data.A transient 2D finite volume method (FVM) was applied to solve the thermal interference under seasonal effects and dynamic heat flux within a vertical coaxial borehole heat exchangers field [22].In consideration of well kick occurrences during the drilling process, the models were solved using the fully implicit finite difference method coupled with the clustering method [23].
For the performance optimization of the CGHE, Bu et al. [24] developed mathematical models to investigate the heat exchange between fluid and rocks.Based on the models, they also examined the effects of critical parameters and determined their optimal values.Yekoladio et al. [25] identified the optimal mass flow of geothermal fluid under minimum pumping power and maximum extracted heat energy.Beier et al. [26] characterized and analyzed the vertical temperature profiles of the working fluid of a CGHE.Nian et al. [27] developed a building and CGHE integrated model that comprehensively takes the heat transfer of wellbores, system formation, and building thermal loads into consideration.This model is successfully used to examine critical system performance parameters, including the production of geothermal energy, room temperature, and outlet water temperature.
Currently, the number of AOGWs worldwide has reached approximately 30 million [16,27].In addition, regions with favorable geothermal parameters are inherently rich in geothermal energy [7].The conversion of those petroleum wells into geothermal heat exchangers can substantially reduce the investment costs of borehole drilling [16,24] and enhance the utilization of geothermal resources, while also mitigating environmental risks [28].Therefore, repurposing old petroleum boreholes for geothermal energy extraction is highly desirable in the pursuit of global carbon neutrality.
Despite the extensive research, a comprehensive performance analysis of a CGHE retrofitted from AOGW, particularly in terms of varied borehole diameters, remains lacking.There is also a lack of a computational model for the CGHE with different diameters Most existing studies merely focus on retrofitting processes, ignoring an in-depth examination of factors that affect the heat extraction performance of such retrofitted heat exchangers.Moreover, there is a notable scarcity of insightful research on optimizing critical parameters during the retrofitting process and analyzing the economic performance of these systems.Given these research gaps, this study aims to develop a numerical heat transfer model of the retrofitted CGHE with varied outer pipe diameters based on the finite difference method (FDM) [20].Based on the developed CGHE model, the influence of four critical factors on the thermal performance of CGHE is evaluated.In addition, the optimal flow rate of CGHE is also investigated from an economic standpoint.The findings of this study can further enhance the adoption of deep CGHE from abandoned oil wells and extend the utilization of deep geothermal energy.

Abandoned Oil Well Structure and Working Principles of a CGHE
The oil well usually consists of a steel pipe with varying pipe diameters along the borehole depth.The drilling process and equipment used for oil well engineering are almost identical to those used for CGHE, which makes retrofitting the oil well into CGHE possible.In this study, the abandoned well, which will be retrofitted to CGHE with coaxial tubes, refers to an oil well that no longer has economic potential after several years of oil production.Compared with a conventional CGHE, the main distinction of the retrofitted CGHE derived from an abandoned well lies in its decreasing diameter of the outer pipe along the borehole depth.Thus, developing a CGHE with varied outer pipe diameters for different depth regions is critical to this transformation.There are various types of retrofitted CGHE, and currently, the three-section CGHE and four-section CGHE are the most commonly used types and are more suitable for heat extraction.Figure 1 describes a cross-section structure of the retrofitted CGHE with four sections along the borehole depth.As shown, the diameter of the outer pipe decreases along the depth, while the diameter of the inner pipe remains constant.

Abandoned Oil Well Structure and Working Principles of a CGHE
The oil well usually consists of a steel pipe with varying pipe diameters along the borehole depth.The drilling process and equipment used for oil well engineering are almost identical to those used for CGHE, which makes retrofitting the oil well into CGHE possible.In this study, the abandoned well, which will be retrofitted to CGHE with coaxial tubes, refers to an oil well that no longer has economic potential after several years of oil production.Compared with a conventional CGHE, the main distinction of the retrofitted CGHE derived from an abandoned well lies in its decreasing diameter of the outer pipe along the borehole depth.Thus, developing a CGHE with varied outer pipe diameters for different depth regions is critical to this transformation.There are various types of retrofitted CGHE, and currently, the three-section CGHE and four-section CGHE are the most commonly used types and are more suitable for heat extraction.Figure 1 describes a crosssection structure of the retrofitted CGHE with four sections along the borehole depth.As shown, the diameter of the outer pipe decreases along the depth, while the diameter of the inner pipe remains constant.
Due to the high ground temperature, the deep CGHE is only used for heat extraction to provide space heating for buildings.Under heat extraction conditions, low-temperature water flows into the CGHE through the outer pipe.During this process, the circulating water is heated by the surrounding ground.As the water reaches the bottom of the CGHE, it reaches its highest temperature and flows into the inner pipe.One advantage of this system is that it operates as a closed circulation system, meaning that it does not extract groundwater directly during the heat extraction process, which can help protect the groundwater environment to some extent.

Model Assumptions
To improve the calculation efficiency while maintaining the prediction accuracy, the following assumptions are adopted: 1.
The thermal properties of ground, pipe, water, and backfilling material are dependent on their specific temperatures.

2.
The properties of the ground are homogeneous in each ground layer.Due to the high ground temperature, the deep CGHE is only used for heat extraction to provide space heating for buildings.Under heat extraction conditions, low-temperature water flows into the CGHE through the outer pipe.During this process, the circulating water is heated by the surrounding ground.As the water reaches the bottom of the CGHE, it reaches its highest temperature and flows into the inner pipe.One advantage of this system is that it operates as a closed circulation system, meaning that it does not extract groundwater directly during the heat extraction process, which can help protect the groundwater environment to some extent.

Model Assumptions
To improve the calculation efficiency while maintaining the prediction accuracy, the following assumptions are adopted: 1.
The thermal properties of ground, pipe, water, and backfilling material are dependent on their specific temperatures.

2.
The properties of the ground are homogeneous in each ground layer.

3.
The penetration of groundwater is neglected.

4.
The lateral and bottom boundaries of the borehole model are set to be temperature constant.

5.
The geothermal heat flux in the ground is assumed to be uniform and constant.

Initial and Boundary Conditions
The initial temperature of the ground has a uniform geothermal gradient along the depth direction in different ground layers.As the geothermal heat flux is assumed to be uniform and constant, the initial ground temperature at any depth in the stratum can be calculated by Equation (1) [20]: where t a is the ambient air temperature ( • C), h a is the convective heat transfer coefficient at the ground surface, (W/(m 2 •K)), H j is the bottom coordinate of the j-layer, q g is the geother- mal heat flux (W/m 2 ), k is the ground thermal conductivity (W/(m•K)); q g h a represents the temperature responses induced by q g at the ground surface, with a convective heat transfer coefficient of h a ; and q g k j represents the temperature responses induced by q g within the j-layer ground, with a conductive heat transfer coefficient of k j .
The bottom boundary of the whole calculation region is set far below the bottom of the borehole, with an extra distance of 200 m beneath the borehole bottom.As no thermal interference occurred at this bottom boundary, the temperature remains constant.
The heat convection boundary condition is set at the ground surface, which can be expressed by Equation (2) [20].
Although the assumption of constant temperature slightly deviates from the actual conditions, which are influenced by variations in ground surface temperature, this impact can be ignored.Considering the drilling depth of more than 2000 m, the model does not need to introduce more complex boundary conditions.
For the two partial differential equations of fluid temperature distribution in the inner and outer tubes, it is considered that the initial temperature of the fluid in the tube is the same as that of the ground at the same depth, which can also be expressed by Equation (1).
The boundary conditions of the partial differential Equation (3) [20] of the flow direction of the tube is where Q (W) is the total heat transfer rate of the CGHE, and c (J/(kg•K) is the specific heat capacity of the circulating fluid.

Heat Transfer Governing Equation
Based on the above assumptions, the heat transfer governing equations of the CGHE are described as follows.Since each layer of ground is axially symmetric, the heat transfer process of each layer of the ground can be regarded as a two-dimensional unsteady heat transfer in the cylindrical coordinate [20]: In the case of heat extraction, the heat flux exchanged between the CGHE and ground peaks at the borehole wall and decreases radially outward.As a result, the ground temperature gradually decreases and can hardly be detected at the radial boundary.Considering heat transfer characteristics in the ground, the varied space steps in the radial direction are employed to improve the computing efficiency, as shown in Figure 2.
In the case of heat extraction, the heat flux exchanged between the CGHE and ground peaks at the borehole wall and decreases radially outward.As a result, the ground temperature gradually decreases and can hardly be detected at the radial boundary.Considering heat transfer characteristics in the ground, the varied space steps in the radial direction are employed to improve the computing efficiency, as shown in Figure 2. Thus, a new radial coordinate is adopted, which can be described as follows: If we apply the new radial coordinate in the control Equation ( 6), it can be rewritten as follows: The flow pattern of the CGHE in this study is the fluid flowing down from the sandwich of two tubes and up from the inner tube.The energy differential equation of the fluid in the inner tube is The energy differential equation for the fluid in the outer tube is where C = Mc is the heat capacity (J/(m•K)) of the circulating fluid;  and  are the density and specific heat capacity of the circulating fluid;  and  are the density and specific heat of the outer pipe wall; and  and  are the density and specific heat of the inner pipe wall.
Considering the inner pipe and circulating water inside the inner pipe, the specific heat capacities should be calculated by Thus, a new radial coordinate is adopted, which can be described as follows: If we apply the new radial coordinate in the control Equation ( 6), it can be rewritten as follows: 1 The flow pattern of the CGHE in this study is the fluid flowing down from the sandwich of two tubes and up from the inner tube.The energy differential equation of the fluid in the inner tube is The energy differential equation for the fluid in the outer tube is where C = Mc is the heat capacity (J/(m•K)) of the circulating fluid; ρ f and c f are the density and specific heat capacity of the circulating fluid; ρ f 1 and c f 1 are the density and specific heat of the outer pipe wall; and ρ f 2 and c f 2 are the density and specific heat of the inner pipe wall.
Considering the inner pipe and circulating water inside the inner pipe, the specific heat capacities should be calculated by where C 1 is the per-length heat capacity of the outer channel, including the outer pipe, circulating water, and backfilling material (J/m•K); and C 2 is the per-length heat capacity of the inner channel.Since the model divides the CGHE into four sections, the thermal resistances of the four sections can be calculated according to different radial boundaries.The thermal resistances between the circulating water in the outer pipe and the borehole wall, R 1 , are given as follows: R 2 is the thermal resistance between the circulating water in the inner and outer pipes, which can be calculated by (11) where k g , k p1 , and k p2 are the thermal conductivity of backfill material, outer pipe, and inner pipe, respectively; d o and d i represent the pipe's outer diameter and inner diameter, respectively; and d b represents the borehole diameter.
Moreover, h 1 and h 2 are the convective heat transfer coefficients of the annular pipe and inner pipe, which can be calculated by Nu is the dimensionless criterion number of the reaction on the strength of the heat transfer process, which can be obtained from [29,30] Nu = 0.023Re 0.8 Pr 0.333 Re ≥ 10 4 (13) Nu = 0.116(Re 0.667 − 125)Pr 0.333 Re > 2200 (15) where L is the length of the borehole in each section, m; and Pr is the dimensionless number of the effect of reactive fluid property on heat transfer.Re can be calculated by where u is the average section velocity, m/s; d is the equivalent diameter, m; and ν is the kinematic viscosity, m 2 /s.

Numerical Discretization
The heat transfer model for the deep CGHE can be considered a two-dimensional unsteady heat conduction problem with specific convective heat transfer boundary conditions on the borehole wall.This study employs the alternative direction finite difference method (ADFDM) [31].The ADFDM combines the forward and backward FDM and uses them alternatively in adjacent time steps.This approach ensures that no more than three unknowns need to be solved for each numerical node, which significantly improves the calculation efficiency.
The model assumes that there is no influence from groundwater flow.Based on the geometric characteristics of the CGHE, the FDM is utilized to establish the computational grid division for the efficient simulation model of the heat exchanger in abandoned oil wells.The grid division is shown in Figure 2.
In the case of converting an abandoned well into a CGHE, there is a variation in the diameter of the outer pipe of the ground heat exchanger along the depth of the borehole.The calculation process differs from the conventional medium-depth CGHE, where the outer tube diameter remains constant.To address this, the fluid temperature in each section of the pipe with the same diameter needs to be calculated separately, as illustrated in Figure 3. unknowns need to be solved for each numerical node, which significantly improves the calculation efficiency.The model assumes that there is no influence from groundwater flow.Based on the geometric characteristics of the CGHE, the FDM is utilized to establish the computational grid division for the efficient simulation model of the heat exchanger in abandoned oil wells.The grid division is shown in Figure 2.
In the case of converting an abandoned well into a CGHE, there is a variation in the diameter of the outer pipe of the ground heat exchanger along the depth of the borehole.The calculation process differs from the conventional medium-depth CGHE, where the outer tube diameter remains constant.To address this, the fluid temperature in each section of the pipe with the same diameter needs to be calculated separately, as illustrated in Figure 3. Based on the ADFDM, the numerical equations for fluid in inner and outer pipes can be obtained.It should be noted that the fluid in the pipe is assumed to be a one-dimensional model, indicating that the water temperature in the pipe on the horizontal level is uniform.The numerical equations of fluid in the outer pipe are as follows: where The numerical equations of fluid in the inner pipe are as follows: Based on the ADFDM, the numerical equations for fluid in inner and outer pipes can be obtained.It should be noted that the fluid in the pipe is assumed to be a one-dimensional model, indicating that the water temperature in the pipe on the horizontal level is uniform.The numerical equations of fluid in the outer pipe are as follows: where . The numerical equations of fluid in the inner pipe are as follows: where where B r = a∆τ (r∆σ) 2 , B z = a∆τ (∆z) 2 .

Equation Solving Algorithm
The complete difference equation can be broken down into several sets of second-order equations, and each set of algebraic equations has a tridiagonal coefficient matrix.These tridiagonal matrices satisfy the solution rule of the catch-up method.The tridiagonal linear algebraic equations can be expressed as follows [20]: . . .
Since the first equation contains only two unknowns, it can be written as Therefore, when i = 1, the . Using this method, in turn, the general expression of the coefficient is where ; and i = 2, 3, . .., n−1 It can be obtained according to the solution process of the catch-up method. where ; and i = 2, 3, . .., n−1.

Model Validation
This study selected the experimental data from Xi'an's non-interference geothermal heat supply project [32] to validate the proposed numerical model.In this project, the drilling depth is 2500 m, the average geothermal gradient is 0.0283 • C/m, the import flow rate of a single well is 13.97 m 3 /h, and the heating period is from 15 November to the 15 March.The inputs include the inlet fluid temperatures, specifications of the CGHE, thermal physical properties of the ground, ground temperature gradient, mass flow rate, and operation time.The output is the outlet water temperature, which was then used for model validation.In Reference [32], the maximum relative error between the simulated and measured outlet temperatures was 4.9%.It can be considered that the simulation results agree with the measured data.
In this study, the outlet water temperature of the project is simulated by the developed numerical model, with the same input parameters as Reference [32].The comparison results among the numerical simulation results of this study and those of Reference [32] and the measured data are shown in Figure 4.It can be seen from Figure 4 that the proposed numerical result in this study follows a trend similar to both the measured data and the simulated data from Reference [32].The maximum relative temperature difference of the numerical simulation data for the referred model versus the developed model in this manuscript is less than 5%.Therefore, the numerical model is considered reliable for further analysis of the performance of the CGHE.
Sustainability 2024, 16, 1603 9 of 21 and the measured data are shown in Figure 4.It can be seen from Figure 4 that the proposed numerical result in this study follows a trend similar to both the measured data and the simulated data from Reference [32].The maximum relative temperature difference of the numerical simulation data for the referred model versus the developed model in this manuscript is less than 5%.Therefore, the numerical model is considered reliable for further analysis of the performance of the CGHE.

Parameter Settings of the Retrofitted CGHE
The heat transfer performance of the deep CGHE is affected by factors such as the subsurface thermal properties, geometric parameters of the abandoned well, and operational parameters.In this study, an abandoned oil well located in Northern China is investigated as a case study, to be retrofitted into a CGHE for space heating.The abandoned well has four sections with different borehole diameters along the borehole depth.According to the onsite measurement, the average thermal properties of each section are determined by sample measurement, using the hot-disk thermal conductivity tester, and the obtained property data are listed in Table 1.The accuracy of the instrument for the thermal properties is within ±10%.The basic geometric parameters of the well are also listed in Table 1.According to the ground temperature gradient tested in the field and the thermal conductivity of each rock layer, the geothermal heat flux can be approximated to be 0.058 W/m 2 , assuming a constant value.Other parameters are set as follows: the outer pipe is the steel pipe with a thermal conductivity of 43 W/(m•K), the circulating medium in the

Thermal Performance Analysis of the Retrofitted CGHE 4.1. Parameter Settings of the Retrofitted CGHE
The heat transfer performance of the deep CGHE is affected by factors such as the subsurface thermal properties, geometric parameters of the abandoned well, and operational parameters.In this study, an abandoned oil well located in Northern China is investigated as a case study, to be retrofitted into a CGHE for space heating.The abandoned well has four sections with different borehole diameters along the borehole depth.According to the onsite measurement, the average thermal properties of each section are determined by sample measurement, using the hot-disk thermal conductivity tester, and the obtained property data are listed in Table 1.The accuracy of the instrument for the thermal properties is within ±10%.The basic geometric parameters of the well are also listed in Table 1.According to the ground temperature gradient tested in the field and the thermal conductivity of each rock layer, the geothermal heat flux can be approximated to be 0.058 W/m 2 , assuming a constant value.Other parameters are set as follows: the outer pipe is the steel pipe with a thermal conductivity of 43 W/(m•K), the circulating medium in the heat exchanger is water, the average temperature of the surface air is 15 • C, and the thermal conductivity of the backfill material is 1.5 W/(m•K).The inner and outer diameters of the inner pipe are set to be 110/90 mm.

Average Heat Extraction Rate in the Heating Season
The operating performance of the CGHE is related to the geometric and physical characteristics and, to a large extent, the working conditions.It is challenging to evaluate the heat transfer performance of CGHE comprehensively.Therefore, a quantitative index of the average heat extraction rate is defined to facilitate smooth communication between engineers and non-professionals.The average heat extraction rate of the CGHE refers to the mean heat extraction of a deep CGHE during the whole operation period, under specific operating conditions.
where ∑ t 0 ∆T is the sum of the difference between the inlet and outlet water temperatures at all times ( • C); and t is time, h.
In this study, the particular condition refers to the fact that the inlet temperature of the CGHE system shall not be lower than 5 • C during the heating period, and 5 • C, 10 • C, 15 • C, and 20 • C are, respectively, used for the simulation.

Orthogonal Test Scheme Design and Test Results
Design parameters, geological conditions, operating conditions, and other factors will significantly impact the heat transfer performance of the CGHE.Especially for abandoned oil wells, many different transformation schemes will influence the value of the average heat extraction rate.The major difference between the abandoned well and the conventional CGHE is the varied outer pipe diameter, along the borehole depth.The smaller diameter of the fourth section leads to increased flow resistance, so whether to consider the fourth section as the heat exchange section needs comprehensive consideration.In this study, four retrofitted configurations of CGHE are considered, as shown in Table 2.It is noticed that the constant d o is selected to be the diameter of the fourth section.Therefore, this study focuses on the impact of four main factors on heat extraction, namely the configuration of the rebuilt abandoned oil well, the inlet temperature, the mass flow rate, and the internal pipe thermal conductivity.The orthogonal test is thus employed to analyze the influence parameters of the average heat extraction rate [33].The orthogonal test method is widely used in scientific research because of its high efficiency and cost-effectiveness.It is necessary to clarify the significance of the influence of the factors.
The total influence factors of the orthogonal test are 4, and the level number of each factor is set to 4. The orthogonal test table, L 16 (4 4 ), is selected according to the factor level.The test scheme and test solution results corresponding to the level values of each factor are shown in Tables 2 and 3, and the average heat extraction is set as the target parameter.

Analysis of Orthogonal Test Results
In the orthogonal test method, a certain factor may cause all the differences in the target parameters when only considering its influence.The factor which causes the most significant change in the average heat extraction rate is considered to be the main influencing factor.The calculation formula is as follows: where n i is the arithmetic average of orthogonal test results when the level i is taken for the column factors; s is the number of tests for these factors under level i; y n is the index value of the k test; and R is the range value in this column.
The average heat extraction rate range analysis of the CGHE in abandoned well reconstruction is shown in Table 4.The range analysis shows that the configurations of CGHEs retrofitted from abandoned oil wells, inlet temperature, mass flow rate, and inner tube thermal conductivity are 43.9, 124, 22.7, and 12.5, respectively.Therefore, the significance of each factor in affecting the heat transfer performance of the CGHE is listed as follows: inlet water temperature > configuration of CGHE retrofitted from abandoned oil wells > flow rate > inner tube thermal conductivity.
As Figure 5 shows, the average heat extraction rate increases with higher flow rates.However, the average heat extraction rate decreases as the inlet temperature rises, which is attributed to the smaller temperature difference between the circulating fluid and the surrounding ground.Moreover, it also shows that the thermal insulation performance of the inner tube plays a vital role in preventing the thermal short circuit phenomenon between the inner tube and the outer tube and also improving the heat extraction performance of the system.The range analysis shows that the configurations of CGHEs retrofitted from abandoned oil wells, inlet temperature, mass flow rate, and inner tube thermal conductivity are 43.9, 124, 22.7, and 12.5, respectively.Therefore, the significance of each factor in affecting the heat transfer performance of the CGHE is listed as follows: inlet water temperature > configuration of CGHE retrofitted from abandoned oil wells > flow rate > inner tube thermal conductivity.
As Figure 5 shows, the average heat extraction rate increases with higher flow rates.However, the average heat extraction rate decreases as the inlet temperature rises, which is attributed to the smaller temperature difference between the circulating fluid and the surrounding ground.Moreover, it also shows that the thermal insulation performance of the inner tube plays a vital role in preventing the thermal short circuit phenomenon between the inner tube and the outer tube and also improving the heat extraction performance of the system.

Case Study and Economic Analysis
In this section, the economic analysis is performed on the retrofitted CGHE from the abandoned oil well, aiming to disclose its practical feasibility and popularization potential of actual applications.

Actual Operation of the Ground Source Heat Pump (GSHP) System with Retrofitted CGHE
In this study, an abandoned oil well located in Jinan, China, is selected as a case to demonstrate the thermal performance of the GSHP system with a retrofitted CGHE.Based on the above discussion, the configuration of a three-section with a variant do of 2626 m is selected to be the optimal design scheme.The thermal conductivity of the inner pipe is set to be 0.17 W/(m•K).According to the real-time heating load of the concerned building in the project, as shown in Figure 6, the appropriate heat pump unit is selected, and its nominal performance under heating conditions is shown in Table 5.

Case Study and Economic Analysis
In this section, the economic analysis is performed on the retrofitted CGHE from the abandoned oil well, aiming to disclose its practical feasibility and popularization potential of actual applications.

Actual Operation of the Ground Source Heat Pump (GSHP) System with Retrofitted CGHE
In this study, an abandoned oil well located in Jinan, China, is selected as a case to demonstrate the thermal performance of the GSHP system with a retrofitted CGHE.Based on the above discussion, the configuration of a three-section with a variant d o of 2626 m is selected to be the optimal design scheme.The thermal conductivity of the inner pipe is set to be 0.17 W/(m•K).According to the real-time heating load of the concerned building in the project, as shown in Figure 6, the appropriate heat pump unit is selected, and its nominal performance under heating conditions is shown in Table 5.The inlet water temperature variations in the CGHE under different water flow rates are simulated to facilitate the identification of the suitable circulation water flow.As shown in Figure 7, in certain conditions, with the circulating water flow rate of 25 m 3 /h and 30 m 3 /h, the inlet water temperature will drop below 5 • C, which is considered impractical for the actual application.However, when readjusting the circulating water flow to 35 or 40 m 3 /h, the inlet water temperature consistently remains above 5 • C during the whole operation period, demonstrating its high reliability in practical applications.The inlet water temperature variations in the CGHE under different water flow rates are simulated to facilitate the identification of the suitable circulation water flow.As shown in Figure 7, in certain conditions, with the circulating water flow rate of 25 m 3 /h and 30 m 3 /h, the inlet water temperature will drop below 5 °C, which is considered impractical for the actual application.However, when readjusting the circulating water flow to 35 or 40 m 3 /h, the inlet water temperature consistently remains above 5 °C during the whole operation period, demonstrating its high reliability in practical applications.It is obvious that a higher fluid temperature can be obtained with a larger flow rate due to its better heat convection.
According to Equations ( 12)-( 16), the convection heat exchange coefficients of the pipe sections can be calculated by the hydraulic diameter, Re, Pr, and Nu.Taking the flow rate of 25 m 3 /h as an example, Table 6 can be obtained.It is obvious that a higher fluid temperature can be obtained with a larger flow rate due to its better heat convection.
According to Equations ( 12)-( 16), the convection heat exchange coefficients of the pipe sections can be calculated by the hydraulic diameter, Re, Pr, and Nu.Taking the flow rate of 25 m 3 /h as an example, Table 6 can be obtained.Similarly, the same calculation can be performed for other flow rates.As shown in Figure 8a, the heat transfer coefficient of the inner tube is significantly higher than that of the outer annular pipe.Furthermore, due to the decreasing diameter of the outer annular pipe at a greater depth, resulting in increased flow velocity, the heat transfer coefficient of the annular wall also exhibits an upward trend along the depth, reaching the highest value in the third (deepest) section.Table 7 is the mathematical relation of convection heat transfer coefficients and flow rates summarized according to Figure 8a.It also can be seen in Figure 8b that the convective heat transfer resistances of the two pipe sides are much more insignificant compared to the heat conductivity resistance.The surrounding ground temperatures after the operation of 20 days (480 h) and the end of the heating season (2880 h) are also simulated in the case of the flow rate of 35 m 3 /h, as shown in Figure 9.It can be seen from the figures that, with the continuous operation of the CGHE, the ground temperature surrounding the CGHE notably decreases, highlighting the necessity of energy rejection during non-heating to maintain its heating capacity.It also can be seen in Figure 8b that the convective heat transfer resistances of the two pipe sides are much more insignificant compared to the heat conductivity resistance.
The surrounding ground temperatures after the operation of 20 days (480 h) and the end of the heating season (2880 h) are also simulated in the case of the flow rate of 35 m 3 /h, as shown in Figure 9.It can be seen from the figures that, with the continuous operation of the CGHE, the ground temperature surrounding the CGHE notably decreases, highlighting the necessity of energy rejection during non-heating to maintain its heating capacity.

Initial System Investment
(1) Drilling cost The drilling cost of the deep CGHE in Shandong province, China, ranges from 500 to 600 CNY/m.This study estimates the drilling cost at 600 CNY/m.The cost of Cb at various drilling depths is shown in Table 8.

Initial System Investment
(1) Drilling cost The drilling cost of the deep CGHE in Shandong province, China, ranges from 500 to 600 CNY/m.This study estimates the drilling cost at 600 CNY/m.The cost of C b at various drilling depths is shown in Table 8.Because the CGHE in this case is modified from the abandoned oil wells, the drilling cost savings can reach up to CNY 1800 thousand.
(2) Pipe material cost According to the actual engineering material selection planning, the outer pipe material is a J55 special steel pipe with a cost of 400 CNY/m, and the inner pipe is a high-density polyethylene (HDPE) pipe with a cost of 40 CNY/m.This amounts to a total cost of 440 CNY/m for the deep casing-buried pipe.The cost of pipe material for various buried depths is calculated and summarized in Table 9. (1) Water pump cost The net power consumption of the water pump in this study is approximated by calculating the frictional pressure drop and local pressure drop of the CGHE.
The frictional pressure drop caused by the viscous force during the fluid flow is called the drag loss, which can be calculated by the following Equation (28) [29] where λ is the friction factor; l is the tube length, m; and ρ is the density of the circulating fluid.The equation for calculating the friction factor (λ) and Reynolds number (Re) along the path is as follows [29]: where K is the height of the equivalent coarse grain, m; and d is the equivalent diameter, m.The weighted average method can calculate the equivalent roughness height of the annular area of the inner and outer tubes of the casing-buried tube heat exchanger.Therefore, the equivalent roughness height of the annular channel should be as follows [34].
where K 1 is the height of equivalent roughness on the inner wall of the outer tube, m; d 1i is the inner diameter of the outer tube, m; K 2 is the equivalent roughness height on the outer wall of the inner tube, m; and d 2o is the outer diameter of the inner tube, m.
According to the research results [31,32], there is a particular relationship between the local pressure drop and pressure drop along the pipeline.In this study, 2% of the pressure drop along the road is used for analysis.
The total pressure drop is the sum of the friction pressure drop and pressure drop, where ∆p z is the total pressure drop of the heat exchanger, Pa; ∆p f is the friction pressure drop, Pa; and ∆p s is the local pressure drop.
Ignoring the resistance loss of the evaporator pipeline and the buried pipeline, the net power consumption of the circulating pump can be approximated to [29] P = ∆p z × q f (32) where P is the net power consumption of the pump, W; and q f is the mass flow rate of circulating water, kg/s.Meanwhile, water pump efficiency is always a necessary concern in the regulation of the GSHP system with retrofitted CGHE.Thus, the water pump efficiency can be calculated by Equation (33) [35]: where g is the acceleration of gravity, m/s 2 ; H pump is the pump lift, m; and P pump is the rated power of the pump, kW, set to 18.5 kW.
Figure 10 shows the change in total pressure drops and pump efficiency in each section pipe under different flow rates.According to Figure 10, the mathematical relationship formula between the flow rate and total pressure drop can be obtained via linear fitting, which is shown in Table 10.(2) Heat pump operation cost The GSHP coefficient of performance (COP) represents the conversion efficiency between energy and heat, specifically the ratio of the heat energy produced by the heat pump to the electrical energy consumed.Figure 11a illustrates the fitted curve of the heat pump inlet water temperature and COP values.Figure 11b indicates the variations in heat pump COP with flow rates, demonstrating that the heat pump has a favorable heat extraction effect with larger flow rates.The GSHP exhibits the highest value of COP when the flow rate is 40 m 3 /h.(2) Heat pump operation cost The GSHP coefficient of performance (COP) represents the conversion efficiency between energy and heat, specifically the ratio of the heat energy produced by the heat pump to the electrical energy consumed.Figure 11a illustrates the fitted curve of the heat pump inlet water temperature and COP values.Figure 11b indicates the variations in heat pump COP with flow rates, demonstrating that the heat pump has a favorable heat extraction effect with larger flow rates.The GSHP exhibits the highest value of COP when the flow rate is 40 m 3 /h.tween energy and heat, specifically the ratio of the heat energy produced by the heat pump to the electrical energy consumed.Figure 11a illustrates the fitted curve of the heat pump inlet water temperature and COP values.Figure 11b indicates the variations in heat pump COP with flow rates, demonstrating that the heat pump has a favorable heat extraction effect with larger flow rates.The GSHP exhibits the highest value of COP when the flow rate is 40 m 3 /h.According to the latest implementation standard of the electricity price in Jinan in the peak, valley, and flat value, the average industrial electricity price is 0.6 CNY/(kW•h).The electricity cost required in one heating season (120 d) can be calculated based on the COP, as shown in Figure 12.As can be seen from the figure, the operation cost of CGHE retrofitted three-section abandoned oil wells.With an increase in the flow rate, the heat pump operation cost decreases rapidly.According to the latest implementation standard of the electricity price in Jinan in the peak, valley, and flat value, the average industrial electricity price is 0.6 CNY/(kW•h).The electricity cost required in one heating season (120 d) can be calculated based on the COP, as shown in Figure 12.As can be seen from the figure, the operation cost of CGHE retrofitted three-section abandoned oil wells.With an increase in the flow rate, the heat pump operation cost decreases rapidly.In summary, Figure 13 enumerates that there are three kinds of total annual operation costs under four different flow rates, namely water pump operation cost, heat pump operation cost, and total operation cost.To guarantee the high operation efficiency of the whole system, it is advisable to maintain a flow rate of no less than 35 m 3 /h.Meanwhile, the larger flow rate can improve the convection heat transfer coefficient, but it also increases the operation cost of the water pump.In summary, Figure 13 enumerates that there are three kinds of total annual operation costs under four different flow rates, namely water pump operation cost, heat pump operation cost, and total operation cost.To guarantee the high operation efficiency of the whole system, it is advisable to maintain a flow rate of no less than 35 m 3 /h.Meanwhile, the larger flow rate can improve the convection heat transfer coefficient, but it also increases the operation cost of the water pump.
From an economic perspective, 35 m 3 /h is selected as the optimal flow rate, which ensures a relatively smaller distribution cost while maintaining a sufficiently high inlet water temperature of the CGHE to avoid system icing.
In summary, Figure 13 enumerates that there are three kinds of total annual operation costs under four different flow rates, namely water pump operation cost, heat pump operation cost, and total operation cost.To guarantee the high operation efficiency of the whole system, it is advisable to maintain a flow rate of no less than 35 m 3 /h.Meanwhile, the larger flow rate can improve the convection heat transfer coefficient, but it also increases the operation cost of the water pump.From an economic perspective, 35 m 3 /h is selected as the optimal flow rate, which ensures a relatively smaller distribution cost while maintaining a sufficiently high inlet water temperature of the CGHE to avoid system icing.

Conclusions
This study demonstrates the practicality of repurposing abandoned wells into deep CGHEs, making use of the existing structures of both three-section and four-section abandoned oil wells.A numerical heat transfer model for CGHEs with varying outer pipe diameters is established, employing a finite difference method to handle complex calculations.The accuracy of the developed model is validated with experiment data, ensuring an accurate analysis of the heat transfer performance of the retrofitted heat exchangers from abandoned wells.
(1) A sensitive analysis was carried out by using the orthogonal experimental method to investigate the influence of four critical variables on the thermal performance of CGHE.The results show that the inlet water temperature plays the most important

Conclusions
This study demonstrates the practicality of repurposing abandoned wells into deep CGHEs, making use of the existing structures of both three-section and four-section abandoned oil wells.A numerical heat transfer model for CGHEs with varying outer pipe diameters is established, employing a finite difference method to handle complex calculations.The accuracy of the developed model is validated with experiment data, ensuring an accurate analysis of the heat transfer performance of the retrofitted heat exchangers from abandoned wells.
(1) A sensitive analysis was carried out by using the orthogonal experimental method to investigate the influence of four critical variables on the thermal performance of CGHE.The results show that the inlet water temperature plays the most important role in heat extraction rates, followed by types of CGHE retrofitted from abandoned wells, flow rate, and inner tube thermal conductivity.(2) An economic analysis of the retrofit, considering the operation costs of the water pump and heat pump, was performed.Given the fact that the CGHE in this study was retrofitted from the abandoned oil wells, the drilling cost can be reduced by up to CNY 1800 thousand.The pump operation cost of CGHE retrofitted from three-section abandoned oil wells is significantly increase with the increase in the flow rate, while the heat pump operation cost decreases because the COP increases.
(3) This study integrates actual project considerations to determine the most appropriate circulation flow based on building load requirements.From an economic perspective, a flow rate of 35 m³/h is selected as the optimal choice, which aligns with both the economic efficiency and operation stability of the middle and deep CGHEs.
In future research, the transformation project of the abandoned wells should be considered, and the interaction between the Wells should be deeply studied.Among them, the heat radius is the problem we must consider.

Figure 1 .
Figure 1.Schematic diagram of CGHE with varied outer pipe diameters.(The colors in the figure represent the water temperature, which is red when it is higher and blue when it is lower).

Figure 1 .
Figure 1.Schematic diagram of CGHE with varied outer pipe diameters.(The colors in the figure represent the water temperature, which is red when it is higher and blue when it is lower).

Figure 2 .
Figure 2. Numerical discretization of the CGHE.(The colors in the figure represent the water temperature, which is red when it is higher and blue when it is lower).

Figure 2 .
Figure 2. Numerical discretization of the CGHE.(The colors in the figure represent the water temperature, which is red when it is higher and blue when it is lower).

Figure 3 .
Figure 3. Flowchart to calculate the circulating fluid temperature of the CGHE with varying outer pipe diameters.

Figure 3 .
Figure 3. Flowchart to calculate the circulating fluid temperature of the CGHE with varying outer pipe diameters.

Figure 4 .
Figure 4. Comparison of simulated and measured outlet temperatures with the [32].

Figure 4 .
Figure 4. Comparison of simulated and measured outlet temperatures with the [32].

Figure 5 .
Figure 5. Analysis of average heat extraction of orthogonal test.

Figure 5 .
Figure 5. Analysis of average heat extraction of orthogonal test.

Figure 6 .
Figure 6.The space-heating load profile during the whole winter period.

Figure 6 .
Figure 6.The space-heating load profile during the whole winter period.

Figure 8 .
Figure 8.Heat transfer coefficients and resistances of the CGHE under different flow rates.(a) Convective heat transfer coefficient.(b) Heat transfer resistances per unit length of CGHE.

Figure 8 .
Figure 8.Heat transfer coefficients and resistances of the CGHE under different flow rates.(a) Convective heat transfer coefficient.(b) Heat transfer resistances per unit length of CGHE.

21 Figure 10 .
Figure 10.Total pressure drops of the CGHE under different flow rates.

Figure 11 .
Figure 11.Variation in heat pump COP with flow rates.(a) Linear fit curve of COP.(b) The COP variations in GSHP with different flow rates.

Figure 11 .
Figure 11.Variation in heat pump COP with flow rates.(a) Linear fit curve of COP.(b) The COP variations in GSHP with different flow rates.

21 Figure 12 .
Figure 12.Variation in heat pump operation cost with flow rates.

Figure 12 .
Figure 12.Variation in heat pump operation cost with flow rates.

Figure 13 .
Figure 13.Operation costs of GSHP with retrofitted 3-section abandoned oil well.

Figure 13 .
Figure 13.Operation costs of GSHP with retrofitted 3-section abandoned oil well.

Table 1 .
Basic parameters of the four different sections.

Table 1 .
Basic parameters of the four different sections.

Table 2 .
Schemes of orthogonal test.

Table 3 .
Results of orthogonal test.

Table 4 .
Analysis of average heat extraction rate of orthogonal test.

Table 5 .
Technical performance parameter table of the heat pump unit.

Table 5 .
Technical performance parameter table of the heat pump unit.

Table 5 .
Technical performance parameter table of the heat pump unit.

Table 7 .
Mathematical relation of convection heat transfer coefficients and flow rates.

Table 8 .
Drilling cost of CGHE with different borehole depths.

Table 8 .
Drilling cost of CGHE with different borehole depths.

Table 9 .
Cost of pipe material of CGHE with different borehole depths.Total Operation Cost of the GSHP System with Retrofitted CGHE This study uses the industrial electricity price in Jinan, China (116.98 • , 36.67 • ), to calculate the operation cost of the water pump and heat pump.