Next Article in Journal
Examining the Effects of Induced Plant Defenses on Spodoptera frugiperda Performance
Previous Article in Journal
Special Issue “Renewable and Sustainable Energy Conversion Systems”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Characteristics of Surrounding Rock and Design of Backfill Material Parameters for Tunnels Passing through Giant Caverns and Underground Rivers

1
School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
2
Beijing Ruiwei Century Geotechnical Engineering Limited Company, Beijing 100080, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(8), 3906; https://doi.org/10.3390/app12083906
Submission received: 5 March 2022 / Revised: 4 April 2022 / Accepted: 6 April 2022 / Published: 13 April 2022

Abstract

:
The influence on tunnel construction of karst and underground rivers is always an important problem in tunnel engineering. In order to demonstrate the rationality of backfill parameters and the effectiveness of supports under the influence of groundwater when a tunnel passes through a large karst cave, the finite element software FLAC3D was used for numerical analysis. Converting the mean values and standard deviations of mechanical response results of the surrounding rock and the supports on multiple sections into the ECULID distance from the origin point on a two-dimensional plane as evaluation indexes, the influence of the reinforcement parameters on the mechanical response of the surrounding rock and the supports was analyzed by orthogonal experiments. Based on fuzzy decision theory, by regarding the ECULID distance between the simulated result of each group and the global optimal value of the multiple evaluation index as the comprehensive evaluation index, a backfill parameter design method was proposed. By comparing the results which used optimal and worst parameters in the FLAC3D, 10 times and 2.5 times differences in longitudinal and horizontal displacement were observed, respectively. Then, the optimal backfill parameters were applied to the actual project for verification. The field monitoring results showed that the optimal backfill parameters can effectively reduce the displacement around the tunnel. After constructing a diversion for the underground river, the water flow in the karst cave did not rise during a rainstorm. This study provides a reference for the design and construction of other projects in the future.

1. Introduction

Limestone and karst are widely distributed in Central and Western China. The distribution area of soluble rocks in China is more than 3.4 × 104 km2, accounting for about one third of the land area, mainly in Yunnan and Guizhou in the southwest [1]. Karst landforms are common geological obstacles in the construction of mountain tunnels and urban underground tunnels. In the complex geological environment, a tunnel needs to pass through strata containing karst caves, and various geological disasters, such as collapses, rock bursts and water and mud inrushes, make for a large number of accidents in tunnel construction [2,3], as shown in Table 1.
A huge risk of karst disaster is water gushing. There have been many studies on water-filled limestone caves that are in front of tunnels [4,5,6,7], including field tests, model tests and numerical simulations.
Li et al. [8], based on the actual conditions of the engineering section from Wangfuzhuang Station to Dayangzhuang Station of the Jinan Metro Line R1, derived a calculation formula of the reserved safe thickness of the tunnel face given the presence of a water-filled limestone cave using COMSOL. Li et al. [9] established a 3D water-filled karst cave and tunnel model, simulated the tunnel’s excavation under different height–span ratios and revealed the variation law of stress and displacement during tunnel excavation. Pan et al. [10], through a model experiment to study the failure process of tunnel faces at the front of water-filled limestone caves, divided the failure processes in the model into the production of group cracks, the formation of water inrush channels and the complete collapse of water-resistant plates. Meng et al. [11] independently developed a model test system for water and mud inrush in a tunnel to study the structural instability characteristics and seepage law of karst pipeline fillers with different permeability coefficients. Gao et al. [12] took Kaiyuansi Tunnel as an example, carried out physical model and numerical simulations and analyzed the water leakage, vertical displacement and water pressure results for the monitored section. Li et al. [13,14,15], based on the highly dangerous karst tunnel in Qiyue Mountain and numerically simulated five tunnel floor water inrush cases using FLUENT software, analyzed the variation characteristics of velocity and pressure and discussed and summarized the water flow characteristics after water inrush at the tunnel floor. However, research on the impacts of underground rivers on tunnels are not enough [16,17].
The existing research mostly focuses on small-size karst caves. Zhao et al. [2] carried out similar simulation experiments and a simulation analysis to study the influence of karst cave size on the stability of surrounding rock. Cheng et al. [18] used the principle of limit analysis to analyze tunnel damage due to collapses caused by the existence of karst caves and obtained the damage range of intermediate surrounding rock. Huang et al. [19], in order to predict the collapse range of a rock mass when excavating a tunnel above a karst cave, used a variational method within the framework of an upper bound theorem to obtain an analytical expression for the collapse surface. Most current research is focused on the location relationship, distance and relative size of karst caves and tunnels.
For super-large tunnels, there is still little research. Chen et al. [20] analyzed the Naqiu Tunnel as an example. The karst cave is 170 m long, which is 10 times the radius of the tunnel. Due to the karst cave not being high, the method of drilling, grouting and bridging has been selected to pass through. Cao et al. [21] took Yangqiaoba Tunnel as an example; the maximum length, width and height of the karst cave are 267, 180 and 120 m respectively. A changing line was selected for the left tunnel and backfilling for the right tunnel. Ma et al. [22] took Longlingong Tunnel as an example; the bottom axial length of the karst cave is 171 m, the width is 65 m, and the distance to the ground surface is 110 m. The roof of the karst cave collapsed to the surface. Finally, the crossing mode of backfilling the cavity was determined and the treatment mode of dynamic compaction plus grouting was adopted to improve the foundation.
Generally, large karst caves are mainly dry. Common crossing methods include route diversion, subgrade backfilling and excavation and bridge construction. The specific method to be used should be determined according to the actual situation and cost. However, the number of projects crossing large karst caves with water are few, and there are few studies that have been made on this subject worldwide. Therefore, it is necessary to study the safety of passing through large karst caves with water.
This paper takes the Yujingshan Tunnel as an example to study the impact of rising groundwater on the project. Firstly, an actual investigation of the surrounding rock properties and groundwater conditions of the Yujingshan Tunnel was carried out; the main parameters of the backfill body include: elastic modulus, density, cohesion, internal friction angle, porosity and softening coefficient of water to surrounding rock. Orthogonal experiments were used to design a test group, and the six parameters of backfill layer settlement, tunnel vault settlement, convergence on both sides and the internal force and bending moment of the initial support were used as evaluation indicators. The finite element software FLAC3D was used for a simulation, and the values of the evaluation indicators in each experimental group were extracted to calculate the means and standard deviations. Then, the range standardization process was performed, and the means and standard deviations for the experimental group were normalized to the [0,1]. d E C U L I D i k is the ECULID distance between ( x i k , σ i k ) and the origin of coordinates. The backfill parameter corresponding to the minimum value of d E C U L I D i k is the optimal parameter. Comparisons with on-site monitoring were made to verify the correctness of the method described in this article. This study provides a reference for engineering design and construction in the future.

2. Tunnel Overview

2.1. Karst Cave Overview

As shown in Figure 1, the tunnel in the Chengdu–Guiyang Railway is located in the section between Xingwen county and Weixin county, within Yunnan Province. The design speed is 250 km per hour. There are double tracks in a single tunnel. The longitudinal slope of the line is a single upward slope of 30‰, with a total length of 6306.28 m and a maximum buried depth of 350 m.
The tunnel mainly passes through coal measure strata and soluble rock strata, of which 71 m of section passes through the coal-bearing stratum of the Permian Changxing Formation (P2C), which is a high gas section; 439 m of section crosses the coal-bearing stratum of the Longtan Formation (P2l), which is a gas outburst section, and 795 m of section crosses the limestone of the Qixia Maokou Formation (P1m + Q) of the Lower Permian system, and karst is strongly developed there.
When the guide gallery was excavated to d3k279 + 948, a huge karst cave was revealed, in the shape of a dome, about 100 m along the line direction and about 200 m perpendicular to the line. The height of the karst cave is 50~120 m. It is high on the right side of the line and low on the left side. The tunnel is located near the top of the karst cave. Most of the extent of the karst cave is stable. However, under the influence of cracks, joints, groundwater and blasting, local block falling from the surface of the cave wall has occurred. The thickness of the bottom accumulation layer is 30~90 m, and its composition is complex, including soft clay, gravelly soil and medium coarse sand. An underground river has developed at the foot of the slope on the left side of the line in the karst cave. The river is 5~15 m wide and flows from the right side of the line to the left. The flow of the underground river during the dry season was observed to be 4~15 m3/s, and the maximum flow during the rainy season is about 20.16 m3/s. According to the drilling and geophysical data, it is very likely to be unstable and damaged due to the influence of branch pipeline water at the bottom and in surrounding areas, the scouring and erosion of the underground river and unbalanced loading on the upper part of the later backfill body. See Figure 2 for details [23].

2.2. Underground River Overview

An underground river has developed at the foot of the slope on the left side of the line within the karst cave, its surface elevation about 961 m, which is about 114 m, below the tunnel rail, and a measurable section length of the underground river is about 739 m. The underground river roughly intersects with the line at 70°, developing in an S shape, and the runoff is discharged from the right side to the left side of the line. The overall runoff direction is N64°E, which is basically developed along the bedding plane of surrounding rock, with a flow of 13.55~20.16 m3/s. The upper reaches of the underground river enter the cave hall 30 m from the left side of d3k279 + 985. The entrance of the underground river is nearly rectangular, about 14 m wide and about 1.6 m high.
Downstream, the underground river flows out of the cave hall at 281 m on the left side of d3k279 + 895, which is a deep pool. From 15 October to 23 October 2016, it rained continuously in the tunnel site. On 26 October 2016 (the third day after the rainfall stopped), the water surface of the pool fell to 30 cm below the roof, with a width of about 8 m.
In order to explore the maximum water inflow of the underground river, the atmospheric precipitation infiltration method was used to predict the water inflow in the catchment area in the section the tunnel passes through. Application of the atmospheric precipitation infiltration method involved selecting areas with similar hydrogeological and meteorological characteristics in the area the tunnel passes through and using multi-year average precipitation data to estimate water inflow in the tunnel. The calculation formula is as follows:
Q s = 2.74 × λ × h × F
where Q s is the calculated water inflow (m3/d); λ is the infiltration coefficient (L/s km2); h is the average rainfall (mm); and F is the distribution area of the aquifer (km2).
According to the investigation of the catchment area of the underground river, the total catchment area is about 85 km2, so F = 85 km2. The infiltration coefficient is taken to be for soluble rock: λ = 0.71 L/s km2. In order to calculate the most unfavorable situation, the historical daily maximum rainfall in Weixin County was taken, so that H = 132 mm (1 July 2008). With these parameters, Equation (1) was used to calculate the maximum water inflow of the underground river: Q s = 8 × 106 m3/d = 92.6 m3/s.
After the underground river was revealed, hydrological observation points were set at the outlet and middle part of the underground river. According to the calculation formula for ditch flow in the tunnel design manual, the calculation formula for these aspects is as follows:
Q = ω c ( R i ) 1 / 2
where Q is the flow of underground (m3/s); V is average velocity (m/s), V = c ( R i ) 1 / 2 ; ω is the cross-sectional area of overflow (m2); R is hydraulic radius (m); R = ω / X ; C is the coefficient, C = ( 1 / n ) R y ; X is the wet perimeter, which refers to the length of the contact surface between the fluid and the trench wall in the section (m); y is the coefficient, y = 2.5 n 1 / 2 0.7 R 1 / 2 ( n 1 / 2 0.1 ) 0.13 ; and n is the roughness coefficient of the trench wall, which can be obtained from the relevant table.
According to the field measurement, the minimum cross-sectional area of the underground river entering the cave hall is about 29 m2, rectangular in shape, 13 m wide and 2.23 m high. According to on-site observation, there are obvious and strong scouring traces on the wall of the intake, indicating that the intake must be filled with an underground river in rainy seasons. The maximum water inflow of the main pipeline intake of the underground river in the Karst hall, Q = 338 m3/s, can be obtained.
The underground water of the cave wall is not developed; there is dripping at the top of the cave. There is seasonal fissure water in some cracks, and many karst pipelines are developed at the bottom of the cave hall.
As shown in Figure 3, on the whole, the karst cave is mainly affected by the underground river. The underground river has been diverted many times in history, and there are many empty pipelines in the lower part of the karst cave. These pipelines also form a mainstream–tributary relationship with the current main underground rivers. These pipelines will have a great adverse impact on the stability of the upper cavity structure. In the rainy season, the groundwater level will rise rapidly, which will have a significant erosion effect on the karst cave. At the top of the giant karst cave, there is a partially permeable part, which is mainly affected by the upper overburden.

2.3. Solution

The section from d3k279 + 855 to d3k279 + 960 of the tunnel passes through the huge karst cave. Combined with the hydrogeological conditions and the development of the karst cave, many expert joint reviews were organized, and it was decided to adopt the overall treatment scheme of “underground river diversion + karst cave backfilling + roof reinforcement + bridge crossing” for passing the karst cave.

2.3.1. Karst Cave Backfilling

In order to provide the excavation construction platform for the tunnel arch, reduce the risk of falling blocks from the tunnel arch and eliminate the influence of the instability of the tunnel wall in the later stage, the karst cave was backfilled. Due to the maximum drop of the karst cave of 120 m, in order to ensure that the backfill body is dense, it is necessary to backfill in layers. As shown in Figure 4, the karst cave was backfilled in five levels. The karst cave was backfilled step by step through channel 4 (988 m elevation), channel 6 (1008 m elevation), channel 7 (1028 m elevation), channel 1 (1048 m elevation) and the upper steps of the tunnel (about 1082 m elevation).
Considering that the karst cave below the elevation of 1008 m is mainly located in the underground river section (close to the bottom of the karst cave), there are many through-cracks in the hall, making for water seepage. In order to ensure drainage fluently at the bottom of the karst cave and prevent the impact of excessive bottom water pressure on the superstructure, the part below the elevation of 1008 m was backfilled with block stones, the block diameter not less than 30 cm. In order to ensure fluent drainage of groundwater such as fissure water and karst pipeline water near the wall of the karst cave, block stone backfilling could also be adopted within 3 m near the wall of the karst cave, the block diameter also not less than 30 cm.

2.3.2. Bridge Crossing

In order to avoid backfilling body settlement during tunnel operation and ensure driving safety, a bridge structure is to be adopted for crossing the karst cave section. In order to prevent the shotcrete of the tunnel wall spalling and falling off, light steel scaffolding shall be set on the bridge for protection to ensure operational safety. The bridge structure crossing the karst cave shall include a 38.85 + 108 + 38.85 m short side span prestressed concrete continuous beam with a beam height of 5–8 m.

2.3.3. Roof Reinforcement

Considering the location relationship between the tunnel and the karst cave, the different stability of surrounding rock and the requirements of bridge construction, overhaul and maintenance, the tunnel crossing the karst cave section has an A-type straight side wall composite lining and a B- and C-type curved side wall composite lining. The B- and C-type sections are located in the scope of the karst cave hall, and anchor mesh shotcrete was used as a permanent protective measure. In order to ensure the operational safety and the stability of the karst cave roof, it was provided with anchor cables to strengthen the protection, as shown in Figure 5.

3. Numerical Simulation

In order to study the influence of groundwater in the karst area on tunnel construction when crossing the karst cave and selecting appropriate backfill materials, the 279 + 923 section with the lowest design elevation of tunnel was selected as the research object, and a model was established using the finite difference software FLAC3D for a simulation. Figure 6a is the grid element division diagram, and Figure 6b is the distribution diagram of surrounding rock in the karst cave. So that the tunnel width on both sides is five times that of the width of the karst cave, the calculation range of the model was set to be 50 m in the longitudinal direction, 1300 m in the transverse direction and 700 m in the vertical direction, including 270 m from the tunnel crown to the model top surface, 120 m from the tunnel bottom plate to the karst cave bottom and 400 m from the model bottom. The simulation range was relatively large, and a total of 1 million elements were generated. Only the upper surface is free; all other surfaces are constrained. Therefore, the surrounding and bottom boundaries were set to be impervious.
Groundwater rise was simulated by using the seepage module that exists in the FLAC3d software. When the groundwater flow rises, the seepage calculation mode is turned on, and a stable pore pressure field is obtained through single seepage calculation, and then the fluid modulus Kf is set to 0 to achieve the mechanical equilibrium calculation [24,25].
During the supplementary survey of karst caves and underground rivers, 11 holes were drilled in the karst cave hall, with a total length of 575.6 m. A total of nine groups of rock samples and two groups of underground river water samples were obtained on site for the experiment. In addition, TEM, geological radar, seismic imaging, the active source surface wave method and the natural source surface wave method were also used to carry out geophysical exploration on the basement of the karst hall and the excavated tunnel. The properties of surrounding rocks inside the karst cave hall were determined by integrating various test methods, as shown in Table 2.
In order to study the impact of the rise of groundwater on the tunnel’s construction, this paper uses orthogonal experiments to study the different effects of the mechanical properties, the degree of compaction and the softening coefficient of water on surrounding rock of backfill materials on the tunnel project. In this paper, six main factors are selected for analysis: the elastic modulus, density, cohesion, internal friction angle, porosity and reduction coefficient. The standard orthogonal table L25 (56) is used to design the experimental team, and each parameter takes five levels, as shown in Table 3. The variation amplitude of the elastic modulus and cohesion is 20%, the variation amplitude of the internal friction angle and porosity is 10% and the variation amplitude of the density and reduction coefficient is 5%.
Counting the numerical calculation results of nine sections which are at intervals of 5 m from Y = 5 m to 45 m, a total of six average parameters were taken as an evaluation index, including the displacement of 0~20 m under the bottom plate, the crown settlement, the convergence of both sides and the internal force and bending moment of initial lining. In order to make the orthogonal experimental results comprehensively reflect the influence of groundwater on tunnel construction, according to fuzzy mathematics, after taking the absolute value of experimental results in each group, the average value x i ¯ and standard deviation σ in the group can be calculated according to Formulas (3) and (4). In practical engineering, it is always required that the average level of each evaluation index is small to ensure the safety of tunnel excavation. The degree of dispersion is low, so there is no local deterioration. Therefore, the smaller the average value and standard deviation in the group, the more favorable the situation is for tunnel excavation. The whole test group is range standardized according to Formula (5), and the average value and standard deviation of the test group are normalized to the interval of [0–1]. According to Formula (6), d E C U L I D i k , which is the distance between ( x i k , σ i k ) and origin, can be calculated as follows:
x ¯ i k = 1 n j = 1 m | x i j k |
σ i k = 1 n j = 1 m ( | x i j k | x ¯ i k ) 2
{ x i k = x ¯ i k x ¯ min x ¯ max x ¯ min σ i k = σ i k σ min σ max σ min
d E C U L I D i k = ( x i k * ) 2 + ( σ i k * ) 2 2
where i is the experimental group, 1 i 25 ; j is the serial number of the tunnel section, 1 j 9 ; k is the sequence number of six surrounding rock stability evaluation indexes: tunnel crown settlement, backfill body displacement under the bottom plate, convergence around the tunnel and internal force and bending moment of initial lining, 1 k 6 ; x ¯ i j k represents the value of k evaluation index of j tunnel section in i test group; x ¯ max , σ max , x ¯ min and σ min is the average value of the k evaluation index and the maximum and minimum value of standard deviations in 25 experimental groups, respectively; and x i k and σ i k are the range standard value of the average value and the standard deviation of k evaluation index in group test i, respectively.
The smaller the d E C U L I D of the experimental group, the higher the degree of optimization. Using this method can not only avoid redundant calculations when predecessors analyze multiple tunnel sections or multiple factors but also avoid the one-sidedness of single section and single factor analysis.
Table 4 illustrates the results of d E C U L I D evaluation indexes corresponding to the different parameters of each test group. The test group of factors E 2 ρ 4 c 5 φ 1 n 2 k 3 represents the elastic modulus, density, cohesion, internal friction angle, porosity and water reduction coefficient for the backfill body, and their values are 800 MPa, 2520 kg/m3, 700 kPa, 35°, 0.18 and 15%, respectively. Based on the negative index characteristics of each evaluation index, the influence of six backfill parameters on different evaluation indexes within the given factor level range is obtained, and the optimal group for different evaluation indexes is obtained, as shown in Table 5. When the crown settlement is taken as the index, the elastic modulus of backfill soil has the greatest influence, followed by the reduction coefficient of water on backfill, cohesion, porosity and density. The internal friction angle has the least influence. When the elastic modulus, density, cohesion, internal friction angle, porosity and water reduction coefficient are 1400 MPa, 2160 kg/m3, 300 kPa, 45°, 0.16 and 5% respectively, the optimization of the tunnel crown settlement is the best.
In order to comprehensively consider the influence of various mechanical parameters of backfill soil on tunnel construction under the influence of groundwater and optimize the scheme design for the backfill soil, the weighted relative deviation distance minimum method [26] was adopted to form a comprehensive evaluation scheme for backfill reinforcement with the whole tunnel as the evaluation scheme, which is based on the principle of fuzzy mathematics and multi-objective comprehensive decision making. The schemes are compared and selected and combined with orthogonal experimental analysis.
According to the principle of fuzzy mathematics, the six evaluation indexes are assigned according to importance degrees. Therefore, the importance degrees of backfill bottom displacement, crown settlement, left-side convergence, right-side convergence, initial lining axial force and initial lining bending moment were determined as 20%, 20%, 15%, 15%, 15% and 15% respectively. Therefore, the fuzzy set of evaluation index weight can be obtained as X = (0.20,0.20,0.15,0.15,0.15,0.15). We employed Equation (7):
{ g i = g i ( ν i , ν 0 ) = 1 x [ k = 1 m ( x k δ k ( i ) ) 2 ] 0.5 x = 1 m k = 1 m x k
where m is the number of evaluation indexes; x is the average value of index weight; x k is the weight value corresponding to k evaluation index; g i * represents the distance between the i test group and the optimal effect within the given factor level range. The smaller the g i * value, the higher the optimization level of the test group.
When the proposed weight fuzzy set is X, the range situation of g i is shown as in Figure 7. The main influencing factor obtained by the comprehensive evaluation method is the elastic modulus of the backfill body and the selected optimal scheme is E 5 ρ 1 c 3 φ 1 n 1 k 1 . If the decision makers adopt different evaluation index weight preferences, the results will be different, but the optimization design of backfill reinforcement parameters can still be carried out according to the method proposed in this paper.
The strength parameter group of the reinforcement optimized using this method is defined as group 26, and the weighted deviation distances of 26 conditions are calculated at the same time. The results are shown in Table 6.
It can be determined that when considering the comprehensive evaluation index, the dECULID distance between this method and the optimal evaluation index group is the smallest, which is 0.0057, indicating that the backfill material scheme designed by this method can effectively reduce the settlement around the tunnel and control the mechanical response of the surrounding rock and lining at a low level, and the mechanical response of each section is evenly distributed. The design goal of backfill reinforcement is achieved.

4. Effect of Comparion

In order to show the impact of different backfill materials on tunnel construction, the largest group E 1 ρ 5 c 5 φ 5 n 5 k 5 and the smallest group E 5 ρ 1 c 4 φ 5 n 4 k 1 of g i in Table 6 are selected for comparison. The two groups are compared to ensure the feasibility of the methods used in this study.

4.1. Backfilling Body Displacement

Figure 8 illustrates the settlement of two different backfill materials under the action of self-weight. The settlement trends of the two different backfill materials are almost the same but the difference in settlement amplitude is large. The main displacement of backfill is concentrated in the upper part of the karst cave. However, due to the principle of coordination deformation between continuous grids in the numerical simulation software, the top part of the backfill which is in contact with limestone is not the largest, and the displacement of this part is from 0 to the maximum, which is slightly different from the actual situation.
When using the optimal backfill material parameters, the maximum settlement of the backfill body is 70 cm, and the settlement of most of the area is below 35 cm. When using the worst backfill material parameters, the maximum settlement of the backfill body reaches 1.18 m and the average settlement is about 55 cm, which is 1.7 times and 1.6 times that of the best material, respectively. At the same time, it can be clearly found that when the worst backfill material is used, the settlement gradient is large, which is not conducive to the compaction of the backfill body and will cause certain damage to the tunnel construction. Therefore, it is necessary to roll the backfill in layers to ensure the compactness of each layer of soil, which can effectively reduce the settlement of the backfill under the influence of self-weight and construction.

4.2. Tunnel Longitudinal Displacement during Tunnel Construction

As shown in Figure 9, the longitudinal displacement of two backfill materials during tunnel construction is quite different after completing backfilling and removing the settlement under the action of their own weight. Due to the tunnel being located near the roof of the karst cave, the mechanical properties of the surrounding rock on both sides of the tunnel vary greatly from the limestone to the backfill body. Therefore, the longitudinal displacement distribution around the tunnel is different. The maximum crown settlement of the backfill body with the optimal parameters is 6 mm, and the maximum uplift of the tunnel bottom is only 2 mm. For the backfill with the worst parameters, the maximum crown settlement and the maximum uplift of the floor are 15 cm and 2 cm, respectively, which are 2.5 times and 10 times more than those under the optimal parameters. Therefore, it can be clearly discovered that the optimal parameters obtained by using the reinforcement parameter design method proposed in this study can effectively reduce the longitudinal displacement of the tunnel, which all meet the national specifications.

4.3. Tunnel Horizontal Displacement during Tunnel Construction

Figure 10 illustrates the tunnel lateral convergence of both sides during tunnel excavation. When using the optimal parameters, the maximum displacements on the left and right sides of the backfill are only 1.8 cm and 2 cm. For the backfill with the worst parameters, the maximum displacements on the left and right sides reach 4.5 cm and 5 cm, which are 2.5 times the values obtained when using the optimal parameters. Compared with the longitudinal displacement, the difference between the two parameters in the transverse displacement is smaller.
Firstly, this is because the distance of backfilling on both tunnel sides is small, and the limestone on both sides ensures the effective horizontal deformation of the backfilling body. Secondly, the time for finishing backfill is too short and horizontal stress is small. At the same time, the displacement of the left side of the tunnel is greater than that of the right side because the right side of the tunnel is closer to the limestone and the deformation is smaller.
In general, the application of the optimal parameters obtained above can effectively reduce the displacement around the tunnel, reduce the surrounding rock stress and provide a safety guarantee for tunnel construction.

5. Field Monitoring

In order to ensure the safety of the tunnel construction process and explore the impact of the underground river on the tunnel, the flow of the underground river, the layered settlement of backfill and the settlement of the tunnel vault were monitored, respectively.

5.1. The Flow of the Underground River Monitoring

In order to achieve a clear understanding of the underground river flow, 130 m to the left of d3k279 + 938 and 70 m to the left of d3k279 + 970 were taken as observation sections in the underground river in the karst cave to observe and measure the flow in the tunnel. The statistics and the amount of rain on that day are plotted in Figure 11.
The two black lines in Figure 11 represent the flow of the underground river in the karst cave and the outlet flow, respectively. It can be seen that the outlet flow is significantly greater than the flow of the underground river in the karst cave, indicating that there are multiple tributaries in Yujing Mountain, and because the limestone belongs to soluble rock, it forms many connecting drainage routes, resulting in the collection of the whole water at the outlet. The change in rainfall and the change in underground river flow in the karst cave are relatively slight, but they have little correlation with the flow at the outlet of the underground river. Therefore, the water system is complex, and so we should pay attention to the monitoring of groundwater to prevent the water level of the underground river rising quickly.
There is a deep pool at the outlet of the karst cave. According to the flow monitoring data for the outlet, the maximum water volume of the outlet is Q = 68.75 m3/s, which is far less than the calculated maximum water inflow and underground river water inflow. Therefore, in the rainy season, the water level in the karst cave must rise. Based on the static water sediment on the surface of the accumulation at the bottom of the cave and the drilling, it is judged that the water level of the underground river in the cave once rose to the elevation of 990 m (an increase of about 30 m). Due to the great uncertainty of the drainage property of the underground river channel, it is very likely that the underground river pipeline will collapse and block the river channel during the tunnel operation, which will further aggravate the rise of the underground river in the karst cave.
In order to avoid the backfill blocking the underground river channel and adversely affecting the operation, the underground river in the karst cave section was diverted before the karst cave treatment, as shown in Figure 12.
The underground river channel mainly uses drainage and blocking to prevent the rise of underground river water level in the rainy season from having a great impact on the safety of the tunnel in the giant karst cave. After the construction of the discharge tunnel, the underground river shall be blocked upstream and downstream. At the same time, the underground river flow in the karst cave section shall be led downstream of the underground river through the NO. 4 and 5 discharge tunnels.
In addition, a new channel is to be constructed at the intersection between the underground river and the tunnel, which will be able to maintain the rapid flow of the underground river. At the foot of the slope of the karst cave, the new channel follows the flow direction of the underground river and is about 30–50 m away from the tunnel wall on the right side, and a circuitous discharge tunnel is set to divert the underground river section so as to direct the upstream water of the underground river downstream through the discharge tunnel. The discharge tunnel is about 450 m long, 7.2 m wide and 6.0 m high.

5.2. Layered Backfill Monitoring

In order to ensure the safety of the tunnel after backfilling, the settlement of layered backfill soil is monitored at d3k279 + 915, as shown in Figure 13. This paper focuses on monitoring the settlement of backfill near the underground river, so the monitoring range is 20~40 m below the tunnel floor, and a monitoring point is set every 2 m.
At the initial stage after backfilling, the backfill body of each layer is gradually consolidated under the action of self-weight and upper construction load, displacement fluctuating frequently within the range of 5~−10 mm (positive value indicates soil uplift).
After 15 January 2019, the backfilling and compaction stage were completed, and the project entered the tunnel construction process. The first step was the tunnel step excavation. There was obvious settlement to backfill of each layer of d3k279 + 915, of which the settlement at 20 m is the largest, reaching 53 mm, and the settlement of each layer surrounding rock was concentrated in 40~50 mm. During this period, the settlement of each layer basically decreased with the increase of depth, which corresponded to the actual situation.
After the completion of tunnel excavation, grouting reinforcement was carried out for the backfill areas at the bottom and both sides, so the settlement of surrounding rock of each layer tended to decrease. Among them, the displacement of deep surrounding rock decreased significantly, and the maximum reduction degree was 10 mm. Therefore, it can be considered that grouting to the bottom of the tunnel can effectively improve the surrounding rock’s mechanical condition.

5.3. Crown Settlement Monitoring

In order to explore the state of the surrounding rock at the top of the tunnel in different construction stages, three sections d3k279 + 875~896 were selected for analysis, as shown in Figure 14a. The final crown settlement of 12 sections in d3k279 + 896~955 in the open excavation period and concealed excavation period was counted, as shown in Figure 14b.
By analyzing the data in the diagram, we can determine that:
(1)
The crown settlement increases gradually with the progress of construction, and the growth trend of each curve is basically consistent. During the open excavation period, the crown settlement gradually increases with the construction process of the step method, and it can be clearly seen in Figure 14a that it has an obvious growth gradient. In the period of concealed excavation, due to the influence of grouting and other reinforcement measures, the crown settlement is fluctuating and continues to rise as part of the whole trend.
(2)
The crown settlement of the d3k279 + 885 section is the largest of the three sections in both open excavation and concealed excavation, and the maximum settlement in the two periods reaches 28.5 mm and 37.3 mm, respectively. In addition, the crown settlement can reach a stable state through reinforcement measures, such as continuous grouting behind the lining.
(3)
Crown settlement is considered to be the main displacement in the tunnel. The final settlement of each section at the end of the open excavation period and concealed excavation period is shown in Figure 14b. The settlement range of the whole monitoring section is 0.3~13.6 mm and 4.4~15.2 mm, respectively. The crown settlement difference between sections is relatively large. Therefore, it is necessary to pay attention to the change in load caused by the change in surrounding rock stability and to strengthen the monitoring of the surrounding rock in different areas and carry out timely supplement grouting for reinforcement.

6. Discussion

Taking the tunnel of Chengdu Guiyang railway crossing the huge karst cave section as project example, this study took the underground river in the karst cave as the research object, used the finite element software FLAC3D and the measured data to explore the adverse impact of underground river level increases on tunnel construction and selected the best parameters for backfill material for practical engineering. The main findings of this study are as follows:
(1)
In this study, the commonly used parameters, including density, elastic modulus, cohesion, internal friction angle, porosity and parameter reduction coefficient in backfill, were selected as the main variables. The six data of layered settlement, tunnel crown settlement, convergence on both sides and internal force and bending moment of initial lining obtained by the finite element software FLAC3D were used as evaluation indexes. The average value and standard deviation of the mechanical response results for the surrounding rock and lining structure in multiple sections were transformed into the ECULID distance from the origin on the two-dimensional plane as the evaluation index. The parameter influence analysis was carried out based on the orthogonal test. It was found that the bulk elastic modulus of backfill and the reduction coefficient of water on materials have the greatest influence on structural mechanics.
(2)
According to the field measured data of 10 different sections, the layered settlement of the backfill body presented three different stages of slow increase, rapid increase and slow increase in the backfilling period, open excavation period and concealed excavation period. The growth range of each section was different but the trend was basically the same.
(3)
According to the measured data of tunnel crown settlement in different sections, in the open excavation period, the crown settlement increased in step with the advance of the construction steps, and the tunnel crown settlement fluctuated continuously in the concealed excavation period. However, due to the different nature of the surrounding rock, the support method and grouting amount, the settlement of each section at the end of the open excavation period and the concealed excavation period varied greatly, but they were all within a controllable range, and the impact of groundwater on tunnel construction was basically eliminated.
In general, the tunnel crossing the giant cavity was well controlled.

7. Conclusions

(1)
Based on the multi-objective comprehensive decision-making theory, the global optimal value of multiple evaluation indexes was used as the optimal backfilling effect, and the multi-dimensional ECULID distance between the test group and the global optimal solution was used as the comprehensive evaluation index. A determination method for reinforcement parameters was proposed. When the weight of the evaluation index was determined, the design group of backfilling reinforcement parameters with the optimal overall mechanical response could be obtained. In this study, the optimal scheme was E 5 ρ 1 c 3 φ 1 n 1 k 1 .
(2)
The optimal parameter scheme obtained by the backfill parameter design method was applied to practical engineering. The monitoring results showed that the displacement of the tunnel mainly occureed in the backfill stage, and after optimizing the backfill parameters, the displacement of the tunnel could be effectively reduced and the stability of the surrounding rock could be improved. At the same time, after the diversion of the underground river, the water flow in the karst cave hall was significantly reduced, and there was no situation of poor drainage.
(3)
The research presented in this paper can provide a research basis and engineering reference for tunnel engineering passing through karst caves. In the next stage, we will introduce the tomographic analysis method in the process of selecting parameters, which provides a firmer basis for the selection of engineering evaluation indicators. In addition, a small tunnel model was established to simulate the changing process of groundwater erosion and damage to the tunnel backfill.

Author Contributions

Conceptualization, J.W. and Y.W.; data curation, J.W. and F.P.; formal analysis, J.W. and Y.W.; investigation, F.P.; validation, J.W; writing—original draft preparation, J.W. and Y.W.; writing—review and editing, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Li, L.; Tu, W.; Shi, S.; Chen, J.; Zhang, Y. Mechanism of water inrush in tunnel construction in karst area. Geomat. Nat. Hazards Risk 2016, 7 (Suppl. S1), 35–46. [Google Scholar] [CrossRef] [Green Version]
  2. Zhao, M.; Ao, J.-H.; Liu, X.-H.; Wang, B. Model test study on the effect of karst size on the stability of surrounding rock in a tunnel. Chin. J. Rock Mech. Eng. 2004, 23, 213–217. [Google Scholar]
  3. Li, K.-L.; Wang, M.-H.; Wang, Q.; Yang, B.; Zhang, Q. Research on treatment measures for large semi-filled karst caves in Yangjiao No. 1 Tunnel. High-Speed Railw. Technol. 2019, 10, 91–95. [Google Scholar]
  4. Cui, Q.-L.; Wu, H.-N.; Shen, S.-L.; Xu, Y.-S.; Ye, G.-L. Chinese karst geology and measures to prevent geohazards during shield tunnelling in karst region with caves. Nat. Hazards 2015, 77, 129–152. [Google Scholar] [CrossRef]
  5. Wang, Y.; Yin, X.; Jing, H.; Liu, R.; Su, H. A novel cloud model for risk analysis of water inrush in karst tunnels. Environ. Earth Sci. 2016, 75, 1450. [Google Scholar] [CrossRef]
  6. Cheng, W.-C.; Cui, Q.-L.; Shen, J.S.-L.; Arulrajah, A.; Yuan, D.-J. Fractal Prediction of Grouting Volume for Treating Karst Caverns along a Shield Tunneling Alignment. Appl. Sci. 2017, 7, 652. [Google Scholar] [CrossRef] [Green Version]
  7. Yang, X. Research on Water Inrush Failure Mode and Catastrophic Evolution Mechanism of Tunnel Hydraulic-Filled Karst Caves; Shandong University: Jinan, China, 2017; p. 92. [Google Scholar]
  8. Li, L.P.; Xiong, Y.F.; Wang, J.; Gao, X.C.; Wang, K.; Sun, H.C.; Fang, Z.D. Comprehensive Influence Analysis of Multiple Parameters on the Safety Thickness against Water Inrush in Shield Tunnel. Int. J. Geomech. 2020, 20, 04020226. [Google Scholar] [CrossRef]
  9. Zhang, L.W.; Fu, H.; Wu, J.; Zhang, X.Y.; Zhao, D.K. Effects of Karst Cave Shape on the Stability and Minimum Safety Thickness of Tunnel Surrounding Rock. Int. J. Geomech. 2021, 21, 04021150. [Google Scholar] [CrossRef]
  10. Pan, D.; Li, S.; Xu, Z.; Lin, P.; Huang, X. Experimental and numerical study of the water inrush mechanisms of underground tunnels due to the proximity of a water-filled karst cavern. Bull. Eng. Geol. Environ. 2019, 78, 6207–6219. [Google Scholar] [CrossRef]
  11. Meng, Y.; Jing, H.; Yin, Q.; Wu, X. Experimental study on seepage characteristics and water inrush of filled karst structure in tunnel. Arab. J. Geosci. 2020, 13, 450. [Google Scholar] [CrossRef]
  12. Gao, C.-L.; Zhou, Z.-Q.; Yang, W.-M.; Lin, C.-J.; Li, L.-P.; Wang, J. Model test and numerical simulation research of water leakage in operating tunnels passing through intersecting faults. Tunn. Undergr. Space Technol. 2019, 94, 103134. [Google Scholar] [CrossRef]
  13. Li, S.; Wu, J.; Xu, Z.; Li, L.; Huang, X.; Xue, Y.; Wang, Z. Numerical analysis of water flow characteristics after inrushing from the tunnel floor in process of karst tunnel excavation. Géoméch. Eng. 2016, 10, 471–526. [Google Scholar] [CrossRef]
  14. Li, S.C.; Wu, J.; Xu, Z.H.; Li, L.P.; Huang, X.; Xue, Y.G.; Wang, Z.C. A Case Study for Escape Route Optimization after Water Inrush in a Backward Excavated Karst Tunnel. In Proceedings of the Fourth Geo-China International Conference, Jinan, China, 25–27 July 2016; pp. 100–108. [Google Scholar]
  15. Wu, J.; Li, S.C.; Xu, Z.H. Numerical analysis of gas-liquid two-phase flow after water inrush from the working face during tunnel excavation in a karst region. Bull. Eng. Geol. Environ. 2019, 78, 2973–3010. [Google Scholar] [CrossRef]
  16. Yang, W.; Yang, X.; Fang, Z.; Shi, S.; Wang, H.; Bu, L.; Li, L.; Zhou, Z. Model test for water inrush caused by karst caves filled with confined water in tunnels. Arab. J. Geosci. 2019, 12, 749. [Google Scholar] [CrossRef]
  17. Kaufmann, G.; Romanov, D. Modelling long-term and short-term evolution of karst in vicinity of tunnels. J. Hydrol. 2020, 581, 124282. [Google Scholar] [CrossRef]
  18. Lyu, C.; Yu, L.; Wang, M.; Xia, P.; Sun, Y. Upper bound analysis of collapse failure of deep tunnel under karst cave considering seismic force. Soil Dyn. Earthq. Eng. 2020, 132, 106003. [Google Scholar] [CrossRef]
  19. Huang, F.; Zhao, L.; Ling, T.; Yang, X. Rock mass collapse mechanism of concealed karst cave beneath deep tunnel. Int. J. Rock Mech. Min. Sci. 2017, 91, 133–138. [Google Scholar] [CrossRef]
  20. Chen, H.; Sha, C. Stability Analysis of Surrounding Rock and Treatment Structures in Superlarge Karst Cave of Naqiu Tunnel. Adv. Civ. Eng. 2018, 2018, 4842308. [Google Scholar] [CrossRef]
  21. Cao, X.; Zhang, L.; Lin, Y. Study on the treatment plan of large dry karst caves in Yangqiaoba Tunnel. Mod. Tunn. Technol. 2014, 51, 185–190. [Google Scholar]
  22. Ma, T. Research on Treatment Technology of Large Karst Cave No.1 Longlingong Tunnel of Yiwan Railway. Railw. Stand. Des. 2010, 33, 125–128. [Google Scholar]
  23. Ye, K. Research on Comprehensive Treatment of Karst Hall in the Newly-Built Yujingshan Tunnel of Guizhou Line; Southwest Jiaotong University: Chengdu, China, 2018; p. 79. [Google Scholar]
  24. Forcellini, D. Analytical fragility curves of shallow-founded structures subjected to Soil-Structure Interaction (SSI) effects. Soil Dyn. Earthq. Eng. 2021, 141, 106487. [Google Scholar] [CrossRef]
  25. Triantafyllaki, A.; Papanastasiou, P.; Loukidis, D. Numerical analysis of the structural response of unburied offshore pipelines crossing active normal and reverse faults. Soil Dyn. Earthq. Eng. 2020, 137, 106296. [Google Scholar] [CrossRef]
  26. Peng, P.; Zhang, D.; Sun, Z. Design method of surrounding rock deformation characteristics and reinforcement parameters for tunnels with weak interlayers. Chin. J. Rock Mech. Eng. 2021, 40, 2260–2272. [Google Scholar]
Figure 1. Map of the Chengdu–Guiyang high-speed railway.
Figure 1. Map of the Chengdu–Guiyang high-speed railway.
Applsci 12 03906 g001
Figure 2. The situation of the karst cave.
Figure 2. The situation of the karst cave.
Applsci 12 03906 g002
Figure 3. Underground river at the foot of the slope in the karst cave.
Figure 3. Underground river at the foot of the slope in the karst cave.
Applsci 12 03906 g003
Figure 4. Backfilling schematic diagram of the cross-sectional treatment of the giant karst cave.
Figure 4. Backfilling schematic diagram of the cross-sectional treatment of the giant karst cave.
Applsci 12 03906 g004
Figure 5. Site picture of tunnel crown reinforcement.
Figure 5. Site picture of tunnel crown reinforcement.
Applsci 12 03906 g005
Figure 6. The numerical model and its construction. (a) Unit division. (b) Stratigraphic distribution.
Figure 6. The numerical model and its construction. (a) Unit division. (b) Stratigraphic distribution.
Applsci 12 03906 g006
Figure 7. Range analysis of comprehensive evaluation method.
Figure 7. Range analysis of comprehensive evaluation method.
Applsci 12 03906 g007
Figure 8. Backfilling body displacement. (a) The worst backfill material parameters. (b) The optimal backfill material parameters.
Figure 8. Backfilling body displacement. (a) The worst backfill material parameters. (b) The optimal backfill material parameters.
Applsci 12 03906 g008
Figure 9. Range analysis of comprehensive evaluation method. (a) The worst backfill material parameters. (b) The optimal backfill material parameters.
Figure 9. Range analysis of comprehensive evaluation method. (a) The worst backfill material parameters. (b) The optimal backfill material parameters.
Applsci 12 03906 g009
Figure 10. Range analysis of the comprehensive evaluation method. (a) The worst backfill material parameters. (b) The optimal backfill material parameters.
Figure 10. Range analysis of the comprehensive evaluation method. (a) The worst backfill material parameters. (b) The optimal backfill material parameters.
Applsci 12 03906 g010
Figure 11. Rainfall and flow monitoring map in the cave.
Figure 11. Rainfall and flow monitoring map in the cave.
Applsci 12 03906 g011
Figure 12. Schematic diagram of underground river diversion.
Figure 12. Schematic diagram of underground river diversion.
Applsci 12 03906 g012
Figure 13. Time-curve of cumulative layered settlement of section d3k279 + 915: (a) 20~30 m; (b) 32~40 m.
Figure 13. Time-curve of cumulative layered settlement of section d3k279 + 915: (a) 20~30 m; (b) 32~40 m.
Applsci 12 03906 g013
Figure 14. Crown settlement of sections. (a) Crown settlement of the d3k279 + 875~895 section. (b) The final vault settlement of each section during open excavation and concealed excavation.
Figure 14. Crown settlement of sections. (a) Crown settlement of the d3k279 + 875~895 section. (b) The final vault settlement of each section during open excavation and concealed excavation.
Applsci 12 03906 g014
Table 1. Social and economic losses caused by karst adverse geological disasters in recent years.
Table 1. Social and economic losses caused by karst adverse geological disasters in recent years.
Unfavorable Geological Phenomena of KarstSocial and Economic Losses
Foundation collapse of the Lizhan railwayEconomic loss reached 1 million yuan
Large water and mud inrush in Yesanguan Tunnel of Yiwan Railway10 deaths and the economic loss was more than 13 million yuan
Water and mud inrush and surface collapse in the Hejiazhai TunnelConstruction period delayed by half a year and the economic loss was 10 million yuan
Water inrush in the Maluqing Tunnel11 deaths
Water and mud inrush in the Dayaoshan TunnelEconomic loss exceeeds 10 million yuan
Table 2. The parameters of the surrounding rock.
Table 2. The parameters of the surrounding rock.
Type of Surrounding RockElasticity Modulus
E/MPa
Poisson’s
Ratio
µ
Density
p/kg m−3
Cohesion
c/kPa
Internal Friction Angle
φ
Porosity
n/%
Permeability
Coefficient
k/(m2/Pa s)
Limestone10,0000.2226502000650.051 × e−11
Backfilling Crushed stone5000.25220030400.301 × e−1
Backfilling Crushed slag10000.252400500450.201 × e−6
Original crushed stone1000.3022005380.305 × e−2
Original crushed clay100.35186010070.201.2 × e−6
Table 3. Levels of sensitivity factors.
Table 3. Levels of sensitivity factors.
Factor NameElastic Modulus
E/MPa
Density p/kg m−3Cohesion
c/kPa
Internal Friction
Angle φ
Porosity
n/%
Reduction Coefficient
k
Level 16002160300350.165%
Level 28002280400400.1810%
Level 310002400500450.2015%
Level 412002520600500.2220%
Level 514002640700550.2425%
Table 4. Standard orthogonal experimental results (indicated by dECULID).
Table 4. Standard orthogonal experimental results (indicated by dECULID).
Serial NumberFactor Level CombinationBackfill Floor DisplacementCrown
Settlement
Left
Convergence
Right
Convergence
Initial
Lining Shaft
Initial Lining Bending
Moment
1 E 1 ρ 1 c 1 φ 1 n 1 k 1 0.7000.3670.5260.2220.1340.321
2 E 1 ρ 2 c 2 φ 2 n 2 k 2 0.6740.4910.6520.5590.4510.413
3 E 1 ρ 3 c 3 φ 3 n 3 k 3 0.6420.7240.7910.8550.7360.535
4 E 1 ρ 4 c 4 φ 4 n 4 k 4 0.6410.9020.6460.8460.9340.730
5 E 1 ρ 5 c 5 φ 5 n 5 k 5 0.70110.809111
6 E 2 ρ 1 c 2 φ 3 n 4 k 5 0.2750.5430.7820.6500.3190.744
7 E 2 ρ 2 c 3 φ 4 n 5 k 1 0.2240.4540.3650.4180.6490.342
8 E 2 ρ 3 c 4 φ 5 n 1 k 2 0.2840.4480.4040.4640.8480.424
9 E 2 ρ 4 c 5 φ 1 n 2 k 3 0.4140.4910.6280.6090.7690.478
10 E 2 ρ 5 c 1 φ 2 n 3 k 4 10.5210.6980.2650.0810.596
11 E 3 ρ 1 c 3 φ 5 n 2 k 4 0.0510.4260.3810.5100.6790.291
12 E 3 ρ 2 c 4 φ 1 n 3 k 5 0.0940.9260.5440.6090.3670.388
13 E 3 ρ 3 c 5 φ 2 n 4 k 1 0.1210.2540.3260.3660.6610.207
14 E 3 ρ 4 c 1 φ 3 n 5 k 2 0.5760.1020.4150.0360.0750.298
15 E 3 ρ 5 c 2 φ 4 n 1 k 3 0.4330.2520.6450.4650.4370.550
16 E 4 ρ 1 c 4 φ 2 n 5 k 3 0.0490.2450.3160.2890.4550.107
17 E 4 ρ 2 c 5 φ 3 n 1 k 4 0.0330.3200.3600.3840.7120.237
18 E 4 ρ 3 c 1 φ 4 n 2 k 5 0.4920.4360.7440.4720.1230.807
19 E 4 ρ 4 c 2 φ 5 n 3 k 1 0.1820.0890.2690.2650.6090.078
20 E 4 ρ 5 c 3 φ 1 n 4 k 3 0.3140.1940.2370.2210.1990.222
21 E 5 ρ 1 c 5 φ 4 n 3 k 2 00.1250.2150.1760.4920.116
22 E 5 ρ 2 c 1 φ 5 n 4 k 3 0.1740.1110.2480.1230.2110.151
23 E 5 ρ 3 c 2 φ 1 n 5 k 4 0.2060.2250.3370.11200.179
24 E 5 ρ 4 c 3 φ 2 n 1 k 5 0.2190.2130.4240.2560.3120.281
25 E 5 ρ 5 c 4 φ 3 n 2 k 1 0.14400.2130.2280.6270
Table 5. Direct analysis results of L25 (56) standard orthogonal experiment.
Table 5. Direct analysis results of L25 (56) standard orthogonal experiment.
Evaluating IndicatorPrimary and Secondary FactorsOptimal Combination
Backfill floor displacement E > c > ρ > φ > k > n E 5 ρ 1 c 4 φ 5 n 4 k 1
Crown settlement E > k > c > n > ρ > φ E 5 ρ 1 c 1 φ 3 n 1 k 1
Left convergence E > k > c > φ > ρ > n E 5 ρ 2 c 4 φ 5 n 5 k 1
Right convergence E > k > c > φ > n > ρ E 5 ρ 1 c 1 φ 2 n 1 k 1
Initial lining shaft c > φ > E > ρ > k > n E 5 ρ 1 c 1 φ 1 n 5 k 5
Initial lining bending moment E > k > φ > ρ > c > n E 5 ρ 2 c 4 φ 1 n 3 k 1
Table 6. Calculation results of weighted relative deviation distance.
Table 6. Calculation results of weighted relative deviation distance.
Item Number g i Item Number g i Item Number g i Item Number g i Item Number g i Item Number g i
10.032520.040130.051440.055550.064560.0411
70.029180.034090.0392100.0465110.0289120.0405
130.0236140.0244150.0335160.0183170.0259180.0397
190.0192200.0168210.0152220.0119230.0149240.0203
250.0178260.0057
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Wang, J.; Peng, F. Study on the Characteristics of Surrounding Rock and Design of Backfill Material Parameters for Tunnels Passing through Giant Caverns and Underground Rivers. Appl. Sci. 2022, 12, 3906. https://doi.org/10.3390/app12083906

AMA Style

Wang Y, Wang J, Peng F. Study on the Characteristics of Surrounding Rock and Design of Backfill Material Parameters for Tunnels Passing through Giant Caverns and Underground Rivers. Applied Sciences. 2022; 12(8):3906. https://doi.org/10.3390/app12083906

Chicago/Turabian Style

Wang, Yu, Jiachen Wang, and Feng Peng. 2022. "Study on the Characteristics of Surrounding Rock and Design of Backfill Material Parameters for Tunnels Passing through Giant Caverns and Underground Rivers" Applied Sciences 12, no. 8: 3906. https://doi.org/10.3390/app12083906

APA Style

Wang, Y., Wang, J., & Peng, F. (2022). Study on the Characteristics of Surrounding Rock and Design of Backfill Material Parameters for Tunnels Passing through Giant Caverns and Underground Rivers. Applied Sciences, 12(8), 3906. https://doi.org/10.3390/app12083906

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop