Dewatering Characteristics and Drawdown Prediction of Suspended Waterproof Curtain Foundation Pit in Soft Soil Areas

: Suspended waterproof curtains combined with pumping wells are the primary method for controlling groundwater levels in foundation pits within soft soil areas. However, there is still a lack of a systematic approach to predict the groundwater drawdown within the foundation pit caused by the influence of these suspended curtains. In order to investigate the variation of ground-water level within the excavation during dewatering processes, the finite difference method is employed to analyze the seepage characteristics of foundation pits with suspended waterproof cur-tains. Basing on the concept of equivalent well, this study examines the coupled effects of aquifer anisotropy (k i ), aquifer thickness (M i ), well screen length (l i ), and the depth of waterproof curtain embedment on the seepage field distortion. A characteristic curve is established for standard conditions, which exposes the blocking effect of the curtain on the amount of groundwater drawdown in the pit. Additionally, correction coefficients are proposed for non-standard conditions, which, in turn, results in a prediction formula with a wider range of applicability. Comparative analysis be-tween the calculated predictions and the field observation data from an actual foundation pit project in Zhuhai City validates the feasibility of the quantitative prediction method proposed in this research, which also provides a 21% safety margin.


Introduction
As urbanization accelerates in southern coastal regions of China, the population density continues to rise, leading to an increasing demand for the construction of large buildings along the coastal areas [1].Foundation pits, as primary underground structures for constructing large buildings, have seen a continuous increase in both quantity and depth, along with an expanding excavation area.In coastal regions of the Pearl River Delta, China (e.g., Zhuhai, Zhongshan, and Jiangmen), with abundant groundwater in aquifers [2], the risks associated with dewatering during construction have correspondingly grown [3][4][5].To ensure a safe construction environment within the foundation pit, the application of dewatering measures in excavation projects has also extended to a broader scope [6].In contrast to inland areas, the soil in coastal regions is primarily composed of alluvial convergent sediments by rivers and marine soft sediments alternately, mainly consisting of silt sand and low-permeability clay.These soil layers are characterized by high compressibility and low permeability [7,8].Due to the large thickness of the submerged layer, the majority of groundwater aquifers are categorized as phreatic aquifers rather than confined aquifers.In addition, the initial groundwater tables are generally higher in coastal areas, and most of the target depths of excavation for many building foundation pit projects are lower than the location of the initial groundwater tables [9], which is more significantly affected by the groundwater level.
The excavation process in these areas does not isolate the hydraulic connection between the inside and outside of the foundation pit, resulting in a difference in the hydraulic head, which allows groundwater to flow into the pit [10].To prevent issues such as sand inflow, sudden floor heave, and water inflowing towards the excavation [11,12], it is crucial to implement necessary dewatering measures.Additionally, improper dewatering can lead to rapid ground settlement, potentially causing damage to nearby structures, especially in densely populated urban areas [13,14].In preventing major financial and material losses resulting from dewatering accidents, it is of utmost importance to prudently control the groundwater level within the foundation pit, ensuring the safety and stability of the construction site, which holds significant practical significance.Dewatering of excavation aims to minimize surrounding ground settlement and pumping volume while ensuring a safe and dry construction environment and avoiding the waste of water resources.The common approach to controlling groundwater pressure involves the use of low-permeability materials such as cement mixing piles, jet-grouted piles, and diaphragm walls [15][16][17] as waterproof curtains, combined with pumping wells for dewatering.Curtain as a waterproof structure of the pit project can effectively reduce the hydraulic connection between the inside and outside of the pit, block the seepage of groundwater in the direction of the pit [18][19][20], and also play the role of supporting structure.Waterproof curtains can be categorized as penetrating curtains and suspended waterproof curtains based on whether they fully penetrate aquifers.In theory, penetrating waterproof curtains can entirely block hydraulic connections between the inside and outside of the excavation.However, they often encounter challenges when attempting to penetrate aquifers, leading to potential defects [21].Thick aquifer layers not only limit the practicality of fully penetrating curtains but also make the exposed curtain above the pit bottom vulnerable to localized cracking and leakage due to the external soil pressure after excavation, posing a risk to construction safety [22].Consequently, suspended waterproof curtains are becoming increasingly popular in foundation pit dewatering projects when compared to full penetrating waterproof curtains.
To investigate the blocking effect induced by suspended waterproof curtains [23], numerical simulation methods such as MODFLOW finite difference software [24][25][26] have been widely adopted by many researchers.Studies have shown that the hydraulic isolation provided by the curtain is manifested in three key aspects: extending the groundwater seepage path, altering the seepage direction, and reducing the seepage area [27].Furthermore, the permeability anisotropy and the thickness of the aquifer [28,29] also have certain impacts on groundwater seepage.Presently, the construction methods for waterproof curtains have become increasingly mature, while the coupled blocking effect of pumping wells and curtains has been studied at a deeper level [30,31].A large number of experiments have been used to verify the dewatering effect of wall-well coupling.The application of numerical simulation methods aids in determining the optimal embedment depth of the curtain [32] and evaluating the feasibility of dewatering methods [33], which provides a reference for the study of the seepage mechanisms and field application of suspended waterproof curtains.Research into the inflow prediction and head differences on both sides of the waterproof curtain for confined aquifer excavation has been extensive [34,35].However, methods for forecasting drawdown in the foundation pit after the operation of dewatering wells often rely heavily on numerical simulations.Such research approaches are hindered by their high modeling complexity, intensive computational requirements, and limited range of applicability.Many of these studies primarily consider the embedding depth of waterproof curtains as a single variable and explore its impact on groundwater seepage without considering the coupled effects of various influencing factors.To make up for the lack of research on groundwater seepage characteristics during the dewatering process of foundation pits in phreatic aquifers, this study aims to quantitatively assess the blocking effects of curtains, taking into account the anisotropy of permeability, aquifer thickness, and the length of the well screen in dewatering wells, to derive a more widely applicable formula for predicting drawdown in unconfined aquifer foundation pits and provide a reference method for the prediction of the effectiveness of the design and construction monitoring of similar foundation pits dewatering in the soft soil areas.To address the problem, this paper first analyzes the seepage mechanisms of suspended waterproof curtains in unconfined aquifers.Subsequently, numerical simulation software MODFLOW is employed to solve standard and variational functions that influence groundwater seepage.Finally, a comparative analysis between the numerical solutions and field observations is conducted to validate the applicability of the proposed equations.

Seepage Analysis
When conducting dewatering operations within an unconfined aquifer foundation pit, groundwater from the surrounding aquifer flows into the dewatering well, resulting in a lowering of groundwater levels in the central region of the well and its immediate vicinity.The influence of the dewatering well on drawdown in the surrounding aquifer becomes less pronounced as the distance from the well center increases.The maximum drawdown occurs at the well center and decreases as the well diameter moves outward from the center of the well.The overall drawdown profile exhibits a continuous funnellike shape, as illustrated in Figure 1a.Since dewatering wells do not fully penetrate the aquifer to the aquitard, groundwater near the dewatering well bottom experiences a complex three-dimensional (3D) seepage field.To simplify the description of groundwater seepage, horizontal and vertical seepage are used as substitutes for 3D seepage.After the installation of the waterproof curtain, the seepage environment around the foundation pit site is influenced, resulting in changes to the groundwater flow field near the curtain.The hydraulic containment effects of the curtain are primarily manifested in three aspects [36,37]: 1. Alteration of seepage direction and elongation of the seepage path: In the absence of the waterproof curtain, groundwater flows horizontally into the well along the radial direction of the pumping well, with vertical seepage occurring around the pumping well, as shown in Figure 1a.However, with the installation of the waterproof curtain during dewatering, there is a change in the direction of groundwater flow outside the curtain at the base of the curtain, as depicted in Figure 1b.Groundwater outside the foundation pit site needs to circumvent the bottom of the waterproof curtain to reach the vicinity of the pumping well, resulting in an extended flow path compared to when no curtains are in place.2. Modification of groundwater drawdown within and outside the foundation pit: In the absence of the curtain, the drawdown profile exhibits a continuous funnel shape.However, when the waterproof curtain is installed, water levels inside and outside the excavation show a head difference on either side of the curtain, leading to a discontinuous, step-like drawdown profile on both sides of the curtain.3. Reduction in the total inflow to the foundation pit: According to Darcy's law [38], when pumping is conducted at a constant flow rate, the elongation of the groundwater flow path due to the presence of the waterproof curtain leads to a reduced seepage area.This decrease in seepage area results in an increased hydraulic gradient, reducing the required dewatering time and improving dewatering efficiency.Consequently, this leads to a reduction in the overall seepage volume within the foundation pit.
All of the above descriptions of groundwater seepage movement patterns in phreatic aquifers are based on the following assumptions: 1.The groundwater flow towards the well is assumed to be approximately horizontal. 2. The aquifer is considered to be homogeneous, isotropic, of uniform thickness, and laterally infinite.3. Before pumping, the natural hydraulic gradient in the aquifer is assumed to be zero.4. Pumping is assumed to be of long duration, and groundwater flow within the aquifer is expected to follow Darcy's law, reaching a state of steady flow.

Drawdown in Equivalent Well
As illustrated in Figure 2a, based on the Dupuit theory of drawdown curves [39], the foundation pit with the suspended waterproof curtain can be analogously represented as a two-dimensional partial pumping well of radius r0 (the radius of the equivalent pumping well for the excavation) in the phreatic aquifer.The calculation of the drawdown in the equivalent well using the Theis formula [40] can be expressed as follows: In the equation, the variables are defined as follows: st represents the drawdown in water level at the location r0 (m); Qw stands for the flow rate of the equivalent dewatering well associated with the excavation (m 3 /d); T is transmissivity of the aquifer (m 2 /d), where T = KH (K: hydraulic conductivity; H: aquifer thickness); S denotes the specific storage coefficient of the aquifer indicates the characteristics of the water outflow when the water level within the aquifer decreases; s corresponds to the observed drawdown in water level (m); r0 is the radius of the equivalent pumping well for the excavation, calculated as r0 =   , where A represents the area of the excavation(m 2 ); and H signifies the thickness of the phreatic aquifer (m).According to the Theis equation, it is evident that the transmissivity (T) and the thickness of the phreatic aquifer (H) directly influence the drawdown in water level after dewatering of the equivalent well associated with the foundation pit.However, this calculation method does not take into account the indirect effects caused by the distortion of the seepage field.The blocking effect of the waterproof curtain induces a complex threedimensional flow rather than a parallel flow, leading to changes in seepage direction, paths, and seepage area, which subsequently deform the seepage field.Hence, the blocking effect of the waterproof curtain cannot be overlooked and requires the introduction of quantifying coefficients.

Numerical Simulation
For the equivalent circular foundation pit described above, a numerical model is established based on typical geological properties in the soft soil areas of the Pearl River Delta, China [41].The foundation pit radius, denoted as r0, is set at 53 m, with a waterproof curtain thickness of 1 m.To stabilize the seepage flow in the model, based on the empirical dewatering scheme, the total dewatering time is set to be 100 days.The groundwater flow in the numerical model is approaching a steady state in this state.

Governing Equation
In order to obtain a conservative estimate for the drawdown design values, the assumptions of Darcy law and continuity are employed, neglecting no-Darcy flow.Considering the effect of unsteady flow, modeled groundwater 3D unsteady flow in the study area is controlled by the following equation.The formulation includes the establishment of a groundwater flow governing equation, initial groundwater head, and boundary conditions within a saturated porous medium [38] In the equation, Kxx, Kyy, and Kzz represent the hydraulic conductivities in the x, y, and z directions, respectively (cm/s); h is the groundwater level at point (x, y, z) (m); Q is the groundwater recharge and discharge (d −1 ); Ss represents the specific storage at point (x, y, z) (m −1 ); T is the elapsed time (h); h0 is the initial groundwater level at point (x, y, z) (m); Γ1 and Γ2 correspond to the first and second class boundary conditions, respectively; nx, ny, and nz are the unit normal vectors along the x, y, and z directions for boundary Γ2; and q is the lateral inflow per unit area on boundary Γ2 (m 3 /d).In this simulation, it is assumed that the soil is isotropic in the horizontal direction, meaning Kxx equals Kyy.

Simplification of Model Parameters
Table 1 outlines the parameters of the simplified calculation model.A phreatic aquifer (which allows the movement of water through it under ordinary conditions) and an aquitard (which has poor permeability and hence will not yield water freely to wells) are considered at 60 m depth.The soil types of the phreatic aquifer are backfill and sandy clay.The soil type of aquitard is silty clay.In the numerical simulation calculations, it is assumed that the soil is isotropic.To investigate the groundwater drawdown under the coupled influence of waterproof curtain and pumping wells, the analysis is conducted under standard conditions, with the waterproof curtain embedding depth (Hu) set at 36 m.The length of the well screen (l) is equal to Hu, the initial groundwater table in the phreatic aquifer is 2 m, and the aquifer thickness (H) is selected as 38 m, as depicted in Figure 3.

Initial Boundary Conditions
To mitigate boundary effects, it is essential to determine appropriate model dimensions.Siechardt's [42] empirical formula is employed to calculate the influence radius of the equivalent well of the phreatic aquifer foundation pit: In the equation, R represents the influence radius of the equivalent well (m); H is the aquifer thickness before pumping (m); K is the aquifer permeability coefficient; and s denotes the drawdown in the equivalent well.
In accordance with the requirements, considering an aquifer thickness of 38 m, the calculated drawdown influence radius is 259 m.Consequently, as illustrated in Figure 4a, the equivalent circular well radius (r0) of the model is set to 54 m, and the finite difference model extended 260 m in all directions.This resulted in the creation of a three-dimensional grid with dimensions measuring 700 m in length and width, and 60 m extending in the zaxis direction.The model is divided into 42 layers, 4 for aquitard and 38 for phreatic aquifer, with 88 columns of values in 96 rows horizontally.The grid is fine in the excavation area of the pit.The grid is coarse at the edge of the excavation.
The boundary conditions of the model are defined with constant head to simulate hydraulic recharge far from the source-sink terms [43], ensuring that the boundary remained at the original head level, unaffected by the drawdown within the foundation pit, as shown in Figure 4b.Referring to hydrogeological survey data, the initial groundwater level for the project is set at 2 m.Within the conceptual model of the excavation established by simulation, the pumping well and cutoff curtain units are projected onto the computational grid.At the well screen, a constant flow boundary condition is established to simulate drainage, while the waterproof curtain is represented by low-permeability units of 1 × 10 −5 .To ensure conservative results, the model employed the parameters listed in Table 1. Figure 4c illustrates the positions of the pumping well and the waterproof curtain configuration.

Drawdown Coefficient
Adjust the pumping volume (Qw) of the equivalent well in the three-dimensional numerical model to solve the groundwater drawdown sc after pumping by the action of the waterproof curtain.The Theis theoretical solution (st) is linearly related to Qw.The coefficient si of the influence of three-dimensional seepage caused by different depths of embedment of the waterproof curtain on the depth of the groundwater drawdown within a certain pumping flow rate can be expressed as follows:

Normalized Form
To ensure non-dimensional computational conditions, dimensionless form parameters are introduced, including the waterproof curtain embedment ratio (bi = H/Hu), the permeability anisotropy (ki = kh/kv), the thickness of phreatic aquifer (Hi = H/38), and the well screen length (li = l/Hu), serving as influencing parameters.For this study, the parameter values are set as ki = 1, Hi = 1, and li = 1, representing the normalized situation for investigating the influence of the waterproof curtain on drawdown within the foundation pit.Under these standard conditions, different waterproof curtain embedment ratios, the Theis analytical solution curves (Qw-st), and the numerical model results (Qw-sc) are presented in Figure 5a. Figure 5 illustrates the impact of the cutoff barrier configuration on the distribution of the seepage field both within and outside the foundation pit.The presence of the cutoff barrier leads to an extension of groundwater flow paths outside the excavation, reducing the seepage area and causing a redistribution of hydraulic heads around the curtain.The simulation results indicate that horizontal groundwater flow predominantly occurs in the lower aquifer region beneath the waterproof curtain.The presence of the waterproof curtain necessitates the dissipation of a greater amount of vertical hydraulic potential energy in groundwater flow, leading to a weakening of the hydraulic connection between the interior and exterior of the foundation pit.
Under the influence of the waterproof curtain, groundwater seepage follows Darcy's law.As observed in Figure 5a, the Qw-s curves exhibit linear distribution for different depths of cutoff wall insertion, indicating that with a constant pumping volume, greater waterproof curtain insertion depths result in a more significant drawdown within the foundation pit, signifying a more pronounced blocking effect on groundwater seepage.However, different waterproof curtain insertion ratios (bi) lead to variations in the slopes of the Qw-s curves.The impact of waterproof curtain insertion depth on the drawdown is quantified by the ratio of the slopes of Qw-sc and Qw-st, representing the influence coefficient (si).The curve in Figure 5b illustrates that, under standard conditions, increasing the waterproof curtain insertion ratio intensifies the drawdown within the foundation pit.In other words, deeper cutoff walls enhance the blocking effect on groundwater seepage.As the depth of cutoff wall insertion increases, the slope of the Qw-s curve continually rises.The bi-si curve demonstrates that the blocking effect of the waterproof curtain grows linearly and then logarithmically with increasing bi.When bi approaches 1, it can be considered that the cutoff walls entirely cut hydraulic connections between the interior and exterior of the foundation pit, resulting in the maximum si value.
In addition, there is a linear and nonlinear turning point in the fitted bi-si curve, which is usually regarded as the blocking ratio in the case of an abrupt change in the groundwater drawdown discipline in the available numerical solutions and experimental studies of the model.In this study, bi = 0.21 is selected as the break point and fits the bi-si curve for the standard aquifer conditions.The fitting results are presented in Equation ( 5) (with correlation coefficients  and  of 0.958 and 0.999 for linear and exponential functions, respectively).

Analysis of Factors Affecting Seepage
The numerical model established above is based on an assumption of isotropic properties in the aquifer.However, in practical engineering, aquifers are often anisotropic, and the influence of this anisotropy on seepage should be considered.Additionally, different projects are located in different geographical areas, resulting in variations in geological conditions and aquifer parameters.Therefore, the aquifer thickness should also be considered as a factor affecting the effectiveness of waterproof curtains.In addition, the actual filter pipe lengths differ from the numerically modeled filter pipe lengths in the standard conditions, depending on the precipitation requirements of a specific project.To make the bi-si curve more applicable to actual geological and construction situations, numerical models are established under different scenarios to simulate and adjust Equation ( 5).The fitted curve under standard conditions from Figure 5b is chosen as the baseline, and the numerical simulations were used to quantify the impact of permeability anisotropy (ki), aquifer thickness (Mi), and well screen length (li) on the blocking effect of the waterproof curtain.

Anisotropy of Permeability
The hydraulic permeability of soil exhibits anisotropy due to factors such as depositional history, initial non-uniform consolidation, and variations in stress, resulting in varying characteristics of groundwater flow within the aquifer [44].The anisotropy ratio of permeability in horizontal and vertical directions, denoted as ki (ki = kh/kv), defines the anisotropy of the aquifer.The vertical flow resistance within the soil becomes more prominent as the anisotropy of hydraulic conductivity increases, and the installation of the curtain exacerbates this effect by inducing vertical flow.The curtain serves to weaken the hydraulic connection between the interior and exterior of the foundation pit, enhancing dewatering efficiency.As the anisotropy of soil permeability increases, the vertical flow velocity decreases, and the blocking effect of the curtain results in the consumption of a greater amount of gravitational potential energy.Keeping all other coefficients constant, only adjusting the ki parameter, the impact of non-standard permeability anisotropy on the bi-si curve is computed, and the simulation results are presented in Figure 6. Figure 6 indicates that as the ki value increases, the bi-si curve does not coincide with the standard curve, providing evidence that ki does have an impact on the blocking effect of the curtain.The value of ki varies from 3 to 9, and with increasing values of ki, the curve deviation becomes more pronounced.The ratio of the numerical solution for the drawdown to the analytical solution increases, suggesting that the permeability anisotropy enhances the blocking effect of the waterproof curtain.To extend Equation ( 5) to different forms of permeability anisotropy, a correction coefficient function denoted as α is introduced.Its expression is given in Equation ( 6):

Thickness of the Aquifer
While maintaining the ratio of the insertion depth of the waterproof curtain (bi) constant, a reduction in the thickness of the aquifer (H) leads to lower piezometric head values in the surrounding area.This results in a decrease in the seepage field area and the vertical hydraulic recharge, ultimately reducing the required dewatering discharge to control drawdown.Changes in the thickness of the aquifer also affect the seepage field distribution, leading to a nonlinear relationship between drawdown and aquifer thickness.Introducing the non-standard aquifer thickness coefficient, by adjusting Hi (Hi = H/38), the impact of non-standard aquifer thickness can be calculated on the bi-si curve, as illustrated in Figure 7.  Figure 7 illustrates that the blocking effect curve exhibits varying vertical displacements with changes in Hi, and the linear patterns remain similar across different Hi values.The thickness of the aquifer significantly influences the blocking effect of the curtain.When Hi > 1, which indicates an increase in aquifer thickness, the curve tends to shift vertically upwards, indicating enhanced curtain barrier effectiveness.Conversely, when Hi < 1, the barrier effect weakens.Similar to the method used for solving permeability anisotropy coefficient, a non-standard aquifer thickness coefficient is determined and represented as β:

Pumping Well Screen Length
The pumping capacity of well screens is continuous, and a smaller proportion of well screen length results in longer drainage paths [45].Under constant discharge conditions, shorter pumping well screen lengths have a more pronounced impact on drawdown levels, leading to greater drawdown depths within the excavation.Therefore, while meeting the fundamental dewatering requirements, it is advisable to select shorter well screen lengths to enhance dewatering efficiency.By adjusting the non-standard well screen length coefficient, li (li = l/Hu), the impact of non-standard well screen lengths on the bi-si curve can be computed, as illustrated in Figure 8.  Figure 8 demonstrates that the blocking effect of the curtain varies with different values of li.While the shape of the bi-si curve remains consistent across different li values, a reduction in the length of the well screen leads to a vertical upward shift in the curve.This shift provides evidence that the well screen length enhances the blocking effect.Similar to the method used for solving aquifer thickness coefficient, a coefficient of variation, γ, for non-standard well screen lengths is determined, and its expression is as follows:

Drawdown Prediction with Suspended Waterproof Curtain
Based on the equivalent well assumption, modifications to the standard curve can be made using variation functions under non-standard conditions.Considering the impact of a single variable on the barrier effect of the curtain, the normalized bi-si curve can be adjusted to a reasonable range.Simultaneously, there are coupling effects among the parameters, leading to a nonlinear increase in si.By extending the range of applicability of the bi-si curve using the variation functions derived above, predictions of dewatering drawdown are brought into alignment with actual observations.This yields non-standard characteristic curves for the phreatic aquifer and curtain impact coefficient, denoted as sni, as expressed below: Combined with the Theis formula, the prediction formula for the dewatering drawdown depth sd of the phreatic aquifer foundation pit under the action of the waterproof curtain in a soft soil area is derived, and its expression is

Project Description
In order to verify the applicability of the calculation method of drawdown depth in the pit under the action of the suspended waterproof curtain proposed above, a largescale foundation pit required for the construction of a building in Zhuhai City is used as a field case for verification.The layout of this large pit is shown in Figure 9a.The project site is located in a land area with flat terrain and relatively open surroundings.Within 50 m on the north side of the site is a residential community with some elevated buildings, surrounded by grassland and accessed by major urban roads.The absolute elevation of the outdoor ground is 5.35 m.
The pit as a whole is an irregular rectangle which oriented northeast-southwest, with a support perimeter of 368 m and an area of about 8920.6 m 2 .As shown in Figure 9b

Hydrogeological Conditions of the Area
The soil profile and soil properties are shown in Figure 10.The soil layers, from top to bottom, are backfill (11), silty clay (12), coarse sand (13), sandy clay (14), fine, chalky sand (15), chalky clay (21), fully weathered granite (31), sandy, gravelly soil (41), and strongly weathered granite (51).The project includes a phreatic aquifer (Aq01), two confined aquifers (Aq I and Aq II), and an aquitard (Ad I and Ad II).Detailed information on the aquifer is as follows: During the investigation, groundwater is encountered in all boreholes.As the thickness of the pressurized water in the coarse sand layer ( 13) is thin and its distribution is discontinuous, all of which have an insignificant impact on the construction of the project, it has been determined that the groundwater is submerged at the proposed site.The initial groundwater level of Aq01 is at a depth of 2.33-3.92m and the source of groundwater is mainly from the upper phreatic aquifer.It is closely connected to atmospheric precipitation and surface water.As a result of the climate and seasons, the water level of the groundwater fluctuates, and when water is abundant, the water level rises.The confined aquifers AqⅠ and AqⅡ are rock fissures filled with pressurized water.It exhibits the following characteristics: it is slightly pressurized, and the distribution of groundwater is greatly affected by the development of fissures in the endowed rock body, with obvious anisotropic characteristics, and it is a non-homogeneous seepage field.In the section where the nodal fissures are more developed, fissure water is more abundant.
The distribution of pumping wells in this project is divided into two parts with 1-1 section, numbered DW1-n(n=1-15) and DW2-n(n=1-12).The length of the pumping wells is 23.0-24.8m, the effective length of the filter pipe is 5-12 m, and the spacing of the precipitation wells is about 14.8-17.5 m.At the same time, the observation wells are set up around the pumping wells.The location of the pumping wells and the observation wells is shown in Figure 11.

Results
In order to ensure the effect of dewatering, the aquifer needs to be generalized at the design stage.Prior to curtain construction, indoor core tests and field pumping tests are carried out to determine the horizontal and vertical permeability coefficients of the aquifer.These coefficients are 6.322 m/d and 4.33 m/d, respectively.
Firstly, from Equations ( 1) and ( 4), the equivalent well radiuses of deep and shallow pits are 41.69 m and 34.3 m, respectively.The total pumping volume of deep and shallow pits are 2400 m 3 /d and 2000 m 3 /d, respectively.The depth of pit reduction under the action of pumping wells without waterproof curtains of deep and shallow pits are 3.77 m and 3.76 m.When the curtain is not taken into account, the depth of groundwater reduction generated by the pumping wells can obviously not meet the demand for pit dewatering.The origin method is unable to accurately predict the actual drawdown of groundwater level in the pit.Blindly increasing the pumping volume will cause unnecessary wastage of water resources.The prediction method described in this paper is based on the simulation of an equivalent circular pit model, which proves that the change in the water level in the pit by the precipitation project is independent of the converging shape of the pit according to the results of the study by Lyu et al. [46].Therefore, according to Equations ( 5)-( 9), the actual parameters of the project are brought in to calculate the curtain-blocking effect coefficients and the value of the variation function.The calculation results are shown in Table 2. Combined with the field observing data, the predicted calculation method case proposed in this paper is verified.The data show that when the total pumping volume of the deep pit on the southwest side and the shallow pit on the northeast side reaches 2400 m 3 /d and 2000 m 3 /d, respectively, the predicted values of groundwater drawdown in the deep and shallow pits in 30 days are 9.2 m and 7.7 m after the correction coefficients of Table 2 are taken into account, and all of them satisfy the project's safe groundwater drawdown requirements [47].
The pumping capacity of the pumping wells DW1-n(n=1,2,3,5,6,7) on the southwest side and DW2-n(n=1,2,3,6,7,8) on the northeast side of the pit are used as reference, with the pumping capacities of the reference pumping wells on the southwest side ranging from 150 to 165 m 3 /d, and the pumping capacities of the pumping wells on the northeast side ranging from 160 to 168 m 3 /d.At the same time, based on the observation data from the observation wells (OW1-1, OW1-2, OW2-1, and OW2-2), the groundwater level changes in the pit are reflected.When the fluctuation of water level in the observation wells is less than 2-3 m, it is considered to reach a stable groundwater level [48].The observation process is accumulated for 30 days.Figure 12 shows the groundwater level observation data in the observation wells.As the working time of the pumping well grows, the groundwater level in the pit decreases in depth.The change in drawdown decreases is first steep and then becomes flat.The groundwater level between the two observation wells in the deep and shallow pits is basically the same, indicating that the water pressure in the pits is uniformly distributed in the process of precipitation, and the maximum depth of the groundwater level in the deep and shallow pits after the final stabilization is −11.2 m and −8.9 m.Accordingly, the study concludes that the predicted value of groundwater drawdown in the pits has a safe redundancy with a value of between 16 and 21%.
Since the waterproof curtain did not completely isolate the hydraulic connection between the inside and outside of the pit, the operation of the pumping wells in the pit has a certain impact on the groundwater outside the pit.The groundwater outside the pit fluctuates, but the head change value is between 1.0 and 1.5 m, which meets the demand for groundwater control outside the pit.As the project is in a coastal area with a strong recharge capacity of the aquifer, the head outside the pit can recover to the initial state similar to the pre-construction state.The foundation pit surface remains stable throughout the excavation and dewatering process, which indicates that the calculation method for predicting the amount of groundwater drawdown in the pit proposed in this study is feasible.

Conclusions
Based on the theory of three-dimensional seepage analysis of phreatic aquifer, a formula is proposed to predict the groundwater drawdown in the foundation pit under the action of the suspended waterproof curtain and pumping wells.The formula is applied to the field foundation pit dewatering project.The comparison of the observed and calculated values shows that the prediction and calculation method is reasonable and the algorithm is simple.In the actual project, the depth of curtain embedding and pumping volume can be determined to determine the depth of groundwater drawdown in the pit.Thus, this calculation method provides a reference for predicting the depth of groundwater drawdown in the pit.Based on the results of this study, the following conclusions are drawn: 1.The waterproof curtain can improve the efficiency of pit dewatering by changing the seepage field distribution and thus reducing the hydraulic connection between the inside and outside of the pit.The Thies formula does not take into account the influence of the change in the flow field morphology on the groundwater seepage, and the influence coefficient must be introduced to correct the calculation results.2. Based on the equivalent well theory, solve the characteristic curve of curtain ratio (bi) and drawdown ratio (si) under each influence coefficient.Quantify the curve to get the curtain-blocking effect coefficient under the normalized state.Considering the coupling effect generated by ki, Mi, li, and the waterproof curtain, the finite difference method is adopted, and the non-standard variability function is obtained by regression analysis of numerical results.The non-standard variability function reflects the coupled influence of the curtain barrier effect and project parameters on the amount of groundwater drawdown in the pit.3.This paper presents a calculation method for predicting the amount of groundwater drawdown in a pit for a certain amount of pumping.The established equation is applied to the Zhuhai field pit pumping design.The predicted value is close to the field wells observing results, and the original design scheme maintains a 16-21% safe margin.In addition, the change of groundwater drawdown outside the pit is less than 1.5 m, which ensures the environment inside and outside the excavation pit is safe.The prediction method is reasonable.

Figure 1 .
Figure 1.Profile view of groundwater seepage and partial pumping well in phreatic aquifer.(a) Without curtain.(b) With curtain.

Figure 2 .
Figure 2. Equivalent pit precipitation well.(a) Top view of equivalent well.(b) Profile view of equivalent well.

Figure 3 .
Figure 3. Profile view of curtain and aquifers.

Figure 5 .
Figure 5. Variation results in the standard state.(a) Qw-s curves with different bi.(b) Normalized curve.

Figure 9 .
Figure 9.The layout of foundation pit.(a) Environment around the pit.(b) Plan view of the foundation pit.

Figure 11 .
Figure 11.Location of the wells.

Table 1 .
Simplification of model parameters.

Table 2 .
Values of the variation function.