A Predictive Model of Mining Collapse Extent and Its Application

Nan, Jia; Liu, Cheng; Liu, Yi

doi:https://doi.org/10.1155/2019/5184287

Advances in Civil Engineering

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 5184287 | https://doi.org/10.1155/2019/5184287

A Predictive Model of Mining Collapse Extent and Its Application

Jia Nan,¹Cheng Liu,¹and Yi Liu¹

Academic Editor: Mario D’Aniello

Received13 Aug 2018

Revised16 Oct 2018

Accepted25 Oct 2018

Published04 Mar 2019

Abstract

To reveal the mechanical behavior mechanism of collapse and to control risks effectively, the instability extent of the collapse area was established through theoretical mechanics and numerical methods, taking one metal mine as a case study; on this basis, a routine reinforcement program was determined, and the effect of the program was evaluated. The results show the following. (1) Analytical formulas of the critical slip angle and the collapse height of the ore body were derived by the mechanics method, and the rock mechanics parameters were obtained by field coring and physical and mechanical experiments. The slipping line angle increases along with uniform force Q and is inversely proportional to the bending stiffness. Meanwhile, the calculation formula for the maximum subsidence of ore body was deduced. (2) Numerical results can be used to determine the basic form of the collapse area, and a “U-shaped” collapse area formed when a plastic area passed completely through, resulting in the overall destruction. (3) The grouting reinforcement program includes “determining the instability region ⟶ roadway temporary support ⟶ improve the water environment and surrounding rock bearing capacity ⟶ mining planning” which were determined on the basis of prediction. (4) The hierarchical structure of the rock body and filling were improved combined with the Delphi method, and the grouting effect evaluation model was constructed and verified using the improved FD-AHP method; the evaluation value indicating that the grouting reinforcement improved the bearing capacity of ore body and filling body in collapse area. The research results provide systematic reference and technical support for the analysis of stope collapse mechanism, prediction of hidden trouble, and the subsequent mining.

1. Introduction

Mining safety is an important part of the public safety, and reasonable and effective use of mineral resources is the responsibility of the whole society [1]. With the rapid development of China’s social economy, the consumption of mineral resources is also increasing, and high-grade easy-mining mineral resources is increasingly depleting [2]. To meet the demand of economic construction for sustained growth of mineral resources and maintain the steady growth of economic development, under the current technical conditions, for the exploitation of mineral, in addition to continuing to expand to deeper regions, it is necessary to fully recycle the hidden resources caused by the difficult mining or abandoned mining in the early stage of mining [3]. The mining collapse problems bring hidden hazards for the subsequent mining [4]. Prediction of the slip boundary and its influence [5] along with the deformation and subsidence of ore body (filling body) [6] will affect the further exploitation of hidden resources directly.

Facing the extremely complicated mining environment, recovering the hidden resources through determining mining collapse extent, and reinforcement are crucial to mining safety. At present, the main methods to study the instability and collapse of stopes include numerical simulation analysis, catastrophe theory, similar material simulation, and uncertainty method [7–9]. Yao [10], followed by Wood [11], employed the finite element method and the boundary element method to study the conditions and heights of overlying rock collapse caused by stoping. Freidin et al. conducted stress distribution and stability analysis studies on goaf roofs with finite element software and obtained the safe range of underground mining span [12]. Peralta adopted the perturbation method to study the stability of fluid-filled roadways under the influence of mining and analyzed the change pattern of impacts exerted by different dynamic waves on existing structures. Grig studied the mechanical behaviors of iron ore and marl based on rock mechanics tests with FLAC 2D software [13]. Malan proposed an elastic-viscoplastic method to study the time dependence of creep deformation of the ore rock fault zone during surface excavation works [14]. With respect to stoping of resources with hidden dangers and the countermeasures to hidden dangers, developed countries studied the exploitation of resources containing hidden dangers in underground mines much earlier than China. With a completely different concept being adopted, residual ore was recovered fully and effectively from San Manuel Copper Mine and White Pine Copper Mine in the United States by using the dissolving leaching mining method, which was safe and had low mining cost [15]. Tan et al. studied the stress distribution pattern of semi-U-shaped valley terrains [16].

Before mining of the hidden resources, the most critical step is to determine the influence range of mining stability, especially the location and angle of the slipping line; subsequently, the hidden area can be managed [16]. Meanwhile, there are no general approaches or methods for references [17, 18]. Therefore, to improve mining safety and reveal the mechanical behavior mechanism of collapse and control risks, combined with a metal ore mine as a case study, in view of the complex mining conditions of the collapsed area, the theoretical mechanics and numerical analysis model are built to calculate the influence extent of collapse area through field coring and laboratory physical and mechanical experiments. Based on the rock pressure and rock mechanics theories to predict the sinkage of collapsible ore body, the numerical calculation model was established to verify the shape and the extent of collapsed area. Next, the grouting reinforcement program can be determined, and the effect of the program can be evaluated by the comprehensive evolution FD-AHP method.

Engineering background of the metal mine is as follows. The range of vertical collapse is from elevation −298.3 m to −232.0 m, north-south coordinates from 290.5 to 420.0 (Figure 1), with the width of each stope of approximately 6.5∼8 m. The mining has been finished and filled (the remaining stakes are interpillar stope), and the second half of the 0# stope is approximately 22.5 m high. The collapse of the goaf caused the filling of two sides into the goaf and the upper ore body subsidence.

2. Theoretical Prediction Model of the Mining Collapse Extent

2.1. Mechanical Calculation Model of the Slipping Line Angle

Firstly, the theoretical mechanics model with the principle of deletion and simplification is established within the collapse range, where the unit length is taken along the circumference and the radial direction is taken to the boundary of the collapse zone, cutting out a long beam along the axial direction. The line displacement and angular displacement on the side connecting with the mine-out area before collapse are constrained; the constraint can be simplified into a fixed end. The line displacement along the axial displacement on the other side is constrained, and the angle is 0; this constraint can be simplified as a siding bracket. Suppose that, during the collapse, excavation unloading is evenly distributed. In the simplified stability analysis of beam-column structure, the axial stress and radial stress are the dominant factors that affect the collapse. The tangential stress is the secondary factor, and the tangential stress will affect the geometry of the beam-column structure and physical structure. The mechanical model of the collapse zone can be abstracted as shown in Figure 2.

Then, the rock mass structure within a certain range of the collapse zone can be simplified as a one-time statically indeterminate structure, and the unknown parameters corresponding to the excess unknown force can be solved using the force method, where F indicates axial pressure of the surrounding rock (or filler) in the collapse zone, and Q indicates radial load. As shown in Figure 3, it is assumed that the constraint force on F side is an extra unknown constraint force. By the principle of superposition, we can see that, for the statically indeterminate structure, the mechanical state is the integrated stacking result in the radial uniform force Q and the axial uniform force F (in Figure 4, this mechanical state can be converted into determinate structure in the radial uniform force Q and axial force F).

(a)

(b)

The bending moment M₀ of any section J can be obtained by the cross section method, as shown in Figure 5 (X-Z coordinate system):

The mechanical model of the statically indeterminate problem can be established from the approximate differential equation and the boundary condition of the deflection curve:

The general solution of the differential equation iswhere H is the maximum deflection. Thus, k, , C, and G are obtained by the following boundary conditions:

When tends to be 0, , the equations of the deflection curve and the corners are as follows:where

Similarly, under the effect of concentration force, M₀ = −M_e0 − F(H − Z), with the abovementioned curve differential equation and boundary condition, the simplified mechanical model of collapse area are as follows:

However, its general solution can be changed intowhereN and M are integral constants, and H is the maximum deflection; N and M can be obtained from the boundary conditions:where

For , the solution to the problem becomes

For the abovementioned statically indeterminate problem, according to the coordination equation, when and , can be obtained by the following equation:

According to the superposition principle, the deflection curve and the rotation equation are as follows:

The slip angle θ₁ is as follows:where the Q value is 10.6 MPa, l represents 40 m, E value is 2.7 MPa, I value is 1000, and x value is 8. Substituting the data into the above equations (14) and (15), the slip angle of the collapse zone can be estimated to be 42°.

To obtain the objective calculation parameters of rock mass, the core is obtained from field drilling, and the mechanical parameters of rock are obtained through mechanical experiments and then are transformed into rock physical and mechanical parameters (Figure 6).

(a)

(b)

It can be concluded that the slipping line angle increases with the increase of uniform force Q and is inversely proportional to the bending stiffness; the slip angle increases with the influence range of the collapse.

2.2. Prediction of Subsidence Displacement of Ore Body

Both sides of mined-out area become unstable due to excavation, and it will enlarge the span of the mined-out area and make the ore body fall. The collapse displacement calculation schematic of mined-out area is shown in Figure 7. The lower half of the shaded area is the falling ore or backfill, the upper half of shaded area is the ore body above the mining area, and the middle part is the hypothetical sinking area. Because the 0# stope in Figure 7 is asymmetrical on both sides, assume the model is tilted. The red square in Figure 7 corresponds to the purple area in Figure 1.

According to the Mohr–Coulomb strength theory [19], A₁ and A₂ are the increased span caused by the instability of the mine-out area, and they can be expressed by the following equation [19]:where D_k is the thickness of ore body and φ is the internal friction angle of the ore body.

The volume V_a formed after the left body instability in the collapse area is

The volume V_b formed after the right body instability in the collapse area is

The volume V_c formed after the collapse of the ore body and filling body in a certain range on the top iswhere k_d is the expansion coefficient of the ore body and A is the width of the original mined-out area.

Because the fallen ore body and filling body piled into the mined-out area, the pile height H can be obtained by the following three equations:where is the oblique width of the top surface of the collapse accumulation body in the mined-out area; H₀ is the height of the collapse accumulation body along the direction of the ore body in the mine-out area; H_z is the total height of the accumulation body in the mine-out area; λ is Poisson’s ratio; and β is the collapse slip angle.

Equation (24) can be obtained by simultaneous formulas (20)–(23):

Assuming that, the collapsed ore bodies fill the entire mined-out area and make the area to be completely blocked, the necessary condition is as follow:

The maximum amount of subsidence of the ore body in the collapse area can be calculated from the above equation:

Therefore, the maximum subsidence of the ore body is related to ore body thickness, internal friction angle, density and dilatancy coefficient, width of the original goaf, ore body inclination angle, and the geometric shape of ore body. The obtained calculation formula can be used as a reference for collapse analysis of similar engineering stopes.

2.3. Numerical Verification Analysis on Collapse Area

The three-dimensional numerical model is established according to the current production map, which is 200 m long, 100 m wide (ore body direction), and 100 m high. In the model, the backfilling waste rock can be simulated using the double yield constitutive model, and the rest are simulated using the Mohr–Coulomb model [20]. The interface between the filling and the backfill is connected by the interface command and uses the complete mesh (187500 units and 199576 nodes were formed in total). Two of the horizontal stresses can be expressed in the form of the pressure coefficient (λ = ν/1 − ν). Meanwhile, rock mass mechanical parameters were obtained by experiment (Table 1).

Based on the results of the inversion calculation, the average elastic modulus is approximately 3.75 GPa, with Poisson’s ratio of 0.25, the uniaxial compressive strength of 24 MPa, the proportional limit of 34.6 MPa for the surrounding rock, ore body and filling body with depth of approximately 300 m, the cohesion of 3.6 MPa, the friction angle 36°, the axial horizontal stress of 6.96 MPa, and the vertical stress of 11.2 MPa.

The excavation of the underground space may lead to the destruction of the ore body, with the characteristics of the size of the ore body plastic zone (to avoid size effect, the model is 3 times the size of the model in Figure 1). Figure 8 shows the distribution of the plastic zone formed after the excavation of the 0# stope. The figure reveals that the plastic zone has a wide range of distribution and a plastic zone with the existence of shear and tensile failure. The red outline is the slip line. The plastic zone in the mined-out area roof and both sides expand to the adjacent stope and result in perforation, leading to large-scale instability of the mined-out area. Therefore, it can be known that rib spalling and collapse exist in both sides of the mined-out area because of transverse tension, resulting in the sheeting bodies and collapse along with caving of the roof of various degrees, eventually leading to the instability destruction of the wide area in the mined-out area.

Figure 9 shows the maximum shear strain increment cloud in the state of collapse, revealing that a “U-shaped” collapse area was formed, and the area with the maximum shear strain having the most serious damage.

3. Grouting Reinforcement Measures

In the subsequent stoping works, guaranteed high-quality filling was the precondition for improvements to the lithology of surrounding rocks in the collapse area, whereas the enhanced bearing capacity of the backfill was a key factor for inhibiting the further sinking and destruction of the backfill. Based on the results of numerical analysis, grouting was adopted to improve the strength of the backfill in the U-shaped collapse area effectively, with the overall stability being enhanced. To ensure successful improvements, it was critical to create a better hydrological environment of surrounding rocks in goafs and enhance the self-bearing capacity of these rocks; safe and efficient stoping could be guaranteed by taking timely and effective measures to improve the overall bearing capacity of the backfill. Therefore, the basic principles for controlling surrounding rocks in the collapse area were proposed as follows: (1) pretreatment should be conducted on the water environment in the collapse area; (2) advance prereinforcement should be implemented on the surrounding rocks in goafs; (3) advance support and temporary support should be provided on the ore removal roadway at the bottom; and (4) improvements should be made to the overall stability of the backfill. The overall plan is as shown in Figure 10.

The plastic zone of the surrounding rock has a large range, and its main failure mode was shear failure or shear and tensile comprehensive damage. The bearing capacity of the surrounding rocks could be increased through timely effective grouting reinforcement. Grouting can enhance strength of surrounding rocks and stability of surrounding rocks in future mining practice. Grouting reinforcement was performed on the largest shearing strain cut through belt within the slip line of the unstable collapse area (including crushing zone and plastic zone). Every time after the excavation works were completed, the surrounding rocks in goafs should be filled with high-strength backfill to ensure the development of sufficient bearing capacity, thus, restraining the further sinking and deformation of the surrounding stopes. With respect to the specific grouting situation, as shown in Figure 11, grouting holes were arranged in those highly safe drilling chambers outside the collapse area, which were used to detect if the adjacent ore bodies had sunk as well as to block and reinforce the voids and cracks of backfill below the ore bodies, ensuring that the backfill was stable without further collapse in the follow-up stoping works.

4. Inspection Grouting Reinforcement Effects

To judge whether the stope in the collapse area could be safely mined or not, detection of the grouting effects was conducted. In this paper, the grade of ore body and backfill were classified as a multiattribution decision-making problem. Using the combination of the RMR method, BQ method [20], and Delphi group decision-making theory, assessment factors influencing ore bodies and backfill were selected according to the practical situations of stopes, making these factors to be better in line with the uncertainties in the real situation.

Comprehensive analysis of the RMR system and the BQ system revealed that the RMR system did not consider the influences of mechanical parameters, such as fracture and elasticity modulus. However, in the backfill, because of the changes in internal structure of ore body during waste filling and collapse, fractures were formed, greatly changing the mechanical properties of ore body. According to results of a field test, a water pressure test could detect the presence of voids and fractures in the ore body and the backfill. Therefore, the mechanical indices of the BQ system and the results of the water pressure test were included in the classification system as evaluation indices. Meanwhile, no joint spacing existed in the backfill, and it was difficult to determine the joint spacing because of the engineering environmental conditions. Given this situation, in this paper, regarding the selection of the hierarchical structure to be built, the joint spacing in the RMR system was ignored, with the state of ground stress being considered.

In summary, the following index layers were used in the ore rock and backfill classification method: ① overall rock mass strength; ② rock quality index (RQD); ③ mechanical indexes in the BQ system; ④ results of the water pressure test; and ⑤ state of ground stress. The hierarchical structure model of the ore bodies and the backfill is shown in Figure 12.

The RMR system did not take into consideration water pressure test and elasticity modulus. Currently, there are no uniform domestic and international standards for classifying the results of water pressure tests. Thus, it is necessary to take the existing research into consideration. The principles for classifying the water pressure test evaluation indices were obtained based on field experience, as shown in Table 2.

The fuzzy Delphi analytic hierarchy process (FD-AHP) is a fuzzy group decision-making method combined with fuzzy mathematics appraisal, analytic hierarchy process, and Delphi group decision-making method [21–23]. FD-AHP was used to calculate the total score of rock mass and the backfill; thus, an improved FD-AHP calculation method for the rock mass and the backfill was formed. The basic steps of improved FD-AHP are as follows:(1)Build a hierarchical structure model to calculate the quality of ore body and backfill, as shown in Figure 12.(2)Build an AHP hierarchical structural model.(3)Create a comparison and judgment matrix. For elements on the criterion level, comparison was made one by one. When an element on the above level is taken as the comparison criterion, a comparison scale b_ij can be used to express the corresponding criticality of element i and element j on one of these levels.(4)When the triangular fuzzy number is used to create a group fuzzy judgment matrix, the two-by-two judgment matrix is as follows:wherewhere N is the total number of experts giving scores and is the corresponding criticality of element i and element j judged by expert K:(5)The group fuzzy weight vector is established based on the group fuzzy judgment matrix b_ij, and the corresponding fuzzy weight vector is established using the geometric mean method. The vector is calculated for any evaluation index i (i = 1, 2, … , m): Through further calculation, the following is obtained:where and are the relationship of triangular fuzzy number multiplication and that of the addition, respectively:where(6)In weight decision-making analysis, let the fuzzy weight vector of each index bewhere , , and are the minimum, median, and maximum values, respectively, of the three components of the fuzzy weight vector. The relative weight value of each evaluation index is calculated using the geometrical average approach, with the decision-making weight being obtained through normalization, and the normalized weight vector is

After normalization, the decision-making weight value of each evaluation index is

Thus, it is available to calculate the total score of the classified layers of the ore rock and the backfill based on the evolved FD-AHP, as shown in the following equation:

Before grouting, the quality improvement score of ore bodies was 81, and that after grouting was 88. Before grouting, the quality improvement score of backfill was 65, and that after grouting was 75. Clearly, the grouting effect was apparent. Based on the above analysis, in the collapse area, with grouting reinforcement, the stress state of the surrounding rock structure was improved; thereby improving the self-bearing capacity of the surrounding rock in the collapse area and the stability of the surrounding rocks to be stoped.

5. Conclusion

(1)According to the structure and force characteristics of the surrounding rock in the collapsed area, the theoretical mechanics model of the influence range of the collapse zone is constructed with the loose circle as the boundary. Meanwhile, based on the force method, the deflection line and the corner equation of a statically indeterminate beam-column of the plastic range are obtained, and the analytic expression of the critical slip angle of the beam-column model is obtained. According to rock physical and mechanical experiments, it can be concluded that the angle of the slip line increased along with the uniform force Q, which is inversely proportional to the bending stiffness; the slip angle increases with the influence range of the collapse.(2)Based on the theory of mine pressure, rock formation control, and rock mechanics, the basic collapse form of collapse area is analyzed, the calculation formula of collapse height of ore body is derived, and the collapse height of the upper ore body is estimated to be ≥4.26 m. Therefore, the maximum subsidence of ore body is related to ore body thickness, internal friction angle, density and dilatancy coefficients, width of the original goaf, ore body inclination angle, and the geometric shape of ore body. The obtained calculation formula can be used as a reference for collapse analysis of similar engineering stopes.(3)With the method of numerical analysis, the mechanism of collapse and its catastrophe process are obtained. The influence range of the collapse area is calculated, and the influence range of the collapse area has the form of a “U shape.”(4)The overall grouting reinforcement scheme was determined according to the above analysis results: determine the instability region ⟶ temporary support for laneway ⟶ improve the water environment and the bearing capacity of the surrounding rock ⟶ make mining planning. The safe existing roadways are utilized to layout grouting holes to implement grouting reinforcement in collapse area.(5)Combined with the field water pressure test, and based on the RMR and BQ methods, the evaluation index of water pressure test results was introduced into the classification method. The hierarchical grading model of ore body and backfill was built, and the improved FD-AHP coupling method was established to evaluate the quality of the filling body before and after grouting. The results show that the quality of the ore body and the filling body was improved obviously, and grouting reinforcement improved the bearing capacity of the ore body and the filling body in the mining collapse area.

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 they have no conflicts of interest.

Acknowledgments

This study was supported by the National Key R&D Program of China (No. 2018YFC0809700) and National Natural Science Foundation of China (Nos. 91646101, 71673158, 91324022, and 91646201).

References

R. Gao, J. Zhang, A. J. S. Spearing, M. Li, B. An, and D. Hao, “Research into stope roof control of compound roof by solid backfilling mining,” International Journal of Mining Science and Technology, vol. 26, no. 4, pp. 609–614, 2016.
View at: Publisher Site | Google Scholar
Z.-q. Luo, C.-y. Xie, J.-m. Zhou, N. Jia, X.-m. Liu, and H. Xu, “Numerical analysis of stability for mined-out area in multi-field coupling,” Journal of Central South University, vol. 22, no. 2, pp. 669–675, 2015.
View at: Publisher Site | Google Scholar
J. C. Saul, P. J. Evans, and B. C. Forster, “Fault zone architecture and permeability structure,” Geology, vol. 24, no. 11, pp. 1025–1028, 1996.
View at: Publisher Site | Google Scholar
T. H. Ma, H. Y. Wang, and G. S. Fu, “Present situation analysis of non-coal mine goaf hazards in China,” Journal of Northeastern University, vol. 30, no. S1, pp. 91–95, 2009.
View at: Google Scholar
K. Guan, W. C. Zhu, L. L. Niu, and Q. Y. Wang, “Three-dimensional upper bound limit analysis of supported cavity roof with arbitrary profile in Hoek-Brown rock mass,” Tunnelling and Underground Space Technology, vol. 69, pp. 147–154, 2017.
View at: Publisher Site | Google Scholar
J. L. Peng, J. X. Yang, and M. X. Sun, “Analysis and control on anomaly water inrush in roof of fully-mechanized mining field,” Mining Science and Technology (China), vol. 21, no. 1, pp. 89–92, 2011.
View at: Publisher Site | Google Scholar
R. Hart, P. A. Cundall, and J. Lemos, “Formulation of a three-dimensional distinct element model-Part II. Mechanical calculations for motion and interaction of a system composed of many polyhedral blocks,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 25, no. 3, pp. 117–125, 1988.
View at: Publisher Site | Google Scholar
A. M. Starfield and P. A. Cundall, “Towards a methodology for rock mechanics modelling,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 25, no. 3, pp. 99–106, 1988.
View at: Publisher Site | Google Scholar
D. Mao, B. Nilsen, and M. Lu, “Numerical analysis of rock fall at Hanekleiv road tunnel,” Bulletin of Engineering Geology and the Environment, vol. 71, no. 4, pp. 783–790, 2012.
View at: Publisher Site | Google Scholar
X. L. Yao, D. J. Reddish, and B. N. Whittaker, “Non-linear finite element analysis of surface subsidence arising from inclined seam extraction,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 30, no. 4, pp. 431–441, 1993.
View at: Publisher Site | Google Scholar
W. L. Wood and A. Calver, “Lumped versus distributed mass matrices in the finite elements solution of subsurface flow,” Water Resources Research, vol. 26, no. 5, pp. 819–825, 1990.
View at: Publisher Site | Google Scholar
A. M. Freidin, S. A. Neverov, A. A. Neverov, and P. A. Filippov, “Mine stability with application of sublevel caving schemes,” Journal of Mining Science, vol. 44, no. 1, pp. 82–91, 2008.
View at: Publisher Site | Google Scholar
L. A. Peralta, G. F. Carrier, and C. C. Mow, “An approximate procedure for the solution of a class of transient-wave diffraction problems,” Journal of Applied Mechanics, vol. 33, no. 1, pp. 168-169, 1966.
View at: Publisher Site | Google Scholar
D. F. Malan, “Time-dependent behaviour of deep level tabular excavations in hard rock,” Rock Mechanics and Rock Engineering, vol. 32, no. 2, pp. 123–155, 1999.
View at: Publisher Site | Google Scholar
W. S. Phillips, D. C. Pearson, C. L. Edwards, and B. W. Stump, “Microseismicity induced by a controlled, mine collapse at white pine, Michigan,” International Journal of Rock Mechanics and Mining Sciences, vol. 34, no. 3‐4, pp. 246.e1–246.e14, 1997.
View at: Publisher Site | Google Scholar
L. Y. Tan, H. F. Yu, and L. Chen, “A new approach for predicting bedding separation of roof strata in underground coalmines,” International Journal of Rock Mechanics and Mining Sciences, vol. 61, pp. 183–188, 2013.
View at: Publisher Site | Google Scholar
F. Q. Liu and J. H. Wang, “A generalized slip line solution to the active earth pressure on circular retaining walls,” Computers and Geotechnics, vol. 35, no. 2, pp. 155–164, 2008.
View at: Publisher Site | Google Scholar
F. Q. Liu, J. H. Wang, and L. L. Zhang, “Axi-symmetric active earth pressure obtained by the slip line method with a general tangential stress coefficient,” Computers and Geotechnics, vol. 36, no. 1-2, pp. 352–358, 2009.
View at: Publisher Site | Google Scholar
A. Hackston and E. Rutter, “The mohr-coulomb criterion for intact rock strength and friction - a re-evaluation and consideration of failure under polyaxial stresses,” Solid Earth, vol. 7, no. 2, pp. 493–508, 2016.
View at: Publisher Site | Google Scholar
H. F. Xu, H. S. Geng, W. D. Li et al., “Theory of strength increment of grouting-reinforced bodies for broken rock mass based on BQ,” Chinese Journal of Geotechnical Engineering, vol. 49, pp. 1147–1151, 2014.
View at: Google Scholar
S.-L. Hsueh, “Assessing the effectiveness of community-promoted environmental protection policy by using a Delphi-fuzzy method: a case study on solar power and plain afforestation in Taiwan,” Renewable and Sustainable Energy Reviews, vol. 49, no. 6, pp. 1286–1295, 2015.
View at: Publisher Site | Google Scholar
X. Ying, G.-M. Zeng, G.-Q. Chen, L. Tang, K.-L. Wang, and D.-Y. Huang, “Combining AHP with GIS in synthetic evaluation of eco-environment quality-A case study of Hunan Province, China,” Ecological Modelling, vol. 209, no. 2, pp. 97–109, 2007.
View at: Publisher Site | Google Scholar
S. Lee and K.-K. Seo, “A hybrid multi-criteria decision-making model for a cloud service selection problem using BSC, fuzzy delphi method and fuzzy AHP,” Wireless Personal Communications, vol. 86, no. 1, pp. 57–75, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Jia Nan 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.

PDF Download Citation

Download other formats

Order printed copies