Influence of shield tunnel construction on ground surface settlement under the condition of upper-soft and lower-hard composite strata

Taking the shield tunnel project of Guangzhou Metro Line 8 from Tongdewei Station to Shangbu Station as the research background, using the research method of finite element simulation and site monitoring, this paper analyses the influence rules of shield tunneling on ground subsidence under the condition of different hard rock height ratios. The research results show that in the process of crossing different hard rock height ratio composite stratum, as the hard rock height ratio decreases, the value of ground settlement decreases and settlement tank becomes shallow. The surface subsidence in different hard rock height ratio strata is obviously different, and the maximum difference is about 8.6 mm; The influence of the hard rock height ratio on the surface longitudinal settlement is mainly reflected in the position change of the beginning and the end of the settlement. With the increase of the hard rock height ratio, the shield construction reduces the amount of the surface longitudinal settlement and its influence range; Through the research, it is found that the hard rock height ratio in the 0-0.2 and 0.5-1 is the sensitive interval, and the settlement value in these two hard rock height ratio interval varies greatly. It is necessary to pay attention to the uneven settlement of the ground surface caused by shield construction in the sensitive hard rock height ratio interval. The research results of this paper can provide reference value for similar shield construction in upper soft and lower hard composite stratum.


Introduction
Traffic construction in major cities has increasingly become prominent with rapid urbanization. Urban rail transit construction (dominated by subways) has been conducted in the infrastructure of large cities in China in recent years to alleviate the increasingly heavy traffic problems. Shield tunneling is widely used in urban underground rail transit construction because of its safety and small impact on the stratum and surrounding environment. However, stratum disturbance is inevitable during shield construction, which causes the uplift and subsidence of disturbed soil and even cracking of buildings or damage of underground pipelines [1][2][3][4][5]. Experts and scholars at home and abroad have conducted extensive research on the impact of shield construction [6][7][8][9][10]. Stratum deformation caused by the construction of urban metro tunnels was reported in 1960s. In 1969, Peck [11] proposed a Gauss equation to predict the stratum deformation caused by shield tunneling. Subsequently, scholars at home and abroad have performed considerable theoretical research, numerical simulation, and model testing related to this problem [12][13][14][15][16]. Liu [17] established the original Peck formula. On this basis, consolidation settlement caused by the dissipation of excess pore water pressure was considered, and a formula for estimating the surface vertical and horizontal settlements through formation loss was presented. In practical engineering, the factors affecting stratum subsidence are extremely complex, and subsidence trough and settlement curve of stratum change accordingly [18][19][20][21].
China has a vast territory and varying geological conditions. Its upper-soft and lower-hard strata are widely distributed in the coastal areas of Guangzhou, Shenzhen, and Fuzhou and they are typical geological conditions in urban rail transit. The engineering characteristics of upper-soft and lower-hard strata have huge difference in physical and mechanical properties between rock and soil and sensitive response to external disturbances [22][23][24]. During the construction of metro shield in such special strata, the upper-soft soil can be easily over excavated, causing immense stratum disturbance; moreover, this soil cannot form a stable pressure arch effect [25][26][27][28], which will cause abnormal movement of the excavated surface soil, large-scale collapse of the soft soil stratum, deviation of the shield driving angle from the predetermined value, and the occurrence of surface subsidence. Large-scale settlement remarkably challenges shield tunnel construction. Therefore, the influence of shield tunneling on strata under the condition of upper-soft and lower-hard composite strata should be investigated [29][30][31][32][33]. Deng et al. [34] took Shenzhen Metro Line 2 as an engineering example, conducted a construction risk analysis under the conditions of upper-soft and lower-hard composite strata, and obtained the specific tunneling parameters through slag soil improvement and tool change with pressure; consequently, they provided valuable data for the theoretical study of shield tunneling under complex geological conditions. Wang et al. [35] conducted indoor model testing and numerical simulation to evaluate the over-excavation phenomenon of soft soil and disturbance characteristics in the upper-soft and lower-hard composite strata of EPB shield. The study provided certain reference value for the construction tools of EPB shield in the upper-soft and lower-hard strata. Wuke et al. [36] took Shenzen Metro Line 7 as an engineering example and established a three-dimensional finite element model of the existing urban road crossing under the condition of upper-soft and lower-hard strata. The numerical simulation results were compared with the field monitoring and measurement results, and the deformation and failure law of the existing urban road crossing under the condition of upper-soft and lower-hard strata were analyzed. The optimum construction method and corresponding support were also summarized. Yang et al. [37] evaluated the influence of upper-soft and lower-hard strata on the stability of tunnel surrounding rock. Physical experiments and numerical simulations were conducted for tunnel excavation in inclined upper-soft and hard strata. Moreover, the evolution law and distribution characteristics of stress, displacement, and failure of surrounding rock were analyzed. Wang et al. [38] performed numerical simulation to simulate the main tunnel construction in the upper-soft and lower-hard strata of a provincial hospital station in Guiyang Metro Line 2. The advantages of main support arch cover in settlement control, plastic zone distribution, and support structure safety were analyzed, and the underlying surrounding rock grade range suitable for the main support arch cover was presented. The aforementioned studies revealed the mechanism of shield tunneling in the upper-soft and lower-hard strata, which provide a certain reference value for follow-up studies. However, the research results have certain limitations due to the difference of geological conditions between the upper-soft and lower-hard composite strata and the complexity of the actual construction. The present study divides the different hard rock height ratios along the upper-soft and lower-hard composite strata of the shield tunnel project of Tongdewei-Shangbu Station of Guangzhou Metro Line 8. The influence law of shield tunnel construction on ground settlement under different hard rock height ratios is obtained through statistical analysis of field settlement data and finite element numerical simulation analysis. Some suggestions are presented on the basis of the actual construction requirements, and the research results can serve as reference for similar projects.

Prediction theory of surface transverse settlement
Peck believed that ground subsidence caused by construction occurs without drainage; thus, the volume of settlement trough should be equal to the volume of stratum loss. On this basis and combined with a large number of measured data, Peck proposed a method for estimating ground subsidence caused by shield tunneling in the construction stage. This method assumes that ground loss is uniformly distributed in tunnel length, the lateral distribution of ground subsidence is similar to a normal distribution curve, and its prediction method is expressed as follows [11]: where is the ground settlement value at the central axis of the tunnel; is the maximum ground settlement value at the central line of the tunnel; is the width coefficient of the ground settlement trough, that is, the distance between the central point of the tunnel and the reverse bending point of the settlement curve; is the bottom layer loss of the unit length of the tunnel caused by excavation, that is, the volume difference between the actual excavation volume of tunnel construction and the volume after tunnel completion; is the loss rate of stratum volume; is the calculated tunnel radius; denotes the tunnel depth; and is the friction angle in soil. As shown in Fig. 1, the ground loss is the difference between the actual excavated soil volume and the completed tunnel volume in shield tunneling. The volume of completed tunnel includes the volume of pressed slurry wrapped around the tunnel. The formation loss rate is expressed as , which is a percentage of the theoretical discharge volume of shield tunneling. The theoretical discharge volume of the circular shield is , where is the advancing length; and the formation loss per unit length is . The surrounding soils move and cause land subsidence to compensate for the stratum loss. As shown in the figure, the area of the settlement trough is equal to formation loss per unit length of the tunnel.

Prediction theory of longitudinal surface settlement
The ground settlement curve along the arch roof of the longitudinal tunnel is shown in Fig. 2.
The estimation formula of settlement at a certain point is expressed as follows [11]: where is the settlement (m) distributed along the arch roof of a longitudinal tunnel, in which negative and positive values denote the uplift and settlement values, respectively; is the distance from the settlement point to origin o of the coordinate axis (m); is the distance (m) from the excavation surface of the shield to origin of the coordinate axis at the starting point of shield driving; is the distance (m) between the shield excavation surface and the coordinate axis origin o; is the shield length (m); ′ − ; ′ − ; represents the ground loss caused by shield excavation (m 8 /m; for example, under-excavation of shield causes negative formation loss, and is negative); and is the formation loss (m 8 /m). Behind the shield excavation face, the formation loss (m8/m) caused by all construction factors, such as insufficient grouting in the shield tail clearance and the change of shield driving direction, is always positive.
When is negative, the results calculated by Eq.
(2) are expressed as follows. In front of the shield tail, a certain range of points on the ground surface uplift shows an uplift curve and becomes a settlement curve in the back of the uplift curve (Fig. 2). Eq. (2) is basically applicable to the normal and uniform construction conditions of shield front and tail ground loss. If the shield front under-excavation varies uncertainly, then the Grouting is delayed, the grouting volume varies uncertainly, and the actual settlement curve does not conform to the calculation curve.

Project survey
The north extension project of Guangzhou Metro Line 8 (Cultural Park-Baiyun Lake) is an underground line that is approximately 16.3 km long. The shield tunneling section between Tongdewei Station and Shangbu Station is located in the middle of the entire line. The section line is spread from south to north along the east side of Xichao Road viaduct after Tongdewei Station. The line is mainly located along the lower side of Xichao Road. Left-line starting and closing mileage is ZDK21+078.150-ZDK21+694.400, right-line starting and closing mileage is YDK21+078.150-YDK21+694.400, and the interval length is 616.25 m. The roof depth of the tunnel structure is approximately 10-16 m, and the elevation is -2.40-7.40 m. The inner diameter of the tunnel is 5400 mm, and the thickness of the segment is 300 mm. The entire ring of the shield segment is divided into six parts, namely, three standard A-shaped segments, two adjacent C-shaped segments, and one K-shaped segment. Shield tunneling is used in this section.

Engineering geological conditions
On the basis of the stratum data presented by boreholes, such as geological age, genetic type, lithological characteristics, weathering degree, and other engineering characteristics of strata, the rock and soil layers along the line are divided into four categories, namely, artificial filling layer, alluvial-diluvial layer, residual layer, and rock weathering zone. The sublayers in each layer are subdivided on the basis of their geological conditions. In this area, the hard rock surface immensely fluctuates, and the soft strata above the rock surface are thick. The main soil and rock layers in the small spacing section of shield tunneling are <1> artificial fill, <4-2> silty soil, <3-1> silty sand, <3-2> medium-coarse sand, <3-3> gravel sand, and <9C-2> slightly weathered limestone. Most of the tunnel bodies have grade VI surrounding rocks and sand, and soft soil is developed above the roof. The local floor has slightly weathered limestone intruding into the tunnel. The strength of slightly weathered limestone is high, and the tunnel section is soft and hard. For shield tunneling, the lower-hard rock has better self-stability, the upper-soft rock has poor self-stability, and the phenomenon of soft-hard unevenness is obvious because the hard rock only appears in the lower half of the excavation face. The cross section of shield tunnel needs to pass through the composite strata with different composite height ratios of hard and soft rocks. The engineering geological conditions in this area are poor. The geological profile of the left tunnel is shown in Fig. 3. The proportion of hard rock in the tunnel face is between 0 and 1.

Concept of hard rock height ratio
The north extension section of Guangzhou Metro Line 8 is composed of upper-soft and lower-hard strata. The properties of the overlying soil layer of shield tunnel are basically the same; however, the hard rock along the tunnel is fluctuating, which leads to different hard rock proportions along different excavation faces. Hard rock height ratio is defined as the ratio of hard rock height to the excavation surface height of shield tunnel along the line: where is the height ratio of hard rock to stratum, ℎ is the height of soft rock, and is the total height of the tunnel excavation face.

Settlement monitoring point arrangement
On the basis of comprehensive considering the surrounding environment grade and the operability of monitoring implementation in Tongdewei-Shangbu Station, ground settlement monitoring in shield tunneling interval can be divided into longitudinal monitoring and lateral monitoring. Longitudinal monitoring points are installed every 5 m along the middle line of the tunnel. Lateral monitoring is distributed from the center of the shield vault to both sides, from near to far, and increases by 3-5 m between the measuring points. The detailed layout are shown in Figs. 7 and 8.  Burying monitoring points on hardened surface: firstly, drill a hole with diameter of 70 mm in hardened road surface with core pulling machine, and then drill a round head steel bar with length of 400-1000 mm and diameter of 12-20 mm into the hole. The head of the steel bar shall be 50 mm lower than the road surface, and the bottom shall be embedded in the undisturbed soil of the subgrade base, and then the sand shall be backfilled in the hole, and the steel bar shall be exposed to the sand surface for 20 mm. The embedding diagram is shown in Fig. 9. Burying monitoring points on the unhardened surface: Directly driving a round head steel bar with length of 400-1000 mm and diameter of 12-20 mm into the undisturbed soil. The steel bar shall be exposed to the surface for 50 mm, and the surrounding of the monitoring point is protected by bricklaying. The field monitoring is shown in Fig. 10.   Transverse surface subsidence is obviously affected by hard rock height ratio. The general trend is that the value of surface subsidence correspondingly decreases and settlement tank becomes shallow with the increase in hard rock height ratio. The surface subsidence in different hard rock height ratio strata is obviously different, and the maximum difference is about 8.6 mm; When the hard rock height ratio is 0.2 and 0.5, the curve of settlement tank does not change much. When the hard rock height ratio is 0.85 and 1, the curve also changes little. The maximum settlement difference of hard rock height ratio interval between 0.2-0.5 and 0.85-1 is only 1.4 mm. But in the hard rock height ratio interval between 0-0.2 and 0.5-1, the minimum difference of surface settlement is 3.2 mm. The settlement value reaches a maximum value of 13.2 when the hard rock height ratio is 0.   Fig. 13 shows the settlement curve, which is the surface longitudinal settlement under different hard rock height ratios during shield tunneling. From the figure, obvious uplift is observed at approximately 15-20 m ahead of the excavation face due to the influence of face pressure of palm, which belongs to negative stratum loss. At approximately 5 m ahead of the excavation face, surface deformation changes from uplift to subsidence with an accelerating trend. At the area of 5-10 m behind the excavation face, the surface subsidence rate is high, and the subsidence amount reaches approximately 83 % of the total subsidence. This subsidence which is influenced by the gap of shield structure belongs normal stratum loss. At approximately 20 m behind the excavation face, surface settlement is no longer affected by shield tunneling. Longitudinal surface subsidence is significantly affected by the hard rock height ratio, which shows that surface subsidence decreases with the increase in hard rock height ratio. The surface longitudinal settlement curves of tunnels with different hard rock height ratios are compared. The results show that with the increase of the hard rock height ratio, the location of surface deformation begins to decrease from approximately 15 m in front of the excavation face to approximately 10 m. Moreover, the surface subsidence rate has a decreasing trend with the increase in hard rock height ratio. Generally, with the increase of the hard rock height ratio, the influence degree and range of shield construction on the surface longitudinal settlement are decreasing.

Establishment of 3D numerical model
This study mainly investigates the law of surface subsidence caused by shield tunneling in the strata with different hard rock height ratios. In order to facilitate comparative analysis, five sections (YD170, YD203, YD200, YD161, and YD162) in Chapter 4 are selected for simulation. The model mesh is divided into tetrahedron + hexahedron mesh, which is a special hybrid mesh of MIDAS GTS. The grid can not only improve the calculation accuracy, but also improve the calculation speed. The grid in the excavation range of the tunnel needs to be divided into smaller ones, and the size is set to 0.5 m. The mesh is roughly divided at the boundary of the model, and the size is set to 3 m. Then, the dimension from the model boundary to the excavation decreases gradually according to the linear gradient. A finite element model with a hard rock height ratio of 0.5 is shown in Fig. 14.

Setting of main parameters of the model
The Shield machine used in tunnel excavation is simulated by plate element, the material is made of steel, the elastic modulus is 250 GPa, the Poisson's ratio is 0.2, and the bulk density is 78 kN/m 3 . The C50 concrete used in shield tunnel segment material is based on the Code for Design of Concrete Structure, which indicates that the elastic modulus of a segment should be 34.5 GPa. In view of the effects of annular and longitudinal splicing of segment on the reduction of segment stiffness, the elastic modulus of the segment in the model calculation is 80 % of the reduction value recommended by the code, that is, 27.6 GPa. The Poisson ratio of the segment is 0.2, and the plate element is used to simulate it. The segment grouting layer is equivalent to the mud water equivalent layer, the equivalent layer thickness is 0.3 m, and the solid element is used for simulation, and the hardening process of the slurry is not considered in the calculation process [39]. All of the above materials only consider the stress and deformation within the elastic range, that is, the elastic constitutive model is adopted. The physical and mechanical parameters of various materials of the model are shown in Table 1. In order to more accurately reflect the situation of shield excavation, the model sets the pressure of tunnel face, cylinder thrust force and grouting pressure. The pressure of tunnel face is 120 kPa according to experience. The cylinder thrust force and grouting pressure shall be determined according to the actual situation, in which the cylinder thrust force is taken as 100 kN/m and the grouting pressure is taken as 150 kPa.

Numerical calculation method
The numerical model in this paper adopts the Mohr Coulomb failure criterion. Mohr proposed the strength formula in 1900: = ( ). According to Mohr's theory, the yield of rock and soil is expressed as follows: where -cohesion of soil; -internal friction angle of soil The Eq. (4) shows that the shear strength of the material is related to the normal stress acting on the plane. It is not the maximum shear stress that causes material failure, but the most dangerous combination of -in a plane. In numerical analysis, Eq. (4) can be expressed by stress invariants as follows: sin sin + cos ( ) + cos = 0, where -equivalent compressive stress; -Mises stress; -lode angle.
In the principal stress space, it appears as an irregular hexahedral vertebral body surface. The section shape on the plane is an irregular hexagon. As shown in Fig. 15(a), the hexagon needs to be determined by the length of and . These two quantities can be expressed as follows: Since the geometric similarity of each plane of the Mohr Coulomb yield surface, the ratio = of any plane is constant. As shown in Fig. 15(b), the slope tan of the straight line does not change with the confining pressure (or hydrostatic pressure), so the criterion is accurate when the confining pressure is in a certain range. The confining pressure of this project will not be large enough to produce compressive failure, so the mohr-coulomb failure criterion can be used to obtain very accurate results. a) Yield surface shape of plane b) Yield surface shape of = − 6 ⁄ meridian plane Fig. 15. Shape of Mohr Coulomb yield surface

Shield excavation simulation
In order to improve the calculation efficiency and get close to the actual construction condition as much as possible, shield driving two ring segments (3 m) as a working condition, a total of 48 rings need to be driven after the completion of construction. Shield Excavation Simulation: (1) In the first step, the initial stress field is analyzed by applying dead weight and soil boundary conditions. The initial displacement is then reset to facilitate the study of the incremental displacement caused by tunnel excavation; (2) The second step is to "passivation" the soil mass of tunnel excavation, applies the pressure of tunnel face to stabilize the excavation face, "activates" the shield shell, simulates the shield tunneling process; (3) The third step "activates" the segment unit to simulate the assembly process of the shield tail segment; (4) The fourth step is to simulate shield tail synchronous grouting. "Activate" the slurry equivalent layer unit. The radial grouting pressure is applied on the grouting layer to consider the influence of grouting pressure on the stratum and segment. At the same time, "activate" cylinder thrust force to provide forward jacking force for the next excavation; (5) Repeat Steps (1)-(4) and continue digging until all segments are assembled.  It can be seen from the settlement nephogram that when the hard rock height ratio is 0.5, the surrounding rock in the lower half of the tunnel excavation face is basically in a stable state, and the impact of the disturbance is very small. There is a certain uplift phenomenon at the bottom of tunnel arch, which is mainly due to the unloading effect caused by tunnel excavation. There is a large settlement at the vault, and the maximum settlement is nearly 40 mm. Due to half of the soft stratum in the tunnel excavation face, it is necessary to strictly control the excavation parameters and adjust the excavation attitude, so as to effectively control the uneven settlement caused by the disturbance of the soft stratum. It can be seen from the calculation nephogram that there is obvious settlement trough phenomenon in the stratum settlement, and the surface settlement conforms to the Peck curve. As the buried depth of the tunnel is = 16.7 m, and the internal friction angle is = 10°-55°, the theoretical value 7.94 m-21.13 m can be obtained by substituting the Eq. (1). The computational nephogram shows = 8.5 m, which is in the theoretical range. Through the prediction of the value of the inflection point i and the range of the settlement groove, the range of the influence of shield construction on the stratum can be determined, and effective control measures can be taken in time to reduce the influence.

Analysis of calculation results
The transverse settlement curves corresponding to the hard strata height ratio with different proportions are plotted with = 10, 20, and 30 m.
For different hard rock height ratio strata, the lateral surface subsidence is immensely affected by the hard rock height ratio, and the surface influence range is between 3D and 3.5D when the values are 0, 0.5, and 1, with corresponding maximum settlement values of 12.8 9.7, and 4.6 mm, respectively (Fig. 18).
As shown in Fig. 18, after the completion of tunnel excavation, the curve of settlement trough is consistent with Peck's prediction, and the law of transverse surface settlement of each monitoring section is basically the same. The sequence of shield excavation has small influence on the law of transverse surface settlement. The depth of settlement tank increases and its width decreases with the increase in hard rock height ratio. The settlement tank is basically shallow and wide when the hard rock height ratio is = 1. The settlement tank gradually exhibits the characteristics of "deep and narrow" with the decrease of . Fig. 19 shows the simulation curve of the longitudinal surface settlement varying with hard rock height ratio. From the figure, the longitudinal settlement value of the tunnel decreases with the increase in hard rock height ratio. When the hard rock height ratio varies from 0 to 0.2, shield tunneling can affect the ground subsidence at approximately 15 m ahead of the excavation face and ground subsidence tends to be stable at approximately 20 m after excavation face. When the ground subsidence tends to be stable, the settlement value reaches 98 % of the total settlement. Shield tunneling only affects the ground subsidence at approximately 10 m before excavation face when the hard rock height ratio varies from 0.5 to 1 and ground subsidence tends to be stable at approximately 20 m after excavation face. Surface subsidence reaches 96 % of the total settlement at approximately 14 m after the excavation surface when the hard rock height ratio is 1. With the increase of the hard rock height ratio, the shield construction reduces the amount of the surface longitudinal settlement and its influence range. Generally, the simulated values fit well with the measured values. However, the simulation of longitudinal surface settlement does not show the phenomenon of surface uplift in front of the excavation face, which may be related to the methods of numerical simulation and the selection of pressure of tunnel face.

Comparative analysis of numerical results and monitoring results
In order to verify the reliability of the numerical model, the monitoring results are compared with the numerical calculation results. As shown in Fig. 20, the relation curve between the maximum surface subsidence above the vault and the hard rock height ratio is shown. The variation of the simulated value and the measured value with the hard rock height ratio is shown in the figure. It can be seen from the figure that the numerical simulation results and the measured results fit well, and the results of both show the same trend of change. The maximum surface subsidence decreases with the increase of the hard rock height ratio. When the hard rock height ratio is in the range of 0.2-0.5, the influence of the hard rock height ratio on the maximum surface settlement value is not significant. When the hard rock height ratio is in the range of 0-0.2 and 0.5-1, the change rate of the maximum subsidence value is greatly affected by the hard rock height ratio. The rate of change between 0-0.2 is the largest. Therefore, special attention should be paid to these two sections in order to prevent uneven ground settlement during the construction of sensitive sections.
From the monitoring results, the numerical calculation results can well predict the settlement trend under different hard rock height ratio, and the monitoring results verify the accuracy of the numerical calculation model.

Conclusions
In this paper, numerical simulation and site monitoring are used to study the influence of shield tunneling on the ground settlement in the upper soft and lower hard strata by introducing the concept of hard rock height ratio. The following conclusions are obtained: 1) The surface subsidence is obviously affected by the hard rock height ratio. With the increase of the height ratio of hard rock, the surface settlement value decreases, and the settlement trough changes from "deep narrow" to "shallow wide". The surface subsidence in different hard rock height ratio strata is obviously different, and the maximum difference is about 8.6 mm. When the height ratio of hard rock is 0, the surface subsidence reaches the maximum value, the maximum value is 13.2 mm.
2) The influence of the hard rock height ratio on the surface longitudinal settlement is mainly reflected in the position change of the beginning and the end of the settlement. With the increase of the hard rock height ratio, the location of surface deformation begins to decrease from approximately 15 m in front of the excavation face to approximately 10 m. The settlement of different hard rock ratio tends to be stable at approximately 20 m after excavation face.
3) The monitoring results verify the accuracy of the numerical results. The hard rock height ratio in the 0-0.2 and 0.5-1 is the sensitive interval, and the settlement value in these two hard rock height ratio interval varie greatly. The maximum settlement difference of hard rock height ratio interval in 0.2-0.5 is only 1.4 mm, but in the hard rock height ratio interval between 0-0.2 and 0.5-1, the minimum difference of surface settlement is 3.2 mm. It is necessary to pay attention to the uneven settlement of the ground surface caused by shield construction in the sensitive hard rock height ratio interval. Jianbing Lv holds a Doctorate from the School of Civil and Transportation Engineering, Central South University. He is now a Professor at the School of Civil and Transportation Engineering, Guangdong University of Technology. At present, it mainly studies the direction of rock, soil and road. In this paper is responsible for writing the full text and refining the central idea.

Xulong Li is a Master of Civil and Transportation Engineering College of Guangdong
University of Technology. The research direction is geotechnical engineering. In this paper is responsible for numerical simulation calculation and analysis.
Helin Fu, Professor of Tunnel and Underground Engineering Department of Central South University, doctoral supervisor, visiting Professor of Freiberg University of Technology and Wuyi University, consultant of Hualu company of Guangdong Communications Group, director of tunnel and underground engineering branch of China Civil Engineering Society, vice chairman of underground space Professional Committee of China Civil Engineering Society. In this paper is responsible for technical support and advice.