Abstract
Modern mining within ecologically fragile areas is under enormous pressure from the ecological and geological environment. Using abandoned coal mines as underground water reservoirs is an effective idea for reducing eco-geological risks in mining districts. Accurate estimation of reservoir capacity is conducive to guiding the construction of an underground water reservoir. In this paper, estimation models of underground water storage capacity were established, which were divided into two kinds: horizontal mining district and inclined mining district. Furthermore, the minimum porosity of broken rock mass was analyzed based on laboratory experiments, and the fitted equation of minimum porosity was developed. A modern coal mine in northwest China was taken as an example, and the feasibility of the proposed model was verified by in-situ measurement. Finally, the effects of porosity of broken rock mass, water storage height and coal seam dip angle on water storage capacity were presented. The coal seam dip angle should be the principal consideration in choosing the site of an underground water reservoir.
HIGHLIGHTS
Underground water reservoir can realize sustainable development in mining areas.
Prediction model of underground water storage capacity was established.
The minimum void ratio of broken rock mass with different gradations was studied.
The coal seam dip angle significantly affects underground water storage capacity.
Graphical Abstract
INTRODUCTION
Coal resources make a great contribution to China's national economic development, and accounted for 56.0% of China's primary energy consumption in 2021 (Liu et al. 2022). It is considered that coal resources will continue to account for more than 50% of energy production and consumption before 2030 (Gu 2015; Xue et al. 2018; Fan et al. 2022). With the gradual exhaustion of coal resources, many collieries have or are facing closure in eastern and central China. The focus of coal mining has rapidly transferred to western mining districts of China, where the ecological environment is fragile, but the proven reserves and buried depth of the coal resource are massive and shallow, respectively (S. L. Liu et al. 2019). High-intensity and large-scale modern coal mining in thick coal seams must result in varying degrees of disturbance to the geological conditions, hydrogeological structure, ecological and natural environment within a mining district (Gunson et al. 2012; Fan et al. 2022). The delicate ecological balance will be destroyed when mining-induced disturbance exceeds the carrying capacity of the eco-geological system within the mining district. In recent years, mining-induced eco-geological environment and public health problems, such as water shortages, water pollution, vegetation death, desertification, and high human health risk, have received more and more attention (Wessman et al. 2014; Xue et al. 2018; Sun et al. 2020). In addition, statistics indicate that about 8 billion tons of mine water resource have been discharged by China's coal mining annually, and the utilization percentage of mine water is only about 25% (Gu 2015; Zhang et al. 2021). It is a huge waste of water resources, especially for a country where the per capita fresh water resource is far below that of the global average. Therefore, two key issues that must be solved for the establishment of ecological civilization and the sustainable development of a modern coal mining industry are water-resource-preserving mining and reasonable mine water treatment, especially within ecologically fragile mining districts (Mhlongo et al. 2018; Fan 2019).
To settle the two problems above, there have been several important studies from many disciplines using comprehensive research methods. Through site observation, Liu et al. (2022), Wei et al. (2017) and Fan et al. (2022) studied the forming process and distribution characteristics of mining-induced horizontal fracture and vertical fracture over the longwall panel in western mining districts of China. According to the affect degree of underground coal mining on the phreatic aquifer water table, four types of environmental engineering geological patterns, i.e., basically unaffected model, gradually restored model after destruction, gradually deteriorated model, and disaster model, were proposed and defined by Liu & Li (2019). Furthermore, according to field investigation of mining-induced hydrogeological conditions and normalized difference vegetation index, Chen et al. (2021) found that the emergence of phreatic water at the land surface in low-lying areas of the working faces may induce the formation of an oasis and wetland system in arid and semi-arid areas. In view of the lack of water resources in arid and semi-arid areas, Chi et al. (2021) presented the scale of coal mining in arid and semi-arid areas under the constraint of the water resources carrying capacity with the aim of realizing mining that conserves the ecological environment. Moreover, many researchers reported some interesting applications of water preservation mining techniques in ecologically fragile mining districts, for example underground water storage (Cairney 1973; Q. Liu et al. 2019; Song et al. 2020), limited-height mining (Miao et al. 2009), backfill mining (Liu et al. 2021), and room-and-pillar mining (Shi & Hou 2006). Among them, the most influential and promising method is underground water storage (Kallioras & Rusinski 2011; Li et al. 2014; Wang et al. 2018).
Sketch map of underground water reservoirs and mine water utilization (Zhang et al. 2021).
Sketch map of underground water reservoirs and mine water utilization (Zhang et al. 2021).
The present situation of underground water reservoirs is still in the early stages. Based on the investigation of the hydrogeological conditions and operation management of the Dagu River underground reservoir, a new assessment system for the operation effects of the underground reservoir was established by Li et al. (2019). Taking Bulianta coal mine as the research background, a discrete element fluid–solid coupling numerical simulation model was constructed to analyze the development characteristics of mining-induced fractures after coal seam mining, and the water replenishment channel of the coal mine underground reservoirs was determined by Chi et al. (2022). How to accurately assess the storage capacity of an abandoned coal mine is considered an urgent problem to settle for extending the technology (Li et al. 2019; Zhang et al. 2021). Mining scope and overburden destruction are important considerations in the storage capacity forecast (Wang et al. 2018; Song et al. 2020), but the effect of the compaction characteristics of broken rock mass and coal seam dip angle are ignored in most predictive models. The result is, in current practice, that most predictive models do not perform well. In this study, prediction models of underground water storage capacity in a flat seam mining district and inclined mining district are established, respectively, and the compaction characteristics of broken rock mass are introduced based on laboratory experiments. Furthermore, the accuracy of the prediction model is verified by engineering practice. Finally, the effects of different factors, i.e., porosity of broken rock mass, water storage height and coal seam dip angle, are discussed. This study can provide helpful guidance for site selection and construction of underground water reservoirs.
MECHANISM OF UNDERGROUND WATER RESERVOIR
Overburden destruction and water migration triggered by underground coal mining.
Overburden destruction and water migration triggered by underground coal mining.
The combined height of the caved zone and the fractured zone is the water flowing fractured zone (WFFZ) in which water can flow through. The water flowing fractured zone contains the majority of the void volume and has most of the water-conducting channels (Palchik 2003). If the height of the WFFZ penetrates the aquifer, i.e., the aquifuge is destroyed, water from the phreatic aquifer and surface flow into the working panel (Kim et al. 1997; Karaman et al. 2001; Liu et al. 2022), causing enormous pressure for the mine water drainage system and the eco-geological environment degradation problem on the surface (Figure 2). The change of abandoned underground void space into underground mine water storage area can effectively ameliorate and even eliminate local eco-geological environment risks, and also saves the costs of mine water drainage (Zhang et al. 2021).
Sketch map of the storage capacity loss of underground water reservoirs in a mining district (modified after Zhang et al. 2021).
Sketch map of the storage capacity loss of underground water reservoirs in a mining district (modified after Zhang et al. 2021).
PREDICTION METHOD OF UNDERGROUND WATER STORAGE CAPACITY
Storage capacity of reservoir in horizontal mining district
Sketch map for calculation of water storage capacity in horizontal mining area.
Storage capacity of goaf in horizontal mining district
The coefficients of Equation (2), i.e., c′ and d, can be obtained through the two data of porosity and thus the porosity of broken strata at other locations in the WFFZ can be predicted according to e(z). In addition, the porosity of the broken strata at the top interface of the WFFZ is 0, and the height of the WFFZ is readily accessible. So, it only takes one other datum of porosity to solve the coefficients of Equation (2).
Storage capacity of WFFZ in horizontal mining district
In addition, the parameters is of hc, hf, α and β are generally determined by field measurements or regional experience.
Storage capacity of reservoir in inclined mining district
The foregoing analysis reveals the prediction method of water storage capacity in a horizontal mining district. However, the majority of coal resources are inclined seams, and the inclined angle has a great impact on the distribution of water, which may decrease the accuracy of the above model. Hence, it is necessary to find a method in which the inclined angle has been taken into account.
Sketch map for calculation of water storage capacity in an inclined mining district.
Sketch map for calculation of water storage capacity in an inclined mining district.
FIELD APPLICATION
Background
The Lijiahao coal mine is located in the southeast part of Ordos city, Inner Mongolia Autonomous Region. The production capacity and service life of this coal mine are 6.0 Mt/a and 80 years, respectively. The eco-geological environment in this area is all too fragile and prone to be destroyed by high-intensity underground mining. To reduce the conflict between efficient coal mining and underground water loss, taking the mining-induced broken zone of No. 31108 and No. 31109 working face as water storage space, an underground water storage reservoir has been constructed. In this area, the coal mining seam has an elevation of about −256.19 m and the original structure is simple. The average mining thickness is 3.3 m, with average dip angle of 1.7°. The width and length of the water storage area are 581.3 m and 3,338.3 m, respectively. The roof of the coal mining seam is mainly sandstone and sandy mudstone The thickness of the caving zone and fractured zone is 19.8 m and 75.9 m (i.e., hc = 19.8 m, hf = 75.9 m) according to field measurements, respectively. In addition, the limit water head value of the artificial dam and coal pillar dam is 13.5 m according to the field test (Ju et al. 2017).
Minimum porosity of broken rock mass
Gradation of the broken rock mass
Testing apparatus
Experimental equipment: (a) the physical diagram; (b) schematic of the internal structure.
Experimental equipment: (a) the physical diagram; (b) schematic of the internal structure.
The test sample was a dry broken rock sample, and the test was carried out at room temperature. To remove moisture from the broken rock sample, the broken rock sample was dried in a drying oven before the test (the internal temperature of the drying oven was 50 °C, and the drying time was 2 h).
Test results and analysis
The p of the broken rock mass was calculated according to Equation (19), and the minimum porosity, emin, of the sample was summarized according to the test results.
The determination R2 coefficient of the fitting formula is 0.939, indicating that the minimum porosity of broken rock mass can be well predicted by Equation (20). The minimum value of emin is 0.121 according to Equation (20), which can be considered the porosity of the broken rock mass at the center of the goaf.
Calculation and validation
Firstly, according to laboratory experiments and Equation (20), the minimum porosity of broken rock mass at the center of the goaf is 0.121, i.e., e(z = 0 m) = 0.121. Furthermore, at the top of the fractured zone, the porosity narrows almost to 0, i.e., e(z = 99 m) = 0. Then the coefficients of Equation (2) are solved out.
Secondly, the height of the mine water is 13.5 m, while the mining thickness is 3.3 m, and the water storage area includes the goaf area and caving zone. In addition, the dip angle is 1.7°. So, the water storage volume Vc of this underground reservoir should be calculated by Equation (16). The water storage volume is 195,922.5 m3.
Mining engineers of the Lijiahao mining district had conducted a drainage experiment after injecting mine water into the goaf. When the height of the reservoir water storage is 13.5 m, the volume of water discharged is 185,759 m3. The relative error of the theoretical calculation and the measured values is only 5.5%, indicating the suitability of the proposed model.
DISCUSSION
In fact, water resources protection in arid and semi-arid regions has always been an important issue that needs to be solved. As a new type of hydraulic structure, underground reservoirs use abandoned goaf to store mine water, which has the characteristics of no land occupation, large storage capacity and small evaporation loss. Underground reservoirs allow for the rational allocation of water resources to meet the needs of coal enterprises in water-scarce areas. At the beginning of the construction of underground reservoirs in abandoned goaf, it is necessary to design the storage capacity to ensure that the water level and storage capacity of the underground reservoir remain within the safe range, so that the underground reservoirs are safe and reliable. The existing research results have been mainly based on empirical formulae, numerical analysis or actual monitoring of mine water inflow (Zhang et al. 2021). These studies often simplify the actual complex and changeable field conditions into a single assessment parameter, and lack of knowledge regarding the choice of important assessment parameter results in difficulties in assessing and designing the storage capacity of underground reservoirs. For example, the parameter of storage coefficient is widely used in the evaluation of reservoir volume in practical engineering applications (Wang et al. 2018; Song et al. 2020; Zhang et al. 2021), which is the comprehensive reflection of reservoir morphology, porosity and water storage height. However, there is no accurate and concrete formula for calculating a storage coefficient, which makes it impossible to determine the storage capacity reasonably (Wang et al. 2018; Song et al. 2020).
A theoretical framework of a storage capacity analytical method for an underground reservoir, which takes a full account of geological and mining conditions, is the focus of the current research. Based on this research gap, this study simplifies the WFFZ into a trapezoidal body, and solves the underground reservoir capacity of horizontal and inclined coal seams with porosity as the intermediate quantity. On this basis, the model in this study organically combines the geological conditions and mining conditions, i.e., mining thickness, mining area, porosity of broken rock mass, coal seam dip angle and the development of the WFFZ. The model is more accurate and credible in the process of practical application due to calculation parameters that are much easier to access and choose. In addition, the model in this study can be used to analyze the impact of different physical parameters on water storage capacity of underground reservoirs from the perspective of quantification.
Effect of porosity on water storage capacity
Water storage capacity of the underground reservoir versus porosity.
Effect of storage height on water storage capacity
Similarly, the relationship between water storage height and water storage capacity has been obtained according to the proposed model. Significantly, the dip angle is 1.7° for all the cases.
Water storage capacity of the underground reservoir versus water storage height.
Water storage capacity of the underground reservoir versus water storage height.
Effect of dip angle on water storage capacity
Water storage capacity of the underground reservoir versus dip angle.
At the same water storage height, the results in Figure 13 indicate that the angle of the underground reservoir has a marked impact on the water storage capacity. The value of water storage capacity decreases with the increase of dip angle value. The water storage capacity attenuates largely before the dip angle reaches 1°, and then tends toward stability with increasing dip angle. The water storage capacity of the underground reservoir reaches a maximum of 2,990,954.4 m3 when the dip angle is 0°, while the value of water storage decreases to 353,260.3 m3 when the dip angle increases to 1°. The change of the water storage capacity of the underground reservoir is 88.19% when the range of the dip angle is from 0° to 1°. That can explain the shape and sectional area of underground reservoir decreasing considerably as dip angle increases. Although the maximum height of the water allowed to be stored in an underground reservoir, which depends on the strength of the dam, can be different at different dipping angles in actual engineering, the results in Figure 13 suggest that the dip angle should be the principal consideration in choosing the site of an underground reservoir.
CONCLUSIONS
Underground water reservoirs provide a new and efficient method of achieving green mining and sustainable development in ecologically fragile mining districts. In this paper, the prediction model of water storage capacity was researched in detail based on laboratory experiments, theoretical calculation and field verification. The main conclusions are summarized as follows.
- (1)
For a specified water storage height, the storage capacity of underground reservoirs depends on the shape of the underground reservoirs and the distribution of porosity. Through simplifying the shape of the WFFZ into a trapezoid, the prediction model of underground water storage capacity in a flat seam mining district and inclined mining district are established through an empirical equation of porosity of the broken strata.
- (2)
Taking Lijiahao underground reservoir as an example, the compaction characteristics of broken rock mass are introduced based on laboratory experiments. The porosity of the broken rock mass (emin) shows a parabolic curve relationship with the characteristic parameter of the gradation curve (p). The emin of the broken rock mass can be well predicted by the fitting curve. The minimum value of emin is 0.121 in Lijiahao coal mine, which is considered to be the porosity of the broken rock at the goaf center. The predicted and experimental results of water storage capacity are 195,922.5 m3 and 185,759 m3 respectively. Measured results indicate that the error between predicted and measured values is only 5.5%. The validation and correctness of the prediction model is demonstrated.
- (3)
The effect of different physical parameters on the water storage capacity of underground reservoirs is researched based on the background of Lijiahao underground reservoir. The storage capacity of underground water storage is positively regression correlated with the porosity of broken rock mass and water storage height, and negatively regression correlated with the coal seam dip angle. The coal seam dip angle should be the principal consideration in choosing the site of an underground reservoir, because a small increase of the dip angle could lead to a significant decrease in the water storage capacity of the underground reservoir.
ACKNOWLEDGEMENT
Financial support for this work is provided by the Central Public-Interest Scientific Institution Basal Research Fund (No. Y320010) and Natural Science Foundation of Jiangsu Province (BK20190646).
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.