Abstract

Many deep mining mines in southwestern Shandong Province of China are covered with thick loose layers. When mining near the loose layers, there is a risk of water and sand inrush, which threatens the personal safety of miners. The prediction of sudden water and sand inrush is difficult due to the comprehensive influence of many factors, and the influencing factors are fuzzy and random. To solve this problem, in this paper, a new risk assessment method of water and sand inrush based on comprehensive weight and cloud model was proposed. Seven factors are selected as indexes: the aquifer thickness, the thickness ratio of sand layer to clay layer, the thickness of bottom clay layer, the coal seam thickness, the percentage of core recovery, the geological structure, and the bedrock thickness. The assessment index system is established, and the index is divided into three grades. A comprehensive weighting method, which combines analytic hierarchy process (AHP), entropy weight method (EWM), and minimum entropy principle, is used to reasonably assign the weight of index. Based on the cloud generator equation, the membership function is obtained. The assessment result of the assessment object is obtained by combining the membership degree and the weight of index. The comprehensive weight-cloud model assessment method is applied to the risk assessment of water and sand inrush in the 6311-2 working face in the sixth mining area of Baodian Coal Mine. According to the assessment results, the following conclusions can be drawn: (1) the bedrock thickness and the coal seam thickness are the main factors of water and sand inrush under loose layer mining; (2) the assessment results obtained by the comprehensive weight-cloud model method are consistent with the actual situation. The assessment method can provide scientific reference for the safe mining under the thick loose layer in the deep mines of southwest Shandong.

1. Introduction

Water and sand inrush is a kind of mine geological disaster that water and sand mixed fluid with high sand content bursts into underground working face and causes property damage and casualties [1]. Some coalfields in North China are covered with thick loose layers, especially in deep mining mines in southwest Shandong Province. When mining near the loose layers, the upper water-rich sand layer is prone to water and sand inrush under the disturbance of mining activities, which affects the normal production of the mine and causes casualties [2, 3]. In order to reduce the occurrence of water and sand inrush disasters and take timely and effective measures, it is necessary to put forward a more accurate assessment method of water and sand inrush disasters under loose layer mining [4, 5].

Experts and scholars studied the problem of water and sand inrush by various methods [613]. Zhong et al. used software PFC3D and software GID to simulate the whole process of water and sand inrush in precast ideal fracture with different opening widths and dip angles in overlying rock strata. Their analysis shows that the opening widths and dip angles of fracture change the flow patterns of water-sand flow inrush and have great impact on the contact force of the fracture channel, flow velocity of water, and the time of mixed water-sand flow [14]. Zhao et al. studied the overlying stratum fracture development and distribution characteristics of water-sand inrush channel through the simulation experiment. They divided the development process of water-inrush channels into three stages: the stage of gradual development, the stage of penetration linking channel formation, and the stage of water-sand inrush, and divided the area of the overlying rock fracture and water-sand inrush into three sections: the zone of overburden fracture gradual development, the zone of water-sand intrusion, and the zone of water-sand-intrusion blocking [15]. Peng et al. made a comprehensive analysis on the mechanism of water and sand inrush disaster from many aspects such as channel, water source, water storage space, power source, and geological structure. They found out that the cause of water and sand inrush disaster under thick overlying bedrock is that the water flowing fractured zone generated by mining causes water to enter the separated cavity between rock formations. The water in the separated cavity penetrates into the loose geological body and breaks into the working face instantly along the concentrated channel generated by cutting the working face, resulting in the occurrence of water and sand inrush disaster [16]. Zhang et al. established the mechanics model of sand inrush in fractures and analyzed the limit equilibrium condition of water and sand inrush in fractures. Through experiments, they quantitatively analyzed the characteristics of water and sand migration and correlation changes of physical parameters in different stages of water and sand inrush and divided the whole process of water and sand inrush in fractures into four stages, namely, start-up stage, continuous outburst stage, silt blockage stage, and outburst equilibrium stage [17]. Ma et al. established a water-sediment flow resistance model in fractures based on the two-phase flow theory and verified it through laboratory-scale test [18]. Using the LBM-DEM coupling simulation method, Pu et al. studied the problem of the water and sand two-phase migration in the single-fracture opening channel model. They compared the changes of section flow rate and sand inrush rate under different boundary pressures, fracture opening widths, and sand layer thickness [19].

As mentioned above, researchers have made many achievements in the mechanism of water and sand inrush. However, the prediction of sudden water and sand inrush is difficult due to the comprehensive influence of many factors, and the influencing factors are fuzzy and random. To solve this problem, a risk assessment method of water and sand inrush is proposed based on comprehensive weight and cloud model in this paper. The membership degree function transformed from the normal cloud generator equation is used to calculate the membership degree of the index. Combine analytic hierarchy process, entropy weight method, and minimum entropy principle to calculate the comprehensive weight of the index. Based on the membership degree and the comprehensive weight, the risk of water and sand inrush of the assessment object is evaluated, hoping it can provide new ideas and methods for the prevention and control of water and sand inrush disasters. The assessment process is shown in Figure 1.

2. Overview of the Study Area

The sixth mining area of Baodian Coal Mine is selected as the study area. The range and location of boreholes near the sixth mining area are shown in Figure 2. Baodian Coal Mine is located in Yanzhou District, Jining City, Shandong Province. The sixth mining area is located in the west of Baodian Coal Mine. The structure is controlled by Yanzhou syncline, and the axial direction is NEE, inclining to the northeast. There is a small south lake syncline in the north of the sixth mining area. The southern development range was Baochang anticline. Fault strike is mostly northeast.

The sixth mining area is a fully concealed North China Carboniferous Permian coalfield. The strata from old to new are Ordovician (O2), Carboniferous (C), Permian (P), Jurassic (J3), and Quaternary (Q). The following is a detailed description: (1)Middle and lower Ordovician (O2, 1): it is the basement of coal measure strata, which is composed of gray and gray-white limestone(2)Carboniferous (C): the Taiyuan formation of upper Carboniferous is composed of dark gray-grayish black mudstone, bauxite mudstone, siltstone, and medium-coarse sandstone, with 0-11 layers of limestone. Among them, the thickness of the tenth lower limestone and the third limestone is large and the horizon is stable, which is the auxiliary marker layer of the sixth mining area(3)Permian (P): Shanxi formation is the main coal-bearing strata in the sixth mining area. It is thick in the north but thin in the south. It is composed of gray-white medium, coarse sandstone, gray siltstone, mudstone, bauxite mudstone, and coal seam. Among them, No. 3 coal seam is the main minable coal seam, which has complete contact with the underlying strata(4)Jurassic (J3): the upper member is gray-green, purple-gray medium-fine sandstone. The middle member is loose red sandstone. The next section is brownish-red siltstone. It is distributed within a very small range in the eastern part of the sixth mining area. It is in angular unconformity contact with underlying coal measures(5)Quaternary (Q): thin in the east and thick in the west, thin in the south, and thick in the north. It is composed of sandy clay, clay sand, clay layer, and medium and coarse sand layers

In all strata, the main coal-bearing area is the Carboniferous strata and Permian strata of Shanxi and Taiyuan formation, which belongs to the type of coal-bearing rock series in North China. The main coal seam is the No. 3 coal seam, with a thickness of about 7.86~10.02 m and an average thickness of about 9.00 m. The thickness of the coal seam is stable, and the buried depth of coal seam is about 200~390 m.

The main aquifers affecting the production of the No. 3 coal seam in the sixth mining area from top to bottom are the gravel aquifer in the Quaternary upper group, the gravel aquifer in the Quaternary lower group, and the sandstone aquifer at the roof and floor of the No. 3 coal seam. Among them, the direct water-filled aquifer of coal seam mining is the sandstone aquifer at the roof and floor of the No. 3 coal seam and the gravel aquifer in the Quaternary lower group, and the indirect water-filled aquifer is the gravel aquifer in the Quaternary lower group (when the sandstone aquifer at the roof of the No. 3 coal seam is the direct water-filled aquifer). Except for the gravel aquifer in the Quaternary upper group, the remaining aquifers are mainly static reserves, and the recharge, runoff, and discharge conditions are poor. With the development of mining activities, the water level of the sand layer in the Quaternary lower group decreased slowly year by year, and the water level of the sandstone at the roof of the No. 3 coal seam decreased significantly.

3. Assessment Methods

3.1. Analytic Hierarchy Process

American operational researcher Saaty put forward the famous analytic hierarchy process (AHP) in the early 1970s. The analytic hierarchy process is a decision-making method which decomposes the elements related to decision-making into objective, criterion, plan, and other levels and, on this basis, makes qualitative and quantitative analysis [20].

As a weight determination method, the analytic hierarchy process is commonly used in the field of mine water disasters, such as water abundance assessment and floor water inrush risk assessment [2123]. In this paper, the improved analytic hierarchy process with three scales is used to calculate the weight.

The traditional AHP needs to check the consistency of the judgment matrix. In this paper, the improved three-scale AHP is used for weight calculation, and the traditional AHP is optimized by using the properties of the optimal transfer matrix, so that it naturally satisfies the consistency. It can greatly reduce the number of iterations and make the subjective factors analytic, thereby reducing the system error.

Supposing there are lower-level indexes under a certain upper-level index, the importance of each index at the same level is compared according to expert consultation, and the comparison matrix is established through the following equation:

Since matrix satisfies and , it is an antisymmetric matrix; then, according to the principle of optimal transfer matrix, the optimal transfer matrix of matrix should conform to Equation (2).

According to the Equations (2) and (3), the judgment matrix is obtained.

Due to the properties of the optimal transfer matrix, no consistency check is required. The equation of the weight of the index is as follows:

Finally, according to Equation (4), the weight of the plan level to the criterion level and the weight of the criterion level to the objective level are calculated, respectively, and then, the weight vector of the plan level to the objective level is obtained, where is the index number in the plan level.

3.2. Entropy Weight Method

As a weight calculation method, entropy weight method (EWM) has been well applied in the weight calculation of multifactor indexes [2427].

Assuming that there are samples and each sample has indexes, the original matrix can be constructed, where represents the data of index of the sample .

Then, according to the original matrix , the normalized matrix is calculated.

For the positive indexes that the greater the better indexes, the calculation equation is as follows:

For negative indexes that the smaller the better indexes, the calculation equation is as follows:

According the normalized matrix , calculate the proportion of the index data of the sample:

Then, calculate the entropy of all index, and the calculation equation for the entropy value of the index is

Calculate the entropy weight of all index, and the calculation equation for the entropy weight of the index is

Finally, get the weight vector .

3.3. Comprehensive Weight

Each weight calculation method has its own scope of application, and sometimes, it is often necessary to use a variety of methods to measure the weight of the same data, so that the comprehensive weight has higher performance and can reflect the real characteristics of the data. According to the minimum entropy principle, this paper processes the weight vector determined by the analytic hierarchy process and the weight vector determined by the entropy weight method to determine the comprehensive weight vector [28, 29]. The calculation process is as follows:

3.4. Cloud Model

Cloud model is an uncertain cognitive model based on fuzzy set theory and probability concept, which was proposed by Liu et al. [27]. The cloud model can be used to deal with the uncertain conversion between qualitative concepts and quantitative description and has been widely used in algorithm improvement, simulation, risk assessment, geological prediction, excavation, and other fields [3033]. In the conversion process from quantitative data (influencing factors data) to qualitative concepts (risk grade), the cloud model can better handle the effects of randomness and ambiguity, thus making the evaluation results more scientific and accurate.

Normal cloud is an important cloud model based on normal distribution and Gaussian membership function. Since the expected value curves of influencing factors in natural science are mostly normal distribution or seminormal distribution [34], the normal cloud model is used in this paper to evaluate the risk of water and sand inrush.

Supposing the set is a domain, the qualitative concept on the domain is defined as . For any belonging to , there exists a random number belonging to . The set of is called the membership degree of belonging to , if satisfies

If satisfies and , the distribution of on is called a normal cloud, and each is called a cloud drop. Ex, En, and He are the numerical characteristics of a qualitative concept, where Ex represents expectation, En represents entropy, and He represents hyperentropy. If the three numerical characteristics of the qualitative concept are known, the normal cloud generator can be used to generate the normal cloud. The process is as follows:

Step 1. Generate a normal random number En with an expected value of En and a standard deviation of He.

Step 2. Generate a normal random number with an expected value of Ex and a standard deviation of .

Step 3. Calculate through Equation (12), and a cloud drop with a membership degree of for the qualitative concept is generated.

Step 4. Repeat Step 1–Step 3 until the number of cloud drops meets the requirements.

The calculation of the numerical characteristics ( and ) and of the index belonging to the risk grade is as follows.

Assuming that the upper and lower boundary values of the index belonging to the risk grade are and , then

Since the boundary value is from one grade to another and should belong to both grades [32], so

The size of is determined according to the fuzziness and randomness of the specific case, and the value is about 0.1 times of [35].

Some scholars use the cloud model to obtain the membership degree by using the cloud generator to randomly generate cloud drops and then obtain the average membership degree [33]. The membership degree obtained by this method has a certain degree of volatility, resulting in the same calculation process which may not be able to obtain the same assessment results. In order to obtain a stable membership degree, the membership function is obtained based on Equation (12). The membership function is shown below: where represents the value of the index and represents the membership degree of the index belonging to the risk grade .

The membership degree matrix of the assessment object is obtained, where is the number of indexes, and is the number of risk grades.

In order to find out the membership degree of the assessment object to a risk grade, it is necessary to multiply the membership degree of index of the assessment object corresponding to the risk grade by the index weight and add them, so as to obtain the membership degree of the assessment object to the risk grade.

The comprehensive weight vector is combined with the membership matrix of the assessment object to obtain the matrix .

In Equation (17), represents the membership of the assessment object to the risk grade . Then, according to the principle of maximum membership degree, the risk grade corresponding to the maximum membership degree is the water inrush risk grade of the assessment object.

4. Preparation of Assessment

4.1. Index Selection

The occurrence of water and sand inrush in the mining under the loose layer depends on the combined effect of various influencing factors. According to the hydrogeological data in the sixth mining area of Baodian Coal Mine, this paper selects seven influencing factors, namely, the aquifer thickness (), the thickness ratio of the sand layer to clay layer (), the thickness of the bottom clay layer (), the coal seam thickness (), the percentage of core recovery (), the geological structure (), and the bedrock thickness (), as index. These factors can be categorized into the characteristics of loose layer and the characteristics of rock layer. The hierarchical structure system of water and sand inrush index for mining under loose layer in the mining under the loose layer is shown in Figure 3. (1)The aquifer thickness (). Generally speaking, the aquifer in the loose layer is mainly sand layer. Thicker aquifer can store more groundwater, and under the influence of mining, there is a greater possibility of water and sand inrush [36](2)The thickness ratio of the sand layer to clay layer (). The sand layer with fractures has strong water storage and water conductivity, while the clay layer has certain water and sand resistance. The thickness ratio of the sand layer to clay layer in the loose layer of the Quaternary lower group determines the risk of water and sand inrush(3)The thickness of the bottom clay layer (). The clay layer at the bottom of the loose layer is a powerful barrier that directly hinders the downward seepage and inrush of water and sand in the upper aquifer. The thicker the clay layer at the bottom is, the more likely it will reduce the possibility of water and sand inrush in the exploitation of coal resources [36](4)The coal seam thickness (). The height of the roof fall zone and fracture zone caused by coal seam mining is related to the cumulative mining thickness of the coal seam. Generally speaking, the greater the cumulative mining thickness of the coal seam, the greater the height of the caving zone and fracture zone of the roof. In mining, it is necessary to set sand-prevention coal and rock pillars to avoid excessive water and sand inrush due to the excessive height of the caving zone and fracture zone. Since No. 3 coal is the main mineable coal seam in the sixth mining area, this paper uses the thickness of the No. 3 coal seam instead of the cumulative mining thickness as the index [37](5)The percentage of core recovery (). The integrity of the core taken during drilling is related to the degree of rock fragmentation. The core of the bedrock is relatively complete, indicating that the bedrock has a low degree of fragmentation and is an effective water-blocking layer. The possibility of water and sand inrush is low when mining [38](6)The geological structure (). The development degree of geological structure can be expressed by density of faults, fault drop, density of joints, etc. The area with developed geological structures has high risk of water and sand inrush. Based on the actual mining experience and geological data, the degree of geological structure development is divided, and the equation is as follows [39]:(7)The bedrock thickness (). The bedrock is thick and stable, the fracture zone cannot develop to the aquifer, and the risk of water and sand inrush is low. If the bedrock is thin or missing, the fracture zone develops to the aquifer, and the risk of water and sand inrush is high [33, 40]

4.2. Risk Grade

The risk of water and sand inrush is divided into three grades, namely, low risk (I), medium risk (II), and high risk (III). The corresponding situation of grade I is that the aquifer at the bottom of the loose layer has little influence on the mining of the working face, and the water and sand inrush will not occur in the mining process. The corresponding situation of grade II is that the aquifer at the bottom of the loose layer has a certain influence on the mining of the working face. For example, the roof of the working face often shows the phenomenon of water leaching, and the water inflow of the working face changes greatly; mining process, sudden water, and sand inrush may occur. The corresponding situation of grade III is that the aquifer at the bottom of the loose layer has a great impact on the mining of the working face. When the roof comes to pressure, the water inflow of the working face changes greatly, and there is a great possibility of sudden water and sand inrush in the mining process. According to the engineering experience, each index is divided into intervals according to the risk grade, as shown in Table 1.

4.3. Weight Calculation

Based on the analysis and statistics of hydrogeological data and borehole data in the study area, 28 borehole data are collected as samples, as shown in Table 2, and the borehole locations are shown in Figure 2.

After consulting and analysis, this paper argues that for the occurrence of water and sand inrush, in the criterion level, the characteristics of the rock layer are greater than the characteristics of loose layer. For the characteristics of loose layer, the index weights in descending order are the aquifer thickness (), the thickness of the bottom clay layer (), and the thickness ratio of the sand layer to clay layer (). For the characteristics of the rock layer, the index weights in descending order are the bedrock thickness (), the coal seam thickness (), the percentage of core recovery (), and the geological structure (). The comparison matrix of criterion level and the comparison matrices and of plan level are obtained as follows:

The weight vector is calculated according to Equation (4). According to the data in Table 2, the weight vector is determined by the entropy weight method, and then, the comprehensive weight vector is obtained by Equation (10). Analyzing the weight vector , the thickness of the bottom clay layer (), the coal seam thickness (), and the bedrock thickness () have a greater impact on the risk of water and sand inrush. Among the weight vector , the aquifer thickness (), the coal seam thickness (), and the bedrock thickness () are larger. Combining the weights determined by the three methods, it can be concluded that the coal seam thickness () and the bedrock thickness () are the main influencing factors of water and sand inrush. The index weights of the three methods are shown in Table 3 and Figure 4.

4.4. Numerical Characteristics

Based on the risk grade range and borehole sample data (Table 2), by using Equations (13) and (14), the numerical characteristics (Ex, En, and He) of each index belonging to risk grades I, II, and III are determined, as shown in Table 4.

According to the normal cloud generation process and the numerical characteristics, the normal cloud of each index belonging to each risk grade is generated (Figure 5). The number of the cloud drops for each normal cloud is 1000. The normal cloud can represent the distribution of the membership degree of the index belonging to a certain risk grade.

As can be seen from Figure 5, Ex determines the position of the center point of the risk grade normal cloud; En determines the range of the risk grade normal cloud. The larger the En, the larger the risk level normal cloud range. He determines the discreteness of the normal cloud of risk grades. When the ratio of He to En is small, the distribution of the normal cloud tends to a curve with a normal distribution. When the ratio of He to En is large, the dispersion of the normal cloud is large.

5. Verification via the Application

The 6311-2 working face is located in the west of Baodian Coal Mine, and the south of the working face is close to the outcrop area of the No. 3 coal seam aeolian oxidation zone in the sixth mining area. The No. 3 coal seam is mined at the working face, the thickness of the coal seam is about 8.12~9.16 m, and the average is about 8.64 m. The Quaternary lower group is composed of gray-green, gray-yellow, and gray-white clay; clay-bearing gravel; and sand. The main aquifers of the lower group are clay gravel and gravel layers.

The coal seam of the working face is a monoclinic structure and belongs to the north wing of the Baojiachang anticline. Small secondary wide and gentle folds are developed in the working face. The maximum water inflow of the 6311-2 working face during mining is 24 m3/h, and it is 18.9 m3/h under normal conditions. The water inflow of the working face is basically the water inflow after mining and the water inflow during production. The working face did not show excessive water inrush locally and at intervals, and no water and sand inrush disaster occurred.

According to the corresponding geological report, the index data of the working face was determined (Table 5), by using Equations (15)–(17) to calculate the membership degree of the working face belonging to each risk grade according to the index data of the working face. According to the calculation results (Table 6), the maximum membership degree of the 6311-2 working face is 0.3610, and the water and sand inrush risk grade corresponding to the maximum membership degree is grade I, which means that water and sand inrush will not occur in mining. The actual mining process of this working face did not appear to have water and sand inrush, which is consistent with the prediction of the comprehensive weight-cloud model assessment method proposed in this paper. This result illustrates the feasibility of this method.

6. Conclusions

To reduce the randomness and ambiguity of the influencing factors in the prediction of the risk of water and sand inrush and to better assess the risk of water and sand inrush in short-distance mining under thick loose layers of coal mines, this paper proposes a new risk assessment method of water and sand inrush based on the comprehensive weight and cloud model. The method is applied to the risk assessment of water and sand inrush in the 6311-2 working face in the sixth mining area of Baodian Coal Mine. The following conclusions are drawn from the research: (1)Analyzing the weights determined by the analytic hierarchy process, the entropy weighting method, and the comprehensive weighting method, the weight of the bedrock thickness and the coal seam thickness are all ranked in the top three of the three weights. It can be considered that among the seven indicators that affect the risk of water inrush and sand inrush, the bedrock thickness and the coal seam thickness have a greater influence on the risk of water inrush and sand inrush and are the main influencing factors(2)The comprehensive weight-cloud model method is applied to assess the risk of water and sand inrush in the working face, and the assessment result is consistent with the actual situation. It shows that the comprehensive weight-cloud model method has good prediction performance and can provide scientific reference for safe mining under thick loose layer in deep mines in southwest Shandong(3)The comprehensive weight-cloud model method is based on the existing sample data. The number of samples, the selection of indexes, and the division of risk grade interval will have a certain impact on the assessment results of the method. In view of the complexity of water and sand inrush in close-distance mining under thick loose layer, in order to obtain more accurate prediction results, it is necessary to collect more engineering examples and sample data, so as to improve the accuracy of the method

Data Availability

All data, models, or codes generated or used during the study are available from the corresponding author by request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors of the paper would like to extend grateful thanks to the National Natural Science Foundation of China (Grant No. 51774199) for its support and Baodian Coal Mine for providing relevant geological data.