# Numerical Simulation of Seepage and Deformation in Excavation of a Deep Foundation Pit under Water-Rich Fractured Intrusive Rock

Authors: Changfeng Yuan; Zhenhui Hu; Zhen Zhu; Zijin Yuan; Yanxiang Fan; Hui Guan; Liang Li
Journal: Geofluids
 (2021)

---

## Abstract

Water is one of the major risk sources in the excavation of deep-large foundation pits in a water-rich area. The presence of intrusive broken diorite porphyrite in the stratum aggravates the risk level of deep foundation pits. Based on a geological survey report and design documents of parameter information, MIDAS/GTS software was used to perform the numerical simulation of an engineering example of a deep foundation pit project of ultradeep and water-rich intrusion into the broken rock station of subway line 4 in a city. The simulation results show the characteristics of seepage path evolution, seepage aggregation areas and points, and the effect of seepage on the deformation of a deep foundation pit during the whole construction of this deep foundation pit. The results show that with the precipitation-excavation of the deep foundation pit, the pore water pressure at the bottom of the foundation pit follows a distribution of three “concave” shapes. High-permeability pressure zones are found around the foundation pit, intruding broken diorite porphyrite zones, and middle coarse sand zones. With further excavation of the foundation pit, the seepage pressure in the middle part of the foundation pit gradually decreases, and the two “concave” distributions in the middle gradually merge together. After excavation to the bottom of the pit, the pore water pressure at the bottom is distributed in two asymmetrical “concave” shapes, and the maximum peak of pore water pressure is found at the intrusion of fractured porphyrites prone to water inrush. The four corners of the foundation pit are prone to form seepage accumulation zones; therefore, suffosion and piping zones are formed. The surface settlement caused by excavation is found to be the largest along the longitudinal axis of the deep foundation pit, whereas the largest deformation is found near the foundation pit side in the horizontal axis direction of the foundation pit. With the excavation of the deep foundation pit, the diaphragm wall converges to the foundation pit with the maximum deformation reaching about 25 mm. After the first precipitation-excavation of the deep foundation pit to the silty clay and the bottom of the pit with the largest uplift, with further precipitation-excavation of the deep foundation pit, the uplift at the bottom of the deep foundation pit changes only slightly.


---

## Body

## 1. Introduction

Disasters of deep foundation pits are mostly related to the improper control of groundwater. In recent years, with the shortage of urban land resources, more and more deep foundation pit projects have been undertaken, and the excavation depth has generally reached more than 20 m. At the same time, the number of deep foundation pits larger than 30 m is also increasing. Water problems are encountered as the foundation pit continues to deepen; in particular, groundwater control is the focus of research in some deep foundation pit projects when located in water-rich fractured rock areas. At present, mainly two methods are available to control confined water, i.e., water stop curtain partition and pumping decompression [1]. Both the abovementioned methods require a water seepage path and evolution trend along with excavation to obtain accurate information about the design to ensure project safety. Regarding the seepage and deformation characteristics of deep foundation pits, studies [2–28] have shown that the seepage velocity of a foundation pit first increases and then decreases with the increase in the waterproof curtain depth and finally becomes stable. The hydraulic gradient is distributed at the bottom of the pit in a “concave” shape. Hydraulic gradient is the largest at the corner of the foundation pit, where it is prone to seepage failure. For the excavation of deep foundation pits under the influence of precipitation and excess pore water pressure, the change in permeability coefficient significantly affects the foundation pit floor heave, and the change in permeability coefficient significantly affects the foundation pit. The optimum bottom heave insertion ratio of the underground diaphragm wall is approximately 0.5; the maximum performance can be reached under this condition. If the diaphragm wall is too deep, the effect on the uplift resistance of the foundation pit will be reduced. A linear relationship is obtained between the maximum horizontal displacement of the underground diaphragm wall and excavation depth. As the depth of phreatic water level increases, the maximum surface settlement and maximum pile displacement on both sides of the pit will decrease linearly. In the bottom layer of foundation pit excavation under seepage condition, the soil around the pile will also rebound due to water pressure. The soil’s rebound amount is higher than the pile’s rebound amount due to the water pressure caused by large displacement rebound of soil around the pile. Although great achievements have been made in this research, the variability of geotechnical engineering in time and space makes each specific deep foundation pit project have its own characteristics; therefore, it is necessary to perform focused research in the abovementioned aspects for deep large foundation pits with different geological conditions.The seepage and deformation evolution trends of deep foundation pit in the whole construction process were studied through a numerical simulation for the excavation of a deep foundation pit of a 37 m deep city subway station with water-rich and locally broken intrusive rocks. The research results provide a reference for the excavation design and construction of deep foundation pits with water-rich and locally invaded broken diorite porphyritic stratum.


## 2. General Situation of the Project

The west side of the deep foundation pit of the line 4 subway station is adjacent to the Daedong River, and a scenic plaza and parking lot are present on the east side. This subway station is a 13 m island platform, double-column three-span station with four underground levels. The station is 158 m long, 26 m wide, and 37 m deep, and it is constructed using the open bottom-up method. The security level of the foundation pit is grade 1.


### 2.1. Engineering Geological Conditions

According to the survey report, the pit area is mainly composed of a Holocene artificial accumulation layer, marine sedimentary layer, and alluvial layer. The underlying bedrock is late Yanshan granite, locally interspersed with diorite porphyrite, and the site area is still tectonically developing.The artificial accumulation layer of the Fourth Series Holocene is filled with vegetation, and the marine sediment layer is mainly made of organic silt. These two layers have a loose structure and very poor homogeneity, poor self-stabilization ability after excavation, poor engineering properties, and strong permeability.The alluvial-diluvial layer mainly consists of silty clay, organic silty clay, coarse sand, coarse gravel sand, and pebble. Silty clay is mainly malleable and partially malleable. Organic silty clay has high compressibility and low bearing capacity. The distribution of coarse sand is discontinuous and saturated. The bearing capacity is relatively low, and the excavation is prone to collapse and water gushing, coarse gravel sand saturation, and high bearing capacity. The bearing capacity of silty clay is relatively low. The pebble and bedrock have direct contact, saturation, high bearing capacity, and low compressibility.The granites in the bedrock lie in the tectonic fracture zone and its influencing range; the borehole mostly exposes the sandy soil-like and massive cataclastic rocks. The bedrock of intrusive broken diorite porphyrite is affected by tectonic action joint development and crushing of rock mass; the joint surface mineral alteration is obvious, so the Aeolian zone rock mass quality classification is IV level.The geological profile along this deep foundation pit is shown in Figure1. It can be seen that at the beginning 0 m, the elevation of medium coarse sand is -2 m. Near 40 m along the deep foundation pit, there are intrusive fractured porphyrites with an elevation of -33 m, and between 100 m and 110 m along the deep foundation pit, the elevation of medium coarse sand is -13 m.Figure 1
Geological profile along the deep foundation pit.


### 2.2. Hydrogeological Condition

The site is located in the lower basin of the Daedong River, and the width of the river is about 80-90 m. The bottom level of the river is about 0.77-1.53 m, and masonry embankments are present on both sides of the river, as shown in Figure2. The river is a rain-sourced mountain stream-type river in the monsoon area, which is prone to flash floods during a heavy precipitation season. The critical daily precipitation of Daedong River flash flood was about 100 mm, and the recurrence interval of flash flood is 2-3 years. Figure 1 shows that the water level of the deep foundation pit is 1.5-2.5 m below the ground.Figure 2
Daedong River.The main types of groundwater in the pit area are quaternary pore water and bedrock fissure water. Quaternary pore water is divided into upper stagnant water and confined water; bedrock fissure water is mainly divided into weathered fissure water and structural fissure water. The upper stagnant water is mainly artificially filled with a soil layer, shallow silt, and sandy soil layer. The confined water aquifer is mainly composed of a sandy soil and gravel soil layer and also includes several clay-water resisting layer. Weathering fissure water is mainly stored in the bedrock’s strong weathering to medium weathering zone rocks, which are sandy-soil, sandy, and gravelly. The weathered fissure development shows a layered distribution. The water surface is uniform, and the underground water level increases with the increase in the terrain. Because the bedrock’s strong weathering zone and the upper quaternary aquifer have no stable continuous water barrier, the weathered fissure water has a certain hydraulic relationship with the pore water of a loose quaternary soil layer. The tectonic fissure water is mainly stored in the tectonic fissure zone of medium-weathered and slightly weathered bedrock and the dense zone of an extrusion fracture of late intruded vein-like rock veins, where it has a large depth of groundwater runoff, complex direction of runoff, and certain bearing pressure.


### 2.3. Design of Deep Foundation Pit Enclosure Structure

According to the design document, the deep foundation pit retaining system consists of three parts, as shown in Figure3. The underground diaphragm wall and internal support are used in the first part of the ground to 15 m underground. The second part uses steel pipe piles and prestressed anchors from 15 m underground to 37 m at the bottom of the pit. In the third part, a concrete uplift pile and concrete cushion are applied at the bottom of the foundation pit. The specific design is as follows:Figure 3
Enclosure system of the deep foundation pit.Part I of the Enclosure System. The first one uses a concrete inner support and crown beam, and the remaining four inner supports use a double 56C I-beam and an 800 mm steel pipe. The foundation pit column uses a 1400mm×1000mm concrete pile. The depth of the underground diaphragm wall is 15 m and penetrates the underground rock mass for about 3-5 m. The waterproof curtain lies on the underground diaphragm wall.Part II of the Enclosure System. After the construction of the underground diaphragm wall was completed, a 25 m long steel pipe pile with a diameter of 200 mm was embedded in the rock with a spacing of 900 mm. After the construction of steel pipe piles was completed, the prestressing anchors were started. Eleven layers of prestressing anchors were installed with a vertical spacing of 2 m and lateral spacing of 1.8 m. Then, 700 kN of prestressing was applied to the 1-4 layers, 500 kN to the 5-8 layers, and 200 kN to the 9-11 layers.Part III of the Enclosure System. After the excavation of the foundation pit was completed, eight rows of concrete uplift piles with a lateral spacing of 3 m, a length of 4 m, and a diameter of 400 mm were placed in the bottom of the foundation pit, finally laying a 200 mm thick concrete bedding layer.


## 2.1. Engineering Geological Conditions

According to the survey report, the pit area is mainly composed of a Holocene artificial accumulation layer, marine sedimentary layer, and alluvial layer. The underlying bedrock is late Yanshan granite, locally interspersed with diorite porphyrite, and the site area is still tectonically developing.The artificial accumulation layer of the Fourth Series Holocene is filled with vegetation, and the marine sediment layer is mainly made of organic silt. These two layers have a loose structure and very poor homogeneity, poor self-stabilization ability after excavation, poor engineering properties, and strong permeability.The alluvial-diluvial layer mainly consists of silty clay, organic silty clay, coarse sand, coarse gravel sand, and pebble. Silty clay is mainly malleable and partially malleable. Organic silty clay has high compressibility and low bearing capacity. The distribution of coarse sand is discontinuous and saturated. The bearing capacity is relatively low, and the excavation is prone to collapse and water gushing, coarse gravel sand saturation, and high bearing capacity. The bearing capacity of silty clay is relatively low. The pebble and bedrock have direct contact, saturation, high bearing capacity, and low compressibility.The granites in the bedrock lie in the tectonic fracture zone and its influencing range; the borehole mostly exposes the sandy soil-like and massive cataclastic rocks. The bedrock of intrusive broken diorite porphyrite is affected by tectonic action joint development and crushing of rock mass; the joint surface mineral alteration is obvious, so the Aeolian zone rock mass quality classification is IV level.The geological profile along this deep foundation pit is shown in Figure1. It can be seen that at the beginning 0 m, the elevation of medium coarse sand is -2 m. Near 40 m along the deep foundation pit, there are intrusive fractured porphyrites with an elevation of -33 m, and between 100 m and 110 m along the deep foundation pit, the elevation of medium coarse sand is -13 m.Figure 1
Geological profile along the deep foundation pit.


## 2.2. Hydrogeological Condition

The site is located in the lower basin of the Daedong River, and the width of the river is about 80-90 m. The bottom level of the river is about 0.77-1.53 m, and masonry embankments are present on both sides of the river, as shown in Figure2. The river is a rain-sourced mountain stream-type river in the monsoon area, which is prone to flash floods during a heavy precipitation season. The critical daily precipitation of Daedong River flash flood was about 100 mm, and the recurrence interval of flash flood is 2-3 years. Figure 1 shows that the water level of the deep foundation pit is 1.5-2.5 m below the ground.Figure 2
Daedong River.The main types of groundwater in the pit area are quaternary pore water and bedrock fissure water. Quaternary pore water is divided into upper stagnant water and confined water; bedrock fissure water is mainly divided into weathered fissure water and structural fissure water. The upper stagnant water is mainly artificially filled with a soil layer, shallow silt, and sandy soil layer. The confined water aquifer is mainly composed of a sandy soil and gravel soil layer and also includes several clay-water resisting layer. Weathering fissure water is mainly stored in the bedrock’s strong weathering to medium weathering zone rocks, which are sandy-soil, sandy, and gravelly. The weathered fissure development shows a layered distribution. The water surface is uniform, and the underground water level increases with the increase in the terrain. Because the bedrock’s strong weathering zone and the upper quaternary aquifer have no stable continuous water barrier, the weathered fissure water has a certain hydraulic relationship with the pore water of a loose quaternary soil layer. The tectonic fissure water is mainly stored in the tectonic fissure zone of medium-weathered and slightly weathered bedrock and the dense zone of an extrusion fracture of late intruded vein-like rock veins, where it has a large depth of groundwater runoff, complex direction of runoff, and certain bearing pressure.


## 2.3. Design of Deep Foundation Pit Enclosure Structure

According to the design document, the deep foundation pit retaining system consists of three parts, as shown in Figure3. The underground diaphragm wall and internal support are used in the first part of the ground to 15 m underground. The second part uses steel pipe piles and prestressed anchors from 15 m underground to 37 m at the bottom of the pit. In the third part, a concrete uplift pile and concrete cushion are applied at the bottom of the foundation pit. The specific design is as follows:Figure 3
Enclosure system of the deep foundation pit.Part I of the Enclosure System. The first one uses a concrete inner support and crown beam, and the remaining four inner supports use a double 56C I-beam and an 800 mm steel pipe. The foundation pit column uses a 1400mm×1000mm concrete pile. The depth of the underground diaphragm wall is 15 m and penetrates the underground rock mass for about 3-5 m. The waterproof curtain lies on the underground diaphragm wall.Part II of the Enclosure System. After the construction of the underground diaphragm wall was completed, a 25 m long steel pipe pile with a diameter of 200 mm was embedded in the rock with a spacing of 900 mm. After the construction of steel pipe piles was completed, the prestressing anchors were started. Eleven layers of prestressing anchors were installed with a vertical spacing of 2 m and lateral spacing of 1.8 m. Then, 700 kN of prestressing was applied to the 1-4 layers, 500 kN to the 5-8 layers, and 200 kN to the 9-11 layers.Part III of the Enclosure System. After the excavation of the foundation pit was completed, eight rows of concrete uplift piles with a lateral spacing of 3 m, a length of 4 m, and a diameter of 400 mm were placed in the bottom of the foundation pit, finally laying a 200 mm thick concrete bedding layer.


## 3. Seepage Field and Deformation Characteristics of the Deep Foundation Pit

Currently, urban underground projects in China are more often designed using MIDAS software, which is developed in the Visual C++ language of the Windows system. It is a software program for the stress and infiltration analysis and design in the construction phase of geotechnical structures. In this study, MIDAS/GTS software was used for simulation analysis, and the geotechnical parameters and material properties used in the simulation were provided from the design file of the deep foundation pit, so that the simulation results can reflect the seepage and deformation characteristics of the deep foundation pit under the design of the deep foundation pit and better guide the design.


### 3.1. MIDAS Numerical Model Simulation

#### 3.1.1. Basic Theory of MIDAS Seepage Analysis

MIDAS seepage analysis includes steady-flow analysis and unsteady-flow analysis [29]. In this study, unsteady-flow analysis was used to establish the model. The flow law used in MIDAS/GTS is Darcy’s law, and the formula of penetration flow per unit area is as follows:
(1)q=ki,where q is the seepage discharge per unit area, k is the permeability coefficient, and i is the hydraulic gradient.It is assumed that the change in the inflow and outflow of a small volume at any position and at any time is the same as the change in the volume moisture content. Simply put, the change of flow in thex, y, and z directions is the same as the change of the sum of the flows and the volumetric moisture content.The basic equation for flow is as follows:(2)∂∂xkx∂H∂x+∂∂yky∂H∂y+∂∂zkz∂H∂z+Q=∂Θ∂t,where H is the total head; kx,ky, and kz are the permeability coefficient in the x, y, and z directions, respectively; Q is the flow rate; Θ is the volumetric water content; and t is the time.According to Darcy’s law, the laminar velocity is very small, so the velocity head can be ignored. The total head can be regarded as the sum of pressure head and position head.(3)H=uwγw+z,where uw is the pore water pressure, γw is the unit weight of water, and z is the elevation.


#### 3.1.2. Establishment of the Numerical Model

The construction phase of deep foundation pits was modelled using the MIDAS/GTS software. The actual excavation part of the deep foundation pit in this project is 158 m long, 26 m wide, and 37 m deep. Considering the boundary conditions, the whole numerical model is 299 m long, 167 m wide, and 55 m deep. There are 125248 nodes and 248745 elements in the whole model. The physical and mechanical parameters of each soil layer are shown in Table1. At present, the pit is under construction and has been excavated to 8 m. Combined with the site monitoring data, the geotechnical parameters shown in Table 1 are consistent with the site monitoring data. The numerical model is shown in Figure 4. The Mohr-Coulomb constitutive model was used in the soil layer of modulus value type. Crown beam, lumbar beam, and steel internal support all used the beam element. This is because the displacement of the uplift pile and pit column in the vertical direction is restrained, and their force characteristics are similar to those of the beam. Therefore, they also use the beam element. The foundation pit floor and underground diaphragm wall were simulated with the plate element. When establishing the waterproof curtain on the underground diaphragm wall because the waterproof curtain belongs to the interface element which is equivalent to cutting the nodes, it will cause the nodes to be uncoupled. Therefore, before the underground diaphragm wall is activated with the waterproof curtain, a rigid connection should be established to bind the cut nodes together to make the model coupled as a whole body. Prestressing anchors were used using the implanted truss element.Table 1
Physical and mechanical parameters of surrounded layers.
Soil layerThickness (m)Elastic modulus (MPa)Poisson’s ratioμγ (kN·m−3)Permeability coefficient (m·d)Cohesion (kPa)Internal friction angleφ (°)Plain fill50.30.1718301020Silty clay4500.33200.0234.916.5Coarse sand7180.320.5150.130Fractured rock10120000.2250.5500050Granite30300000.1260.0013000060Figure 4
MIDAS/GTS numerical simulation diagram.(1) Boundary Conditions of Model Displacement. The normal displacement around the model and surface is constrained, and the top surface is the free surface. The seepage boundary condition is the impermeable boundary at the place where the continuous underground diaphragm wall is precipitated around the foundation pit. The permeability coefficient is 0 m/d at the place where the waterproof curtain is set. However, in the excavation of the pit to 15 m below, the excavation was carried out deep into the rock layer. Because the permeability coefficient is very small, the actual project did not set the waterproof curtain. However, the project section has the possibility of infiltration of groundwater outside the pit. For example, under about 10 m underground, a section of sandy fractured rock with a large permeability coefficient invaded the granite interior. This section is very likely to allow underground water from outside the pit flow into the pit.


#### 3.1.3. Deep Foundation Pit Construction Phase

The simulation of the deep foundation pit construction phase is shown in Figure5. The excavation of this deep foundation pit is divided into six steps. The precipitation is divided into six steps. Five supports were applied. Prestressed anchors were applied after the excavation of the 4th step of precipitation, and finally, the uplift piles and the foundation slab were applied.Figure 5
Deep foundation pit construction phase.


### 3.2. Seepage Field Characteristics

The seepage path and corresponding pore water pressure of the deep foundation pit after six precipitation events were extracted, as shown in Figure6. The seepage path of the deep foundation pit after each precipitation is different. Intrusive broken diorite porphyrite and medium coarse sand are the key points and belts for waterproofing the deep foundation pit. With the precipitation-excavation of the deep foundation pit, the pore water pressure at the bottom of the foundation pit shows a distribution of three “concave” shapes. A high-permeability pressure zone occurs around the foundation pit, the intrusive broken diorite porphyrite zone, and the middle coarse sand zone. With further excavation of the foundation pit, the seepage pressure in the middle of the foundation pit gradually decreases, and the two “concave” shapes in the middle of the foundation pit gradually merge together. After excavation to the bottom of the foundation pit, the pore water pressure at the bottom of the foundation pit shows two asymmetrical “concave” shapes, and the peak value is located at the place where broken diorite porphyrite is intruded, which is prone to water inrush. The four corners of the foundation pit are prone to seepage accumulation and subduction corrosion and finally form the piping zone.Figure 6
Seepage path and pore water pressure evolution of the deep foundation pit after precipitation.
(a)
Seepage path and pore water pressure distribution at the bottom of the pit of the 1st precipitation(b)
Seepage path and pore water pressure distribution at the bottom of the pit of the 2nd precipitation(c)
Seepage path and pore water pressure distribution at the bottom of the pit of the 3rd precipitation(d)
Seepage path and pore water pressure distribution at the bottom of the pit of the 4th precipitation(e)
Seepage path and pore water pressure distribution at the bottom of the pit of the 5th precipitation(f)
Seepage path and pore water pressure distribution at the bottom of the pit of the 6th precipitationFigure6(a) shows the seepage path of the 1st precipitation and the distribution of pore water pressure at the bottom of the pit. Three seepage gathering points are present in the deep foundation pit. Two of them are at the corners of the pit on both sides of the pit, and the stratum above them is medium-coarse sand. The other one is intrusive broken diorite porphyrite area. Although there are three seepage accumulation points, the height of the foundation pit seepage path at this time only reaches the silty clay layer, and the pore water pressure of the foundation pit floor shows negative pressure.Figure6(b) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. The height of the seepage path of the deep pit only reaches the silted clay layer, and the pore water pressure at the bottom of the pit slightly varies with the positive pore water pressure at the bottom of the pit.Figure6(c) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. When the deep foundation pit is excavated to about 15 m underground, i.e., when the diaphragm wall and steel pipe pile are converted, the maximum pore water pressure at the bottom of the pit is 135 kN/m2. It is highly recommended to make waterproof work before continuing the construction.Figure6(d) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. The seepage path of the deep foundation pit became more and more concentrated, and the pore water pressure is distributed in 3 “concave” shapes at the bottom of excavation. High-permeability pressure zones are present around the foundation pit, the intrusive broken diorite porphyrite zone, and the middle coarse sand zone.Figure6(e) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. Although the pore water pressure is still distributed in three “concave” shapes at the bottom of the pit at the excavation, the two “concave” shapes between the intrusive broken diorite porphyrite and the terminal mile of pit tend to merge together.Figure6(f) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. After the end of precipitation, the pore water pressure at the bottom of the pit showed two very obvious “concave” shapes, mainly due to the presence of intrusive broken diorite porphyrite in the deep pit, which changed the groundwater flow path of pit excavation and showed multiple “concave” patterns.


### 3.3. Deformation Characteristics

Figure7 shows the deformation evolution of deep foundation pit precipitation. The surface settlement caused by excavation is the largest along the longitudinal axis of the deep foundation pit in the area invaded by the broken zone, and the largest deformation along the horizontal axis of the foundation pit was found near the foundation pit side. After the excavation of the foundation pit entered the moderately weathered rock stratum, the surface settlement became stable. With the excavation of the deep foundation pit, the convergence deformation of the underground diaphragm wall to the pit occurred, with a maximum value of about 25 mm. The change in the underground diaphragm wall was only above the construction of the medium weathering rock layer. With the construction of steel pipe piles, the deformation of the underground diaphragm wall became stable. After the first precipitation excavation of the deep foundation pit to the powder clay, the bottom of the pit had an uplift to the maximum. With further precipitation-excavation, the uplift of the bottom had only a small change.Figure 7
Deformation evolution of deep foundation pit precipitation excavation.
(a)
Deformation characteristics of the 1st precipitation excavation(b)
Deformation characteristics of the 1st precipitation excavation(c)
Deformation characteristics of the 1st precipitation excavation(d)
Deformation characteristics of the 1st precipitation excavation(e)
Deformation characteristics of the 1st precipitation excavation(f)
Deformation characteristics of the 1st precipitation excavation


## 3.1. MIDAS Numerical Model Simulation

### 3.1.1. Basic Theory of MIDAS Seepage Analysis

MIDAS seepage analysis includes steady-flow analysis and unsteady-flow analysis [29]. In this study, unsteady-flow analysis was used to establish the model. The flow law used in MIDAS/GTS is Darcy’s law, and the formula of penetration flow per unit area is as follows:
(1)q=ki,where q is the seepage discharge per unit area, k is the permeability coefficient, and i is the hydraulic gradient.It is assumed that the change in the inflow and outflow of a small volume at any position and at any time is the same as the change in the volume moisture content. Simply put, the change of flow in thex, y, and z directions is the same as the change of the sum of the flows and the volumetric moisture content.The basic equation for flow is as follows:(2)∂∂xkx∂H∂x+∂∂yky∂H∂y+∂∂zkz∂H∂z+Q=∂Θ∂t,where H is the total head; kx,ky, and kz are the permeability coefficient in the x, y, and z directions, respectively; Q is the flow rate; Θ is the volumetric water content; and t is the time.According to Darcy’s law, the laminar velocity is very small, so the velocity head can be ignored. The total head can be regarded as the sum of pressure head and position head.(3)H=uwγw+z,where uw is the pore water pressure, γw is the unit weight of water, and z is the elevation.


### 3.1.2. Establishment of the Numerical Model

The construction phase of deep foundation pits was modelled using the MIDAS/GTS software. The actual excavation part of the deep foundation pit in this project is 158 m long, 26 m wide, and 37 m deep. Considering the boundary conditions, the whole numerical model is 299 m long, 167 m wide, and 55 m deep. There are 125248 nodes and 248745 elements in the whole model. The physical and mechanical parameters of each soil layer are shown in Table1. At present, the pit is under construction and has been excavated to 8 m. Combined with the site monitoring data, the geotechnical parameters shown in Table 1 are consistent with the site monitoring data. The numerical model is shown in Figure 4. The Mohr-Coulomb constitutive model was used in the soil layer of modulus value type. Crown beam, lumbar beam, and steel internal support all used the beam element. This is because the displacement of the uplift pile and pit column in the vertical direction is restrained, and their force characteristics are similar to those of the beam. Therefore, they also use the beam element. The foundation pit floor and underground diaphragm wall were simulated with the plate element. When establishing the waterproof curtain on the underground diaphragm wall because the waterproof curtain belongs to the interface element which is equivalent to cutting the nodes, it will cause the nodes to be uncoupled. Therefore, before the underground diaphragm wall is activated with the waterproof curtain, a rigid connection should be established to bind the cut nodes together to make the model coupled as a whole body. Prestressing anchors were used using the implanted truss element.Table 1
Physical and mechanical parameters of surrounded layers.
Soil layerThickness (m)Elastic modulus (MPa)Poisson’s ratioμγ (kN·m−3)Permeability coefficient (m·d)Cohesion (kPa)Internal friction angleφ (°)Plain fill50.30.1718301020Silty clay4500.33200.0234.916.5Coarse sand7180.320.5150.130Fractured rock10120000.2250.5500050Granite30300000.1260.0013000060Figure 4
MIDAS/GTS numerical simulation diagram.(1) Boundary Conditions of Model Displacement. The normal displacement around the model and surface is constrained, and the top surface is the free surface. The seepage boundary condition is the impermeable boundary at the place where the continuous underground diaphragm wall is precipitated around the foundation pit. The permeability coefficient is 0 m/d at the place where the waterproof curtain is set. However, in the excavation of the pit to 15 m below, the excavation was carried out deep into the rock layer. Because the permeability coefficient is very small, the actual project did not set the waterproof curtain. However, the project section has the possibility of infiltration of groundwater outside the pit. For example, under about 10 m underground, a section of sandy fractured rock with a large permeability coefficient invaded the granite interior. This section is very likely to allow underground water from outside the pit flow into the pit.


### 3.1.3. Deep Foundation Pit Construction Phase

The simulation of the deep foundation pit construction phase is shown in Figure5. The excavation of this deep foundation pit is divided into six steps. The precipitation is divided into six steps. Five supports were applied. Prestressed anchors were applied after the excavation of the 4th step of precipitation, and finally, the uplift piles and the foundation slab were applied.Figure 5
Deep foundation pit construction phase.


## 3.1.1. Basic Theory of MIDAS Seepage Analysis

MIDAS seepage analysis includes steady-flow analysis and unsteady-flow analysis [29]. In this study, unsteady-flow analysis was used to establish the model. The flow law used in MIDAS/GTS is Darcy’s law, and the formula of penetration flow per unit area is as follows:
(1)q=ki,where q is the seepage discharge per unit area, k is the permeability coefficient, and i is the hydraulic gradient.It is assumed that the change in the inflow and outflow of a small volume at any position and at any time is the same as the change in the volume moisture content. Simply put, the change of flow in thex, y, and z directions is the same as the change of the sum of the flows and the volumetric moisture content.The basic equation for flow is as follows:(2)∂∂xkx∂H∂x+∂∂yky∂H∂y+∂∂zkz∂H∂z+Q=∂Θ∂t,where H is the total head; kx,ky, and kz are the permeability coefficient in the x, y, and z directions, respectively; Q is the flow rate; Θ is the volumetric water content; and t is the time.According to Darcy’s law, the laminar velocity is very small, so the velocity head can be ignored. The total head can be regarded as the sum of pressure head and position head.(3)H=uwγw+z,where uw is the pore water pressure, γw is the unit weight of water, and z is the elevation.


## 3.1.2. Establishment of the Numerical Model

The construction phase of deep foundation pits was modelled using the MIDAS/GTS software. The actual excavation part of the deep foundation pit in this project is 158 m long, 26 m wide, and 37 m deep. Considering the boundary conditions, the whole numerical model is 299 m long, 167 m wide, and 55 m deep. There are 125248 nodes and 248745 elements in the whole model. The physical and mechanical parameters of each soil layer are shown in Table1. At present, the pit is under construction and has been excavated to 8 m. Combined with the site monitoring data, the geotechnical parameters shown in Table 1 are consistent with the site monitoring data. The numerical model is shown in Figure 4. The Mohr-Coulomb constitutive model was used in the soil layer of modulus value type. Crown beam, lumbar beam, and steel internal support all used the beam element. This is because the displacement of the uplift pile and pit column in the vertical direction is restrained, and their force characteristics are similar to those of the beam. Therefore, they also use the beam element. The foundation pit floor and underground diaphragm wall were simulated with the plate element. When establishing the waterproof curtain on the underground diaphragm wall because the waterproof curtain belongs to the interface element which is equivalent to cutting the nodes, it will cause the nodes to be uncoupled. Therefore, before the underground diaphragm wall is activated with the waterproof curtain, a rigid connection should be established to bind the cut nodes together to make the model coupled as a whole body. Prestressing anchors were used using the implanted truss element.Table 1
Physical and mechanical parameters of surrounded layers.
Soil layerThickness (m)Elastic modulus (MPa)Poisson’s ratioμγ (kN·m−3)Permeability coefficient (m·d)Cohesion (kPa)Internal friction angleφ (°)Plain fill50.30.1718301020Silty clay4500.33200.0234.916.5Coarse sand7180.320.5150.130Fractured rock10120000.2250.5500050Granite30300000.1260.0013000060Figure 4
MIDAS/GTS numerical simulation diagram.(1) Boundary Conditions of Model Displacement. The normal displacement around the model and surface is constrained, and the top surface is the free surface. The seepage boundary condition is the impermeable boundary at the place where the continuous underground diaphragm wall is precipitated around the foundation pit. The permeability coefficient is 0 m/d at the place where the waterproof curtain is set. However, in the excavation of the pit to 15 m below, the excavation was carried out deep into the rock layer. Because the permeability coefficient is very small, the actual project did not set the waterproof curtain. However, the project section has the possibility of infiltration of groundwater outside the pit. For example, under about 10 m underground, a section of sandy fractured rock with a large permeability coefficient invaded the granite interior. This section is very likely to allow underground water from outside the pit flow into the pit.


## 3.1.3. Deep Foundation Pit Construction Phase

The simulation of the deep foundation pit construction phase is shown in Figure5. The excavation of this deep foundation pit is divided into six steps. The precipitation is divided into six steps. Five supports were applied. Prestressed anchors were applied after the excavation of the 4th step of precipitation, and finally, the uplift piles and the foundation slab were applied.Figure 5
Deep foundation pit construction phase.


## 3.2. Seepage Field Characteristics

The seepage path and corresponding pore water pressure of the deep foundation pit after six precipitation events were extracted, as shown in Figure6. The seepage path of the deep foundation pit after each precipitation is different. Intrusive broken diorite porphyrite and medium coarse sand are the key points and belts for waterproofing the deep foundation pit. With the precipitation-excavation of the deep foundation pit, the pore water pressure at the bottom of the foundation pit shows a distribution of three “concave” shapes. A high-permeability pressure zone occurs around the foundation pit, the intrusive broken diorite porphyrite zone, and the middle coarse sand zone. With further excavation of the foundation pit, the seepage pressure in the middle of the foundation pit gradually decreases, and the two “concave” shapes in the middle of the foundation pit gradually merge together. After excavation to the bottom of the foundation pit, the pore water pressure at the bottom of the foundation pit shows two asymmetrical “concave” shapes, and the peak value is located at the place where broken diorite porphyrite is intruded, which is prone to water inrush. The four corners of the foundation pit are prone to seepage accumulation and subduction corrosion and finally form the piping zone.Figure 6
Seepage path and pore water pressure evolution of the deep foundation pit after precipitation.
(a)
Seepage path and pore water pressure distribution at the bottom of the pit of the 1st precipitation(b)
Seepage path and pore water pressure distribution at the bottom of the pit of the 2nd precipitation(c)
Seepage path and pore water pressure distribution at the bottom of the pit of the 3rd precipitation(d)
Seepage path and pore water pressure distribution at the bottom of the pit of the 4th precipitation(e)
Seepage path and pore water pressure distribution at the bottom of the pit of the 5th precipitation(f)
Seepage path and pore water pressure distribution at the bottom of the pit of the 6th precipitationFigure6(a) shows the seepage path of the 1st precipitation and the distribution of pore water pressure at the bottom of the pit. Three seepage gathering points are present in the deep foundation pit. Two of them are at the corners of the pit on both sides of the pit, and the stratum above them is medium-coarse sand. The other one is intrusive broken diorite porphyrite area. Although there are three seepage accumulation points, the height of the foundation pit seepage path at this time only reaches the silty clay layer, and the pore water pressure of the foundation pit floor shows negative pressure.Figure6(b) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. The height of the seepage path of the deep pit only reaches the silted clay layer, and the pore water pressure at the bottom of the pit slightly varies with the positive pore water pressure at the bottom of the pit.Figure6(c) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. When the deep foundation pit is excavated to about 15 m underground, i.e., when the diaphragm wall and steel pipe pile are converted, the maximum pore water pressure at the bottom of the pit is 135 kN/m2. It is highly recommended to make waterproof work before continuing the construction.Figure6(d) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. The seepage path of the deep foundation pit became more and more concentrated, and the pore water pressure is distributed in 3 “concave” shapes at the bottom of excavation. High-permeability pressure zones are present around the foundation pit, the intrusive broken diorite porphyrite zone, and the middle coarse sand zone.Figure6(e) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. Although the pore water pressure is still distributed in three “concave” shapes at the bottom of the pit at the excavation, the two “concave” shapes between the intrusive broken diorite porphyrite and the terminal mile of pit tend to merge together.Figure6(f) shows the seepage path of the 2nd precipitation and the distribution of pore water pressure at the bottom of the pit. After the end of precipitation, the pore water pressure at the bottom of the pit showed two very obvious “concave” shapes, mainly due to the presence of intrusive broken diorite porphyrite in the deep pit, which changed the groundwater flow path of pit excavation and showed multiple “concave” patterns.


## 3.3. Deformation Characteristics

Figure7 shows the deformation evolution of deep foundation pit precipitation. The surface settlement caused by excavation is the largest along the longitudinal axis of the deep foundation pit in the area invaded by the broken zone, and the largest deformation along the horizontal axis of the foundation pit was found near the foundation pit side. After the excavation of the foundation pit entered the moderately weathered rock stratum, the surface settlement became stable. With the excavation of the deep foundation pit, the convergence deformation of the underground diaphragm wall to the pit occurred, with a maximum value of about 25 mm. The change in the underground diaphragm wall was only above the construction of the medium weathering rock layer. With the construction of steel pipe piles, the deformation of the underground diaphragm wall became stable. After the first precipitation excavation of the deep foundation pit to the powder clay, the bottom of the pit had an uplift to the maximum. With further precipitation-excavation, the uplift of the bottom had only a small change.Figure 7
Deformation evolution of deep foundation pit precipitation excavation.
(a)
Deformation characteristics of the 1st precipitation excavation(b)
Deformation characteristics of the 1st precipitation excavation(c)
Deformation characteristics of the 1st precipitation excavation(d)
Deformation characteristics of the 1st precipitation excavation(e)
Deformation characteristics of the 1st precipitation excavation(f)
Deformation characteristics of the 1st precipitation excavation


## 4. Conclusion

Based on a project of a deep foundation pit in a subway station, the characteristics of groundwater flow path evolution and deformation field during the whole process of deep foundation pit precipitation excavation with water-rich intrusive broken diorite porphyrite were simulated using a numerical simulation method, and the conclusions are as follows:(1)
Intrusive broken diorite porphyrite and medium coarse sand are the key points and belts for waterproofing a deep foundation pit. High-permeability pressure zones are present around the foundation pit, the intrusive broken diorite porphyrite zone, and the middle coarse sand zone. Water inrush easily occurs in the place of broken diorite porphyrite. The four corners of the foundation pit are prone to seepage accumulation and subduction corrosion and finally form the piping zone(2)
The pore water pressure at the bottom of the deep foundation pit gradually shows a distribution of three “concave” shapes with the excavation of precipitation. With further excavation, the seepage pressure in the middle part of the pit gradually decreases, and the two “concave” shapes in the middle gradually merge together. After excavation to the bottom of the pit, the pore water pressure in the bottom of the pit is asymmetrically distributed in two “concave” shapes(3)
The surface settlement caused by deep foundation pit excavation is the largest along the longitudinal axis of the deep foundation pit in the intruding and crushing zone. In the horizontal axis direction of the foundation pit, it is the largest near the side of the foundation pit, and the surface settlement at one time of the excavation depth is approximately 0. The surface settlement becomes stable after the foundation pit excavation reaches the moderately weathered rock stratum(4)
With the excavation of the deep foundation pit, the convergence deformation of the diaphragm wall towards the interior of the foundation pit occurs, and the maximum value of convergence deformation is about 25 mm. The change in the underground diaphragm wall was only above the construction of medium weathering rock layer. With the construction of steel pipe piles, the convergence deformation of the underground diaphragm wall becomes stable(5)
After the first precipitation excavation of the deep foundation pit to silty clay, the uplift at the bottom of the pit is the largest, but with further precipitation-excavation, the uplift only changes slightly



---

*Source: 6628882-2021-02-13.xml*