Influencing Factors’ Analysis of Disturbance Degree of Key Stratum Based on the Response Surface Method

)e key stratum controls the activities of the overlying strata or the whole strata up to the surface, which is one of the important research objects in the coal seam mining. Based on the analysis of several geological factors affecting on the key stratum, the definition of “disturbance degree of key stratum” (KSDD) was proposed. And, the KSDD is quantified by the value among 0 to 10. )rough the response surface method, experiments of three factors (mining height, buried depth, and interlayer spacing) with three different lithology types (soft, medium, and hard) between key stratum and coal seam are signed. And, the KSDD of each scheme is calculated by the developed calculation system. )e response surface regression models of KSDD with three lithology types are established. And, the single influence and interactive influences of the three factors on the KSDD with different lithology types are studied. )e results show that the following. (1) Mining height and buried depth are positively correlated with the value of KSDD, and the interlayer spacing is negatively correlated with KSDD. However, when the value of interlayer spacing exceeds 30m, the change of the KSDD tends to be gentle. (2) )e value of KSDD is not only affected by a single factor but also affected by the interaction of various factors. With the increase of burial depth, the decrease of interlayer spacing and the impact of mining height on key stratum are more severe. (3) )e influence order of each factor on KSDD is as follows: the interlayer spacing>mining height> buried depth. (4) Although the three factors interact with each other, the three factors decrease with the increase of the lithology proportional coefficient. According to the above research results, based on the calculation results of KSDD on five mines, the variation laws of KSDD with actual situation are analysed. And, the calculation results further verify the above experimental rules, which provide a certain reference and theoretical basis for the design of backfilling parameters and the management of the roof.


Introduction
Rock strata with different thickness and strength were formed above the coal seam because of different diagenetic time periods. According to the key stratum theory, the stratum which controls the whole or part of rock mass activity are called key stratum [1,2]. A large number of engineering practice show that the key stratum was hard. So, a new structure still can be formed to support the overlying strata when key stratum was broken [3][4][5]. erefore, the studies of key stratum, such as the breaking position, structural stability, and integrity of key stratum, are of great significance to the safety and production of coal seam mining [6][7][8][9][10].
A large number of scholars have enriched and expanded the theory of key stratum in order to control strata movement. Based on the medium-thick plate theory, the influences of the thickness on the fracture mode of the key layer have been studied. e criterion of fracture mode of key stratum with different thickness has been obtained [11,12]. Based on the theory of double key strata, the roof structure models of double key strata have been established, and the mining methods of graded water conservation have been obtained [13,14]. Based on the model of the key stratum structure after sliding and instability, the restability condition of the key stratum has been obtained [15]. e degree of subsidence deflection has been proposed through the establishment of the surface deflection subsidence model. Based on the key stratum theory, the influences of the thickness of loose strata and the bulking factor on the surface deflection subsidence have been analysed [16]. e stability and periodic pressure mechanism of single and double key strata in overlying strata with large mining height have been revealed. e calculation formula of support resistance has been obtained [17]. Based on numerical simulation, the influences of different distances of double seams' mining on overlying water resisting key stratum have been studied [18]. e above studies have analysed various failure laws of key stratum under different geological conditions. A series of relatively complete theoretical and surrounding rock control methods are established. However, the current research studies mainly focus on the influences of key stratum on the working face. And, there are few studies on the influences of mining geological conditions on the disturbance of key stratum. In addition, there is no corresponding evaluation method and theoretical basis. erefore, in order to explore the influences of different geological conditions on key stratum, this paper focuses on the following three research objectives. (1) Several geological factors affecting the key stratum and the proposal of the "disturbance degree of key stratum" are analysed (hereinafter referred to as KSDD). (2) e Design-Expert software will be used to design the response surface [19][20][21] experimental scheme of three factors (mining height, buried depth, and interlayer spacing) with three different lithology types (soft, medium, and hard) between key stratum and coal seam.
e response surface regression model of KSDD will be established. (3) e single influence and interactive influences of the three factors on the KSDD with different lithology types will be explored. ese studies will provide a certain reference and theoretical basis for the stability identification of key water resisting stratum, the design of backfilling parameters, and the control of surface subsidence.

e Definition of KSDD.
Underground coal mining breaks the original stress state of surrounding rock [22][23][24]. Due to the coal seam mining, the degrees of disturbance on the stratum above the coal seam are different. erefore, in order to explore the change laws of disturbance degree on key stratum with different mining geological conditions. e degree of disturbance caused by coal seam mining on key stratum is defined as the "Disturbance Degree of Key Stratum" (KSDD), which is expressed by the symbol φ.
According to the definition of the KSDD, which means the key stratum is greatly affected by the disturbance with coal seam mining with large value of KSDD, the key stratum under the same conditions is more likely to be destroyed at this time. Conversely, the key stratum is relatively affected by the disturbance with coal seam mining with a small value of KSDD. And, the key stratum under the same conditions is easier to keep intact at this time.

e Calculation Formula of KSDD.
e related factors are considered comprehensively, including quantitative discrimination method of the feasibility of upward mining in close multiple coal seams [25], mining height, buried depth [26], distance between key stratum and coal seam (hereinafter referred to as interlayer spacing), tensile strength, and rock hulking coefficient. en, according to the definition of KSDD, the disturbance intensity of the key stratum caused by coal mining is quantified by the value from 0 to 10. Formula (1) [25] is obtained: where φ is the disturbance degree of key stratum (KSDD), M is the mining height, m, H is the interlayer spacing, m, D is the buried depth of the coal seam, m, H 1 , H 2 . . . H i are the thicknesses of the ith (i � 1, . . ., n) rock above the mining coal seam, m, R t1 , R t2 . . . R ti are the tensile strength of the ith(i � 1, . . ., n) rock above the mining coal seam, MPa, K P1 , K P2 . . . K Pi are the rock bulking factors of the ith (i � 1, . . ., n) rock above the mining coal seam, n is the number of all rock strata between the Kay stratum and coal seam, λ is the influence factor of the disturbance degree of the key stratum, and C is the geological constant. e value of KSDD is closer to 10, and the disturbance effect of coal mining on the key stratum is greater. In contrast, the value of KSDD is closer to 0, and the disturbance effect of coal mining on the key stratum is smaller.

e Influencing Factors of KSDD.
e location distribution of the various factors included in formula (1) are shown in Figure 1: In order to simplify formula (1), according to the location relationship of various factors, the ratio of the interlayer spacing and the mining height (H/M) in formula (1) is defined as "disturbance influence coefficient." e ratio of the interlayer spacing and the buried depth (H/D) is defined as "interburden-to-overburden," And, the K p i · (H i /H) · R t i in formula (1) is defined as "lithology proportion coefficient," which is expressed by the symbol ξ.
In order to simplify the experimental research and analysis, based on formula (1), the lithology proportion coefficient (ξ) is expressed as follows: where ξ is the lithology proportional coefficient, α is the lithology strength coefficient, and β is the lithology ratio: (3) Figure 2 shows the relationship between the defined above and the factors in formula (1) with the KSDD.

Calculation System of KSDD.
As the geological conditions of coal seam are complex and various, the calculation process of KSDD is cumbersome. erefore, based on 2 Advances in Civil Engineering formula (1), the calculation system of KSDD is developed with C# language on ". NET" platform. e home page of the system is shown in Figure 3. According to Figure 4, the calculation of KSDD can be divided into the following four steps. (1) Click "File"-"New" to enter the calculation interface. (2) Click "Add Rock Stratum" to enter the rock information input interface. (3) Input parameters including the rock strata, the thickness of the single stratum, the cumulated thickness of the rock strata, the tensile strength, the bulking factor, the mining height, the interlayer spacing, and the buried depth. (4) Finally, click to "Calculate the KSDD" to get calculation results.
In addition, batch calculation can be performed according to Figure 5, and the calculation steps are as follows. (1) Save relevant data as ".csv" file in fixed "Excel" format. (2) Click "File"-"Import" and select the ".csv" data file to be imported. (3) After the data are imported, the

Advances in Civil Engineering
KSDD can be calculated. (4) Wait until the data calculation is completed, click "File"-"Export" to get a lot of data.

Experimental Design
In order to simplify the subsequent data analysis, the lithology between key stratum and coal seam is set as soft, medium, and hard, respectively. According to the actual geological conditions, the tensile strength of three lithology is set as 1 MPa, 3 MPa, and 5 MPa, respectively. Similarly, the bulking factor is set as 1.2, 1.3, and 1.4, respectively, as shown in Table 1.
In order to obtain the lithology proportion coefficient and show the superimposed influences between strata, three experimental schemes are set up. As shown in Table 2, the soft rock in Scheme 1 reaches 80%, which means the overall lithology between the key stratum and the coal seam is relatively soft. However, the hard rock in Scheme 3 reaches 80%, which means the overall lithology between the key stratum and the coal seam is relatively hard.
According to Table 2, the influences of various factors on KSDD are studied in three schemes. Each scheme selected mining height (M), buried depth (D), and interlayer spacing (H) as the research object. en, three levels are obtained. e final three factors and three levels' experimental scheme is shown in Table 3.

e Experimental Results.
e Box-Behnken in Design-Export [27] is used to obtain the experimental design table.
e KSDD of each scheme is calculated by the calculation system, and the results are shown in Table 4. In addition, the response surface functions of KSDD with the three kinds of lithology matching schemes are obtained as follows.

4.2.
e Analysis of Variance. e results of response surface variance are shown in Table 5.
In response surface experiment, the significance of each factor in the regression model is related to P valued. If P < 0.05, it means that the factor is significant, otherwise, it is not significant. In addition, if the F valued is larger, it means that the factor has a greater impact on the dependent variable [28,29].
As shown in Table 5, it can be seen that the three factors are significant because P valued of three factors are more than 0.05. However, when the lithology proportional coefficient increases, the significance of the model decreases. Moreover, by comparing the F valued with three schemes, the influence order of each factor on KSDD can be obtained as follows: interlayer spacing > mining height > buried depth.

4.3.
e Analysis of Single Factor. In order to more intuitively analyse the influence of each factor on KSDD, the relationships between each factor and the value of KSDD are obtained as shown in Figures 6-8. Figure 6, the values of KSDD are linearly positively correlated with the mining height. at is to say, the increase of mining height increases the disturbance effect on the key stratum. In addition, the range of KSDD between 1 m and 5 m is 0.238 in (a), the range 1 m and 5 m is 0.147 in (b), and the range between 1 m and 5 m is 0.043 in (c). It can be inferred that, as the lithology proportion coefficient increases, the impacts' range of mining height changes on the key stratum decreases.    Figure 7, the buried depth of coal seam has a positive linear correlation with KSDD. To be specific, the increase of buried depth increases the disturbance effect on the key stratum. In addition, the range of KSDD between 300 m and 700 m is 0.143 in (a), the range 300 m and 700 m is 0.088 in (b), and the range between 300 m and 700 m is 0.026 in (c). It can be inferred that, as the lithology proportion coefficient increases, the impacts of buried depth changes on the key stratum decrease. Figure 8, the values of KSDD are negatively correlated with the interlayer spacing. at is to say, the increase of interlayer spacing decreases the disturbance effect on the key stratum. en, when the value of interlayer spacing is less than 30 m, the influence is more severe. However, when the value of interlayer spacing is more than 50 m, the disturbance effects on the key stratum tend to be gentle. In addition, as the lithology proportion coefficient increases, the impacts of interlayer spacing changes on the key stratum decrease.

Influences of Disturbance Influence Coefficient.
As shown in Figure 9, with the increase of the disturbance influence coefficient, the impact on the key stratum decreases. However, when the disturbance influence coefficient increases to a certain value, the decrease of KSDD tends to be gentle. And, when the disturbance influence coefficient is less than 10, the impact on the key stratum is more serious,  but when the disturbance influence coefficient is more than 10, the impact is gentle. In addition, as the lithology proportion coefficient increases, the curve gradually becomes gentle. at is to say, the increase of lithology proportion coefficient reduces the influence range of the disturbance influence coefficient.

Influences of Interburden-to-Overburden.
Similarly, with the increasing values of interburden-to-overburden, the impacts on the key stratum decrease, as shown in Figure 10. However, when the value of interburden-to-overburden increases to a certain value, the decrease of KSDD tends to be gentle. And, when the value of interburden-to-overburden is less than 0.06, the impact on the key stratum is more serious, but when the disturbance influence coefficient is more than 0.06, the impact is gentle. In addition, as the lithology proportion coefficient increases, the curve gradually becomes gentle. at is to say, the increase of lithology proportion coefficient reduces the influence range of interburden-to-overburden. In general, the impacts of disturbance influence coefficient and interburden-to-overburden on KSDD are roughly similar. However, compared with the disturbance influence coefficient, the impact of interburden-to-overburden on the key stratum is relatively small. It further shows that the impact of buried depth on the key stratum is less than mining height.

e Analysis of Multifactors.
e KSDD is affected not only by the single factor (mining height, buried depth, interlayer spacing, and lithology proportional coefficient) but also by the interaction of factors. erefore, in order to   explore the impacts of the interaction between factors on the key stratum, based on the calculation results of KSDD, the 3D response surfaces of each factor interaction are obtained (Figures 11-13). Figure 11, with the increase of mining height and burial depth, the response surface becomes steeper, which indicates that the impacts on the key stratum increase. However, with the increase of the lithology proportion coefficient, the response surface gradually tends to be gentle. at is to say, the increase of lithology proportion coefficient decreases the range of impacts by mining height and buried depth on the key stratum. In addition, according to the distribution trend of contour, with the increase of buried depth, the disturbance effect of mining height on the key stratum is strengthened.

Interaction between Mining Height and Interlayer
Spacing on Key Stratum. As shown in Figure 12, with the increase of interlayer spacing and the decrease of mining height, the response surface gradually becomes gentle, which indicates that the impacts on the key stratum decrease. However, with the increase of the lithology proportion coefficient, the response surface gradually tends to be gentle. at is to say, the increase of the lithology proportion coefficient decreases the range of impacts by interlayer spacing and mining height on the key stratum. In addition, according to the distribution trend of contour, the impacts on the key stratum is the most severe when the interlayer spacing is within 30 m. erefore, when the key stratum is close to the coal seam, the mining height can be appropriately reduced to slow down the disturbance effect of coal mining on the key stratum.

Interaction between Buried Depth and Interlayer
Spacing on Key Stratum. As shown in Figure 13, with the increase of interlayer spacing and the decrease of coal seam buried depth, the response surface gradually becomes gentle, which indicates that the impacts on the key stratum decrease. However, with the increase of lithology proportion coefficient, the response surface gradually tends to become gentle. at is to say, the increase of lithology proportion coefficient decreases the range of impacts by interlayer spacing and buried depth on the key stratum. In addition, according to the distribution trend of contour, it can be seen that, with the increase of buried depth, the range of impacts by interlayer spacing on the key stratum increases. erefore, when coal seam is in the deep area or the coal seam is close to the key stratum, the key stratum is more seriously affected by coal seam mining. And the method of backfilling can be used to reduce the disturbance effect of coal seam mining on the key stratum.

Engineering Verification
Based on the geological data of Longde Coal Mine [30], Yuhua Coal Mine [31], Da'anshan Coal Mine [32], Qianjiaying Coal Mine [25], and Fangezhuang Coal Mine [33], the values of KSDD in each mining area is calculated by the calculation system of KSDD. e results are shown in Table 6.
As shown in Table 6, when the buried depth and interlayer spacing are basically the same, the values of KSDD are affected by mining height and lithology proportion coefficient. erefore, roof management should be strengthened in areas with large mining height to prevent roof fall and other disasters caused by mining height increase. When the mining height and interlayer spacing are basically the same, the values of KSDD are mainly affected by the buried depth and the lithology proportion coefficient. erefore, when mining in Yuhua Coal Mine and other areas with large buried depth, the key stratum is more likely to be destroyed, which makes the overlying strata more difficult to control and easily bring disasters such as rock burst. When the lithology proportion coefficient decreases, the values of KSDD are greatly affected by the changes of mining height, buried depth, and interlayer spacing. at is to say, when there is  Advances in Civil Engineering less intermediate hard rock, the changes of mining height, buried depth, and interlayer spacing need to be paid attention to, which will lead to more damage to the key stratum. In general, the above rules are basically consistent with the actual situation, which further verifies the above experimental results.

Conclusion
Synthesizing the experiment designed of response surface, calculation results of KSDD, and analyses of single factor and multifactors, the following main conclusions are drawn: (1) ough the analysis of the key stratum affected by several geological factors, the definition of "disturbance degree of key stratum" (KSDD) was proposed.
To be specific, KSDD is the degree of disturbance caused by coal seam mining on the key stratum. Based on the calculation formula of KSDD, the calculation system is developed. And, the value of KSDD is quantified by the value among 0 to 10. In addition, the response plane method is used to design the experiment of three kinds of lithology proportional coefficient with a different mining height, buried depth and interlayer spacing. e response surface regression models of three lithology proportional coefficients are established.
(2) By analysing the disturbance effect of a single factor on the key stratum and its significance, the following are found. ① e order of influence intensity of each factor on the KSDD is as follows: interlayer spacing > mining height > buried depth. ② Mining height and buried depth are positively correlated with the value of KSDD. ③ e interlayer spacing is negatively correlated with KSDD, but when the interlayer spacing increases to a certain value, the influence range of KSDD tends to be small. (3) e KSDD is affected not only by a single factor but also by the interaction of various factors. at is to say, with the increase of buried depth, the impacts of mining height on the key stratum increase. And, with the decrease of interlayer spacing, the disturbance effects of mining height on the key stratum are more severe. In addition, with the increase of buried depth, the change of interlayer spacing increases the influence range of KSDD. However, with the increase of lithology proportional coefficient, the influence range of KSDD decreases with the three factors.
(4) e actual geological data of five mines are selected to calculate the value of KSDD. And, the variation law of KSDD in each mine is analysed. e calculation rules of KSDD of five mines further verified the above experimental results, which provide some reference and theoretical basis for the management of the roof and the design of backfilling parameters.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper. Advances in Civil Engineering 9