This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
With increased mining depth, the failure of protective pillars in deep high-stress roadways has received much attention [1]. The influence of the protective pillar design on the stability of roadways has a bearing on worker safety and the effective mining of coal resources. Protective pillar around the gob provides a safe and stable production environment for the roadway and decreases disturbance to the roadway from high external stress. An improper protective pillar design can cause frequent roof sagging, pillar rib bulging, and severe coal bump to occur in a roadway, which will eventually cause the collapse of the surrounding rock structure [2, 3]. Therefore, protective pillar width design is a necessity for maintaining the stability of the surrounding rock. Many methods of designing an appropriate protective pillar width have been proposed, including theoretical analysis calculations and numerical simulations.
A variety of calculation models have been established to analyse the proper width of a protective pillar and decrease the risk of pillar failure [4, 5]. Ghasemi et al. [6] used limit equilibrium analysis to explore the relationship between the range of elastic area and the stability of the protective pillar. Cao and Zhou [7] investigated the influence of movement in key blocks of broken strata above the roadway on the stability of the protective pillar and defined a reasonable width for the pillar. Yang et al. [8] showed the scope of stress distribution in the upper and lower strata of a protective pillar by analysing the influence of roof failure on the bearing capacity of the pillar.
In reality, many factors, such as the protective pillar’s form, the rock mechanical properties, the stress state, and the panel mining methods, should be considered in designing the width of the protective pillar. Since some of these factors are not considered in theoretical analysis calculations, they are limited in application. Various numerical simulations have been adopted to analyse the effect of the protective pillar width on the roadway stability because of their low cost and high efficiency [9]. Moreover, as numerical simulation includes numerous factors during modeling, the results are often more reliable than those of theoretical calculations. Shabanimashcool and Li [10] established local and global numerical models to study the influence of falling strata on the stress and breaking range of a protective pillar during panel retreat. Bai et al. [11] investigated roof failure characteristics and stress change rules of gob-side entry driving by stages, noting that a rational protective pillar width could enhance the roadway stability. Ma et al. [12] analysed the shear failure behavior for protective pillars and found that the failure mode varied with the pillar width and inclination of the coal seam. Wang et al. [13] adopted FLAC3D software to study the dynamic response and failure mechanism of protective pillars; the results showed that enlarging the elastic core in the pillar increased the risk of coal bump. However, based on previous studies, we found that the numerical simulations seldom considered the following: (1) When the overlying strata in a gob are cut off along a preexisting fissure, the caving materials will fill the gob and provide supporting resistance to the roof strata, which will relieve stress on the protective pillar and decrease its width. (2) The energy stored in the rock mass is closely related to its stress state. The stored energy will vary as the stress changes. Thus, the evolution laws of energy stored in a rock mass can be inferred by the change in stress, decreasing or avoiding the deformation failure of the surrounding rock and the occurrence of coal bump. In view of these limitations in current simulations and considering the supporting feature of gob caving materials and the relationship between the stress and energy of the rock mass, a numerical modeling is developed to analyse failure mechanism for protective pillars.
In this paper, the relationship between roadway stability and protective pillar width is investigated based on numerical simulations and field tests. First, the failure characteristics of high-stress roadways are studied when the width of protective pillar is 20 m. Next, a novel modeling approach is proposed. In this modeling, a strain-softening model is adopted to describe the mechanical behaviour of protective pillars, and a double-yield model is used to simulate the gob materials mechanical behaviour. The stress changes and energy density distribution characteristics of roadways with five different widths are analysed. Finally, the effect of the optimal protective pillar on controlling the deformation of the surrounding rock is evaluated by a field test. The modeling approach and design principle presented in this paper can be used to analyse the protective pillar design of high-stress roadways at other similar sites.
2. Engineering Background
2.1. Geology Conditions
As shown in Figure 1, the selected study site is in the Wangzhuang coal mine, Shanxi Province, China. The average depth and thickness of coal seam are 300 m and 4 m. Roadway roof and floor mainly include mudstone, siltstone, and sandy mudstone. A detailed stratigraphic column is shown in Figure 2.
[figure omitted; refer to PDF]
Panels 8101, 8102, and 8103 are located in mining area #8 of the test coal mine, with a strike length of 1088 m and dip length of 210 m, as shown in Figure 3. Panel 8101 is located to the north of panel 8102, where the coal seam has been exploited. Panel 8103 is located to the south of panel 8102, which is unexploited. At present, the width of the protective pillar between the 8101 tailgate and 8102 headgate is 20 m (Figure 3).
[figure omitted; refer to PDF]
For determining a rational pillar width for a high stress roadway, the following considerations could be made:
(1) The width of the protective pillar could not be too thin, such as the left side of point A in Figure 8. Because of its weak bearing capacity, the protective pillar could not maintain roadway stability in this case.
(2) The protective pillar width could not be too wide, indicated by the curve between points B and D. On account of the strain energy accumulation, the protective pillar may lead to coal bumps in this situation.
(3) The width of the protective pillar could not reach the critical width, such as point B. Under this condition, the peak stress would be shifted from the solid coal region to the protective pillar. The bearing strength of the protective pillar could not resist the peak stress effectively, and the roadway was in a state of high stress.
A suitable width for protective pillars should provide enough capacity to bear roof strata loads while maintaining the pillar integrity, reducing the pillar width, and increasing coal recovery rate. Therefore, a thin pillar could be chosen, corresponding to the left side of point B (Figure 8). The roadway could undertake a certain degree of large deformation in this condition, and its integrity could be maintained.
3.3. Determination of the Cutting Height for the Overlying Strata in the Gob
During mining, as the working face constantly moved forward, the strata behind the working face were cut off along the preexisting fissure. The gob caving material had a certain bearing capacity when compressed by the upper strata. The gob caving material could support the overlying strata and reduce the pressure on the pillars. Therefore, when designing protective pillars, it is necessary to account for the supporting characteristics of the caving materials.
To effectively utilize the supporting characteristics of the caving materials, it is important to determine the cutting height of the overlying strata in the gob. After roof stratum caving, the volume of caving materials expands. In other words, the pile height of the caving rock mass will be larger than the original height. The volume expansion of the caving materials can be fully described with the bulking factor
According to equation (1), when Δ = 0, the cutting height of the roof strata is calculated by the following equation [15, 16]:
Based on a number of field observations, the caving height is 4 to 8 times of the coal seam height [17, 18]. For panel 8102, the coal seam height is 4 m, and the bulking factor is determined to be 1.25 [19, 20]. Therefore, the cutting height of the strata above the gob is 16 m.
4. Establishment of Numerical Model
4.1. Modeling Scheme
To analyse the relationship between the roadway stability and protective pillar width, a numerical model was established based on FLAC3D software. The dimensions of the model were 160 m × 59.1 m × 2 m. The model consisted of panels 8101 and 8102 and their roadway system (Figure 9). A vertical stress of 7.5 MPa, representing the loads of overlying rock, is applied to the model top boundary. Based on the in situ stress monitoring data, the ratio of horizontal to vertical stress was determined to be 1.1. For the four vertical model planes, the horizontal displacement was constrained. The displacements in the horizontal and vertical directions were restrained in bottom boundary. The Mohr–Coulomb model was adopted to simulate the rock mass, except for the protective pillar and gob caving materials. Table 1 shows the mechanical properties of the rock mass.
[figure omitted; refer to PDF]
Table 1
Mechanical properties of the rock mass.
Lithology | Density (kg/m3) | Compressive strength (MPa) | Elasticity modulus (GPa) | Poisson's ration ( | Cohesion (MPa) | Friction (°) |
Fine sandstone | 2750 | 57.2 | 12.5 | 0.22 | 2.6 | 29 |
Mudstone | 1900 | 13.1 | 5.0 | 0.29 | 1.4 | 25 |
Sandy mudstone | 2450 | 27.8 | 7.1 | 0.26 | 1.8 | 27 |
Siltstone | 2680 | 47.6 | 9.8 | 0.24 | 2.3 | 31 |
Mudstone | 1900 | 13.1 | 5.0 | 0.29 | 1.4 | 25 |
Coal seam | 1600 | 9.8 | 1.5 | 0.32 | 0.8 | 22 |
Mudstone | 1900 | 13.1 | 5.0 | 0.29 | 1.4 | 25 |
Fine sandstone | 2750 | 57.2 | 12.5 | 0.22 | 2.6 | 29 |
Siltstone | 2680 | 47.6 | 9.8 | 0.24 | 2.3 | 31 |
The numerical simulation comprised five steps: (1) model establishment and application of in situ stress, (2) initial stress balance in the model, (3) 8101 tailgate excavation, (4) retreat of panel 8101 and equivalent gob materials simulation, and (5) 8102 headgate excavation with five different pillar widths. On the basis of the 8102 headgate conditions, the protective pillar height was maintained at 3.5 m, and the protective pillar width was simulated at 5, 10, 15, 20, and 25 m, as shown in Figure 10.
[figure omitted; refer to PDF]
Many strength formulas of pillar have been put forward in the recent decades [23–25]. It is noted that Salamon-Munro proposed the following empirical formula, with pillar W/H (width-to-height) ratios from 2 to 20 [23].
The empirical strength formula of the pillar has good consistency with developed average strength formulas. Therefore, the Salamon-Munro strength formula was adopted in this study to verify the pillar numerical model. A commonly used iterative method is adopted [21, 26]. The parameters of the strain-softening model are determined by matching the strength of the pillar acquired by numerical simulation with that determined by the Salamon-Munro formula. Table 2 shows the input parameters for calibrating the pillar strain-softening model. The simulated and calculated strengths of pillar are displayed in Figure 12. The results show that the simulated strength of the pillar matches with the calculated strength of pillar from the Salamon-Munro empirical formula very well, suggesting that the parameters can be adopted to accurately simulate the mechanical behaviour of a protective pillar.
Table 2
Strain-softening properties of pillar with plastic strain.
Strain | 0 | 0.0025 | 0.005 | 0.0075 | 0.01 |
Cohesion (MPa) | 0.8 | 0.68 | 0.54 | 0.4 | 0.28 |
Friction angle (°) | 24 | 23 | 22 | 21 | 21 |
4.3. Modeling of the Caving Materials in the Gob
After panel retreat, the roof strata behind the panel collapse, and the caving rock in the gob is compacted and consolidated. Regarding the influence of the gob caving rock on the protective pillars, the mechanical behaviour of the gob caving materials was simulated by the double-yield model [21, 27]. The input parameters of the double-yield model included material properties and cap pressure [28]. The cap pressure parameters were determined by the following equation proposed by Salamon:
Table 3
The double-yield model cap pressures.
Strain | 0 | 0.01 | 0.02 | 0.03 | 0.04 | 0.05 | 0.06 | 0.07 | 0.08 | 0.09 |
Stress (MPa) | 0 | 0.26 | 0.55 | 0.88 | 1.24 | 1.65 | 2.13 | 2.67 | 3.31 | 4.06 |
Strain | 0.10 | 0.11 | 0.12 | 0.13 | 0.14 | 0.15 | 0.16 | 0.17 | 0.18 | 0.19 |
Stress (MPa) | 4.96 | 6.06 | 7.44 | 9.21 | 11.57 | 14.88 | 19.84 | 28.10 | 44.64 | 94.24 |
To obtain the material parameters of the double-yield model (bulk modulus, shear modulus, friction angle, and dilation angle), a gob element (Figure 13) was selected. A 10−5 m/s vertical velocity was used to represent the load and apply to the element top. The horizontal deformation of four vertical faces of the element was restrained, while the vertical deformation at the element bottom was set to zero. A fitting method in previous related research was adopted to match the gob element stress-strain curve with Salamon’s model [21, 29]. Figure 14 compares the stress-strain curves of two models, and they match very well. Based on the comparison result, the final parameters of the gob material in the double-yield model were obtained and used in FLAC3D (Table 4).
[figure omitted; refer to PDF]
On the other hand, it can be considered that stresses
Under the action of the average stress
Under the action of the deviatoric stresses
The element strain energy density includes the volume-changed and shape-changed strain energy density. Therefore, the strain energy density of the element can be expressed [32, 33]:
According to equation (11), it can be clearly seen that energy stored in the surrounding rock element is closely related to its own stress state. A coupled program of the strain energy density and stress of a surrounding rock element was developed in this study, which was used to investigate the relationship between the internal energy and stress of the rock mass.
4.5. Verification of the Gob Caving Materials Supporting Capacity
To validate the effectiveness of the supporting capacity of caving materials in gob 8101, Figure 16 displays the vertical stress evolution law in the caving area. There is 0.42 MPa of vertical stress at the gob edge. Meanwhile, the vertical stress increases gradually with distance increase from the gob edge. The maximal vertical stress reaches 7.45 MPa at an approximately 73.6 m distance from the gob edge and remains relatively steady. The caving materials in the gob bear 99% of the initial vertical stress (7.45/7.5 MPa). With increasing distance from the gob edge, the caving materials are compressed and the vertical stress in the caving area increases gradually. The vertical stress of the gob reaches the initial stress (7.5 MPa) at approximately 73.6 m from the gob edge, that is, 25% (73.6/300) of the overlying depth. Through extensive field tests, Wilson and Carr [34] and Campoli et al. [35] proposed that the gob vertical stress could reach the initial stress of the rock mass at a distance of 0.2–0.38 times the buried depth of coal seam. In this study, the obtained vertical stress distribution agrees well with the related conclusions [14, 25, 36], indicating that the simulation result of the supporting capacity of the gob caving materials is reliable.
[figure omitted; refer to PDF]
The results (Figures 17 and 18) show that the protective pillar width influences greatly the pillar load-bearing capacity and roadway deformation. Most surprisingly, when the protective pillar widths are 10 m and 25 m, the difference in roadway deformation is not obvious. This demonstrates that the existence of an elastic intact zone enhances the stability of the protective pillar. Consequently, increasing the protective pillar width could enhance the reliability of the pillar load-bearing capacity and improve the roadway stability.
5.2. Discussion of the Rational Width for a Protective Pillar
The original balanced state of the main roof above the roadway is strongly disturbed during the panel retreat period. The main roof is broken [37–39], as shown in Figure 19(a). Rock masses A, B, and C interact, which influences the stability of the 8102 headgate. During the roof fracture period, rock mass B above the roadway begins to rotate and subside. The rotation and subsidence have a great impact on the stability of the protective pillar. The rock B fracture position and rotating speed are closely related to the pillar width. According to the ultimate balance theory, the geometry size
[figures omitted; refer to PDF]
From the results of the numerical simulation, in the case of a 5 m wide protective pillar, the peak vertical stress of the virgin coal rib occurs 14.3 m from the rib edge. The length between the peak vertical stress location in the virgin rib and the rib edge plus the width of the protective pillar and the roadway equals the geometry size of rock B, indicating that rock mass B will fracture at the location where the peak vertical stress of the virgin coal rib occurs (Figure 19(a)). A 5 m wide protective pillar cannot withstand the loads of the overlying strata. Rock mass B rotates violently and the protective pillar is crushed. With a 20 m wide protective pillar, the vertical stress peak in the virgin rib is transferred to the pillar. Based on the geometry size of rock B, it could be concluded that its fracture position is located at the junction of the roof and the virgin coal rib (Figure 19(b)), which is consistent with the field result showing that roof sagging and falling often occur at the place close to the junction between the roof and the virgin coal rib. In the panel retreat period, the superimposed effect of mining stress and high in situ stress surpasses the bearing strength of the protective pillar, and the continued pillar deformation leads to the fracture and rotation movement of rock mass B in the main roof. Finally, the protective pillar is crushed and fails. With a 10 m wide protective pillar, the vertical stress peak is still located in the virgin coal rib; that is, the area of virgin rib is subjected to the main load of the overlying strata. The pillar bears relatively less load, which puts the roadway in a state of low stress. Meanwhile, the protective pillar could provide support to the fractured rock mass B and prevent it from rotary movement, which maintains the stability of the roadway. When the width of protective pillar is 25 m, the bearing capacity is improved due to the expansion of the elastic area in the pillar. The protective pillar can not only maintain its own stability but also prevent the fractured rock mass B of the main roof from rotating, transitioning the roadway to a stable state. However, a greater pillar width wastes more coal resources.
A rational pillar width should have the capacity to withstand abutment pressure and maintain a stable roadway, while maximizing the recovery of coal resources. Based on the above discussion, the recommended 10 m wide protective pillar may be more conducive to improving the roadway stability and coal recovery rate, while reducing economic loss due to large deformation and coal bump failure.
6. Field Test
A protective pillar width of 10 m was applied in field to maintain the stability of the surrounding rock 8103 headgate. The protective pillar stress distribution and deformation of the 8103 headgate were measured to validate the reliability of the simulation results.
6.1. Stress Measurement of the Protective Pillar in the Field
The vertical stress at different distances from the protective pillar rib of the 8103 headgate is plotted in Figure 20. The field data are also displayed for comparison with the modeling results. The measured results and simulated results have good agreement, verifying the numerical modeling accuracy.
[figures omitted; refer to PDF]
6.2. Roadway Deformation Laws
The measured roadway deformation results at different stages are shown in Figure 21. The deformation of the roadway became stable after 50 days of excavation. The convergences of the roof, virgin coal rib, protective pillar rib, and floor were 115 mm, 74 mm, 87 mm, and 16 mm, respectively. During the panel 8103 retreat period, the roadway deformation mainly occurred 60 m in front of the working face; the total convergences at the roof, virgin coal rib, protective pillar rib, and floor were 301 mm, 162 mm, 214 mm, and 24 mm, respectively, which represented reductions of 70%, 59%, 68%, and 77% (Figure 22), compared with the deformation that occurred under the previous protective pillar width. The stability of the surrounding 8103 headgate is improved, as shown in Figure 23. The field observation results indicate that the proposed protective pillar width is feasible for maintaining the stability of the roadway.
[figures omitted; refer to PDF]
[figure omitted; refer to PDF][figure omitted; refer to PDF]7. Conclusions
This study was mainly focused on analysing the effect of protective pillar width on the roadway stability to identify a pillar design principle based on field tests and numerical simulations. The main conclusions in this paper are summarized as follows:
(1) A novel numerical model was established to analyse the failure mechanism of the protective pillar. As an innovative design method of protective pillars, the supporting features of the gob caving materials on overlying strata and the relationship between the internal stress and stored energy of the rock mass were considered in the modeling process. The stress change and distribution characteristics of the energy density are regarded as important conditions in designing the width of a protective pillar and evaluating the stability of the roadway.
(2) The modeling results showed that, with a 20 m wide pillar, the peak vertical stress and energy density in the protective pillar were 18.5 MPa and 563.7 kJ/m3, respectively. Excessive stress and elastic energy resulted in considerable deformation and coal bump failure in headgate 8102. With a pillar width of 10 m, the location of the peak vertical stress and energy density moved from the protective pillar to the virgin rib. The main loads of the overlaying strata were borne by the virgin rib, and the pillar was subjected to a relatively low load. The roadway was in a state of low stress and could maintain its stability.
(3) The results of the field measurements showed that a 10 m wide protective pillar was able to effectively control the roadway deformation and release most of the storing energy in protective pillar. Meanwhile, the proposed modeling approach and protective pillar design principles in this study can provide a useful basis for application in similar coal mines.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (no. 51674250) and the State Key Program of the National Natural Science Foundation of China (no. 50834005). In addition, the authors are very grateful for the linguistic assistance of AJE (American Journal Experts) in the manuscript preparation period.
[1] S.-Q. Yang, M. Chen, H.-W. Jing, K.-F. Chen, B. Meng, "A case study on large deformation failure mechanism of deep soft rock roadway in Xin'An coal mine, China," Engineering Geology, vol. 217, pp. 89-101, DOI: 10.1016/j.enggeo.2016.12.012, 2017.
[2] R. Gu, U. Ozbay, "Numerical investigation of unstable rock failure in underground mining condition," Computers and Geotechnics, vol. 63, pp. 171-182, DOI: 10.1016/j.compgeo.2014.08.013, 2015.
[3] Z. Zhang, M. Deng, X. Wang, W. Yu, F. Zhang, V. D. Dao, "Field and numerical investigations on the lower coal seam entry failure analysis under the remnant pillar," Engineering Failure Analysis, vol. 115,DOI: 10.1016/j.engfailanal.2020.104638, 2020.
[4] W. Gao, M. Ge, "Stability of a coal pillar for strip mining based on an elastic-plastic analysis," International Journal of Rock Mechanics and Mining Sciences, vol. 87, pp. 23-28, DOI: 10.1016/j.ijrmms.2016.05.009, 2016.
[5] H. Wu, X. Wang, W. Wang, G. Peng, Z. Zhang, "Deformation characteristics and mechanism of deep subsize coal pillar of the tilted stratum," Energy Science & Engineering, vol. 8 no. 2, pp. 544-561, DOI: 10.1002/ese3.546, 2020.
[6] E. Ghasemi, M. Ataei, K. Shahriar, "Prediction of global stability in room and pillar coal mines," Natural Hazards, vol. 72 no. 2, pp. 405-422, DOI: 10.1007/s11069-013-1014-2, 2014.
[7] Z. Z. Cao, Y. J. Zhou, "Research on coal pillar width in roadway driving along goaf based on the stability of key block," CMC-Computers, Materials Continua, vol. 48, pp. 77-90, 2015.
[8] J. X. Yang, C. Y. Liu, B. Yu, F. F. Wu, "The effect of a multi-gob, pier-type roof structure on coal pillar load-bearing capacity and stress distribution," Bulletin of Engineering Geology and the Environment, vol. 74 no. 4, pp. 1267-1273, DOI: 10.1007/s10064-014-0685-6, 2015.
[9] M. Chen, S.-Q. Yang, Y.-C. Zhang, C.-W. Zang, "Analysis of the failure mechanism and support technology for the Dongtan deep coal roadway," Geomechanics and Engineering, vol. 11 no. 3, pp. 401-420, DOI: 10.12989/gae.2016.11.3.401, 2016.
[10] M. Shabanimashcool, C. C. Li, "A numerical study of stress changes in barrier pillars and a border area in a longwall coal mine," International Journal of Coal Geology, vol. 106, pp. 39-47, DOI: 10.1016/j.coal.2012.12.008, 2013.
[11] J.-B. Bai, W.-L. Shen, G.-L. Guo, X.-Y. Wang, Y. Yu, "Roof deformation, failure characteristics, and preventive techniques of gob-side entry driving heading adjacent to the advancing working face," Rock Mechanics and Rock Engineering, vol. 48 no. 6, pp. 2447-2458, DOI: 10.1007/s00603-015-0713-2, 2015.
[12] T. Ma, L. Wang, F. T. Suorineni, C. Tang, "Numerical analysis on failure modes and mechanisms of mine pillars under shear loading," Shock and Vibration, vol. 2016,DOI: 10.1155/2016/6195482, 2016.
[13] H. Wang, Y. Jiang, Y. Zhao, J. Zhu, S. Liu, "Numerical investigation of the dynamic mechanical state of a coal pillar during longwall mining panel extraction," Rock Mechanics and Rock Engineering, vol. 46 no. 5, pp. 1211-1221, DOI: 10.1007/s00603-012-0337-8, 2013.
[14] Z. Zhang, J. Bai, Y. Chen, S. Yan, "An innovative approach for gob-side entry retaining in highly gassy fully-mechanized longwall top-coal caving," International Journal of Rock Mechanics and Mining Sciences, vol. 80,DOI: 10.1016/j.ijrmms.2015.09.001, 2015.
[15] V. Palchik, "Bulking factors and extents of caved zones in weathered overburden of shallow abandoned underground workings," International Journal of Rock Mechanics and Mining Sciences, vol. 79, pp. 227-240, DOI: 10.1016/j.ijrmms.2015.07.005, 2015.
[16] F. Z. Qi, Z. G. Ma, "Investigation of the roof presplitting and rock mass filling approach on controlling large deformations and coal bumps in deep high-stress roadways," Latin American Journal of Solids and Structures, vol. 190, 2019.
[17] S. S. Peng, Longwall Mining, pp. 46-53, 2006.
[18] M. G. Qian, W. C. Shi, J. L. Xu, Ground Pressure and Strata Control, 2010.
[19] B. X. Huang, Y. Z. Wang, S. G. Cao, "Cavability control by hydraulic fracturing for top coal caving in hard thick coal seams," International Journal of Rock Mechanics and Mining Sciences, vol. 74, pp. 75-57, DOI: 10.1016/j.ijrmms.2014.10.011, 2015.
[20] S. Yan, T. Liu, J. Bai, W. Wu, "Key parameters of gob-side entry retaining in a gassy and thin coal seam with hard roof," Processes, vol. 6 no. 5, pp. 51-65, DOI: 10.3390/pr6050051, 2018.
[21] W. Li, J. Bai, S. Peng, X. Wang, Y. Xu, "Numerical modeling for yield pillar design: a case study," Rock Mechanics and Rock Engineering, vol. 48 no. 1, pp. 305-318, DOI: 10.1007/s00603-013-0539-8, 2015.
[22] Z. Z. Zhang, W. J. Wang, S. Q. Li, "An innovative approach for gob-side entry retaining with thick and hard roof: a case study," Technical Gazette, vol. 25 no. 4, pp. 1028-1036, 2018.
[23] M. D. G. Salamon, A. H. Munro, "A study of the strength of coal pillars," Journal of the Southern African Institute of Mining and Metallurgy, vol. 68, pp. 55-67, 1967.
[24] X. Du, J. Lu, K. Morsy, S. S. Peng, "Coal pillar design formulae review and analysis, " , pp. 153-160, .
[25] G. Feng, P. Wang, Y. P. Chugh, "Stability of gate roads next to an irregular yield pillar: a case study," Rock Mechanics and Rock Engineering, vol. 52 no. 8, pp. 2741-2760, DOI: 10.1007/s00603-018-1533-y, 2019.
[26] G. Zhang, S. Liang, Y. Tan, F. Xie, S. Chen, H. Jia, "Corrigendum to: numerical modeling for longwall pillar design: a case study from a typical longwall panel in China," Journal of Geophysics and Engineering, vol. 16 no. 3,DOI: 10.1093/jge/gxz042, 2019.
[27] M. Shabanimashcool, C. C. Li, "Numerical modelling of longwall mining and stability analysis of the gates in a coal mine," International Journal of Rock Mechanics and Mining Sciences, vol. 51, pp. 24-34, DOI: 10.1016/j.ijrmms.2012.02.002, 2012.
[28] M. Salamon, "Mechanism of caving in longwall coal mining,” Rock mechanics contributions and challenges," Proceedings of the 31st US Symposium on Rock Mechanics, pp. 161-168, .
[29] G.-E. Zhang, F.-L. He, H.-G. Jia, Y.-H. Lai, "Analysis of gateroad stability in relation to yield pillar size: a case study," Rock Mechanics and Rock Engineering, vol. 50 no. 5, pp. 1263-1278, DOI: 10.1007/s00603-016-1155-1, 2017.
[30] Z. Z. Zhang, F. Gao, "Research on nonlinear characteristics of rock energy evolution under uniaxial compression," China Journal of Rock Mechanics and Engineering, vol. 31, pp. 1198-1207, 2012.
[31] D. Song, E. Wang, Z. Liu, X. Liu, R. Shen, "Numerical simulation of rock-burst relief and prevention by water-jet cutting," International Journal of Rock Mechanics and Mining Sciences, vol. 70, pp. 318-331, DOI: 10.1016/j.ijrmms.2014.05.015, 2014.
[32] Q. B. Meng, M. W. Zhang, Han, "Effects of acoustic emission and energy evolution of rock specimens under the uniaxial cyclic loading and unloading compression," Rock Mechanics and Rock Engineering, vol. 49, pp. 3873-3886, 2016.
[33] C. Q. Zhu, L. Wang, X. G. Han, "State judgement model of the coal and rock medium and its engineering application," Advances in Civil Engineering, vol. 2020,DOI: 10.1155/2020/4670876, 2020.
[34] A. H. Wilson, F. Carr, "A new approach to the design of multientry developments of retreat longwall mining, " Proceedings of the 2nd International Conference on Ground Control in Mining, .
[35] A. A. Campoli, T. M. Barton, F. C. Dyke, M. Gauna, "Gob and gate road reaction to longwall mining in bump-prone strata," Bureau of Mines, vol. 9445, pp. 48-56, 1993.
[36] Z. Z. Zhang, M. Deng, Bai, "Stability control of gob-side entry retained under the gob with close distance coal seams," International Journal of Mining Science and Technology, vol. 11, pp. 234-249, 2020.
[37] J. N. Liu, M. C. He, Y. J. Wang, "Stability analysis and monitoring method for the key block structure of the basic roof of non-coal pillar mining with automatically formed gob-side entry," Advances in Civil Engineering, vol. 2019,DOI: 10.1155/2019/5347683, 2019.
[38] Z. Li, S. Yu, W. Zhu, "Dynamic loading induced by the instability of voussoir beam structure during mining below the slope," International Journal of Rock Mechanics and Mining Sciences, vol. 132,DOI: 10.1016/j.ijrmms.2020.104343, 2020.
[39] Z. Tian, Z. Zhang, M. Deng, S. Yan, J. Bai, "Gob-side entry retained with soft roof, floor, and seam in thin coal seams: a case study," Sustainability, vol. 12 no. 3, pp. 1197-1214, DOI: 10.3390/su12031197, 2020.
[40] S. Yan, J. Bai, X. Wang, L. Huo, "An innovative approach for gateroad layout in highly gassy longwall top coal caving," International Journal of Rock Mechanics and Mining Sciences, vol. 59, pp. 33-41, DOI: 10.1016/j.ijrmms.2012.11.007, 2013.
[41] H. Wu, X. K. Wang, Yu, "Analysis of influence law of burial depth on surrounding rock deformation of roadway," Advances in Civil Engineering, vol. 2020,DOI: 10.1155/2020/8870800, 2020.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2021 FuZhou Qi et al. This is an open access article distributed under the Creative Commons Attribution License (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License. https://creativecommons.org/licenses/by/4.0/
Abstract
The width design of protective pillars is an important factor affecting the stability of high-stress roadways. In this study, a novel numerical modeling approach was developed to investigate the relationship between protective pillar width and roadway stability. With the 20 m protective pillar width adopted in the field test, large deformation of roadways and serious damage to surrounding rocks occurred. According to the case study at the Wangzhuang coal mine in China, the stress changes and energy density distribution characteristics in protective pillars with various widths were analysed by numerical simulation. The modeling results indicate that, with a 20 m wide protective pillar, the peak vertical stress and energy density in the pillar are 18.5 MPa and 563.7 kJ/m3, respectively. The phenomena of stress concentration and energy accumulation were clearly observed in the simulation results, and the roadway is in a state of high stress. Under the condition of a 10 m wide protective pillar, the peak vertical stress and energy density are shifted from the pillar to roadway virgin coal region, with peak values of 9.5 MPa and 208.3 kJ/m3, respectively. The decrease in vertical stress and energy density improves the stability of the protective pillar, resulting in the roadway being in a state of low stress. Field monitoring suggested that the proposed 10 m protective pillar width can effectively control the large deformation of the surrounding rock and reduce coal bump risk. The novel numerical modeling approach and design principle of protective pillars presented in this paper can provide useful references for application in similar coal mines.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details


1 School of Civil & Architecture Engineering, Zhongyuan University of Technology, Zhengzhou 450007, Henan, China
2 State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, Jiangsu, China
3 Earthquake Administration of Henan, Zhengzhou 450016, Henan, China