Abstract
Protective seam mining is an effective gas pressure relief method in deep mining. Effective theoretical calculation methods in the current studies on the prediction of pressure relief protection effect of interbed coal and rock masses and their distribution laws are lacking. Thus, the evaluation and research with respect to pressure relief effect in protective seam mining relatively lag behind. This situation restricts the engineering feasibility evaluation and decision making in the protective seam mining. Therefore, the influence of upper protective seam mining on the pressure relief protection effect of coal and rock mass between underlying beds was investigated in this study. On the basis of an analysis of concrete engineering projects, a mechanical model was constructed for the pressure relief protection effect of upper protective seam mining on the coal and rock mass between underlying beds. The distribution equation of pressure relief expansion ratio in the underlying protected seam was also derived. The influence laws of main influencing factors on the pressure relief protection effect of the protected seam were revealed as well. In the end, the pressure relief effect was analyzed and verified for the protected seam before and after mining through numerical simulation and similarity simulation test. The pressure relief effect of upper protective seam mining on the coal and rock mass between underlying beds and the distribution characteristics were deeply explored in this study, which could provide a theoretical reference for the decision making in the gas extraction engineering design and pre-evaluation of extraction effect. Results show that the effective pressure relief zone (expansion rate>0.3%) of the protected seam beneath the goaf is located within the range of approximately 40 m from the coal wall to the rear part. It also presents an approximate “Λ-shaped distribution characteristic,” that is, it experiences migration and evolution with the advancement in the working face. Moreover, the peak pressure relief lags behind the coal wall on the working face by nearly 10–20 m. In the numerical simulation, the expansion ratio in the goaf also presents an approximate “Λ-shaped distribution.” Its effective pressure relief zone is the 50 m range from the coal wall to the rear part of the goaf, and the peak value lags behind the coal wall by around 15 m. The theoretical results and numerical simulation results are basically consistent with the physical experiment results. The expansion rates are 1.25%, 1.268%, and 1.32%, respectively. The elastic modulus E of coal seam and interbed spacing H are the main influencing factors of the swelling deformation and are negatively correlated with the expansion ratio. In the actual mining process, E and H of the protected layer can be measured to infer the expansion deformation of the protected layer.
1. Introduction
The deep high ground stress and coal and gas outburst disaster become increasingly severe as the coal resource development and utilization extend to deep parts. During the deep mining process, the mining disturbance behaviors have a great bearing on the changes in the mechanical state of deep coal body and seepage field. Protective seam mining is an effective regional gas pressure relief method, and a large quantity of pressure relief gas at the protected seam is emitted into the working face and goaf via the floor cracks in the mining process; thus, extracting the pressure relief gas in the protected seam is necessary [1–5].
Domestic (Chinese) and foreign scholars have always been dedicated to the studies on protective seam pressure relief mining and coal and gas disaster. Wang et al. [6] studied the evolution laws of floor cracks in the short-distance protective seam mining. They demonstrated that the coal mass in the pressure relief and permeability improving zone had large swelling deformation and high permeability, and this zone was also an efficient pressure relief gas intercept and extraction zone. Ren et al. [7] acquired the load-carrying condition of the protected seam through the similar material simulation test and determined the stress state after the coal seam deformation through the force analysis of the protected coal seam. They concluded that greater mining disturbance borne by the protected seam meant greater peak breaking stress of the coal mass and higher volumetric strain during the upper protective seam mining process under the same test conditions. Xue et al. [8] obtained the expression of internal displacement of the rock mass using a semi-infinite body model. They also constructed a “two-zone” crack distribution model and its simplified mechanical model to probe the pressure relief mechanism of protected seam. On the basis of a study of gas emission laws on the short-distance protective seam mining face, Wang et al. [9] optimized the pressure relief gas extraction parameters of protective seam. They obtained that the gas emissions on the protective seam mining face mostly came from the pressure relief gases at the protected seam. Zhang et al. [10] discussed the local stress concentration and rock breaking behavioral mechanisms in rock intercalation after the protective seam mining from longitudinal and transverse perspectives via numerical simulation and mechanical analysis. Under the engineering background of 10# mine in Pingdingshan, Zhang et al. [11] studied the crack evolution laws and distribution characteristics in the coal and rock mass between two coal seams. They found that the cracks were expanded with the increase in intensity of mining activity, and most cracks developed at low angle or became parallel to the stratum. Through physical model and numerical simulation, Wang et al. [12] investigated the stress distribution and crack evolution laws in surrounding rocks during the protective seam mining process under certain mining conditions. They stated that the displacement of surrounding rock and the evolution of cracks were affected by the support pattern.
At present, extensive research has been conducted at home and abroad on the development of overburden fissure zones and gas pressure relief caused by traditional protective coal seam mining, and important research results have been obtained [13–15]. Good research progress has also been made in the expansion trend area below the mined-out area. However, with the continuous advancement of the protective layer working face, there is currently no fixed theoretical formula for the expansion rate and expansion range of coal and rock masses with different buried depths under the goaf. It is impossible to grasp the expansion and deformation at different positions of the coal and rock mass underneath the goaf during the advancement of the working face, nor can it accurately grasp the swelling area of the coal and rock mass underneath. The previous engineering experience or similar working face treatment methods are often used to determine the expansion changes of the underlying coal body, often with large deviations. Although there are related mechanical models to analyze the interaction between the protective layer and the protected layer, it still fails to reveal the specific factors that affect the expansion and deformation of the protected layer. Therefore, it is difficult to understand the real-time change law of the expanded range of the mined-out area. This paper takes the 11129 working face of Zhangji Coal Mine of Huainan Mining Group as the engineering background and explores the changing law and range of the expansion ratio of the protective layer under the goaf. The stress evolution process and the effective expansion zone change of the protected coal seam are theoretically deduced. The analytical solution of the expansion rate of the protected layer and its effective pressure relief range in the mining mode of the protective layer is given. It is verified by a combination of numerical simulation and similar simulation experiment. The research is expected to guide the engineering design and pre-evaluation decision making of protective layer mining and gas drainage pressure relief effects.
2. Analysis of Project Profile
The 8# coal in east (1) mining area of Zhangji Coal Mine has complicated geological conditions and high risk of coal and gas outburst. Thus, it does not meet the direct mining conditions. Upper protective seam mining is adopted to realize pressure relief and gas extraction of 8# coal. Meanwhile, 9# coal is mined at protective seam, where the average thickness of 9-1# coal is 1.9 m and that of 9-2# coal is 0.9 m. The average thickness of dirt band is 1.4 m, and the upper roof consists of siltstone with an average thickness of 8.0 m. The immediate roof is composed of quartz sandstone with an average thickness of 7.6 m. The lithology of immediate floor is mudstone with an average thickness of 2.5 m. The lithology of lower floor is siltstone with an average thickness of 5.2 m. The average thickness of 8# coal at the underlying protected seam is 4 m, and it is 9 m away from 9# coal. The recoverable strike and inclined lengths of 11129 working face, which is the first mining face of 9# coal, are 1,200 and 240 m, respectively. Its burial depth from the ground is 710 m. The general occurrence state of coal seam is a monoclinal structure, which is locally developed with secondary folds. The coal seam is generally high in the west and low in the east, and the dip angle ranges from 2° to 6°. The comprehensive histogram of the working face is shown in Figure 1.

8# coal at the protected seam is of low permeability, low gas drainage rate, and complicated geological conditions. Thus, the upper protective seam mining should be conducted for 8# coal to realize gas extraction and pressure relief. (1) The coal pillars reserved in the 9# coal goaf form displacement and stress fields in the mining process of 8# coal, which brings difficulty for the supporting and roadway layout of 8# coal working face and impedes the safety mining and gas extraction of 8# coal. (2) During the 8# coal mining process at protected seam, the existence of 9# coal goaf will accelerate the fracture of bearing bed. As a result, the 9# coal goaf runs through 8# coal working face. Consequently, gas, water, and gangues in the goaf at upper coal seam are poured into 8# coal working face, and this condition harms the safety production of this coal mine [16]. (3) The pressure relief of 8# coal is good for gas extraction. Thus, the pressure relief zone and gas extraction engineering are closely combined to realize safe and efficient mining of 8# coal and maximize the economic benefit.
3. Mechanical Model Analysis of Pressure Relief Protection Effect of Upper Protective Seam Mining on Coal and Rock Masses between Underlying Beds
3.1. Mechanical Model of Mining-Induced Stress Distribution on the Floor of Working Face
Coal seam is usually under elastic deformation state before mining; in the mining process, the initial stress field of rock on the floor is affected by the mining, the floor in front of the working face is under the action of concentrated stress, the floor stress field is changed and redistributed, and the floor stress distribution depends on the transfer of concentrated stress in front of the working face to the lower rock mass of the coal seam floor [17, 18]. The following assumptions are considered: the coal and rock mass under the initial stress state of rock do not influence the stress redistribution of underlying coal and rock mass, the abutment pressure in front of the working face can be regarded as the triangular banding load from the coal wall to stress peak point and the trapezoid banding load in front of the stress peak point, the floor is regarded as a uniform elastomer, and the problem is solved by a plain strain problem. The underlying rock and coal stress is then theoretically calculated [18].
To facilitate the computational analysis, the stress change in the coal mass and goaf is expressed in the form of increment:
The initial stress of rock is deducted from formula (1), and the distribution law of the stress increment at the coal side can be obtained. The maximum value of its stress increment is (k − 1)P, while that in the goaf is −P. Thus, the distribution law of abutment pressure increment in coal body on the working face and the goaf can be acquired [19], as shown in Figures 2 and 3.


The expression of vertical stress is [19]
According to the equilibrium of elastic mechanics, the equilibrium can be reached only when the following condition is satisfied:
The following condition is then solved:
To solve the stress at one point M within the half planar body, the coordinate axes are taken as shown in Figure 3, and the coordinates of point M are (x, y). A minimum length dξ is taken at the place with a distance of from the origin of coordinates O on the straight line Oy. The force dP = qdξ it bears is regarded as a minimal concentrated force, and the vertical and horizontal distances from point M to this minimal concentrated force dP are x and y-ξ, respectively. Thus, the stress caused by the minimal concentrated force dF = qdξ at point M is as observed in the following formula:
The stresses caused by all minimal concentrated forces are superposed to solve the vertical and horizontal stress distributions in front of the working face. Therefore, the integral equation of formula (5) is obtained as follows:
The vertical and horizontal stresses are solved using formula (6).
The equation for the mining-induced stress distribution in the surrounding rocks on the floor perpendicular to the working face is
The following condition is then solved:
The equation for the mining-induced stress distribution in the surrounding rocks on the floor in the advancement direction of the working face is
The following condition can be solved:
3.2. Mechanical Relation of Upper Protective Seam Mining with the Pressure Relief Protection Effect of Coal and Rock Masses between Underlying Beds
According to the physical and geometric equations in the plane strain problem of the elastic mechanics, the vertical displacement can be solved as follows:where m and n are stress coefficients of P(ξ) corresponding to d1, d2, d3, and d4; E is the elastic modulus; μ is Poisson’s ratio; and K is a constant.
The floor stress distribution in the advancement direction of the working face is divided into four stress zones: AB, BO, OC, and CD. The integrals of the abovementioned formulas are solved. Then, complete stress distribution, strain, and displacement equations in vertical direction of the floor can be obtained.
In the protective seam mining process, the displacements of 8# coal roof and floor are calculated using the aforementioned formula, and the expansion ratio Q is further calculated aswhere Δu is the displacement difference between roof and floor of the protected seam, m; h is the average thickness of protected seam, m; and a is the distance from protective seam to the roof of protected seam, m.
The fluctuation of expansion ratio Q can reflect the evolution law of pressure relief at the protected seam, and Q>0.3% is taken as an effective pressure relief zone.
4. Influence Factor Analysis of Upper Protective Seam Mining on Pressure Relief Protection Effect of Protected Seam
The pressure relief effect of the protected seam is mainly related to stress concentration coefficient k, elastic modulus E of coal seam, coal seam thickness h, and interbed spacing H. According to the mine ground pressure laws and geological structural conditions, the engineering parameter k is taken as 2.5, μ value is 0.3, E = 1 GPa, h = 4 m, H = 9 m, d1 = 30 m, d2 = 10 m, d3 = 10 m, and d4 = 40 m. After external and human factors are excluded, different parameters are taken to comparatively analyze the changes in the roof stress and expansion ratio of 8# coal, as shown in Table 1 and Figures 4–11.








As shown above, the concentrated stress coefficient k mainly influences the abutment pressure in front of the coal wall, and it is in direct proportion to the abutment pressure. With the continuous increase in the coefficient k, 8# coal beneath the abutment pressure goes through swelling deformation, which has a very minor influence on the stress and swelling deformation in the rear goaf. Under k = 2.5, the maximum increment of abutment pressure in front of the coal wall is 18.5 MPa. As the elastic modulus is continuously increased, no influence is generated on the stress change of the protected seam, and the expansion ratio is gradually reduced. Moreover, the elastic modulus becomes inversely proportional to expansion ratio, and it is the result of joint action of coal seam properties and surrounding environment. Before the protective seam mining, the elastic modulus of the protected seam should be measured in advance to estimate the overall swelling deformation of the protected seam during the protective seam mining process. During the upper protective seam mining process, the change in the thickness of overlying coal seam has no obvious influence on the overall stress change or swelling deformation. As a result, the thickness of coal seam is not the primary influencing factor of the swelling deformation. With the increase in the burial depth of the protected seam, the downward stress propagation is continuously reduced with the increase in interbed spacing, the swelling deformation is gradually mitigated, and the peak swelling deformation is gradually distant from the coal wall. The scopes of influence at two sides are also slightly enlarged, and both approach 0 in the end.
In summary, the main influencing factors of the swelling deformation of the protected seam are elastic modulus E of the coal seam and interbed spacing H. Other factors exert minor influences. In the actual mining process, the layer spacing H can be estimated, and the underlying coal rock mass is sampled by drilling holes into the protected layer. The mechanical properties of the rock samples were tested using the MTS rock mechanical performance testing machine, and the stress-strain curves of the rock samples were obtained, and then the elastic modulus E of the underlying coal rock mass was obtained.
5. Instance Analysis of Upper Protective Seam Mining Project in Zhangji Coal Mine
5.1. Example Calculation of Protective Seam Mining on 11129 Working Face of Zhangji Coal Mine
According to the geological conditions of 11129 working face and mine ground pressure law, the distance from the coal wall of the working face to the peak stress concentration is d2 = 10 m, the distance from the peak stress concentration to the initial stress zone of rock in deep coal and rock mass in the front is d1 = 30 m, the spatial width of the working face is d3 = 10 m, the distance from the goaf to the initial stress zone of rock in the rear deep compacted goaf is solved as d4 = 40 m, the stress concentration coefficient k is taken as 2.5, the initial stress of rock is 19 MPa, and the distance from the working face to the origin is 40 m. The aforementioned formula was used to calculate the displacement change of the roof and floor of the protected seam beneath the floor in the advancement direction of the working face. Mathematical software Wolfram Mathematica was used to solve the OA, AB, BC, and CD segments. The total vertical stress was obtained through the superposition, and the integral was taken to solve the vertical displacement.
The expression of vertical displacement in different phases is
The vertical stress of the protected seam iswhere
The roof and floor displacement of the protected seam iswhere
The change in the roof and floor displacement is calculated. Then, the formula of expansion ratio Q at the protected seam is further solved as
According to the aforementioned calculation results, the change charts of the roof and floor displacement of the protected seam and the expansion ratio are drawn as shown in Figures 12 and 13.


The distribution environment of roof and floor displacement changes in 8# coal was analyzed, as shown in Figure 12. As the upper protective seam mining was influenced by the abutment pressure concentration on the front coal seal, 8# coal undertook increasing vertical pressure. Its roof and floor displacements reached the maximum values of 0.32 and 0.26 m, respectively. With the stress relief of surrounding rocks in the goaf behind the coal wall, underlying 8# coal experienced upward swelling deformation. The maximum deformation of the top plate is 0.34 m at 15 m behind the coal wall, and the maximum deformation of the bottom plate is 0.289 m at 18 m behind the coal wall. The expansion deformation value of the top and bottom plates of the protected layer gradually attenuates as it approaches and moves away from the coal wall and eventually approaches zero.
On the basis of Figure 13, the change curve chart of expansion ratio at the protective seam beneath the floor of the working face was analyzed. The protective seam mining facilitated the stress relief of surrounding rocks on the floor behind the coal wall and was influenced by the downward transfer of abutment pressure on the front coal wall. The coal and rock mass on the working face and that beneath the goaf went through swelling deformation, and the expansion ratio of underlying coal and rock mass reached the maximum value of 1.25% at the place 12 m behind the coal wall. With the advancement in the working face, the expansion ratio of overlying coal and rock mass was gradually reduced and approached 0. The effective scope of expansion (expansion ratio>0.3%) was from the coal wall to the place 40 m behind the coal wall, and the floor heaving was the most obvious in this scope. Meanwhile, the pressure relief effect of the protected seam was the best.
5.2. Simulation Analysis of Protective Seam Mining Effect on 11129 Working Face in Zhangji Coal Mine
A numerical model was established via FLAC3D numerical simulation software, with dimensions of 540 m × 1000 m × 169 m (Figure 14). The mechanical properties of rock samples were tested by MTS rock mechanical property testing machine, and the stress-strain curve of rock samples was obtained. The elastic modulus can be solved by slope. In order to ensure the reliability of data, the average value of three lithologic experiments in each group is taken. The mechanical parameters of the roof and floor are shown in Table 2. Fixed horizontal displacement constraints were set at the bottom edge, front, back, left, and right boundaries of the model. The Mohr–Coulomb failure criterion was also adopted to analyze the mechanical characteristics of coal and rock mass in this model. Domestic and foreign scholars have conducted several experiments and studies on the abutment effect of equivalent filled rock mass in goafs [20–24]. In the present study, the goaf was simulated with double-yield model. The stress-strain relation of the rock mass in the goaf is presented in Table 3, and the mechanical parameters of the rock mass in the double-yield goaf are presented in Table 4.

The advancement in 11129 working face from the open-off cut along the strike (positive direction of axis y) was simulated, the open-off cut was nearly 200 m from the boundary, and the influence of boundary effect was eliminated. The excavation step size was 25 m, the goaf was filled, balanced, and advanced by 600 m along the strike, and the mining was stopped at y = 800 m. The displacement changes of 8# coal roof and floor at 9 m beneath it during the mining process are shown in Figures 15–19.





Figures 15–17 shows that with the advancement in the working face, a goaf was formed, the stress was redistributed, and an abutment pressure zone higher than the initial stress of rock was formed in the coal and rock mass around the goaf. The surrounding coal and rock masses subsided, and the overlying coal and rock mass in the goaf experienced upward swelling deformation. The overall displacement was large at two sides and small in the middle, and that in the middle presented a waved change. The maximum floor heaving amounts of underlying 8# coal roof and floor in the goaf were 0.268 and 0.215 m, respectively. The maximum abutment pressure borne by the coal wall around 8# coal beneath the goaf was 37.4 MPa, which was higher than the initial stress (18.4 MPa) of rock. The marginal position of the goaf carried the minimum stress, and the stress in the middle of the goaf was gradually increased [25–27].
As shown in Figures 18 and 19, the overall central displacements of the roof and floor fluctuated within 0.1–0.2 m and 0.07–0.15 m, the overall average expansion ratio was 1.286% when the mining face was advanced by 600 m, and the scope of influence of average expansion was from the coal wall to 50 m behind the coal wall. The peak pressure relief lagged behind the coal wall by nearly 20 m, and it presented an “Λ-shaped distribution” on the whole. With the advancement in the working face, the goaf was gradually compacted, and local rebounding and recompaction process existed within a small scope.
5.3. Similarity Simulation Test Analysis of Protective Seam Mining Effect of 11129 Working Face in Zhangji Coal Mine
The similarity simulation test was conducted under the engineering background of 9# coal at the upper protective seam and 8# coal at the protected seam in Zhangji Coal Mine. The model testbed was made of steel frame with dimensions of 3.0 m (length) × 1.2 m (height) × 0.3 m (width). The model test should satisfy the following conditions: geometric similarity, kinematic similarity, dynamic similarity, similarity of boundary conditions, and proportionating physical quantities. Therefore, the linear ratio was taken as 1 : 100, volume weight ratio was 3 : 5, and time ratio was 1 : 12. The ratio table of the main layers of the experimental model is shown in Table 5. The model and the distribution of measuring points are shown in Figures 20 and 21.


Strain gauges were laid on the protected seam and deformation monitoring points were arranged to perform real-time monitoring of 8# coal pressure and deformation of coal and rock mass. A total of 8 strain gauges are arranged in the horizontal direction, the middle 2#–7# strain gauges are spaced 30 cm apart, and the two ends of the strain gauges 1# and 8# are spaced 30 cm apart from the end of the steel frame. The strain change curves at measuring points 1#–8# are shown in Figures 22 and 23.


As shown in Figures 22 and 23, the protective seam mining resulted in the stress relief of the underlying protected seam, and the roof and floor of the protected seam experienced upward swelling deformation. As the working face was advanced from the open-off cut to 40 m before the mining stopping line, the maximum roof and floor displacements were 0.267 and 0.216 m, respectively. The maximum expansion ratio (1.32%) was reached at 14.5 m behind the coal wall. The effective pressure relief zone was from the coal wall to 47 m behind the coal wall, which was approximate to the theoretical analysis and numerical simulation results.
5.4. Guidance and Suggestions for the Gas Extraction Design Scheme of Underlying Protected Seam
The 9# coal was close to the underlying 8# coal seam, and pressure relief gas of 8# coal (8# coal gas pressure was 1.52 MPa, and gas content was 5.5 m3/t) would surge into the goaf of 9# coal working face during the mining period of 9# coal working face. The extraction and drainage were simultaneously implemented in the upper protective seam mining, drilling was done toward upper strata in the floor roadway in the form of meshes, and the gas drainage of 8# coal seam was conducted within the effective pressure relief zone. Thus, the gas in the pressure relief zone was effectively prevented from permeating into the upper goaf, and the safe and efficient mining of the working face could be guaranteed.
6. Conclusions
(1)A mechanical model was constructed for the pressure relief protection effect of upper protective seam mining on the coal and rock mass between underlying beds. The equation for the expansion ratio distribution of underlying protected seam was also acquired. The pressure relief evolution laws of the protected seam before and after mining were further analyzed and verified through numerical simulation and similarity simulation test. The theoretical results are basically consistent with the numerical simulation results. The expansion rates are 1.25% and 1.268%, respectively, and the physical experimental results are 1.32%. The major reason for the results may be man-made reasons in the process of model making and mining.(2)The effective pressure relief zone of the protected seam was located within the range from the coal wall to nearly 40 m behind it, and it presented an approximate “Λ-shaped distribution.” With the advancement, migration, and evolution of the working face, the peak pressure relief lagged behind the coal wall by around 10–20 m. The numerically simulated expansion ratio in the goaf also presented an approximate “Λ-shaped distribution.” The effective pressure relief zone was from the coal wall to approximately 50 m behind the goaf, and the peak value lagged behind the coal wall of the working face by around 15 m.(3)The elastic modulus E of coal seam and interbed spacing H were the main influencing factors of the swelling deformation, and they were negatively correlated with the expansion ratio. Meanwhile, the influence of coal seam thickness on the overall selling deformation could be neglected. The interbed spacing was one of the primary influencing factors of the stress change of the protected seam. The change in the stress concentration coefficient k had an obvious effect on bearing the pressure of the coal and rock mass in front of the coal wall. However, it had no obvious pressure relief effect on the underlying protected seam in the rear goaf.(4)The expansion ratio and expansion zone of the protected seam were mastered in a real-time way, and the pressure relief degree and its scope of influence were reflected. With the advancement in the working face of protective layer, the peak expansion ratio in the goaf was also followed up timely. A good pressure relief zone was also provided for the gas extraction from the highly gassy coal seam. This study provides a theoretical reference for the decision making in the gas extraction engineering design and pre-evaluation of extraction effect.Data Availability
The data generated or analyzed during this study are included within this article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Jiaxin Dang and Qingwei Bu conducted theoretical calculation and analysis; Jiaxin Dang, Min Tu and Xiangyang Zhang conducted field data collection; Jiaxin Dang and Qingwei Bu conducted numerical simulation analysis; Jiaxin Dang wrote the manuscript.
Acknowledgments
This study was supported by the NSFC Project (52074008), Anhui Collaborative University Innovation Project (GXXT-2020-056), Inner Mongolia Natural Science Foundation Project (2019MS05055), and Open Research Fund Project of Key Laboratory of Safety and High-Efficiency Coal Mining of Ministry of Education (JYBSYS2019208).