Next Article in Journal
Assessing Potential of Dor River as Small Hydro Project for Lessening Energy Crisis and Enhancing Tarbela Reservoir Life in Khyber Pakhtunkhwa, Pakistan
Next Article in Special Issue
Analytical Method for Groundwater Seepage through and Beneath a Fully Penetrating Cut-off Wall Considering Effects of Wall Permeability and Thickness
Previous Article in Journal
Experimental Investigation on Bragg Resonant Reflection of Waves by Porous Submerged Breakwaters on a Horizontal Seabed
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Groundwater Flow Changes Associated with Drainage within Multilayer Aquifers in a Semiarid Area

1
State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Nanjing Hydraulic Research Institute, Nanjing 210029, China
2
Research Center for Climate Change, Ministry of Water Resources, Nanjing 210029, China
3
Hetaoyu Coal Mine, Huaneng Qingyang Coal and Power Co., Ltd., Qingyang 745000, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(17), 2679; https://doi.org/10.3390/w14172679
Submission received: 28 June 2022 / Revised: 23 August 2022 / Accepted: 26 August 2022 / Published: 29 August 2022
(This article belongs to the Special Issue Flow and Transport Processes in Groundwater Systems)

Abstract

:
In order to evaluate the impact of groundwater drainage on groundwater flow, the Hetaoyu coal field was taken as a case study in the Longdong area, China, where the coal seam was covered with multilayer aquifers. A three-dimensional unsteady groundwater flow model and a one-dimensional fracture water flow model were calculated by joint equations for changing hydrogeological structures under coal mining. According to the results, mine construction had greatly affected groundwater reserves in the Quaternary phreatic aquifer, Cretaceous Huanhe confined aquifer, and Luohe confined aquifer. The groundwater drainage was mainly from the Cretaceous aquifer, in which the aquifer reserves of the Luohe Formation decreased by 30,861.8 m3/m, accounting for about 92% of the total changes in local groundwater reserves. A drop funnel with an area of about 2.3 km2 would be formed under the groundwater discharge of 187.6 m3/h for the main inclined shaft excavation of the Hetaoyu coal mine. With the continuation of mining activities, the mine water flow will reach 806.83 m3/h and would result in descending funnel area of about 4.5 km2, the groundwater level drawdown at least 16 m, which would exceed the limited value regulated by the government. Therefore, in order to ensure the safety of coal mining and protect groundwater resources, the Hetaoyu Coal Mine departments should take some water loss prevention and control projects to reduce the drawdown of groundwater.

1. Introduction

The groundwater resource is critical for human survival and economic development in the semiarid area of Northwest China, where the coal seam and oil-bearing formation are usually overlaid by multilayer aquifers [1]. However, for coal and oil safety production, a large amount of groundwater should be discharged to the surface [2], which would influence local groundwater storage [3]. In order to protect the groundwater resource, the groundwater flow was usually simulated before coal mining [4]. Owing to the complexity of hydrogeological conditions in the mining area, the hydrogeological parameters vary strongly with time and space [5], and even the geological structure could change with mining [6]. So, it was difficult to estimate groundwater flow for each aquifer accurately.
There were several methods to estimate groundwater flow, such as the numerical model, analytical method, water balance method, and hydrochemistry method [7].
The numerical model is based on observed geological data [8] and can adjust the parameters conveniently when hydrogeological conditions change. Because of the groundwater drainage, the balance was disturbed, which resulted in the hydraulic connection changes between previously separate aquifers [9]. Based on the traditional groundwater numerical model, some scholars attempted to adjust parameters to solve the discontinuity in the groundwater model when the new water-connected zone formed. Sun et al. (1989) first introduced two correction factors to the hydrology conductivity model when groundwater went through a fault [10]. Liang et al. (2012) encrypted the numerical model grids for the fault and then assigned different values to each node [11]. Based on the previous study, Karimzade et al. (2017) established the probability distribution functions for the fracture size, fracture orientation, and fracture density, respectively. According to these functions, flow paths in the rock were identified [12]. Mu et al. (2020) established a numerical couple model between fluid and solid matter to simulate the fault-zone deformation and water inflow under natural and disturbed conditions [13]. Song and Liang (2021) applied the numerical simulation method to analyze the process of fracture activation driven by mining, and the results showed good agreement between the theoretical calculation and the numerical simulation [14]. Xue et al. (2022) established a coupling model of gas drainage and groundwater loss to analyze the groundwater loss [15]. Liu et al. (2020) used the mixing calculations to identify sources of groundwater recharge in a Coastal Gold Mine [16]. Therefore, the numerical model has been widely used in groundwater inflow simulation, but there is much less research on coupling the 3D groundwater flow and 1D fracture water flow in a mining area.
Compared with the numerical model, the water balance method is based on the principle of mass conservation, but it requires more hypotheses for groundwater prediction [17]. Based on groundwater dynamics in a certain equilibrium period, the relationship between charge and discharge items was established, so the equilibrium equation can be used to predict the mine water inflow, which is suitable for the water-filled deposit with nonseepage groundwater inflow and simple water equilibrium conditions [18]. Elmarami et al. (2017) pointed out that the water balance method was suitable for large-scale calculation. The smaller the scale, the greater the uncertainty of each equilibrium term [19]. The groundwater flow calculated by this method can only be used as a reference value, not as an accurate result. In addition, climate change would have a significant effect on groundwater recharge [20]. If the recharge amount is smaller than other equilibrium terms, it will be very greatly influenced by evaporation. So, the water balance method is not suitable for fine calculation in the mining area with a multilayer aquifer.
On the other hand, the hydrochemistry method was useful for identifying the source of groundwater drainage by hydrochemical tracer test [21]. Qiao et al. (2019) took 76 groundwater samples for hydrogeochemical analysis, which indicated that coal mining caused the rock strata to sink, obstructing the water connection between aquifers [22]. Duan et al. (2019) identified the source of the coastal mine water inrush using hydrogeochemical analysis and inferred the potential water flow channels and the prevailing supply patterns [23]. Pilla et al. (2006) adopted hydrochemistry and isotope geochemistry as tools to define the recharge mechanisms, the origin of groundwater, and the hydraulic confinement of deep aquifers [24]. Powell et al. (2003) assessed the microbial quality of groundwater within the Permo-Triassic Sherwood Sandstone aquifers underlying the cities of Nottingham and Birmingham, and the groundwater sources were detected by microbial quality [25]. However, the hydrochemical method could not reveal the mechanism of groundwater flow, and it is difficult to couple the model with groundwater flow and fracture water flow in changing environment.
Therefore, the numerical model is closely associated with hydrogeological conditions, which could simulate the groundwater flow more accurately. So, in this paper, the coupled model of three-dimensional unsteady groundwater flow and one-dimensional fracture water flow was established based on hydrogeological borehole data, groundwater level, and so on. The hydrogeological structure model would be modified when the new water-conducted zone formed. The results would provide some guidance on the utilization and protection of local groundwater resources.

2. Materials and Methods

2.1. Study Area

The Hetaoyu coal field is located in Zhengning county in the southeastern city of Qiangyang, Gansu Province. It belongs to the Longdong Loess Plateau, and the landform is mainly composed of the Loess plateau area, Loess ridge, and valley area. As shown in Figure 1, the altitude is generally 889 to 1310 m. The plateau height decreases gradually from 1300 m to 1170 m with higher elevation in the north and east. The valley area mainly consists of Jinghe River and its tributaries, including Wuritiangou River and Silang River. Jinghe River Valley is about 210 m lower than the Loess Plateau.
It belongs to a semihumid and semiarid climate, with an annual average precipitation of 598.44 mm, annual average evaporation of 1496.66 mm, and average relative humidity of 56–75%. Precipitation is mainly concentrated from July to September. The frost-free period is 101–117 days, and the maximum permafrost depth is 89 cm.
There are three exploitable coal seams in this coalfield, named coal-2 layer, coal-5 layer, and coal-8 layer. The design production capacity of the coal mine is 12.0 Mt/a. The main inclined shaft for coal mining is along the southwest direction, and the vertical angle is 7°. The Quaternary aquifer, Huanhe Confined aquifer, and Luohe Confined aquifer would be exposed in sequence, and then a large amount of the groundwater would be discharged to the surface. So, the groundwater around the main inclined shaft would be affected seriously, but further analysis of its impact is necessary for evaluating groundwater flow changes.
Therefore, in order to analyze the effect of the main inclined shaft excavation on the groundwater and the groundwater dynamics under coal-8 mining in the future, the study area took the main inclined shaft as the center; other boundaries were parallel or vertical to the equal-head line of the confined aquifer of Luohe Formation, which is situated at 107°59′2″–108°04′32″ E and 35°15′59″–35°20′3″ N, with an area of 32.3 km2 within the Hetaoyu coalfield. There were four geological boreholes and three hydrological holes in this area, named #1, #2, and #3, as shown in Figure 2. The #1 recorded the groundwater level in Luohe Confined aquifer from November 2019 to May 2021. The #3 recorded the groundwater level in Huanhe Confined aquifer from January 2019 to November 2019. The #2 recorded the mixed groundwater level of the Quaternary aquifer, Huanhe Confined aquifer, and Luohe Confined aquifer from November 2019 to May 2021. The isolevel map of the Luohe Formation aquifer is shown in Figure 2.
The vertical geological structure was divided into seven layers, including three aquifers, three aquicludes, and a coal seam, by integrating similar and thin lithologies. The aquifers are Quaternary phreatic aquifers, Cretaceous Huanhe Formation confined aquifers, and Cretaceous Luohe Formation confined aquifers. The hydrogeological section is shown in Figure 3. Because of the discontinuous distribution of aquifers in the Huanhe Formation and Luohe Formation, there is a need to consolidate the aquifers within the same layer for convenient calculation. The Quaternary phreatic aquifers are composed of loess, and confined aquifers in the Cretaceous are composed of siltstone, fine sandstone, coarse sandstone, medium sandstone, and medium conglomerate. The aquicludes are mainly silty mudstone, mudstone, and clay.

2.2. Model Description

There are several methods to estimate groundwater dynamics, such as the numerical model, analytical method, water balance method, mathematical, statistical method, and so on. However, the numerical model can not only describe the aquifer system and boundary conditions more realistically but also couple with different models more flexibly by writing computational codes. So, a three-dimensional groundwater dynamics model was adopted in this paper. Using a Cartesian coordinate system, the general equation for water flow with saturated, heterogeneous, and anisotropic porous media is expressed as follows
x ( K x H x ) + y ( K y H y ) + z ( K z H z ) + W = S s H t
where H is the total head, W is the recharge rate, Ss is the specific storage of the porous medium, Kx, Ky, Kz are the saturated hydraulic conductivity in x , y , and z directions, respectively.
According to the pumping test results from the Huanhe Formation confined aquifers and Luohe Formation confined aquifers, the initial values of hydrological parameters are listed in Table 1.
A one-dimensional fracture water flow model under the dominant water-conducted direction is usually used to represent the flow within the water-conducted zone. The numerical model is as follows, based on Navier–Stokes Equation
q f = γ 4 μ J f 0 b / 2 ( b 2 4 y 2 ) d y = g δ 3 12 μ J f
v f = g δ 2 12 μ J f = K f J f
where q f is the unit discharge within the water-conducted zone, v f is the average flow velocity within the water-conducted zone, g is the groundwater specific gravity, is the width of fracture, μ is the kinetic viscosity coefficient of water, K f is the conductivity of fracture, J f is the hydraulic gradient of flow.
The permeability changes are associated with mining-induced stress and water pressure. It is usually a negative relationship between hydraulic conductivity and the principal stress of fracture zones described by a negative exponential function [26]. It is a positive relationship between hydraulic conductivity and pore water pressure of fracture zones. The relationship could be expressed as follow [13]
K f = K 0 e α ( σ m p )
where σ m is the mean principal stress, α is a fitting parameter, p is the pore water pressure, and K 0 is the hydraulic conductivity.
The maximum height of the water-conducted fracture zone can be calculated by Equation (5)
h c o n d = m n = 1 k 1 h n ( γ n 1 )
where h c o n d is the maximum height of water conducted zone, m is the thickness of the coal seam, n is the number of rock layers, k is the total number of rock layers, h n is rock thickness at n layer, and γ n is the coefficient of rock-breaking expansion.
Equation (1) can be discretized by Equation (6) with a finite difference scheme for isotropic porous media. The calculation nodes are shown in Figure 4.
K x Δ x 2 ( H i , j 1 , k n 2 H i , j , k n + H i , j + 1 , k n ) + K y Δ y 2 ( H i , j , k 1 n 2 H i , j , k n + H i , j , k + 1 n ) + K z Δ z 2 ( H i 1 , j , k n 2 H i , j , k n + H i + 1 , j , k n ) + W n = S s Δ t ( H i , j , k n H i , j , k n 1 )
where Kx, Ky, Kz, and Ss are constant in each aquifer, n and n − 1 are the time at n and n − 1, H i , j , k n is the total head of node , i , j , k at n time, and W n is the total recharge rate at time n. In this paper, W is the mine drainage varied with time during main inclined shaft excavation, as shown in Figure 5.

3. Results and Discussion

3.1. Spatial and Temporal Variation of Groundwater Level for Different Aquifers

The variations of groundwater levels for different aquifers are shown in Figure 6a,b. It could be noticed that the trends of groundwater level varied from each other, especially in #3 the groundwater level had been rising with the increase in groundwater drainage. That phenomenon indicated that although there was a large amount of groundwater discharge in the coal mine area, the influence on regional groundwater flow was limited. The distance between the main inclined shaft and #3 was about 3 km, which might be the threshold of the drainage influence.
The Luohe aquifer in #1 presented a downward trend accompanied by groundwater drainage since November 2020, but the groundwater level was relatively steady with time in #2 at about 1 km east of #1. That meant that the influence area of drainage might only reach near #1 and #2.
Therefore, although groundwater drainage caused the groundwater level to drop, the groundwater level at the location about 3 km away from the drainage area was not affected by the drainage. A similar phenomenon would occur in coastal areas, where the groundwater replenishes at a rapid rate. The groundwater level dropped by 29 m to 32 m after 7 days of pumping. After stopping pumping, the groundwater level would recover quickly, so the change in water level in a short time would not lead to large land subsidence and have little impact on the surrounding environment [27].

3.2. The Relationship between Groundwater Levels of Various Aquifers

It is well known that the impact of groundwater drainage is greatest in the mining area and gradually decreases with increasing distance. According to the above analysis, the variation of groundwater level in #2 was lagging behind that of #1. In order to quantify the lag time, the correlation coefficient was calculated between #1 and #2. The result is shown in Figure 7.
As shown in Figure 7, the lag time maximum correlation coefficient appeared at 28d, indicating that under the groundwater drainage, each aquifer did not present a synchronous downward trend; after the first aquifer dropped, the near aquifer began to decline about 28d later. This meant that the groundwater flow was exchanged frequently in this area, and the losses of water would be recharged in a short time. This result also indicated that the groundwater drainage in the mining area had little influence on the regional groundwater flow.

3.3. Influence on Local Groundwater Dynamics under Groundwater Drainage

Based on the above results, taking the groundwater drainage as a discharge item in Equation (6), the groundwater dynamic was analyzed, and the source of groundwater drainage was identified. The groundwater discharge of 187.6 m3/h was used as an example in Figure 5. The driving length of the main inclined shaft was drilled to the aquifer of the Luohe Formation in June 2019, and the impact of drainage on groundwater flow is shown in Figure 8.
From Figure 8, with a large amount of groundwater discharging to the surface, the groundwater level of Luohe Formation dropped by about 14 m, and a descending funnel of about 2.3 km2 was formed in the drainage area, which affected the characteristics of local groundwater flow. However, there was no change in the direction of groundwater flow on the regional scale, where the regional groundwater flow was discharged from the northeastern mountainous area to the Malian River. Nevertheless, a large amount of groundwater discharged into the Malian River was reduced as a result of groundwater drainage.
According to the drop depth of water level, the reserve changes of each aquifer are shown in Figure 9.
It showed that the Quaternary phreatic aquifer reserves and the Cretaceous Huanhe Formation aquifer reserves decreased by 249.3 m3/month and 3270.4 m3/month, respectively. When the driving length of the main inclined shaft was excavated into the Luohe confined aquifer, the groundwater drainage increased sharply, and the Cretaceous Luohe Formation aquifer reserves decreased by 30,861.8 m3/month, which accounted for about 92% of the total reserve reduction. Therefore, the Luohe Formation aquifer was the main source of groundwater drainage. In the process of coal mining, it is necessary to pay more attention to the influence of water flow from the confined aquifer of Luohe Formation and take corresponding water control projects. By plugging the groundwater recharge channel on fracture, the water yield in a coal mine could be reduced by 600 m3/h in a typical North China coal mine [28].
When the groundwater discharge was up to 250 m3/h, the water level of the Huanhe Formation aquifer and Luohe Formation aquifer dropped by about 11 m and 16 m, respectively. A descending funnel of about 2.4 km2 was formed in the drainage area, which affected the characteristics of local groundwater flow. Based on the numerical calculation, the descending trend of groundwater level for Cretaceous Huanhe Formation and Luohe Formation aquifers under different groundwater discharges is shown in Figure 10.
Compared with the groundwater level of the Luohe Formation in Tingnan Coal Mine near Hetaoyu Coal Mine, the groundwater level of the Luohe Formation dropped by 52.5 m from 2004 to 2015 [29]. However, with the improvement of water control works in each coal mine recently, the groundwater level has declined more and more slowly. Zhang (2015) obtained similar results using the large well method [30].
Therefore, it can be inferred that the groundwater level would decrease with the increase in drainage, but the groundwater drawdown rate decreased gradually. The downstream groundwater level suffered more influence.
The results are of great significance to guide the construction of mine water-control projects, such as curtain grouting, to protect local groundwater resources. According to the groundwater level control requirements of the government, the mine water inflow could be adjusted by engineering technologies so that the groundwater level depth of each aquifer was in a controllable range.

3.4. Prediction Groundwater Inflow Changes for Groundwater Drainage

The average coal seam thickness of the coal-2 layer, coal-5 layer, and coal-8 layer was 1.87 m, 0.95 m, and 8 m, respectively. Based on geological boreholes and survey data, the results of the maximum height of the water-conducted fracture zone in the coal-2 layer, coal-5 layer, and coal-8 layer were calculated by Formula (5). The results are shown in Table 2.
As the driving length of the main inclined shaft increases, the new water conducting would appear until the height of the water-conducted fracture zone was larger than the thickness of the above aquicludes. Then, the above geological structure model and the hydraulic conductivity would change accordingly. According to the trends of the driving length in Figure 5, the future groundwater drainage was predicted by Equation (6). The simulation result is shown in Figure 11.
With a linear increase in the driving length, when it reached 1410.4 m in the Hetaoyu coal field, the water-conducted fracture zone appeared, resulting in the groundwater drainage of 517.42 m3/h. It would generate another, new water-conducted fracture zone when the driving length reached 3083.2 m, and groundwater drainage increased to 806.83 m3/h. So, in order to ensure the safety of coal production, the drainage capacity of the pumps should be greater than 806.83 m3/h in future coal mining.
Taking the groundwater drainage with 806.83 m3/h as an example, the results are shown in Figure 12. The groundwater level of the Huanhe Formation aquifer and Luohe Formation aquifer would drop by about 16 m and 23 m, respectively. A descending funnel of about 4.5 km2 would be formed in the drainage area, which would affect the characteristics of local groundwater flow. The groundwater discharging from the Luohe Formation aquifer would account for 86%, about 693.88 m3/h.
However, when the groundwater drainage rate increased to 806.83 m3/h, the regional groundwater level would not drop drastically. It meant that the local groundwater drainage had little influence on the groundwater flow of the whole basin. Therefore, 4.5 km2 around the drainage area was the key groundwater protection area, and the coal mine departments should strengthen the monitoring of groundwater level and fracture zone in this area.

3.5. Uncertainty of Numerical Model

The numerical model presents a good simulation of groundwater flow, but it usually simplifies the hydrogeological structure, parameters, and recharges for convenience calculation. So, there were some uncertainties in the coupled numerical model.
Firstly, as shown in Figure 3, the stratigraphic structure was divided into seven sections, but, in natural conditions, each lithology is discontinuous horizontally or vertically, and the aquifers are distributed in the form of multilayer nesting. However, in this paper, the aquifers were uniformly divided, and the nearby aquifers were consolidated. So, the geological structure had some uncertainties, but the total thickness of each aquifer in this paper did not change in order to reduce its influence.
Secondly, the strata are heterogeneous in the mining area, and the hydrogeological parameters vary with time and space. In this paper, the values of specific storage and hydraulic conductivity were derived from the pumping test at five holes, but there would be a scale effect on applying point data to the whole area. In order to avoid the influence of uncertainty of hydrogeological parameters, the parameters were verified by comparing them with the water table of the hydrological hole based on the observed values obtained by the pumping test. Moreover, according to the pumping test results at different layers, the parameter partitions were refined to ensure the reliability of the water inflow model.
In addition, the recharge item in Formula (1) is only the mine drainage. The rainfall infiltration, evapotranspiration, and irrigation pumping were not calculated; this paper focused on the impact of mine drainage on the groundwater system, thus simplifying other factors. Considering the thickness of the Loess layer of about 200 m in the mining area, rainfall infiltration and evapotranspiration had little influence on the aquifers. Based on the results of this paper, the irrigation pumping also had little influence on the Huanhe Formation aquifers and Luohe Formation aquifers. Generally, the uncertainty of recharge had a limited impact on the results.
Therefore, there were some uncertainties in the numerical model. Although the corresponding treatment was carried out, its influence on the results still needs further study.

4. Conclusions

The numerical model was used to estimate the groundwater dynamics with groundwater drainage in the Longdong area, where the groundwater was the main source for residents and agriculture. Based on the observation data and simulation results, the spatial and temporal variations in the groundwater level for different aquifers were analyzed, and the sources of groundwater drainage were quantified. Moreover, the groundwater flow under the influence of the water-conducted fracture zone was predicted for future coal mining.
(1) The groundwater drainage changed irregularly, which was not linearly related to coal production. It was mainly affected by the hydrogeological structure and lithological parameters. The groundwater level at the location about 3 km away from the drainage area was not affected by the drainage;
(2) When the driving length of the main inclined shaft was excavated in all aquifers, the Luohe Formation aquifer was the main source of the mine water inflow, which accounted for about 92% of the total water inflow. So, the Luohe Formation aquifer was the main object to be protected and monitored by the local government;
(3) During the coal mining period, when the new water-conducting formed, the mine water inflow would increase to 806.83 m3/h under the influence of the water-conducted fracture zone, which would cause 23 m groundwater drawdown and descending funnel of about 4.5 km2. This drawdown might exceed the limit value regulated by the Qingyang Government. The Hetaoyu Coal mine departments should take some corresponding water prevention and control projects to reduce the drawdown of groundwater and pay more attention to monitoring the groundwater level within 4.5 km2.

Author Contributions

Conceptualization, L.C., X.W., G.L. and H.Z.; methodology, L.C. and X.W.; formal analysis, L.C. and X.W.; writing—original draft preparation L.C.; writing—review and editing, L.C., X.W., G.L. and H.Z.; visualization, L.C. and X.W.; supervision, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51722905), Youth Fund Program of Basic Research Operating Expenses of Central Public Welfare Research Institute for Nanjing Hydraulic Research Institute (Grant No. Y520024), and Young Top-Notch Talent Support Program of National High-level Talents Special Support Plan and China Water Resource Conservation and Protection Project (Grant No. 126302001000150005).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable.

Acknowledgments

The authors are grateful to the students who participated in this study, Jiaqi Sun for collaborating on data collection and question development, and Wenke Wang for suggestions on an early version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zeng, Y.F.; Wu, Q.; Liu, S.Q.; Zhai, Y.L.; Lian, H.Q.; Zhang, W. Evaluation of a coal seam roof water inrush: Case study in the Wangjialing coal mine, China. Mine Water Environ. 2018, 37, 174–184. [Google Scholar] [CrossRef]
  2. Husam, M.B.; Christopher, S.L. Applied Groundwater Modelling for Water Resource Management and Protection. Water 2022, 14, 1142. [Google Scholar] [CrossRef]
  3. Fan, C.; Du, Y.F.; Ni, J.Y.; Zhao, Z.R.; Peng, S.Q. Effect of Shallow-Buried High-Intensity Mining on Soil Water Content in Ningtiaota Minefield. Water 2021, 13, 361. [Google Scholar] [CrossRef]
  4. Li, P.; Wu, J.; Tian, R. Geochemistry, Hydraulic Connectivity and Quality Appraisal of Multilayered Groundwater in the Hongdunzi Coal Mine, Northwest China. Mine Water Environ. 2018, 37, 222–237. [Google Scholar] [CrossRef]
  5. Li, J.; Wang, W.K.; Cheng, D.W.; Li, Y.; Wu, P.; Huang, X.Q. Hydrogeological structure modelling based on an integrated approach using multi-source data. J. Hydrol. 2021, 600, 126435. [Google Scholar] [CrossRef]
  6. Zhang, Q.; López, D.L. Use of time series analysis to evaluate the impacts of underground mining on hydraulic properties in groundwater of Dysart Woods, Ohio. Mine Water Environ. 2019, 38, 566–580. [Google Scholar] [CrossRef]
  7. Roger, D.G.; Elias, H.T.; Hung, K.C. Regional Groundwater Modeling of the Guarani Aquifer System. Water 2020, 12, 2323. [Google Scholar] [CrossRef]
  8. Aalianvari, A. Combination of engineering geological data and numerical modeling results to classify the tunnel route based on the groundwater seepage. Geomech. Eng. A 2017, 13, 671–683. [Google Scholar] [CrossRef]
  9. de la Hera-Portillo, Á.; López Gutiérrez, J.; Mayor, B.; López-Gunn, E.; Henriksen, H.J.; Gejl, R.N.; Zorrilla-Miras, P.; Martínez-Santos, P. An Initial Framework for Understanding the Resilience of Aquifers to Groundwater Pumping. Water 2021, 13, 519. [Google Scholar] [CrossRef]
  10. Sun, N.Z.; Qi, G.M.; Li, J.; Liu, M.Q.; Sun, C. A New Approach for Predicting Mine Water Inflow from Inundation Coalfiled with Fracture and Faults. J. China Coal Soc. 1989, 2, 1–10. (In Chinese) [Google Scholar] [CrossRef]
  11. Liang, S.C.; Zhao, X.S.; Li, S.; Feng, L.H.; Liu, X.L.; Qin, X.M. The Application of GMS in Numerical Simulation of Groundwater and Faults Disposing in Gaizi River Source. Adv. Mater. Res. 2012, 518–523, 4047–4056. [Google Scholar] [CrossRef]
  12. Karimzade, E.; Sharifzadeh, M.; Zarei, H.R.; Shahriar, K.; Seifabab, M.C. Prediction of water inflow into underground excavations in fractured rocks using a 3D discrete fracture network (DFN) model. Arab. J. Geosci. 2017, 10, 206. [Google Scholar] [CrossRef]
  13. Mu, W.; Wu, X.; Deng, R.C.; Hao, Q.; Qian, C. Mechanism of Water Inrush Through Fault Zones Using a Coupled Fluid–solid Numerical Model: A Case Study in the Beiyangzhuang Coal Mine, Northern China. Mine Water Environ. 2020, 39, 380–396. [Google Scholar] [CrossRef]
  14. Song, W.; Liang, Z. Theoretical and numerical investigations on mining-induced fault activation and groundwater outburst of coal seam floor. Bull. Eng. Geol. Environ. 2021, 80, 5757–5768. [Google Scholar] [CrossRef]
  15. Xue, Y.; Liu, J.; Liang, X.; Wang, S.; Ma, Z.Y. Ecological risk assessment of soil and water loss by thermal enhanced methane recovery: Numerical study using two-phase flow simulation. J. Clean. Prod. 2022, 334, 130183. [Google Scholar] [CrossRef]
  16. Liu, G.W.; Ma, F.S.; Liu, G.; Guo, J.; Duan, X.L.; Gu, H.Y. Quantification of Water Sources in a Coastal Gold Mine through an End-Member Mixing Analysis Combining Multivariate Statistical Methods. Water 2020, 12, 580. [Google Scholar] [CrossRef]
  17. Anatoliy, F.; María del Refugio, B.-L.; Iryna, T.; Lydia, B.L.; Diego, A.P.; Javier, A.; Velázquez-Muñoz, F.Á. Thermal Regime and Water Balance of Two Tropical High-Mountain Lakes in the Nevado de Toluca Volcano, Mexico. Water 2022, 14, 1104. [Google Scholar] [CrossRef]
  18. Byrne, C.F.; Stormont, J.C.; Stone, M.C. Soil water balance dynamics on reclaimed mine land in the southwestern United States. J. Arid. Environ. 2017, 136, 28–37. [Google Scholar] [CrossRef]
  19. Elmarami, H.; Meyer, H.; Massmann, G. Combined approach of isotope mass balance and hydrological water balance methods to constrain the sources of lake water as exemplified on the small dimictic lake Silbersee, northern Germany. Isot. Environ. Health Stud. 2017, 53, 184–197. [Google Scholar] [CrossRef]
  20. Adeline, U.; Xi, C.; Tie, L.; Li, Z.Y.; Solange, U.; Richard, M.J.; Edovia, D.U.; Romaine, I.; Aline, U. Future Climate Change Impact on the Nyabugogo Catchment Water Balance in Rwanda. Water 2021, 13, 3636. [Google Scholar] [CrossRef]
  21. Wang, Q.; Wang, X.Y.; Liu, X.M.; Zhen, X.G.; Guo, J.W.; Chen, G.S.; Zhang, B. Prevention of Groundwater Disasters in Coal Seam Floors Based on TEM of Cambrian Limestone. Mine Water Environ. 2018, 37, 300–311. [Google Scholar] [CrossRef]
  22. Qiao, W.; Li, W.P.; Zhang, S.C.; Niu, Y.F. Effects of coal mining on the evolution of groundwater hydrogeochemistry. Hydrogeol. J. 2019, 27, 2245–2262. [Google Scholar] [CrossRef]
  23. Duan, X.L.; Ma, F.S.; Guo, J.; Zhao, H.J.; Gu, H.Y.; Liu, S.Q.; Sun, Q.H. Source Identification and Quantification of Seepage Water in a Coastal Mine, in China. Water 2019, 11, 1862. [Google Scholar] [CrossRef]
  24. Pilla, G.; Sacchi, E.; Zuppi, G.; Braga, G.; Ciancetti, G. Hydrochemistry and isotope geochemistry as tools for groundwater hydrodynamic investigation in multilayer aquifers: A case study from Lomellina, Po plain, South-Western Lombardy, Italy. Hydrogeol. J. 2007, 14, 795–808. [Google Scholar] [CrossRef]
  25. Powell, K.L.; Taylor, R.G.; Cronin, A.A.; Barrett, M.H.; Pedley, S.; Sellwood, J.; Trowsdale, S.A.; Lerner, D.N. Microbial contamination of two urban sandstone aquifers in the UK. Water Res. 2003, 37, 339–352. [Google Scholar] [CrossRef]
  26. Chen, Y.; Selvadurai, A.; Liang, W. Computational Modelling of Groundwater Inflow During a Longwall Coal Mining Advance: A Case Study from the Shanxi Province, China. Rock Mech. Rock Eng. 2018, 52, 917–934. [Google Scholar] [CrossRef]
  27. Cai, J.B.; Wang, P.; Shen, H.; Su, Y.; Huang, Y. Water Level Prediction of Emergency Groundwater Source and Its Impact on the Surrounding Environment in Nantong City, China. Water 2020, 12, 3529. [Google Scholar] [CrossRef]
  28. Lin, Y.; Wu, Y.Z.; Pan, G.Y.; Qin, Y.T.; Chen, G.S. Determining and plugging the groundwater recharge channel with comprehensive approach in Siwan coal mine, North China coal basin. Arab. J. Geosci. 2014, 8, 6759–6770. [Google Scholar] [CrossRef]
  29. Sun, K.; Fan, L.M.; Xia, Y.C.; Chen, J.; Gao, S.; Wu, B.; Peng, J.; Ji, J. Impact of coal mining on groundwater of Luohe Formation in Binchang mining area. Int. J. Coal Sci. Technol. 2021, 8, 88–102. [Google Scholar] [CrossRef]
  30. Zhang, Z. Estimation and Analysis of Mine Water Inflow in Hetaoyu Coalmine, Gansu. Coal Geol. China 2015, 27, 32–36, 56. (In Chinese) [Google Scholar] [CrossRef]
Figure 1. The location of the Hetaoyu coal field.
Figure 1. The location of the Hetaoyu coal field.
Water 14 02679 g001
Figure 2. The schematic diagram of the study area and the isolevel map of the Luohe Formation aquifer.
Figure 2. The schematic diagram of the study area and the isolevel map of the Luohe Formation aquifer.
Water 14 02679 g002
Figure 3. Diagram of the hydrogeological section of the mining area.
Figure 3. Diagram of the hydrogeological section of the mining area.
Water 14 02679 g003
Figure 4. Schematic diagram of calculation nodes of the groundwater flow model.
Figure 4. Schematic diagram of calculation nodes of the groundwater flow model.
Water 14 02679 g004
Figure 5. Time variations of groundwater drainage with coal mining in Hetaoyu Coal Mine.
Figure 5. Time variations of groundwater drainage with coal mining in Hetaoyu Coal Mine.
Water 14 02679 g005
Figure 6. Temporal variations of groundwater level under groundwater drainage in Hetaoyu Coal Mine. (a) Luohe aquifer (1#) and mixed groundwater level (2#). (b) Huanhe aquifer (3#).
Figure 6. Temporal variations of groundwater level under groundwater drainage in Hetaoyu Coal Mine. (a) Luohe aquifer (1#) and mixed groundwater level (2#). (b) Huanhe aquifer (3#).
Water 14 02679 g006
Figure 7. The correlation under different lag-time for #1 and #1.
Figure 7. The correlation under different lag-time for #1 and #1.
Water 14 02679 g007
Figure 8. The water table contour map of the Cretaceous Luohe Formation aquifers under the groundwater discharge was 187.6 m3/h.
Figure 8. The water table contour map of the Cretaceous Luohe Formation aquifers under the groundwater discharge was 187.6 m3/h.
Water 14 02679 g008
Figure 9. The time variation of aquifer reserve reduction under coal mining.
Figure 9. The time variation of aquifer reserve reduction under coal mining.
Water 14 02679 g009
Figure 10. Groundwater level drawdown under the different groundwater drainage.
Figure 10. Groundwater level drawdown under the different groundwater drainage.
Water 14 02679 g010
Figure 11. Prediction of the mine water inflow caused by the water-conducting fracture zone under coal mining.
Figure 11. Prediction of the mine water inflow caused by the water-conducting fracture zone under coal mining.
Water 14 02679 g011
Figure 12. The groundwater level contour map under the groundwater drainage rate was 806.83 m3/h: (a) the Huanhe Formation aquifers; (b) the Luohe Formation aquifer.
Figure 12. The groundwater level contour map under the groundwater drainage rate was 806.83 m3/h: (a) the Huanhe Formation aquifers; (b) the Luohe Formation aquifer.
Water 14 02679 g012
Table 1. Hydrogeological parameters of the main strata lithology in the mining area.
Table 1. Hydrogeological parameters of the main strata lithology in the mining area.
LithologySaturated Hydraulic ConductivityPorositySpecific YieldStorage Coefficient
m/dm−1
Loess0.0360.490.180.02
Silty mudstone0.0017280.210.120.000003
Siltstone0.000008640.050.210.000003
Fine sandstone0.00017280.10.230.000003
Mudstone0.0017280.210.120.000003
Coarse sandstone0.00087280.30.270.000003
Medium sandstone0.00157280.20.250.000003
Medium conglomerate00.10.130.000003
Clay0.0008640.450.030.001
Table 2. The maximum height of the water-conducted fracture zone of each coal seam in Haoheyu Coal Mine.
Table 2. The maximum height of the water-conducted fracture zone of each coal seam in Haoheyu Coal Mine.
Coal SeamDistance from Quaternary Bottom Plate (m)Normal Distance from the Bottom of the Lower Cretaceous Conglomerate
(m)
Dip Angle of Coal Seam
(°)
Coal Seam Thickness
(m)
Height of Water Conducted Fracture Zone
(m)
coal-2 layer690.9379.0465.0564.34
783.54104.7452.750.12
764.8588.7962.2128.35
coal-5 layer734.9890.1868.8686.78
840.86134.7495.2882.13
799.47123.4162.435.17
coal-8 layer748.49121.89723.56186.94
657.34104.43623.47185.47
832.02155.96613182.63
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, L.; Wang, X.; Liang, G.; Zhang, H. Evaluation of Groundwater Flow Changes Associated with Drainage within Multilayer Aquifers in a Semiarid Area. Water 2022, 14, 2679. https://doi.org/10.3390/w14172679

AMA Style

Chen L, Wang X, Liang G, Zhang H. Evaluation of Groundwater Flow Changes Associated with Drainage within Multilayer Aquifers in a Semiarid Area. Water. 2022; 14(17):2679. https://doi.org/10.3390/w14172679

Chicago/Turabian Style

Chen, Li, Xiaojun Wang, Gelong Liang, and Haicheng Zhang. 2022. "Evaluation of Groundwater Flow Changes Associated with Drainage within Multilayer Aquifers in a Semiarid Area" Water 14, no. 17: 2679. https://doi.org/10.3390/w14172679

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