Next Article in Journal
Future Projection of Water Resources of Ruzizi River Basin: What Are the Challenges for Management Strategy?
Previous Article in Journal
Influence of Marine Currents, Waves, and Shipping Traffic on Sulina Channel Fairway at the Mouth of the Black Sea
Previous Article in Special Issue
Comprehensive Study on the Electrical Characteristics and Full-Spectrum Tracing of Water Sources in Water-Rich Coal Mines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Solution for Characterizing Mine Water Recharge of Water Source Heat Pump in Abandoned Coal Mines

1
Inner Mongolia Research Institute, China University of Mining and Technology-Beijing, Ordos 017010, China
2
National Engineering Research Center of Coal Mine Water Hazard Controlling, China University of Mining and Technology, Beijing 100083, China
3
Key Laboratory of Shallow Geothermal Energy, Ministry of Natural Resources of the People’s Republic of China, Beijing 100195, China
4
Key Laboratory of Mine Water Hazard Controlling, National Mine Safety Administration, Beijing 100083, China
5
College of Geoscience and Surveying Engineering, China University of Mining and Technology, Beijing 100083, China
6
School of Energy and Mining Engineering, China University of Mining and Technology, Beijing 100083, China
*
Authors to whom correspondence should be addressed.
Water 2024, 16(19), 2781; https://doi.org/10.3390/w16192781
Submission received: 23 August 2024 / Revised: 23 September 2024 / Accepted: 25 September 2024 / Published: 30 September 2024
(This article belongs to the Special Issue Innovative Technologies for Mine Water Treatment)

Abstract

:
Due to tremendous mining operations, large quantities of abandoned mines with considerable underground excavated space have formed in China during the past decades. This provides huge potential for geothermal energy production from mine water in abandoned coal mines to supply clean heating and cooling for buildings using heat pump technologies. In this study, an analytical model describing the injection pressure of mine water recharge for water source heat pumps in abandoned coal mines is developed. The analytical solution in the Laplace domain for the injection pressure is derived and the influences of different parameters on the injection pressure are investigated. This study indicates that a smaller pumping rate results in a smaller injection pressure, while smaller values of the hydraulic conductivity and the thickness of equivalent aquifer induce larger injection pressures. The well distance has insignificantly influenced the injection pressure at the beginning, but a smaller well distance leads to a larger injection pressure at later times. Additionally, the sensitivity analysis, conducted to assess the behavior of injection pressure with concerning changes in each input parameter, shows that the pumping rate and the hydraulic conductivity have a large influence on injection pressure compared with other parameters.

1. Introduction

Over the past several decades, the massive consumption of fossil fuels has caused excessive greenhouse gas emissions around the whole world, leading to severe environmental problems such as global warming and the frequent occurrence of extreme weather [1]. In order to tackle climate change and its negative impacts, it is critical to take firm action to increase the proportion of clean energy while weakening the share of traditional fossil fuels [2,3,4,5]. A long-term goal is that the proportion of renewable energy will increase to 32% by the year of 2030, as set by the European Union [6]. The shares of electricity generated by renewable energy and non-water renewable energy have accounted for 19.35% and 10.68%, respectively, which far surpassed the prediction proportion for 2035 performed by the U.S. Energy Information Administration [7]. For a long time, coal has been China’s main energy source, accounting for over 50% of the primary energy consumption and production structure, but the utilization of enormous quantities of coal has induced many problems, including large amounts of greenhouse gas emissions and various environmental issues. According to the Paris Climate Agreement, national governments are required to increase the share of clean energy and reduce reliance on fossil fuels, particularly coal [8]. Under these new circumstances, it is estimated that the total energy consumption of China will reach a peak with 5.8 billion tons of standard coal and the CO2 emissions will reach 10.58 billion tons by the year of 2030 [9].
With the proposal of a carbon emissions peak and carbon neutrality in China, the adoption of clean energy from abandoned coal mines has gained tremendous interest in recent years [8,10,11,12,13], especially when extracting geothermal energy from the mine water in abandoned coal mines, so as to provide heating or cooling for buildings by adopting heat pump technologies (as shown in Figure 1) [14,15,16,17]. During coal exploitation, huge amounts of underground excavated voids (goafs, tunnels, shafts, roadways, etc.) are formed. The closure of coal mines results in the termination of the typical drainage work of mine water. The end of the draining of mine water will gradually induce so-called groundwater rebound [18,19]. These underground excavated voids will be filled with large quantities of mine water, which will assimilate and store the heat released from the interior layer of earth and provide great opportunities for geothermal energy extraction [20,21].
In China, around one-third of mines are water-abundant and most coal mines are located below 800–1000 m, which is greatly compatible with the depth of hydrothermal geothermal energy [19]. Additionally, China closed 5300 coal mines from 2016 to 2020 [19] and it has been estimated that the number of closed coal mines will reach 15,000 by the year 2030 [10]. According to incomplete statistics, the closure of coal mines provided approximately 8 × 107 m3 of underground excavated space resources between 2016 and 2020 [22]. Therefore, tremendous underground mine water storage space from abandoned coal mines in China will provide huge development and utilization possibilities for geothermal energy extraction. Most importantly, the exploitation of geothermal energy from abandoned coal mines not only can avoid problems of safety and mine water environmental pollution, but can also reduce the consumption of fossil fuels, which will lead to the achievement of carbon emission reductions [15,19].
As early as the 1980s, the innovative concept of extracting geothermal energy from mine water in abandoned mines was proposed in Springhill, Canada [23]. In the past decades, the exploitation of geothermal energy from abandoned coal mines has been widely used [14,24,25,26,27,28]. Madiseh and Abbasy analyzed the feasibility of closed geothermal heat exchangers installed in underground mine stopes and evaluated the performance of heat transport for such systems [29]. Chudy investigated the geothermal energy possibilities of mine water in the Nowa Ruda region of Poland, employing an open-loop heat pump system with a rejection system. Their findings indicate that this system has the capacity to generate 10 GWh of geothermal energy, effectively reducing carbon dioxide emissions by 4.05 tons [30]. Meanwhile, Wang et al. introduced a new case of an actual abandoned coal mine in Shandong, China to investigate the intricate flow and heat transfer dynamics within the interconnected network of roadways [31]. Also, further analysis was conducted on heat recovery potential under varying operational conditions. Furthermore, Jardon et al. performed an investigation into the potential recovery of geothermal energy from mine water’s thermal value by leveraging water-to-water heat pumps. Their findings indicated that mine water reservoirs in Central Asturias possess an annual energy supply capacity approximating 260,000 thermal MWh [32]. Furthermore, in order to gain a comprehensive overview of geothermal recovery from abandoned mines, Ramos et al. illustrated the implementation procedures and established minimum requirements for heat pump installations in abandoned mines within the Upper Harz region based on 18 projects worldwide [15]. Moreover, based on unique non-isothermal and non-iso-solutal hydrodynamics, Bao and Liu formulated a scientific framework for modeling transient heat extraction from mine water in a single mine shaft through an open-loop system [33]. The simulations explain that the impact of heat extraction on thermohaline stratification stability is negligible under normal pumping rates. Al-Habaibeh et al. investigated the feasibility of extracting energy from an abandoned coal mine by employing a single mine shaft ground source heat pump (GSHP) system [17]. The performance analysis demonstrated that this new system is an effective means of providing sustainable energy for the heating of buildings and can maintain a high and stable coefficient of performance (COP).
Additionally, a considerable amount of theoretical models have been developed to address the issues of flow and heat transfer in many studies. Based on the recurrent calculation algorithm, Rodriguez and Diaz presented a semi-empirical analytical model which can be employed to assess the temperature increment of a fluid flowing through mine galleries [34]. The numerical modellings of fluid flow and heat transfer to investigate the hydraulic and thermal behavior in an open vertical shaft associated with a flooded mine were carried out by Hamm and Sabet [35]. Raymond and Therrien developed a three-dimensional numerical model for simulating heat advection, heat conduction, and mechanical heat dispersion in the underground tunnels and mining shafts of a flooded mine [36]. Guo et al. presented a numerical model incorporating key parameters that can be used to evaluate geothermal resources within abandoned coal mines [16]. The calculation results suggested the importance of fracture permeability and fractured zone height in assessing the geothermal energy potential of similar mine reservoirs. Sun et al. investigated the effects of heat transfer on the heat exchange system of the pipes at an abandoned mine by establishing a flow-heat coupling numerical model. The results showed that the angle of the single-row borehole has a remarkable influence on the heat exchange performance of the system [37]. A 3D numerical calculation model of multi-physical filed coupling for an abandoned mine water source heat pump system based on COMSOL 6.2 Multiphysics software was developed by Zhang et al. [38]. Based on the established model, thermal breakthrough for 14 distinct pumping-recharge well layout schemes was numerically simulated. The results showed that the pumping flow rate, well spacing, hydraulics, and aquifer thickness have significant impacts on thermal breakthrough.
To the best of our knowledge, previous studies concerning geothermal energy extraction from mine water in abandoned coal mines predominantly focused on the suitability of certain locations, technical feasibility, heat transfer modeling, and design optimization. However, studies on the injection pressure of mine water recharge for water source heat pumps in abandoned coal mines remain quite limited. Therefore, the objective of this research is to propose an analytical model for describing the injection pressure of mine water recharge for geothermal energy extraction systems in an abandoned coal mine. This could have theoretical significance for the implementation of geothermal energy extraction from mine water in these abandoned coal mines.

2. Mathematical Model and Analytical Solution

2.1. Mathematical Model

Figure 2 depicts a schematic of a mathematical model for an open-loop geothermal energy extraction system in an abandoned coal mine. To simplify the analysis of the recharge of mine water for such a system, several hypotheses in this study are made as below: (1) the goaf space can be considered as the porous media (equivalent to an aquifer), by the same assumption as in Guo et al. [16]; (2) the equivalent aquifers are regarded as homogeneous and are of a uniform thickness; (3) mine water flow in the equivalent aquifer obeys Darcy’s equation; (4) the release and storage of water induced by changes in hydraulic heads are instantaneous.

2.2. Derivation of the Analytical Solution

Sen [39] proposed that the continuity equation in an annular between radius r and r + Δ r in a time increment Δ t can be represented as follows:
Δ t Q r , t Q r + Δ r , t = 2 S s π d r Δ r s r , t s r , t + Δ t
in which Q is the pumping rate [L3T−1], s refers to the drawdown [L], S s is the specific storage [L−1], r is the radical coordinate [L], and t denotes time [T].
Equation (1) as Δ t 0 and Δ r 0 becomes the following:
Q r , t r = 2 S s π d r s r , t t
where d is the thickness of an equivalent aquifer [L].
The pumping rate for any concentric cylindrical surface with radius r can be defined as follows [39]:
Q = 2 π r d q
in which q refers to the specific discharge [LT−1].
Then, substituting Equation (3) into Equation (2), the final form of the continuity equation can be given as follows:
q r , t r + 1 r q r , t = S S s r , t t
Recalling the property of q r , t = s r and substituting it into Equation (4), the governing equation for characterizing groundwater flow toward the production or injection well in the underground water reservoir can be expressed as follows:
K 2 s r , t r 2 + 1 r s r , t r = S s s r , t t
in which K is the hydraulic conductivity [LT−1].
The boundary conditions can be described as follows, respectively:
s r , 0 = 0
and
s , t = 0
For an open-loop geothermal energy extraction system in an abandoned coal mine, the well radius is relatively large in practical engineering implementation and the wellbore storage cannot be ignored. Thus, similarly as in Papadopulos and Cooper [40], the wellbore storage is considered in this study and the well boundary condition is given as follows:
lim r r w r s r , t r = r w 2 2 K d s w r , t t Q 2 π K d
where rw denotes the well radius [L].
Applying the Laplace transformation to Equation (5), one can obtain the following:
d 2 s ¯ r , p d r 2 + 1 r d s ¯ r , p d r S s p K s ¯ r , p = 0
in which and s ¯ and p refer to the transient drawdown in the Laplace domain and the Laplace variable, respectively.
Then, in order to simplify Equation (9) with defining α = S s p K , Equation (9) can be reduced to:
d 2 s ¯ r , p d r 2 + 1 r d s ¯ r , p d r α s ¯ r , p = 0
Equation (10) is a linear differential Bessel equation with second order, and the general solution of Equation (10) can be given as follows:
s ¯ r , p = C 1 I 0 α r + C 2 K 0 α r
where I 0 · and K 0 · denote the modified Bessel functions of the first and second kind with the order of 0, respectively. The integration constants C1 and C2 can be obtained based on boundary conditions Equations (7) and (8).
Similarly, when the Laplace transformation is applied to boundary conditions Equation (7), one can obtain:
s ¯ , p = 0
Substituting Equation (12) to Equation (11), one obtains:
C 1 = 0
Equation (11) can be then rewritten as:
s ¯ r , p = C 2 K 0 α r
Also, by employing the Laplace transformation to Equation (8), one obtains:
r w d s ¯ r w , p d r = r w 2 2 K d p s ¯ w r w , p Q 2 π K d 1 p
Taking the derivative of Equation (14), one obtains:
d s ¯ r , p d r = C 2 α K 1 α r
The integration constant C2 can then be determined by combining Equations (15) and (16) as follows:
C 2 = Q π p p r w 2 K 0 α r w + 2 d r w K α K 1 α r w
Substituting Equation (17) to Equation (14), one can obtain the analytical solution of transient drawdown for the production well in the Laplace domain:
s ¯ 1 r 1 , p = Q K 0 α r 1 π p p r w 2 K 0 α r w + 2 d r w K α K 1 α r w
where r1 denotes the radial distance from a certain location in the aquifer to axis of the production well [L] and s ¯ 1 is the transient drawdown of the production well in the Laplace domain.
Similarly, the according analytical solution of transient drawdown for the injection well in the Laplace domain can also be given as follows:
s ¯ 2 r 2 , p = Q K 0 α r 2 π p p r w 2 K 0 α r w + 2 d r w K α K 1 α r w
where r2 is the radial distance from the axis of injection well to a certain location in the equivalent aquifer. s ¯ 2 refers to the transient drawdown of the injection well in the Laplace domain.
Then, applying the superposition principle and combining Equations (18) and (19), the transient drawdowns induced by the production and injection well in the equivalent aquifer can be given as follows:
s r , t = s 1 r 1 , t + s 2 r 2 , t
Applying the Laplace transform to Equation (20), one can obtain:
s ¯ r , p = s ¯ 1 r 1 , p + s ¯ 2 r 2 , p
Substituting Equations (18) and (19) to Equation (21), one obtains:
s ¯ r , p = Q K 0 α r 1 K 0 α r 2 π p p r w 2 K 0 α r w + 2 d r w K α K 1 α r w
The transient drawdown, same as in Liu et al. [41], can be expressed with pressure as:
s r , t = f 0 f r , t ρ g
where f and f0 are the pressure and initial pressure in the equivalent aquifer [ML−1T−2], respectively. ρ refers to the density of water [ML−3], and g denotes the gravitational acceleration [LT−2].
Similarly, applying the Laplace transformation to Equation (23), one obtains:
s ¯ r , p = f 0 p f ¯ r , p ρ g p
where f ¯ denotes to the injection pressure in the Laplace domain.
Substituting Equation (22) to Equation (24), the analytical solution of pressure in the Laplace domain can be described as follows:
f ¯ r , p = f 0 p ρ g Q K 0 α r 1 K 0 α r 2 π p p r w 2 K 0 α r w + 2 d r w K α K 1 α r w
In order to calculate the pressure in the injection well, one has r1 = l and r2 = rw in Equation (25). Then, Equation (25) becomes:
f ¯ r , p = f 0 p ρ g Q K 0 α l K 0 α r w π p p r w 2 K 0 α r w + 2 d r w K α K 1 α r w
where l represents the distance between the production and injection well [L].
As presented in Equation (26), the term of the Bessel function K 0 · makes it difficult to obtain the time domain–analytical solution of Equation (26). However, several numerical inverse methods, such as Stehfest [42,43]; Crump [44]; and De Hoog et al. [45], can be applied to address this problem. In this study, the Stehfest method, which has been revealed to tackle such inverse problems accurately, is applied to invert the Laplace domain solution numerically. In this study, we have developed a MATLAB program to calculate the injection pressure in the time domain by combining the Stehfest method and Equation (26). It is found that accurate and stable results can be obtained by choosing the number of terms of the infinite series n ≥ 18 after performing multiple numerical tests. Thus, in order to reduce the computational cost, n = 18 is used in the following calculations.

3. Results

3.1. Comparison of This Study with the Existing Analytical Solution

To verify the accuracy of the analytical solution in this study (Equation (26)), the analytical solution for injection pressure is validated by comparing it with the analytical solution proposed by Ma et al. [46]. These two analytical solutions (Equation (26) in this study and Ma et al. [46]) for water recharge of a water source heat pump are employed to calculate injection pressure with parameters as follows: Q = 40 m3/h, K = 1.25 m/h, d = 60 m, l = 27 m, rw = 0.1 m, and S = 0.007 m−1. The values of the aforementioned parameters are referenced from the work of Ma et al. [46]. As shown in Figure 3, the results of injection pressure calculated using the analytical solution in this study are consistent with the results presented in the previous study (e.g., Ma et al. [46]). This illustrates that the analytical solution of injection pressure proposed in this work is reliable and can be used to analyze the variations of injection pressure under different conditions for geothermal energy extraction systems in an abandoned coal mine.

3.2. Injection Pressure versus Time with the Pumping Rate

Figure 4 presents the transient injection pressure distributions for different pumping rates Q that range from 40 to 70 m3/h. Other parameters in this scenario are given as follows: K = 1.25 m/h, d = 10 m, l = 20 m, rw = 0.1 m, and S = 0.007 m−1. As depicted in Figure 4, the changes in injection pressure are significant when the pumping rate is reduced from 70 to 40 m3/h. Larger pumping rate values result in greater injection pressure. This can be explained by the fact that, given the equivalent aquifer in goaf and the length of pumping screen, which are constant, operations at higher pumping rates will inherently possess a correspondingly higher pumping rate for every unit length of its well screen. This indicates that it needs more power to transmit water for larger pumping rates. Thus, one can see in Figure 4 that greater injection pressures are obtained for larger pumping rates.

3.3. Injection Pressure versus Time with the Well Distance

Figure 5 illustrates the transient injection pressure distributions for different well distances between the production and injection wells from 15 to 30 m. In addition, other parameters used to generate these curves include Q = 40 m3/h, K = 1.25 m/h, d = 10 m, rw = 0.1 m, and S = 0.007 m−1. It can be found in Figure 5 that injection pressure curves merge for different well distances at initial stages, but a smaller well distance leads to a reduced injection pressure at later times. This behavior is due to the fact that the well interference between the production and injection well is not obvious at early times, while it gradually occurs with the increase in time. Also, shorter well distances will cause stronger well interference and have better hydraulic connection. This means that shorter well distances are beneficial for recharging water. Consequently, a smaller injection pressure can be found at late times for a shorter well distance.

3.4. Injection Pressure versus Time with the Specific Storage

Furthermore, it is also important to analyze the impact of specific storage on injection pressure, as indicated in Figure 6. The parameters in this scenario are as follows: Q = 40 m3/h, K = 1.25 m/h, d = 10 m; l = 20 m, rw = 0.1 m, S = 0.001 m−1, 0.002 m−1, 0.005 m−1, and 0.01 m−1. It can be found in Figure 6 that a larger value of specific storage leads to a reduced injection pressure at the beginning, while at later times, all injection pressure curves for different specific storages approach the same value. This is also understandable because an equivalent aquifer with larger specific storage has a greater capacity to release water compared with one with a smaller storage coefficient. Meanwhile, since the process of water release from the equivalent aquifer nears completion at late times (t > 1 h), the curves of injection pressure approach the same asymptotic value for different specific storages, indicating that the effect of specific storage on injection pressure can be neglected.

3.5. Injection Pressure versus Time with the Hydraulic Conductivity

The transient injection pressure distributions for different hydraulic conductivity K are depicted in Figure 7. In this case, other relevant parameters remain unchanged as before: Q = 40 m3/h, d = 10 m, l = 20 m, rw = 0.1 m, and S = 0.007 m−1. The results shown in Figure 7 indicate that changes in injection pressure are remarkable when the hydraulic conductivity K is increased from 0.01 to 0.1 m/h. One can observe that a larger value of hydraulic conductivity leads to a reduced injection pressure. As a consequence, it can be seen that a weaker injection pressure correlated to larger hydraulic conductivity, as illustrated in Figure 7.

3.6. Injection Pressure versus Time with the Thickness of Equivalent Aquifer in Goaf

The distribution curves of transient injection pressure for different thicknesses of equivalent aquifer d from 4 to 16 m are revealed in Figure 8. The values of base parameters chosen for this scenario are as follows: Q = 40 m3/h; K = 1.25 m/h; l = 20 m; rw = 0.1 m; and S = 0.007 m−1. One can observe from Figure 8 that the injection pressure for the abandoned coal mine water source heat pump system is perceptibly influenced by the thickness of the equivalent aquifer in goaf. The thickness of the equivalent aquifer with a larger value induces a smaller injection pressure. Additionally, it is noted that the changes in injection pressure with smaller thickness of equivalent aquifer (from 4 to 8 m) are notably larger than those for the larger thickness of equivalent aquifer, which is increased from 8 to 16 m. This occurs because the fully penetrating well with a larger thickness of equivalent in goaf has a correspondingly larger well screen, indicating that the pumping rate for a unit length of the well screen is larger for a smaller thickness of equivalent aquifer while the pumping rate remains constant. Also, differences in injection pressure are smaller for larger thicknesses of equivalent aquifer (e.g., from 8 to 16 m) than for smaller thicknesses of equivalent aquifer, as demonstrated in Figure 8.

3.7. Sensitivity Analysis

Sensitivity analyses are conducted to observe the effect of various parameters, including Q, l, S, and K, on the injection pressure for the geothermal energy extraction system in an abandoned coal mine. In this study, the widespread normalized parameter sensitivity analysis proposed by Kabala [47] is implemented, and the coefficient of normalized sensitivity can take the following form:
Y m , n = R n X m R n
where Ym,n represents the normalized sensitivity due to variations in the m th parameter Rn at the nth time and Xm refers to a function of Rn and n. A finite-difference equation is employed to approximate the partial derivative term on the right side of Equation (27) [48]:
X m R n = X m R n + Δ R n X m R n Δ R n
where the small increment Δ Rn is normally taken as 0.01Rn, similarly as in Yeh [48], Wen et al. [49], and Tu et al. [50,51,52]. In addition, the following parameters are employed in this scenario to perform the sensitivity analysis: Q = 40 m3/h, K = 1.25 m/h, l = 20 m, d = 10 m, rw = 0.1 m, and S = 0.007 m−1.
Figure 9 presents temporal distribution curves of the normalized sensitivity for selected parameters. The sensitivity analysis reveals that two parameters, the hydraulic conductivity K and the pumping rate Q, have notable impacts on the injection pressure, but the injection pressure is most sensitive to the hydraulic conductivity, as shown in Figure 9. It is interesting to note that the hydraulic conductivity K shows a negative impact on injection pressure, and this effect approaches a steady state at late times (t > 1 h). Also, the overall behavior of the normalized sensitivity for the pumping rate Q is obviously similar to that of the hydraulic conductivity K, except that the pumping rate Q has a positive impact on injection pressure. Additionally, as illustrated in Figure 9, the influences of the well distance l and the specific storage S on the injection pressure are minimal and can be neglected.

4. Application of the Proposed Model

The application of the proposed analytical model in this study provides a theoretical base for further analyzing the injection pressure of an open-loop geothermal energy extraction system in an abandoned coal mine, considering the effects of non-Darian, well skin, and different aquifer systems. Additionally, the presented analytical solution would be of significance in further estimating the hydraulic parameters of such a system in an abandoned coal mine.
In practical engineering applications, many factors such as well configuration, hydrogeological conditions, thermal geological conditions, costs, and energy consumption need to be considered when designing an open-loop geothermal energy extraction system in an abandoned coal mine. As revealed in Figure 9, the pumping rate in real engineering applications needs to be adjusted according to the efficiency of the geothermal energy extraction system. As indicated in Figure 4, larger pumping rates result in greater injection pressure, which means that, although larger pumping rates could make the system extract more energy, the energy consumption and costs for injection will increase. Therefore, factors such as cost, energy consumption, and efficiency should be comprehensively considered while designing a proper pumping rate for a geothermal extraction system in an abandoned coal mine. Similarly, the well distance is also an important factor to optimize the efficiency of such a geothermal energy extraction system in an abandoned coal mine. It can be seen in Figure 5 that the injection pressure for a smaller well distance is relatively small compared with that for a larger well distance. This suggests that the well interference for a smaller well distance is more significant, which can be very beneficial for water injection, but it will cause thermal breakthrough more quickly and decrease the efficiency of this system. Given that hydraulic conductivity can remarkably affect the injection pressure, it is necessary to investigate hydrological conditions of an abandoned coal mine, especially related to hydraulic conductivity, by employing pumping tests, geophysical surveys, etc.

5. Conclusions

In this study, an analytical model for characterizing mine water recharge for geothermal energy extraction system in an abandoned coal mine has been proposed. The semi-analytical solution for injection pressure is developed leveraging a Laplace transform. With the presented model, the effects of different parameters on injection pressure are analyzed. The specific conclusions of this study can be drawn are as below:
(1)
The injection pressure varies significantly with pumping rate, hydraulic conductivity, and the thickness of equivalent aquifer. A larger pumping rate results in a greater injection pressure, while a larger value of the hydraulic conductivity and the thickness of equivalent aquifer led to a reduced injection pressure.
(2)
The well distance has little impact on the early-stage injection pressure, although it is significant for the late-stage injection pressure. In contrast, larger values of specific storage result in lessened injection pressure at the beginning, but the effects of specific storage on injection pressure can be neglected at a late stage.
(3)
The limitations of the proposed analytical model in this study could also be clarified. The analytical model for injection pressure in an open-loop geothermal energy extraction system in an abandoned coal mine is a relatively idealized model, which do not consider the well skin, inhomogeneity of equivalent aquifer, non-Darcian flow, and different aquifer systems. Additionally, it is noteworthy that although the proposed analytical solution in this study is verified with the existing analytical solution, it would be better in the future to use actual data fromthe field or laboratory for further investigation. In future research, it is necessary to establish a new analytical model considering the effects of non-Darcian flow in such a system. This is more in line with the field reality and will help in improving our understanding of the complex water flow regime in an open-loop geothermal energy extraction system in an abandoned coal mine.

Author Contributions

K.T.: project supervision, methodology, original draft writing, data curation, formal analysis. X.P.: data curation, drawing, writing, formal analysis. H.Z. (Hongwei Zhang): project supervision, writing—review and editing, data curation. X.L.: supervision, formal analysis. H.Z. (Hongyi Zhao): validation, investigation. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (Grant No. 42302278), the China Postdoctoral Science Foundation (Grant No. 2024T171007), the Major Research Instrument Development Project of the National Natural Science Foundation of China (Grant No. 42027801), the Key Laboratory of Shallow Geothermal Energy, Ministry of Natural Resources of the People’s Republic of China (Grant No. KLSGE202401-02), the Beijing Natural Science Foundation (Grant No. 3232026), the Ordos City Landmark Team Project (2022).

Data Availability Statement

The datasets and computer codes used in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ling, Z.; Huang, T.; Li, J.; Zhou, S.; Lian, L.; Wang, J.; Zhao, Y.; Mao, X.; Gao, H.; Ma, J. Sulfur dioxide pollution and energy justice in Northwestern China embodied in West-East Energy Transmission of China. Appl. Energy 2019, 238, 547–560. [Google Scholar] [CrossRef]
  2. Xie, H.; Ren, S.; Xie, Y.; Jiao, X. Development opportunities of the coal industry towards the goal of carbon neutrality. J. China Coal Soc. 2021, 46, 2197–2211. (In Chinese) [Google Scholar] [CrossRef]
  3. Li, N.; Chen, W. Energy-water nexus in China’s energy bases: From the Paris agreement to the Well Below 2 Degrees target. Energy 2019, 166, 277–286. [Google Scholar] [CrossRef]
  4. Pu, H.; Bian, Z.; Zhang, J.; Xu, J. Research on a reuse mode of geothermal resources in abandoned coal mines. J. China Coal Soc. 2021, 46, 677–687. [Google Scholar] [CrossRef]
  5. Pan, X.; Wang, H.; Wang, L.; Chen, W. Decarbonization of China’s transportation sector: In light of national mitigation toward the Paris Agreement goals. Energy 2018, 155, 853–864. [Google Scholar] [CrossRef]
  6. Bórawski, P.; Bełdycka-Bórawska, A.; Szymańska, E.J.; Jankowski, K.J.; Dubis, B.; Dunn, J.W. Development of renewable energy sources market and biofuels in The European Union. J. Clean. Prod. 2019, 228, 467–484. [Google Scholar] [CrossRef]
  7. Chen, J.; Shi, Q.; Shen, L.; Huang, Y.; Wu, Y. What makes the difference in construction carbon emissions between China and USA? Sustain. Cities Soc. 2019, 44, 604–613. [Google Scholar] [CrossRef]
  8. Fan, J.; Xie, H.; Jie, C.; Jiang, D.; Li, C.; Tiedeu, W.N.; Ambre, J. Preliminary feasibility analysis of a hybrid pumped-hydro energy storage system using abandoned coal mine goafs. Appl. Energy 2020, 258, 114007. [Google Scholar] [CrossRef]
  9. Yang, X.; Wan, H.; Zhang, Q.; Zhou, J.-C.; Chen, S.-Y. A scenario analysis of oil and gas consumption in China to 2030 considering the peak CO2 emission constraint. Pet. Sci. 2016, 13, 370–383. [Google Scholar] [CrossRef]
  10. Yuan, L.; Jiang, Y.; Wang, K.; Zhao, Y.; Hao, X.; Xu, C. Precision exploitation and utilization of closed/abandoned mine resources in China. J. China Coal Soc. 2018, 43, 14–20. (In Chinese) [Google Scholar] [CrossRef]
  11. Lyu, X.; Zhang, T.; Yuan, L.; Fang, J. Prospects for the transformation and development of carbon storage in abandoned mines of coal enterprises from the perspective of carbon neutrality. Int. J. Coal Sci. Technol. 2023, 10, 36. [Google Scholar] [CrossRef]
  12. Liu, F.; Cao, W.; Zhang, J.; Cao, G.; Guo, L. Current technological innovation and development direction of the 14th Five-Year Plan period in China coal industry. J. China Coal Soc. 2021, 46, 1–15. (In Chinese) [Google Scholar] [CrossRef]
  13. Ge, Z.; Deng, K.; Zhang, L.; Zuo, S. Development potential evaluation of CO2-ECBM in abandoned coal mines. Greenh. Gases Sci. Technol. 2020, 10, 643–658. [Google Scholar] [CrossRef]
  14. Watzlaf, G.R.; Ackman, T.E. Underground Mine Water for Heating and Cooling using Geothermal Heat Pump Systems. Mine Water Environ. 2006, 25, 1–14. [Google Scholar] [CrossRef]
  15. Ramos, E.P.; Breede, K.; Falcone, G. Geothermal heat recovery from abandoned mines: A systematic review of projects implemented worldwide and a methodology for screening new projects. Environ. Earth Sci. 2015, 73, 6783–6795. [Google Scholar] [CrossRef]
  16. Guo, P.; Zheng, L.; Sun, X.; He, M.; Wang, Y.; Shang, J. Sustainability evaluation model of geothermal resources in abandoned coal mine. Appl. Therm. Eng. 2018, 144, 804–811. [Google Scholar] [CrossRef]
  17. Al-Habaibeh, A.; Athresh, A.P.; Parker, K. Performance analysis of using mine water from an abandoned coal mine for heating of buildings using an open loop based single shaft GSHP system. Appl. Energy 2018, 211, 393–402. [Google Scholar] [CrossRef]
  18. Menendez, J.; Ordonez, A.; Fernandez-Oro, J.M.; Loredo, J.; Diaz-Aguado, M.B. Feasibility analysis of using mine water from abandoned coal mines in Spain for heating and cooling of buildings. Renew. Energy 2020, 146, 1166–1176. [Google Scholar] [CrossRef]
  19. Guo, P.; Wang, M.; Dang, G.; Zhu, T.; Wang, J.; He, M. Evaluation method of underground water storage space and thermal reservoir model in abandoned mine. Rock Mech. Bull. 2023, 2, 100044. [Google Scholar] [CrossRef]
  20. Hall, A.; Scott, J.A.; Shang, H. Geothermal energy recovery from underground mines. Renew. Sustain. Energy Rev. 2011, 15, 916–924. [Google Scholar] [CrossRef]
  21. Ramos, E.P.; Falcone, G. Recovery of the geothermal energy stored in abandoned mines. In Clean Energy Systems in the Subsurface: Production, Storage and Conversion; Springer: Berlin/Heidelberg, Germany, 2013; pp. 143–155. [Google Scholar] [CrossRef]
  22. Xie, H.; Zhao, J.W.; Zhou, H.W.; Ren, S.H.; Zhang, R.X. Secondary utilizations and perspectives of mined underground space. Tunn. Undergr. Space Technol. 2020, 96, 103129. [Google Scholar] [CrossRef]
  23. Jessop, A.M.; Macdonald, J.K. Clean energy from abandoned mines at Springhill, Nova Scotia. Energy Sources 1995, 17, 93–106. [Google Scholar] [CrossRef]
  24. Kranz, K.; Dillenardt, J. Mine Water Utilization for Geothermal Purposes in Freiberg, Germany: Determination of Hydrogeological and Thermophysical Rock Parameters. Mine Water Environ. 2010, 29, 68–76. [Google Scholar] [CrossRef]
  25. Qi, T.; Zhang, F.; Pei, X.; Feng, G.; Wei, H. Simulation research and application on response characteristics of detecting water-filled goaf by transient electromagnetic method. Int. J. Coal Sci. Technol. 2022, 9, 17. [Google Scholar] [CrossRef]
  26. Neves, R.; Cho, H.; Zhang, J. Techno-economic analysis of geothermal system in residential building in Memphis, Tennessee. J. Build. Eng. 2020, 27, 100993. [Google Scholar] [CrossRef]
  27. Guo, P.; Wang, M.; Sun, X.; He, M. Study on off-season cyclic energy storage in underground space of abandoned space. J. China Coal Soc. 2022, 47, 2193–2206. [Google Scholar] [CrossRef]
  28. Zhang, C.; Luo, B.; Xu, Z.; Sun, Y.; Feng, L. Research on the Capacity of Underground Reservoirs in Coal Mines to Protect the Groundwater Resources: A Case of Zhangshuanglou Coal Mine in Xuzhou, China. Water 2023, 15, 1468. [Google Scholar] [CrossRef]
  29. Ghoreishi-Madiseh, S.A.; Hassani, F.; Abbasy, F. Numerical and experimental study of geothermal heat extraction from backfilled mine stopes. Appl. Therm. Eng. 2015, 90, 1119–1130. [Google Scholar] [CrossRef]
  30. Chudy, K. Mine Water as Geothermal Resource in Nowa Ruda Region (SW Poland). Water 2022, 14, 136. [Google Scholar] [CrossRef]
  31. Wang, H.; Xu, Y.; Yuan, L.; Sun, Y.; Cai, Y. Analysis of geothermal heat recovery from abandoned coal mine water for clean heating and cooling: A case from Shandong, China. Renew. Energy 2024, 228, 120659. [Google Scholar] [CrossRef]
  32. Jardon, S.; Ordonez, A.; Alvarez, R.; Cienfuegos, P.; Loredo, J. Mine Water for Energy and Water Supply in the Central Coal Basin of Asturias (Spain). Mine Water Environ. 2013, 32, 139–151. [Google Scholar] [CrossRef]
  33. Bao, T.; Liu, Z. Geothermal energy from flooded mines: Modeling of transient energy recovery with thermohaline stratification. Energy Convers. Manag. 2019, 199, 111956. [Google Scholar] [CrossRef]
  34. Rodríguez, R.; Díaz, M.B. Analysis of the utilization of mine galleries as geothermal heat exchangers by means a semi-empirical prediction method. Renew. Energy 2009, 34, 1716–1725. [Google Scholar] [CrossRef]
  35. Hamm, V.; Sabet, B.B. Modelling of fluid flow and heat transfer to assess the geothermal potential of a flooded coal mine in Lorraine, France. Geothermics 2010, 39, 177–186. [Google Scholar] [CrossRef]
  36. Raymond, J.; Therrien, R. Optimizing the design of a geothermal district heating and cooling system located at a flooded mine in Canada. Hydrogeol. J. 2014, 22, 217–231. [Google Scholar] [CrossRef]
  37. Sun, Y.; Zhang, X.; Li, X.; Duan, C. A geothermal energy heat exchange system suitable for abandoned mines and its optimization. Heat Mass Transf. 2023, 59, 1749–1766. [Google Scholar] [CrossRef]
  38. Zhang, Z.; Zu, W.; Zhang, W.; Wang, K.; Ma, X.; Cui, P. Investigation of Theoretical Models, Pumping-Recharge Well Arrangements and System Performance of Abandoned Mine Water Source Heat Pump. Energy Built Environ. 2024, in press. [Google Scholar] [CrossRef]
  39. Şen, Z. Type Curves for Two-Regime Well Flow. J. Hydrol. Eng. 1988, 114, 1431–1542. [Google Scholar] [CrossRef]
  40. Papadopulos, I.S.; Cooper, H.H., Jr. Drawdown in a well of large diameter. Water Resour. Res. 1967, 3, 241–244. [Google Scholar] [CrossRef]
  41. Liu, G.; Wu, S.; Fan, Z.; Zhou, Z.; Xie, C.; Wu, J.; Liu, Y. Analytical derivation on recharge and periodic backwashing process and the variation or recharge pressure. J. Jilin Univ. (Earth Sci. Ed.) 2016, 46, 1799–1807. [Google Scholar] [CrossRef]
  42. Stehfest, H. Algorithm 368: Numerical inversion of Laplace transforms. Commun. ACM 1970, 13, 47–49. [Google Scholar] [CrossRef]
  43. Stehfest, H. Remark on Algorithm 368: Numerical inversion of Laplace transforms. Commun. ACM 1970, 13, 624. [Google Scholar] [CrossRef]
  44. Crump, K.S. Numerical inversion of Laplace transforms using a Fourier-series approximation. J. ACM 1976, 23, 89–96. [Google Scholar] [CrossRef]
  45. DeHoog, F.R.; Knight, J.H.; Stokes, A.N. An improved method for numerical inversion of Laplace transforms. SIAM J. Sci. Stat. Comput. 1982, 3, 357–366. [Google Scholar] [CrossRef]
  46. Ma, C.; Zhou, W.; Yun, T.; Liu, X. Study of recharge ability for groundwater heat pump with pumping and recharging in two wells. Acta Energiae Solaris Sin. 2015, 36, 593–598. [Google Scholar] [CrossRef]
  47. Kabala, Z.J. Sensitivity analysis of a pumping test on a well with wellbore storage and skin. Adv. Water Resour. 2001, 24, 483–504. [Google Scholar] [CrossRef]
  48. Yeh, W. Review of parameter identification procedures in groundwater hydrology: The inverse problem. Water Resour. Res. 1986, 22, 95–108. [Google Scholar] [CrossRef]
  49. Wen, Z.; Liu, K.; Chen, X. Approximate analytical solution for non-Darcian flow toward a partially penetrating well in a confined aquifer. J. Hydrol. 2013, 498, 124–131. [Google Scholar] [CrossRef]
  50. Tu, K.; Wu, Q.; Simunek, J.; Chen, C.; Zhu, K.; Zeng, Y.; Xu, S.; Wang, Y. An Analytical Solution of Groundwater Flow in a Confined Aquifer with a Single-Well Circulation System. Water Resour. Res. 2020, 56, e2020WR027529. [Google Scholar] [CrossRef]
  51. Tu, K.; Wu, Q.; Simunek, J.; Zhu, K.; Chen, C.; Zheng, W.; Zeng, Y.; Xu, S. An approximate analytical solution for non-Darcian flow in a confined aquifer with a single well circulation groundwater heat pump system. Adv. Water Resour. 2020, 145, 103740. [Google Scholar] [CrossRef]
  52. Tu, K.; Wu, Q.; Zhang, N.; Liu, X. An analytical model for a single-well circulation system in a confined aquifer with a finite thickness well skin. International J. Numer. Anal. Methods Geomech. 2023, 47, 1289–1309. [Google Scholar] [CrossRef]
Figure 1. An overview of geothermal energy extraction with heat pump technology from an abandoned coal mine.
Figure 1. An overview of geothermal energy extraction with heat pump technology from an abandoned coal mine.
Water 16 02781 g001
Figure 2. The schematic showing a geothermal energy extraction system in an abandoned coal mine.
Figure 2. The schematic showing a geothermal energy extraction system in an abandoned coal mine.
Water 16 02781 g002
Figure 3. Comparison of the analytical solution in this study with the results of Ma et al. [46].
Figure 3. Comparison of the analytical solution in this study with the results of Ma et al. [46].
Water 16 02781 g003
Figure 4. The variations of injection pressure for different pumping rates.
Figure 4. The variations of injection pressure for different pumping rates.
Water 16 02781 g004
Figure 5. The variations of injection pressure for different well distances.
Figure 5. The variations of injection pressure for different well distances.
Water 16 02781 g005
Figure 6. The variations of injection pressure for different specific storages.
Figure 6. The variations of injection pressure for different specific storages.
Water 16 02781 g006
Figure 7. The variations of injection pressure for different hydraulic conductivities.
Figure 7. The variations of injection pressure for different hydraulic conductivities.
Water 16 02781 g007
Figure 8. The variations of injection pressure for different thicknesses of equivalent aquifer in goaf.
Figure 8. The variations of injection pressure for different thicknesses of equivalent aquifer in goaf.
Water 16 02781 g008
Figure 9. The normalized sensitivity of injection pressure to selected parameters.
Figure 9. The normalized sensitivity of injection pressure to selected parameters.
Water 16 02781 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tu, K.; Pan, X.; Zhang, H.; Li, X.; Zhao, H. An Analytical Solution for Characterizing Mine Water Recharge of Water Source Heat Pump in Abandoned Coal Mines. Water 2024, 16, 2781. https://doi.org/10.3390/w16192781

AMA Style

Tu K, Pan X, Zhang H, Li X, Zhao H. An Analytical Solution for Characterizing Mine Water Recharge of Water Source Heat Pump in Abandoned Coal Mines. Water. 2024; 16(19):2781. https://doi.org/10.3390/w16192781

Chicago/Turabian Style

Tu, Kun, Xiaoqiang Pan, Hongwei Zhang, Xiang Li, and Hongyi Zhao. 2024. "An Analytical Solution for Characterizing Mine Water Recharge of Water Source Heat Pump in Abandoned Coal Mines" Water 16, no. 19: 2781. https://doi.org/10.3390/w16192781

APA Style

Tu, K., Pan, X., Zhang, H., Li, X., & Zhao, H. (2024). An Analytical Solution for Characterizing Mine Water Recharge of Water Source Heat Pump in Abandoned Coal Mines. Water, 16(19), 2781. https://doi.org/10.3390/w16192781

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop