# Numerical Simulation of the Volumetric Strain Distribution of a Protected Coal Seam

Authors: Hengyi Jia
Journal: Geofluids
 (2021)

---

## Abstract

To investigate the deformation characteristics of protected coal seams, the numerical simulation of the mining of an upper protective coal seam was carried out in the present study. Based on the basic definition of strain, a method for the extraction of the strain data of the protected coal seam was proposed, and the strain distribution characteristics were obtained. It was found that thex-direction strain is mainly distributed near the coal pillars on both sides and inside the goaf, the y-direction strain is mainly distributed at the working face, the initial mining line, and inside the goaf, and the z-direction strain is mainly distributed at the working face, the initial mining line, the coal pillars on both sides, and inside the goaf. The distribution characteristics and the value of volumetric strain were found to be basically consistent with the z-direction strain. As the working face advances, the protected coal seam undergoes compression and damage expansion in turn. The turning point between compression and damage expansion is approximately 15 m in front of the working face. The variation law of gas drainage in the boreholes of the protected coal seam is closely related to the distribution characteristics of volumetric strain. The results of this research are of great significance for the comprehensive investigation of the effects of pressure relief and the permeability enhancement of protective coal seam mining.


---

## Body

## 1. Introduction

Coal and gas outburst accidents have long caused substantial casualties and property losses to coal enterprises in China [1, 2]. With the increase of the mining depth, both the in situ stress and gas pressure increase, and the risk of coal and gas outburst is further enhanced [3, 4]. Both theoretical research and production practice demonstrate that the mining of protective coal seams can significantly reduce the in situ stress, expand cracks, discharge gas, and ultimately eliminate the outburst risk of protected coal seams [5–9].The elimination of the outburst risk of protected coal seams depends on the degree of pressure relief after mining. Scholars have carried out extensive research on the deformation laws of protected coal seams from the perspectives of the relative distance between protective and protected coal seams [10, 11], double protective coal seam mining [12], the mining depth [13], and the impacts of protective coal seam mining on the roadways of protected coal seams [14]. However, while these studies investigated the deformation of protected coal seams in the vertical direction, the deformation of protected coal seams should be reflected in three directions (the vertical, mining, and tilt directions). Shi and Liu [15] obtained the deformation laws of protected coal seams in both the vertical and mining directions based on the theory of rock fracture and damage, they did not obtain the deformation law along the tilt direction. He et al. [16] proposed a “four fixed point enclosed area” analysis method for the deformation of protected coal seams based on similar simulation experiments and obtained the volumetric strain distribution characteristics. However, it is time-consuming, laborious, and inconvenient to change the simulation conditions to carry out similar simulation experiments. In contrast, in the numerical simulation method, the simulation conditions can easily be changed by modifying the model parameters. Therefore, obtaining the distribution characteristics of the volumetric strain via numerical simulation is important for the comprehensive investigation of protected coal seam deformation.In this research, FLAC3D numerical simulation software was employed to study the pressure relief effect of a protected coal seam. According to the definition of strain, a method for extracting the strain data of the protected coal seam was proposed, and the strain distribution characteristics were obtained.


## 2. Numerical Simulation of Protective Coal Seam Mining

### 2.1. General Situation of the Coal Mine under Study

The Zhongxing coal mine is located in Shanxi Province, China. The main coal seam of the coal mine is the #2 coal seam, which has low permeability and a severe risk of coal and gas outburst. The locally mineable #02 coal seam is located about 10 m above the #2 coal seam and has a lower risk of outburst and a reduced gas content as compared to the #2 coal seam. To eliminate the outburst risk of the #2 coal seam, the #02 coal seam was mined first as an upper protective coal seam. Because the #02 coal seam is locally mineable, the size of its working face is less than that of the protected coal seam. The sizes of the working faces of the protective and protected coal seam are, respectively,1160×75m and 1595×190m (length×width). The layouts of the working faces are presented in Figure 1, and their basic parameters are reported in Table 1.Figure 1
Working faces layouts.Table 1
Basic parameters of the coal seams.
ParameterProtective coal seamProtected coal seamThickness of the coal seam0.51 m (in average)2 m (in average)Tilt angle of the coal seam7° (in average)7° (in average)Original gas pressure—0.82 MPaOriginal gas content3.37 m3/t10.4 m3/t


### 2.2. Establishment of the Numerical Model

Based on the strata conditions of the Zhongxing coal mine, the numerical model shown in Figure2 was established by FLAC3D numerical simulation software. The strata conditions and mechanical parameters of the coal mine are reported in Table 2. Mining a coal seam will cause the deformation of the rock mass and a reduction in rock strength. Thus, the rock mass can be regarded as elastic-plastic material that exhibits strain-softening characteristics after it is damaged [17]. Therefore, the Mohr-Coulomb strain-softening model was selected, and the tilt angle of the rock strata in the model was set to zero for simplicity.Figure 2
Numerical calculation model.Table 2
Mechanical parameters of coal and rock.
LithologyDensity/kg/m3Bulk modulus/GPaShear modulus/GPaFriction angle/°Cohesion/MPaDilatancy angle/°Uniaxial tensile strength/MPaSandy mudstone27403.73.1352.6102.4Mudstone26902.01.5311.5101.0Fine sandstone27403.52.2302.0102.0Coal14501.00.8200.6100.5Limestone30803.93.5434.8103.2According to the stress condition of the Zhongxing coal mine, compressive stress of 20 MPa was imposed at the top of the model, and a fixed boundary was imposed at the bottom of the model. The other sides of the model were set as rolling boundaries with an imposed stress of 20 MPa. The length, width, and height of the model were, respectively, 300, 300, and 131 m. The mining length of the working face was 200 m, which extends from 50 to 250 m along they-direction. Moreover, the mining width was 100 m, which extends from 100 to 200 m along the x-direction, and the mining height was 1 m. The model contained 349,200 grids and 400,243 zones (as shown in Figure 2).


### 2.3. Strain Data Extraction Method

Under the influence of the mining stress, a protected coal seam will be compressed or expand in three directions (the vertical, mining, and tilt directions). In other words, the deformation is mainly reflected in the change of the coal seam volume. The volumetric strain of a coal seam is equal to the sum of its strains in the three directions, and the extraction method of the strain data is presented as follows.The strain in thez-direction (vertical direction) is
(1)ηz=dfloor‐droofm,where ηz is the strain of the protected coal seam in the z-direction, dfloor is the vertical displacement of the coal seam floor, droof is the vertical displacement of the coal seam roof, and m is the initial thickness of the coal seam.The length and width of the grids in the #2 coal seam were both 2.5 m. When extracting the vertical strain, the vertical displacement of the grid boundary line between2≤x≤3 along the coal seam floor was extracted first. The distance was then increased by 10 m along the x-direction, and the vertical displacement of the grid boundary line between 12≤x≤13 was extracted. By analogy, a vertical displacement matrix of the coal seam floor with 120rows×30columns was obtained. Via the same method, the vertical displacement matrix of the coal seam roof was also obtained. According to Eq. (1), the vertical strain matrix of the #2 coal seam was calculated.The strain in thex-direction (tilt direction) is
(2)ηx=dxi‐dxjlxi−xj,where ηx is the strain of the protected coal seam in the x-direction, dxi is the displacement of the coal seam in the x-direction at position xi, dxj is the displacement of the coal seam in the x-direction at position xj, and lxi∼xj is the initial width from xi to xj.When extracting the strain in thex-direction, the displacement in the x-direction of the grid boundary line between 2≤x≤3 was first obtained, which corresponds to the displacement dxi. Then, the displacement of the grid boundary line between 4.5≤x≤5.5 was obtained, which corresponds to the displacement dxj. The initial width (2.5 m) between the two grid boundary lines corresponds to the lxi−xj. According to Eq. (2), the x-direction strain between the two grid boundary lines can be calculated. The distance was then increased by 10 m along the x-direction, and the strain of the grid boundary lines between 12≤x≤13 and that between 14.5≤x≤15.5 were obtained. By analogy, a strain matrix of the coal seam in the x-direction with 120rows×30columns was obtained.The strain in they-direction (mining direction) is
(3)ηy=dyi‐dyjlyi∼yj,where ηy is the strain of the protected coal seam in the y-direction, dyi is the displacement of the coal seam in the y-direction at position yi, dyj is the displacement of the coal seam in the y-direction at position yj, and lyi∼yj is the initial length from yi to yj.The grid boundary line between2≤x≤3 was first determined, after which the displacements in the y-direction of the grids along the grid boundary line were extracted. For any grid, the side with a smaller y-coordinate corresponds to dyi, the side with a larger y-coordinate corresponds to dyj, and the initial length of the grid corresponds to lyi∼yj. The strain of each grid can be calculated by substituting the extracted data into Eq. (3). The distance was then increased by 10 m along the x-direction, and the strains of the grids along the grid boundary line between 12≤x≤13 were obtained by the same method. By analogy, a strain matrix of the coal seam in the y-direction with 120rows×30columns was obtained.The volumetric strain matrix of the protected coal seam can be obtained by directly adding the elements of the three strain matrices. The calculation formula is as follows:(4)ηv=ηx+ηy+ηz,where ηv is the volumetric strain of the protected coal seam.


## 2.1. General Situation of the Coal Mine under Study

The Zhongxing coal mine is located in Shanxi Province, China. The main coal seam of the coal mine is the #2 coal seam, which has low permeability and a severe risk of coal and gas outburst. The locally mineable #02 coal seam is located about 10 m above the #2 coal seam and has a lower risk of outburst and a reduced gas content as compared to the #2 coal seam. To eliminate the outburst risk of the #2 coal seam, the #02 coal seam was mined first as an upper protective coal seam. Because the #02 coal seam is locally mineable, the size of its working face is less than that of the protected coal seam. The sizes of the working faces of the protective and protected coal seam are, respectively,1160×75m and 1595×190m (length×width). The layouts of the working faces are presented in Figure 1, and their basic parameters are reported in Table 1.Figure 1
Working faces layouts.Table 1
Basic parameters of the coal seams.
ParameterProtective coal seamProtected coal seamThickness of the coal seam0.51 m (in average)2 m (in average)Tilt angle of the coal seam7° (in average)7° (in average)Original gas pressure—0.82 MPaOriginal gas content3.37 m3/t10.4 m3/t


## 2.2. Establishment of the Numerical Model

Based on the strata conditions of the Zhongxing coal mine, the numerical model shown in Figure2 was established by FLAC3D numerical simulation software. The strata conditions and mechanical parameters of the coal mine are reported in Table 2. Mining a coal seam will cause the deformation of the rock mass and a reduction in rock strength. Thus, the rock mass can be regarded as elastic-plastic material that exhibits strain-softening characteristics after it is damaged [17]. Therefore, the Mohr-Coulomb strain-softening model was selected, and the tilt angle of the rock strata in the model was set to zero for simplicity.Figure 2
Numerical calculation model.Table 2
Mechanical parameters of coal and rock.
LithologyDensity/kg/m3Bulk modulus/GPaShear modulus/GPaFriction angle/°Cohesion/MPaDilatancy angle/°Uniaxial tensile strength/MPaSandy mudstone27403.73.1352.6102.4Mudstone26902.01.5311.5101.0Fine sandstone27403.52.2302.0102.0Coal14501.00.8200.6100.5Limestone30803.93.5434.8103.2According to the stress condition of the Zhongxing coal mine, compressive stress of 20 MPa was imposed at the top of the model, and a fixed boundary was imposed at the bottom of the model. The other sides of the model were set as rolling boundaries with an imposed stress of 20 MPa. The length, width, and height of the model were, respectively, 300, 300, and 131 m. The mining length of the working face was 200 m, which extends from 50 to 250 m along they-direction. Moreover, the mining width was 100 m, which extends from 100 to 200 m along the x-direction, and the mining height was 1 m. The model contained 349,200 grids and 400,243 zones (as shown in Figure 2).


## 2.3. Strain Data Extraction Method

Under the influence of the mining stress, a protected coal seam will be compressed or expand in three directions (the vertical, mining, and tilt directions). In other words, the deformation is mainly reflected in the change of the coal seam volume. The volumetric strain of a coal seam is equal to the sum of its strains in the three directions, and the extraction method of the strain data is presented as follows.The strain in thez-direction (vertical direction) is
(1)ηz=dfloor‐droofm,where ηz is the strain of the protected coal seam in the z-direction, dfloor is the vertical displacement of the coal seam floor, droof is the vertical displacement of the coal seam roof, and m is the initial thickness of the coal seam.The length and width of the grids in the #2 coal seam were both 2.5 m. When extracting the vertical strain, the vertical displacement of the grid boundary line between2≤x≤3 along the coal seam floor was extracted first. The distance was then increased by 10 m along the x-direction, and the vertical displacement of the grid boundary line between 12≤x≤13 was extracted. By analogy, a vertical displacement matrix of the coal seam floor with 120rows×30columns was obtained. Via the same method, the vertical displacement matrix of the coal seam roof was also obtained. According to Eq. (1), the vertical strain matrix of the #2 coal seam was calculated.The strain in thex-direction (tilt direction) is
(2)ηx=dxi‐dxjlxi−xj,where ηx is the strain of the protected coal seam in the x-direction, dxi is the displacement of the coal seam in the x-direction at position xi, dxj is the displacement of the coal seam in the x-direction at position xj, and lxi∼xj is the initial width from xi to xj.When extracting the strain in thex-direction, the displacement in the x-direction of the grid boundary line between 2≤x≤3 was first obtained, which corresponds to the displacement dxi. Then, the displacement of the grid boundary line between 4.5≤x≤5.5 was obtained, which corresponds to the displacement dxj. The initial width (2.5 m) between the two grid boundary lines corresponds to the lxi−xj. According to Eq. (2), the x-direction strain between the two grid boundary lines can be calculated. The distance was then increased by 10 m along the x-direction, and the strain of the grid boundary lines between 12≤x≤13 and that between 14.5≤x≤15.5 were obtained. By analogy, a strain matrix of the coal seam in the x-direction with 120rows×30columns was obtained.The strain in they-direction (mining direction) is
(3)ηy=dyi‐dyjlyi∼yj,where ηy is the strain of the protected coal seam in the y-direction, dyi is the displacement of the coal seam in the y-direction at position yi, dyj is the displacement of the coal seam in the y-direction at position yj, and lyi∼yj is the initial length from yi to yj.The grid boundary line between2≤x≤3 was first determined, after which the displacements in the y-direction of the grids along the grid boundary line were extracted. For any grid, the side with a smaller y-coordinate corresponds to dyi, the side with a larger y-coordinate corresponds to dyj, and the initial length of the grid corresponds to lyi∼yj. The strain of each grid can be calculated by substituting the extracted data into Eq. (3). The distance was then increased by 10 m along the x-direction, and the strains of the grids along the grid boundary line between 12≤x≤13 were obtained by the same method. By analogy, a strain matrix of the coal seam in the y-direction with 120rows×30columns was obtained.The volumetric strain matrix of the protected coal seam can be obtained by directly adding the elements of the three strain matrices. The calculation formula is as follows:(4)ηv=ηx+ηy+ηz,where ηv is the volumetric strain of the protected coal seam.


## 3. Strain Distribution of the Protected Coal Seam

According to the method discussed in the previous section, the strain data of the protected coal seam were extracted. The data were processed by using the “surf” function and “meshgrid” command in MATLAB software, and the strain distributions were obtained. Figures3–7 present the strain distribution diagrams of the protected coal seam when the working face of the protective coal seam was mined to different distances. In the figures, positive values of strain indicate that the coal seam underwent compression, while negative values indicate that the coal seam underwent expansion.Figure 3
X-direction strain of protected coal seam.
(a)
Mined 50 m(b)
Mined 100 m(c)
Mined 150 m(d)
Mined 200 mFigure 4
Y-direction strain of protected coal seam.
(a)
Mined 50 m(b)
Mined 100 m(c)
Mined 150 m(d)
Mined 200 mFigure 5
Z-direction strain of protected coal seam.
(a)
Mined 50 m(b)
Mined 100 m(c)
Mined 150 m(d)
Mined 200 mFigure 6
Volumetric strain of protected coal seam.
(a)
Mined 50 m(b)
Mined 100 m(c)
Mined 150 m(d)
Mined 200 mFigure 7
The main view of the volumetric strain.Figure3 presents the x-direction strain of the protected coal seam, which was found to be mainly distributed near the coal pillars on both sides and inside the goaf, while the strain at the working face and the initial mining line was not obvious. The coal seam directly below the coal pillars was in a stress concentration zone, and the significantly increased vertical stress compressed the coal seam in the vertical direction and elongated it in the x-direction. The coal seam at the inner side of the coal pillars was located in a pressure relief zone, in which the vertical stress decreased the most and the coal seam was compressed in the x-direction under the increased horizontal deviatoric stress.Figure4 presents the y-direction strain of the protected coal seam, which was found to be mainly distributed at the working face, the initial mining line, and inside the goaf. The strain near the coal pillars on both sides was not obvious. The coal seam directly below the working face was found to expand in the y-direction, which was due to the vertical stress concentration in this area. The coal seam at the inner side of the initial mining line was compressed in the y-direction significantly; it was located in a pressure relief zone, and the degree of unloading of vertical stress in this zone was greater than that of horizontal stress. After the protective coal seam was mined, the floor strata tended to move to the goaf and stop in this zone. Under the dual action of the floor strata movement and vertical stress unloading, the coal seam in this zone was compressed in the y-direction.Figure5 exhibits the z-direction strain of the protected coal seam, which was found to be distributed at the working face, the initial mining line, the coal pillars on both sides, and inside the goaf. The coal seam directly below the working face, the initial mining line, and the coal pillars was compressed along the z-direction, whereas the coal seam in the goaf expanded significantly. In the goaf, the vertical strain was greater than the horizontal strain, which was determined by vertical stress unloading. With the advance of the working face, the stress in the central area of the goaf was gradually restored, and the degree of expansion of the coal seam exhibited a decreasing trend.Figure6 displays the volumetric strain of the protected coal seam. The distribution characteristics and the value of volumetric strain were found to be basically consistent with those of the z-direction strain. With the exploitation of the protective coal seam, the scope of the goaf expanded, and the volumetric strain presented an inverted saddle shape. To further investigate the distribution characteristics of the volumetric strain, the main view is shown in Figure 7. It can be seen from the figure that the volumetric strain experienced the stages of increasing, decreasing, and then increasing again from the boundary to the interior of the model. The volumetric strain reached the maximum value in the stress concentration zone and decreased to the minimum value near the working face. The turning point of volumetric strain from increasing to decreasing was found to be about 15 m in front of the working face, which indicates that this was the starting point of the protected coal seam changing from elastic compression to damage expansion.


## 4. Field Observation of the Gas Drainage Effect of the Protected Coal Seam

The changes in the permeability and volumetric strain of protected coal seam both reflect the degree of pressure relief of the coal seam; thus, the permeability and volumetric strain have internal unity. Assuming that the gas pressure of the protected coal seam remains constant during protective coal seam mining, the gas drainage quantity of the coal seam is proportional to the permeability. Therefore, the change law of the gas drainage quantity can be used to verify the correctness of the distribution characteristics of the volumetric strain.Two gas boreholes were preconstructed in the transportation roadway of the protected coal seam prior to the mining of the protective coal seam. The boreholes were arranged along the coal seam inclination direction and were used to measure the gas drainage quantity and gas drainage concentration with the advance of the working face (as shown in Figure1). The initial horizontal distances between boreholes #1 and #2 and the working face were 120 and 144 m, respectively. The length of the boreholes was 120 m, and the diameter was 94 mm.Figure8 presents the variation laws of the gas drainage quantity and gas concentration of the boreholes. The abscissa in the figure indicates the relative position between the boreholes and the working face. For example, x=100 means that the borehole was located 100 m in front of the working face, and x=−50 means that the borehole was located 50 m behind the working face. It can be seen from the figure that the mixed gas drainage quantity was relatively high at the initial stage of mining and decreased in a fluctuating manner with the advance of the working face. From about 15 m in front of the working face, the mixed gas drainage quantity, the pure gas drainage quantity, and the gas concentration began to increase rapidly. When the boreholes entered the goaf, these gas drainage values reached the maximum and remained stable for a period of time, after which they ultimately fluctuated inside the goaf. The variation laws of the gas drainage quantity and the gas concentration were found to be closely related to the distribution characteristics of the volumetric strain; the turning point of the rapid increase of gas drainage coincided with the starting point of damage expansion. These results prove the correctness of the distribution characteristics of the volumetric strain.Figure 8
Statistics of the gas drainage data [18].
(a)
#1 borehole(b)
#2 borehole


## 5. Discussion

Based on a numerical simulation, a method for the extraction of the strain data of a protected coal seam was proposed, and the strain distribution characteristics were obtained. Thex-direction strain was found to be mainly distributed near the coal pillars on both sides and inside the goaf, while the y-direction strain was found to be mainly distributed at the working face, the initial mining line, and inside the goaf. The deformation of the coal seam was affected by the change of stress and the deformable conditions, both of which caused the distribution difference of the horizontal strains. A stress concentration zone was formed around the stope, at which the significantly increased vertical stress caused the compression of the coal seam along the z-direction. The stress inside the goaf was greatly unloaded in all directions, and the coal seam swelled. The simulation results show that the distribution characteristics and the value of the volumetric strain were basically consistent with the z-direction strain, indicating that the z-direction strain has a decisive influence on the volumetric strain. Due to the restrictions of field measurement conditions and technology, scholars have often used vertical strain (z-direction strain in this paper) to express the pressure relief effect of protected coal seams; although this operation is not perfect in theory, it is acceptable in fact.With the advance of the working face, the protected coal seam entered the stress concentration zone from the original stress zone, and the coal seam experienced compression and damage expansion in turn. The turning point from compression to damage expansion was found to be about 15 m in front of the working face. With the further advance of the working face, the protected coal seam expanded rapidly after entering the goaf and then remained stable. Based on the results of gas drainage, it was determined that the mixed gas drainage quantity was higher in the initial stage and then gradually decreased with the advance of the working face, which was caused by the compression of the coal seam. The mixed gas drainage quantity, pure gas drainage quantity, and gas concentration began to increase rapidly at 15 m in front of the working face and then remained stable after entering the goaf. The variation law of gas drainage data was found to be closely related to the distribution characteristics of the volumetric strain, which proves the correctness of the numerical simulation results. Thus, it is feasible to verify the distribution characteristics of the volumetric strain of a protected coal seam via the gas drainage quantity and gas concentration in the boreholes.The volumetric strain of a protected coal seam reveals only the degree of pressure relief of the coal seam and does not reveal the degree of permeability increase. If an evolution model of coal permeability with volumetric strain as a variable is established, and if the model is combined with coal seam volumetric strain, the distribution characteristics of permeability can be obtained. Xie et al. [19] considered the volumetric change of damaged coal with cracks induced by the coupling of the abutment pressure under different mining conditions, pore pressure, and expansion caused by gas adsorption, and proposed a new term for mining-enhanced permeability that reflects the change in permeability by the volumetric change of coal. Cheng et al. [20] put forward a new model with strain as a variable for unloading fractured coal by considering the deformation effects of the effective stress and gas adsorption/desorption. These studies laid a theoretical foundation for obtaining the permeability enhancement effect of protected coal seams. However, there is a lack of research on the distribution characteristics of the volumetric strain of protected coal seams. Because coal and gas outburst accidents are becoming increasingly more serious [21–24], the present study both supplements and improves the previous research, and the results are of great significance to revealing the effect of the permeability enhancement of protected coal seams.


## 6. Conclusion

(1)
The deformation of coal seams is affected by both the change in stress and the deformation conditions. The vertical strain of a coal seam dominates the change law and the value of volumetric strain. Although it is not comprehensive to use the vertical strain to express the effect of pressure relief of protective coal seam mining in theory, it is acceptable in fact(2)
The change laws of the gas drainage quantity and gas concentration of a protected coal seam are closely related to the distribution characteristics of the volumetric strain. The turning point of the rapid increase of the gas drainage quantity and gas concentration coincides with the point of damage expansion. Thus, it is feasible to verify the distribution characteristics of the volumetric strain of a protected coal seam via the use of gas drainage data(3)
The volumetric strain only reveals the degree of pressure relief of a protected coal seam, but does not reveal the degree of permeability increase. The results of this study are of great significance for the comprehensive investigation of the effect of the permeability enhancement of protective coal seam mining



---

*Source: 3679862-2021-08-13.xml*