Research Article  Open Access
Ziguang Zhang, Tao Xu, Guangyong Cao, Xiaofeng Liu, "Distribution Features of the Minimum Rock Cover Thickness of the Surrounding Rock SelfStability of the Metro Tunnel in the SoilRock Dualistic Stratum", Advances in Civil Engineering, vol. 2020, Article ID 9120983, 12 pages, 2020. https://doi.org/10.1155/2020/9120983
Distribution Features of the Minimum Rock Cover Thickness of the Surrounding Rock SelfStability of the Metro Tunnel in the SoilRock Dualistic Stratum
Abstract
Based on the maximum value of the selfstability of surrounding rock, the key issues affecting the distribution of the minimum mock cover thickness of the surrounding rock selfstability (MRCTSRS) of the metro tunnel in the soilrock dualistic stratum are discussed. The MRCTSRS was studied for different soft stratum thicknesses Hs and excavation spans D based on cusp catastrophe theory. Moreover, the 3D spatial distribution was determined, and the fitted curves were constructed according to the MRCTSRS, Hs, and D. This facilitated the assessment of the stability of the surrounding rock in the excavation of a metro tunnel. Combined with the construction practices used for the Qingdao metro project, the validity of the obtained MRCTSRS has been verified by the inverse analysis of the monitored data from the metro tunnel. Furthermore, the application of the MRCTSRS in determining the reasonable buried depth of a metro line has been explored. The research results have provided theoretical support and scientific basis for the preparation, revision, and improvement of the relevant codes, standards, and guidelines for metro tunnel planning and design.
1. Introduction
The fundamental goal of an excellent metro tunnel design scheme is to make full use of the selfstability of the surrounding rock. Recently, many subsurface metro tunnels have been carried out in the soilrock dualistic stratum in cities, such as Qingdao [1], Dalian [2], Guangzhou [3], Chongqing [4], and Guiyang [5]. When a metro tunnel is excavated in the urban area where the buried depth of the bedrock is shallow and the soft soil stratum thickness is thin, the tunnel should be placed in the rock stratum and maintained such that the minimum rock cover thickness meets the requirements of the selfstability of the rock surrounding. In this case, the construction of a metro tunnel can greatly reduce the tunnel support measures, save construction costs, reduce the tunnel construction safety risk and construction difficulty, and reduce the impact on the surrounding environment.
Many scholars have studied the relationship between the stability of the surrounding rock and the rock cover thickness. Sun et al. [1] used an example of a metro station in Qingdao to determine the safety factor of a tunnel in surrounding rock at different buried depths by using the strength reduction method. They also considered the safety factor equal to 1 as the basis for judging the surrounding rock selfstability of the metro tunnel and determined the relationship between the minimum rock cover thickness and the soil stratum thickness. Yang and Xiao [6] established the cusp catastrophe model of the safety thickness of a karst tunnel floor and derived the expression for the critical safety thickness. Xu et al. [7] designed a study based on the Qingdao Jiaozhou Bay Subsea Tunnel using two mutually verified discrimination methods: the minimum vertical displacement of the tunnel vault and the maximum safety factor of the tunnel surrounding rock. They then studied the minimum rock cover thickness of the subsea tunnel and considered that the minimum rock cover thickness of a subsea tunnel is the reasonable and effective thickness of the tunnel. Wang et al. [8] pointed out that the reasonable buried depth of the metro station in the soilrock dualistic stratum should be determined by using the double index system of the overburden thicknessspan ratio and the rock cover thicknessspan ratio and noted that the rock cover thicknessspan ratio is more sensitive to the safety of tunnel construction than the overburden thicknessspan ratio. Qiu et al. [9] studied the relationship between the minimum rock cover thickness and the surrounding rock stability of the subsea tunnel by the numerical method. Zhang et al. [10] discussed and expounded upon the issues of the surrounding rock selfstability of the metro tunnel in the soilrock dualistic stratum from the three aspects of the stress, the displacement, and the safety factor. Wang et al. [11] studied the conversion timing of tunnel excavation method in upper soft and lower hard strata based on displacement direction angle theory under different working conditions.
The catastrophe theory was put forward by Thom, a French scholar, in the 1970s [12] and was then expanded into a perfect system by Trotman and Zeeman [13], Zeeman et al. [14], and Poston and Stewart [15]. Catastrophe theory is a mathematical method used to study the phenomenon of discontinuity in nature. The cusp catastrophe model is one of seven primary catastrophe models in catastrophe theory, which has been widely used in the study of tunnel surrounding rock stability, and it has achieved ideal results [16–22]. Xia et al. [16] used cusp catastrophe theory and the discontinuous deformation analysis method to study the stability of the tunnel surrounding rock and obtain the safety factor. Based on the measured deformation data of the surrounding rock, Ren et al. [17] established the catastrophe model of the tunnel surrounding rock cusp by using catastrophe theory and deduced the instability criterion of the surrounding rock. Wang et al. [18] used the cusp catastrophe model to study the stability of gas storage pillars in layered salt caverns. Zhang and Han [19] used catastrophe theory to study the collapse mechanism and possible collapse block shape of a shallow, unlined tunnel. Zhou [20] et al. established a model based on the cusp catastrophe theory and analyzed the stability of the support system in the goaf of room pillar gypsum. Zhang et al. [21] and Huang and Zhang [22] used catastrophe theory to study the instability mechanism of a shallow tunnel.
The problem of the minimum mock cover thickness of the surrounding rock selfstability (MRCTSRS) of the tunnel has been studied by many scholars, who have used numerical calculations to analyze the stress, the deformation, and the safety factor. However, based on cusp catastrophe theory, rarely no relevant report has been found on the MRCTSRS of the metro tunnel in the soilrock dualistic stratum. This paper based on the calculation and analysis of the MRCTSRS under different soft stratum thickness Hs and excavation span D was studied based on cusp catastrophe theory and numerical calculation, the 3D spatial distribution was drawn, and the fitted curves was constructed among MRCTSRS, Hs, and D, which facilitated the assessment on the stability of the surrounding rock in excavation of metro tunnel. Combined with the construction practice of Qingdao metro project, the validity of the obtained MRCTSRS has been verified by the inverse analysis on the data monitored for metro tunnel, and the application of MRCTSRS in determining the reasonable buried depth of metro line has been expounded. The research results can provide theoretical support for the design of reasonable depth of the vertical section of metro tunnel in soilrock dualistic stratum and further offer a scientific basis for the preparation, revision, and improvement of relevant codes, standards, and guidelines.
2. Methodology and Research Process
2.1. Cusp Catastrophe Theory
The cusp catastrophe model uses the simplest catastrophe mode with lag, divergence, and other behaviors, which is unique to the highorder catastrophe model. The critical interface is easy to construct and has a strong geometric rationale. Of the seven primary catastrophe models used in catastrophe theory, the cusp catastrophe model is the one with the most extensive application [23–29].
The potential function V of the cusp catastrophe model is a 3D space (x, u, ) composed of the control variables u and and the state variable x as follows:
An equation for the cusp catastrophe model for the balanced curved surface has been obtained after first deriving the potential function (1) with respect to x as follows:
The cusp catastrophe model uses a continuous balanced surface with a 3D phase space, the interior of which is folded into upper, middle, and lower pages, as shown in Figure 1. Different areas have different numbers of balance points. On the middle page, the balance point falls on the maximum value of the potential function, so it is unstable. However, on the upper and lower pages, the potential function takes its minimum value, so the balance point is stable. The boundary between the upper and lower pages, on the one hand, and the middle page, on the other hand, is composed of points with vertical tangents, which is a set of mutation points as follows:
Combining formulas (2) and (3) results in a bifurcation set, which is the projection of the fold of the balancecurved surface in the uv plane as follows:
When applying catastrophe theory, scholars have used formula (4) as a direct criterion to determine the cusp mutation. When Δ > 0, the system is in a stable state; when Δ = 0, the system is in a critical equilibrium state; when Δ < 0, the system is in an unstable state [6].
2.2. Analysis Model
See Figure 2, for the analysis model of the MRCTSRS of the metro tunnel in the soilrock dualistic stratum. In the figure, H represents the buried tunnel depth; D, the excavation span; h, the excavation height; and Hs and Hr, respectively, the soft soil stratum thickness and the rock cover thickness.
A statistical analysis was performed of the geotechnical test results of 3227 geological boreholes in nearly 2000 sets of soilrock strata along the Qingdao metro tunnel. Accordingly, the values of the physical and mechanical parameters of the soft soil stratum and the rock stratum in the soilrock dualistic stratum structure are shown in Table 1.

2.3. Target
The vault is the point at which the deflection is the largest around the cavern of the shallow metro tunnel. Therefore, the vertical displacement of the vault is the key indicator when evaluating the stability of the surrounding rock of the tunnel [30, 31]. If the rock cover thickness is large over the vault of the tunnel in the soilrock dualistic stratum and it is gradually reduced, the vertical displacement of the vault at first is slowly reduced. Then, the displacement gradually increases after reaching its minimum as the thickness continues to decrease. When the thickness decreases to a certain value, the displacement sharply rises [7]. Currently, the vertical displacement of the tunnel vault is selected as the target for evaluating the selfstability of the surrounding rock of the tunnel. Furthermore, the corresponding rock cover thickness at the mutation point is used as the MRCTSRS of the tunnel.
2.4. Process
For the MRCTSRS of the metro tunnel in the soilrock dualistic stratum based on featured point displacement in cusp catastrophe theory, the calculation is performed as follows:(1)First, a numerical calculation model is constructed of the metro tunnel in the soilrock dualistic stratum. In the numerical calculation, the stratum is considered according to the plane strain of ideal elasticplastic materials, which should conform to the Mohr–Coulomb strength criterion and the large deformation hypothesis. ANSYS numerical calculation software and the inscribed circle DP4 yielded criterion based on the Mohr–Coulomb conditions, which were adopted for the onetime excavation of the whole section without considering the influence of groundwater and other factors. The upper boundary of the calculation model is set at the ground surface, and the lower boundary is set from the vault bottom of the excavation section down to 3.5 times the excavation height. Moreover, the distance between the left and right boundaries is more than 7 times the excavation span. The initial geostress incorporates not only the gravity stress of the rock and soil but also a uniformly distributed load of 20 kPa imposed on the ground surface. A large deformation mode is used during the calculation. The failure of the surrounding rock by large deformation can be reproduced by updating the coordinates. The vertical and horizontal ranges of the unit grid around the tunnel are 0.5 m, and those of the more peripheral unit grid are set at 1.0 m.(2)By continuous calculations, the rock cover thickness is determined with respect to the smallest vault vertical displacement of the tunnel surrounding rock as the thickness gradually increases. Then, the thickness is reduced gradually at a specific interval (e.g., 1.0 m or 0.5 m) for multiple calculations. The vertical displacements Si of the tunnel vault in each calculation are recorded to derive a sequential SiHr curve between the displacement Si and the rock cover thickness Hr. The number of calculations should increase in the mutation part of the SiHr curve.(3)The calculation should begin starting with the rock cover thickness with respect to the smallest vault vertical displacement of the tunnel surrounding rock as the thickness gradually increases. The end interpolation method is used to fit the data of vertical displacement S and overburden rock thickness Hr, producing a quartic function as follows: where are undetermined coefficients.(4)The fitting equation is transformed, and the standard potential function of the cusp mutation is derived as , for which and .(5)The bifurcation set equation is derived as according to cusp catastrophe theory.
In this research, Hr with respect to is selected as the MRCTSRS. Its calculation is shown in Figure 3.
2.5. Example Analysis
The MRCTSRS calculation is demonstrated in detail in the following example of a metro tunnel with soft soil stratum thickness Hs of 12 m. The shape and dimensions of the tunnel section are shown in Figure 4. The calculation results of the vertical displacement Si under different Hr values are shown in Table 2.

The data of S and Hr are used, so Hr is set to 11 m to perform the polynomial fitting. When Hr is 1 m, the fitting is carried out on 11 sets of data with Hr vales from 1 m to 11 m, producing the result as follows:
Substitute the values of the four parameters of ∼ from formula (6) into formulas (1) and (4). The result is and .
, indicating that Hr does not meet the MRCTSRS requirements.
When Hr is 2 m, , indicating that it meets the requirements. Therefore, the MRCTSRS should be between 1 m and 2 m.
To improve further the calculation accuracy of the MRCTSRS, Hr is set as 1.5 m. The result shows , indicating that the Hr meets the requirements; when Hr is 1.4, , indicating that it does not meet the requirements. Therefore, the MRCTSRS should be between 1.4 m and 1.5 m. The MRCTSRS is finally set at 1.45 m when the soft soil stratum thickness is 12 m in this paper.
3. Distribution Features of the MRCTSRS
3.1. MRCTSRSs of the Metro Station Tunnel and the Running Tunnel
The example of the engineering of the Qingdao metro tunnel is used in this research to investigate the MRCTSRS of the station tunnel and the running tunnel with two main types of sections in a soilrock dualistic stratum (see Figure 5).
The results of the MRCTSRS calculations for the metro station tunnel and the running tunnel under different soft soil stratum thicknesses Hs are shown in Table 3 and Figure 6. The MRCTSRS of the metro tunnel approximately linearly increases or decreases with respect to Hs. The linearity fitting is used to predict the MRCTSRSs of the metro station and the running tunnel. The fitted equation for the station tunnel MRCTSRS and Hs is shown in formula (7). The adjusted Rsquare value is 0.9991, and the MRCTSRS is approximately linear along with Hs. The fitted equation of the running tunnel MRCTSRS and Hs is shown in formula (8), and the adjusted Rsquare value is 0.99258. With formulas (7) and (8), the MRCTSRSs for the station tunnel and the running tunnel can be derived with respect to different Hs.

(a)
(b)
The MRCTSRS of the metro tunnel increases approximately linearly with the increase in Hs. This result occurs because, when excavating the metro tunnel in the soilrock dualistic stratum under the requirements of the MRCTSRS, the overlying soft soil stratum bears a higher load, while the rock cover stratum bears more of the structure. The rock mechanics performance around the tunnel is similar to the force mode of a composite arch or a composite beam. Therefore, when the Hs is greater, the corresponding load and the MRCTSRS become larger.
3.2. Discussion of How MRCTSRS Works for the Metro Tunnel
The stability characteristics of the tunnel surrounding rock are closely related to the shape and size of the tunnel section. By focusing on the metro station tunnel, the general rule of how the MRCTSRS varies with the excavation span D is discussed in the present study while maintaining the original section shape and highspan ratio but using only the scaling tunnel section size. The tunnel section size is shown in Figure 7. The MRCTSRS calculation results of the metro tunnel are shown in Table 4 when excavation spans D are 10.4 m and 31.2 m with respect to Hs.

The 3D spatial distribution of MRCTSRS ∼ Hs ∼ D among the metro tunnel MRCTSRS, soft soil stratum thickness Hs, and excavation span D is shown in Figure 8. The figure overly shows the distribution characteristics of the metro tunnel MRCTSRS with respect to different Hs and D. At the same time, according to the figure, the surrounding rock stability and supplementary measures in the excavation for the metro tunnel in the soilrock dualistic stratum can be determined as follows:(1)When Hr falls below the interface plane, the excavation for the surrounding rock cannot meet the selfstability condition. Strong supplementary engineering measures are needed to prevent its instability to ensure normal tunnel excavation.(2)When Hr falls on the interface plane, the excavation of the surrounding rock just meets the selfstability condition. Some supplementary engineering measures are needed to prevent its instability to ensure normal tunnel excavation.(3)When Hr falls above the interface plane, the excavation for the surrounding rock can meet the selfstability condition. Only simple supplementary engineering measures are needed to prevent the geometric instability of rocks to ensure normal tunnel excavation.
The fitting equation of MRCTSRS ∼ Hs ∼ Dis shown as formula (9), and the adjusted Rsquare is 0.99597. By this formula, the MRCTSRS with respect to different Hs and D can be derived:
4. Case Analysis and Engineering Application
4.1. Project Overview
4.1.1. Project Descriptions
The planned metro lines for Qingdao’s rail transit are approximately 807 km long and occur in 16 lines. Qingdao Metro Line 3, the first urban rail transit construction, is approximately 24.8 km long and includes 22 stations. Its construction is followed by Line 2, covering a length of approximately 25.2 km with 22 stations. The two lines have 16 underground excavated stations. Except for Dunhua Road Station in the towercolumn structure, the other excavated stations all adopt a singlearch largespan structure with an excavation span of 19.2 m∼23.8 m (see Figure 9 and Table 5). The metro running tunnel mainly adopts the horseshoeshaped subterranean excavation form with a single span of 5.8 m∼6.2 m for a oneway single lane.

4.1.2. Stratigraphic Features
The stratum distribution along Qingdao Metro from top to bottom includes the Quaternary stratum as well as the strongly weathered, the moderately weathered, and the microweathered (or nonweathered) granite strata. To simplify the analysis and highlight the main points, the Quaternary stratum and the strong weathering stratum were considered soft soil stratum, and the strata in the middle weathering strata and below are considered rock stratum, thus constructing a distribution model for the soilrock dualistic stratum. The structural characteristics of the dualistic stratum along the two metro lines in Qingdao are shown in Figure 10.
According to the statistical analysis of 3227 geological borehole samples along the two metro lines, the average thickness of the soft soil stratum is 10.2 m; the maximum is 38.0 m; the median is 10.0 m; and the lower and upper quartiles are, respectively, 5.0 m and 14.8 m (as shown in Figure 11 and Table 6). Further analysis shows that the number of borehole samples with soft soil stratum thicknesses of 9 m or less along the metro lines accounts for 45.6% of the samples; those with less than 12 m, 61.4%; those with less than 15 m, 77.5%; and those with less than 18 m, 89.7% (as shown in Table 7).


4.2. Analysis of Metro Tunnel Monitoring Data
4.2.1. Data Resource
The monitoring data come from six station tunnels of the two metro lines, which have approximately the same section size, excavation depth and technique, support measures, etc. The number of monitoring points, the soft soil stratum thickness Hs, the rock cover thickness Hr, and the MRCTSRS are shown in Table 8. The statistical data are from a third party. The monitoring points for the ground surface settlement are selected within the main area of the tunnel, among which the abnormal points should be excluded in the statistics.

4.3. Inversion Analysis
Table 9 shows that, with regard to the overburden thicknesses exceeding the MRCTSRS requirements for the stations of Zhiquan Road, Haichuan Road, and Junfeng Road, the average surface settlement is between 2.29 mm and 6.37 mm, with a maximum value between 6.22 mm and 12.59 mm, which are small in general. Regarding those for the stations of Zhongshan Park, Jiangxi Road, and Zhanshan Road, the average surface settlement is between 37.5 mm and 54.75 mm, with a maximum value between 71.54 mm and 80.96 mm, which are large in general.

Surface settlement can comprehensively reflect the effect of metro construction quality control and the degree of construction influence on the surrounding rock as well as the soil and neighboring environment. For a metro tunnel with shallow buried depth under the approximately same conditions of section size, excavation depth and technique, support measures, etc., the rock cover thickness on the vault meets the MRCTSRS requirements or not, the surface settlement caused by the tunnel excavation is different by an order of magnitude. This situation fully demonstrates that the selfstability of the surrounding rock has a good correlation with the rock cover thickness (as shown in Figure 12). The rationality of the conclusions obtained in this study has been fully verified.
4.4. Engineering Application
The fundamental goal of an excellent metro line design is to make full use of the selfstability of the surrounding rock and to minimize supplementary engineering measures. Costs should be considered only when engineering security is guaranteed. Furthermore, geological risks and construction difficulty should be reduced. The degree to which the surrounding rock can be selfstabilized depends mainly on the distribution conditions of the strata along the metro lines, the understanding of these conditions, and the ability to execute the required tasks. Therefore, the Geologically Easily Stable Area (GSEA) in this study refers to the area that meets the MRCTSRS requirements for the excavation of the tunnel along the metro in the soilrock dualistic stratum, otherwise referred to as the Geologically Difficultly Stable Area (GSDA).
The metro station and the running tunnel are considered, respectively, as a singlearch largespan excavation structure with the excavation span of 20.8 m and a horseshoeshaped structure with the excavation span of 6.2 m. The buried depth is set at equal intervals of 3 m falling within the scope of 6 m∼30 m. Table 10 and Figure 13 show the calculation results of the proportions of the GSEA lengths of the station (considering only 16 subterranean excavation stations) and the running tunnels, respectively, as well as the lengths of the station and the running tunnel of Qingdao’s two metro lines under different buried depths.

As the buried depth increases, the proportion of the GSEA in the station tunnel rises first rapidly but then slowly. When the buried depth of the station tunnel is at 6 m, the proportion is 49.2%. When the depth increases to 12 m, the proportion increases rapidly to 76.4%. However, when the depth reaches 18 m, the proportion increase to 90.7%. The change of GSEA with buried depth in running tunnel is similar to that in station tunnel. When the buried depth of the running tunnel is 6 m, the proportion is only 26.3%. However, when the depth increases to 15 m, the proportion rapidly increases to 71.6%. As the depth continues to rise to 18 m, the proportion increases to 87.5%. When the depth reaches 21 m, the proportion slowly increases to 93.9%.
The metro line is mainly composed of the station and the running tunnel. The buried depth of the line depends on the depth of the station tunnel. The depth of the running tunnel is mainly determined by factors such as the buried depth of the station tunnel and the longitudinal slope. The height difference between the station and the running tunnels is set as 12 m. Therefore, when the buried depth of the Qingdao metro station falls within 12m∼15 m, the GSEA of the station tunnel accounts for 76.4%∼84.1%, and the GSEA of the running tunnel accounts for 96.6%∼98.5%.
5. Conclusion
(1)The MRCTSRS of the metro tunnel approximately linearly increases or decreases with respect to the soft stratum thickness Hs in the soilrock dualistic stratum. When the Hs is larger, so is the corresponding MRCTSRS.(2)The 3D spatial distribution was determined, and the fitted curves were constructed according to the MRCTSRS, Hs, and D. This facilitated the assessment of the stability of the surrounding rock in the excavation of a metro tunnel.(3)The variation rule of the Geologically Easily Stable Area with respect to buried depth has been revealed by examining the distribution features of the strata along the Qingdao metro lines. Moreover, the application of the MRCTSRS in engineering has been expanded to determine the reasonable buried depth of the metro line.
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.
Authors’ Contributions
Ziguang Zhang was in charge of the paper and completed the construction of thesis framework system and thesis writing. Tao Xu and Xiaofeng Liu were responsible for the collection, sorting, and analysis of relevant project case data in this paper . Guangyong Cao completed the related work of numerical calculation in this paper.
Acknowledgments
This study was funded by the Scientific Research Project of Anhui Jianzhu University (no. 2019QDZ24) and Science and Technology Support Project of Anhui Sijian Holding Group Co., Ltd. (no. HYB20190126).
Supplementary Materials
Figures and Tables of Survey Data in Measurements. (Supplementary Materials)
References
 K. G. Sun, W. P. Xu, W. G. Qiu, H. B. Li, Q. Y. Xianand, and T. Li, “Study on the characteristics of safety distribution changing with buried depth for metro station in uppersoft and lowerhard stratum,” Advances in Civil Engineering, vol. 2018, Article ID 6047919, 14 pages, 2018. View at: Publisher Site  Google Scholar
 Y. L. Shang, S. J. Du, T. Y. Han, and B. Shao, “Case study on deformation control of uppersoft and lowerhard large span tunnel station using combined control technology and monitoring demonstration,” Sains Malaysiana, vol. 46, no. 11, pp. 2091–2099, 2017. View at: Publisher Site  Google Scholar
 T. J. Liu, S. Y. Zheng, X. W. Tang, and W. X. Xu, “Experimental study on seismic response of a largespan and columnfree subway station in composite strata,” Shock and Vibration, vol. 2019, Article ID 8412390, 15 pages, 2019. View at: Publisher Site  Google Scholar
 R. H. Lv, “On adaptability assesment of a TBM used in urban subways,” Modern Tunnelling Technology, vol. 54, no. 1, pp. 31–39, 2017. View at: Google Scholar
 J. B. Wang, Q. Huo, Z. P. Song, and Y. W. Zhang, “Study on adaptability of primary support arch cover method for largespan embedded tunnels in the uppersoft lowerhard stratum,” Advances in Mechanical Engineering, vol. 11, no. 1, pp. 1–15, 2019. View at: Publisher Site  Google Scholar
 X. L. Yang and H. B. Xiao, “Safety thickness analysis of tunnel floor in karst region based on catastrophe theory,” Journal of Central South University, vol. 23, no. 9, pp. 2361–2367, 2016. View at: Publisher Site  Google Scholar
 B. S. Xu, S. C. Li, R. C. Liu, and C. L. Zhao, “Study on the reasonable cover thickness of a subsea tunnel with the numerical calculation criterion method,” in Proceedings of the 4th geoChina International Conference: Emerging Technologies in Tunnel Engineering, Modeling, Design, Construction, Repair, and Rehabilitation, pp. 81–91, Shandong, China, July 2016. View at: Publisher Site  Google Scholar
 X. D. Wang, Y. Yuan, X. Wu, and Z. Du, “Study on method to determining the ration depth of metro station in composed of soil and weathered rock stratum,” in Proceeding of the 6th European Congress on Computational Methods in Applied Sciences and Engineering, Vienna, Austria, September 2012. View at: Google Scholar
 X. B. Qiu, S. C. Li, and W. Z. Chen, “Numerical analysis of rock cover thickness of subsea tunnel,” Tunnelling and Underground Space Technology, vol. 19, p. 426, 2004. View at: Google Scholar
 Z. G. Zhang, K. G. Sun, and W. G. Qiu, “Reasonable thickness of metro station with single arch largespan in uppersoft and lowerhard stratum,” International Conference on Mechanics and Architectural Design, vol. 2016, no. 5, pp. 72–82, 2016. View at: Google Scholar
 W. Wang, M. J. Hu, H. W. Zhang, and X. M. Zhang, “Conversion timing of tunnel excavation methods in uppersoft and lowerhard stratum based on displacement direction angle theory,” International Journal of Geomechanics, vol. 19, no. 4, Article ID 05019003, 2019. View at: Publisher Site  Google Scholar
 R. Thom and D. H. Fowler, Structural Stability and Morphogenesis, BenjaminAddison Wesley, New York, NY, USA, 1975.
 J. M. T. Thompson and E. C. Zeeman, “Classification of elementary catastrophes of codimension ≤ 5, structural stability , the theory of catastrophes and applications in the sciences,” Lecture Notes in Mathematics, vol. 525, SprigerVerlag, Berlin, Germany, 1976. View at: Google Scholar
 E. C. Zeeman, “Bifurcation,catastrophes and turbulence,” New Directions in Applied Mathematics, SprigerVerlag, New York, NY, USA, 1982. View at: Google Scholar
 T. Poston and I. Stewart, Catastrophe Theory and its Application, Pitman, London, UK, 1978.
 C. C. Xia, C. B. Xu, and X. Zhao, “Study on the strength reduction DDA method and its application to mountain tunnel,” International Journal of Computational Methods, vol. 9, no. 3, Article ID 1250041, 2012. View at: Publisher Site  Google Scholar
 S. Ren, Z. Wang, D. Y. Jiang, and C. Liu, “Study on the catastrophe model of the surrounding rock and simulating the constructing process by DEM in gonghe tunnel,” in Proceedings of the 2013 Third International Conference on Intelligent System Design and Engineering Applications, pp. 1355–1357, Hong Kong, China, January 2013. View at: Publisher Site  Google Scholar
 T. T. Wang, X. Z. Yan, H. L. Yang, and X. J. Yang, “Stability analysis of the pillars between bedded salt cavern gas storages by cusp catastrophe model,” Science China.Technological Sciences, vol. 54, no. 6, pp. 1615–1623, 2011. View at: Publisher Site  Google Scholar
 C. P. Zhang and K. H. Han, “Collapsed shape of shallow unlined tunnels based onFunctional catastrophe theory,” Mathematical Problems in Engineering, vol. 2015, Article ID 681257, 13 pages, 2015. View at: Publisher Site  Google Scholar
 Y. J. Zhou, X. D. Xu, X. T. Li, M. Li, and Y. G. Yang, “Study on catastrophe instability of support system in gypsum goaf based on energy dissipation theory,” Advances in Civil Engineering, vol. 2018, Article ID 4293584, 9 pages, 2018. View at: Publisher Site  Google Scholar
 R. Zhang, H. B. Xiao, and W. T. Li, “Functional catastrophe analysis of collapse mechanism for shallow tunnels with considering settlement,” Mathematical Problems in Engineering, vol. 2016, Article ID 4820716, 11 pages, 2016. View at: Publisher Site  Google Scholar
 X. L. Huang and R. Zhang, “Catastrophe stability analysis for shallow tunnels considering settlement,” Journal of Central South University, vol. 25, no. 4, pp. 946–960, 2018. View at: Publisher Site  Google Scholar
 M. A. Ghorbani, R. Khatibi, B. Sivakumar, and L. Cobb, “Study of discontinuities in hydrological data using catastrophe theory,” Hydrological Sciences Journal—Journal des Sciences Hydrologiques, vol. 57, no. 7, pp. 1137–1151, 2010. View at: Publisher Site  Google Scholar
 X. P. Zheng, J. H. Sun, and Y. Cheng, “Analysis of crowd jam in public buildings based on cuspcatastrophe theory,” Building and Environment, vol. 45, no. 8, pp. 1755–1761, 2010. View at: Publisher Site  Google Scholar
 Y. Tao, J. Cao, J. M. Hu, and Z. C. Dai, “A cusp catastrophe model of mid–longterm landslide evolution over low latitude highlands of China,” Geomorphology, vol. 187, pp. 80–85, 2013. View at: Publisher Site  Google Scholar
 B. Hu and N. Xia, “Cusp catastrophe model for sudden changes in a person’s behavior,” Information Sciences, vol. 294, pp. 489–512, 2015. View at: Publisher Site  Google Scholar
 J. C. Xu and Y. D. Ni, “Prediction of greycatastrophe destabilization time of a granite residual soil slope under rainfall,” Bulletin of Engineering Geology and the Environment, vol. 78, no. 8, pp. 5687–5693, 2019. View at: Publisher Site  Google Scholar
 Y. M. Tian, B. Zheng, H. L. Shen, S. N. Zhang, and Y. R. Wang, “A novel index based on the cusp catastrophe theory for predicting harmful algae blooms,” Ecological Indicators, vol. 102, pp. 746–751, 2019. View at: Publisher Site  Google Scholar
 S. W. Ma, Z. Q. Luo, J. H. Hu, Q. F. Ren, Y. G. Qin, and L. Wen, “Determination of intervening pillar thickness based on the cusp catastrophe model,” Advances in Civil Engineering, vol. 2019, Article ID 8253589, 11 pages, 2019. View at: Publisher Site  Google Scholar
 J. L. Qiu, Y. L. Xie, and H. B. Fan, “Centrifuge modelling of twintunnelling induced ground movements in loess strata,” Arablan Journal of Geosciences, vol. 10, no. 22, pp. 12517–13017, 2017. View at: Publisher Site  Google Scholar
 J. X. Lai, X. L. Wang, J. L. Qiu, J. X. Chen, Z. N. Hu, and H. Wang, “Extreme deformation characteristics and countermeasures for a tunnel in difficult grounds in southern Shaanxi, China,” Environmental Earth Sciences, vol. 77, no. 19, p. 14, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Ziguang Zhang 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.