Abstract

The composite roof structure of coal roadway is complex, and its stability is related to lithology of strata, strata thickness, number of strata, strata location, and interlayer cohesive force. Based on a simplified model of the composite roof structure, simulation test, analytic hierarchy process (AHP), and Matlab, programs are employed to comprehensively analyze the structural stability of roofs. The structural stability of the composite roof is demonstrated to decrease with the lithology of strata, strata thickness, and interlayer cohesive force and increase in number of strata and distance of hard-and-thick strata from the roadway. The decreasing order of influence of these factors on the composite roof structure is as follows: lithology of strata, strata thickness, number of strata, strata location, and interlayer cohesive force. Lithology of strata is the main factor affecting stability. The bending strength decreases with the increase in number of strata, and the influence of strata position on stability decreases from the first layer to the fourth layer. According to the AHP, an expression for the comprehensive influence coefficient (k) for the stability of the composite roof structure is proposed, and the surrounding rocks are divided into three levels using this equation: grade I with 0.7 < k < 1 (stable), grade II with 0.4 < k < 0.7 (slightly stable), and grade III with 0 < k < 0.4 (unstable). A scientific basis to evaluate the stability and control of the composite roof of a complex structure is thus provided.

1. Introduction

The composite roof structure of a coal seam is complex and changeable, and its stability seriously affects the stability of the surrounding rock of a roadway; this frequently affects the surrounding rock support of the roadway and leads to roof caving accidents [14]. Therefore, research on the structural stability of the composite roof is critical to the stability analysis of the surrounding rock of the roadway.

The composite roof generally comprises rock that is layered by soft or hard rock interbred or containing weak interlayers and has weak plane development characteristics for the stratum structure, many weak interlayers, small strata thickness, and low interlayer cohesion force [57]. Many scholars have studied and discussed the stability of composite roof structures of coal seams [812]. The main factors affecting the stability of the composite roof structure are considered to be the lithology of strata, strata thickness, number of strata, strata location, and interlayer cohesive force. The comprehensive strength index of the composite roof is introduced to analyze stability; however, the equation does not consider the influence of strata location, number of strata, and interlayer cohesion force [1315]. The application of this equation has certain limitations. In this paper, the theoretical analysis, laboratory simulation test, and analytic hierarchy process (AHP) [1619] are applied to propose an expression for the comprehensive influence coefficient (k) of the composite roof according to the influence of each main factor on the stability. This work provides support for more comprehensive and systematic analyses and evaluation of the structural stability of the composite roof.

2. Theoretical Analysis of the Stability of Composite Roof

For composite roof structures, the ratio of the thickness of each stratum to the span of coal roadway must satisfy the thin plate hypothesis [20, 21] (i.e., the thickness of the rock stratum is less than 1/5 of the span of coal roadway). In theoretical analysis, the composite roof is generally simplified for simply supported beam mechanical models [22]. As shown in Figure 1, the flexural strengths of the composite roof under different strata lithology, strata thickness, and number of strata are calculated to measure the stability of the composite roof structure. The maximum stress in the span of rock beam in each stratum is as follows:where σi represents the maximum stress in the rock beam span (MPa); i = 1, 2, 3; q stands for the vertical uniform load (kN); N displays the horizontal axial force (kN); l is the rock beam span (m); and b is the width of the rock beam. Taking b as 1 m, E1, E2, and E3 are the elastic moduli of the rock; I1, I2, and I3 are moments of inertia; and h1, h2, and h3 are the thickness of the rock formation (m).

The above equation indicates that the tensile strength of the composite roof structure decreases with the elastic modulus and increase in the number of strata when the total thickness of the composite roof structure remains constant under certain conditions of q, N, and l, and the weak interlayer with lesser thickness easily forms a stress concentration area compared with other strata. Although it is concluded that the stability of the composite roof is related to the lithology, thickness, and number of strata, the mechanical theory does not consider the influence of strata location and interlayer cohesive force and cannot fully reflect the influence of the structural characteristics on stability. Therefore, based on the mechanics theory, the influence of structural characteristics of the composite roof on stability is studied by a simulation test.

3. Experimental Study on Stability of Composite Roof for Various Factors

3.1. Test Scheme

Under the action of vertical ground and horizontal tectonic stress, the latter will weaken the tensile stress of the lower edge of each rock layer of the composite roof structure of the coal roadway when the bending and sinking of the roof structure is low. Furthermore, the latter will produce a second-order effect when the bending and sinking of the roof structure is high, aggravating the flexural deformation of the rock stratum and easily leading to the instability and damage of the composite roof coal roadway. Therefore, the stability of the composite roof structure can be analysed by measuring the bending strength of the composite roof structure.

Based on the simplified composite beam mechanics theory of the composite roof structure [13], two different lithology test materials were selected to establish the experimental model ((a) in Figure 2): yellow and black sandstone. The physical and mechanical properties of these materials are listed in Table 1. In the course of the test, by varying the number of strata, position of the strata, and interlayer cohesive force, concentrated loads were applied to different composite beam structures (Figure 2), and the bending strength was tested (Table 2) to analyze the stability of the composite roof structure. Because the interlayer cohesion force of the composite roof is weak in practical engineering, low strength cement mortar is used to simulate the interlayer cohesion force of rock strata ((b) in Figure 2).

3.2. Analysis of Test Results
3.2.1. Influence of Lithology and Interlayer Cohesion on the Stability of Composite Roof

To study the influence of the composite roof structure on the stability of the roadway surrounding rock, a comprehensive strength index [13] was applied, which is defined as follows:where qk displays the comprehensive strength index of the composite beam, qi is the strength of each rock layer of the composite beam, hi represents the thickness of each rock layer of the composite beam, and h is the total thickness of the composite beam.

The following observations can be drawn from the comparison of the bending strengths of the combined models of different lithologies (Figure 3):(1)For the two-to-four-layer composite models, the flexural strengths decrease with the comprehensive strength index, ranging from black to yellow sandstones, and the average reduction rates of the flexural strengths of the combined models with interlayer cohesion force, but no interlayer cohesion force is 0.38. Thus, the lithology is indicated to have a great influence on the stability of the coal seam composite roof under vertical ground stress.(2)The flexural strengths of the composite models without interlayer cohesive forces are lower than those of the composite models with interlayer cohesive forces, with an average reduction rate of 0.26, which indicates that interlayer cohesive force is one of the prominent factors affecting the stability of the composite roof of coal roadway.(3)To consider the influence of interlayer cohesive force on the stability of composite roof, according to the linear relationships between the compressive strengths of the composite models without and with interlayer cohesive forces, as shown in Figure 4, the coefficient A is introduced to improve the expression of the comprehensive strength index; the coefficient A is 1 when there is no interlayer cohesive force and is 1.3333 when there is interlayer cohesive force. The comprehensive strength index qk1 is given by the following equation:

3.2.2. Influence of the Number of Rock Layers on the Stability of Composite Roof

By comparing the bending strengths of different combinations of layers (Figure 5), the following observations are made:(1)Under the vertical load, the bending strengths of the combined models with the same thickness decrease with increase in the number of layers. The average bending strength ratio of the one-layer, two-layer, three-layer, and four-layer interlayer cohesive force models is 1.00 : 0.67 : 0.53 : 0.47, and the average reduction rates are 0.00, 0.34, 0.47, and 0.54, respectively. When the number of layers in the combined model are more, the rate of reduction is lower, and with the increase of the number of layers, the stability of the roof is reduced by the influence of the variation of number of layers; furthermore, the stability of the roof stratum for more than four layers is stabilized by the influence of the number of layers. Table 2 indicates that the combined strength indexes of the one-to-four-layer combined models are the same for identical thicknesses, indicating that the comprehensive strength index does not consider the influence of the number of layers on the stability of the composite roof.(2)To analyze the influence of strata number, the coefficient Bi is introduced to improve the equation of the comprehensive strength index, whose value is shown in Table 3. The value of Bi is determined using the ratio of the average bending strength of the combined model of different layers of the same thickness of the composite beam (one layer :  two layers :  three layers :  four layers = 1 : 0.67 : 0.53 : 0.47) are considered. An increase in the number of layers in the rock layer reduces the stability of the roof. When the rock layer reaches more than four layers, the composite roof structure tends to stabilize, and the value of Bi above the four layers remains constant. The comprehensive strength index qk2 is given by the following equation:

3.2.3. Influence of Position on the Stability of Composite Roof

Figure 6 is a comparison of the flexural strengths of the combined models at different locations. Furthermore, combining with the data of Figure 6 and Table 2, the following are noted:(1)By changing the position of each rock stratum, the farther the black sandstone with better lithology and the thicker the rock stratum from the roadway, the higher the strength of the composite roof structure. However, with the increase in the number of strata, the influence of strata location on the strength of the composite roof is smaller. Because the horizontal load is not applied to the composite beam during the test, the influence of rock strata location on the composite beams is contrary to the deformation and failure law of the roof of the actual coal roadway [23]. Thus, in conclusion, the closer the hard thick rock layer is to the roadway, the greater the bending strength of the composite roof, and the influence of each rock layer on the overall strength decreased from the first layer to the fourth layer.(2)From the two-layer to four-layer combination models, the comprehensive strength indexes are equal; however, the flexural strengths are different. This indicates that the comprehensive strength index does not consider the influence of rock strata location on the stability of the composite roof structure. The comprehensive strength index is further improved as given below.

The analysis demonstrates that, under the action of a vertical load, when the thickness of the composite roof structure of the coal roadway is constant, the stability decreases with strata lithology, strata thickness, and interlayer cohesive force and the increase in strata number and distance of the hard-and-thick-strata from the roadway. The factors affecting the composite roof structure can thus be arranged in the decreasing order as follows: lithology of strata, strata thickness, number of strata, strata location, and interlayer cohesive force.

4. Comprehensive Analysis of Various Factors Based on AHP

4.1. Principles of AHP

AHP is a method of quantitative analysis based on qualitative analysis. The purpose of this method is to decompose complex problems into several levels and several factors and to compare the importance of these two categories to determine their weights [24, 25]. The steps are illustrated as follows:(1)Determine the factors for analysis.(2)Construct a judgment matrix according to the principles of AHP. The relevant indicators and factors of each layer are compared using the 9-level scale method (summarized in Table 4) to establish the judgment matrix F. Here, bij is the element in the judgment matrix, F =(bij) n × n, where bij > 0, bji = 1/bij, and bii = 1.(3)The maximum eigenvalues of the judgment matrix and the corresponding eigenvectors were calculated using Matlab software, and the relative weights of each factor are obtained.(4)Judgment matrix consistency test: owing to the subjective judgment of the person involved in the establishment of the judgment matrix, the judgment may be inconsistent. To ensure the validity of the evaluation, the consistency ratio CR is required to check the consistency of the judgment matrix:where CR represents the consistency ratio. When CR < 0.1, the degree of inconsistency of F is within the allowable range, and the judgment matrix is feasible. Here, CI displays a consistency indicator, CI = (λ − n)/(n − 1), where λ is the maximum eigenvalue of the judgment matrix F. n is the sum of the eigen roots of F. RI is the average random consistency test index [26]. The value of RI for n = 1, 2, 3, ..., 10, is shown in Table 5.

4.2. Calculation of Rock Layer Position Weight
4.2.1. Constructing a Judgment Matrix

According to a qualitative analysis of the factors influencing the stability of the composite roof structure by simulation tests combined with the principle of AHP, the relative importance degree between the two rock layers is analysed. According to the 9-scale method, the judgment matrices F1, F2, and F3 of the 2, 3, and 4 layers can be constructed as follows:

4.2.2. Calculation Results and Analysis

The contribution rate of the composite roof rock stratum (the weight of each rock layer) is summarized in Table 6. The maximum eigenvalues and eigenvectors of the matrices F1, F2, and F3 are calculated using Matlab software. The maximum eigenvalues are λ1 = 2.0000, λ2 = 3.0092, and λ3 = 4.0310, and the eigenvectors are normalized to calculate the weight ratios of each rock layer as , , and .

By substituting each of the rock stratum eigenvalues into equation (5), the consistency ratios CR1 = 0.0000, CR2 = 0.0089, and CR3 = 0.0116 are calculated, all of which are less than 0.1000, indicating that F1, F2, and F3 are within the allowable range of inconsistency and meet the consistency test.

According to the AHP, the contribution rates () of each rock layer of the two-layer, three-layer, and four-layer composite roof are obtained, and the comprehensive strength index expression qk3 of the composite roof is improved as follows:

4.3. Calculation of Weights of Influencing Factors of Composite Roof Structure
4.3.1. Constructing a Judgment Matrix

According to the research of the experiment and inquiry with experts, the degree of influence of the formation lithology, thickness of the stratum, number of strata, position of the stratum, and cohesion of the interlayer on the stability of the composite roof is qualitatively given. The judgment matrix F4 can be constructed according to the 9-scale method as follows:

4.3.2. Calculation Results and Analysis

Matlab software was used to calculate the maximum eigenvalue and eigenvector of matrix F4, in which the largest eigenvalue is λ4 = 5.0681. The eigenvectors were normalized to calculate the weight ratio of each rock layer as . By substituting the eigenvalues λ4 of each rock layer into equation (5), the consistency ratio CR = 0.0150, which is less than 0.1000, indicating that F4 is within the allowable range of inconsistency and satisfies the consistency test. Therefore, the influence of stratum lithology, rock thickness, stratum layer, rock stratum location, and interlayer cohesive force on the stability of the composite roof structure is 0.4185, 0.2625, 0.1600, 0.0972, and 0.0618, respectively.

5. Classification of Composite Roof Stability Based on the Comprehensive Influence Coefficient

5.1. Composite Roof Comprehensive Influence Coefficient

According to the analysis of the factors affecting the stability of the composite roof structure, the influence of strata lithology, rock thickness, strata layer, interlayer cohesive force, and the comprehensive influence coefficient of the composite roof structure is determined. Because equations (3), (4), and (6) are derived based on equation (2), the influence of lithology and thickness of strata has been considered in equation (2); the degrees of influence of strata lithology, strata thickness, number of strata, location, and interlayer cohesive force on the stability of composite roof structure are 0.4185, 0.2625, 0.1600, 0.0972, and 0.0618, respectively. When the influence of number of strata, strata location, and interlayer cohesive force is further considered, the degrees of influence of number of strata, strata location, and interlayer cohesive force based on the strata lithology and strata thickness are then 0.23, 0.14, and 0.09, respectively. Considering the comprehensive equations (2)–(4), and (6) and the ratio of the influence of the interlayer cohesive force, strata location, and number of strata location on the stability of the composite roof structure (0.23 : 0.14 : 0.09) along with the lithology of the rock stratum, the expression k of the comprehensive influence coefficient of the composite roof structure is derived:where k is the comprehensive influence coefficient of the composite roof structure, is the strength contribution rate of each stratum of composite roof (), qi is the strength of each weak layer of composite roof, hi is the thickness of each weak layer of composite roof, h is the total thickness of composite roof strata, and A is the coefficient of adhesion between layers. Taking A as 1 when there is no interlayer cohesive force and A as 1.3333 when there is interlayer cohesive force, qimax is the best rock layer strength for composite roof stability, Bi is the layer coefficient of each rock layer of the composite roof, and these values are summarized in listed Table 3.

5.2. Composite Roof Stability Classification

Because the coal seam composite roof has many weak interlayers, different rock layers with different strengths, and large thickness differences, the stability is the result of the combination of lithology of strata, strata thickness, strata location, number of strata, and interlayer cohesive force. The stability of the composite roof structure is further classified according to the results of the comprehensive influence coefficient calculation, as summarized in Table 7.

6. Conclusion

In this study, based on the simplified theoretical theory of composite beam mechanics, combined with the simulation test, AHP, and Matlab software, the influence of the rock layer strength, thickness of the rock stratum, number of strata, position of the stratum, and interlaminar adhesion force is considered comprehensively; the comprehensive influence coefficient of the composite roof structure is proposed, and its stability is further analysed. The following detailed conclusions are drawn:(1)Under the vertical load, the stability of the composite roof structure of the coal seam decreases with the lithology of the stratum, interlayer cohesive force, and thickness of the stratum as well as the increase of the stratum and distance between the hard thick stratum and the roadway; the influence degree of each factor on the composite roof structure from large to small is as follows: lithology, thickness of strata, number of strata, location of strata, and interlayer cohesion force. After changing the position of the weak strata, the closer the hard and thick strata are to the roadway, the greater is the flexural strength of the composite roof, and the influence of the number of strata on the overall strength decreases from the first layer to the fourth layer. After introducing the comprehensive strength coefficient, it is found that the index does not consider the influence of the number of layers, location of the strata, and interlayer cohesive force on the stability of the composite roof structure, which requires further improvement based on future investigations.(2)Considering the influence of the location of the rock stratum, the number of strata, and the interlayer cohesive force on the composite roof structure, the comprehensive strength index equation is improved. Furthermore, the weight between the influencing factors of the composite roof structure stability and the weight of the position of each stratum on the stability of composite roof are obtained quantitatively using the AHP. Based on this, the equation of the comprehensive influence coefficient is proposed, which provides a basis for the stability analysis of the composite roof structure.(3)The stability of the composite roof structure can be divided into three levels using the equation of the comprehensive influence coefficient: stable (0.7 < k < 1), slightly stable (0.4 < k < 0.7), and unstable (0 < k < 0.4). The smaller the comprehensive influence coefficient is, the worse the stability of the composite roof structure is in the coal seam.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (51174159 and 59974019), the China Postdoctoral Science Foundation (2015M572580), and the Scientific Research Plan Projects of Education Department of Shanxi Province of China (15JK1471).