Abstract
The plastic zone range is an important parameter in the analysis of damage characteristics and the degree of damage to the rock surrounding a roadway. Based on the establishment of a plastic zone calculation model considering the creep effect, this paper obtains the characteristics of the change in the plastic zone damage range with time by solving the model. Additionally, the validity of the model is verified by field experiments. The research results can provide guidance for gas pressure measurement and gas drainage in coal mines.
Similar content being viewed by others
Introduction
Gas disasters constitute the primary threat to safe production in many coal mines1,2. On the one hand, determining the damage range of the plastic zone of the roadway can improve the accuracy of measuring the gas pressure (an important index in predicting the risk of coal and gas outbursts) in a coal mine; on the other hand, it can improve the sealing effect of the gas drainage borehole, thereby improving the gas drainage rate of the coal seam3,4,5,6. Therefore, clearly establishing the damage range of the plastic zone of the roadway can indirectly prevent the occurrence of coal mine gas disasters.
Many scholars have carried out research on the plastic zone of a roadway. Galin7 considered that the shape of the plastic zone around the roadway in the uniform stress field is circular, and the formula for calculating the circular plastic zone was derived. Additionally, Chen et al.8 established the boundary implicit equation for the plastic zone of a circular roadway. Abdel et al.9, Leitman and Villaggio10 considered that the shape of the plastic zone of roadways located in non-uniform stress fields should be elliptical. Guo et al.11 and Wang et al.12 pointed out that the shape of the plastic zone of a roadway could be circular, elliptical or butterfly-shaped under different in situ stress conditions. Zhou and Wu13, Yang et al.14 used a numerical simulation method to study the plastic zone range of a roadway under different stress conditions.
The practice of mining shows that at constant stress, the creep of coal and rock will occur as time goes on; that is, the deformation of coal and rock will continue to grow with time15,16,17,18,19. Even under low stress conditions, the creep of coal and rock would also occur. Due to the creep effect, the plastic zone of a roadway will increase with time and finally approach a certain value. It can be seen that the previous research on the plastic zone of a roadway did not consider the creep effect of the surrounding rock of the roadway and thus cannot better guide engineering practice. In this paper, a plastic zone calculation model considering the creep effect was established to analyze the variation in the plastic zone damage range with time. The validity of the model is also verified by field experiments.
Calculation model of the plastic zone of a roadway considering the creep effect
Model assumptions
To obtain the specific range of the plastic zone of a roadway, some necessary practical assumptions are made to simplify the solution. The basic assumptions are as follows:
-
(1)
Since the axial length of the roadway is much larger than the radial width, the stress problem of the roadway section can be regarded as a plane strain problem.
-
(2)
The stress around the roadway is assumed to be uniform and the shape of the plastic zone of the roadway is assumed to be circular.
-
(3)
The failure of the plastic zone of the rock mass meets the Mohr–Coulomb yield criterion.
-
(4)
The stress distribution of any roadway section is shown in Fig. 1. Figure 1 shows the plastic zone, viscoelastic zone and initial stress zone.
Calculation model of the plastic zone of a roadway
-
(1)
Stress distribution in the plastic zone
When the roadway is in the plane strain state, the static equilibrium equation of the plastic zone can be expressed as follows.
where \(\sigma_{r}^{p}\) is the radial stress of the plastic zone of the roadway; \(\sigma_{\theta }^{p}\) is the tangential stress of the plastic zone of the roadway; and \(r\) is the radius of the surrounding rock from the center of the roadway.
According to Mohr–Coulomb plastic conditions, Eq. (2) can be obtained.
where \(\varphi\) is the internal friction angle, and \(C\) is the cohesion.
By substituting Eq. (2) into Eq. (1), we obtain
By integrating Eq. (3), Eq. (4) is obtained.
On the wall of the roadway, \(r = r_{0}\) and \(\sigma_{r}^{p} = p_{i} (t)\). \(r_{0}\) is the radius of the roadway and \(p_{i} (t)\) is the support resistance changing with time. According to this boundary condition, we can obtain
Substituting Eq. (5) back into Eqs. (4), (6) can be obtained from Eq. (2).
-
(2)
Stress distribution in the viscoelastic zone
The constitutive equation of viscoelastic zone adopts the Poynting–Thomson standard linear solid model. The physical equation of the viscoelastic zone then can be defined by
where \(S_{ij}\) is the partial stress tensor; \(\mathop {S_{ij} }\limits^{ \bullet }\) is the time derivative of the stress partial tensor; \(\tau\) is the relaxation time; \(G_{\infty }\) and \(G_{0}\) are the long-term and instantaneous shear moduli respectively; \(e_{ij}\) is the partial strain tensor; and \(\mathop e\limits^{ \bullet }_{ij}\) is the time derivative of the strain partial tensor.
where \(\sigma_{ij}\) is the stress tensor, \(\sigma_{m}\) is the average stress, \(\sigma_{kk}\) is the body stress, \(\varepsilon_{ij}\) is the strain tensor, \(\varepsilon_{m}\) is the average strain, and \(\delta_{ij}\) is the Labradoroperator.
Regardless of the size of the plastic zone and whether or not there is support resistance, the stress at the interface of the viscoelastic zone and the plastic zone is constant. Let the calculated point \(r = aR_{0} (t)\), \(a\) be the scale factor, and \(R_{0} (t)\) the radius of the plastic zone changing with time. For the point changing with the radius of the plastic zone, the stress state will not change with time. Therefore, the term for stress rate in Eq. (7) is 0, and the average strain and average strain rate are considered to be 0. Equation 7 then can be reduced to
where \(\eta_{ret}\) is the delay time, \(\eta_{ret} = \tau G_{0} /G_{\infty }\), and \(P\) is the original rock stress.
In addition, in the case of axisymmetry, the volumetric deformation of the viscoelastic region is 0; according to the foregoing assumption, that is
Then, according to Eq. (10), we can obtain
where \(u\) is the displacement, and \(A(t)\) is a function of time.
Since the plastic zone changes with time, let \(r = R_{0} (t)\) in Eq. (11). We then obtain
By substituting the first and third formulas of Eq. (12) into the first formula of Eq. (9), and considering that \(\sigma_{r} = \sigma_{{R_{0} }} = P(1 - \sin \varphi ) - C\cos \varphi\), Eq. (13) can be obtained.
When t = 0, \(A(0) = \frac{{MR_{0}^{2} (0)}}{{4G_{0} }}\) and substituting this into Eq. (13), we can obtain the general solution of Eq. (13).
Substituting Eq. (14) into Eqs. (9), (11), and (12), respectively, the displacement and stress distribution of the viscoelastic region can be obtained.
where \(u_{r}^{c}\) is the displacement of the viscoelastic zone; \(\sigma_{r}^{c}\) and \(\sigma_{\theta }^{c}\) are the radial stress and tangential stress of the viscoelastic zone, respectively, and the letter \(c\) represents the abbreviation for creep.
Since Eq. (15) calculates the displacement and stress distribution of the viscoelastic region, r > \(R_{0} (t)\) in Eq. (15).
-
(3)
Plastic zone radius under creep conditions
According to the assumption of Eq. (10), the volume strain of the plastic zone is also 0, and the displacement of plastic zone can be calculated by reference to Eq. (11).
When \(r = R_{0} (t)\), we can obtain
Therefore,
Substituting Eq. (18) into Eq. (16), we can obtain
Equation 19 is similar to Eq. (15) but is different here in that r < \(R_{0} (t)\).
According to Eq. (19), the displacement of the roadway wall can be calculated. At this time, \(r = r_{0}\). Thus, we obtain
From Eq. (20), the range of plastic zones varying with time can be solved.
When calculating the radius of the plastic zone as a function of time according to Eq. (21), it is necessary to measure the displacement value of the wall of the roadway with time.
Equation (21) represents the case of no support. If the support resistance is considered, the support resistance is expressed by Eq. (22).
where \(p_{i} (t)\) is the support resistance.
The radius of the plastic zone as a function of time can be obtained from Eq. (22) and can be expressed by Eq. (23).
To calculate the radius of the plastic zone in Eq. (23), the relationship between support resistance and time needs to be measured. The support resistance increases with time and the increase in support resistance can prevent the deformation of the plastic zone. This paper considers the most unfavorable situation, that is, the existence of unsupported resistance. If the support is not considered, Eq. (23) can be replaced by the following equation when there is no measured displacement of the roadway wall.
where \(\delta\) is the long-term strength reduction coefficient of rock, \(K_{p} = 1.5P\) and \(k\) is the roadway relaxation time.
Equation (24) is a time-dependent model of the plastic zone radius under creep conditions.
Model application
Application background
The gas pressure of a coal seam is one of the important indexes for the prediction of the risk of coal and gas outbursts in coal mines. The bag-type grouting sealing method (Fig. 2) is a commonly used method to measure the gas pressure of coal seams in China. The test principle is as follows. When the capsular bags 1 and 2 are delivered to a predetermined position, the grouting pipe is used to grout the bag. After the capsular bags 1 and 2 expand, the bags and the borehole wall are in close contact. Then, with the grouting pressure exceeding the opening pressure of the grouting nozzle, the space between the two capsular bags is injected with grout to form a pressure chamber.
The key to accurately determining the gas pressure of the coal seam by the capsular bag grouting method is finding the range of the plastic zone of the roadway. If capsular bag 1 and capsular bag 2 are located in the plastic zone of roadway, the gas in the pressure-measuring chamber may leak out along the borehole wall and through cracks in the plastic zone during the pressure measurement, and the measured gas pressure will thus be smaller than the actual value or no indication on the pressure gauge. Conversely, when the positions of capsular bag 1 and capsular bag 2 are outside the plastic zone of roadway, the pressure-measuring chamber is ensured to be airtight and the correct gas pressure value will be obtained. Therefore, the accuracy of the calculated plastic zone of roadway can be approximately reflected by the measured value of the gas pressure.
Plastic zone range of a roadway at the experimental coal mine
The 11111 workface of the No. 13 Coal Mine of Pingdingshan Tianan Coal Mining Corporation Limited in the Henna province of China was selected as the research background. This specific coal mine is classified as a coal and gas outburst coal mine. The workface has an average thickness of 5.85 m and a coal seam inclination of 10°–19°. Considering the danger of coal and gas outburst, a floor rock roadway was used to construct the boreholes to drain the gas from the 11111 workface. The height and width of the roadway are 3 m and 4 m, respectively. The lithology of the floor rock roadway is marlstone and the roadway is 12–19 m away from the coal floor.
The buried depth of the floor rock roadway is approximately 600 m. The stress of the floor rock roadway is 15 MPa. The long-term strength reduction coefficient of marlstone is 0.3723. According to laboratory measurement, the relaxation time of marlstone is 1/24 d. By substituting the above parameters into Eq. (24), the relationship between the plastic zone radius of the floor rock roadway and time can be obtained, as shown in Fig. 3. Before the 100th day of the formation of the roadway, the plastic zone develops rapidly. During the 100th–150th day of roadway formation, the plastic zone development begins to slow down. After the 200th day of roadway formation, the plastic zone is basically stable and no longer expanding, and the radius of the plastic zone is approximately 7 m.
Therefore, during the sealing of the gas-pressure-measuring borehole for the 11111 workface, the sealing position of the starting capsular bag (capsular bag 1) should be at least 7 m away from the borehole orifice.
Gas pressure measuring at the experimental coal mine
The specific position of the gas pressure measuring borehole at the 11111 workface is shown in Fig. 4. There are 4 boreholes, and the boreholes are named borehole 1–borehole 4.
According to the calculation result of the plastic zone radius of the floor rock roadway, considering the certain margin coefficient, the starting capsular bag (capsular bag 1) is designed to be 8 m away from the borehole orifice when the pressure borehole is sealed. The gas pressure recovery curve measured by the pressure measurement process is shown in Figs. 5, 6, 7, 8. The final gas pressures of boreholes 1–4 are 2.71 MPa, 2.54 MPa, 2.30 MPa and 2.02 MPa, respectively.
In the process of coal seam gas pressure measurement, the direct gas content determination method was used to measure the gas content of the coal seam and the measurement results are shown in Table 1.
In addition, there is the following relationship between gas pressure and gas content24.
where \(W\) is the gas content, m3/t; \(a\) is the adsorption constant, 33.333 m3/t; \(b\) is also the adsorption constant, 0.767 MPa−1; \(P\) is the absolute gas pressure of the coal seam, MPa; \(Ad\) is the ash content of the coal, 9.00%; \(Mad\) is the moisture content of the coal, 1.56%; \(\pi\) is the porosity of the coal, 0.0637 m3/m3; and \(\gamma\) is the density of the coal, 1.42 t/m3.
Combined with the gas pressure, the gas content can be inversely calculated according to Eq. (25). The inverse calculation results are shown in Table 1.
It can be seen from Table 1 that the gas content calculated by the coal seam gas pressure is consistent with the measured gas content with a maximum phase difference of only 0.87 m3/t. This indicates that the measured value of the gas pressure is accurate and can correctly reflect the gas occurrence characteristics of the 11111 workface. Therefore, the validity of the plastic zone calculation model is also verified.
Conclusions
In this paper, a plastic zone model considering the creep effect is proposed. The model has good applicability and can be used to calculate the plastic zone under different in situ stress conditions and different roadway sizes.
According to the calculation model for the plastic zone, the radius of the plastic zone of a roadway in the No. 13 Coal Mine of Pingdingshan Tianan Coal Mining Corporation Limited in China is 7 m. Based on the radius of the plastic zone of the roadway, the starting capsular bag is designed to be 8 m away from the borehole orifice when the gas pressure measurement borehole is sealed in the roadway. The final measured gas pressure is 2.02–2.70 MPa. In the process of coal seam gas pressure measurement, the direct gas content determination method was used to measure the gas content of the coal seam. In addition, combined with gas pressure, the gas content can be inversely calculated. The gas content thus calculated from the coal seam gas pressure is consistent with the measured gas content, indicating that the gas pressure measurement value is accurate. The gas pressure measurement results verify the validity of the plastic zone calculation model.
Data availability
All data used to support the findings of this study are available from the corresponding author upon request.
References
Chen, H., Cheng, Y., Ren, T., Zhou, H. & Liu, Q. Permeability distribution characteristics of protected coal seams during unloading of the coal body. Int. J. Rock. Mech. Min. 71, 105–116 (2014).
Chen, H., Wang, Z., Chen, X., Chen, X. & Wang, L. Increasing permeability of coal seams using the phase energy of liquid carbon dioxide. J. CO2. Util. 19, 112–119 (2017).
Wang, L. et al. An improvement sealing method for gas pressure measurement by drilling a reducing-nipple borehole. Energy Sci. Eng. 6, 595–606 (2018).
Lu, S. et al. Numerical assessment of the energy instability of gas outburst of deformed and normal coal combinations during mining. Process. Saf. Environ. 132, 351–366 (2019).
Lu, S., Zhang, Y., Sa, Z. & Si, S. Evaluation of the effect of adsorbed gas and free gas on mechanical properties of coal. Environ. Earth Sci. 78, 218 (2019).
Lu, S. et al. Damage-induced permeability model of coal and its application to gas pre-drainage in combination of soft coal and hard coal. Energy Sci. Eng. 7, 1352–1367 (2019).
Galin, L. The plane elasto-plastic problem. J. Appl. Math. Mech. 10, 367–386 (1946).
Chen, L., Peng, J., Fan, W. & Sun, P. Analysis based on unified strength theory of non-uniform stress field circular roadway for plastic zone. J. China Coal. Soc. 32, 20–23 (2007).
Abdel-Meguid, M., Rowe, R. & Lo, K. Three-dimensional analysis of unlined tunnels in rock subjected to high horizontal stress. Can. Geotech. J. 32, 1208–1224 (2003).
Leitman, M. & Villaggio, P. Plastic zone around circular holes. J. Eng. Mech. 135, 1467–1471 (2009).
Guo, X., Zhao, Z., Gao, X., Wu, X. & Ma, N. Analytical solutions for characteristic radii of circular roadway surrounding rock plastic zone and their application. Int. J. Min. Sci. Technol. 29, 263–272 (2019).
Wang, P. et al. Evolution law of plastic zone of surrounding rock on deep soft rock roadway and its control. J. Hunan Univ. Sci. Technol. 34, 1–10 (2019).
Zhou, R. & Wu, D. Thickness analysis of soft coal rock loose circle based on FLAC3D numerical simulation. J. Qiqihar Univ. 34, 1–4 (2017).
Yang, S., Chen, M., Jing, H., Chen, K. & Meng, B. A case study on large deformation failure mechanism of deep soft rock roadway in Xin’An coal mine, China. Eng. Geol. 217, 89–101 (2017).
Zhang, H., Wang, Z., Zheng, Y., Duan, P. & Ding, S. Study on tri-axial creep experiment and constitutive relation of different rock salt. Saf. Sci. 20, 801–805 (2012).
Mishra, B. & Verma, P. Uniaxial and triaxial single and multistage creep tests on coal-measure shale rocks. Int. J. Coal. Geol. 137, 55–65 (2015).
Yang, S., Xu, P. & Ranjith, P. Damage model of coal under creep and triaxial compression. Int. J. Rock. Mech. Min. 80, 337–345 (2015).
Mansouri, H. & Ajalloeian, R. Mechanical behavior of salt rock under uniaxial compression and creep tests. Int. J. Rock. Mech. Min. 110, 19–27 (2018).
Xu, G., He, C., Yang, Q. & Wang, B. Progressive failure process of secondary lining of a tunnel under creep effect of surrounding rock. Tunn. Undergr. Space Technol. 90, 76–98 (2019).
Guo, Z. & Dong, F. Loose circle of surrounding rock and roadway support. J. Min. Saf. Eng. 3, 111–114 (1995).
Hu, S. & Liu, H. Leakage mechanism of coal seam gas drainage borehole and its application research progress. Saf. Coal Mines 47, 170–173 (2016).
Sun, Y., Lu, W., Yang, K., Zhou, C. & Huang, X. Study on disposal technology of air leakage around borehole by three-pouch closure device. J. Saf. Sci. Technol. 11, 67–72 (2015).
Construction Ministry of the People's Republic of China. Code for Design of Building Foundation. Beijing (2002).
National Coal Mine Safety Supervision Bureau. Code for Coal Mine Gas Drainage. Beijing (2006).
Acknowledgements
This work was supported by the National Natural Science Foundation of China (No. 51804102 and No. 51874122), the Henan Key Laboratory of Coal Green Conversion (CGCF201809, Henan Polytechnic University), the State Key Laboratory Cultivation Base for Gas Geology and Gas Control of Henan Province (WS2019B12, Henan Polytechnic University), the Program for Innovative Research Team of Henan Polytechnic University, and the Fundamental Research Funds for the Universities of Henan Province (NSFRF170901).
Author information
Authors and Affiliations
Contributions
H.C. and X.C. wrote the main manuscript; H.C. prepared all the figures; Z.W. carried out the data processing; Z.L. and F.A. carried out formula derivation. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chen, H., Chen, X., Wang, Z. et al. Plastic zone range of a roadway considering the creep effect. Sci Rep 10, 20341 (2020). https://doi.org/10.1038/s41598-020-77384-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-77384-5
This article is cited by
-
An innovative destressing technology and key parameters determination in both sides of a deep roadway
Scientific Reports (2023)
-
Numerical analysis and field tests of self-relief pressure zone ahead of mining workface: implications for methane extraction during underground mining
Environmental Earth Sciences (2023)
-
A calculation method of gas emission zone in a coal mine considering main controlling factors
Scientific Reports (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.