A Numerical Investigation of the Stress Relief Zones Around a Longwall Face in the Lower Seam for Gas Drainage Considerations

Chunlei Zhang,Y.P.Chugh,Ruimin Feng,Yong Zhang,Wei Shen,Jingke Wu and Yushun Yang

1Faculty of Architecture and Civil Engineering,Huaiyin Institute of Technology,Huai’an,223001,China

2Mining and Mineral Resources Engineering,Southern Illinois University,Carbondale,Illinois,62901,USA

3School of Mechanics and Civil Engineering,China University of Mining&Technology(Beijing),Beijing,100083,China

4Department of Civil Engineering,University of Arkansas,Fayetteville,AR 72701,USA

5School of Energy and Mining Engineering,China University of Mining&Technology(Beijing),Beijing,100083,China

ABSTRACT Extraction of a protective coal seam(PVCS)-below or above a coal seam to be mined with the potential of coal and gas outburst risk-plays an important role not only in decreasing the stress field in the surrounding rock mass but also in increasing the gas desorption capacity and gas flow permeability in the protected coal seam(PTCS).The PVCS is mined to guarantee the safe mining of the PTCS.This study has numerically evaluated the stress redistribution effects using FLAC3D model for a longwall face in Shanxi Province.The effects of mining depth,mining height and inter-burden rock mass properties were evaluated using the stress relief angle and stress relief coefficient.Vertical stress distribution,stress relief angle and stress relief coefficient in the PTCS were analyzed as the face advanced in the PVCS.The results showed that the stress relief achieved in different locations of the PTCS varied as the face advanced.Sensitivity analyses on the pertinent variables indicate that the stress relief in the PTCS is affected most by the mining depth followed by the inter-burden lithology and the mining height.Furthermore,the elastic moduli of different layers within the inter-burden rock mass are more important than their uniaxial compressive strength(UCS)and Poisson’s ratio.These observations can guide gas drainage borehole design to minimize the accidents of coal and gas outbursts.

KEYWORDS Methane drainage;longwall mining;protective/protected coal seam;multiple seams;FLAC3D soft ware;stress relief zones

1 Background and Problem Statement

Over 50% of the coal mines in China mine coal seams with high gas content.The definition of high gas content coal mine is very complex but the commonly used indicators are the absolute gas emission (AGE) and the relative gas emission (RGE).When the AGE is greater than or equal to 40 m3/min,or the RGE is greater than or equal to 10 m3/t,the mine can be regarded as having high gas content [1].Similarly,about 44% of all coal mines in China are also known to have high coal and gas outburst risk [2].The coal and gas outburst risks seriously threaten coal mine safety and reduces operations productivity [3].The National Bureau of Statistics reported(2017) that 2,171 workers died from gas-related accidents in 2005 alone,and gas-related accidents accounted for about 64% of all coal mine accidents from 2006 to 2016 [4].Implementation of various comprehensive gas prevention and control measures in coal mines has greatly decreased fatal injuries in recent years.Since the gas content and gas pressure in coal seams potentially increases with the mining depth,gas and coal outburst accidents continue to threaten the safety of coal mines.

Most coal seams in China are characterized by low gas pressure,low gas saturation,low permeability due to their high degree of metamorphism.This has resulted in low success in gas drainage before coal mining [4–7].Many researchers have suggested that achieving stress relief in coal mass can improve gas permeability of coal and rock mass,which can enhance gas extraction,and reduce gas content and pressure [2,8–14].After decades of theoretical studies and field trials,the concept of excavating the PVCS has proved to be an effective method to reduce coal and gas outburst risk in the PTCS [13,15].The 139th article in the “Safety Regulations for Coal Mines”stipulates:“In the case of mining coal seam group with coal and gas outburst risk,mining protective layers must be adopted as a measure to prevent coal and gas outburst” [1].

2 Description of the Case Study Mine

The mine is located in NW China around Gaoping City of Shanxi Province (Fig.1).The current mine production is 5 Mt of raw coal annually.There are 16 coal seams in the area but #3 and #15 are the two primary minable seams with some portions of #2 and #8 seams also minable.The coal seams depths range from 300 to 800 m with strike in the E-W.Initial longwall panels were 400 m deep.The gas content in the coal field varies and it increases from about 0.4 m3/t in the eastern part to about 25 m3/t in the western part.The mine is currently extracting #3 low permeability coal seam with methane content of 6.5–15.09 m3/t.

Tab.1 shows the selected parameters of the major minable coal seams in the study area.The thickness of #2 and #3 coal seams varies 0–2.95 m and 4.67–6.58 m,respectively.Selected parameters of #3 and #8 coal seams are listed in Tab.2 [13].Seams are characterized by low permeability and low hardness.The gas content of #3 coal seam can reach up to 15.09 m3/t.Gas outburst and dynamic gas emissions frequently occur during the excavation of the roadways for#4 306 longwall faces.The #8 coal seam is a non-outburst prone coal seam with gas content of only 2–3 m3/t.In order to eliminate coal and gas outburst risk in #3 coal seam,the lower #8 coal seam was selected as the PVCS.The LW84306 (Longwall 84306) face (Fig.2) is located in#8 seam with panel width and length of 300 and 1128 m.The average mining thickness is 1.22 m and ranges from 1–3 m.The LW4306 (Longwall 4306) working face mines #3 seam of 5.23 m with an average dip angle of 0◦that ranges from 0–6◦.The panel length and width are about 1100 and 270 m respectively and mining depth ranges from 500 to 550 m.With the mining moving toward the west,the stress relief effect of the PVCS evolves with the changes in its mining height,and the properties of the inter-burden rock mass (rock mass between seam #3 and seam #8) and the mining depth.

Figure 1:Location of Changping coal mine in Shanxi Province

Table 1:Selected parameters of the major minable coal seams

Table 2:Selected parameters of No.3 coal seam

Figure 2:Layout of spatial locations of working faces

3 Literature Review

It is desired to have a good understanding of stress redistribution around a longwall face,and associated fracture development,and gas flow behavior during the mining process and then effectively apply the PVCS concept in the field.A few studies on these topics have been conducted under the condition of extracting single seam and multiple seams [13–15].Formation and evolution of gas-flow channels have been investigated around the front and side abutment pressure areas,floor and gob behind the coal face area during longwall mining [16–21].Recently,research studies have been conducted on the mechanism of strata deformation and stress distribution under PVCS mining [22–26].Huang et al.[26] used a numerical method to model the stress redistribution in the PTCS during the excavation of an upper PVCS.The results showed that the increase in the mining height of the upper PVCS promotes pressure release in the PTCS while the hard rock floor of PVCS reduces the pressure release.Yang et al.[27] used 2D and 3D numerical models to investigate mining-induced stress and the extraction effect on pressure relief for large mining height with upward mining sequence.The results showed that high stiffness and thick strata in the main roof are important in delaying the caving time and decreasing void ratio of the caved rock.The compressive and tensile zones around the coal seam being mined are controlled by the number,thickness,and strength of the stiff and thick layers.In order to control coal and gas outbursts in a pressure-relieved boundary areas,Liu et al.[28] experimented with a mining method based on laboratory experiments and field measurements in Qingdong mine.The results showed that coal mass permeability has a negative exponential relationship with the vertical stress in the pressure-relief boundary areas.Additionally,coal permeability in these areas is significantly different with that around the center of the pressure-relief area.Coal and gas outburst accidents are likely to be triggered in these areas.Effective precautions against coal and gas outburst accidents can be implemented based on stress re-distribution characteristics.Cao et al.[29] obtained the protection range of the lower PVCS in a thin coal seam through FLAC3D simulation in order to promote gas extraction.Cheng et al.[30] conducted theoretical analyses and physical experiments to study the characteristics of the pressure-relief gas source and accumulation zone,and the pressure relief gas extraction in the PTCS by mining PVCS with multiple weak inter-burden layers.The results indicated that drilling drainage boreholes in the gas accumulation zone (in combination with omni-directional stereo pressure-relief gas extraction technology) can improve gas drainage.Yuan et al.[31] studied permeability increase and pressure relief for Auger Mining (AM) through numerical simulation.They developed a permeability law as a function of stress,which showed that the intermediate coal pillar (ICP) width and the AM height and length are the main factors controlling permeability evolution.These studies suggested that the PVCS analyses must consider dynamic changes in stress redistribution for gas drainage.

Rock mass permeability is influenced by the two important factors- rock mass fracturing and stress redistribution.Stress relief and mining-induced fractures can increase the permeability of coal and rock mass.Therefore,identification of the stress relief zones in the upper PTCS above a lower PVCS is an important objective of this paper.The stress relief area in PTCS depends on several parameters such as mining depth,panel width,face advance rate,gob loading behavior and mining height of the PTCS,the thickness and inter-burden strata properties between the PVCS and PTCS.All these factors should be considered to analyze the effect of mining the PVCS.However,in this study three main factors- mining depth,mining height of the PVCS and inter-burden strata properties were chosen since their effects on the stress relief of PTCS during the excavation of PVCS are not well documented [32].“Stress relief coefficient” and “stress relief angle” are used to evaluate the stress relief influence of the PVCS excavation on PTCS.Through these two indicators,the stress relief effect can be quantified.Three-level and fourfactor orthogonal experimental design was used to study the effect of selected variables.Thirteen FLAC3D models were run for the three selected variables.Three typical inter-burden lithologies between PTCS and PVCS were analyzed-sandy mudstone,fine-grained sandstone and limestone.

3.1 Stress Relief Effects of PVCS Mining

The concept of gas drainage before mining and PVCS mining has been used extensively to eliminate coal and gas outburst risk [2].The PVCS is mined out before the excavation of the PTCS coal seam with coal and gas outburst risk.The PTCS could be located above or below the PVCS.In this study,#8 coal seam is the lower PVCS.The relief area represents the area where the coal and gas outburst risks for PTCS are reduced.Relief anglesαandβare the major parameters used to describe the relief effect after the excavation of a PVCS (Fig.3).These are expressed asα=arctan(b2/a2)andβ=arctan(b1/a1)[13].The stress relief zone of the PTCS and the surrounding rock mass represents the area where the coal and gas outburst risk is reduced and/or eliminated.Theb1andb2are the shortest distances between the two coal seams.The termsa1anda2are determined by the stress relief coefficient,which is defined by Eq.(1) below:

whereσis the vertical stress in the target PTCS during the face advance of the PVCS (MPa),σ0is the initial stress before the excavation (MPa),andis stress relief coefficient.Whenσ0<σ,then0,which means that pressure concentration occurs in the target PTCS after stress redistribution.Whenσ0>σ,then0,and the pressure relief occurs in the target coal seam.Whenσ0=σ,then=0,and there is no stress change in the target coal seam and the PTCS is not affected by the mining activity.

The pressure relief area is controlled by geology,gas pressure,deformations associated with mining activity,and the vertical stress in PTCS [13].The outburst-risk in the PTCS is reduced when stress relief occurs in the PTCS [33].Therefore,this paper adopts stress relief as the protection criteria.The critical stress relief value is given by the following formula,as given in Eq.(2).

whereσzcis the vertical stress for the coal seam;α1is the coal seam angle;λis the lateral pressure coefficient;γis the bulk density;andHis the initial depth of the PTCS coal seam.For the case study mine,γ=25kN/m3,H=470 m,andα1=0,then |σzc|≤2500×470=11.7 MPa.Therefore,when the stress decreases to 11.7 MPa,the gas outburst should not occur and it can be taken as the critical stress value.Substituting initial vertical stress of #3 coal seam of 13 MPa into Eq.(1),the critical value of stress relief coefficient of 0.1((13−11.7)/13=0.1)is obtained [13].

Figure 3:Illustration of the concept of protective coal seam mining effects in the upper and lower coal seams

Selection of the PVCS should meet two conditions:(1) Its excavation should not destroy the ability to mine PTCS (which means the PTCS should be located outside of the caved zone of the PVCS);and (2) It should achieve the best pressure relief effect.Upon excavation of the lower PVCS,the overlying strata may be divided into three relatively distinct zones based on fracture development:the caved zone,the fractured zone,and the continuous deformation zone [34,35].The height of these three zones [34] can be expressed below in terms of mining height:

HereHmis the height of the caved zone,Hdis the height of the fractured zone,and M is the mining height.These formulas were developed based on statistical analysis of a large field database [34].These are applicable where mining height is less 3 m and overburden strata are of medium strength (like sandstone,argillaceous limestone,sandy shale,mudstone,and sandy mudstone).The different mining heights of the PVCS will affect the stress relief in the PTCS coal seams (Tab.3).

3.2 Stress Relief and Permeability

Coal and rock mass have naturally-existing geological discontinuities (joints,cracks,and faults).Gas is stored in the pores and fractures mainly in an adsorbed state.The mining of the PVCS increases coal deformations and redistributes stress around adjacent coal and rock mass with stress relief or concentration.The permeability of coal and rock mass in some areas may increase through additional bedding plane fractures and cracks.The gas in both the PTCS and PVCS may desorb and flow into high-permeability regions and the newly formed fracture network [36].Wang et al.[37] investigated changes in the permeability of sedimentary rocks (mudstone,medium-sandstone etc.) using the confined compression test and related it to their complete stress-strain history (Fig.4).The results showed that the permeability varied with the strain states.Prior to reaching the peak strength,the permeability decreased with increased loading.However,permeability increased beyond the peak strength during the strain–softening behavior.

Table 3:Pressure relief effect in #3 coal seam for different mining heights of #8 coal seam

Figure 4:Permeability and complete stress-strain curves for selected rock types under confined compression [37] (a) mudstone (b) medium-grain sandstone

Li et al.[20] studied permeability changes using a stress-based permeability model.The established model for both the fracture and matrix systems was given in Eq.(5):

whereKis the bulk modulus of the coal-fracture assemblage mass;Δpf,Δpmare the pore pressure differences in the fracture and coal matrix;Δσis the stress difference;αis the Biot’s coefficient;kf0,km0are respectively the initial permeability values of the fracture and coal matrix;Φf0,Φm0are the initial porosities of the fracture and coal matrix;andεL,pLare Langmuir-type constants.The∊Lis the sorption-induced strain constant,which represents the volumetric strain when the pore pressure approaches infinity;Rmis the modulus reduction ratio,ranging between 0 and 1;pLis the pore pressure at which the measured sorption-induced strain is equal to half ofεL.

The gas migration rate is related to coal permeability,which is directly determined by the opening and closure of fractures.Therefore,the relationship between stress and permeability can be studied indirectly by studying the relationship between changes in stress,cracks formation and deformation,and rock mass permeability.Fig.5 illustrates relationship between crack orientation and applied stress.The relationship between normal stress and crack opening can be expressed as follows (Eqs.(6)–(9)) [24,38–41]:

whereσnis the normal stress,kn0is the initial stiffness of the fracture,δis the degree of fracture closure,andδmis the maximum degree of fracture closure.

Figure 5:Crack deformation as normal stress decreases (the symbol ↓suggests that the stress decreases [24].(a) Horiontal crack (b) vertical crack

The degree of fracture closure is expressed as follows:

The permeability of a single crack is expressed as:

In the above equation,b=b0−δ,when the crack is completely closed,b0=δm,andb=δm−δ.

Permeability can then be expressed in a dimensionless form:

whereσn0=kn0δmis the initial normal stress;kfis the dimensionless permeability;σn/σn0is the ratio of the normal stress to its initial value.The relationship betweenkfandσn/σn0in Eq.(9)implies that the permeability of coal and rock mass increases with the decrease in normal stress.Thus,the influence of mining PVCS on the permeability changes in the PTCS can be related to changes in stress and fracture permeability.This paper has attempted to analyze stress relief zones as the face advances in the PVCS mining and sensitivity analyses of selected variables using an orthogonal experimental design.The results should assist in designing the drainage borehole layouts in the PTCS.

4 Numerical Model Development

4.1 Description of the Numerical Model

Numerical models were developed to simulate rock mass behavior during the mining of the lower PVCS at the study mine [42].Z-coordinate represents the vertical direction of the model,and X-and Y-coordinates are the two horizontal directions.The boundary conditions for the model were applied:the bottom horizontal boundary is constrained in the Z-direction while the top boundary is free to deform;and the surrounding displacement boundaries are horizontally constrained.The dimensions of the model are 200 m×300 m×96 m (Fig.6).Element size increases with increasing distance from the areas of interest where the element size is 1 m×1 m×1 m.The #8 PVCS coal seam is mined for #3 PTCS coal seam,and the range(along Z direction) of these two coal seams are 31–34 m and 70–76 m from the bottom of the model.There are six lithologies with different engineering properties (Tab.4) between the two coal seams.The model overburden consists of 500 m of rock and topsoil,and it is simulated by applying a constant vertical stress (12.5 MPa) on the upper boundary of the numerical model,and the other strata with a thickness of 20 m (0.5 MPa) above the PTCS is included in the model.The pre-mining horizontal stresses included in the model are based on the data measured in the field.Their values in the Y and X directions were approximately 13.26 and 14.56 MPa [4].A linear gob loading model behind the face was incorporated in the numerical models based on our previous research [13].The “gob stress recovery distance” (the cover pressure re-establishment distance) from the working face is estimated to be 120 m,which is named as Gob1 in this paper.The pressure relief angle for the original model is calculated when the pressure relief coefficient is larger than 0.1 [13].All monitoring points for analysis within the model are located at the #3 coal seam level along the mining direction.The longwall face movement is simulated with an advance of 2 m for each cut.Tab.4 lists the mechanical properties of the rock mass for each layer were developed based on GSI (Geological Strength Index) and Hoek–Brown failure criterion [43].

The generalized Hoek–Brown criterion (Eqs.(10)–(16)) is expressed as [43]:

whereσciis the uniaxial compressive strength (UCS) of the intact rock,σ′1andσ′3are the major and minor effective principal stresses at failure,andmb,sandaare material constants.mbis a reduced value of the material constantmiand is given by:

where s and a are constants for the rock mass,and we have the following relationships:

whereDis a factor depending upon the degree of disturbance to which the rock mass has been subjected by blast damage and stress relaxation.It varies from zero for undisturbed subsurface rock mass to one for very disturbed rock masses.D=0 is assumed in our study.

Figure 6:FLAC3D numerical model and simulated lithologies [13]

Table 4:Hoek–Brown parameters,GSI and Rock mass properties used in numerical modeling

The estimation of the rock mass strength is:

4.2 Data Analyses Methodology

Two panels of the PVCS (#8 coal seam) are simulated in the FLAC3D model,but only 84306 longwall face is excavated in this study.As shown in Fig.7,the length and width of the longwall face are 160 and 270 m,respectively.The location of headgate and tailgate entries are shown in Fig.7a,and a coal pillar with 20 m width is setup between these two gate entries.To improve the calculation speed of the model,only half the width of the longwall face is modeled.A set up room with 7 m width is excavated first from the starting line (20 m from the edge of model).Shields are setup in this area and are advanced with the coal face advance.The face is advanced in increments of 2 m along the Y direction.The stress relief in different locations of the PTCS is analyzed using 12 monitoring lines along the Y-directions at the Z coordinate of 73,which is set up in the #3 coal seam.

Figure 7:Elevation and plan views of modeled faces.(a) Plan view (b) Elevation view

Fig.8 shows that the cross section with Z coordinate equals to 73 cuts through the FLAC3D model in the middle of #3 coal seam (Fig.8a).The distance from the left model edge to the coal pillar is 140 m.Therefore,12 monitoring lines were established along the cutting plane from Y=0 to Y=200 with total of 1000 points selected for analysis with the interval between adjacent monitoring lines and two adjacent points being 10 and 0.2 m,respectively (Fig.8b).The vertical stress change for different face advance values (lower #8 coal seam) is collected for every single point along the monitoring lines.The data above are used to calculate the stress relief coefficient while the PVCS face advances.

Figure 8:Data monitoring lines layout (a) cross section in FLAC3D model (b) the monitoring lines location

5 Results and Discussion

5.1 Stress Redistribution and Stress Relief Coefficients

Vertical stress redistribution and pressure relief coefficient in the PTCS are shown in Figs.9 and 10 for different face advance positions in #8 coal seam for (60,0,73)–(60,200,73)monitoring lines.

Figure 9:Vertical stress distribution for different face advance values in the PTCS

Figure 10:Pressure relief coefficient for different face advance values in the PTCS

1) For 30 m face advance in the PVCS,the minimum vertical stress in the PTCS is located about 18 m behind the coal face.The stress distribution curve in the PTCS is approximately V-shaped,and the maximum pressure relief area is located around the middle of the gob.The vertical stress is about 10.5 MPa,and the pressure relief coefficient is about 0.19.The vertical stress increased ahead of the longwall face and behind in the set up room areas (Figs.9 and 11).The peak stress values in these regions are 14.5 and 14.9 MPa,respectively,with stress concentration factor values (SCF) of 1.1 and 1.14.Since the immediate roof in the PVCS has not caved,the pressure relief and stress concentrations in the PTCS are relatively small.

Figure 11:Vertical stress contours along the face advance direction in the PTCS for 30 m PVCS face advance

2) For 50 m face advance in the PVCS (Fig.12),the minimum vertical stress in the PTCS is about 27 m behind the longwall face with the vertical stress values of about 9.5 MPa,and the stress relief coefficient of 0.28 (Figs.9 and 12).The V-shape depression extends over a larger area as compared to 30 m face advance.The maximum pressure relief zone is located around the center of the gob about 10 to 30 m behind the working face.The PTCS has lower vertical stresses with a stress relief coefficient of more than 0.15.The peak stress values in the PTCS are 15 and 16.1 MPa around the front of the longwall face and behind the open-off cut area of the PVCS with the SCF values of 1.15 and 1.24.These values are higher than for 30 m face advance.The main roof collapses,and part of the overburden weight is transferred to the coal face and the open mined-out area.

3) For 100 m face advance in the PVCS,the vertical stress contour in the PTCS is shown in Fig.13.In conjunction with Figs.9 and 10,it should be noted that the minimum vertical stress in the PTCS is about 35 m behind the coal face.The vertical stress value in the PTCS is about 9.8 MPa,and the stress relief coefficient is about 0.26 similar to face advance of 50 m.The vertical stress in the protected layer is generally U-shaped,and the maximum pressure relief in the protected layer is located on the gob side (Fig.10).The vertical stress in the PTCS significantly reduces around 15 to 75 m behind the coal face,and the stress relief coefficient increases up to about 0.15.The vertical stress distribution in the PTCS is asymmetrical,and the degree of pressure relief is slightly lower around the peak value area,indicating that the caving rock mass in the gob is gradually compacted with face advance,and its ability to carry more load.The peak stress value in the PTCS both in the face area and mined-out area increases further,reaching values of 16 and 17 MPa with the SCF values of 1.23 and 1.31.The pressure relief angle is calculated using the critical values of stress relief coefficients (Section 3.1).The pressure relief angles around the face and around the back of the mined-out areas are 75.1◦and 53.7◦(Fig.13).

Figure 12:Vertical stress contours along face advance in the PTCS for 50 m PVCS face advance

Figure 13:Vertical stress contour along face advance direction of the PTCS for 100 m face advance of the PVCS

In summary,the degree and extent of pressure relief in the PTCS increases gradually with the face advance in the PVCS.With collapse of the immediate roof and main roof behind the face,the overburden stress is transferred ahead of the face and on to the caved gob.With the increase in roof deformation and its vertical movement,the broken rock mass in the gob is gradually compacted,which gradually increases the supporting effect for the main roof,thus reducing to some degree the stress relief of the PVCS in the middle and rear of the gob.Assuming the central axis of the gob as the center line,stress relief around the front of the face and back in the gob area are asymmetrical.Most of pressure relief is concentrated ahead of the face.

The above analyses show that the stress relief degree and extent in the PTCS changes with different face advance distance of the PVCS.The stress relief zone in PTCS for 100 m face advance of PVCS (cutting plane Z-coordinate=73) is shown in Fig.14.To quantify the stress relief,the vertical stress distribution and the pressure relief coefficient for different locations of PTCS along the x-direction,stress data were collected and analyzed for twelve monitoring lines(Figs.8 and 9) in the PTCS (Z-coordinate=73) for 100 m face advance of PVCS,which is 130,120,110,...and 20 m away from the coal pillar.

Figure 14:Stress relief zones in the PTCS for 100 m PVCS face advance for Z-coordinate=73

For 100m face advance of PVCS,the vertical stress in the PTCS (Figs.15 and 16) is gradually reduced,and the pressure relief coefficient increases over a distance from 130 m (X-coordinate=10) away from the coal pillar to 30 m (X-coordinate=110) near the coal pillar.This is true for the pressure relief area above the middle of the PVCS gob or the high-pressure area ahead of the coal face and behind the mined-out area in the gob.

Figure 15:Vertical stress values along x-direction in the PTCS for 100 m face advance in PVCS

Figure 16:Pressure relief coefficient for different locations in the PTCS along x-direction for 100 m face advance

When the distance from the coal pillar is 130 m which is close to the central axis(X-coordinate=0) of the working face,the vertical stress behind the mined-out area and coal face area are 16.8 and 15.8 MPa with SCF values of 1.29 and 1.22,and the maximum pressure relief coefficient of about 0.22.

When the distance from the coal pillar is 30 m (blue lines in Figs.15 and 16),the vertical stress in the mined-out and the coal face areas are 15.5 and 14.2 MPa with SCF values of 1.19 and 1.09 with the maximum pressure relief coefficient of about 0.27.

The maximum pressure relief coefficient in the PTCS is similar for different distances away from the coal pillar.However,as the monitoring lines get farther and farther from the coal pillar,the location of the maximum pressure relief coefficient gets closer to the coal face.Therefore,the PTCS goes through several stages along the direction of face advance:Compression,stress relief,and expansion,increased stress relief and expansion,stabilized stress relief and expansion,decrease in stress relief and expansion,compression,and stable stress relief and expansion.

Fig.16 also shows that the pressure relief angle variation along the direction of the face advance for different distances away from the coal pillar.Tab.5 includes the pressure relief angle for the pressure relief coefficient of 0.1.Overall,the pressure relief angle in the mined-out area is generally smaller than that ahead of the coal face.As the distance to the pillar decreases,the pressure relief angle on both sides increases first and then decreases,and the maximum value is 67◦in the mined-out area and 80◦ahead of the coal face side.

Table 5:Relief angle in the PTCS along X-direction

In order to analyze the vertical stress changes along the strike of the PTCS for different positions in the gob during the mining of the PVCS,four cross-section planes were analyzed perpendicular to the Y-axis for different positions behind the face in the gob during the PVCS mining:100 m (Fig.17a),60 m (Fig.17b),20 m (Fig.17c),and−20 m (Fig.17d).The pressure relief extent was also analyzed when the PVCS is advanced 120 m.With the increase in distance to the working face,the extent and degree of stress relief gradually extend and reach the maximum at about 40 m behind the face.If the distance from the working face decreases further,both the range and the degree of stress relief gradually decrease.

Figure 17:Vertical stress contour along the strike of the PTCS for 120 m face advance of the PVCS 100 m behind the coal face (b) 60 m behind the coal face (c) 20 m behind the coal face(d) 20 m ahead of the coal face

5.2 Effects of Mining Height on Stress Relief

Fig.18 shows stress relief coefficient in the PTCS for different mining heights as the working face in the PVCS advances.These data are for the data monitoring line (60,0,72)–(60,200,72).

For 30 m face advance,the stress relief extent and the coefficient are almost the same.The maximum values of the stress relief coefficient for 1 and 3 m mining heights are 0.12 and 0.14.The location center of the relieved area is about 17 m behind the working face.

For 50 m face advance,the above coefficients are 0.21 and 0.28.Similar data for 70 and 120 m face advances vary 0.23 to 0.24.The above results show that the influence of the mining height on the extent and the degree of stress relief in the PTCS vary only in the initial stage of mining.

Additional analyses were performed to quantify the influence of the mining height on the stress relief coefficients and stress relief angles along the strike direction (Tab.6).The stress relief angle for two mining heights decrease with values of 82◦for 30 m face advance to 41.3◦for 120 m face advance on the mined-out area side.The angle for 3 m mining height is slightly larger than for 1 m mining height.The difference decreases with face advance and are almost the same for face advance of 120 m.The values for 3 m mining height are,however,always larger than for 1m mining height on the coal face side.

Figure 18:Variation of the stress relief coefficient in the PTCS for different mining heights for different face positions

Table 6:Relief angles along strike (Y-coordinate direction) for different mining heights for different face advance

5.3 Effect of the Inter-Burden Lithology between the PVCS and PTCS on Stress Relief

The rock mass properties within the interburden layers affect the stress relief in the PTCS.Rock mass properties of interest include elastic moduli,cohesion,and bulk density.The general consensus is that the disturbance and damage within the overlying strata will increase with the decrease in strength of the inter-burden layers,and that should increase the extent and degree of stress relief.To quantify these effects,three interburden lithologies with engineering properties in Tab.4 were simulated.

Fig.19 shows the stress relief coefficients for different interburden lithologies.They decrease with the higher stiffness or elastic modulus of the interburden layers (limestone

Figure 19:Pressure relief coefficients in the PTCS with different properties of inter-burden rock mass for different face advances of the PVCS (a) 30 m (b) 70 m (c) 120 m (d) 150 m

5.4 Sensitivity Analysis of Selected Variables

Nine numerical models were run based on the three-level and four-factor orthogonal experiment design without a consideration of interaction among those three factors (Table L9,34) [44].The models and results are shown in Tab.7.

Table 7:Parameters of FLAC3D models for sensitivity analyses

The F value reflects the degree of influence of each factor on the experimental results.It is obtained by the ratio of the squared sum of the mean deviations of the influencing factors to the sum of the squared deviations of the error [44].Tab.8 shows the analysis of variance results.The table shows that the order of importance for the selected variables is mining depth,inter-burden lithology,and mining height.

Table 8:Analysis of variance results for selected variables

6 Concluding Remarks

Mining a PVCS has been widely used to reduce the coal and gas outburst risk where multiple coal seams are present near.This paper has numerically analyzed for a case study mine the progressive effects of this practice on the extent (stress relief angle) and the degree of stress relief as the face in a lower PVCS seam is advanced.The relative importance of mining depth,inter-burden lithology between the two seams,and mining height on the above variables are also analyzed.

(1) With face advance in the PVCS,each point in the protected seam (PTCS) undergoes thru varying stress redistribution in different stages of mining that can be beneficially used from gas drainage point of view.The 3-D shapes of these zones outlined in the results section can aid gas drainage boreholes design.

(2) Both the extent of stress relief and the coefficient of stress relief in the PTCS vary ahead of and behind the face in the mined-out areas.These are affected by the mining face position in relation to caving behind the face in the PVCS.

(3) For this case study,the stress relief angle behind the face in the mined-out area is smaller than ahead of the face.The pressure relief angle along the gob side and the face side increases first and then decreases with face advance to an equilibrium value.The maximum relief angle value is 67◦on the gob side and 80.8◦on the coal face side.

(4) Mining height affects the stress relief angle during the early stages of mining on the gob side.Ahead of the face,the stress relief angle increases with the mining height.

(5) The pressure relief in the PTCS is affected most by the mining depth,followed by the inter-burden lithology engineering properties and the mining height.The elastic modulus of the inter-burden layers has the most effect on the stress relief in the PTCS among considered engineering properties.

Acknowledgement:The authors would also like to express special thanks to the editor and anonymous reviewers for their professional and constructive suggestion.

Funding Statement:This paper was supported by the Natural Science Foundation of Jiangsu Higher Education Institutions (No.20KJB440002),the National Natural Science Foundation of China (Project Nos.51804129,51808246 and 51904112),China Postdoctoral Science Foundation(No.2020M671301),the Postdoctoral Science Foundation of Jiangsu Province (Nos.2019K139 and 2019Z107),the Huai’an Science and Technology Plan project (No.HAB201836),the Industry Education Research Cooperation Projects in Jiangsu Province (No.BY2020007),Undergraduate Innovation and Entrepreneurship Training Program (No.202011049111XJ) and the Foundation of Huaiyin Institute of Technology (No.Z301B20530).

Conflicts of Interest:The authors declare that they have no conflicts of interest in reporting the present study.