# Study on the Retention of Large Mining Height and Small Coal Pillar under Thick and Hard Roof of Bayangaole Coal

Authors: Shankun Zhao; Jianping Zuo; Lei Liu; Kaijun Wu
Journal: Advances in Civil Engineering
 (2021)

---

## Abstract

The coal pillar stress distribution at the 311102 working face in the Bayangaole Mine is analyzed and revealed. In addition, borehole stressmeter, PASSAT monitoring system, and numerical modelling are fully utilized. Based on the patterns of acoustic wave velocity distribution, it is discovered that the impact created by mining activity can expand into the working face around 40 m, where the peak stress concentration is found about 15 m ahead. According to borehole stressmeter readings, mildly impacted, ordinarily impacted, and severely impacted zones are distinguished. The equilibrium theory and corresponding calculation indicated that the coal body in front of the working face has a plastic zone width of 4.96 m. The stress-displacement analysis based on numerical simulation showed that the relationship between peak vertical stress and pillar width is unimodal and bimodal. Specifically, both 5 and 10 m wide pillars showed a unimodal stress-width correlation and the peak vertical stresses are all located at the pillar center, whereas 15 m wide pillar has a bimodal stress-width relationship. In comparison, 10 m wide pillar holds the maximum in-situ stress. In consideration of site conditions and economic influences, 6 m wide coal strip coal pillar is designed at the working face 311102. As a result, stopping was successfully completed, and remarkable economic benefits were achieved.


---

## Body

## 1. Introduction

In China, coal is the major component of nonrenewable energy resources. Optimizing recovery rate and mining safety is of great significance in both China’s energy security system establishment and sustainable development strategy implementation. In fact, China has abundant reserves of thick coal seam, where thick coal seam mines make up 41% of the total coal mines. And, the reserves of thick coal mines, which are currently in production, accounted for 45% of the total coal reserves in China. Many mining methods were developed for excavation of extremely thick coal seams based on various geological conditions. However, it is still common that most coal mines adopted relatively large pillars for maintaining the preparation roadway, and it does waste a substantial amount of coal resources. Therefore, the research for nonpillar excavation becomes the key technology for the sake of safety and effectiveness.The traditional design width of the strip coal pillar is in-between 10 and 40 m with coal loss rate ranging from 10% to 30%, and this design approach mainly focuses on ventilation, transportation, and safety. When its roadway stability is impacted by stopping activities, the stress concentration always causes maintaining difficulties. Therefore, some mining experts proposed nonpillar (small coal pillar) goaf-side entrydriving method, which only leaves 3–7 m thick strip coal pillar on the gob side, to ensure the roadway stability and improve recovery rate after the overburden is stabilized. Goaf-side entrydriving is normally started after several months or if the working face has already advanced to a certain distance in the upper panel. The excavation advances along with the goaf side and only leaves a relatively thin coal pillar to isolate and protect the roadway. According to the mine pressure distribution, goaf-side entrydriving is normally longitudinally located in the stable bearing pressure zone, and transversally located in the coal wall plastic zone. Therefore, gob-side entrydriving can not only improve the recovery rate sufficiently but also reduce both maintaining cost and workload as the excavation is conducted in the plastic zone, where the confining pressure is relatively low.In 1911, the first empirical formula for coal pillar strength was presented by Bunting via testing the strength of anthracite blocks of different sizes and shapes [1]. Then, Gaddy [2] presented the relation between the strength of test blocks decreasing and the size increasing, and simultaneously the calculation formula for coal pillar strength with consideration of laboratory test values after testing the strength of coals with different sizes was also proposed. Iannacchione and Mark [3] evaluated several theories of coal pillar strength after testing the actual coal pillar stress. Based on the investigation and analysis of a large number of unstable coal pillars in South Africa, Salamon et al. [4, 5] summarized that most successful coal pillar designs adopted the inversion method. Chen et al. [6] studied the influence of confining pressure and confining pressure unloading rate on crack growth characteristics. At the same time, after analyzing 23 unstable coal pillars and 20 stable pillars in India, a strength formula suitable for coal pillars with different ratios of width to height was presented by Loui and Sheorey [7].In recent years, China’s coal mining technicians conducted a large number of research works on the controlling techniques of country-rock under the direct influence of stopping activities. Li et al. [8] collaborated six classical coal mining cases from the East to the West of China, and the key factors that induced deformative failure on thin coal pillars were analyzed and collaborated. By applying numerical analysis, Bai et al. [9] investigated how coal pillar width, coal seam mechanical characteristics, and rock bolt supporting strength affects both country-rock deformation and coal pillar strength during the gob-side entrydriving. Sun et al. [10] came up with a novel method for predicting movement and damage of inner burden caused by coal mining. Based on FLAC3D numerical analyses and engineering site experiments, Zheng et al. [11] obtained a series of stress distribution behaviors in coal pillars during the entire stopping process, and specifically investigated the influences, which are created by pillar width variation, on both stopping pillar and the evolution of coal rib stress; Zhang et al. [12] further researched and ascertained the width of coal pillars based on five aspects, which are the distribution of upper goaf rib-bearing stress, upper goaf stress field distribution, upper goaf displacement field distribution, the correlations between roadway country-rock and pillar width, as well as the equilibrium theory of small pillar width. Based on macro and micro failure experiments and numerical analysis, Zuo et al.[13] established the theoretical model of gradient failure of the surrounding rock mass of deep roadway, and presented the relationship between relative stress gradient and the average tangential stress level provided by the surrounding rock mass.Wei [14] measured the deformation and stress distribution of room and pillar-mined coal pillars. He et al. [15] established the global constitutive double yield and strain-softening model and studied the abutment pressure and plastic failure range of coal seams with different thickness and pillar width. Frith and Reed [16] explore the reasons why empirically derived coal pillar strength equations tend to be problem-specific and should be considered as providing no more than a pillar strength ‘‘index.” In order to reveal the mechanism of coal pillar impact from the perspective of energy evolution, Xue et al. [17] conducted an experiment on the bursting liability of coal samples during the third period of tunneling in no. 6 roadway of a certain mine. Li [18] studied the influence of coal pillar size on stress direction, plastic zone distribution, and roof stability of deep mining roadway. In addition, the study reveals the mechanical nature of the influence of coal pillar size on the shape of surrounding rock plastic zone in deep roadway. Based on the downhole pressure control theory, Shi et al. [19] expound the distribution characteristics of supporting stress in coal seam before work and its damage to coal seam. And, Chen et al. [20] studied the deformation of strip coal pillars under long-term stress by means of a borehole pressure gauge and displacement monitor. Chen et al. [21] proposed to use the natural strain described in Hooke’s law to establish the constitutive relation of coal and rock mass so as to accurately simulate the deformation of coal and rock mass.In general, a strip coal pillar always has certain strength and is hence capable of supporting overlying strata in the earlier stages. However, some coal pillars which were once stable may become unstable and collapse, leading to surface subsidence. American coal mines have suffered surface collapse due to coal pillar instability more than 100 years after construction. An abandoned mine in Scotland caused surface subsidence 118 years after it had been scrapped. The Wieliczka salt mine in Poland suddenly exhibited surface subsidence after 140 years of mining, resulting in the destruction of all surface buildings. And most unstable coal pillars in South Africa have suffered from surface subsidence after many years of mining. These facts indicate that although the deformation of pillars becomes stable in the short term, they may deteriorate after a long period of damage accumulation, leading to unstable failure in the coal pillar or seam roof. Therefore, the long-term bearing capacity of the strip coal pillar is a key factor in determining whether or not the safety of surface structures (buildings) in the mining area could be maintained for a long period of time.In this study, the strip mining of Shandong Daizhuang Coal Mine was used as the experimental site to measure the bearing capacity of a strip coal pillar at different depths of a pillar during or after mining. In this way, the complicated effects of various factors on the strength and bearing capacity of the coal pillar are skillfully avoided. In addition, the study is also attempted to determine the actual bearing capacity of the strip coal pillar at different pillar depths. The results are of practical significance for the optimal design of strip coal pillars for long-term stability.Bayangaole mine is in Wushen county, Ordos city, Inner Mongolia. Coal pillar supporting technique was adopted, and its ventilation roadway is the auxiliary haulage gate of the upper working face. As the roadway is doubly impacted by mining activities, palpable floor heaves, country-rock deformation, and even local rock burst were triggered during the stopping process. Zuo et al. [22] analyzed strength and deformation and failure characteristics of samples under different combination modes, and discussed the influence of different combination modes on bursting liability and the interaction between coal and rock. Regarding coal pillar width design, theoretical analysis, numerical simulation, and engineering analogy are fully utilized in order to reveal the coal pillar mechanical characteristics, failure mode, overburden movement behaviors, and abutment pressure distribution. Zuo et al. [23] studied the influence of temperature, stress, composition, and other factors on rock permeability, and analyzed the changing law of sandstone samples under the changing rate of confining pressure and pore pressure after high temperature treatment. Meanwhile, this research is significant to the recovery rate improvement; roadway supporting optimization; geological hazard prevention; especially for rock burst; economic benefit improvement; and safety assurance. Sun et al. [24] studied the dynamic fracturing process of coal seam hard rock coal pillar.


## 2. Bayangaole Mine Geological Conditions and Stress Distribution

Bayangaole mine is located 3 km to the northeast of Ulan Tolgoi, Wushen county, Ordos city, Inner Mongolia. This mine has a regular polygon planar configuration, in which the East-West length is 10.64 km, North-South width is 7.79 km, and the area is 65.2734 km2.Bayangaole mine is in the southeast of Dongsheng coal field, and most coal seams are deeply buried with burial depth greater than 660 m. Three major seams have been discovered in this mine, where 2-1 and 2-2 coal seams are partially extractable, and 3-1 coal seam is fully extractable. Based on the existing geological investigations, 3-1 coal seam only interacts with simple geological structures, and rock parting was barely found. The roof consists of silty/sandy mudstone with some minor siltstone and medium to small grained sandstone, whereas the floor is mainly constituted of sandy mudstone.


### 2.1. General Situation of Working Face 311102

Working face 311102 is the second face of the 11th panel zone in 3-1 coal seam at Bayangaole mine, with a trend length of 3578 m, a dip length of 260 m, an area of 930280 m2, and a burial depth of 610–629 m. Figure 1 illustrates the surroundings of this working face, where the goaf area caused by the 31110 working face is on the east side, the unexcavated working face 311103 is on the west side, whereas the water source protecting the pillar of coal and the 3-1coal auxiliary roadway are on the north and south side, respectively. The average coal seam thickness of this working face is 5.42 m and the coal seam has an average dip angle of 1.5 degrees.Figure 1
Schematic diagram of four-neighborhood relationships on working face of 311102.Figure2 shows the bore-log profile of coal seam 3-1, where the borehole directly passed through the coal panel. The immediate roof in this coal seam is sandy mudstone with a thickness variation of 0.98–5.97 m, dark grey in color, and contains phytolite pieces and sliding surfaces with shale-like fractures. The main roof is small grained sandstone with a thickness variation of 9.9–36.83 m, pile grey in color, interbedded with thin calcareous sandstone. Moreover, it has medium-high hardness. The main floor is medium-grained sandstone, 8.45 m thick, grey to white in color with medium harness. The physical-mechanical parameters of the rocks mentioned above are all listed in Table 1.Figure 2
3-1 coal roof and floor lithology column chart.Table 1
Table of physical and mechanical parameters of coal strata.
Coal bed nameThickness (m)Bulk modulus (Pa)Shear modulus (Pa)Angle of internal friction (°)Cohesion (Pa)Tensile strength (Pa)Density (kg/m3)Siltstone136.03e95.30e928.081.06e72.07e62516Medium sandstone222.26e101.91e1037.086.20e62.48e62304Sandy mudstone31.89e109.73e931.461.02e71.11e62503Coal51.20e98.60e8278.00e61.85e61400Sandy mudstone41.89e109.73e931.461.02e71.11e62503Medium sandstone82.73e92.31e932.548.00e61.17e62108Siltstone186.03e95.30e928.081.06e72.07e62516


### 2.2. Ground Stress Characteristics

In-situ stress is one of the most significant parameters that are required for the study on coal pillar and roadway stability. Overcorying stress relief method was adopted for in-situ stress testing, and a site testing approach, which utilizes a refined Hollow Inclusion Strain Gauge, was applied. The refined Hollow Inclusion Strain Gauge has developed a new collection board that provides the functions of simultaneous collection, collection without external power supply, and automatic drifting correction. In addition, temperature compensation circuit was also designed and developed. ADuC847 microprocessor was integrated into this circuit with 5 V battery power supply; voltage regulation chip was used instead of the traditional voltage regulation model and each channel has an individual analog switch. The application of simultaneous collection with temperature compensation technology can ensure the signal integrity. This circuit has 16 signal input channels, in which 1–12 are the strain signal channels, 13–14 are the temperature compensation channels, and 15–16 are the signal distortion compensation channels. The temperature compensation algorithm built in this refined Hollow Inclusion Strain Gauge can reach up to 15% correction for measured data, and it only has 8% maximum error band with theoretical data.Testing locations are determined based on following principles. A representative area of formation; In a complete or as complete rock mass as possible, it is generally required to stay away from faults and avoid rock fracture zones and fault development zones. Zuo et al. [25] used LS-DYNA numerical simulation software to analyze and compare the variation trend of the mises stress evolution and effective stress strain at the observation points around the hole under three layout modes in Bayangaole coal mine. People should stay away from or as far as possible from larger excavated bodies, such as large goaf, large adit chamber, etc. Liu et al. [26] carried out multi-scale hierarchical analysis of rock mass and predicted its mechanical and hydraulic characteristics. Yu et al. [27] conducted tests on mechanical properties of sandstone and coal bolts under conventional triaxial compression at different angles. Wu et al. conducted experimental and theoretical research on the mechanical properties of coal rock bolt (RCB) composite system [28]. The stress concentration areas such as bend, fork, bend and top of roadway and stope should be avoided, and the stress measuring point must be located in the stress area of the original rock—the area where the original stress state is not disturbed by the engineering. In order to study the variation of stress state with depth, measurements should be made at three or more levels as far as possible. In addition to the above principles, the actual conditions on-site should also be considered. According to the above principles, combined with detailed scene investigation, it is decided to arrange three geostress measuring points in Bayangaole mine. The relevant parameters of each measuring point are shown in Table 2, and the positions of each measuring point are shown in Figure 3.Table 2
Borehole parameters table at measuring point.
DotLocationElevation (m)Depth (m)Hole depth (m)13-1 high fructose 6 conduction point of coal return ventilation roadway65062082Guide point of S4 in the main channel 31110365162083Guide point J42 of the auxiliary channel on the 311102 working face6466368Figure 3
Schematic diagram of ground stress measuring point layout. (a) Point 1. (b) Point 2. (c) Point 3.
(a)(b)(c)Most of the stress relief curves obtained from the three measuring points in Bayangaole mine change regularly, indicating that the working state of most strain gauges is normal. Stress relief curves of the three measuring points are shown in Figure4.Figure 4
Stress relief curve at measuring point. (a) 1# measuring point. (b) 2# measuring point. (c) 3# measuring point.
(a)(b)(c)According to the tested strain readings, rock mechanical parameters, and the bore-hole geometric parameters, the principal in-situ stresses with their magnitude and direction are analyzed and calculated. The in-situ stress measurement results of Bayangaole mine are shown in Table3.Table 3
Calculation results of principal stress at each measuring point of Bayangaol mine.
PointDepth (m)Maximum principal stress (σmax)Intermediate principal stress (σmid)Minimum principal stress (σmin)Numerical (MPa)Direction (°)Dip angle (°)Numerical (MPa)Direction (°)Dip angle (°)Numerical (MPa)Direction (°)Dip angle (°)162022.95112.5817.7619.8419.1110.6815.98259.5269.08262029.45102.65−13.418.05−10.15−10.3415.53116.32−72.9363624.87269.70−2.8016.410.31−12.1115.71166.8977.56Table3 shows the testing results of in-situ stresses at three locations in Bayangaole mine, and the characteristics of in-situ stress field are discovered. The major principal stress is almost in a horizontal direction, and the ratio of maximum principal stress to its gravitational stress is 1.43–1.89 with an average value of 1.63. It is indicated that the in-situ stress filed in Bayangaole mine is dominated by its horizontal tectonic stress. Vertical stress is found to be equal to or slightly greater than the overburden weight. The average values of both major principal stresses and vertical stress are found to be almost directly proportional to the burial depth. As a result, the average values of σhmax = 25.76 MPa, σhmin = 18.1 MPa, and σv = 15.74 MPa are used for numerical simulation purposes.


## 2.1. General Situation of Working Face 311102

Working face 311102 is the second face of the 11th panel zone in 3-1 coal seam at Bayangaole mine, with a trend length of 3578 m, a dip length of 260 m, an area of 930280 m2, and a burial depth of 610–629 m. Figure 1 illustrates the surroundings of this working face, where the goaf area caused by the 31110 working face is on the east side, the unexcavated working face 311103 is on the west side, whereas the water source protecting the pillar of coal and the 3-1coal auxiliary roadway are on the north and south side, respectively. The average coal seam thickness of this working face is 5.42 m and the coal seam has an average dip angle of 1.5 degrees.Figure 1
Schematic diagram of four-neighborhood relationships on working face of 311102.Figure2 shows the bore-log profile of coal seam 3-1, where the borehole directly passed through the coal panel. The immediate roof in this coal seam is sandy mudstone with a thickness variation of 0.98–5.97 m, dark grey in color, and contains phytolite pieces and sliding surfaces with shale-like fractures. The main roof is small grained sandstone with a thickness variation of 9.9–36.83 m, pile grey in color, interbedded with thin calcareous sandstone. Moreover, it has medium-high hardness. The main floor is medium-grained sandstone, 8.45 m thick, grey to white in color with medium harness. The physical-mechanical parameters of the rocks mentioned above are all listed in Table 1.Figure 2
3-1 coal roof and floor lithology column chart.Table 1
Table of physical and mechanical parameters of coal strata.
Coal bed nameThickness (m)Bulk modulus (Pa)Shear modulus (Pa)Angle of internal friction (°)Cohesion (Pa)Tensile strength (Pa)Density (kg/m3)Siltstone136.03e95.30e928.081.06e72.07e62516Medium sandstone222.26e101.91e1037.086.20e62.48e62304Sandy mudstone31.89e109.73e931.461.02e71.11e62503Coal51.20e98.60e8278.00e61.85e61400Sandy mudstone41.89e109.73e931.461.02e71.11e62503Medium sandstone82.73e92.31e932.548.00e61.17e62108Siltstone186.03e95.30e928.081.06e72.07e62516


## 2.2. Ground Stress Characteristics

In-situ stress is one of the most significant parameters that are required for the study on coal pillar and roadway stability. Overcorying stress relief method was adopted for in-situ stress testing, and a site testing approach, which utilizes a refined Hollow Inclusion Strain Gauge, was applied. The refined Hollow Inclusion Strain Gauge has developed a new collection board that provides the functions of simultaneous collection, collection without external power supply, and automatic drifting correction. In addition, temperature compensation circuit was also designed and developed. ADuC847 microprocessor was integrated into this circuit with 5 V battery power supply; voltage regulation chip was used instead of the traditional voltage regulation model and each channel has an individual analog switch. The application of simultaneous collection with temperature compensation technology can ensure the signal integrity. This circuit has 16 signal input channels, in which 1–12 are the strain signal channels, 13–14 are the temperature compensation channels, and 15–16 are the signal distortion compensation channels. The temperature compensation algorithm built in this refined Hollow Inclusion Strain Gauge can reach up to 15% correction for measured data, and it only has 8% maximum error band with theoretical data.Testing locations are determined based on following principles. A representative area of formation; In a complete or as complete rock mass as possible, it is generally required to stay away from faults and avoid rock fracture zones and fault development zones. Zuo et al. [25] used LS-DYNA numerical simulation software to analyze and compare the variation trend of the mises stress evolution and effective stress strain at the observation points around the hole under three layout modes in Bayangaole coal mine. People should stay away from or as far as possible from larger excavated bodies, such as large goaf, large adit chamber, etc. Liu et al. [26] carried out multi-scale hierarchical analysis of rock mass and predicted its mechanical and hydraulic characteristics. Yu et al. [27] conducted tests on mechanical properties of sandstone and coal bolts under conventional triaxial compression at different angles. Wu et al. conducted experimental and theoretical research on the mechanical properties of coal rock bolt (RCB) composite system [28]. The stress concentration areas such as bend, fork, bend and top of roadway and stope should be avoided, and the stress measuring point must be located in the stress area of the original rock—the area where the original stress state is not disturbed by the engineering. In order to study the variation of stress state with depth, measurements should be made at three or more levels as far as possible. In addition to the above principles, the actual conditions on-site should also be considered. According to the above principles, combined with detailed scene investigation, it is decided to arrange three geostress measuring points in Bayangaole mine. The relevant parameters of each measuring point are shown in Table 2, and the positions of each measuring point are shown in Figure 3.Table 2
Borehole parameters table at measuring point.
DotLocationElevation (m)Depth (m)Hole depth (m)13-1 high fructose 6 conduction point of coal return ventilation roadway65062082Guide point of S4 in the main channel 31110365162083Guide point J42 of the auxiliary channel on the 311102 working face6466368Figure 3
Schematic diagram of ground stress measuring point layout. (a) Point 1. (b) Point 2. (c) Point 3.
(a)(b)(c)Most of the stress relief curves obtained from the three measuring points in Bayangaole mine change regularly, indicating that the working state of most strain gauges is normal. Stress relief curves of the three measuring points are shown in Figure4.Figure 4
Stress relief curve at measuring point. (a) 1# measuring point. (b) 2# measuring point. (c) 3# measuring point.
(a)(b)(c)According to the tested strain readings, rock mechanical parameters, and the bore-hole geometric parameters, the principal in-situ stresses with their magnitude and direction are analyzed and calculated. The in-situ stress measurement results of Bayangaole mine are shown in Table3.Table 3
Calculation results of principal stress at each measuring point of Bayangaol mine.
PointDepth (m)Maximum principal stress (σmax)Intermediate principal stress (σmid)Minimum principal stress (σmin)Numerical (MPa)Direction (°)Dip angle (°)Numerical (MPa)Direction (°)Dip angle (°)Numerical (MPa)Direction (°)Dip angle (°)162022.95112.5817.7619.8419.1110.6815.98259.5269.08262029.45102.65−13.418.05−10.15−10.3415.53116.32−72.9363624.87269.70−2.8016.410.31−12.1115.71166.8977.56Table3 shows the testing results of in-situ stresses at three locations in Bayangaole mine, and the characteristics of in-situ stress field are discovered. The major principal stress is almost in a horizontal direction, and the ratio of maximum principal stress to its gravitational stress is 1.43–1.89 with an average value of 1.63. It is indicated that the in-situ stress filed in Bayangaole mine is dominated by its horizontal tectonic stress. Vertical stress is found to be equal to or slightly greater than the overburden weight. The average values of both major principal stresses and vertical stress are found to be almost directly proportional to the burial depth. As a result, the average values of σhmax = 25.76 MPa, σhmin = 18.1 MPa, and σv = 15.74 MPa are used for numerical simulation purposes.


## 3. Monitoring the Stress Distribution of Mining Coal Pillar in Working Face 311102

During the mining, a series of ore pressure shows up in the working face 311102. In particular, serious floor heave and slab appear in the advance range of 200 m in the return air roadway of 02 working face. After the working face advanced over 200 m, serious floor heaving in the auxiliary haulage gate was also triggered by stopping. Hence, the mining stress in 02 working face was closely monitored for further pillar design purposes.


### 3.1. Wave Velocity Field of Coal Body in Pillar

Due to the stress field complexity and variation in coal pillar, the distribution of sonic wave velocity was tested as the stress detector can only acquire the point value. PASAT-M portable microseismic detecting technique as a geophysical method is designed to detect the structure of underground medium. Both sonic wave velocity distribution and corresponding elastic parameters were also obtained based on the analysis of seismic wave characteristics, which were acquired by PASAT-M. For the prevention of severe mining disasters, PASAT-M is used to monitor the failure motion in country-rock during excavation, and the PASAT-M testing system is sketched out in Figure5(a).Figure 5
PASAT-MPortable microseismic detection system and coal pillar stress field. (a) Seismic CT technical schematic diagram of the coal mining face. (b) Seismic CT technical schematic diagram of the coal mining face. (c) 311102 working face pillar of the coal PASAT—M monitoring programme. (d) 311102 coal pillar wave velocity field between main transport and auxiliary transport on the working face.
(a)(b)(c)(d)In the case of seismic travel time imaging, it is assumed that seismic waves propagate in the medium inside the detection area in the form of rays, as shown in Figure5(b). For image reconstruction of velocity field, common methods include backprojection (BPT), algebraic reconstruction (ART), joint iterative reconstruction (SIRT), conjugate gradient least square method, etc. SIRT was first proposed by Gilbert. Its formula is(1)f0the initial value is arbitrary,fjk+1=fjk+∑i=1Mdijti−di,fk/∑j=1Ndijμ+∑i=1Mdij.In the formulai = 1, 2, …, M is the ray symbol, j = 1, 2, …, N is the unit number, k = 0, 1, …, are the number of iterations.Whenfk+1−fk∞<ε, take f^=fk+1. Where is the allowable error of a given iteration?Compared with the traditional BPT and ART methods, the SIRT algorithm has the advantages of good convergence, high resolution of reconstructed image, and low dependence on the accuracy of initial value selection. Most of the results are inversion results obtained by the SIRT method for geological interpretation of seismic wave CT images.The detection range is 130 m, from 5 to 135 m in front of the 311102 working face: the firing guns are arranged in the auxiliary transportation channel of the 311102 working face. There are a total of 27 guns, with a spacing of 5 m, blasting hole depth of 2 m, and charge of 100 g per hole. The receiving ends are arranged in the main canal groove, with a total of 11 (red short lines in Figure5(c)), and the spacing is 13 m. It is necessary to connect the signal line at 600 m, and reserve a round-shift (i.e., 24 hours) stopping distance with the working face for the first blasting hole in construction, usually at 20 m, to avoid the stopping of the hole. The specific layout of the detection system diagram is shown in Figure 5(c). On December 27, 2015, the stress field of coal pillars between the main transport and auxiliary transport channel of the 311102 working face was detected. The field design fired guns 27 times; 25 effective guns were received, of which the 4th and 13th guns were invalid. As we can see from Figure 5(c), the distribution of velocity lines at the excitation end and the receiving end. The denser the intersecting lines in the figure, the more accurate the effect of the velocity field. In the figure, because the linear density at both ends of detection is small, it is easy to cause boundary effect, so the analysis is omitted. Figure 5(b) shows the detection results of PASAT-M portable microseismic mode for the main and auxiliary coal pillars in the 02 working face. The horizontal axis is the main strike distance of the coal pillars, and the vertical axis is the width of the coal pillars.As shown in Figure5, the peak pillar stress is close to the working face (Red and Yellow), and pillar stress in 15–45 m zone is generally smaller than that of 15–45 m zone but greater than the other zones. The relatively higher pillar stress near the haulage gate indicates that the shifted symmetrical pillar stress distribution was clearly affected by mining activities on the 02 working face. In the 15–50 m zone, peak pillar stress is still at the pillar center, but the higher side stress is shifted to the main gate. Pillar stress in the 45–90 m zone presents a saddle shape distribution, where side-stress is relatively higher and inner-stress is relatively lower with an 8 m diameter zone. This phenomenon revealed the conservativity of the 30 m thick coal pillar design in Bayangaole mine and showed optimization opportunities. In a stress zone greater than 90 m, much lowered pillar stress is presented in blue color, which proved that the advance stress range is smaller than 90 m, and it matches with the borehole stressmeter readings.


### 3.2. Observation of 3# Coal Pillar Stress at the Location of the Roadway

For the estimation of stress distribution behaviors, borehole stressmeters were applied to discover the local point stress, and form correlations with the sonic wave velocities that occur across entire coal pillar. In total, 9 borehole stressmeters were placed in the coal pillar between the 02 haulage gate and 03 ventilation gate with a drilling depth of 10 m and a drilling height of 1.5–1.8 m. In consideration of pillar stress variation on the 02 working face, boreholes were designed with inconsistent distances along with the pillar wall between the 02 haulage gate and 03 ventilation gate as shown in Figure6.Figure 6
3# schematic diagram of borehole stressmeter layout in the roadway.Stressmeters were installed simultaneously; the stress readings were recorded for each 30 min when the stress variation is greater than 0.5 MPa, and all the stress data were automatically saved from January to June in 2016. As the amount of recorded stress data is massive, minor stress variation and outliers, which were created by the relatively far distance from the working face and the instability of monitoring system, have been neglected and the mean stress was taken every day.It can be seen from Figure7 that the pillar stress increases first and drops suddenly overall along with the working face advancement. The tested pillar stress at monitoring point No.2 was increased and then dropped; this phenomenon matches with the previous testing result, and it proved that the roof was broken at 6 m from the working face which is in the small tectonic stress dropping zone.Figure 7
The variation of stress with distance at different measuring points (ΔP-L) graph. A: 2 m, B: 4.44 m, C: 6.44 m, D: 8.44 m, E: 10.94 m, and F: 13.44 m.According to the characteristics of pillar stress variation as shown in Figure7, coal body along with the working face direction that are affected by excavation activities can be distinguished into 3 zones as shown below.(1)
Mildly Impacted zone is in the coal body which is more than 24 m ahead of the working face. In this zone, bearing pressure and country-rock pressurization are both minor, hence the coal body is barely affected by the excavation activities.(2)
Ordinarily Impacted zone is in the coal body which is between 8 and 24 m ahead of the working face. In this zone, the increment of bearing pressure and the corresponding country-rock deformation are more intensive compared with the mildly impacted zone.(3)
Severely Impacted zone is in the coal body which is less than 8 m ahead of the working face. In this zone, the excavation caused a significant bearing pressure increment, where the maximum bearing pressure is located around 5 m ahead of the working face, whereas the bearing pressure at the working face is much smaller.It can be seen from Figure8 that the stress is high at the measuring point 1# (2 m), low at the measuring points 2# and 3# (4.44 m and 6.44#), rising at the measuring points 4# and 5#, peak at the measuring point 6# (13.44), and then decreases. According to the above results, the coal pillar stress is maximum at 13.44 m, which is consistent with the simulation result that the coal pillar stress is superimposed at 15 m. At 4.44 m and 6.44 m, the stress is low, indicating a small structure.Figure 8
Internal stress curve of coal pillar.Therefore, the initial excavation has a relatively low impact on coal pillar, and a bimodal shape relationship was discovered based on the stress distribution curve.When the working face is within 13.94 m from the ventilation gate, coal pillar was greatly impacted, and the pillar stress distribution showed a similar bimodal correlation with both peak stresses at 2 and 8 m from the coal rib. The stress crest appeared in-between 4 and 6 m.By referring to the previous loose circle test results, the loss circle diameter of the doubly impacted pillar is determined as 3.6 m. In consideration of water proofing and ventilation efficiency, this long wall mining excavation is design with a 4–6 m coal pillar thickness, tunnel centerline at 8 m, and a roadway width within 4 m.


## 3.1. Wave Velocity Field of Coal Body in Pillar

Due to the stress field complexity and variation in coal pillar, the distribution of sonic wave velocity was tested as the stress detector can only acquire the point value. PASAT-M portable microseismic detecting technique as a geophysical method is designed to detect the structure of underground medium. Both sonic wave velocity distribution and corresponding elastic parameters were also obtained based on the analysis of seismic wave characteristics, which were acquired by PASAT-M. For the prevention of severe mining disasters, PASAT-M is used to monitor the failure motion in country-rock during excavation, and the PASAT-M testing system is sketched out in Figure5(a).Figure 5
PASAT-MPortable microseismic detection system and coal pillar stress field. (a) Seismic CT technical schematic diagram of the coal mining face. (b) Seismic CT technical schematic diagram of the coal mining face. (c) 311102 working face pillar of the coal PASAT—M monitoring programme. (d) 311102 coal pillar wave velocity field between main transport and auxiliary transport on the working face.
(a)(b)(c)(d)In the case of seismic travel time imaging, it is assumed that seismic waves propagate in the medium inside the detection area in the form of rays, as shown in Figure5(b). For image reconstruction of velocity field, common methods include backprojection (BPT), algebraic reconstruction (ART), joint iterative reconstruction (SIRT), conjugate gradient least square method, etc. SIRT was first proposed by Gilbert. Its formula is(1)f0the initial value is arbitrary,fjk+1=fjk+∑i=1Mdijti−di,fk/∑j=1Ndijμ+∑i=1Mdij.In the formulai = 1, 2, …, M is the ray symbol, j = 1, 2, …, N is the unit number, k = 0, 1, …, are the number of iterations.Whenfk+1−fk∞<ε, take f^=fk+1. Where is the allowable error of a given iteration?Compared with the traditional BPT and ART methods, the SIRT algorithm has the advantages of good convergence, high resolution of reconstructed image, and low dependence on the accuracy of initial value selection. Most of the results are inversion results obtained by the SIRT method for geological interpretation of seismic wave CT images.The detection range is 130 m, from 5 to 135 m in front of the 311102 working face: the firing guns are arranged in the auxiliary transportation channel of the 311102 working face. There are a total of 27 guns, with a spacing of 5 m, blasting hole depth of 2 m, and charge of 100 g per hole. The receiving ends are arranged in the main canal groove, with a total of 11 (red short lines in Figure5(c)), and the spacing is 13 m. It is necessary to connect the signal line at 600 m, and reserve a round-shift (i.e., 24 hours) stopping distance with the working face for the first blasting hole in construction, usually at 20 m, to avoid the stopping of the hole. The specific layout of the detection system diagram is shown in Figure 5(c). On December 27, 2015, the stress field of coal pillars between the main transport and auxiliary transport channel of the 311102 working face was detected. The field design fired guns 27 times; 25 effective guns were received, of which the 4th and 13th guns were invalid. As we can see from Figure 5(c), the distribution of velocity lines at the excitation end and the receiving end. The denser the intersecting lines in the figure, the more accurate the effect of the velocity field. In the figure, because the linear density at both ends of detection is small, it is easy to cause boundary effect, so the analysis is omitted. Figure 5(b) shows the detection results of PASAT-M portable microseismic mode for the main and auxiliary coal pillars in the 02 working face. The horizontal axis is the main strike distance of the coal pillars, and the vertical axis is the width of the coal pillars.As shown in Figure5, the peak pillar stress is close to the working face (Red and Yellow), and pillar stress in 15–45 m zone is generally smaller than that of 15–45 m zone but greater than the other zones. The relatively higher pillar stress near the haulage gate indicates that the shifted symmetrical pillar stress distribution was clearly affected by mining activities on the 02 working face. In the 15–50 m zone, peak pillar stress is still at the pillar center, but the higher side stress is shifted to the main gate. Pillar stress in the 45–90 m zone presents a saddle shape distribution, where side-stress is relatively higher and inner-stress is relatively lower with an 8 m diameter zone. This phenomenon revealed the conservativity of the 30 m thick coal pillar design in Bayangaole mine and showed optimization opportunities. In a stress zone greater than 90 m, much lowered pillar stress is presented in blue color, which proved that the advance stress range is smaller than 90 m, and it matches with the borehole stressmeter readings.


## 3.2. Observation of 3# Coal Pillar Stress at the Location of the Roadway

For the estimation of stress distribution behaviors, borehole stressmeters were applied to discover the local point stress, and form correlations with the sonic wave velocities that occur across entire coal pillar. In total, 9 borehole stressmeters were placed in the coal pillar between the 02 haulage gate and 03 ventilation gate with a drilling depth of 10 m and a drilling height of 1.5–1.8 m. In consideration of pillar stress variation on the 02 working face, boreholes were designed with inconsistent distances along with the pillar wall between the 02 haulage gate and 03 ventilation gate as shown in Figure6.Figure 6
3# schematic diagram of borehole stressmeter layout in the roadway.Stressmeters were installed simultaneously; the stress readings were recorded for each 30 min when the stress variation is greater than 0.5 MPa, and all the stress data were automatically saved from January to June in 2016. As the amount of recorded stress data is massive, minor stress variation and outliers, which were created by the relatively far distance from the working face and the instability of monitoring system, have been neglected and the mean stress was taken every day.It can be seen from Figure7 that the pillar stress increases first and drops suddenly overall along with the working face advancement. The tested pillar stress at monitoring point No.2 was increased and then dropped; this phenomenon matches with the previous testing result, and it proved that the roof was broken at 6 m from the working face which is in the small tectonic stress dropping zone.Figure 7
The variation of stress with distance at different measuring points (ΔP-L) graph. A: 2 m, B: 4.44 m, C: 6.44 m, D: 8.44 m, E: 10.94 m, and F: 13.44 m.According to the characteristics of pillar stress variation as shown in Figure7, coal body along with the working face direction that are affected by excavation activities can be distinguished into 3 zones as shown below.(1)
Mildly Impacted zone is in the coal body which is more than 24 m ahead of the working face. In this zone, bearing pressure and country-rock pressurization are both minor, hence the coal body is barely affected by the excavation activities.(2)
Ordinarily Impacted zone is in the coal body which is between 8 and 24 m ahead of the working face. In this zone, the increment of bearing pressure and the corresponding country-rock deformation are more intensive compared with the mildly impacted zone.(3)
Severely Impacted zone is in the coal body which is less than 8 m ahead of the working face. In this zone, the excavation caused a significant bearing pressure increment, where the maximum bearing pressure is located around 5 m ahead of the working face, whereas the bearing pressure at the working face is much smaller.It can be seen from Figure8 that the stress is high at the measuring point 1# (2 m), low at the measuring points 2# and 3# (4.44 m and 6.44#), rising at the measuring points 4# and 5#, peak at the measuring point 6# (13.44), and then decreases. According to the above results, the coal pillar stress is maximum at 13.44 m, which is consistent with the simulation result that the coal pillar stress is superimposed at 15 m. At 4.44 m and 6.44 m, the stress is low, indicating a small structure.Figure 8
Internal stress curve of coal pillar.Therefore, the initial excavation has a relatively low impact on coal pillar, and a bimodal shape relationship was discovered based on the stress distribution curve.When the working face is within 13.94 m from the ventilation gate, coal pillar was greatly impacted, and the pillar stress distribution showed a similar bimodal correlation with both peak stresses at 2 and 8 m from the coal rib. The stress crest appeared in-between 4 and 6 m.By referring to the previous loose circle test results, the loss circle diameter of the doubly impacted pillar is determined as 3.6 m. In consideration of water proofing and ventilation efficiency, this long wall mining excavation is design with a 4–6 m coal pillar thickness, tunnel centerline at 8 m, and a roadway width within 4 m.


## 4. Research on Reasonable Theoretical Width of Coal Pillar Retention

Strip coal pillar not only isolates stopping roadway from goaf but also maintains roadway stability. As the coal pillar is the main body that supports the roadway, its deformation and failure behaviors have a great influence on the controlling of roof strata. Based on the research, the extent of the influence that coal pillar geometry has created on the stopping roadway failure behaviors is one of the main factor relating to the stability of stopping road way country-rock. Therefore, reasonable coal pillar geometry can not only reduce both roadway deformation and maintaining workload but also improve the recovery rate compared with the overly designed pillar geometry.The traditional section pillar design leaves the coal pillar to support the roof between upper haulage drift and lower ventilation drift; hence, it allows the drifts in lower panels to avoid the peak abutment pressure zone. It is true that two-way advancing method in section drifts is beneficial for technical management, ventilation, transportation, drainage, and mining safety. However, the loss of coal pillar is as high as 10%∼30%, and the return air roadway is affected by secondary mining. Therefore, the maintenance of roadway is difficult, and the support cost is high. Not only this, the abandonment pressure that spreads to the floor through coal pillar can not only affect both surrounding excavation and floor stability but also produce risks on inducing rock burst.The most common failure modes, when the coal pillar stress and deformation reach its ultimate bearing capacity of deformation limit, are shear failure, splitting failure, rheodestruction failure, sliding failure along with weak plane and bedding, and heaving failure when coal pillar is punched into the floor.Coal pillar width is the major factor that affects its stability and roadway maintenance. The width of coal pillar determines the horizontal distance between roadway and mining space and affects the influence degree of abutment pressure caused by mining on roadway and the load of coal pillar. Therefore, coal pillar’s ultimate bearing capacity not only depends on its boundary conditions and mechanical properties but also relies on its geometry.Goaf in the upper panel is on one side of the strip coal pillar, whereas the preparation roadway in the lower panel is on another side. Plastic deformation zones are formed by both upper goaf and preparation roadway on both pillar sides with plastic zone widths ofx0 and x2, respectively, as shown in Figure 9. Therefore, the basic conditions for the stability of the coal pillar for roadway protection are as follows: after plastic deformation occurs on both sides of the coal pillar, there is an elastic core of a certain width in the center of the coal pillar, and the width of the elastic core should not be less than 2 times of the height of the coal pillar. On the contrary, when the elastic core width is smaller than 2 times the pillar height such as excavating a very narrow roadway in coal pillar, the redistributed pillar stress will reduce the effective bearing area and finally result in a dramatic drop in the pillar bearing capacity.Figure 9
Elastic-plastic deformation zone and stress distribution of coal pillar.There are two ways to avoid stress concentration based on the fact that peak rib bearing stress is located inside of the coal pillar as shown in Figure10. One approach is that the coal pillar must have a sufficient width; hence, the stress concentration cannot be released easily since it is confined by relatively thick coal body. Another approach is that the coal pillar must be small enough; hence, the peak stress will not fall into the coal pillar.Figure 10
Schematic diagram of mining stress.The alignment of Bayangaole mine 311102 large height working face is designed along with the 3–1 coal seam floor. After the face excavation started, the redistributed country-rock in-situ stress first caused rib coal failure, and then this failure expanded and shifted to the boundary of the elastic zone. Rock parting caused a lower cohesionC0 and internal friction angle ϕ0 when compared with that of pure coal mass; the controlling object is focused on the upper coal body over the large roadway height.A substantial amount of research outputs proved that the redistribution of bearing stress in country-rock is caused by the movement in overlaying strata and the behaviors of bearing stress distribution is mainly depending on the mechanical properties of coal mass and the overlaying strata movement. It is somehow similar to,(2)σy=C0tanϕ0e2tanϕ0/Mλx−C0tanϕ0,R0=Mλ2tanϕ0lnKγH+C0/tanϕ0C0/tanϕ0.In the formula:λ—Coefficient of lateral pressure. The formula λ=μ1−μ can get λ = 0.75, K—Coefficient of stress concentration due to stopping, take 3; γ—The average bulk weight of the overlying strata, take 2.5 t/m3; H—Mining depth, take 610 m; M—One-off mining thickness of coal seam, take 5.4 m; C0—The cohesive force of the interface such as the bedding plane of the coal and rock mass itself, from the laboratory test results, take 8 MPa; ϕ0—The internal friction angle of the interface such as the bedding plane of the coal and rock itself, from the laboratory test results, take 27°. Plug the data into the formula, calculated R0 = 4.96 m.Therefore, the theatrical plastic zone width of the section pillar (ultimate stress equilibrium width) in Bayangaole mine is determined as 4.96 m. If the small pillar mining approach is adopted, the peak stress zone can be effectively avoided, hence its stability is ensured. This theatrical estimation verifies the stress monitoring results from Section3. For further verifications, numerical simulation analysis was adopted to investigate the pillar stability.


## 5. Numerical Simulation of Coal Pillars with Different Widths

### 5.1. Simulation Scheme

FLAC3D model was established based on the geological features collaborated from 01, 02, 03 working faces, the 11th panel zone, Bayangaole mine. And, the small pillar excavation was simulated at 03 working face with trial pillar widths of 5, 8, 10, 12, and 15 m. By simulating the stopping activities of 02 and 03 working faces, the reasonable strip coal pillar width was determined based on pillar stress distribution, overlaying strata deformation, and failure behaviors.The FLAC3D model contains 19314 elements, 208050 nodes, and it has a geometry of 645–655 m (length)∗ 360 m (width) ∗ 74 m (height) with a 50 m extended outer boundary. The trial coal pillar widths are 3, 5, 8, 10, 12, 15, 20, 25, and 30 m. As shown in Figure 11, by simulating the stopping activities of 02 and 03 working faces, the reasonable strip coal pillar width was determined based on pillar stress distribution, overlaying strata deformation, and its failure behavior.Figure 11
FLAC3D model.Vertical stressV at 620 m = 15.74 MPa. Considering that the height of the whole model is 70 m, the overburden of the coal seam is 38 m, and the thickness of the coal seam is 5 m, the stress test point is generally in the middle of the coal seam. Therefore, a load of 14.7 MPa is applied vertically to the upper part of the model, which represents the weight of the upper part of the rock mass. The boundary conditions of the model include displacement boundary and stress boundary, where the boundary conditions of the model are set as follows: the left and right boundaries only constrain the displacement in the X direction, and the front and rear boundaries only constrain the displacement in the Y direction, that is, the single constraint boundary. The lower boundary is a fully constrained boundary, and the upper boundary is unconstrained. The Mohr-Coulomb yield criterion is used in the failure criterion of coal and rock mass.Fissure structures, including joints, have great influences on rock mass physical-mechanical parameters, and it consequently affects the results of numerical simulation. Therefore, rock mass properties were optimized based on the intact rock parameters obtained from the laboratory test, and the corresponding rock mass parameters are all listed in Table4.Table 4
Table of physical and mechanical parameters of coal strata.
Coal bed nameThick (m)Bulk modulus (GPa)Shear modulus (GPa)Internal friction (°)Cohesion (MPa)Tensile strength (MPa)Apparent density (kg/m3)Siltstone136.035.3028.0810.62.072516Medium sandstone2222.619.137.086.202.482304Sandy mudstone318.99.7331.4610.21.112503Coal51.200.86278.001.851400Sandy mudstone418.99.7331.4610.21.112503Medium sandstone82.732.3132.548.001.172108Siltstone186.035.3028.0810.62.072516


### 5.2. The Stress Variation Rule of Coal Pillar with Different Widths during Stopping of Working Face

As the strip coal pillar is located between 02 and 03 working face, influences created by mining actives are doubly accumulated. The entire simulation is started with 02 working face excavation, and the various trail coal pillar widths were set after the stress calculation is in the equilibrium state. And then, the stopping of the 03 working face was initiated, hence the coal pillar stress variation behaviors can be analyzed when stopping is simulated at both working faces. Figure12 shows the stress distribution variations in coal body during the stopping in both 02 and 03 working faces.Figure 12
Vertical stress distribution diagram of coal pillars with different widths of 02# and 03# working faces during mining.
(a)(b)As shown in Figure12, the coal pillar width variation does not affect the stress distribution much during the 02 working face excavation. When coal pillar width increased from 3 to 30 m, the corresponding peak vertical stresses are always in-between 40 and 50 MPa, but the pillar stress distribution changed as the peak stress zone shifted to the 02 goaf zone.In comparison, the coal pillar width variation has a relatively big influence on stress distribution at the 03 working face. The 3 m wide coal pillar was crushed, and its maximum stress concentration shifted to the roadway. Whereas, when coal pillar widths are 5, 8, 12, 15, 20, 25, and 30 m, the corresponding maximum stress and stress concentration coefficients are 71 MPa, 4.4; 82 MPa, 5.1; 145 MPa, 9; 116 MPa, 7.2; 89 MPa, 5.5; 73.6 MPa, 4.6; 69.3 MPa, 4.3; and 66.8 MPa, 4.1, respectively. An increasing trend is observed in the stress concentration factor when the pillar width increases from 5 to 10 m. When coal pillar width increases from 10 to 12 m, its concentration coefficient remains over 7, and it poses extreme danger to roadway stability; On the contrary, the stress concentration coefficient drops when pillar width increases from 10 to 30 m. Especially for coal pillar greater than 20 m, its corresponding stress concentration remains below 5, which is beneficial to the pillar stability.It is found that when coal pillar thickness is 15 m, the corresponding pillar stress is smaller than that of 10 and 12 m thick pillars, but still at a relatively high level. Therefore, the numerical simulation proved that pillar width between 5 and 8 m benefits its stability the most due to the relatively small corresponding peak vertical stress.


### 5.3. Stress Characteristics of Coal Pillar with Different Widths

The pillar stress characteristics were analyzed based on a representative scenario when the 03 working face is advanced by 130 m, as well as the strip pillar is influenced by both 02 goaf and 03 working face excavations. The vertical stress vs pillar width curves are shown in Figure13. It can be seen from Figure 13 that all the pillar stress distributions present a unimodal correlation and all peak stresses are shifted to the goaf side with a shifting range of 0.5–3 m. Meanwhile, vertical stress is always the biggest and the difference in-between both horizontal stresses is barely noticed. However, the horizontal stress, which is parallel to the roadway alignment, is still slightly greater than that of perpendicular to the roadway alignment. When coal pillar width is in-between 3 and 5 m, it was crushed and depressurized. When pillar width increases from 5 to 15 m, both vertical stress and alignment-parallel horizontal stress are concentrated. When pillar width further increases from 20 to 25 m, a bimodal correlation appeared and the peak stress close to the 02 working face is slightly higher than that of the 03 working face. It is also found that the stress at the pillar center is higher than the original rock stress, and the alignment-perpendicular horizontal stresses throughout the entire strip coal pillar are all smaller than the original rock stress.Figure 13
Stress distribution in coal pillar with different widths. (a) Vertical stress in pillar. (b) Horizontal stress in pillar (axial direction of parallel roadway). (c) Horizontal stress in pillar (vertical roadway axial direction).
(a)(b)(c)


### 5.4. Characteristics of Roadway Stress Field with Different Pillar Widths

The internal stress of the coal pillar is monitored by installing borehole stressmeter, and then the maximum stress in the coal pillar and solid coal of roadway side is extracted. The relationship between the characteristics of the roadway peak vertical stress and pillar width when the 03 working face was advanced by 130 m is shown in Figure14. It can be seen from Figure 14 that the pillar width variation has certain influence on the vertical stresses in both coal pillar and coal rib. Figure 13 reveals that the curves of both peak vertical stress vs coal pillar width present a bimodal trend. With the increase of coal pillar width, the vertical pillar stress increases sharply and then drops mildly, whereas the lower coal rib stress slightly increases and then drops. Figure 14 also indicates that coal pillar bearing capacity is positively correlated to the pillar width, and the lower coal rib stress has a pillar-side shifting trend with the increases of pillar width. When coal pillar width reaches 8 m, the upper rib coal pillar stress becomes much closer to that of lower coal rib, and until then all the lower coal rib stresses are slightly greater than pillar stress. With a further increment in coal pillar (width ≥ 10 m), the critical pillar stress remains much greater than that of lower coal rib stress, and it becomes the major load bearing object.Figure 14
Stress distribution of roadway under different pillar widths.Overall, peak pillar stress and peak lower-rib stress appear together when the coal pillar is 10 m wide, and the difference in-between both upper and lower-rib peak stresses is minor due to the side supporting effect created by 02 goaf.


### 5.5. Characteristics of Roadway Displacement Field with Different Pillar Widths

The analyses of roadway country-rock deformation at 200 m ahead of the open-off cut focus on roof-floor convergence and rib-rib convenience. As shown in Figure15, the accumulation of roof-floor convergence is inversely correlated to the working face distance, and it rapidly builds up when the working face approaches from 50 m away. The relationship between coal pillar width and relatively roof-floor displacement also presents a bimodal correlation. When coal pillar width increases from 5 to 10 m, the relatively roof-floor displacement also increases and reaches its peak value. And then, it gradually drops when coal pillar width increases from 10 to 30 m.Figure 15
Relative displacement curve of roadway roof and floor at different pillar widths. (a) Vertical displacement curve. (b) Horizontal displacement curve of two sides.
(a)(b)Regarding the rib-rib deformation, the results also reveal that its relative displacement is negatively correlated to the working face distance, and this trend intensifies further when the working face is within 40 m.Finally, based on the analyses of stress-displacement distributions in both coal pillar and coal rib, the strip pillar width is therefore recommended as 5–8 m.


## 5.1. Simulation Scheme

FLAC3D model was established based on the geological features collaborated from 01, 02, 03 working faces, the 11th panel zone, Bayangaole mine. And, the small pillar excavation was simulated at 03 working face with trial pillar widths of 5, 8, 10, 12, and 15 m. By simulating the stopping activities of 02 and 03 working faces, the reasonable strip coal pillar width was determined based on pillar stress distribution, overlaying strata deformation, and failure behaviors.The FLAC3D model contains 19314 elements, 208050 nodes, and it has a geometry of 645–655 m (length)∗ 360 m (width) ∗ 74 m (height) with a 50 m extended outer boundary. The trial coal pillar widths are 3, 5, 8, 10, 12, 15, 20, 25, and 30 m. As shown in Figure 11, by simulating the stopping activities of 02 and 03 working faces, the reasonable strip coal pillar width was determined based on pillar stress distribution, overlaying strata deformation, and its failure behavior.Figure 11
FLAC3D model.Vertical stressV at 620 m = 15.74 MPa. Considering that the height of the whole model is 70 m, the overburden of the coal seam is 38 m, and the thickness of the coal seam is 5 m, the stress test point is generally in the middle of the coal seam. Therefore, a load of 14.7 MPa is applied vertically to the upper part of the model, which represents the weight of the upper part of the rock mass. The boundary conditions of the model include displacement boundary and stress boundary, where the boundary conditions of the model are set as follows: the left and right boundaries only constrain the displacement in the X direction, and the front and rear boundaries only constrain the displacement in the Y direction, that is, the single constraint boundary. The lower boundary is a fully constrained boundary, and the upper boundary is unconstrained. The Mohr-Coulomb yield criterion is used in the failure criterion of coal and rock mass.Fissure structures, including joints, have great influences on rock mass physical-mechanical parameters, and it consequently affects the results of numerical simulation. Therefore, rock mass properties were optimized based on the intact rock parameters obtained from the laboratory test, and the corresponding rock mass parameters are all listed in Table4.Table 4
Table of physical and mechanical parameters of coal strata.
Coal bed nameThick (m)Bulk modulus (GPa)Shear modulus (GPa)Internal friction (°)Cohesion (MPa)Tensile strength (MPa)Apparent density (kg/m3)Siltstone136.035.3028.0810.62.072516Medium sandstone2222.619.137.086.202.482304Sandy mudstone318.99.7331.4610.21.112503Coal51.200.86278.001.851400Sandy mudstone418.99.7331.4610.21.112503Medium sandstone82.732.3132.548.001.172108Siltstone186.035.3028.0810.62.072516


## 5.2. The Stress Variation Rule of Coal Pillar with Different Widths during Stopping of Working Face

As the strip coal pillar is located between 02 and 03 working face, influences created by mining actives are doubly accumulated. The entire simulation is started with 02 working face excavation, and the various trail coal pillar widths were set after the stress calculation is in the equilibrium state. And then, the stopping of the 03 working face was initiated, hence the coal pillar stress variation behaviors can be analyzed when stopping is simulated at both working faces. Figure12 shows the stress distribution variations in coal body during the stopping in both 02 and 03 working faces.Figure 12
Vertical stress distribution diagram of coal pillars with different widths of 02# and 03# working faces during mining.
(a)(b)As shown in Figure12, the coal pillar width variation does not affect the stress distribution much during the 02 working face excavation. When coal pillar width increased from 3 to 30 m, the corresponding peak vertical stresses are always in-between 40 and 50 MPa, but the pillar stress distribution changed as the peak stress zone shifted to the 02 goaf zone.In comparison, the coal pillar width variation has a relatively big influence on stress distribution at the 03 working face. The 3 m wide coal pillar was crushed, and its maximum stress concentration shifted to the roadway. Whereas, when coal pillar widths are 5, 8, 12, 15, 20, 25, and 30 m, the corresponding maximum stress and stress concentration coefficients are 71 MPa, 4.4; 82 MPa, 5.1; 145 MPa, 9; 116 MPa, 7.2; 89 MPa, 5.5; 73.6 MPa, 4.6; 69.3 MPa, 4.3; and 66.8 MPa, 4.1, respectively. An increasing trend is observed in the stress concentration factor when the pillar width increases from 5 to 10 m. When coal pillar width increases from 10 to 12 m, its concentration coefficient remains over 7, and it poses extreme danger to roadway stability; On the contrary, the stress concentration coefficient drops when pillar width increases from 10 to 30 m. Especially for coal pillar greater than 20 m, its corresponding stress concentration remains below 5, which is beneficial to the pillar stability.It is found that when coal pillar thickness is 15 m, the corresponding pillar stress is smaller than that of 10 and 12 m thick pillars, but still at a relatively high level. Therefore, the numerical simulation proved that pillar width between 5 and 8 m benefits its stability the most due to the relatively small corresponding peak vertical stress.


## 5.3. Stress Characteristics of Coal Pillar with Different Widths

The pillar stress characteristics were analyzed based on a representative scenario when the 03 working face is advanced by 130 m, as well as the strip pillar is influenced by both 02 goaf and 03 working face excavations. The vertical stress vs pillar width curves are shown in Figure13. It can be seen from Figure 13 that all the pillar stress distributions present a unimodal correlation and all peak stresses are shifted to the goaf side with a shifting range of 0.5–3 m. Meanwhile, vertical stress is always the biggest and the difference in-between both horizontal stresses is barely noticed. However, the horizontal stress, which is parallel to the roadway alignment, is still slightly greater than that of perpendicular to the roadway alignment. When coal pillar width is in-between 3 and 5 m, it was crushed and depressurized. When pillar width increases from 5 to 15 m, both vertical stress and alignment-parallel horizontal stress are concentrated. When pillar width further increases from 20 to 25 m, a bimodal correlation appeared and the peak stress close to the 02 working face is slightly higher than that of the 03 working face. It is also found that the stress at the pillar center is higher than the original rock stress, and the alignment-perpendicular horizontal stresses throughout the entire strip coal pillar are all smaller than the original rock stress.Figure 13
Stress distribution in coal pillar with different widths. (a) Vertical stress in pillar. (b) Horizontal stress in pillar (axial direction of parallel roadway). (c) Horizontal stress in pillar (vertical roadway axial direction).
(a)(b)(c)


## 5.4. Characteristics of Roadway Stress Field with Different Pillar Widths

The internal stress of the coal pillar is monitored by installing borehole stressmeter, and then the maximum stress in the coal pillar and solid coal of roadway side is extracted. The relationship between the characteristics of the roadway peak vertical stress and pillar width when the 03 working face was advanced by 130 m is shown in Figure14. It can be seen from Figure 14 that the pillar width variation has certain influence on the vertical stresses in both coal pillar and coal rib. Figure 13 reveals that the curves of both peak vertical stress vs coal pillar width present a bimodal trend. With the increase of coal pillar width, the vertical pillar stress increases sharply and then drops mildly, whereas the lower coal rib stress slightly increases and then drops. Figure 14 also indicates that coal pillar bearing capacity is positively correlated to the pillar width, and the lower coal rib stress has a pillar-side shifting trend with the increases of pillar width. When coal pillar width reaches 8 m, the upper rib coal pillar stress becomes much closer to that of lower coal rib, and until then all the lower coal rib stresses are slightly greater than pillar stress. With a further increment in coal pillar (width ≥ 10 m), the critical pillar stress remains much greater than that of lower coal rib stress, and it becomes the major load bearing object.Figure 14
Stress distribution of roadway under different pillar widths.Overall, peak pillar stress and peak lower-rib stress appear together when the coal pillar is 10 m wide, and the difference in-between both upper and lower-rib peak stresses is minor due to the side supporting effect created by 02 goaf.


## 5.5. Characteristics of Roadway Displacement Field with Different Pillar Widths

The analyses of roadway country-rock deformation at 200 m ahead of the open-off cut focus on roof-floor convergence and rib-rib convenience. As shown in Figure15, the accumulation of roof-floor convergence is inversely correlated to the working face distance, and it rapidly builds up when the working face approaches from 50 m away. The relationship between coal pillar width and relatively roof-floor displacement also presents a bimodal correlation. When coal pillar width increases from 5 to 10 m, the relatively roof-floor displacement also increases and reaches its peak value. And then, it gradually drops when coal pillar width increases from 10 to 30 m.Figure 15
Relative displacement curve of roadway roof and floor at different pillar widths. (a) Vertical displacement curve. (b) Horizontal displacement curve of two sides.
(a)(b)Regarding the rib-rib deformation, the results also reveal that its relative displacement is negatively correlated to the working face distance, and this trend intensifies further when the working face is within 40 m.Finally, based on the analyses of stress-displacement distributions in both coal pillar and coal rib, the strip pillar width is therefore recommended as 5–8 m.


## 6. Analysis of Stress Distribution of Coal Pillars with Different Widths When Excavating Goaf in the 03 Working Face

### 6.1. Stress Variation Rule of Coal Pillar in Working Face with Widths of 5 m, 10 m, and 15 m

In order to see the stress distribution of the working face more comprehensively, when the coal pillar widths were 5 m, 10 m, and 15 m, the stress distribution diagrams were selected. In addition, when the working face of 03 was stopping at 130 m, it can be seen from Figure16 that the stress distribution forms of coal pillars in the mining area under three widths are different. The goaf formed by working face 02 will still affect the stress distribution of coal pillar to a certain extent. Under the influence of mining in working faces 02 and 03, concentrated stress will be formed, and the maximum stress is located in a section behind the goaf. At the same time, under the action of advanced mining stress and lateral stress of the working face, the superposition of stress at the intersection of solid coal side and the pillar of the working face forms a high stress zone. With the decrease of the width of coal pillar, the stress concentration in the middle of coal pillar and the end of the working face on one side of coal pillar is increasing, and the distance between the stress concentration area and the roadway side is also getting closer.Figure 16
Stress nephogram at 130 m stope. (a) Width of pillar is 5 m. (b) Width of pillar is 10 m. (c) Width of pillar is 15 m.
(a)(b)(c)


### 6.2. Stress Characteristics of Coal Pillar with Different Widths

Figure17 shows the stress characteristics of coal pillar when working face pushes 50 m. As shown in Figure 17, the stress in 5 m and 10 m coal pillars presents a single peak distribution; the peak position is basically in the center of the coal pillar, and the vertical stress in 15 m coal pillar presents a bimodal distribution. The vertical stress in the coal pillar is the largest, the difference between the two horizontal stresses is not big, but the overall horizontal stress parallel to the roadway axis is slightly larger than the horizontal stress perpendicular to the roadway axis. With the increase of the width of the pillar, the concentrated stress peak in the pillar first increases and then decreases, and the stress value is the largest in the pillar with a width of 10 m.Figure 17
Stress distribution in coal pillar with different widths. (a) Vertical stress in pillar. (b) Horizontal stress in pillar (axial direction of parallel roadway). (c) Horizontal stress in pillar (vertical roadway axial direction).
(a)(b)(c)


### 6.3. Characteristics of Roadway Displacement Field with Different Pillar Widths

Fixed the position of 200 m in front of the 03 working face cut, and studied the top and bottom movement and convergence of the two sides in the mining When the mining distance reaches 200 meters in the front of the open-off cut, both the roof-floor's movement displacement and two sidewalls' movement displacement of different roadways for different coal pillar widths have been studied, respectively, in mining process. Figure18 shows the relative vertical displacement curves of roadway roof and floor at different pillar widths, and Figure 19 shows the relative horizontal displacement curves of two sides of roadway at different pillar widths.Figure 18
Relative vertical displacement curve of roadway roof and floor at different pillar widths.Figure 19
Horizontal displacement curves of two sides of roadway at different pillar widths.It can be seen from Figure18 that with the decrease of distance from the working face, the moving amount of top and bottom increases and the trend accelerates. When the width of coal pillar increases from 5 m to 15 m, the relative displacement of the roof and floor first increases and then decreases. When the width of coal pillar changes from 5 m to 10 m, the relative displacement of roof and floor gradually increases. When the width of coal pillar changes from 10 m to 15 m, the displacement reaches the maximum.As can be seen from Figure19, along the strike, with the decrease of distance from the working face, the displacement value gradually increases and the curve is concave. When the width of coal pillar increases from 5 m to 10 m, the relative horizontal displacement of two sides of the roadway increases. When the width of coal pillar increases from 10 m to 15 m, the horizontal displacement of both sides of roadway decreases. However, the deformation of the roadway with a width of 5 m is slightly larger than that of the roadway with a width of 15 m. Considering that the coal pillar with a width of 10 m can be mined, it is suggested that the reasonable width of the coal pillar at the side of goaf should be less than 10 m. Combined with the previous field monitoring and referring to the small pillar retention situation of similar mines, the reasonable pillar width of BayanGaole coal mine is finally determined to be 6 m. At present, the working face has been successfully stopped and produced obvious economic benefits.


## 6.1. Stress Variation Rule of Coal Pillar in Working Face with Widths of 5 m, 10 m, and 15 m

In order to see the stress distribution of the working face more comprehensively, when the coal pillar widths were 5 m, 10 m, and 15 m, the stress distribution diagrams were selected. In addition, when the working face of 03 was stopping at 130 m, it can be seen from Figure16 that the stress distribution forms of coal pillars in the mining area under three widths are different. The goaf formed by working face 02 will still affect the stress distribution of coal pillar to a certain extent. Under the influence of mining in working faces 02 and 03, concentrated stress will be formed, and the maximum stress is located in a section behind the goaf. At the same time, under the action of advanced mining stress and lateral stress of the working face, the superposition of stress at the intersection of solid coal side and the pillar of the working face forms a high stress zone. With the decrease of the width of coal pillar, the stress concentration in the middle of coal pillar and the end of the working face on one side of coal pillar is increasing, and the distance between the stress concentration area and the roadway side is also getting closer.Figure 16
Stress nephogram at 130 m stope. (a) Width of pillar is 5 m. (b) Width of pillar is 10 m. (c) Width of pillar is 15 m.
(a)(b)(c)


## 6.2. Stress Characteristics of Coal Pillar with Different Widths

Figure17 shows the stress characteristics of coal pillar when working face pushes 50 m. As shown in Figure 17, the stress in 5 m and 10 m coal pillars presents a single peak distribution; the peak position is basically in the center of the coal pillar, and the vertical stress in 15 m coal pillar presents a bimodal distribution. The vertical stress in the coal pillar is the largest, the difference between the two horizontal stresses is not big, but the overall horizontal stress parallel to the roadway axis is slightly larger than the horizontal stress perpendicular to the roadway axis. With the increase of the width of the pillar, the concentrated stress peak in the pillar first increases and then decreases, and the stress value is the largest in the pillar with a width of 10 m.Figure 17
Stress distribution in coal pillar with different widths. (a) Vertical stress in pillar. (b) Horizontal stress in pillar (axial direction of parallel roadway). (c) Horizontal stress in pillar (vertical roadway axial direction).
(a)(b)(c)


## 6.3. Characteristics of Roadway Displacement Field with Different Pillar Widths

Fixed the position of 200 m in front of the 03 working face cut, and studied the top and bottom movement and convergence of the two sides in the mining When the mining distance reaches 200 meters in the front of the open-off cut, both the roof-floor's movement displacement and two sidewalls' movement displacement of different roadways for different coal pillar widths have been studied, respectively, in mining process. Figure18 shows the relative vertical displacement curves of roadway roof and floor at different pillar widths, and Figure 19 shows the relative horizontal displacement curves of two sides of roadway at different pillar widths.Figure 18
Relative vertical displacement curve of roadway roof and floor at different pillar widths.Figure 19
Horizontal displacement curves of two sides of roadway at different pillar widths.It can be seen from Figure18 that with the decrease of distance from the working face, the moving amount of top and bottom increases and the trend accelerates. When the width of coal pillar increases from 5 m to 15 m, the relative displacement of the roof and floor first increases and then decreases. When the width of coal pillar changes from 5 m to 10 m, the relative displacement of roof and floor gradually increases. When the width of coal pillar changes from 10 m to 15 m, the displacement reaches the maximum.As can be seen from Figure19, along the strike, with the decrease of distance from the working face, the displacement value gradually increases and the curve is concave. When the width of coal pillar increases from 5 m to 10 m, the relative horizontal displacement of two sides of the roadway increases. When the width of coal pillar increases from 10 m to 15 m, the horizontal displacement of both sides of roadway decreases. However, the deformation of the roadway with a width of 5 m is slightly larger than that of the roadway with a width of 15 m. Considering that the coal pillar with a width of 10 m can be mined, it is suggested that the reasonable width of the coal pillar at the side of goaf should be less than 10 m. Combined with the previous field monitoring and referring to the small pillar retention situation of similar mines, the reasonable pillar width of BayanGaole coal mine is finally determined to be 6 m. At present, the working face has been successfully stopped and produced obvious economic benefits.


## 7. Conclusions

Coal pillar width is the major factor that affects pillar stability and roadway maintenance. It dominates the horizontal distance between roadway and stopping space, and affects the bearing pressure under excavation and coal pillar bearing capacity. The coal pillar’s ultimate bearing capacity not only depends on its boundary conditions and mechanical properties but also relies on its geometry. This paper revealed and analyzed the in-situ stress distribution of strip coal pillar in Bayangaole mine by utilizing rock mechanics testing, borehole stressmeter, PASSAT monitoring system, and numerical simulation. The small pillar excavation was practiced in Bayangaole mine and 5 essential conclusions are drawn.(1)
Based on the distribution patterns of acoustic wave velocity throughout coal pillar at the 311102 working face, it is discovered that the impact created by mining activity expands into the working face around 80 m, where the peak stress concentration is found about 15 m ahead of the working face.(2)
Mildly Impacted, Ordinarily Impacted, and Severely Impacted zones are distinguished based on the borehole stressmeter testing results in Bayangaole mine. Mildly impacted zone is more than 24 m away from the working face. Ordinarily impacted zone is 8–24 m away from the working face, and the severely impacted zone is within 8 m from the working face. The peak coal body in-situ stress appears around 5 m ahead of the working face, whereas the in-situ stress is extremely low at the working face.(3)
The most common failure modes, when the coal pillar stress and deformation reach their ultimate bearing capacity of deformation limit are shear failure, splitting failure, rheodestruction failure, sliding failure along with weak plane and bedding, and heaving failure when coal pillar is punched into the floor. Based on the ultimate equilibrium theory, the plastic zone width of the 311102 working face coal body in Bayangaole mine is determined as 4.96 m.(4)
The numerical simulation conducted both stress and displacement analyses on coal pillars with various widths. The results indicated that 5 m and 10 m wide pillars showed a unimodal stress-width correlation and the peak vertical stresses are all located at the pillar center, whereas 15 m wide pillar has a bimodal stress-width relationship. In comparison, 10 m wide pillar holds the maximum in-situ stress.(5)
In consideration of in-situ stress test results, sonic wave velocity test results, theatrical analysis, numerical simulation analysis, actual site conditions, and economic effects, 6 m wide pillar is designed at the working face 311102. As a result, stopping was successfully completed and remarkable economic benefits were achieved.



---

*Source: 8837189-2021-08-10.xml*