Abstract
As mining height increases, the influence of coal wall caving on safety production becomes stronger. There is no systematic and effective method to analyse the risk of coal wall caving and its slip caving depth. First, this paper established the Bishop mechanical model of sliding instability of coal wall, and then it deduced the general equation of a safety factor for every slip surface, which can be used to judge the stability of the coal body on the slip surface. Moreover, taking the 8102 workface in the Wulonghu Mine, China, as an example, this paper evaluated the calculation method of slip surface safety factor in detail and obtained the critical slip surface position and the maximum slip depth of a coal wall. Overall, the results showed that the maximum slip depth based on the Bishop strip method is more consistent with the measured data compared with other methods and thus has strong significance and practical engineering value for selecting the most suitable method and its parameters of regulating coal wall caving.
1. Introduction
Thick coal seams account for a large proportion of China’s coal resources, and fully mechanized mining with large mining height is the main mining method for 3.5∼8.2 m thick coal seams [1–3]. However, as mining height increases, the pressure on the workface is more and more intense, such that coal wall caving becomes one of the main factors restricting safe and efficient production [4–6]. Severe coal wall caving not only affects the normal production of the workface but also threatens the life and safety of miners [7]. Therefore, controlling coal wall caving has become a key technical issue for fully mechanized mining with large mining heights that has attracted the attention of scholars and engineers. Fu et al. [8] established the coal wall excavation unloading effect model by using the unloading rock mass mechanics theory and the fracture mechanics theory to analyse the coal wall instability mechanism. The research results show that cracks in the coal wall would change into mixed mode cracks under the effects of the unloading stress field. Pang and Wang [9] divided the process of hard coal seam rib spalling into tensile cracking damage and sliding instability based on analysing the effect of the stress path in the process of rib spalling and established the tensile cracking-sliding mechanical model of rib spalling for hard and thick coal seams. Hua and Xie [10] found that the stress state of the coal wall after excavation changes from a three-direction stress state to a two-direction stress or one-direction stress state. Yan et al. [11, 12] used the theory of the “pressure bar” to analyse the deflection characteristics of a coal wall with good integrity. As a result, it is concluded that the part that is most prone to cave is in the middle and upper parts of the coal wall. Yuan et al. [13] built a mechanical model of wedge-shaped sliding coal body of a coal wall workface and analysed the key factors of the sliding instability of a wedge-shaped sliding coal body. Wang et al. [14, 15] simplified the trace of coal wall caving from the arc type to the oblique straight line type and considered the caving coal body as a triangular block. Lastly, combined with the safety margin criterion, an equation for judging whether coal wall caving is generated is obtained. Yao et al. [5, 16–19] analysed the factors influencing coal wall caving, including coalface depth, mining height, panel advance velocity, coal strength, and coal seam dip angle. According to these factors, a series of methods for regulating coal wall caving are proposed, including optimizing the production process, optimizing the hydraulic support system, water injection, and wood-bolt reinforcement.
In summary, the existing research results are mostly focused on analysing the mechanisms of the coal wall caving, the forms of coal wall caving, the factors affecting coal wall caving, and the technical measures to prevent and control coal wall caving. However, an effective method for analysing the coal wall slip risk and the coal wall caving depth has not been found. This status directly limits the accurate determination of the method and its parameters for controlling the coal wall. In view of the above engineering problems, this paper draws on Bishop’s strip method, which is commonly used in slope stability analysis, and applies the Bishop method to the analysis of the coal wall stability of a workface. Considering the roof pressure, the safety factor of any slip surface of a coal wall of the workface can be calculated by the Bishop method, and the exact value of the maximum caving depth can be obtained. This research can provide a theoretical basis for the determination of the methods and their parameters for controlling coal wall caving in a large mining height workface.
2. Method Introduction
The Bishop strip method [20, 21] is a slope stability analysis method proposed by Professor A. W. Bishop of Imperial College London in the 1950s. Based on the Mohr–Coulomb criterion, the Bishop method assumes that the sliding surface of the slope is a circular arc surface and the sliding body is vertically divided into several vertical soil strips. At the same time, the method also assumes that the antisliding safety factors on the sliding surface of each soil strip are the same and are equal to the average safety factor of entire sliding surface. The Bishop strip method is widely used in solving the safety factor of slope stability engineering because the calculation results are more realistic [22–25].
The coal body near the coal wall of the workface is ruptured due to the abutment pressure, and the shear strength is greatly reduced. Coupled with the cutting disturbance of the shearer drum, the coal body easily slips along a “sliding surface” that is inclined or approximates a circular arc, as shown in Figure 1. In addition, field test results show that this form of coal wall caving accounted for more than 75% of the total number of working face cavings [13, 26]. Therefore, coal wall caving can be regarded as the mechanical process of the sliding coal body sliding along an arc-shaped sliding surface through the foot of the coal wall. Considering the roof pressure, the Bishop method is used to establish an improved mechanical model of the coal wall caving, and the slip surface safety factor is calculated. According to the calculation result of the safety factor, the stability of the coal body on the slip surface can be judged.

3. Calculation of the Safety Factor of a Coal Wall Slip Surface
3.1. Calculation of Safety Factor of Coal Wall Slip Surface during Horizontal Mining
When the length of the workface is much larger than the caving depth, the workface per unit length can be treated according to the plane strain problem. Using the Bishop strip method principle, a simplified mechanical model of the coal wall slip caving is established, as shown in Figure 2(a). Suppose that the slip surface AB is a circular arc surface with O as the centre, the radius of the arc is R, the distance between the centre, O, the roof is y at the top of the roof, point A is at the foot of the coal wall, point B is on the contact surface between the coal seam and the roof, AC is the height of the workface, and the distance d of BC is defined as the slip depth of the slip surface, AB. The slip caving body, ABC, is divided into k vertical coal bars.

(a)

(b)
According to the force characteristics of the coal wall, the m-th coal strip is taken for force analysis, as shown in Figure 2(b). Fm is the normal force on the contact surface between the roof and the coal strip; xm is the horizontal distance from the centre, O, to Fm, that is, the force arm of Fm to O; sm is the horizontal distance from the centre, O, to the centre line of the coal strip; fm is the frictional resistance of the roof to the coal strip; y is the vertical distance from the centre, O, to fm, that is, the force arm of fm to O; Ym, and Xm, are the tangential force and normal force between adjacent coal strips, respectively; Wm is the self-weight of the coal strip; Nm is the normal force on the slip surface of the bottom of the coal strip; Tm is the tangential force on the slip surface; is the angle between the horizontal direction and the tangent line of the slip surface at the bottom of the coal strip; lm is the length of the bottom of the coal strip along the tangential direction; and hm is the height of the coal strip.
Use K to indicate the safety factor on the slip surface AB, and thenwhere is the average shear strength along the entire slip surface and is the shear stress along the entire slip surface.
When K > 1, the coal body is stable; when K < 1, the coal body is unstable; and when K = 1, the coal body is in a critical instability state.
According to the Mohr–Coulomb criterion, the average shear strength along the entire slip surface AB iswhere c is the coal body cohesion, is the normal stress on the shear plane, is the pore water pressure, and is the internal friction angle of the coal body.
The pore stress, , is assumed to be zero because the coal body is on the workface; thus, the above equation is simplified to
The Bishop strip method assumes that is balanced with the tangential force Tm on the slip surface, i.e,
According to the vertical balance condition of the m-th coal strip, equation (5) can be obtained:
Substituting equation (4) into equation (5),where .
When the sliding body reaches the limit equilibrium state, the sum of the torque of the coal strip to the centre, O, of the slip surface should be zero. At this point, the force between the adjacent coal strip should be regarded as the internal force and will cancel each other out:
Equations (4) and (6) can be substituted into equation (7):
Equation (8) is a general expression of the safety factor of the slip surface calculated by Bishop’s strip method. In this equation, Ym and are unknowns. To find K, the value must be estimated, but Bishop proved that if , the error caused is only 1%; therefore, it is assumed that the vertical shear forces on both sides of the coal strip are equal, i.e., . Substituting into equation (8), the simplified Bishop safety factor K expression is obtained:where , in which is the vertical stress acting on the coal strip by the roof; , in which is the bulk density of the coal strip; and , in which is the frictional resistance coefficient.
3.2. Calculation of the Safety Factor of a Coal Wall Slip Surface during Up-Dip/Down-Dip Mining
The Bishop mechanical model of up-dip and down-dip mining, when the coal seam dip angle is , is shown in Figures 3(a) and 3(b).

(a)

(b)
At this time, the safety factor, K, expression of the slip surface AB iswhere .
When the workface uses up-dip mining, takes a positive value; when the workface uses down-dip mining, takes a negative value; and when , equation (10) is reduced to equation (9).
Equation (10) is a general equation for calculating the safety factor, K, of a coal wall slip surface by the simplified Bishop strip method. Since the right-hand side of the equation also contains K, equation (10) needs to be solved by an iterative method. First, assume a K value to find , and then substitute into equation (10) to find K. If this K value does not match the assumed value, re-substitute the K value into and re-substitute into equation (10) to obtain a new K value. Repeat these steps until the assumed K value is close to the obtained K value. Generally, K = 1 is assumed first. For the convenience of the calculations, the value can be directly found by the values of and , as shown in Figure 4 [20, 21].

4. Field Application
The 8102 workface of the Wulonghu Coal Mine is a monoclinic structure with N-S-trending and E-W-tendency. The dip angle is 5°∼25°, with an average value of 18°, which is a moderately inclined coal seam. The dip length of the working face is 156–233 m, the strike length is 623–802 m, the mining height is 4.1 m, and the burial depth is approximately 600 m. Fully mechanized mining with a large mining height is used for full thickness mining, and the method of all roof caving is used to deal with the roof. The bulk density, , of the coal body is 14.0 kN/m3, the cohesion, c, is 1.0 MPa, and the internal friction angle, , is 25°.
4.1. Determination of the Slip Surface Location
To find the position of the critical slip surface with the safety factor K = 1, several slip surfaces of several different positions should be assumed first. Calculating every safety factor for each slip surface by equation (9), the critical slip surface is the slip surface whose K value is equal to or closest to 1. It can be seen in Figure 2 and equation (9) that the arc surface of the crossing point, A, is determined by the two parameters of the sliding arc radius, R, and the centre height, y, so that different values of R and y are combined to obtain a plurality of slip surfaces.
Through a field measurement, the maximum caving depth is d ≤ 1.9 m, and the minimum caving depth in the workface that affects the normal production is d ≥ 0.7 m. The radius of the slip surface is estimated to be 12∼20 m from the caving depth and the height of the coal seam (4.1 m); therefore, the radius Rj values are 12 m, 14 m, 16 m, 18 m, and 20 m, respectively. When the radius R of the slip surface is determined, the height of the centre of the circle is larger, and the depth of the determined slip surface is larger. When the height of the centre of the circle is 6.5 m, the minimum depth of the determined slip surface is 1.954 m, which is larger than the maximum caving depth of the field measured value. Therefore, the height yi values of the centre of the circle are 0.5 m, 2 m, 3.5 m, 5 m, and 6.5 m, respectively. Therefore, 25 slip surfaces of ABij are obtained from the combination of Rj and yi, and the slip depths of each slip surface are shown in Table 1.
4.2. Calculation of Safety Factor Kij
According to the mining height and the laws of strata behavior, it is estimated that the maximum slip depth of the coal wall is 1.4∼1.9 m. According to the engineering characteristics of the coal wall caving and the reinforcement methods of the coal wall, the prediction accuracy of the maximum slip depth is 0.2 m. The Bishop strip method is used to divide the coal body of each slip surface into vertical coal strips with widths of bm = 0.2 m. Figure 5 shows the coal strip division corresponding to the slip surface determined by R3 = 18 and y3 = 3.5 m. A total of 7 coal strips are divided, and their geometric parameters are shown in Table 2.

Coal wall caving generally occurs in the plastic zone, and the values of the normal stress, qm, on the contact surface between the roof and the coal strip can be directly calculated by equation (12) [27], and the results are shown in Table 2.where c is the coal body cohesion, MPa; is the friction coefficient of the contact surface of the coal seam and the roof; p is the horizontal thrust of the sidewall guard plate, which is not considered here; M is the thickness of the coal seam, m; is the distance of the coal strip from the edge of the coal wall, m; and and are generalized von Mises criterion coefficient.
Substituting the geometric parameters and mechanical parameters of the coal strip into equation (9), the K33 value is solved iteratively, and the safety factor on the slip surface is 0.896. Similarly, the safety factors of other slip surfaces can be calculated. Figure 6 shows the iterative curve of the safety factor Ki3 values of the five slip surfaces determined by radius R3 = 18 m and yi = 0.5 m, 2 m, 3.5 m, 5 m, and 6.5 m. It can be seen from the figure that when the centre height, y, is 5 m, or 6.5 m, the safety factors of the determined slip surfaces AB43 and AB53 are all greater than 1, which indicates that the possibility of the coal wall sliding down along these slip surfaces is small; when y is 0.5 m, 2 m, or 3.5 m, the safety factors of the determined slip surfaces AB13, AB23, and AB33 are less than 1, which indicates that the possibility of the coal wall sliding down along the slip surface is great.

The safety factors of the 25 slip surfaces of 12 m, 14 m, 16 m, 18 m, and 20 m and yi = 0.5 m, 2 m, 3.5 m, 5 m, and 6.5 m combinations are shown in Table 3.
The change contours of the safety factor, K, with the slip arc radius, R, and the centre height, y, are shown in Figure 7. It can be seen that the contour line of K = 1 is a continuous and smooth curve, which indicates that the safety factors of a series of slip surfaces are equal to 1. The K value is generally small in the lower part and large in the upper part, and the slip surface safety factor determined by R and y corresponding to the coordinates of any point at the bottom of the curve K = 1 is less than 1. According to the positional characteristics of the slip surfaces determined by R and y shown in Figure 7, it is shown that the safety factor increases with the increase of the slip surface slip depth in an up-dip mining workface.

4.3. Determination of the Maximum Caving Depth
It is known from Figure 7 that there are many slip surfaces with a safety factor of K = 1, so it is necessary to determine the location of the critical slip surface with the maximum slip depths, and then the maximum caving depth can be obtained. The change rules of the safety factor, K, with slip arc radius, R, and centre height, y, are fitted to the K surface in Figure 8 by using MATLAB software. The change rules of the slip depth with slip arc radius, R, and centre height, y, are fitted to the d surface. A continuous and smooth Kcr curve is formed by intersecting the plane of safety factor K = 1 with the K curved surface, projecting the Kcr curve onto the d surface to obtain a spatial dcr curve, finding the maximum F point on the dcr curve and obtaining F coordinates (17.96, 4.20). When the radius of the coal wall slip surface is 17.96 m and the centre height is 4.20 m, the slip surface is the farthest critical slip surface from the coal wall.

(a)

(b)
From Figures 7 and 8, the overall trend is that the smaller the coal wall caving depth is, the smaller the coal wall safety factor is and the more unstable the coal wall is. So, the farthest critical slip surfaces divide the coal body in front of the coal wall into two areas, as shown in Figure 9. The safety factors of the slip surfaces in area ① are less than or equal to 1, which indicates that the coal body may easily slip along the slip surface of a certain K < 1 in this area. Additionally, the coal body in this area may easily be disturbed by a coal cutter, so it is prone to local caving. The safety factors of the slip surfaces in area ② are greater than 1, which indicates that the coal body will not slide along a slip surface in this area. Therefore, the slip surface AB is the critical slip surface of the coal wall, and the maximum coal wall slip depth is 1.63 m. Based on this, the control method and parameters of the coal wall can be designed accurately.

4.4. Comparative Analysis of Prediction Accuracy
In references [10, 28], based on the Mohr–Coulomb strength theory, the equations for calculating the maximum coal wall caving depths for compression-shear-type failure and sliding-type failure are (13) and (14), respectively:where M is the thickness of the coal seam, m; is the coal seam dip angle, and is the internal friction angle of the coal body.
Based on the aforementioned 8102 workface of the Wulonghu Coal Mine as background, the maximum depth is calculated by using equations (13) and (14), respectively, and compared with the actual field measurement, and the accuracy of each method is analysed. The results are shown in Table 4.
Table 4 shows that the maximum coal wall caving depth obtained according to formula (13) is 4.96 m, which is 3.14 m from the measured value. The maximum coal wall caving depth obtained according to formula (14) is 33.37 m, which is 31.55 m different from the measured value, indicating that the formula is not suitable for predicting the coal wall caving depth of the large mining height; the maximum coal wall caving depth obtained by the Bishop strip method proposed in this paper is 1.63 m, which is only 0.19 m away from the measured value, indicating that the formula predicts accurately.
5. Conclusions and Discussion
(1)According to the mechanical characteristics of a coal wall with a large mining height workface, the mechanical model of the coal wall caving is established using the Bishop strip method, and the solution equation of the safety factor K of the coal body on the slip surface is obtained.(2)In up-dip workface, the smaller the coal wall caving depth is, the smaller the coal wall safety factor is and the more unstable the coal wall is.(3)Taking workface 8102 of Wulonghu Coal Mine as an example, the safety factor, K, of the assumed slip surface is calculated and the critical slip surface position of slip coal caving is obtained. As a result, the maximum slip caving depth is 2.2 m, which is closer to the measured results than the calculation results of the existing compression-shear-type failure equation and sliding-type failure equation. This paper mainly addresses the slip-type coal caving that is common in fully mechanized mining having large mining heights, so these research results cannot explain the phenomenon of swallow-tail-type coal caving and so on. Otherwise, in the derivation of the safety factor equation of the slip surface, the influences of the initial support force and the thrust of the sidewall guard are not taken into account. In the future, research in these areas needs to be further strengthened.
Data Availability
Our data are from the filed case. The calculation parameters are based on geological data of the 8102 workface in the Wulonghu Coal Mine of China. All initial data of the workface are described in detail in Section 4 of this paper.
Conflicts of Interest
The authors have no conflicts of interest.
Authors’ Contributions
C. L. completed the establishment of the Bishop mechanical model of coal wall sliding instability and the derivation of the coal wall safety factor, K. In the case analysis, X. L. calculated the depth of the coal wall caving depth and compared it with the traditional calculation method. T. K. put forward the innovation and research direction of the article and played a corrective role in the whole process of the paper. All authors gave final approval for publication.
Acknowledgments
This study was supported financially by the National Natural Science Foundation of China (U1810102).