Prediction and Maintenance of Water Resources Carrying Capacity in Mining Area—A Case Study in the Yu-Shen Mining Area

: The problem of water resources damage caused by coal mining has restricted the sustainable development of Yu-Shen mining area. Illustrating the relationship between mining and water resources carrying capacity is of great signiﬁcance to solve this problem. In this study, the authors proposed an appraisal and prediction model of water resource carrying capacity in the mining area (WRCCMA) based on the analytic hierarchy process (AHP)-fuzzy comprehensive evaluation method. A triple-leveled structure model was developed, and the main inﬂuencing factors of the WRCCMA and the membership functions were analyzed. The prediction model was applied to Yubujie colliery to test its validity by investigating the changes of vegetation coverage and the ground deformation of the colliery and its adjacent coal mine before and after mining. Subsequently, we obtained the WRCCMA of the study area and zoning map of di ﬀ erent grades of WRCCMA in the mining area by applying this model to the whole Yu-Shen mining area. Furthermore, three countermeasures to maintain the WRCCMA and realize water conservation coal mining (WCCM) were provided to collieries with di ﬀ erent WRCCMA grades, including mining methods selection, mine water reutilization, and water-resisting layer reconstruction. Reasonable mining methods and water-resisting layer reconstruction can reduce the development of water conductive fractures and thus prevent groundwater from penetrating into the goaf. Mine water reutilization provides a source of water demand for collieries and families, contributing to the reduction of abstraction of water resources. These three countermeasures can help to maintain the WRCCMA. This paper successfully combines the fuzzy theory with mining engineering and provides theoretical and practical guidance for other mining areas in arid and semi-arid regions of Northwest China.


Introduction
With the depletion of coal resources in central and eastern China, the focus of coal mining has gradually shifted to arid and semi-arid areas in Northwest China, where the proven reserves make up to more than two thirds of the country's total, and contributes to over one third of China's annual coal production [1,2]. However, Northwest China has a dry climate, sparse vegetation, and water resources are scarce, accounting for only 3.9% of the entire country [3]. Water resources are likely to be

Geological Characteristics
The comprehensive stratigraphic column of Yu-Shen mining area, as shown in Figure 2. The aquifer is mainly Salawusu Formation, usually forming a complete aquifer with overlying Aeolian sand Formation. The aquifuge has consisted of Lishi Formation and Baode Formation. Yan'an Formation is the coal-bearing Formation, which has five main mineable coal seams [4,45].

Underground Aquifer Characteristics
From top to bottom, there are total four types of aquifers in the Yu-Shen mining area, including: the unconsolidated porous phreatic aquifer, Salawusu phreatic aquifer, burnt rock phreatic aquifer, and porous bedrock confined aquifer, as shown in Figure 3. The unconsolidated porous phreatic aquifer is thin and is mainly recharged from precipitation, usually forming a complete aquifer with the underlying Salawusu Formation aquifer. The porous bedrock confined aquifer is characterized by a minor dimension of pore, high permeability, extremely low water abundance, and limited distribution range. The burnt rock phreatic aquifer itself cannot form a water-storing structure and is primarily recharged from the Salawusu phreatic aquifer. The Salawusu Formation aquifer is widely distributed in the mining area, and its thickness ranges from 0 m to 67.3 m. The buried depth of water table is generally less than 10 m [4,45]. As a result, it is the major aquifer in this area which needs to be considered during the evaluation of WRCCMA.

Identification of Influencing Factors of WRCCMA
The influence factors of WRCCMA are fuzzy and complex, which is consistent with the characteristics of fuzzy mathematics capable of solving many multi-leveled and fuzzy problems [11]. The fuzzy comprehensive evaluation method is based on fuzzy mathematics. It can transform qualitative evaluation into quantitative evaluation using the membership degree theory of fuzzy mathematics and is used to make an overall evaluation of objects restricted by many factors. The analytic hierarchy process (AHP)-fuzzy comprehensive evaluation method was employed and a triple-leveled structure model was constructed to identify the influencing factors of WRCCMA and calculate their weights. In the AHP Sustainability 2020, 12, 7782 5 of 27 model, six factors were selected as secondary indicators, i.e., overburden strata system, geological structure system, groundwater system, coal mining system, ecological system, and social system. Meanwhile, 15 factors, including the thickness of aquifuge, fault strength, mining method, vegetation coverage, and water resources per capita, were chose as tertiary indicators, as shown in Figure 4. By analyzing the results of field investigation and the existing research, we constructed the membership function of each influencing factor, which was the basis of evaluating the WRCCMA [11,49].

Identification of Influencing Factors of WRCCMA
The influence factors of WRCCMA are fuzzy and complex, which is consistent with the characteristics of fuzzy mathematics capable of solving many multi-leveled and fuzzy problems [11]. The fuzzy comprehensive evaluation method is based on fuzzy mathematics. It can transform qualitative evaluation into quantitative evaluation using the membership degree theory of fuzzy mathematics and is used to make an overall evaluation of objects restricted by many factors. The analytic hierarchy process (AHP)-fuzzy comprehensive evaluation method was employed and a triple-leveled structure model was constructed to identify the influencing factors of WRCCMA and calculate their weights. In the AHP model, six factors were selected as secondary indicators, i.e., overburden strata system, geological structure system, groundwater system, coal mining system, ecological system, and social system. Meanwhile, 15 factors, including the thickness of aquifuge, fault strength, mining method, vegetation coverage, and water resources per capita, were chose as tertiary indicators, as shown in Figure 4. By analyzing the results of field investigation and the existing research, we constructed the membership function of each influencing factor, which was the basis of evaluating the WRCCMA [11,49].

Overburden Strata System
1. The distance from aquifer to coal seam b With the working face advancing, the overlying strata above the goaf will move and produce water conducted fractures because of losing the support from coal body. When the face reaches critical mining, if b is greater than the height of the WCFZ, the groundwater will not penetrate into the mined-out area, and a water inrush accident will not occur. However, when it is less than the height of the WCFZ and greater than the height of the caved zone, the fractures will develop and reach the aquifer. The groundwater table lowering and water resources loss will happen, but sand gushing in mines will be impossible. With the continuous advancement of the face after critical mining, the water flowing fissures will naturally be compacted and closed, and the water level may return to its original level on account of the recharge from precipitation, lateral runoff, and irrigation return flow. When b is less than the caved zone height, mining is prone to causing water and sand inrush. The groundwater enters a vicious cycle and water level will not recover [9][10][11]. In this situation, the grade of WRCCMA is the lowest. Therefore, the membership function is as follows: where Hk refers to the height of the caved zone (m); Hd refers to the height of WCFZ (m); b refers to the distance from an aquifer to coal seam (m).

Aquifuge thickness
The water-resisting property of the aquifuge is also an important indicator to evaluate the WRCCMA. The stronger the water-resistant ability, the lower risk of the groundwater will be lost due to mining. The scholars and experts engaged in mining generally believe that clay has a strong water-resistant property, and the clay stratum is usually thought of as an impermeable layer that the WRCCMA  With the working face advancing, the overlying strata above the goaf will move and produce water conducted fractures because of losing the support from coal body. When the face reaches critical mining, if b is greater than the height of the WCFZ, the groundwater will not penetrate into the mined-out area, and a water inrush accident will not occur. However, when it is less than the height of the WCFZ and greater than the height of the caved zone, the fractures will develop and reach the aquifer. The groundwater table lowering and water resources loss will happen, but sand gushing in mines will be impossible. With the continuous advancement of the face after critical mining, the water flowing fissures will naturally be compacted and closed, and the water level may return to its original level on account of the recharge from precipitation, lateral runoff, and irrigation return flow. When b is less than the caved zone height, mining is prone to causing water and sand inrush. The groundwater enters a vicious cycle and water level will not recover [9][10][11]. In this situation, the grade of WRCCMA is the lowest. Therefore, the membership function is as follows: where H k refers to the height of the caved zone (m); H d refers to the height of WCFZ (m); b refers to the distance from an aquifer to coal seam (m).

Aquifuge thickness
The water-resisting property of the aquifuge is also an important indicator to evaluate the WRCCMA. The stronger the water-resistant ability, the lower risk of the groundwater will be lost due to mining. The scholars and experts engaged in mining generally believe that clay has a strong water-resistant property, and the clay stratum is usually thought of as an impermeable layer that the water bodies cannot flow through by gravity [31]. The authors found that the thickness of the aquifuge is the most important influencing factor of the water-resisting property. If it is less than 12 m, groundwater has a tendency to infiltrate into gob through the aquifuge, while if it is greater than 40 m, it can effectively prevent water seepage [4,32,57].

Geological Structure System
The fault lowers the integrity and mechanical properties of the rock mass. The strata near the fault are usually broken, loose, and fractures are well-developed [33]. The faults can easily be connected with mining-induced water flowing fissures and form a complete water flowing channel for groundwater penetrating into the goaf. Moreover, those faults with high fault throw and long extension length also provide a large space for groundwater storage [34]. When the face advances ahead and the water conductive fractures reach the fault, the water bursting into the mine will take place. The greater the fault strength, the more the fractures are developed, which is unfavorable to WCCM, and therefore the lower the grade of WRCCMA.
The fault strength was selected as the evaluation indicator of WRCCMA. The formula of fault strength is as follows [35]: The membership function of fault strength is [36]:

Buried depth of groundwater table
The arid and semi-arid areas in Northwest China are ecologically fragile, and water resources are incredibly scarce [4]. Surface vegetation needs to absorb water from the aquifer through the root system to meet its own growth needs. The buried depth of aquifer can directly affect the luxuriance, degradation of plants, as well as the type, composition, and quantity of vegetation [58]. When the depth is too shallow, the water contained in the soil will flood the vegetation root system, which is not conducive to the respiration of the root system, thereby inhibiting the growth of vegetation. Furthermore, the shallow depth can easily lead to soil salinization and is harmful to the growth of vegetation. On the other hand, when the depth exceeds the critical depth from which the vegetation root can absorb groundwater, the vegetation will die due to water shortage [59]. The change of water level (D) before and after coal exploitation was selected as an evaluation indicator.

Recharge of groundwater
In arid and semi-arid areas, surface runoff is limited and the groundwater provides the important water supply for production and living. The groundwater is mainly recharged from precipitation, followed by lateral flow and irrigation infiltration. Thus, the coefficient of recharge from precipitation was selected as the evaluation standard and its formula is as follows [60]: where α refers to the coefficient of recharge from precipitation; µ refers to the specific yield; ∆h refers to the rise of water level caused by cumulative recharges from precipitation in one year (mm); R refers to annual precipitation (mm). The membership function is:

Groundwater mineralization (M)
When the buried depth of groundwater is within the suitable growth range of plants, the groundwater salinity has a significant impact on the growth of plants. The critical salinity of different plants' withering and death vary. The suitable range for most vegetation is 0.5-3 g/L, and 3-5 g/L will not inhibit the vegetation growth. When it is greater than 8 g/L, most of the plants will perish and die [59,60]. The water abundance of an aquifer is related to its thickness and permeability. The greater the thickness and permeability of the aquifer, the stronger the ability for the aquifer to store water and receive replenishment and thereby the higher the water abundance. The water yield property has a positive correlation with the grade of WRCCMA. Based on the grading standard [4], the water-richness is divided into three grades, low, medium, and high, the corresponding uniform drawdown unit flow is q < 0.1 L·s −1 ·m −1 , 0.1 L·s −1 ·m −1 < q ≤ 5.0 L·s −1 ·m −1 , and q > 5.0 L·s −1 ·m −1 , respectively.
High occurrence 0.7, Medium occurrence 0.3, Low occurrence 0, No occurrence (9) 5. Groundwater quality In the course of coal exploitation, the coal mine solid waste is accumulated to form a waste heap on the surface. The leaching water formed after the solid waste is washed, leached, or soaked will seriously pollute the shallow water and leads to water quality deterioration. This phenomenon will not only threaten the health of local citizens, but also do damage to the regional ecological environment conservation, contributing to the lowering of the WRCCMA grade. The comprehensive pollution index p was selected as the evaluation standard [11]:

Mining methods
The WCFZ heights of various mining methods are different [26]. As a consequence, the degrees of mining's influence on overlying aquifer are varying. Longwall full-seam mining will only be implemented when a coal seam is thin while the overlying strata are thick, and the height of mining-induced WCFZ is less than the distance between the aquifer and the coal seam. Among all mining methods, the longwall full-seam mining has the lowest ability to preserve water bodies above colliery [27]. In arid and semi-arid mining areas, coal mining under aquifers usually adopts WCCM methods, including backfill mining, partial mining, and slice mining [28]. When other conditions are the same, compared to full-seam mining, slice mining can reduce the height of the WCFZ to a limited extent [29]. Partial mining, i.e., height-restricted mining, strip mining, and room and pillar mining, can reduce movement and deformation of the overlying strata and thereby the height of WCFZ. However, all of these partial mining are realized by abandoning coal pillars and restricting the mining thickness of coal seam, wasting large amounts of good coal resources [30]. The backfill mining replaces coal pillars by filling bodies to fill up the gob and support the roof and overlying strata. It is an effective mining method to control strata movement together with the development of water conducted fractures, and thus can better preserve overlying aquifer.

Mining parameters
Mining parameters, including mining height, advancing speed, width, and advancing length of longwall panel, are important factors of the WCFZ height [4]. If the water flowing fractures penetrate the overlying aquifer, groundwater water loss will happen and thereby lower the grade of WRCCMA. With the rapid advancement of the panel, the caved zone is gradually compacted [4]. Hence, the water flowing channel can be completely closed, and the water level above the goaf may recover progressively. The authors studied the relationship between the mining speed and the change of the groundwater table. The extent of crack development and penetrance under the condition of mining speed greater than 15 m/d is weaker than that less than 10 m/d. When the mining speed is less than 10 m/d, the water level drops significantly, while when it is greater than 15 m/d, its effect on the water level is not obvious. Furthermore, the greater the width and the advancing length of the coal face, the easier it is for the overlying strata to sink as a whole, and the less likely it is for water conductive fissures to develop. According to the results of on-site investigation and combined with the existing research, the multivariate membership function of mining parameters is [4,11]: where x 1 refers to the mining height; x 2 refers to the mining speed; x 3 refers to the width of face; x 4 refers to the advancing length of face.

Surface deformation
The surface deformation will exert influence on the WRCCMA. Mining-induced land subsidence will alter the landform and the direction of rivers and affect the normal growth of plants, resulting in vegetation stunting and even vegetation coverage reduction. Moreover, land cracks developing down to the aquifer will change the coefficient of recharge from precipitation and affect the recharge, runoff, and discharge of groundwater. Hence, the original balance and circulation state of groundwater will be broken. The subsidence coefficient q was used as evaluation standard and its formula is as follows [5,10]: where W max refers to the maximum subsidence value after critical mining (m); M refers to the thickness of coal seam (m); α refers to the dip angle of coal seam ( • ). In light of the subsidence factors of different lithological strata measured in some coal mines, the membership function is: 3.5. Ecological System

Vegetation coverage (β)
Vegetation plays a significant role in maintaining ecological balance and curbing land desertification. There are three types of water supplies for vegetation, i.e., atmospheric precipitation, surface water, and groundwater, among which the groundwater is the main source of water supply because rainfall and surface runoff are minimal. The circulation and balance of water resources control the growth and types of vegetation, alteration in circulation and balance will unavoidably lead to changes in vegetation coverage. Vegetation coverage is closely related to the WRCCMA and is an intuitive expression of it. Up to now, there has been no acknowledged classification standard for vegetation coverage in the Yu-Shen mining area. Based on the ecological characteristics of arid and semi-arid areas and referring to the classification standards adopted in neighboring areas, its grade was divided into five levels [59]:

Total water resources
The total water resources are closely correlated to the eco-environment and human social activities and are an important factor affecting the WRCCMA. The draining depressurization, which uses artificial drainage to lower the water level and head pressure, makes the amounts of water resources continuously decrease for a long period of time [60]. Besides this, as the evaporative capacity in arid and semi-arid area is high, the discharged mine water will suffer from evaporation loss. Taking the ratio of total water resources after and before mining (η) as the evaluation standard: 3.6. Social System

Water utilization rate (υ)
The utilization rate of water resources is a main indicator that reflects whether the WRCCMA has the potential for improvement. Taking into account the requirements of eco-environment conservation and biodiversity, the development and utilization rate of water resources should not exceed 30%. If it exceeds 40%, the sustainable development of water resources will be affected [17].
2. Water resources per capita (ϕ) Water resources per capita are an intuitive reflection of WRCCMA. A low level means that water resources cannot meet the water demand of humans and thus, the grade of WRCCMA is low, and vice versa. Water resources per capita of extremely severe, severe, medium, and slight water shortage areas are <500 m 3 , 500-1000 m 3 , 1000-2000 m 3 , 2000-3000 m 3 , respectively [23].

Mathematical Modeling
In the AHP multilevel model, level A represents the overall goal, level B denotes the level of evaluation criteria, and level C is comprised of sub-criteria. Judgment matrices were developed by this multilevel structure to realize the numerical expression and quantification of WRCCMA [61]. The evaluation factors are as follows: Let the fuzzy subset V denote the discourse domain of WRCCMA. The Φ value, indicating the level of WRCCMA, is the membership degree of U in V. The domain of V is defined by Equation (20): The following equation can calculate the Φ value: where u i (u i ) refers to the membership degree of the i-th factor and w i refers to the weight of the i-th factor. The grade of WRCCMA can be classified as carrying surplus, capable of carrying, moderate carrying, slightly over-carrying, and severely over-carrying [4,11], as shown in Table 1. In order to keep the WRCCMA and preserve the eco-environment from continuous deteriorating, mining must be stopped, since its impacts on the WRCCMA have already exceeded the critical state.

Weight Determination
Researchers and scholars engaged in coal mining, hydrogeology, eco-environment, and WRCC have been invited to assign the weights of the factors of this AHP model of WRCCMA. In terms of the relative weight evaluation table of the factors from the experts, the following judgment matrices were obtained: The relative weight evaluation of sub-factors of the groundwater system from the experts was taken as an example. There are three steps totally to calculate the weights of these sub-factors of matrix W B 3 ∼C .
Step 1: Calculating the largest eigenvalue and its eigenvector The largest eigenvalue λ max of the matrix W B 3 ∼C is 5.0059 and the eigenvector W = (0.4085, 0.0977, 0.0854, 0.2042, 0.2042).
Step 2: Testing the random coincidence coefficient C.R.
Equations (22) and (23) were used to conduct the consistency test: where C.I. refers to the consistency indicator; λ max refers to the largest eigenvalue; n refers to the number of the influencing factors of this matrix.
where C.R. refers to the consistency ratio and R.I. refers to the average consistency index. If C.R. < 0.1, the relative weights are reasonable. Otherwise, the judgment matrix needs to be adjusted by reassigning the weights of these sub-factors. The C.R. of the matrix W B 3 ∼C is 0.0013 < 0.1, indicating the relative weight from the experts is reasonable and acceptable.
Step 3: Calculating the weights of the influencing factors in the premise of C.R. passing the check.
The weight of the groundwater system is 0.3795. Then, 0.3795 × W = (0.1550, 0.0371, 0.0324, 0.0775, 0.0775) is the weight of the five sub-factors of groundwater system.
The weights of the other five matrixes were also calculated using this method. The weight distribution is as shown in Table 2.

Numerical Simulation of Groundwater Level Decline
The mathematical evaluation model was applied to Yubujie colliery to verify its reliability and rationality. The narrow strip WCCM method was used in Yubujie colliery. The narrow strip that was mined was 12 m wide, and the strip coal pillars left unmined was 8 m wide [4], as shown in Figure 5. Since there were no data on water level drop in this colliery, fluid-solid coupling module of Fast Lagrangian Analysis of Continua (FLAC 3D ) finite element software was employed to simulate groundwater table lowering under narrow strip mining. The numerical simulation model was as shown in Figure 6. Mined-out narrow strip   The Mohr-Coulomb model was used during the calculation. The "CONGIG fluid" code was used to enter the seepage mode, and the "Initial pp" code was employed to set the pore pressure and pore water pressure gradient. The seepage mode was set as the isotropic. The permeability coefficient of the clay (aquifuge) was 10 -13 m 2 /Pa-sec. Other hard rocks, including sandstone, fine sandstone, medium sandstone siltstone, gritstone, and mudstone had a larger permeability coefficient than clay. The porosity was set to default value 0.5, and the tensile strength of fluid was set to 10 15 Pa. The Biot coefficient was set to 1, and the saturation was set to 1. The top sandstone (aquifer) was set as a free surface, and the surrounding boundary and bottom boundary were set as permeable boundaries by using the "Fix pp" code. Fluid could flow into or out of the model along the permeable boundary. The mechanical parameters of the rock were as shown in Table 3. The groundwater depression cone and the contour map of water table drop were as shown in Figure 7. The numerical simulation results show that the water level dropped by 1.32 m, suggesting a narrow variation in groundwater levels.

Narrow strip face
Mined-out narrow strip The front and rear sides of the model were fixed in the X direction, the left and right sides were fixed in the Y direction, and the bottom boundary was fixed in the Z direction. The Mohr-Coulomb model was used during the calculation. The "CONGIG fluid" code was used to enter the seepage mode, and the "Initial pp" code was employed to set the pore pressure and pore water pressure gradient. The seepage mode was set as the isotropic. The permeability coefficient of the clay (aquifuge) was 10 -13 m 2 /Pa-sec. Other hard rocks, including sandstone, fine sandstone, medium sandstone siltstone, gritstone, and mudstone had a larger permeability coefficient than clay. The porosity was set to default value 0.5, and the tensile strength of fluid was set to 10 15 Pa. The Biot coefficient was set to 1, and the saturation was set to 1. The top sandstone (aquifer) was set as a free surface, and the surrounding boundary and bottom boundary were set as permeable boundaries by using the "Fix pp" code. Fluid could flow into or out of the model along the permeable boundary. The mechanical parameters of the rock were as shown in Table 3. The groundwater depression cone and the contour map of water table drop were as shown in Figure 7. The numerical simulation results show that the water level dropped by 1.32 m, suggesting a narrow variation in groundwater levels.

The Evaluation of the WRCCMA for Yubujie Colliery
The distance from the aquifer to the coal seam of Yubujie colliery is 220 m on average, and the thickness of the aquifuge varies from 20 to 90 m. The height of the caved zone and the height of WCFZ are 33.9 m and 158.8 m, respectively. There is no large fault in the Yubujie colliery and the fault strength is less than 0.22. After mining, the water level has decreased by 1.32 m, and the coefficient of recharge from precipitation is 18.4%. Moreover, the salinity ranges from 270.2 to 675.4 mg/L and the aquifer is highly water-rich. The comprehensive water pollution index p is less than 0.2. The thickness of the coal seam varies from 4.1 to 7.2 m, with an average of 6.0 m, and the subsidence coefficient is 0.18. The vegetation coverage of Yubujie colliery is 64%, and the ratio of total water resources after mining and before mining is 90.4%. The water resources utilization is 34.04%. The per capita water resources vary from 1000 to 2000 m 3 /person [57]. The comprehensive evaluation value Φ was calculated to be 0.84, indicating the WRCCMA is capable of carrying.

The Changes in Vegetation Coverage and Surface Deformation after Mining
Vegetation coverage and landform alteration is the direct reflection of the WRCCMA. An on-site investigation was conducted in Yubujie colliery, as shown in Figure 8. The area framed with a red dash line is a part of Yubujie colliery, which used the narrow strip WCCM method. The vegetation grew luxuriantly after mining in this area, showing that the grade of the WRCCMA is high, and water resources are sufficient to meet the demand of the surface vegetation. By contrast, the area framed in the yellow dash line is the adjacent coal mine of Yubujie colliery, which used the traditional longwall

The Evaluation of the WRCCMA for Yubujie Colliery
The distance from the aquifer to the coal seam of Yubujie colliery is 220 m on average, and the thickness of the aquifuge varies from 20 to 90 m. The height of the caved zone and the height of WCFZ are 33.9 m and 158.8 m, respectively. There is no large fault in the Yubujie colliery and the fault strength is less than 0.22. After mining, the water level has decreased by 1.32 m, and the coefficient of recharge from precipitation is 18.4%. Moreover, the salinity ranges from 270.2 to 675.4 mg/L and the aquifer is highly water-rich. The comprehensive water pollution index p is less than 0.2. The thickness of the coal seam varies from 4.1 to 7.2 m, with an average of 6.0 m, and the subsidence coefficient is 0.18. The vegetation coverage of Yubujie colliery is 64%, and the ratio of total water resources after mining and before mining is 90.4%. The water resources utilization is 34.04%. The per capita water resources vary from 1000 to 2000 m 3 /person [57]. The comprehensive evaluation value Φ was calculated to be 0.84, indicating the WRCCMA is capable of carrying.

The Changes in Vegetation Coverage and Surface Deformation after Mining
Vegetation coverage and landform alteration is the direct reflection of the WRCCMA. An on-site investigation was conducted in Yubujie colliery, as shown in Figure 8. The area framed with a red dash line is a part of Yubujie colliery, which used the narrow strip WCCM method. The vegetation grew luxuriantly after mining in this area, showing that the grade of the WRCCMA is high, and water resources are sufficient to meet the demand of the surface vegetation. By contrast, the area framed in the yellow dash line is the adjacent coal mine of Yubujie colliery, which used the traditional longwall full-seam mining method. It can be seen that after mining, a small subsidence basin appeared above the coal face, and the trend of vegetation deterioration and land desertification is conspicuous. Hence, from the perspective of the changes in the vegetation coverage and the ground deformation after mining, the impacts of different mining methods on the WRCCMA are also obviously different. This also confirms that it is reasonable for experts to give high scores to mining methods in the fuzzy mathematical evaluation model. full-seam mining method. It can be seen that after mining, a small subsidence basin appeared above the coal face, and the trend of vegetation deterioration and land desertification is conspicuous. Hence, from the perspective of the changes in the vegetation coverage and the ground deformation after mining, the impacts of different mining methods on the WRCCMA are also obviously different. This also confirms that it is reasonable for experts to give high scores to mining methods in the fuzzy mathematical evaluation model.

WRCCMA Prediction for Yu-Shen Mining Area
By analyzing the collected data as well as the existing research, the membership degree of each evaluation indicator was determined. Subsequently, combined with the weight distribution, we calculated the comprehensive evaluation value Φ of each drill hole under longwall full-seam mining in the Yu-Shen mining area.

Overburden strata system
Based on the drilling data, we obtained the distance from the aquifer to the first-mined coal seam by subtracting the buried depth of the Salawusu Formation aquifer from the buried depth of the firstmined coal seam, and its contour map was plotted using Kriging method, as shown in Figure 9.
Many scholars have done a lot of theoretical and practical research on the height of WCFZ in East China. The empirical formulae for calculating the height of WCFZ, which are suitable for different overburden lithology and dip angles of coal seams, were given when adopting longwall mining method [27][28][29]. However, most of the main mineable coal seams in eastern mining area belong to Permian system, while the Yu-Shen mining area is mainly Jurassic coal seam [4,30]. The hydrogeological and geological engineering conditions in this area are significantly different from those in the eastern mining area [60]. Thus, the traditional prediction formulae of WCFZ height are inapplicable to Yu-Shen mining area. In this paper, field measurement and numerical simulation were employed to determine the heights of the caved zone and the WCFZ. The caved zone height and WCFZ height of Yuyang colliery, Yushuwan colliery, Jinjitan colliery, and Hanglaiwan colliery were obtained by drilling fluid consumption observation assisted with the core engineering geological logging, water pressure test in the borehole, television logging, and geophysical logging. The mining height of these coal mines varies from 3.5 to 5.5 m, which cannot represent the whole Yu-Shen mining area. Therefore, the universal distinct element code (UDEC) was used to establish

WRCCMA Prediction for Yu-Shen Mining Area
By analyzing the collected data as well as the existing research, the membership degree of each evaluation indicator was determined. Subsequently, combined with the weight distribution, we calculated the comprehensive evaluation value Φ of each drill hole under longwall full-seam mining in the Yu-Shen mining area.

Overburden strata system
Based on the drilling data, we obtained the distance from the aquifer to the first-mined coal seam by subtracting the buried depth of the Salawusu Formation aquifer from the buried depth of the first-mined coal seam, and its contour map was plotted using Kriging method, as shown in Figure 9.
Sustainability 2020, 12, x FOR PEER REVIEW 17 of 28 different numerical simulation models of various stratigraphic structures (sand soil bedrock, sand bedrock, bedrock, soil-bedrock) and different mining height (2-12 m). During the numerical simulation, the physical and mechanical parameters of these models were continuously optimized until the simulation results were consistent with the measured results. Then, the formulae for predicting the height of the caved zone and WCFZ in the Yu-Shen mining area were obtained by using linear regression, as shown in Formulas (24) and (25). On the basis of the formulae and the coal seam thickness of different drill holes, we calculated the height of the caved zone and the WCFZ of boreholes in the Yu-Shen mining area, and plotted their contour maps, as shown in Figure 9. The subordinate degree of the distance from aquifer to coal seam of each borehole was then determined.
where Hc is the height of caved zone (m) and M is the mining height (m).  (25) where Hf is the height of WCFZ (m) and M is the mining height (m).
The aquifuge in the mining area mainly consists of the loess in Quaternary Middle Pleistocene Lishi Formation and the red soil in Neogene Pliocene Baode Formation. The Lishi Formation is a loess formation composed of grayish-yellow sandy loam and silty clay, and its thickness ranges from 0 to 109.5 m, with an average of 23.0 m. The Baode Formation is a red soil formation comprised of brownred clay and loam, with a thickness of 30 m on average [4,60]. The contour map of the thickness of aquifuge is shown in Figure 9. 2. Geological structure system.
According to the geological exploration and seismic exploration results, the geological structure of the Yu-Shen mining area is simple, and there are seven normal faults identified in total. The fault characteristics are shown in Table 4.

Groundwater system
Sarawusu Formation aquifer is the primary supply source of water demand for local residents and is decisive for the eco-environment of the study area. Thus, it is the main aquifer for WRCCMA research in the Yu-Shen mining area. The high water-rich regions of the aquifer cover a total area of 504 km 2 , and the salinity in these areas ranges from 0.188 to 0.355 g/L. The total area of medium waterrich areas is 1911 km 2 , with a mineralization degree of 0.170 to 0.561 g/L. The low water-rich regions of the aquifer is 1919 km 2 in area, with mineralization varying from 0.214 to 0.588 g/L. The comprehensive pollution index p of the Sarawusu Formation aquifer is from 0.21 to 0.7, and the groundwater quality is generally II−III grade. For most areas in the Yu-Shen mining area, the water table drop after mining is less than 1 m. However, in the southeast and northeast boundaries of this mining area, the water level declines drastically. The coefficient of recharge from precipitation in Yu-Shen mining area varies from 0.2 to 0.54 [4,57]. The water-rich zoning map of Sarawusu Formation aquifer is shown in Figure 10. The zoning map of mining-inudced groundwater level drop of Sarawusu Formation aquifer is shown in Figure 10. The zoning map of the coefficient of recharge from precipitation is shown in Figure 10.

Many scholars have done a lot of theoretical and practical research on the height of WCFZ in East
China. The empirical formulae for calculating the height of WCFZ, which are suitable for different overburden lithology and dip angles of coal seams, were given when adopting longwall mining method [27][28][29]. However, most of the main mineable coal seams in eastern mining area belong to Permian system, while the Yu-Shen mining area is mainly Jurassic coal seam [4,30]. The hydrogeological and geological engineering conditions in this area are significantly different from those in the eastern mining area [60]. Thus, the traditional prediction formulae of WCFZ height are inapplicable to Yu-Shen mining area. In this paper, field measurement and numerical simulation were employed to determine the heights of the caved zone and the WCFZ. The caved zone height and WCFZ height of Yuyang colliery, Yushuwan colliery, Jinjitan colliery, and Hanglaiwan colliery were obtained by drilling fluid consumption observation assisted with the core engineering geological logging, water pressure test in the borehole, television logging, and geophysical logging. The mining height of these coal mines varies from 3.5 to 5.5 m, which cannot represent the whole Yu-Shen mining area. Therefore, the universal distinct element code (UDEC) was used to establish different numerical simulation models of various stratigraphic structures (sand soil bedrock, sand bedrock, bedrock, soil-bedrock) and different mining height (2-12 m). During the numerical simulation, the physical and mechanical parameters of these models were continuously optimized until the simulation results were consistent with the measured results. Then, the formulae for predicting the height of the caved zone and WCFZ in the Yu-Shen mining area were obtained by using linear regression, as shown in Formulas (24) and (25). On the basis of the formulae and the coal seam thickness of different drill holes, we calculated the height of the caved zone and the WCFZ of boreholes in the Yu-Shen mining area, and plotted their contour maps, as shown in Figure 9. The subordinate degree of the distance from aquifer to coal seam of each borehole was then determined.
where H c is the height of caved zone (m) and M is the mining height (m).
where H f is the height of WCFZ (m) and M is the mining height (m). The aquifuge in the mining area mainly consists of the loess in Quaternary Middle Pleistocene Lishi Formation and the red soil in Neogene Pliocene Baode Formation. The Lishi Formation is a loess formation composed of grayish-yellow sandy loam and silty clay, and its thickness ranges from 0 to 109.5 m, with an average of 23.0 m. The Baode Formation is a red soil formation comprised of brown-red clay and loam, with a thickness of 30 m on average [4,60]. The contour map of the thickness of aquifuge is shown in Figure 9.
According to the geological exploration and seismic exploration results, the geological structure of the Yu-Shen mining area is simple, and there are seven normal faults identified in total. The fault characteristics are shown in Table 4.

Groundwater system
Sarawusu Formation aquifer is the primary supply source of water demand for local residents and is decisive for the eco-environment of the study area. Thus, it is the main aquifer for WRCCMA research in the Yu-Shen mining area. The high water-rich regions of the aquifer cover a total area of 504 km 2 , and the salinity in these areas ranges from 0.188 to 0.355 g/L. The total area of medium water-rich areas is 1911 km 2 , with a mineralization degree of 0.170 to 0.561 g/L. The low water-rich regions of the aquifer is 1919 km 2 in area, with mineralization varying from 0.214 to 0.588 g/L. The comprehensive pollution index p of the Sarawusu Formation aquifer is from 0.21 to 0.7, and the groundwater quality is generally II−III grade. For most areas in the Yu-Shen mining area, the water table drop after mining is less than 1 m. However, in the southeast and northeast boundaries of this mining area, the water level declines drastically. The coefficient of recharge from precipitation in Yu-Shen mining area varies from 0.2 to 0.54 [4,57]. The water-rich zoning map of Sarawusu Formation aquifer is shown in Figure 10. The zoning map of mining-inudced groundwater level drop of Sarawusu Formation aquifer is shown in Figure 10. The zoning map of the coefficient of recharge from precipitation is shown in Figure 10.

Coal mining system
The Jurassic Yan'an Formation is the coal-bearing strata in the mining area. There are five main mineable coal seams, i.e., 1 −2 , 2 −2 , 3 −1 , 4 −2 , and 5 −2 coal seams. These seams are generally flat with small dip angles of 1 • to 3 • , while their thicknesses vary significantly, ranging from 0 to 12 m. The contour map of the first-mined coal seam thickness is shown in Figure 11. Through on-site investigation, the author found that when the mining height is 4 m, 5 m, 6 m, and 7 m, the corresponding maximum subsidence values are 2.2 m, 2.8 m, 3.3 m, and 3.8 m, respectively, in the Yu-Shen mining area. The formula of maximum subsidence value was obtained by linear fitting, as shown in Formula (26).
Based on the formula and the coal seam thickness of different boreholes, we calculated the maximum subsidence values of different drill holes under longwall full-seam mining in the Yu-Shen mining area, and the contour map of them was plotted, as shown in Figure 11.

Coal mining system
The Jurassic Yan'an Formation is the coal-bearing strata in the mining area. There are five main mineable coal seams, i.e., 1 −2 , 2 −2 , 3 −1 , 4 −2 , and 5 −2 coal seams. These seams are generally flat with small dip angles of 1° to 3°, while their thicknesses vary significantly, ranging from 0 to 12 m. The contour map of the first-mined coal seam thickness is shown in Figure 11. Through on-site investigation, the author found that when the mining height is 4 m, 5 m, 6 m, and 7 m, the corresponding maximum subsidence values are 2.2 m, 2.8 m, 3.3 m, and 3.8 m, respectively, in the Yu-Shen mining area. The formula of maximum subsidence value was obtained by linear fitting, as shown in formula (26).
Based on the formula and the coal seam thickness of different boreholes, we calculated the maximum subsidence values of different drill holes under longwall full-seam mining in the Yu-Shen mining area, and the contour map of them was plotted, as shown in Figure 11.

Ecological system
The vegetation coverage in Yulin city is 45.82% [62]. The total water resources of Yulin city in 2018 was 2.432 billion m 3 , 90.4% of that in 2017, which is 2.69 billion m 3 [63].

Social system
The water utilization rate of Yulin city from 2014 to 2018, published by the Shaanxi Provincial Water Resources Department, is shown in Table 5. It can be seen that the water utilization rate of Yu-Shen mining area is 34.04% on average [63]. In addition, in light of the distribution of water resources per capita from 2014 to 2018 released by the National Bureau of Statistics [64], the water resources per capita of Yu-Shen mining area are 1000-2000 m 3 , revealing that it is a moderate water scarcity area.

Results of WRCCMA Prediction
Based on the comprehensive evaluation value Φ of the WRCCMA of each borehole, the Kriging method was employed to draw the contour map of the Φ under longwall full-seam mining in the Yu-Shen mining area, as shown in Figure 12. The evaluation results in this paper are consistent with other evaluation methods, indicating this model is scientific, reasonable, and reliable [60]. In addition, this evaluation system can not only accurately predict and evaluate the WRCCMA in the Yu-Shen mining area but also can be employed in other arid and semi-arid regions in Northwest of China, since they share the arid and semi-arid climate, fragile ecology, and water scarcity in common. Based on the prediction results of WRCCMA grades of different collieries (boreholes) in the Yu-Shen mining area, three countermeasures, namely, mining methods selection, mine water reutilization, and water-resisting layer reconstruction, were proposed to achieve WCCM.

Mining methods selection
For areas where there is no water distribution, or the WRCCMA is carrying surplus, the As a whole, the comprehensive evaluation values of drill holes in the northwest area are generally larger than those in the southeast. The areas of grade I are geographically dispersed, and are distributed in the central west, southwest, and northeast area in Yu-Shen mining area, where the longwall full-seam mining can realize the coordinated development among the high-efficient mining, water resources preservation, and ecology protection. In the western border of the study area, the WRCCMA belongs to grade II and grade III, which are generally lower than that in the adjacent eastern boreholes. The main reason is that the areas near the west border are eolian sand landform and are characterized by low vegetation coverage and serious land desertification. The areas of grade II and grade III of WRCCMA are distributed in the middle part, and a small part of northeast border and WCCM methods should be used in these areas to maintain the WRCCMA. The areas of grade IV are distributed in the northeast border and southeast border of the study area. The areas of grade V are distributed in the southeast of Yu-Shen coal area, in which the coal mining should be strictly prohibited because the WRCCMA is severly over-carrying.
The evaluation results in this paper are consistent with other evaluation methods, indicating this model is scientific, reasonable, and reliable [60]. In addition, this evaluation system can not only accurately predict and evaluate the WRCCMA in the Yu-Shen mining area but also can be employed in other arid and semi-arid regions in Northwest of China, since they share the arid and semi-arid climate, fragile ecology, and water scarcity in common.
Based on the prediction results of WRCCMA grades of different collieries (boreholes) in the Yu-Shen mining area, three countermeasures, namely, mining methods selection, mine water reutilization, and water-resisting layer reconstruction, were proposed to achieve WCCM.

Mining methods selection
For areas where there is no water distribution, or the WRCCMA is carrying surplus, the traditional longwall full-seam mining method can be used to extract coal resources. However, when the WRCCMA is capable of carrying, the coal mines should adopt harmonious mining or slice mining methods. By adjusting the mining sequence of the panels, these two methods can reduce and counteract the movement and deformation of overlying strata. Hence, the WCFZ height will be restricted, and the impact of coal mining on water bodies will decrease. Moreover, coal mines with medium WRCCMA should employ partial mining methods, including height-restricted mining, room and pillar mining, and narrow strip mining, since the water flowing fractures generated by traditional longwall full-seam mining in these collieries are prone to reaching the overlying aquifer, resulting in water level drop and groundwater loss. These three WCCM methods can prevent groundwater from seeping and thus maintain the WRCCMA by discarding part of coal body to support the overlying strata and restrict their movement. If the WRCCMA is severely over-carrying, mining activity should be completely banned in order to preserve the precious water resources.
For collieries whose WRCCMA is slightly over-carrying, backfill mining must be employed, because the WRCCMA has reached the warning value. This method substitutes extracted coal body with filling bodies and allows all coal resources in an area to be extracted, and is recognized as an ideal approach to protecting water from mining-induced damage. However, traditional paste filling material and high water swelling material are confronted with several problems, such as insufficient sources, small filling scales, and high filling costs, and it is difficult for them to be popularized and applied in a larger range. On the other hand, coal mine solid waste was no longer needed in the process of production and living activities and was therefore accumulated on the ground for a long time. The harmful contaminants contained in the solid waste are prone to infiltrating into the shallow water by leaching process, contributing to the water quality deteriorating and the WRCCMA declining. In addition, as there is no coal body to support the overburden strata after mining, the water conductive fractures can easily reach the aquifer and thus lead to water table lowering and water bursting in mine, which not only jeopardizes the mine safety, but also severely affects the domestic water of local residents. Therefore, backfilling the solid waste into the goaf to restrict the movement of overlying strata and the development of the water flowing fissures can not only preserve shallow water resources from being polluted, but prevent water seepage. As an effective method of disposing of mining wastes, solid waste backfill mining method allows longwall full-seam mining under water bodies. It will make significant contributions to the harmless treatment of solid waste and sustainable development of the coal area.

Mine water reutilization
Seventy percent of coal mines in arid and semi-arid areas of Northwest China are short of water, over 40% of which suffer from serious water shortage. Collieries throughout China discharge about 4.2 billion m 3 of mine water every year, while the reutilization rate is only 40% [58][59][60][61]. In the process of coal mining, aquifer dewatering has caused large-scale water level drop and water loss. Also, due to the large evaporative capacity in arid and semi-arid areas, the discharged mine water on the surface can quickly be evaporated and lost. If the WRCCMA of the colliery is carrying surplus or capable of carrying, it should purify the polluted mine water while making a profit. The treated and purified water can provide a source of water for coal mines and households. The over-exploitation of water resources can be avoided, and thereby the WRCCMA can be maintained.
In addition, mining areas where the WRCCMA is moderate carrying or slightly over-carrying are usually the ones with a short distance from the aquifer to the coal seam and without effective water-resistant strata (aquifuge). Once the coal seam is mined, the water flowing fractures will connect the overlying aquifer and the goaf, which inevitably leads to a large quantity of groundwater infiltration into gob, meaning the scale of mine drainage treatment can no longer meet the requirements. Therefore, it is necessary to construct underground reservoirs to store and utilize mine water. The subsurface reservoir can make full use of the function of environmental self-purification of rock mass in goaf. It is unnecessary for the mine water to be discharged to the ground, thus reducing the evaporation loss. In recent years, more than 32 underground coal mine reservoirs have been built in Shen-Dong mining area, with a water storage capacity of more than 31 million m 3 . The water stored in a reservoir can be used for colliery production and local residents' living, and excessive abstraction of groundwater can be avoided.

Water-resisting layer reconstruction
The mining-induced water conductive fissures of slightly over-carrying colliery are likely to break the balance of original stress of the aquifer and change the hydrogeological conditions and the recharge, runoff, discharge of the groundwater. More seriously, the groundwater entrained with sand or mud flows into the mined-out area, causing water bursting and sand gushing in mines and threatening the mine safety. Under this background, the impermeable layer reconstruction was proposed to reduce the permeability of rock mass by grout injection in the layer under the overlying aquifer, so as to block the water flowing channels between the aquifer and the goaf. This technology can enable the aquifer to be recharged and to restore water, realizing water level recovery and thereby helping to maintain the WRCCMA. At the moment, this technology is utilized more frequently to prevent the water inrush from the floor in underground coal mining. The impermeable layer reconstruction in overlying strata of coal seam has been applied in Dashui colliery and Shangwan colliery, and good economic and environmental benefits have been achieved [2].

Discussion
The prediction results of the model have good agreement with the ones of other studies. Overall, the grade of WRCCMA in the Yu-Shen mining area decreases from west to east. The main reasons are that the areas near the west boundary of the study area feature high or medium water abundance and thick bedrock, while the east areas are usually covered by thin bedrock and water yield property in these areas is usually low or even without Salawusu Formation aquifer distribution. Other similar research usually treats the mining area as a whole and evaluates the WRCC for the entire coal area.
The evaluation results of mining area are usually used as the instructions and guidance for the mining planning and scientific mining scale [2,11]. By contrast, we evaluated and predicted the WRCC based on each borehole distributed in the whole mining area. Therefore, the prediction results of our article can provide instructions for each colliery with different grades of WRCC. Compared to other studies, our findings can better instruct coal production of each colliery and help coal mines to maintain their WRCC, making significant contributions to the sustainable and green development of the mining area. To sum up, our findings have advantages over findings of other studies in terms of instructing each colliery's production in the premise of preserving water resources.
WCFZ height is decisive to achieve WCCM and thus is an important factor affecting the WRCCMA. We obtained the formula for predicting the WCFZ height of four different stratigraphic structures (sand-soil-bedrock, sand-bedrock, bedrock, and soil-bedrock) by numerical simulation and field measurements, which is only related to mining height. However, the coefficient in the proposed formula also takes the influence of other factors on WCFZ height into consideration. Therefore, it is reasonable to use mining height as the only variable in the equation. The equation can contribute to accurately evaluating the WRCCMA of the Yu-Shen mining area.
The Yu-Shen mining area is a part of Yulin city. Therefore, the membership degree of some evaluation indicators of WRCCMA, including vegetation coverage, total water resources, water utilization rate, and water resources per capita, were determined by the data of Yulin city, since there are no relevant collected data of different boreholes in the Yu-Shen mining area up to now. However, the membership degrees of different collieries are not always the same. Further field investigation is thus needed to make the evaluation of the WRCCMA more accurate, so as to better instruct WCCM.
Although these three countermeasures offered to maintain the WRCCMA can provide guidance for WCCM, it is necessary for planners and engineers to take into account the specific engineering and geological conditions of collieries during implementation. Taking mining methods, for example, coordinated mining should consider the layout of the coal face and mining sequence to partly reduce and offset the tensile and compressional strains of overlying strata and thereby the development of the water flowing fissures. For narrow strip mining, the width of the narrow strip and strip coal pillar need to be optimized to maximize the coal extraction in the premise of water conservation. In addition, the strip coal pillars left unmined may fail owing to weathering after years of mining, which suggests that narrow strip mining is only suitable under certain geological conditions. Considering the current technological and economic conditions and equipment available, it is suitable for this WCCM method to be adopted at low-production coal mines. Furthermore, for room and pillar mining, in order to improve the recovery ratio of coal resources in the premise of WCCM, the sizes of room and pillar should be analyzed in detail. Both of these two mining methods can be particularly used in collieries where the shapes coal seams are irregular and longwall panels are unsuitable to arrange. Besides, for backfill mining, filling rate, filling materials, and filling pipelines need to be determined and designed based on the different conditions. In order to preserve water resources while maximizing economic benefits, two or more WCCM methods can be employed together in a colliery. For example, after narrow strip mining, the backfill mining method can be utilized to extract coal pillars left unmined, so as to maintain the WRCCMA while improving the recovery rate. For specific conditions, the construction of underground coal mine reservoir requires the determination of water source prediction, site selection, storage capacity design, dam construction, and water quality guarantee. As for impermeable layer reconstruction, it is necessary to select the specific grouting materials, grouting parameters, and grouting pressure, in terms of the physical and mechanical properties of overlying strata [65].

Conclusions
The following conclusions are drawn from the research: 1. Based on the AHP-fuzzy comprehensive evaluation method, a three-level AHP model, with 6 sub-factors and 15 tertiary factors, was developed to identify the influencing factors of the WRCCMA Sustainability 2020, 12, 7782 24 of 27 and determine the weight distribution. The subordinate functions of these influencing factors were constructed. According to the weights distribution, among the secondary indicators, the coal mining system is the most important one, with a weight of 0.3319. The mining method and buried depth of aquifer are the two most important tertiary indicators, with weights of 0.1791 and 0.1550, respectively.
2. The fluid-solid coupling module of FLAC 3D finite element software was employed to simulate the seepage of the groundwater in Yubujie coal mine, and its water level has dropped 1.32 m. The prediction model of WRCCMA was applied in the Yubujie colliery, and the comprehensive evaluation value Φ was calculated to be 0.84, indicating that the WRCCMA is capable of carrying. The ground deformation and the vegetation coverage alteration after mining in the colliery were analyzed through field investigation, and the validity of this evaluation model was verified.
3. The comprehensive evaluation value Φ of each borehole was obtained in Yu-Shen mining area. The Kriging method was employed to plot the contour map of Φ in this study area. In view of different grades of WRCCMA of collieries, three countermeasures, i.e., selecting reasonable mining methods, reutilizing the mine water, and reconstructing the water-resisting layer, have been provided to help them achieve WCCM and maintain the WRCCMA. Reasonable mining methods can control the deformation and break of overlying strata and thus can prevent groundwater from infiltrating into the mined-out area. Mine water reutilization can provide a source of water for coal mines and households and avoid the over-exploitation of water resources. The water-resisting layer reconstruction can reduce the permeability of rock mass by grout injection in the layer under the overlying aquifer, so as to block the water flowing channels between the aquifer and the goaf. These three countermeasures can not only be used in the Yu-Shen mining area, but also can be popularized and applied in a larger range, such as Xinjiang, Ningxia, and Inner Mongolia coal production base in the Northwest of China. More broadly, these three countermeasures can be used in non-coal mines under water bodies to protect the previous water resources. There are severe conflicts between coal mining and water preservation. Only by predicting the WRCCMA accurately can water preservation and coal resources extraction develop in a coordinated way. Without the precise evaluation of WRCCMA, water loss on a larger scale and a water crisis will occur in the coming years. This article combines fuzzy mathematics with mining engineering to construct a fuzzy comprehensive evaluation model to predict the WRCCMA. The research results can provide a theoretical basis and practical reference for coal mining planning in ecologically fragile arid and semi-arid mining areas. It is conducive to promoting the coordinated development of coal resources extraction, water resources preservation and ecology protection, contributing to the realization of green mining and sustainable development of the mining area in Northwest China.