Research Article  Open Access
Xiaosheng Liu, Weijun Wang, Yizhong Luo, Jin Zhan, "Study on the Allowable Exposure Time of Temporary Goaf in the Open Stope and Backfill Mining Method with Large Structure Parameters", Advances in Civil Engineering, vol. 2020, Article ID 8853410, 9 pages, 2020. https://doi.org/10.1155/2020/8853410
Study on the Allowable Exposure Time of Temporary Goaf in the Open Stope and Backfill Mining Method with Large Structure Parameters
Abstract
Based on the uniaxial rheological experiment data of iron ore rock mass and filling body in largescale stope obtained by using the rock servo controlling creep equipment, the characteristics of the creep curve were analyzed. The creep models of iron ore and filling body which can show the attenuation creep of rock were constructed. And a nonlinear creep model was obtained. According to the rheological data of the iron ore and the filling body, the parameters of the new nonlinear creep model were identified to obtain the material parameters of the creep model. Then, the creep model parameters were fitted and reduced to calculate the reasonable exposure time of temporary goaf in the largescale stope. The results show that the reasonable exposure time of stope is 520 ～ 650 days, and the reasonable exposure time of filling body is 410 ～ 520 days. The model can well describe the initial attenuation creep stage and steady creep stage in the creep curve, which proves the correctness and rationality of the model. The study provides a reference for mining design and safety production of similar mines with largescale stope structure.
1. Introduction
In rock engineering, a large number of instability phenomena are related to the creep characteristics of rock. Rock creep is one of the important reasons for the deformation and instability of surrounding rock in rock engineering [1–3]. The rheological property of rock means that the rock presents deformation, fluidity, failure, and other properties with time under the conditions of external load and temperature, mainly presenting the properties of creep, stress relaxation, elastic aftereffect, strain rate effect, aging strength, rheological damage, and fracture [4–6]. In the past few decades, many scholars have carried out experiments on rock creep characteristics and obtained the conclusion that the rock creep continuously adjusts and recombines with time, resulting in rock stress and strain which also continuously increase and change with time [7–11].
In the process of ductile creep, the microcrack growth with time is also one of the main reasons for the failure of hard rock. In the process of brittle and ductile creep, rock samples show different timedependent crack growth patterns in the microcosmic mechanism, leading to significant differences in macroscopic deformation behavior [12, 13]. A series of uniaxial compressive strength tests and multistage creep tests are carried out on rock specimens under different immersion conditions. The results show that the peak strength and elasticity modulus of red sandstone not only decrease with water content but also decrease with immersion time, which can be better expressed by the negative exponential function. On the basis of studying the stressstrain curve of rock, many creep models of rock have been put forward by combining with various simple creep elements, such as Maxwell model, Bingham model, Burgers model, generalized Kelvin model, and KB model. But these models cannot describe the whole process of rock creep very well, especially for the accelerated creep stage. Among these models, the KB model and Burgers model are widely used and are improved to obtain good research results [13–18]. It should be pointed out that, because of the test instrument, the rock creep models under high temperature, high pressure, and anisotropy are seldom discussed, which is also the difficulty of rock creep research in the future.
In recent years, a large number of studies have been carried out on the cemented filling body strength, filling material proportion, backfill mechanical properties, and stability analysis; for example, the factors affecting the strength of tailings filling body are studied through different proportion tests [19], the damage model of tailings filling body is established through tests [20], and the mechanical properties in the roof of tailings cemented filling body are studied [21, 22]. However, there are few papers on the creep research for the cemented filling body in largescale stope mined by the backfill mining method. Therefore, the establishing creep model of filling body to determine the reasonable exposure time of goaf is not only of theoretical value but also of great necessity to calculate and analyze the stability of filling body in engineering practice. In the whole creep process of rocks and filling body, accelerated creep is relatively difficult to simulate. The current theoretical research focuses mainly on simulating accelerated creep characteristics. In this study, through a series of uniaxial compression creep tests, an ideal nonlinear material element which can simulate accelerated creep characteristics is constructed and then combined with basic elements to simulate the whole creep process. The uniaxial creep tests are carried out on the iron ore rock and filling body under the loading conditions to study the deformation characteristics of rock and filling body under different stress states. On the basis of summarizing the damage research results of different materials, statistical damage is introduced into the composite element model, and the corresponding creep model is established. And the corresponding creep parameters are obtained by model parameter identification. Then, the creep equation is used to predict the timeliness of largescale stope structure, and the reasonable stope exposure time is discussed from the safety.
2. Mining Conditions and Requirements of a Large Iron Mine
An iron mine in Anhui Province, China, is located about 18 km in the south of Huaihe River. The mining area is basically flat. The surface of the mining area is basically farmland and is part of the residential buildings. In order to safeguard the surface from collapse, strict protection is needed. The deposit is covered by Quaternary topsoil and strongly weathered rock stratum; its average thickness is 120 m, and the thickness of upper clay stratum is tens of meters. The integrity and stability of strongly weathered rock are poor. The occurrence of the ore body is steep, the local structural fracture zone is developed, and the rock and ore are stable. The buried depth of the main ore body is less than 150 ～ 900 m. The ore body trend is nearly northsouth and inclines to west, with a dip angle of 68° ～88°. The thickness of the ore body is 25.5 ～ 96.5 m, the average thickness is 49.65 m, the strike length is about 2000 m, and the average geological grade is 30.69%.
The required production scale of the mine is 7.5 million t/a, the mining capacity of the stope is 2000 t/d, the ore dilution rate is less than or equal to 10%, and the ore recovery rate is more than or equal to 90%. In order to protect the earth’s surface and the structures and to achieve the above technical and economic indicators, it is necessary to adopt the openstope and backfill mining method with large structure parameter. According to research, the method is divided into three parts: drilling in sublevel, ore drawing in level, and filling in level. The space is open during mining. After mining, the ground pressure is controlled and managed by the filling body and surrounding rock. The stope is vertical to the ore body trend, with a width of 40 m, a length of about 40–80 m, and a height of 100 m. The stope is divided into a room and pillar for mining. The widths of the room and pillar both are 20 m. There is no pillar between the stope and the stope. The sublevel height is 25 m, and the whole stope is divided into 4 sublevels to drill the upward fanshaped mediumdeep holes. The mining method is shown in Figure 1.
Because the room and pillar are open during mining, the key to successfully realize the mining method is to ensure that the stope does not collapse and largescale displacement before filling. The demonstration work is carried out through the creep principle and test.
3. Creep Experiment
There are many methods for rock creep experiment, and the most common ones are uniaxial compression experiment, triaxial compression experiment, and shear creep experiment. In order to determine the creep parameters of iron ore and filling body in the underground stope of an iron mine in a short period of time through the test, so as to study the time effect of rock stability, uniaxial compression creep experiment is chosen for the study.
3.1. Experiment Object and Uniaxial Compression Experiment
The iron ore belongs to quartz specularite. The filling body has three kinds of cementsand ratio of 1 : 4, 1 : 6, and 1 : 8. Due to the long time of creep experiment, this study chooses the quartz specularite and the filling body with cementsand ratio of 1 : 6 as the experiment objects.
Firstly, the conventional uniaxial compression experiment is carried out on the specimens, and the experiment results are shown in Table 1. Figure 2 shows the photos of creep experiment specimens.

(a)
(b)
3.2. Creep Experiment Instrument and Loading Mode
The strengths of the two kinds of specimens in this creep experiment are quite different. The uniaxial compressive strength of iron ore is relatively high, with an average value of 163.59 MPa, while that of the filling body is extremely low, and the average uniaxial compressive strength of the filling body with a cementsand ratio of 1 : 6 is 3.34 MPa in 28 d. So it is difficult to complete the test work on the same instrument for two kinds of specimens with great strength difference. Therefore, high requirements are put forward for the experiment instrument, and it is necessary to use different tonnage testing machine for testing.
For the iron ore with high uniaxial compressive strength, an Instron 1346 servo rigid material testing machine is used for testing. For the filling samples with very low uniaxial compressive strength, the Instron 1342 servo rigid material testing machine is used for testing, with the maximum axial pressure of 250 kN and the sensor of 10 kN, and its accuracy can meet the observation value of 0.001 mm. The load control mode is adopted in the test, the test data are recorded automatically, and the test curve is drawn synchronously.
The loading mode is the constant loading mode. A group of specimens with the same rock are used to bear different constant loads. Because of the different properties of the rock specimens, the test results are discrete. However, the test results are reliable and not affected by the loading state. Therefore, the uniaxial creep experiment adopts constant load mode.
4. Creep Experiment and the Model of Ore and Rock
The uniaxial creep experiment of iron ore is carried out by using the Instron 1346 servo rigid material testing machine. Before loading, the specimens are placed indoor for 20 ～ 30 days in the natural environment. The temperature of the experiment environment shall be controlled at 20°C ± 2°C. The constant load mode is adopted. Due to the high strength of iron ore, the loading rate is 90 kN/min, which is about 0.75 MPa/s. The axial displacement is recorded by using a displacement extensometer with a gauge distance of 5 mm. The control mode is load control. The creep experiment of iron ore is shown in Figure 3.
After loading the predetermined constant load, the strain is read immediately as the initial strain value of the load. The strain values are recorded at the intervals of 5, 10, 15, and 30 min in the first 1 h and then are recorded at intervals of half an hour.
Under the action of constant load (81.80 MPa) with a stress level of 50% R, specimen Y26 does not damage in 192 hours. Specimen Y27 also is not damaged in 192 hours under the action of constant load (98.15 MPa) with the stress level of 50% R. Under the action of constant load (114.51 MPa) of a stress level of 70% R, specimen Y28 is not damaged in 192 hours but presents a sign of possible damage. Under the action of constant load (130.87 MPa) with a stress level of 80% R, specimen Y29 gets damaged in 175 hours. The creep of each test specimen is shown in Table 2. The test curve of the relationship between axial strain and time of each test specimen is shown in Figure 4.

The iron ore belongs to hard rock, its creep is relatively small at low stress level, and its deformation tends to be stable after a short period of initial creep. But at high stress level, the creep is relatively obvious, which is also proved by this experiment. In the experiment, due to a good control on the load, the accelerated creep process is observed in two specimens besides the attenuation creep process and steady creep process.
The five element model used in rock mechanics is called the BK model, also known as the Bingham–Vogt model, which is composed of a Bingham body and a Kelvin (Vogt) body in series.
When σ < σ_{s}, the friction plate is a rigid body. This model has the same characteristics as the generalized Kelvin model. The constitutive equation and creep equation of the model are as follows:
When σ < σ_{s}, the creep of the BK model is in a stable creep. According to formula (2), when t = 0, only the Bingham body works, and when t ⟶ ∞, the damper cylinder does not work, which is equivalent to connecting the spring of the Bingham body and Kelvin (Vogt) body in series.
When σ ≥ σ_{s}, the performance of the BK model is similar to that of the Burgers model, but the difference is that the stress overcoming the friction resistance σ_{s} in the model should be removed in the BK model. Therefore, by directly substituting (σ−σ_{s}) for σ in the constitutive equation of the Burgers model, the constitutive equation of the BK model is obtained as follows:
The creep and unloading equations of the BK model can be obtained from the creep and unloading equations of the Kelvin body and Bingham body:
The relaxation equation of the BK model is too complex to be listed here. From equation (4), when t = 0, only the spring works, and . When t ⟶ ∞, ⟶ ∞, but the change rate tends to be constant, and the creep belongs to metastable creep. Similarly, the unloading equation of the BK model can be obtained by the unloading equation superposition of the Maxwell body and Kelvin body. When t ⟶ ∞, , and some permanent deformation remains after unloading.
5. Creep Experiment and the Model of Filling Body
The uniaxial creep experiment of filling body is carried out by using the Instron 1342 servo rigid material testing machine. Before loading, the specimens are placed in the indoor for 20 ～ 30 days in the natural environment. The temperature of the experiment environment shall be controlled at 20°C ± 2°C. The constant load mode is adopted, and the loading rate is 2 kN/min, which is about 0.015 MPa/s. The axial displacement is recorded by using an LVDT displacement sensor. The control mode is load control. The creep experiment of filling body is shown in Figure 5.
After loading the predetermined constant load, the axial displacement is read immediately as the initial displacement value of the load. In the first 1 h, the axial displacement values are recorded at the intervals of 5, 10, 15, and 30 min and then are recorded every half an hour.
Under the action of constant load (1.00 MPa) with a stress level of 30% R, specimen Yc67 is not damaged in 200 hours. Specimen Yc68 also is not damaged in 200 hours under the action of constant load (1.34 MPa) with a stress level of 40% R. Under the action of constant load (1.67 MPa) with the stress level of 50% R, specimen Yc69 is not damaged in 200 hours but presents a sign of possible damage. Under the action of constant load (2.00 MPa) with a stress level of 60% R, specimen Yc610 gets damaged in 196 hours. The creep of each test specimen is shown in Table 3. The test curve of the relationship between axial strain and time of each specimen is shown in Figure 6.

The strength of filling body is low, and its creep deformation is large. When the stress level is more than 50% R, the creep is obvious, which experiences three stages: attenuation creep, steady creep, and accelerated creep. There are obvious signs before the material fracture. In the creep curve, there is an upwarping section in rapid increase of change rate.
The Burgers model can well describe the creep curve with initial creep and steady creep, and the model is simple and practical, which has been widely used. Therefore, the Burgers model is selected to describe the compression creep characteristics of the filling body.
The Burgers model is also called the MK model, which is composed of a Maxwell body (M body) and a Kelvin body (K body) in series.
The constitutive equation and creep equation of the Burgers model are as follows:
According to equation (7), the creep equation of the Burgers model can be obtained by the creep equation superposition of the Maxwell body and Kelvin body. When t = 0, only the spring works, and . When t ⟶ ∞, approaches infinity, but its change rate tends to be constant and the creep belongs to metastable creep.
6. Creep Model Parameters and Exposure Time Calculation
The creep constitutive equations of iron ore and backfill in the following sections under uniaxial compression belong to a nonlinear function. In order to determine the parameters of the constitutive equations, the Marquardt method in the least squares method is used to fit the parameters of the creep model.
6.1. Fitting Value of Creep Parameters of Iron Ore
The creep constitutive model of iron ore uses the KB model for fitting. In the KB model, there are 5 rheological parameters need to be calculated, in which EB is the elastic modulus. When σ acts on the specimen, the elastic deformation caused by EB is completed in an instant. Therefore, when t = 0, the corresponding ratio of σ to ε is equal to the EB value. The initial flow stress can be read directly from the creep curve. The parameters to be determined by the least squares method are EK, ηB, and ηK. The rheological parameters are related to the stress level. It is necessary to fit the creep parameters under different stress levels. During fitting the creep curves of specimens Y26 ∼ Y210, the correlation coefficients are 0.92, 0.90, 0.93, 0.91, and 0.92 respectively. Table 4 shows the calculation results of creep parameters of iron ore.

6.2. Fitting Value of Creep Parameters of the Filling Body
The creep constitutive model of filling body adopts the Burgers model for fitting. In the Burgers model, there are four rheological parameters need to be calculated, in which EB is the elastic modulus. When σ acts on the specimen, the elastic deformation caused by EB is completed in an instant. Therefore, when t = 0, the corresponding ratio of σ to ε is equal to the EB value. Comparing the Burgers model with the KB model, it is found that the elastic modulus EB of the two models is exactly the same. The other three rheological parameters (EK, ηB, and ηK) are needed to be determined by the least squares method. The rheological parameters are related to the stress level. Similarly, it is necessary to fit the creep parameters of filling body under different stress levels. The correlation coefficients of the creep curves of specimens Yc67 ∼ Yc611 are 0.95, 0.93, 0.93, 0.92, and 0.94, respectively. Table 5 shows the calculation results of creep parameters of the filling body.

6.3. Recommended Creep Parameters
Because the time of indoor creep experiment is limited, the creep parameters of rock and filling body obtained by tests must have differences with the real creep parameters of rock and backfill onsite. Therefore, the fitting value of creep parameters of iron ore and filling body should be modified. The general method is to make corresponding reduction through reference to similar projects and combining with the engineering geological conditions or viscoelastic displacement inversion results.
Considering the difference between the onsite rock mass and the indoor rock specimens, referring to the research results of similar projects at home and abroad, and combined with the engineering geological conditions of the iron ore and the filling body in the mine, the fitting values of creep parameters in Tables 4 and 5 are reduced. As the joint fissures in iron ore are relatively developed, the reduction coefficient for the iron ore is selected as 0.7. For the filling body, the development degree of the joint fissures is relatively low, and the reduction coefficient is selected as 0.8. The creep parameters of onsite iron ore body and filling body after modification are shown in Tables 6 and 7, respectively.


6.4. Reasonable Exposure Time of LargeScale Stope Structure
On the basis of the indoor uniaxial compression creep experiment, the creep models of iron ore (quartz specularite) and filling body with a cementsand ratio of 1 : 6 are established by analyzing the characteristics of the creep curve. Then, the model parameters are fitted by the least squares method, and the corresponding creep parameters are obtained. Then, the fitting value of creep parameters is modified for the engineering application. So the creep equation can be used to calculate the stability time of iron ore and filling body. This stability time is considered as reasonable exposure time.
The reasonable exposure time for the iron ore and the filling body with the cementsand ratio of 1 : 6 is from the start of stope blasting to the end of roofcontacted filling. The initial setting time of the filling slurry is about 6–10 hours, and the final setting time is about 2 days. During this time, the filling body gradually develops strength and supports the roof of the stope. Because the solidification time of filling slurry is relatively short, it can be neglected in this study.
The results of the uniaxial compression creep tests are obtained as follows: the ratio of axial ultimate strain to conventional ultimate strain is 93% ∼ 95% in the case of creep failure of iron ore and is 94% ∼ 96% in the case of creep failure of the filling body.
The creep deformation of iron ore and filling body is related to the stress level. The finite element numerical simulation shows that the stress level of iron ore is generally lower than 100.0 MPa, namely, lower than 60% of the uniaxial compressive strength of iron ore. And the stress level of filling body is generally lower than 2.0 MPa, namely, lower than 60% of the uniaxial compressive strength of filling body.
The axial ultimate strain is 0.49 in the creep failure of iron ore, and the ratio of the axial ultimate strain to the conventional ultimate strain is about 94%. The axial ultimate strain is 0.55 in the creep failure of filling body, and the ratio of the axial ultimate strain to the conventional ultimate strain is 95%. The stress level is 50% R ∼ 60% R. According to formulas (2) and (7), the stability time of iron ore and filling body is calculated, respectively.
6.4.1. Stability Time of Iron Ore
The equation is established according to the parameters in Table 6 and the relevant initial values. When the stress level is 50% R, the equation is as follows:
When the stress level is 60% R, the equation is as follows:
The stability time of iron ore is calculated as 658.6 days through solving equation (8) and as 527.8 days through solving equation (9). So it is suggested that the reasonable exposure time of iron ore in the stope is 520 ～ 650 days.
The strength of quartz specularite in this mine is lower than that of magnetite in other mines. If the exposure time is determined based on the quartz specularite in the mining area, the magnetite in other mines will be more stable in the same exposure time.
6.4.2. Stability Time of the Filling Body
The equation is established according to the parameters in Table 7 and the relevant initial values. When the stress level is 50% R, the equation is as follows:
When the stress level is 60% R, the equation is as follows:
The stability time of filling body with the cementsand ratio of 1 : 6 is calculated as 523.2 days through solving equation (10) and as 417.5 days through solving equation (11). So it is suggested that the reasonable exposure time of filling body in the stope is 410 ～ 520 days.
There are three different cementsand ratios for the filling slurry, including 1 : 4, 1 : 6, and 1 : 8. Among them, the filling body with the cementsand ratio of 1 : 6 is widely used, with a reasonable strength, and its research results provide a good reference for the filling body with other two kinds of cementsand ratio.
7. Conclusion
(1)Uniaxial compression creep tests are carried out on iron ore and filling body, respectively. The creep time of iron ore is 795 hours and that of filling body is 889 hours. The total time of two kinds of samples is 1684 hours.(2)Through the analysis on the characteristics of creep curves, the creep models of rocks (quartz specularite) and filling body with a cementsand ratio of 1 : 6 in the largescale stope are established, respectively. The parameters of the models are fitted by the least squares method to obtain the corresponding creep parameters. Then, the fitting value of creep parameters is modified for the engineering application. According to the stress environment of iron ore and filling body, the creep equations are used to calculate the stability time of iron ore and filling body. The reasonable exposure time of largescale stope is determined from the aspect of safety.(3)The stability time of quartz specularite at the stress level of 50% R ～ 60% R (R is the uniaxial compressive strength) is 527.8 ～ 658.6 days, so it is suggested that the reasonable exposure time of the largescale stope is 520 ～ 650 days. The stability time of the filling body with the cementsand ratio of 1 : 6 at the stress level of 50% R ～ 60% R is 417.5 ～ 523.2 days, so it is suggested that the reasonable exposure time of the filling body is 410 ～ 520 days.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
References
 N. Brantut, M. J. Heap, P. G. Meredith, and P. Baud, “Timedependent cracking and brittle creep in crustal rocks: a review,” Journal of Structural Geology, vol. 52, pp. 17–43, 2013. View at: Publisher Site  Google Scholar
 X. Feng, B. Chen, S. Li et al., “Studies on the evolution process of rockbursts in deep tunnels,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 4, no. 4, pp. 289–295, 2012. View at: Publisher Site  Google Scholar
 X.T. Feng, S.F. Pei, Q. Jiang, Y.Y. Zhou, S.J. Li, and Z.B. Yao, “Deep fracturing of the hard rock surrounding a large underground cavern subjected to high geostress: in situ observation and mechanism analysis,” Rock Mechanics and Rock Engineering, vol. 50, no. 8, pp. 2155–2175, 2017. View at: Publisher Site  Google Scholar
 A. Nicolas, J. Fortin, J. B. Regnet et al., “Brittle and semibrittle creep of tavel limestone deformed at room temperature,” Journal of Geophysical Research: Solid Earth, vol. 122, no. 6, pp. 4436–4459, 2017. View at: Publisher Site  Google Scholar
 J. Renner, B. Evans, and G. Siddiqi, “Dislocation creep of calcite,” Journal of Geophysical Research: Solid Earth, vol. 107, no. B12, pp. 6–16, 2002. View at: Publisher Site  Google Scholar
 N. Brantut, P. Baud, M. J. Heap, and P. G. Meredith, “Micromechanics of brittle creep in rocks,” Journal of Geophysical Research: Solid Earth, vol. 117, no. B8, 2012. View at: Publisher Site  Google Scholar
 C. H. Scholz, “Mechanism of creep in brittle rock,” Journal of Geophysical Research, vol. 73, no. 10, pp. 3295–3302, 1968. View at: Publisher Site  Google Scholar
 D. M. Cruden, “A theory of brittle creep in rock under uniaxial compression,” Journal of Geophysical Research, vol. 75, no. 17, pp. 3431–3442, 1970. View at: Publisher Site  Google Scholar
 R. J. Martin, “Timedependent crack growth in quartz and its application to the creep of rocks,” Journal of Geophysical Research, vol. 77, no. 8, pp. 1406–1419, 1972. View at: Publisher Site  Google Scholar
 B. K. Atkinson, “Subcritical crack propagation in rocks: theory, experimental results and applications,” Journal of Structural Geology, vol. 4, no. 1, pp. 41–56, 1982. View at: Publisher Site  Google Scholar
 N. L. Carter and F. D. Hansen, “Creep of rocksalt,” Tectonophysics, vol. 92, no. 4, pp. 275–333, 1983. View at: Publisher Site  Google Scholar
 J. Zhao, X. T. Feng, X. Zhang, and C. Yang, “Brittle and ductile creep behavior of Jinping marble under true triaxial stress,” Engineering Geology, vol. 258, 2019. View at: Publisher Site  Google Scholar
 C. Yu, S. Tang, C. A. Tang et al., “The effect of water on the creep behavior of red sandstone,” Engineering Geology, vol. 253, pp. 64–74, 2019. View at: Publisher Site  Google Scholar
 N. Cristescu and U. Hunsche, Time effects in rock mechanics, vol. 350, Wiley, New York, NY, USA, 1997.
 Ö. Aydan, T. Ito, U. Özbay et al., “ISRM suggested methods for determining the creep characteristics of rock,” Rock Mechanics and Rock Engineering, vol. 47, no. 1, pp. 275–290, 2014. View at: Publisher Site  Google Scholar
 C. F. Yang and Y. J. Zeng, “Investigation on creep damage constitutive theory of salt rock,” Chinese Journal of Rock Mechanics and Engineering, vol. 21, no. 11, pp. 1602–1604, 2002. View at: Google Scholar
 L. S. Tsai, Y. M. Hsieh, M. C. Weng, T. H. Huang, and F. S. Jeng, “Timedependent deformation behaviors of weak sandstones,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 2, pp. 144–154, 2008. View at: Publisher Site  Google Scholar
 Z. Tomanovic, “Rheological model of soft rock creep based on the tests on marl,” Mechanics of Timedependent Materials, vol. 10, no. 2, pp. 135–154, 2006. View at: Publisher Site  Google Scholar
 K. Zhao, X. Yu, and S. T. Zhu, “Acoustic emission investigation of cemented paste backfill prepared with tantalum–niobium tailings,” Construction and Building Materials, vol. 237, Article ID 117523, 2020. View at: Publisher Site  Google Scholar
 K. Zhao, Q. Wang, S. Gu et al., “Mining scheme optimization and stope structural mechanic characteristics for a deep and large ore body,” Journal of Minerals, vol. 71, no. 11, pp. 4180–4190, 2019. View at: Publisher Site  Google Scholar
 K. Zhao, H. B. Yan, and X. Feng, “Stability analysis of pillar based on energy law,” Chinese Journal of Theoretical and Applied Mechanics, vol. 48, no. 4, pp. 976–983, 2016. View at: Google Scholar
 K. Zhao, Q. Wang, Y. Yan et al., “Geological structural surface evaluation model based on unascertained measure,” Geofluids, vol. 2019, pp. 1–7, 2019. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Xiaosheng Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.