Stability Analysis of Surrounding Rock and Treatment Structures in Superlarge Karst Cave of Naqiu Tunnel

*is study relied on the treatment and construction process of the corridor hall-type large karst cave of Naqiu Tunnel, and field tests and numerical simulation were carried out to investigate the stability of surrounding rock and treatment structures during the tunnel excavation process. From the results of different excavation footages, both the displacement characteristics between the backfill pile-bearing platform-retaining wall compound treatment structures with surrounding rock and the stress changes of the initial support were analyzed. Combining with the field tests and simulation results, it can be obtained that the compound treatment structure under the bias voltage has unequal leftward displacement and settlement. When the distance to excavation surface gets closer, the bias effect will become greater on the horizontal displacement than on the vertical displacement. Both the displacement and stress changes tend to be stable after the excavation is completed, which proves that the treatment structure of the foundation, bearing platform, retaining wall, and lining structure on the steep slope section is stable and effectively suppresses the effect of the surrounding rock bias voltage. All above conclusions can provide theoretical basis and construction reference for similar projects.


Introduction
With the continuous development of infrastructure and economic construction in Central China and Western China, the transportation infrastructure has been developing rapidly.In the construction of highways and railways in hilly areas, it is necessary to build tunnels [1,2].It should be noted that karst disasters are frequent in the central and western regions, and the karst environment is various in different regions.e tunnel structure will be affected by both the natural geological form and artificial excavation, which will make the construction process more complicated and difficult [3,4].
In recent years, research studies on the influence of the karst cave on tunnel engineering have been carried out.Mengjun et al. [5] used model tests and numerical simulations to study the influence of the diameter ratio of the karst cave and tunnel and the relationship between the spatial position of the karst cave and tunnel geometry on the displacement, stress, and stability of the surrounding rock.Based on the karst tunnel project of the high-speed railway segment, Yangchun and Xiaojun [6] used FLAC3D to simulate the excavation stability of the tunnel with horizontal karst caves.Yuanhai et al. [7] studied the effect of karst caves on the stability of the surrounding rock for the Dayaoshan Tunnel and compared the filling and grouting reinforcement schemes when tunneling through the karst caves.Based on the elastic-plastic theory, Xianju and Mengtao [8] analyzed the stability of the surrounding rock and support of the karst tunnel using finite element software.Peiguo and Huajun [9] analyzed the large deformation of the surrounding rock of tunnels through filled karst caves and introduced the mechanism of the initial arching foot to cause the overall deformation of the tunnel and proposed the main treatment measures.Shubin et al. [10] conducted a model test for the tunnel excavation near the pressure-bearing chamber and obtained the stress-strain characteristics of the key locations in the tunnel excavation process at a certain cavity size and internal pressure.
e above research is concentrated on medium and small karst tunnels.However, the tunnel engineering research for the superlarge corridor hall cave is not yet thorough.It takes into account the development of a large cavity invading the tunnel profile, backfill foundation in the steep slope geological section, complex treatment structure, and other special factors.e stability analysis of complex structures such as pile foundation on steep slope-retaining wall-primary support and surrounding rock is still a major problem to be studied.erefore, the stability of wall rock and complex treatment structures during the excavation of extremely large karst tunnels is studied in depth in this paper through numerical simulation and field testing, taking the superlarge karst cave in the left line of Naqiu Tunnel as the engineering background.On the one hand, it verifies the rationality of the actual design and construction.On the other hand, it can provide guidance and theoretical basis for the treatment and excavation of similar karst tunnels.

Engineering Background
Naqiu Tunnel is located in Furong Town, Yongshun County, Xiangxi, Hunan Province.It is a one-way doubletrack tunnel with a total length of 2860 m. e superlarge karst cave of the corridor hall type is located on the left line K14 + 350-520.
e cave is 170 m long and is basically coincided with the direction of the tunnel.e location of the tunnel and cave is shown in Figure 1(a).
e surrounding rock of the tunnel is composed of microweathered limestone.It is a medium-thick layered structure and belongs to Class III surrounding rock.Its destruction form is dominated by falling blocks.e karst cave is basically developed and formed, and there are local surrounding rock structures with poor stability.ey belong to Grade IV surrounding rocks and may produce smallscale rockfalls.e original landform of the karst cave is shown in Figure 1(b).
According to the detailed geological survey report of the karst section, three schemes for cave treatment were put forward: bridge, line change, and backfill.From the aspects of reasonableness, safety, environmental protection, easy construction, and low cost, technical comparisons are considered.Specific backfilling geology and hydrological conditions are considered as well.Finally, the backfilling scheme is selected [11,12], as shown in Figure 1(a).
For the special steep slope section of K14 + 385−K14 + 415, it is also the focus of this study.e average height difference of the cave is 30 m, and the width is relatively large.e impact of steep slopes on the backfills in different backfill areas may produce large inhomogeneous deformation over time and reduce the stability of the backfilling foundation.In order to avoid the influence of vehicle load and structural load and bias load of surrounding rock after excavation, the backfill pile-retaining wall joint treatment scheme was uniquely adopted on the basis of the backfilling scheme.e implementation of this program can enhance the stability of the foundation.At the same time, it can form an integral structure with the surrounding rock and tunnel structure to better play the role of transmission load.e specific scheme is as follows.
As one of the supporting structures of the retaining wall and platform, the pile foundation adopts a 1.5 m × 1.2 m rectangular pile structure and is arranged at the steep slope of the bottom of the cave.e total number is 9, and the length of the pile is between 20 and 30 m. e pile body is applied as a retaining wall, and the layered backfilling of the hole slag is carried out together.After reaching the elevation of the pile top, the pile body is poured to form a lower backfill-piled foundation roadbed; the vertical and horizontal frame girder structure is laid below the C30 concrete platform, and it is connected with the rock-socketed pile foundation to form an integral structure.e C30 reinforced concrete retaining wall structure is laid on the left side of the tunnel lining structure to support the bias effect caused by the cave roof and the irregular surrounding rock excavation.At the later stage of construction, the small cavity on the right side of the tunnel lining is filled with C20 concrete pumping, as shown in Figure 2.

Calculation Model and Theoretical Parameters
As shown in Figure 2(b), the cross section of K14 + 385 is selected as the study object in this paper.Based on the backfill pile-retaining wall and initial support joint treatment scheme, the section is modeled and calculated by FLAC3D.According to the actual situation, the depth of the tunnel is about 70 m.e model uses the center of the tunnel as the origin of coordinates.e overall model is 120 m wide (x-axis, horizontal), 30 m long (z-axis, axial), and 146.27 m high (y-axis, vertically radial).e left and right model boundaries are 50 m away from the tunnel center construction site.e computational model is divided into 94080 elements and 99231 nodes.e mesh is all hexahedral elements, as shown in Figure 3.
e Mohr-Coulomb constitutive model is used for the material of surrounding rock and the backfill, and the elastic model is used for the pile, bearing platform, and retaining wall [13,14].According to the treatment scheme, geological survey, and rock mechanics test, the material mechanics parameters were determined, as shown in Table 1.
e boundary condition used in this model is the displacement constraint.
e bottom and four sides of the model are fixed by the normal displacement.e top surface uses the free boundary.e gravity acceleration is 10 N/Kg, and the lateral pressure coefficient is λ � 0.6.e excavation process was simulated according to the sequence of the initial ground stress field → displacement clearing → fullsection excavation → initial support → next excavation, which is carried out 10 times by using the "null" model in the FLAC3D, and the loop footage is 3 m.e results are used to analyze the stress and strain of the pile foundation, platform, 2 Advances in Civil Engineering retaining wall, surrounding rock, and initial lining structure under di erent processes.(1) e horizontal and vertical displacements of the three piles are all negative value; it proves that, after the excavation of the tunnel, there is a bias pressure from the surrounding rock, which makes the top of the pile arise leftward horizontal displacement and vertical settlement.

Displacement Analysis of the Pile Foundation
(2) e displacement on the left side of the pile top is slightly larger than that on the right side, the horizontal displacement value of no. 3 pile after the end of excavation reaches the maximum 1.79 mm, the vertical displacement of the pile foundation is much smaller than the horizontal displacement, and the maximum vertical displacement di erence between three pile foundations is 0.05 mm.erefore, the e ect of excavation on the horizontal displacement of the pile foundation is much greater than that on the vertical displacement.

Advances in Civil Engineering
(3) As the excavation progresses, the horizontal displacement of the pile foundation near the excavation face is larger.As shown in Figure 5(a), the horizontal displacement of no. 1 pile increased with construction progresses firstly.As the excavation proceeds, the leftward displacement of nos. 2 and 3 piles starts to increase, and when the excavation of the face is completed, the horizontal displacement of the pile top reaches its maximum.As shown in Figure 5(b), the vertical displacement of the three piles increases quickly in early construction and gradually begins to reduce after the excavation of 9 m, tending to the fixed value of 0.13 mm after the end of the excavation.e influence scope of overall excavation on the pile foundation is 9∼10 m (Figure 6).

Displacement Analysis of the Bearing Platform.
Figure 7 shows the layout of the monitoring point of the bearing platform, where nos.7 and 9-10 feature points are located below the retaining wall.According to the simulation results, the displacement of the monitoring point of the bearing platform is shown in Figure 8, and the negative value indicates horizontally left (or vertically downward), which is analyzed from the graph: (1) e displacements of nos.1-11 monitoring points are all negative value, which means that the bearing platform has a corresponding displacement under the superstructure load and surrounding rock bias pressure.In the later construction, it is necessary to backfill the left cavity of the bearing platform.

Advances in Civil Engineering
(2) e vertical displacement of the bearing platform is much smaller than the horizontal displacement, the horizontal displacement of top nos.4∼6 monitoring points is larger than that of bottom nos.1∼3, and the horizontal displacements of the monitoring points 4, 7, and 8 in the same section are successively decreased.e maximum horizontal displacement of the bearing platform is 2.7 mm, which appears on the top right corner in no. 4 point.(3) After the excavation, the in uence of the surrounding rock on the right side on the bearing platform is more signi cant than the structure of the retaining wall.When the distance to excavation surface gets closer, the horizontal displacement of the bearing platform will become greater.e vertical displacement of points 9 to 11 of the bearing platform increased rapidly during the early stage of excavation; the displacement peak is reached when the excavation reaches 9 m, and then it began to gradually decrease and tends to the value of 0.22 mm with construction progress.It can be seen that the impact of the construction on the vertical displacement of the bearing platform is 9 m.

Analysis of Retaining Wall Stress and Strain
4.3.1.Analysis of Retaining Wall Displacement.Figure 9 shows the layout of the retaining wall feature points.Nos.1-3 points are located at the intersection of the retaining wall and surrounding rock, and nos.4-6 points are located at the middle of the right side of the retaining wall.According to the results of simulation calculation, the displacement of the monitoring point of the retaining wall is shown in Figure 10.A negative value indicates that the displacement is to the left (or vertically downwards), as analyzed from the gure: (1) e displacements of the feature points are all negative value, which re ecting the excavation construction produces a bias pressure e ect.Meanwhile, it re ects that the retaining wall will play a supporting role under the bias of the surrounding rock.
(2) e vertical displacement of the retaining wall is less than the horizontal displacement, the horizontal displacement of the upper part is less than that of the middle part of the retaining wall, and the maximum value is 1.45 mm on no.6 point; the vertical displacement of the left side of the retaining wall is smaller than that of the right side, and the maximum value is 0.62 mm, which appears at the intersection of surrounding rock and retaining wall.Meanwhile, the vertical displacement of the bottom is also less than that of the other upper points.(3) As the excavation progresses, the overall horizontal displacement of the retaining wall also gradually increases, and the displacement rate of the lower part is greater than that of the upper part; at the same excavation section, when the distance to excavation surface gets smaller, the horizontal displacement of the retaining wall will become greater, and the peak 6 Advances in Civil Engineering value is reached after the completion of excavation.e vertical displacement rapidly increases during the initial period of excavation and reaches a maximum at 12 m, and then it decreases and tends to a constant value.It can be seen that the range of in uence of the construction on the vertical displacement of the retaining wall is about 12 m.(1) e maximum value of vertical compressive stress appears at the lower right corner of the retaining wall, with a maximum value of 3 MPa.At the same time, the vertical compressive stress changes from the initial 0.2 MPa to the maximum 0.75 MPa at the top of the retaining wall.(2) From the view of dynamic construction, the vertical maximum compressive stress in the lower right corner of the retaining wall is gradually moved forward along the excavation direction, and the stress value also increases as the excavation distance increases and reaches the peak value when the excavation is completed.(3) For the maximum compressive stress of the right lower corner of the retaining wall, the construction quality should be strengthened, and the stability of the connection between the bottom of the retaining wall and the bearing platform should be maintained.Meanwhile, for the trend of displacement to the left karst cave at the retaining wall, it is necessary to back ll the left karst cave with a backpressure to keep the overall structure stable after the completion of the construction.

Displacement Analysis of Tunnel Surrounding Rock.
Figure 12 shows the displacement change of the surrounding rock of the tunnel at the excavation section Y 0, in which the positive value indicates the horizontal direction to the right (vertically upward) and the negative value indicates the horizontal direction to the left (vertically down), as analyzed from the gure: (1) After excavation, the horizontal displacement of the right side wall and right arch of the tunnel appears convergent to the left with the maximum of 2.56 mm. e rest of the monitoring points are horizontally o set to the right, and the maximum rightward displacement appears in the vault of 2.8 mm.
(2) e vertical displacement of surrounding rock is far greater than the horizontal displacement, and there is a large upward displacement of 4.5 mm at the right arch foot of the tunnel.e vertical displacement of the rest of the feature points is all settlement, and the maximum of vertical settlement of 4.6 mm appears in the vault.e settlement of the right arch waist is larger than that of the left arch waist, which proves that the supporting structure of the left retaining wall plays a key role in suppressing the vertical displacement.(3) As the excavation progresses, the vertical and horizontal displacements of each monitoring point increase with the increase of excavation distance.

Advances in Civil Engineering
After excavation of 15 m, the change rate decreased and gradually tended to stable.It can be seen that the range of in uence of tunnel excavation on the front section is 15 m.

Stress Analysis of the Initial Support.
Figure 13 shows the vertical stress of the initial support structure under di erent tunneling footage.e negative value represents the compressive stress, and the positive value represents the tensile stress, as analyzed from the gure: (1) During the excavation process, the stress of the initial lining support gradually increases with the increase of the excavation distance.e overall structure is mainly presented as a compressive stress, and the maximum compressive stress of 5 MPa appears at the support point on the right wall of the tunnel and the right arch foot.e compressive stress value of the whole support is 1∼5 MPa, which is far less than the ultimate compressive strength of 15.4 MPa of C20 shotcrete, which indicates that the initial support structure is safe and stable.(2) From the results of the nal stress, it can be seen that there is a tensile stress of 1 to 2 MPa at the right arch position of the initial support structure, and at the same time, there is a stress di erence between each section of support, so the quality of the connection between the supporting structures in each cycle tunneling footage should be enhanced.Meanwhile, the quality and thickness of the initial support should be increased in the stress concentration area of the right side wall.

Stability Test and Analysis of the Pile Foundation.
e pile top strain test adopted the JMZX-215A intelligent embedded concrete strain gauge, and the strain gauge was placed on the top of the main reinforcement of the pile to test the concrete strain at the top of the pile before the pile was poured.Figure 14 shows the test conditions of no. 2 pile foundation below the eld-retaining wall foundation: It can be seen from Figure 14 that a negative strain value was produced at the pile head, indicating that the superstructure loads and the excavation have produced compressive stress.In 1 to 10 days, small vertical strains of the pile head were generated due to the concrete pouring of the bearing platform.e next 11∼25 days, the vertical strains of the pile head increased slowly as the concrete layered pouring of the retaining wall.And the vertical strain rate of the pile increased rapidly after 29 days, indicating that the upper retaining wall structure was a ected by the pressure subjected to the surrounding rocks, and thus, the pressure of the bearing platform and the pile head was also produced, which made the strain of the pile increase rapidly.At the end of the excavation, the vertical strain of the pile began to decrease gradually and became stable with the completion of the construction, indicating that the pile foundation structure was stable.

Stability Test and Analysis of the Retaining Wall.
e stress test at the junction of the retaining wall and the surrounding rock of the karst cave was based on the JMZX-5010am intelligent pressure cell and the JMDL-3210A singlepoint displacement meter.Figure 15(  In Figure 15(a), when the excavation was carried out, the surrounding rock pressure acting on no. 1 pressure cell was increased rapidly, and the rate of change reached 0.15 MPa/d and increased to 0.27 MPa within 2 to 3 days.In contrast, the pressure change of nos. 2 and 3 pressure cells was lagging behind, and an obvious increment in the 8th and 20th days started to generate.e pressure values of all the above three pressure cells tended to be stable nally, indicating that the con ning pressure was basically stable after excavation.In Figure 15(b), the horizontal displacements of the retaining wall increased totally as the excavation distance increased.e horizontal displacements at the midpoint of the retaining wall structure during the initial excavation were small but increased at a rate of 0.25 mm/d as the excavation surface was approaching,nally slowed down, and tended to a stable value in the 18th day.In summary, the maximum pressure of the retaining wall of 0.33 MPa and the measured leftward horizontal displacement value of midpoint of 1.25 mm are all smaller than the simulation results.It can be seen that the stress released by the stress redistribution of the surrounding rock due to the surface excavation was small.

Structure Test and Analysis of the Initial Support.
JMZX-5040am dual-mode pressure cells were, respectively, installed on the tunnel vault, tunnel arch (right), and right side wall of the K14 + 400 section to test the pressure change between the surrounding rock and the initial support.
As shown in Figure 16, after the initial support was constructed, the pressure value of the lining increased with the following excavation, and the change rate at the characteristic point of the right side wall increased fastest and reached a maximum value of 0.95 MPa. e pressure of the tunnel vault acting on the lining structure had the minimum value of 0.39 MPa, which indicated that the right con ning pressure of surrounding rock had great in uence on the right side wall.e bias pressures of the three measuring points were far less than the ultimate compressive strength of the initial lining, so we concluded the lining structure was stable.Furthermore, the pressure change rate was almost close to zero in the 28th day, indicating that the deformation of the surrounding rock was stopped.From the above analysis, the stability period of surrounding rock in this tunnel section was inferred in 28∼30 days.

Displacement Test and Analysis of Surrounding Rock.
e internal displacement of the surrounding rock test of the K14 + 385 section was conducted using a JMDL-3210A single-point displacement meter, and its speci cations were 1.5 m, 2.5 m, and 3.5 m, respectively.e settlement test and convergence test of the tunnel vault were performed with level measurement and a JSS30A aggregator, respectively.As shown in Figure 17, a positive value represents the displacement of the surrounding rock was biased toward the inside of the tunnel, and a negative value indicated the tunnel vault had a settlement and tended to convergence.
(1) Both the internal displacement and convergence of the surrounding rock had a tendency of coming toward the inside of the tunnel.e displacement of the surrounding rock at internal 1.5 m was larger than that at other two measuring points, which is 1.4 mm, while the internal displacement of the surrounding rock at 3.5 m is the smallest, of 0.79 mm, which indicates the closer the measuring point near the loose circle, the bigger the displacement of the measuring point would be, and vice versa.the change of displacement of the whole surrounding rock was basically zero, which indicated that surrounding rock reached a stable state and the change period of the surrounding rock displacement in this section is 22∼25 days.

Comparison of Field Test and Calculation Simulation.
Comparison between eld tests and numerical simulation results is given in Table 2, and it shows the following: (1) Compared with the results of eld tests and numerical simulations, the stress and strain change laws of each structure during the excavation stage were basically the same.e displacement error was less than 2 mm, and the stress error was less than 1 MPa.ere is a great agreement between the eld test results and simulation results.
(2) e strain test results of the pile head were slightly smaller than the simulation results.Meanwhile, there was no large o set and settlement in the overall structure of the site, which proved that the design of the pile foundation was reasonable and applicable.
(3) e pressure at the top of the retaining wall on the site as well as the displacement at the midpoints was smaller than the simulation results, indicating that the actual surrounding rock pressure is lower than the design value.At the same time, according to the onsite observation, the roof of the karst cave had not collapsed and become unstable after the treatment structure construction, which veri ed that the supporting e ect of the back ll pile foundation-bearing platform-retaining wall-combined supporting structure was good.(4) Due to the relative lagging nature of the on-site support process, the settlement and convergence of the tunnel vault on-site were larger than the simulation results by 1 to 2 mm.Combined with the results of the eld test and numerical simulation, there were no large displacement changes; thus, the stability of the surrounding rock of the superlarge karst cave was re ected as good from another point of view.(5) ere may be no su cient contact between the pressure cell and the surrounding rock, which led to  the result that both the pressure of the initial support and the pressure of the surrounding rock were less than the simulation results.e maximum surrounding rock pressure stress appeared on the right wall on the site, which was far less than the ultimate compressive strength of the design of the support.Also, the deformation of the arch foot did not appear in the field, which reflected the supporting effect was good.

Conclusion
Both the characteristics and change laws of stress and strain of the surrounding rock and treatment structures in the superlarge karst cave of Naqiu Tunnel were discussed under different excavation distances situation.
e following conclusions can be drawn from the comparison between field tests and numerical simulation: (1) During the excavation, the top of the pile foundation produced a larger leftward displacement under the bias pressure of surrounding rock, which increases with the excavation distance.e vertical displacement is less than the horizontal displacement and gradually decreases after the excavation of 9 m. e influence scope of overall excavation on the pile foundation is 9-10 m. (2) e bearing platform and the retaining wall also produced a leftward displacement under the bias pressure.e effect of excavation on the horizontal displacement is much larger than that on the vertical displacement.A large stress concentration appears in the lower right corner of the retaining wall, and the stability of the connection between the bottom structure of the retaining wall and the bearing platform should be maintained during construction.At the same time, to ensure the stability of the operation in the future, backfill is needed for the left karst cave in the later stage of construction.(3) Compared with the horizontal displacement, the vertical displacement of the surrounding rock around the tunnel was greater, and the maximum vertical settlement of 4.6 mm appeared at the tunnel vault.e settlement of the tunnel arch (right) was larger than that of the left arch, and it proved that the left retaining wall support structure contributed a restraining effect to the vertical displacement.With the construction progresses, the horizontal and vertical displacements became greater with the increase of the excavation distance, and the rate of change decreased and gradually became stable after excavation of 15 m; it can be referred the influence scope of the surrounding rock excavation was 15 m.(4) e maximum compressive stress value of the supporting structure is far less than its ultimate compressive strength, and the supporting effect is proved to be good based on on-site testing and observation.
(5) e results of numerical simulation and field tests proved that the combined treatment structure of backfill pile foundation-bearing platform-retaining wall initial lining structure in the steep slope was stable and reasonable, and it can effectively restrain the bias pressure effect of surrounding rock, which can provide theoretical basis and construction reference for similar projects in the future.
Figure 4 shows the position of the plane pile foundation.According to the simulation calculation, the displacement of the pile top under di erent excavation distances is shown in Figure 5.A negative value indicates horizontally to the left (or vertically downward), which is analyzed from the graph:

Figure 2 :Figure 1 :
Figure 2: Treatment structure: (a) the longitudinal section of the treatment structure; (b) the cross section of K14 + 385.

Figure 4 :
Figure 4: Position of the pile foundation.
Figure 11  shows the stress changes of the retaining wall under four working conditions.Negative values indicate compressive stress, and positive values indicate tensile stress.e following can be seen from the gure:

Figure 10 :
Figure 10: Displacement of the retaining wall at di erent excavation distances: (a) horizontal displacement; (b) vertical displacement.

Figure 11 :Figure 12 :
Figure 11: Vertical stress change of the retaining wall under di erent working conditions: (a) excavation of 3 m; (b) excavation of 12 m; (c) excavation of 21 m; (d) excavation of 30 m.
a) shows the test results of 3 pressure cells in the front, middle, and rear of the top of the retaining wall, while Figure 15(b) shows the test result of the horizontal displacement meter in the middle point of the retaining wall.Information available from Figure 15 is given below.

Figure 14 :Figure 15 :
Figure 14: Strain test of the pile.(a) Vertical strain of no. 2 pile top.(b) Site installation.

Figure 16 :
Figure 16: Displacement test of surrounding rock.(a) Change in surrounding rock and initial lining pressure of the K14 + 400 section.(b) Site layout.

Figure 17 :
Figure 17: Displacement test of surrounding rock.(a) Internal displacement of surrounding rock of the K14 + 385 section.(b) Dome settlement and peripheral convergence of the K14 + 385 section.

Table 2 :
Comparison between eld test results and numerical simulation results.