Abstract

Coal entry heading is one of the most hazardous activities in coal mine operations because a certain area of an unsupported roof inevitably forms and poses a significant threat to the safety of miners. In order to accelerate the coal entry heading, a simplified method including theoretical analysis and laboratory and in situ tests was developed to predict the influence of heading rate on the stability of the unsupported immediate roof. The results demonstrate that the deflection of the unsupported immediate roof at the heading face is on a scale of millimetre; hence, monitoring the deformation by conventional observation methods is difficult. The proposed model shows that, within the unsupported immediate roof, the peak values of normal stresses σx (perpendicular to the direction of excavation) and σy (parallel to the direction of excavation) and shear stresses τxz (perpendicular to the direction of excavation) and τyz (parallel to the direction of excavation) have different changing trends. The peak values of σx and σy both rise with the increasing advancing distance; however, σy reaches the tensile strength within a shorter range than σx. Moreover, the peak values of τxz and τyz initially increase with the increasing advancing distance and then stabilize or decline. The major threat to roof stability at the heading face is tensile failure parallel to the heading direction. According to the industry practices, it is proved that our method can make a good prediction of the mechanical state of the unsupported immediate roof, further deriving the heading rate with a considerable safety margin.

1. Introduction

As the biggest coal-producing and -consuming country in the world, more than 70% of China’s energy is provided by burning coal [1, 2]. Longwall mining is the main coal-extracting method used in China, with over 85% of the coal produced by longwall mining. Since the 1990s, with the development of mechanized coal mining, the “super longwall panel” that can produce 1–10 million tonnes of raw coal per year has been commonly implemented in China. With this technology, the advancing rate of a longwall panel can exceed 3000 m/year [3, 4] and will continue to increase. The rapidly advancing longwall panel will inevitably destroy a large part of the coal entry. Therefore, to avoid expensive production interruptions, Chinese coal mine managers and engineers are examining ways to apply the highly efficient heading of coal entry.

Highly efficient heading can be achieved by balancing speed with safety. In China, the most common heading technique utilizes a roadheader combined with a roof bolter (Figure 1). After excavating an advancing distance L0, the roadheader retreats about 5 m from the coal face; then, the roof bolter is moved close to the coal face, and the miners insert bolts into the rock to support the newly exposed roof. This frequent position exchange is time-consuming and reduces the heading speed greatly. Although increasing L0 can lower the frequency of the operation, it increases the risk that the longer unsupported roof will trigger roof collapse which poses a hazard to miners and equipment. Therefore, for highly efficient coal entry heading, we need to know the safe advancing distance L0 that can maintain roof stability as the roadheader moves forward.

Extensive investigation was conducted to prevent roof collapses in underground coal mines. Several geotechnical variables that affect roof stability have been identified. These include the geology, mine opening geometry, horizontal and vertical stress regime, abutment pressure, and support [512]. Most studies of mine roof stability try to determine the relationship between roof failures and the abovementioned variables by using statistical analysis [1319]. The detailed measurements help us understand the mechanics of roof instability and failure; however, the stratigraphy around coal mines can be complex, and the mining parameters vary between mines. Hence, the results of the above studies may be site-specific and not applicable to other mines.

Although the stability of the unsupported roof near the heading face is a vital factor for mining engineers when considering increasing the heading rate and preventing roof collapse, few studies have investigated this issue in depth. Hence, in this study, we develop an accurate model to predict the mechanical behavior of the unsupported roof. To improve the performance and effectiveness of the constructed model, a series of tests are conducted to obtain the parameters used in the model. The model is applied to a real coal excavation case, verifying its reliability as a tool for increasing mining safety and achieving highly efficient heading operation.

2. Problem Definition

The bearing structure of a typical heading face is illustrated in Figure 2. After excavation, the strong main roof supports the overlaying strata and maintains self-stability. Below the main roof, the weak strata (usually mudstone or siltstone) bend downward, causing separation of the weak strata from the main roof. At this stage, most of the weak strata are supported by the immediate roof. Hence, excavation-induced roof failure will initially begin at the immediate roof surface where the stress is concentrated and then will propagate upward into the deeper roof. This causes a huge loss of bearing capacity at the immediate roof and increases the likelihood of triggering roof collapse. Thus, the roof stability near the heading face is directly related to the stress regime and the engineering characteristics of the immediate roof.

The unsupported immediate roof can be simplified as a rectangular plate model (Figure 3) whose three edges are hinged at the rock stratum and one edge is hinged at the row of bolts next to the heading face. Previous studies generally regarded the area of the plate as the area of the exposed roof [2024]. This conception assumes that the coal ribs are rigid. However, in fact, the coal seam is much weaker than the roof and floor. The entry excavation causes the initial failure of the shallow coal rib, and it will gradually extend to the deep coal body. Hou and Ma [25] proposed a mechanical model (Figure 4) to calculate the normal stress distribution in the coal-roof interface of a coal rib. The model results indicate that the shallow coal rib loses most of its bearing capacity due to failures. This increases the area of the unsupported immediate roof. Therefore, in this work, the mechanical model (Figure 3) takes into account the roof over the broken coal rib.

The normal stress σz is obtained bywhere is the normal stress of the limit equilibrium zone in the coal rib; is the cohesion of the coal rib; is the internal friction angle of the coal rib; is the support resistance at the coal-rib edge; is the ratio of horizontal stress to vertical stress that can be obtained by in situ measurements; is the thickness of the coal seam. In Figure 4, α0 is the width of the limit equilibrium zone; H is the mining depth; is the average unit weight of the overlying strata which can be obtained by in situ measurements; and is the stress concentration factor, and the proposed value is 2 [25].

A realistic value for the width of the broken coal rib (l in Figure 4) is important for ensuring that the analysis returns valid results, but there is almost no reported in situ compressive test on coal ribs. Nevertheless, available results of extensive in situ compressive tests on coal pillars can provide valuable information for our modeling. Three sets of existing data on large-scale coal tests collected from References [2628] are depicted in Figure 5. After the peak strength, the curves eventually switch to a stable stage which represents the residual strength. The tests conducted are reviewed and listed in Table 1. The ratio of residual strength to peak strength ranges from 0.15 to 0.38 with the average value of 0.25. Therefore, here, a value equal to a quarter of the peak strength of the coal rib is used as the critical stress σb to determine the width of the broken coal rib l. By substituting σb into Equation (1), the width l can be expressed aswherein which is the peak strength of the small coal specimen and is the coefficient of the scale effect. Thus, the lengths and of the model shown in Figure 3 can be expressed aswhere denotes the width of the heading entry.

3. Model Solution

The elastic thin plate theory [29] is used to solve the rectangular model proposed in Section 2. It avoids many complex settings and produces results that meet practical engineering requirements. The deflection of the model in Figure 3 can be expressed aswherein which denotes the load concentration. represents the bending rigidity and can be expressed aswhere is the elastic modulus, is the plate thickness, and is Poisson’s ratio.

The bending moments Mx and My and the shear forces Qx and Qy are given as

The maximum values of the normal stresses σx and σy appear on the surface of the plate, while the maximum values of the shear stresses τxz and τxz appear on the middle plane of the plate. They are given as

The loads applied on the mechanical model depicted in Figure 3 are the overburden pressure and the coal-rib support pressure. Solving the model using the above equations directly is complicated and difficult. However, according to the superposition principle, the overall load q is equivalent to the superposition of the loads illustrated in Figure 6.

The overburden pressure q0 is evenly distributed and can be expressed aswhere is denoted as the average unit weight of the overlying strata and denotes the thickness of the immediate roof.

Compared to the overburden pressure q0, the distribution range of the coal-rib support pressures qL, qR, and qF is narrow. Hence, to simplify the calculation, the expressions of qL, qR, and qF are simplified as linear relations, yielding

Substituting Equations (11)–(14) into Equation (7) yields

Then, inserting Equations (15)–(18) into Equations (5) and (6) gives , , , and . Finally, the deflection of the model in Figure 3 can be expressed as

4. Case Study

A real coal entry was selected as a case study to test and verify the theoretical model presented above. The example is tailgate entry no. 12311 serving for coal seam no. 11-2 extraction in the Panyi Mine, Huainan Mining Industry Group Co., Ltd., Anhui Province, China. The site was chosen because the Panyi Mine is suffering the extremely low-efficiency heading of coal entry, with an average heading rate of 220 m/month. The main reason for the slow rate is the short advancing distance (an average of 1 m) and the corresponding frequent position exchange of the machinery.

4.1. Geological Setting

Coal seam no. 11-2 is buried at an average depth of 760 m, and its average thickness is 2.4 m. Tailgate entry no. 12311 is cut as 2.1–2.4 m high and 5 m wide and is oriented N83°W. A geologic column of the mine roof obtained from a vertical core hole drilled at the site is shown in Figure 7. The roof may be roughly divided into three units:(i)A slightly stronger siltstone layer in the lowest 1.1 m, comprising the immediate roof(ii)A sequence of coal, weak claystone, and slickensided shale layers from 1.1 m to 14.4 m(iii)A significantly stronger sandstone above 14.4 m, comprising the main roof

The properties of the coal ribs for the model in Section 3 were obtained from previous tests on coal from the adjacent Xieyi Mine [30]. The coal specimens taken from the Xieyi Mine were subjected to multistage triaxial testing, and the results are presented in Table 2. By using the linear Mohr–Coulomb criterion, the cohesion and friction angle of the coal specimens were determined as 7.92 MPa and 27°, respectively. The uniaxial strength and cohesion of the coal specimens cannot be directly recognized as the in situ properties of the coal ribs. However, by in situ large-scale coal testing, Bieniawski [31] found that the coefficient of the scale effect can be set as 0.17. Accordingly, the in situ compressive strength and cohesion of the coal ribs used in the model were 3.91 MPa and 1.35 MPa, respectively.

As illustrated in Figure 2, our model ignored the bearing capacity of weak formations between the main roof and immediate roof and regarded them as the load applied on the immediate roof. By doing this, the calculation process has been greatly simplified; meanwhile, the safety margin of results can be increased. In order to understand the mechanical properties of the immediate roof of tailgate entry no. 12311, a series of tests were conducted on siltstone specimens obtained from the roof core (Figure 8): uniaxial compressive tests, Brazilian disc tests, and shear-compression tests. The results are presented in Tables 35, respectively. The relationship between the average shear strength and normal stress is illustrated in Figure 9. Based on the Mohr–Coulomb criterion, the cohesion and friction angle of the siltstone were derived as 4.86 MPa and 35.6°, respectively.

The stress field was measured in two locations in haulage entry no. 12311 and adjacent recovery room no. 12521 (Figure 10). The results are shown in Table 6. By decomposing the principal stresses and extracting the stress vector in the direction perpendicular to tailgate entry no. 12311, the coefficient of horizontal stress A can be obtained as 0.9. The average unit weight of the overlying strata γ can be expressed aswhere is the average vertical stress and is the mining depth. Substituting the measurements in Table 5 into Equation (20), is obtained as 26836 N/m3.

4.2. Results and Discussion

In Section 4.1, the detailed parameters describing the geological conditions of tailgate entry no. 12311 were determined by experimental and in situ tests. To study the roof stability at the heading face in tailgate entry no. 12311, we substituted these parameters into the model proposed in Section 3. Given the complexity of solving the equations, the commercial math software MATLAB R2012a (MathWorks, Inc., Natick, USA) was used.

As the first step to explore the possibility of increasing the advancing distance, the stress state of the current unsupported immediate roof (1 m long and 5 m wide) was analyzed. The distributions of the deflection, maximum normal stresses σx and σy, and maximum shear stresses τxz and τyz in the immediate roof are illustrated in Figure 11. The deflection is on a scale of decimillimetre; hence, it cannot be perceived by miners working in the heading face. The distribution of σx forms a hump-shaped pattern, with the two peaks located on the midline of the unsupported roof (Figure 11(b)). But the maximum value is only 0.137 MPa which is far below the tensile strength of the immediate roof (2.18 MPa). At this low stress level, the lateral broken coal ribs support the unsupported roof to some extent, which produces the bilateral compressive areas presented in Figure 11(b). Compared to the maximum normal stress σx, σy (Figure 11(c)) is relatively large, reaching a peak value of 0.486 MPa at the center of the unsupported roof. As the stress level rises, the supporting effect of the broken ribs weakens. The distribution of the maximum shear stress τxz (Figure 11(d)) shows two high-stress areas which are asymmetric, both along the lateral coal ribs. The maximum shear stress τyz (Figure 11(e)) has two high-stress areas at the front and back edges. The stress along the back edge which is corresponding to the row of bolts next to the heading face is more concentrated. The peak values of τxz and τyz are 0.152 MPa and 0.373 MPa, respectively. They are both small compared with the shear strength of the immediate roof. Overall, the results indicate that the midline, center, bilateral edges, and back edge of the unsupported immediate roof are possible failure areas, and more importantly, they confirm that the stress state of the current immediate roof is below the threshold level needed to trigger the roof failure; thus, there is a large potential to increase the advancing distance.

To investigate the stress evolution characteristics of the unsupported immediate roof with the increasing advancing distance, we conducted a serial of calculations. The results were plotted against advancing distance and are shown in Figure 12. The peak value of deflection shows a nonlinear increase with the increasing advancing distance, but it is still on a scale of millimetre which is hard to measure or observe. The normal stress and shear stress show clearly different trends. Specifically, the peak values of σx and σy both rise with the increasing advancing distance, but σy reaches the tensile strength within a shorter range than σx. Moreover, the peak values of τxz and τyz initially increase with the increasing advancing distance and then become stable or decline. The simulation was designed to test whether the shear stress can reach the minimum shear strength determined as 4.86 MPa in Section 4.1. However, as depicted in Figures 12(d) and 12(e), it is far below the minimum shear strength even at an unrealistic length of advancing distance (29 m).

According to the above analysis, it is apparent that the major threat to roof stability at the heading face is tensile failure parallel to the heading direction. However, a problem arises from the determination of the tensile strength of the immediate roof. Since the coal mine roof is commonly disturbed by faults, beddings, and joints, its strength cannot be equated to the strength derived from the small rock sample. Therefore, a reduction parameter is required to relate the roof strength to the rock strength in the laboratory. So far, extensive works have been devoted to this issue; several rock mass classification systems (for e.g., RMR [32], MRMR [33], Q-system [34], and GSI [35]) have been established and successfully applied in engineering practice, but the fairly high time cost is a major shortcoming of these classification systems blocking their application to this study. As mentioned, the main intention of the establishment of the foregoing analytical method is to increase the efficiency of coal entry heading. However, the previous rock mass classification systems generally need a variety of input parameters (joint spacing, joint roughness, water reduction factor, etc.), some of which are based on the time-consuming in situ observation and borehole logging. Besides, these classification systems are not designed to rate the roof impacted by a variable geological condition; hence, parameter adjustments are required to be continuously conducted with the advance of heading, which would greatly slow down the heading rate. In view of the problem, a simplified but competent relationship between roof strength and laboratory rock strength is necessary to allow the further analysis.

The existing empirical relations of in situ deformation modulus and rock quality designation (RQD) paved a way for overcoming the problem because of the following reasons:(1)The positive correlation between rock mass strength and deformation modulus has been proven by numerous researchers [3638]; therefore, it becomes convenient to relate the rock mass strength to RQD.(2)Compared to the rock mass strength, the dataset of the deformation modulus of rock mass is much larger, which can be attributed to the utilization of the seismic technique in its determination. A larger dataset indicates a more reliable relation.(3)RQD is a single parameter and can be easily calculated from core drilling; this simplicity makes it possible to evaluate the variable condition of the immediate roof.

Zhang and Einstein [39] collected about 120 sets of published data covering mudstone, siltstone, sandstone, shale, dolerite, granite, limestone, greywacke, gneiss, and granite gneiss and derived an empirical relation between the modulus reduction ratio Em/Ei and the RQD:where and are the deformation moduli of the rock mass and the intact rock, respectively. The reduction ratio of deformation modulus is employed to relate the roof strength to the laboratory-scale rock strength. Based on RQD, the immediate roof conditions are classified into four types and presented in Table 7. To increase the safety margin, the reduction ratio () of each type is calculated by substituting the lower value of RQD into Equation (21). Subsequently, is introduced into the calculation of the advancing distance in heading of tailgate entry no. 12311. The critical tensile stress below which the immediate roof can maintain self-stability is given aswhere is the tensile strength of the small sample drilled from the immediate roof and has been determined as 2.18 MPa in Table 4. The determined of each roof condition type is presented in Table 7.

It can be observed in Figure 12 that the tensile stress parallel to the heading direction (σy) exhibits the fastest increase with the increasing advancing distance. Therefore, the advancing distance of each type of roof condition is determined by substituting their critical tensile stress into the linear relationship presented in Figure 12(c). It has to be noted that the advancing distance of the roof condition with RQD lower than 70% is kept as the original value (1 m). Before applying the advancing distances to the entry heading, the expected mechanical state of the unsupported immediate roof was studied (an example of 1.9 m advancing distance is shown in Figure 13). The results confirm that the stress distributions are all within a safe level and indicate that the advancing distance in the modeled coal seam of the Panyi Mine can be confidently increased.

The expected length of tailgate entry no. 12311 is 1370.5 m. Before employing the newly designed advancing distance, the heading face had advanced 370 m taking 55 days. Thus, the average heading rate was 202 m/month. After adopting the newly designed advancing distance, the remaining 1000.5 m was excavated in 74 days with no roof falls. Thus, even with a considerable safety margin, the heading rate was improved from 202 m to 406 m per month. The frequency of position changes between the roadheader and the roof bolter decreased from 6.7/day to 5.4/day, which reduced the labor intensity. Finally, the proposed method was verified as a reliable method for predicting the mechanical state of the unsupported immediate roof, providing a safe and effective tool to help achieve highly efficient heading of coal entry.

5. Conclusions

This study focuses on developing a systematic method to predict the mechanical behavior of the unsupported immediate roof and preventing roof collapse while pursuing the highly efficient heading of coal entry. The theoretical model was applied to data collected from a Chinese coal mine. Based on the results, the following conclusions were drawn:(1)The deflection of the unsupported immediate roof at the heading face was on a scale of millimetre; hence, it is difficult to monitor the deformation of the unsupported immediate roof using conventional observation methods. It is thus extremely dangerous to determine the extent of safe advancing distance by observing the deformation of the unsupported immediate roof.(2)The proposed model shows that the peak values of normal stress and shear stress in the unsupported immediate roof show different changing trends. The peak values of σx and σy both rise with the increasing advancing distance, while σy reaches the tensile strength within a shorter range than σx. Moreover the peak values of τxz and τyz initially increase with the increasing advancing distance and then stabilize or decline.(3)The peak values of τxz and τyz cannot reach the minimum shear strength over the range of advancing distance from 1 m to 29 m. Therefore, shear failure is unlikely to occur in the unsupported immediate roof.(4)The major cause of roof instability at the heading face is tensile failure parallel to the heading direction.(5)Based on engineering practice, the theoretical method was verified as a reliable model for predicting the mechanical state of the unsupported immediate roof and can be used confidently to achieve safe and highly efficient heading of coal entry.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities (2018BSCXC31) and Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX18_1963).