Numerical and Field Measurement Investigation on Foundation Pit Excavation Blasting of Anchor in Suspension Bridge

The foundation pit of a suspension bridge project in the Three Gorges Reservoir area is investigated in this paper. The pit is located under an unstable rock mass and landslide body; its base lithology is mudstone. The bridge foundation pit project mainly adopts blasting excavation to accelerate construction progress. However, as a hazardous technique to engineering safety, the explosion vibration easily causes deterioration of the surrounding strata, thereby inducing slope instability and rock mass collapse. Besides, three major challenges should be considered: complex terrain conditions, difficulties in the blasting excavation of anchors, and the extremely high risk of construction. Therefore, comprehensive risk control measures using the methods of hierarchical excavation and minimum charge blasting are put forward. After the measures were verified to be feasible through finite element simulation, it was successfully applied to actual construction. In addition, this paper proposes using fiber concrete to reinforce slope retaining walls, and simulates the reinforced effect based on the research above. The results indicate that the risk control scheme is reasonable, which not only ensures the construction process but also guarantees the stability of the slope and unstable rock body. At the same time, the slope is reinforced with fiber concrete, which effectively decreases the protection wall thickness. Finally, the article can provide a valuable reference for similar engineering projects around the world.


Introduction
Currently, a slope protection wall is usually used, which is poured using conventional Portland cement concrete. The construction safety risk is greatly increased when blasting excavation must be carried out for an engineering project, where it takes place in certain types of terrain, such as a wide, unstable rock belt. However, conventional preventive measures are faced with a series of challenges, such as difficulties in construction, long project periods, high budget expenses and low prediction accuracy. In this regard, many scholars have tried to introduce interdisciplinary approaches to eliminate intractable engineering issues and have made effective progress. For example, by introducing new materials, the dangerous area is accurately set up, and the structural performance is also improved with regard to dangerous area reinforcement [1]. Second, an early warning system with a machine learning algorithm or other computer-aided techniques can be applied. Chen Z et al. [2] have made an early attempt in this field. They decomposed the original measurement signal using the ensemble empirical mode decomposition (EEMD) method, and further combined deep neural networks (DNNs), gated recurrent units (GRUs), and set displacement monitoring points at crucial positions of the foundation pit, there was still uncertainty in determining the key points of the foundation pit, due to the lack of preconstruction scheme simulation. Yang Jichao et al. [22] considered the actual engineering situation of the Qingshuihe Bridge, and further proposed the safety control of the gravity anchorage deep foundation pit of the suspension bridge utilizing layered excavation, bolt hanging net injection, and prestressed anchor cable protection. Although the proposed method effectively prevents the collapse of rock debris in the foundation pit slope, it does not include risk monitoring of the whole safety control process. Results show that there are deficiencies if we only rely on construction personnel to inspect hidden dangers. Xiao Anbin et al. [23] used FLAC3D to conduct a finite element simulation of the pre-construction scheme of the deep foundation pit slope project of the north bank of the Baiyang Yangtze River Highway Bridge. The numerical model established by the finite difference program can obtain the dynamic deformation and stability evolution law of the rock mass of the deep foundation pit slope project, which is conducive to the construction and the determination of the key safety points. However, the limitation was that it only used the station to monitor the displacement in the field, which lagged behind the judgment of the stability of the slope.
To sum up, many scholars have solved the problem of slope instability in the excavation process by taking preventive measures before construction, or the timely detection of safety risks through monitoring during construction. However, few have introduced the methods of other disciplines, using both the finite element simulation of the prevention program, and the real-time monitoring of the construction process of the program to conduct a dynamic stability analysis of the slope [24], and even to eliminate the safety hazards.
The anchorage foundation pit of the suspension bridge is currently facing a series of technical difficulties, such as the consideration of special terrain conditions, blasting difficulties of anchorage excavation, and high risks in construction. Therefore, how to carry out effective risk control is of great significance. Based on the studies of construction challenges and detailed risk control measures, this paper puts forward several comprehensive risk control measures using the method of hierarchical excavation and minimum charge blasting. In addition, the paper verifies the feasibility and influence of the measures, and summarizes their advantages and disadvantages, which can further provide a reference for similar engineering projects around the world. The methodology of this study flow of the full text is shown in Figure 1.
project. In view of similar engineering situations, scholars have conducted some exploration studies. For example, Zhang Hui et al. [21], adopted the method of anchor bolt support and net jetting C20 concrete to protect the slope in the Qipanzhou Yangtze River Highway Bridge. They aimed at resolving the technical difficulties of the excavation of a gravity rock-socketed anchorage deep-foundation pit in complex terrain. Although they set displacement monitoring points at crucial positions of the foundation pit, there was still uncertainty in determining the key points of the foundation pit, due to the lack of preconstruction scheme simulation. Yang Jichao et al. [22] considered the actual engineering situation of the Qingshuihe Bridge, and further proposed the safety control of the gravity anchorage deep foundation pit of the suspension bridge utilizing layered excavation, bolt hanging net injection, and prestressed anchor cable protection. Although the proposed method effectively prevents the collapse of rock debris in the foundation pit slope, it does not include risk monitoring of the whole safety control process. Results show that there are deficiencies if we only rely on construction personnel to inspect hidden dangers. Xiao Anbin et al. [23] used FLAC3D to conduct a finite element simulation of the pre-construction scheme of the deep foundation pit slope project of the north bank of the Baiyang Yangtze River Highway Bridge. The numerical model established by the finite difference program can obtain the dynamic deformation and stability evolution law of the rock mass of the deep foundation pit slope project, which is conducive to the construction and the determination of the key safety points. However, the limitation was that it only used the station to monitor the displacement in the field, which lagged behind the judgment of the stability of the slope.
To sum up, many scholars have solved the problem of slope instability in the excavation process by taking preventive measures before construction, or the timely detection of safety risks through monitoring during construction. However, few have introduced the methods of other disciplines, using both the finite element simulation of the prevention program, and the real-time monitoring of the construction process of the program to conduct a dynamic stability analysis of the slope [24], and even to eliminate the safety hazards.
The anchorage foundation pit of the suspension bridge is currently facing a series of technical difficulties, such as the consideration of special terrain conditions, blasting difficulties of anchorage excavation, and high risks in construction. Therefore, how to carry out effective risk control is of great significance. Based on the studies of construction challenges and detailed risk control measures, this paper puts forward several comprehensive risk control measures using the method of hierarchical excavation and minimum charge blasting. In addition, the paper verifies the feasibility and influence of the measures, and summarizes their advantages and disadvantages, which can further provide a reference for similar engineering projects around the world. The methodology of this study flow of the full text is shown in Figure 1.

Risk Control Measures
The proposed project is the suspension bridge anchor foundation pit. One side is adjacent to the high and steep slope, and the other side is adjacent to the Yangtze River. Its elevation is shown in Figure 2. The suspension bridge in the project uses gravity anchors on both banks. The anchorage area is located in the mid-slope, and the altitude of the ground surface distributes between 210 m and 250 m. The designed height for the foundation pit bottom is 198 m (from the bottom of the anchor blocks). In addition, the area of anchorage foundation pit excavation is 5958 m 2 , the area on the bottom of the foundation pit is 1026 m 2 , and the slope area is 5633 m 2 . Therefore, the total volume of the anchorage foundation pit excavation is 110,400 m 3 . Futhermore, the basement of the anchorage foundation pit is located on a moderately weathered mudstone stratum, and its allowable bearing capacity is not less than 1.35 MPa. After excavating, the height of rock slopes are between 12.0 m and 52.0 m, mainly consisting of mudstone and sandstone, covered with a small amount of rocky soil. The maximum grade of the foundation pit slope is 1:0.3, while the minimum is 1:1, and the highest level for the slope is 6. The schematic design of the foundation pit excavation is shown in Figure 3.
The proposed project is the suspension bridge anchor foundation pit. One side i jacent to the high and steep slope, and the other side is adjacent to the Yangtze Rive elevation is shown in Figure 2. The suspension bridge in the project uses gravity anc on both banks. The anchorage area is located in the mid-slope, and the altitude o ground surface distributes between 210 m and 250 m. The designed height for the f dation pit bottom is 198 m (from the bottom of the anchor blocks). In addition, the ar anchorage foundation pit excavation is 5958 m 2 , the area on the bottom of the founda pit is 1026 m 2 , and the slope area is 5633 m 2 . Therefore, the total volume of the ancho foundation pit excavation is 110,400 m 3 . Futhermore, the basement of the anchorage f dation pit is located on a moderately weathered mudstone stratum, and its allowable b ing capacity is not less than 1.35 MPa. After excavating, the height of rock slopes ar tween 12.0 m and 52.0 m, mainly consisting of mudstone and sandstone, covered w small amount of rocky soil. The maximum grade of the foundation pit slope is 1:0.3, w the minimum is 1:1, and the highest level for the slope is 6. The schematic design o foundation pit excavation is shown in Figure 3.  The bridge sits in complex geographical conditions. The cable support tower ba located not only in a cliff area, but also across the riverside highway. Therefore, ther five main engineering challenges, as follows: (1) Challenges for the oversized anchorage foundation pit area and excessive excava volume; (2) A dangerous rock zone is above the construction area of the anchorage founda pit. Also, solitary boulders are observed above the landslide body. Therefore, ser risks exist in the construction; (3) Vibration-induced influence on the slope stability due to the blasting excavation not be ignored; The proposed project is the suspension bridge anchor foundation pit. One side is adjacent to the high and steep slope, and the other side is adjacent to the Yangtze River. Its elevation is shown in Figure 2. The suspension bridge in the project uses gravity anchors on both banks. The anchorage area is located in the mid-slope, and the altitude of the ground surface distributes between 210 m and 250 m. The designed height for the foundation pit bottom is 198 m (from the bottom of the anchor blocks). In addition, the area of anchorage foundation pit excavation is 5958 m 2 , the area on the bottom of the foundation pit is 1026 m 2 , and the slope area is 5633 m 2 . Therefore, the total volume of the anchorage foundation pit excavation is 110,400 m 3 . Futhermore, the basement of the anchorage foundation pit is located on a moderately weathered mudstone stratum, and its allowable bearing capacity is not less than 1.35 MPa. After excavating, the height of rock slopes are between 12.0 m and 52.0 m, mainly consisting of mudstone and sandstone, covered with a small amount of rocky soil. The maximum grade of the foundation pit slope is 1:0.3, while the minimum is 1:1, and the highest level for the slope is 6. The schematic design of the foundation pit excavation is shown in Figure 3.  The bridge sits in complex geographical conditions. The cable support tower base is located not only in a cliff area, but also across the riverside highway. Therefore, there are five main engineering challenges, as follows: (1) Challenges for the oversized anchorage foundation pit area and excessive excavation volume; (2) A dangerous rock zone is above the construction area of the anchorage foundation pit. Also, solitary boulders are observed above the landslide body. Therefore, serious risks exist in the construction; (3) Vibration-induced influence on the slope stability due to the blasting excavation cannot be ignored; The bridge sits in complex geographical conditions. The cable support tower base is located not only in a cliff area, but also across the riverside highway. Therefore, there are five main engineering challenges, as follows: (1) Challenges for the oversized anchorage foundation pit area and excessive excavation volume; (2) A dangerous rock zone is above the construction area of the anchorage foundation pit. Also, solitary boulders are observed above the landslide body. Therefore, serious risks exist in the construction; (3) Vibration-induced influence on the slope stability due to the blasting excavation cannot be ignored; (4) Significant elevation differences on the ground surface lead to difficulties in selecting the channel during excavation; (5) The anchorage site has a large grade of terrain slope and a deep foundation pit.
According to the above technical challenges, a series of comprehensive risk control measures using layered and graded excavation and minimum charge blasting methods are put forward, to meet the safety requirements of the actual engineering project.

Layered and Graded Excavation
The vertical slope of the anchorage foundation pit is divided into six layers to control the risks caused by difficulties such as the excessive area and volume of the excavation. Meanwhile, the top-bottom technique is applied for earth rock excavation according to the ground elevation. Each main layer of the slope is cut into multi-thickness pieces based on various geological conditions for easy excavation, and the control value for the maximum excavation depth is limited to 3 m. The sequence of the designed excavation is shown in Figures 3 and 4. The different colors (purple to yellowish brown) from bottom to top in Figure 5 represent the stratified excavation, corresponding to layers 1-6.
(4) Significant elevation differences on the ground surface lead to difficulties in selecting the channel during excavation; (5) The anchorage site has a large grade of terrain slope and a deep foundation pit.
According to the above technical challenges, a series of comprehensive risk control measures using layered and graded excavation and minimum charge blasting methods are put forward, to meet the safety requirements of the actual engineering project.

Layered and Graded Excavation
The vertical slope of the anchorage foundation pit is divided into six layers to control the risks caused by difficulties such as the excessive area and volume of the excavation. Meanwhile, the top-bottom technique is applied for earth rock excavation according to the ground elevation. Each main layer of the slope is cut into multi-thickness pieces based on various geological conditions for easy excavation, and the control value for the maximum excavation depth is limited to 3 m. The sequence of the designed excavation is shown in Figures 3 and 4. The different colors (purple to yellowish brown) from bottom to top in Figure 5 represent the stratified excavation, corresponding to layers 1-6.   It is essential to use some method to protect all layers and slopes, such as anchorageshotcrete support, bolt-faced walls, or anchor wire-faced walls. The height of rock slope is between 12.0 m and 52.0 m, mainly consisting of mudstone and sandstone, covered with a small amount of rocky soil. The maximum grade of the foundation pit slope is 1:0.3, while the minimum ratio is 1; the highest level for the slope is 6, and the lowest level is 1, respectively. Furthermore, a hand of height slope separately is 8 m and 10 m. It is essential to use some method to protect all layers and slopes, such as anchorageshotcrete support, bolt-faced walls, or anchor wire-faced walls. The height of rock slope is between 12.0 m and 52.0 m, mainly consisting of mudstone and sandstone, covered with a small amount of rocky soil. The maximum grade of the foundation pit slope is 1:0.3, while the minimum ratio is 1; the highest level for the slope is 6, and the lowest level is 1, respectively. Furthermore, a hand of height slope separately is 8 m and 10 m.

Minimum Charge Blasting Method
Blasting excavation was necessary for the foundation pit due to the high-strength rock body. Therefore, the minimum charge blasting method was conducted for the excavation volume, to reduce the blasting influence on the stability of the slope and unstable rock zone. The minimum dosage explosive method is a method to obtain the total amount of blasting charge by calculating the required blasting load and the maximum amount of blasting charge in a single section. The total amount of blasting charge required by the project is then determined according to the relevant specifications of the excavation of the foundation pit. Thus, the amount of blasting excavation charge will have the least influence on the stability of the foundation pit, that is, the minimum amount of blasting charge required for the excavation project.

Establishment of the Finite Element Model
In the project, conventional mechanical excavation was combined with drilling and blasting techniques for the anchorage foundation pit excavation. The software MADIS/GTS (Geotechnical & Tunnel Analysis System) was applied to conduct the finite element analysis. Firstly, the Mohr-Coulomb elastoplastic material model [25] was used for a direct integral solution in the geotechnical calculation, and the relevant parameters are shown in Table 1. Secondly, in model building, the soil mass was discretized using high-precision 10-node tetrahedral elements, and the grid around the structural elements was encrypted to obtain better convergence and calculation accuracy. The wall protection was simulated by a plate element, and the anchor cable support structure was simulated by an implantable truss and beam element. The size range of the model was X[−130, 130], Y[−130, 130], Z[−70, 120], which was divided into 32,463 nodes and 57,080 elements. Moreover, a linear time history analysis was applied, and the initial geo-stress equilibrium of the analysis model only considered the dead weight stress and ignored the influence of tectonic stress. Since there is no rock mass parameter under saturated water and no slope wetting line under rainstorm conditions in the geological prospecting report, the effects of groundwater and rainfall were not considered in the slope blasting excavation. Finally, the dynamic absorption boundary was adopted in this model, and the viscous boundary was set to consider the influence of energy loss caused by wave energy escape on the dynamic properties of the surrounding rock. The bottom boundary of the model limits the horizontal and vertical displacements of elements, while the side boundary limits the horizontal displacements. A combination of contact element and collision element was used to connect different soil layers and model components. The blasting vibration loading was carried out using the reaction spectrum method. The seven physical parameters of natural weight, saturated compressive strength, Poisson's ratio, saturated compressive strength, shear angle, internal friction angle and elastic modulus of the three rock strata soil samples were obtained through geological exploration data and indoor physical experiments, as shown in Table 1. However, it was necessary to further calculate the relevant parameters of the Mohr-Coulomb constitutive model adopted in the finite element simulation. In this regard, referring to the calculation method of concrete Mohr-Coulomb parameters by Selimir, the parameters required by five constitutive models, namely: cohesion, internal friction angle, volume modulus, shear modulus, and tensile strength were determined according to the soil samples of the three strata. . The soil layers were established, and were divided into 32,463 joints and 57,080 elements. The constraint location was fixed in the bottom area and the other areas were free. The following hypothesis was proposed in terms of modeling: (1) The initial ground stress balance merely considers self-weight stress, but ignores tectonic stress; (2) Current slope blasting excavation overlooked the effect of groundwater and rainfall, due to the lack of rock parameters under saturated moisture conditions, and the slope's phreatic lines as a result of rainstorms. The overall model is shown in Figure 5. The different colors (purple to yellowish brown) from bottom to top in Figure 5 represent the stratified excavation, corresponding to layers 1-6. The sequence of the designed excavation pit is shown in Figure 4. Figure 4 is composed of six figures (a) to (f), representing the layer-by-layer excavation design plan of the foundation pit, from the sixth level (a) to the first level (f), respectively. Table 2 lists the design parameters of the supporting structure. Figure 6 shows the anchor bar and anchor cable of the foundation pit system. Figure 7 illustrates the retaining walls of the foundation pit.     (1) Calculation of the equivalent explosive diameter for the maximum single-blow blasting charge

The Design Parameters of the Supporting Structure
The calculation of the explosive diameter was performed as follows: the total blasting charge is Qmax = 18 kg, the total length is L2 = 1.5 m, the blasting density is ρ = 1000 kg/m 3 , the maximum single-blow blasting charge is 18 ÷ 5 = 3.6 kg; the equal explosive diameter is D = 2 × (3.6 ÷ 1000 ÷ π ÷ 1.5)0.5 = 0.055 m = 55 mm, and the explosive diameter is approximately 55 mm. In practice, three 32 mm charges were lashed horizontally and placed into a borehole with a cross-sectional area consistent with the equal explosive diameter of 55 mm.

Blasting Load
(1) Calculation of the equivalent explosive diameter for the maximum single-blow blasting charge The calculation of the explosive diameter was performed as follows: the total blasting charge is Q max = 18 kg, the total length is L2 = 1.5 m, the blasting density is ρ = 1000 kg/m 3 , the maximum single-blow blasting charge is 18 ÷ 5 = 3.6 kg; the equal explosive diameter is D = 2 × (3.6 ÷ 1000 ÷ π ÷ 1.5)0.5 = 0.055 m = 55 mm, and the explosive diameter is approximately 55 mm. In practice, three 32 mm charges were lashed horizontally and placed into a borehole with a cross-sectional area consistent with the equal explosive diameter of 55 mm.
(2) Blast loading The blast loading can be simplified into a triangular load curve based on the explosive vibration theory [26]. The performance indices of 2 rock emulsion explosives are shown in Table 3.
The applicable calculation method on peak blasting value is as follows.
(1) Detonation pressure The calculation of the coupled charging condition is shown as Equation (1).
The equation of the uncoupled charging condition is shown as Equation (2).  (2) The initial pressure peaks of the shock wave The initial pressure peaks of the detonation wave on the borehole wall under the influence of a single borehole are shown in Equation (3).
P r -the initial pressure peaks of the shock wave in rock body (Pa), ρ r -density of rock body (kg/m 3 ); C er -compression wave velocity of the rock body (m/s).
Equation (3) indicates that pressure peaks on the borehole wall differ when the same explosive explodes in various rocks; P r is enlarged with increasing wave impedance in rock (ρ r C er ). The pressure increases with the enhancement of surrounding rock (smaller level) and vice versa. However, the wave transmission process between the borehole wall and bench slope wall follows the wave attenuation law, and the pressure values of the two are not equal.
(3) Blasting equivalent load The equivalent pressure equation is shown below, indicating the conversion from the pressure peaks on the borehole wall to the vibration circle is (4).
The relative parameters can be calculated according to the characteristics of the blasting simulation, as shown in Table 4. The equivalent pressure on the hypothetical boundary can be approximately calculated by the above calculation process. It also represents the explosive peak value in simulation. The explosive peak value of the anchorage foundation pit excavation can be further calculated by the above result, and the relative parameters are shown in Table 5. The maximum charge of a single section is 18 kg, and the calculated load is 9.88 Mpa, resulting in an equivalent value Pmax = 9.88 MPa of elastic vibration. In addition, the rise and fall times are 3 ms and 12 ms, respectively, and the total blasting time is 1 s.

Damping Calculation
Damping is a main factor that causes energy consumption and reduces the amplitude of structural vibration. Furthermore, as a critical factor in wave propagation, damping is beneficial to the attenuation of shock wave energy. Moreover, the Rayleigh [27] damping coefficient, a common damping factor in explosive calculations, assumes the damping matrix as a linear system. It is the combination between mass matrix and stiffness matrix, as shown in Equation (5): where damping matrix C represents the linear combination between mass matrix M and stiffness matrix K; α is the mass ratio damping coefficient; and β is the stiffness ratio damping coefficient. α and β are calculated by Equation (6): ω i and ω j are two separate reference frequencies of stage i and j, ξ indicates the structural damping ratio.
The value of ξ depends on the structural type, material property, and loading wave (ξ is 0.01 due to the lack of damping ratio in the experiment). By applying eigenvalue analysis for the shock wave, two main vibration periods can be obtained: the first main vibration period is 1.06 s; the second main vibration period is 0.87 s.

Boundary Conditions
In dynamic analysis, conventional boundary conditions would cause enormous deviation by reflecting the wave. Therefore, the viscous boundary conditions proposed by Lysmer and Wassis in 1972 [28] were used. The calculation of the soil damping ratio in the x, y and z directions defines the viscous boundaries, as shown in Equations (7) and (8).

Simulation on the Blasting Vibration Field
Based on the designed blasting requirements, conventional loose control blasting (ϕ = 70 mm) method should be applied, with a distance of 50 m from the dangerous rock belt. Working conditions for the simulation are as follows: the horizontal distance between the first-stage slope (close to the side of the slope top) and the unstable rock body is 50 m, which is consistent with the field measurement data. The maximum single-blow blasting charge is 18 kg.
The Mohr-Coulomb elastoplastic material model with a fixed bottom is selected for geotechnical calculation. The surrounding is set as a viscous boundary to deduce the influence of the reflecting wave. The total explosive charge is controlled by stages 1, 3 and 5, with a maximum one-time blasting of 18 kg. As shown in Figure 8, the biggest detonation pressure is 9.88 Mpa; the rising time is 3 ms, the falling time is 12 ms, and the total duration is 1 s. Figures 9 and 10 show the test point on the foundation pit, and the loading area. The blasting load is located in the first stage of the slope (near the side of the slope), while the monitoring point is located at the midpoint of W1 and W2 on the dangerous rock. Therefore, the calculating process for the horizontal distance between them can be as follows:

Blasting Vibration Velocity on the Basement of Perilous Rock
According to Blasting Safety Regulations (GB 6722-2014) [29], the Sadaovsk formula is used to calculate the allowable range of blasting vibration: R-the allowable safety distance for blasting (m); V-the allowable safety vibration velocity for blasting(cm/s); Q-blasting charge, total charge of simultaneous blasting (kg); K,α-parameters and attenuation coefficients related to terrain and geological conditions between the explosive location and protected object.
The blasting safety procedures take the peak particles vibration velocity and main vibration frequency of the protected object as the blasting vibration criterion, and blasting

Blasting Vibration Velocity on the Basement of Perilous Rock
According to Blasting Safety Regulations (GB 6722-2014) [29], the Sadaovsk formula is used to calculate the allowable range of blasting vibration: R-the allowable safety distance for blasting (m); V-the allowable safety vibration velocity for blasting(cm/s); Q-blasting charge, total charge of simultaneous blasting (kg); K,α-parameters and attenuation coefficients related to terrain and geological conditions between the explosive location and protected object.
The blasting safety procedures take the peak particles vibration velocity and main vibration frequency of the protected object as the blasting vibration criterion, and blasting

Blasting Vibration Velocity on the Basement of Perilous Rock
According to Blasting Safety Regulations (GB 6722-2014) [29], the Sadaovsk formula is used to calculate the allowable range of blasting vibration: R-the allowable safety distance for blasting (m); V-the allowable safety vibration velocity for blasting(cm/s); Q-blasting charge, total charge of simultaneous blasting (kg); K,α-parameters and attenuation coefficients related to terrain and geological conditions between the explosive location and protected object.
The blasting safety procedures take the peak particles vibration velocity and main vibration frequency of the protected object as the blasting vibration criterion, and blasting (1) Gravity balance and displacement clearing should be carried out before slope excavation; (2) The static calculation of excavation and support is essential for stages 6, 5, 4, 3 and 2; (3) Blasting is applied on stages 1, 3 and 5 to retrieve the vibration field.

Blasting Vibration Velocity on the Basement of Perilous Rock
According to Blasting Safety Regulations (GB 6722-2014) [29], the Sadaovsk formula is used to calculate the allowable range of blasting vibration: R-the allowable safety distance for blasting (m); V-the allowable safety vibration velocity for blasting (cm/s); Q-blasting charge, total charge of simultaneous blasting (kg); K,α-parameters and attenuation coefficients related to terrain and geological conditions between the explosive location and protected object.
The blasting safety procedures take the peak particles vibration velocity and main vibration frequency of the protected object as the blasting vibration criterion, and blasting vibration velocity control is carried out in the area that may affect the surrounding houses. From a safety perspective, the vibration velocity of a dangerous rock burst is controlled within 1.0 cm/s, and the blasting vibration speed of civil buildings should be controlled within 0.5 cm/s, to not affect the surrounding buildings or disturb the people. Therefore, the parameters, such as K = 200 and α = 1.8, are determined based on the blasting scheme and experience from the project in the Three Gorges Reservoir area. At the same time, the explosion point is located 50 m outside the unstable rock, the maximum single blasting charge is 18 kg, and the vibration velocity of the dangerous rock base is 1.0 cm/s, according to the Sadaovsk formula.

The Velocity Curve of the Blasting Vibration on the Basement under Perilous Rock
With the finite element simulation, Figures 11-14 show the velocity curve of blasting vibration on the basement under the perilous rock.

The Velocity Curve of the Blasting Vibration on the Basement under Perilous Rock
With the finite element simulation, Figures 11-14 show the velocity curve of blasting vibration on the basement under the perilous rock.
The location of blasting leading to the movement of the vibration curve on the basement of perilous rock at T1 and T2 is obviously opposite to T3, and the total vibration curve is coupled by T1 to T3. Furthermore, in terms of single velocity, T1 and T2 appear to be the peak; meanwhile, T3 is the minimum within 0.21, and the vibration velocity of T3 attains the maximum of 0.72 cm/s within 0.35 s. The total vibration velocity on the basement of the perilous rock reaches a maximum value of 1.19 cm/s within 0.36 s. The vibration velocity gradually decays after attaining the maximum value.

The Velocity Curve of the Blasting Vibration on the Basement under Perilous Rock
With the finite element simulation, Figures 11-14 show the velocity curve of blasting vibration on the basement under the perilous rock.
The location of blasting leading to the movement of the vibration curve on the basement of perilous rock at T1 and T2 is obviously opposite to T3, and the total vibration curve is coupled by T1 to T3. Furthermore, in terms of single velocity, T1 and T2 appear to be the peak; meanwhile, T3 is the minimum within 0.21, and the vibration velocity of T3 attains the maximum of 0.72 cm/s within 0.35 s. The total vibration velocity on the basement of the perilous rock reaches a maximum value of 1.19 cm/s within 0.36 s. The vibration velocity gradually decays after attaining the maximum value.     Figures 15-21 show the vibration velocities of the retaining walls at 200 ms, 300 ms, 500 ms, 800 ms, and 1000 ms throughout the blasting blasting vibration mainly impacts the first-stage slope at 50 ms. The vib attain the largest value of 24.05 cm/s on the first-stage slope at 100 ms. T impacts the edge of the retaining walls of the second-stage slope. Moreov vibration has an influence on all retaining walls after 300 ms. The blasting rising tendency at first, and then decreases, with the values of the vib ranging from 50 ms to 1000 ms. The variation trend for the vibration velo also similar, but slightly lagged compared to the blasting load.   Figures 15-21 show the vibration velocities of the retaining walls at 200 ms, 300 ms, 500 ms, 800 ms, and 1000 ms throughout the blasting s blasting vibration mainly impacts the first-stage slope at 50 ms. The vibr attain the largest value of 24.05 cm/s on the first-stage slope at 100 ms. Th impacts the edge of the retaining walls of the second-stage slope. Moreov vibration has an influence on all retaining walls after 300 ms. The blasting rising tendency at first, and then decreases, with the values of the vibr ranging from 50 ms to 1000 ms. The variation trend for the vibration velo also similar, but slightly lagged compared to the blasting load. The location of blasting leading to the movement of the vibration curve on the basement of perilous rock at T1 and T2 is obviously opposite to T3, and the total vibration curve is coupled by T1 to T3. Furthermore, in terms of single velocity, T1 and T2 appear to be the peak; meanwhile, T3 is the minimum within 0.21, and the vibration velocity of T3 attains the maximum of 0.72 cm/s within 0.35 s. The total vibration velocity on the basement of the perilous rock reaches a maximum value of 1.19 cm/s within 0.36 s. The vibration velocity gradually decays after attaining the maximum value. Figures 15-21 show the vibration velocities of the retaining walls at 50 ms, 100 ms, 200 ms, 300 ms, 500 ms, 800 ms, and 1000 ms throughout the blasting simulation. The blasting vibration mainly impacts the first-stage slope at 50 ms. The vibration velocities attain the largest value of 24.05 cm/s on the first-stage slope at 100 ms. This significantly impacts the edge of the retaining walls of the second-stage slope. Moreover, the blasting vibration has an influence on all retaining walls after 300 ms. The blasting load exhibits a rising tendency at first, and then decreases, with the values of the vibration velocities ranging from 50 ms to 1000 ms. The variation trend for the vibration velocity and load is also similar, but slightly lagged compared to the blasting load.

Blasting Vibration Velocity near the Pit
attain the largest value of 24.05 cm/s on the first-stage slope at 100 ms. This significan impacts the edge of the retaining walls of the second-stage slope. Moreover, the blasti vibration has an influence on all retaining walls after 300 ms. The blasting load exhibit rising tendency at first, and then decreases, with the values of the vibration velocit ranging from 50 ms to 1000 ms. The variation trend for the vibration velocity and load also similar, but slightly lagged compared to the blasting load.

Stability Analysis of the Slope
According to the geological prospecting reports, the ground surface bedroc face of the higher bank mainly consists of mudstone. It is a relatively water-resista tum, with a steep slope of sandstone. Therefore, the bedrock fissure water is co tively poor, and the influence of rainfall is temporarily ignored. The nonlinear time +SRM method in the model was used to calculate the slope under complex terr geological conditions, which has the important advantage of considering the soil tutive relationship and the deformation influence on stress. Furthermore, the wor ing condition on the first-stage slope (near the side of the slope) is simulated, with zontal distance of 50 m away from the perilous rock, and the maximum single-blow ing charge of 18 kg.

Stability Analysis of the Slope
According to the geological prospecting reports, the ground surface bedrock interface of the higher bank mainly consists of mudstone. It is a relatively water-resistant stratum, with a steep slope of sandstone. Therefore, the bedrock fissure water is comparatively poor, and the influence of rainfall is temporarily ignored. The nonlinear time history +SRM method in the model was used to calculate the slope under complex terrain and geological conditions, which has the important advantage of considering the soil constitutive relationship and the deformation influence on stress. Furthermore, the worst blasting condition on the first-stage slope (near the side of the slope) is simulated, with a horizontal distance of 50 m away from the perilous rock, and the maximum single-blow blasting charge of 18 kg.
As shown in Figures 22-26, the potential slip surface of the anchorage foundation pit is judged by the metrics, including overall displacement, slope plastic zones, effective plastic strain and equivalent strain. In addition, the potential slip surface on the rock mass behind the pit is sheared from the foundation pit of the slope foot. The results show that the overall instability of the slope will occur in the weak zone and the stress concentration

Stability Analysis of the Slope
According to the geological prospecting reports, the ground surface bedrock interface of the higher bank mainly consists of mudstone. It is a relatively water-resistant stratum, with a steep slope of sandstone. Therefore, the bedrock fissure water is comparatively poor, and the influence of rainfall is temporarily ignored. The nonlinear time history +SRM method in the model was used to calculate the slope under complex terrain and geological conditions, which has the important advantage of considering the soil constitutive relationship and the deformation influence on stress. Furthermore, the worst blasting condition on the first-stage slope (near the side of the slope) is simulated, with a horizontal distance of 50 m away from the perilous rock, and the maximum single-blow blasting charge of 18 kg.
As shown in Figures 22-26, the potential slip surface of the anchorage foundation pit is judged by the metrics, including overall displacement, slope plastic zones, effective plastic strain and equivalent strain. In addition, the potential slip surface on the rock mass behind the pit is sheared from the foundation pit of the slope foot. The results show that the overall instability of the slope will occur in the weak zone and the stress concentration area (plastic zone), and its soil element will undergo permanent deformation in varying degrees. When the above parts are connected with each other, the overall instability of the slope will appear in an interconnected shear fracture surface.

Experiment on Blasting Vibration
The finite element model data based on risk control measures indicate that the vibration in the unstable rock belt is comparatively small, and slope stability is good. According to the results, construction safety can be ensured. Therefore, the above measures were applied in the actual project, and the field data was monitored in real time. Figure 27 shows the two testing devices for the blasting vibration monitoring experiment. The vibration monitoring test must satisfy professional standards and regulations, such as the safety regulations for blasting and Code for measurement methods of dynamic properties of subsoil. The GNSS slope deformation monitor station is used in monitoring slope deformation, and the GNSS relay station is installed close to the site. Therefore, it is convenient to transmit signals to a central system, which can achieve persistent and unattended data monitoring.

Measuring Point Arrangement
(1) Blasting vibration monitoring point In the blasting experiment, three test lines were arranged with three measuring points spaced on the steps, as shown in Figures 28 and 29. The distance between test points and blasting sources is shown in Table 7, and the overall number of blasting times is 6, of which 5 are for the penultimate level of the foundation pit, and 1 time for the lowest level. Based on the strength reduction method, the 2D dynamic safety factor of the slope is 1.38, and the 3D value is 1.78, as shown in Table 6. Table 6. Safety factor of slope stability on anchorage foundation pit.

Calculation Condition Two-Dimension Overall Stability Three-Dimension Overall Stability Normal Safety Factor
Foundation pit dynamic excavation 1.38 1.78 1.25

Experiment on Blasting Vibration
The finite element model data based on risk control measures indicate that the vibration in the unstable rock belt is comparatively small, and slope stability is good. According to the results, construction safety can be ensured. Therefore, the above measures were applied in the actual project, and the field data was monitored in real time. Figure 27 shows the two testing devices for the blasting vibration monitoring experiment. The vibration monitoring test must satisfy professional standards and regulations, such as the safety regulations for blasting and Code for measurement methods of dynamic properties of subsoil. The GNSS slope deformation monitor station is used in monitoring slope deformation, and the GNSS relay station is installed close to the site. Therefore, it is convenient to transmit signals to a central system, which can achieve persistent and unattended data monitoring.

Experiment on Blasting Vibration
The finite element model data based on risk control measures indicate that the vibration in the unstable rock belt is comparatively small, and slope stability is good. According to the results, construction safety can be ensured. Therefore, the above measures were applied in the actual project, and the field data was monitored in real time. Figure 27 shows the two testing devices for the blasting vibration monitoring experiment. The vibration monitoring test must satisfy professional standards and regulations, such as the safety regulations for blasting and Code for measurement methods of dynamic properties of subsoil. The GNSS slope deformation monitor station is used in monitoring slope deformation, and the GNSS relay station is installed close to the site. Therefore, it is convenient to transmit signals to a central system, which can achieve persistent and unattended data monitoring. (1) Blasting vibration monitoring point In the blasting experiment, three test lines were arranged with three measuring points spaced on the steps, as shown in Figures 28 and 29. The distance between test points

Measuring Point Arrangement
(1) Blasting vibration monitoring point In the blasting experiment, three test lines were arranged with three measuring points spaced on the steps, as shown in Figures 28 and 29. The distance between test points and blasting sources is shown in Table 7, and the overall number of blasting times is 6, of which 5 are for the penultimate level of the foundation pit, and 1 time for the lowest level.   (2) Slope deformation monitoring point The slope deformation monitoring points were arranged on all levels of the slope platforms, and the distance between each point varied from 15 m to 30 m. According to the observation marks on the stable stones, the basic monitor point was arranged in a stable area, far away from the monitoring slope. In addition, 9 surface horizontal displacement monitoring points were set for the slope. According to Blasting Safety Regulations (GB 6722-2014) [29], the sampling frequency used to measure the vibration in the field test was 50 Hz, and referring to the allowable range of blasting vibration calculated by the Sadaovsk formula, the maximum change ratio of horizontal displacement on the slope deformation monitoring point was 10 mm/d, and the cumulative amount of change warning value is 50 mm. According to the evaluation method of monitoring and alert value   (2) Slope deformation monitoring point The slope deformation monitoring points were arranged on all levels of the slope platforms, and the distance between each point varied from 15 m to 30 m. According to the observation marks on the stable stones, the basic monitor point was arranged in a stable area, far away from the monitoring slope. In addition, 9 surface horizontal displacement monitoring points were set for the slope. According to Blasting Safety Regulations (GB 6722-2014) [29], the sampling frequency used to measure the vibration in the field test was 50 Hz, and referring to the allowable range of blasting vibration calculated by the Sadaovsk formula, the maximum change ratio of horizontal displacement on the slope deformation monitoring point was 10 mm/d, and the cumulative amount of change warn-  (2) Slope deformation monitoring point The slope deformation monitoring points were arranged on all levels of the slope platforms, and the distance between each point varied from 15 m to 30 m. According to the observation marks on the stable stones, the basic monitor point was arranged in a stable area, far away from the monitoring slope. In addition, 9 surface horizontal displacement monitoring points were set for the slope. According to Blasting Safety Regulations (GB 6722-2014) [29], the sampling frequency used to measure the vibration in the field test was 50 Hz, and referring to the allowable range of blasting vibration calculated by the Sadaovsk formula, the maximum change ratio of horizontal displacement on the slope deformation monitoring point was 10 mm/d, and the cumulative amount of change warning value is 50 mm. According to the evaluation method of monitoring and alert value stipulated in "Technical Standard for Monitoring Construction Foundation Pit Engineering" (GB50497-2019) (pp. 33-35) [30]; also combined with the characteristics of this foundation pit excavation project, the displacement monitoring, and the alert value range is 40-50 mm. However, the displacement alert value of the current project is the maximum accumulated change of horizontal displacement at the monitoring point of the foundation pit slope. Therefore, the final displacement alert value is 50 mm, and the alert value of the maximum displacement change rate is 10 mm/d. The arrangement of the test points is shown in Figure 30, and the construction site is shown in Figure 31. nsors 2022, 22, x FOR PEER REVIEW 40-50 mm. However, the displacement alert value of the current projec accumulated change of horizontal displacement at the monitoring poin pit slope. Therefore, the final displacement alert value is 50 mm, and th maximum displacement change rate is 10 mm/d. The arrangement o shown in Figure 30, and the construction site is shown in Figure 31.   40-50 mm. However, the displacement alert value of the current projec accumulated change of horizontal displacement at the monitoring point pit slope. Therefore, the final displacement alert value is 50 mm, and the maximum displacement change rate is 10 mm/d. The arrangement of shown in Figure 30, and the construction site is shown in Figure 31.

Monitoring Data of Blasting Vibration Velocity
In total, the blasting experiment monitored 17 groups of effective data that contain vibration amplitude and dominant frequency from three directions. The first blasting experiment comprised four blasting times within two days and collected five sets of effective monitoring data on the penultimate level of the foundation pit. The second test obtained a set of valid monitoring data on the lowest level within another two days. The third explosive work using a blasting charge of 352 kg was completed in one day and obtained all the data for points 2 and 3, as shown in Table 8. During the test, the measured vibration was filtered to eliminate the obvious interference vibration (more than three times the designed vibration speed, or less than two times the average minimum vibration speed). Horizontal slope displacement monitoring on the surface is the main task of the project, which can last for one year. The monitoring data show that the displacement of the slope surface had a steady growth trend, and it increased with the foundation pit slope excavation depth. Furthermore, the peak of the horizontal displacement in monitoring point 4 was approximately 12 mm, which was 24% of the control value. Nevertheless, the cumulant displacement index was less than the alert value, and the maximum ratio of displacement change was about 4 mm/d, which was less than the alert value of 10 mm/d. Moreover, the monitor data for other points were far lower than the alert value, which meets the requirements of the designed construction. The horizontal displacement monitoring curve is shown in Figure 32

Comparison between Experiment and Simulation
The data of the finite element simulation are shown in Section 3.3, and the data of the field measurements are shown in Section 4.2, see Table 9. In the field measurement data, the maximum vibration velocity of the dangerous rock foundation is 54.45% of that of the finite element simulation, while at the location of the foundation pit retaining structure, the ratio value is 74.25%. The results of the finite element simulation show that the construction data can ensure the construction safety of the site. However, the data obtained from the simulation has higher values than the field measurement data. The reason is that not all conditions can be fully simulated, thus the most unfavorable conditions are considered, and the numerical calculation is more conservative. In conclusion, the finite element simulation in this paper is effective and is useful to the actual engineering project.

Conclusions
This paper aims at the challenges of the foundation pit excavation of a suspension bridge project in the Three Gorges Reservoir area. A series of comprehensive risk control measures have been put forward by applying a hierarchical excavation and minimum charge blasting method. Meanwhile, the feasibility of the scheme is also verified via a finite element simulation, and further applied to the actual construction. There are five main conclusions as follows: (1) The layered and graded excavation scheme can effectively reduce slope deformation, which is critical for construction safety. Therefore, excavation should follow the construction sequence, and the minimum explosive blasting method should be adopted in the excavation of the foundation pit. Thus, construction safety could be ensured below the dangerous rock belt. (2) The blasting is located on the first slope, close to the side of the slope top, with an approximate horizontal distance of 50 m from the blasting location to the unstable rock mass. The maximum single-blow blasting charge is 18 kg. Meanwhile, the Horizontal slope displacement monitoring is the main work, which lasts for one year. The monitoring data show that the slope displacement increases with the increase of depth, but the increasing trend is stable. In addition, the peak horizontal displacement of monitoring point 4 was about 12 mm, and when the cumulant displacement index was less than the alarm value, it reached 24% of the control value, and the maximum ratio of displacement change was about 4 mm/d, less than the alarm value 10 mm/d. The monitoring data of other points are much lower than the warning value, which meets the requirements of design and construction. The horizontal displacement monitoring curves are shown in Figure 32

Comparison between Experiment and Simulation
The data of the finite element simulation are shown in Section 3.3, and the data of the field measurements are shown in Section 4.2, see Table 9. In the field measurement data, the maximum vibration velocity of the dangerous rock foundation is 54.45% of that of the finite element simulation, while at the location of the foundation pit retaining structure, the ratio value is 74.25%. The results of the finite element simulation show that the construction data can ensure the construction safety of the site. However, the data obtained from the simulation has higher values than the field measurement data. The reason is that not all conditions can be fully simulated, thus the most unfavorable conditions are considered, and the numerical calculation is more conservative. In conclusion, the finite element simulation in this paper is effective and is useful to the actual engineering project.

Conclusions
This paper aims at the challenges of the foundation pit excavation of a suspension bridge project in the Three Gorges Reservoir area. A series of comprehensive risk control measures have been put forward by applying a hierarchical excavation and minimum charge blasting method. Meanwhile, the feasibility of the scheme is also verified via a finite element simulation, and further applied to the actual construction. There are five main conclusions as follows: (1) The layered and graded excavation scheme can effectively reduce slope deformation, which is critical for construction safety. Therefore, excavation should follow the construction sequence, and the minimum explosive blasting method should be adopted in the excavation of the foundation pit. Thus, construction safety could be ensured below the dangerous rock belt. (2) The blasting is located on the first slope, close to the side of the slope top, with an approximate horizontal distance of 50 m from the blasting location to the unstable rock mass. The maximum single-blow blasting charge is 18 kg. Meanwhile, the nonlinear time-history analysis method is used to simulate the vibration field of the anchorage foundation pit slope. As a result, the maximum value of the total velocity can reach 1.19 cm/s within 0.36 s, and the vibration velocity gradually decays after attaining the maximum value. (3) Blasting of first-stage slope has a significant influence on the edge retaining walls of second-stage slope, and the vibration velocities attain 24.05 cm/s in 100 ms. Furthermore, the protecting wall should be comprehensively determined by safety regulations for blasting (concrete), numerical simulation results, and field measurement data. (4) The blasting is located on the first slope (close to the side of the slope top), with a maximum single blasting charge of 18 kg. The two-dimensional dynamic safety factor of slope is 1.38, while the three-dimensional factor is 1.78 based on the strength reduction method. Both safety factors of the excavation slope meet the design requirements of the blasting excavation process. (5) The foundation pit project of the suspension bridge mentioned in this paper has been successfully completed. The results show that risk control measures can guarantee engineering qualities and ensure construction safety. The engineering experience can be set as a reference for similar projects around the world.

Institutional Review Board Statement:
The study did not require ethical approval.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.