#### Abstract

Great practical significance and engineering application value can be achieved when the large amount of discharged soil produced by EPB shield tunnels is recycled and comprehensively utilized. As one of the key processes of shield construction, synchronous grouting needs a large amount of bentonite, cement, fly ash, sand, and other materials. The research on the reuse of shield muck as synchronous grouting material is carried out based on Zhengzhou subway project. The physical properties and phase of the discharged soil from EPB shield tunnels are studied by using laboratory tests and XRD. The statistics show that the shield muck meets the performance requirements of bentonite and fine sand in synchronous grouting materials. The optimal grout ratio of the reused muck is obtained based on the optimization idea of multiobjective programming by MATLAB. Considering the combined effect of seepage field, stress field, and the timeliness of the grout, the influences of grouting pressure and the filling rate of synchronous grouting on surface settlement, plastic zone of strata, and segment deformation are analyzed by using finite difference method. The results prove that the surface settlement and segment deformation can be better controlled when the grouting pressure is at 0.18 MPa and the grouting rate is at 120%–150%.

#### 1. Introduction

EPB shield method has become the most widely used construction method in urban subway construction in China. A large amount of muck is generated during shield construction. The shield muck not only occupies a large amount of land, reduces soil quality, and affects air and water quality but also has certain security risks. Besides, it is expensive to transport the muck directly to waste disposal area. Therefore, the research on the reuse of shield muck is of great practical significance and practical value.

Presently, some researches on the reuse of shield muck have been carried out, including the following aspects: (1) engineering grout treatment technology. Zha et al. [1] studied a formula for converting wastewater-based drilling fluid into environmentally friendly soil. The treatment process can achieve harmless treatment and effective utilization of drilling waste. Through data analysis of the unconfined compressive strength, triaxial compressive strength, indirect tensile strength, humidity sensitivity, and seismic modulus, Shon et al. [2] conclude that it is feasible to use modified drilling waste as a basic material for road construction. Whitaker et al. [3] studied the effects of drilling mud application rate on soil salt accumulation and leaching, and the effects of applying time and rate of the drilling mud on wheat production are discussed. Dong [4] studied the treatment technology on drilling mud and explored the solidification treatment technology on drilling waste through laboratory research. Peng et al. [5] systematically analyzed and studied several main harmless treatment methods (solidification treatment, secondary utilization of waste grout) in China. (2) Reuse of shield muck: through the treatment process of grout in Weser river tunnel excavation, Grohs [6] found that the separated sand can be used as the load-bearing layer, sound insulation wall’s cover layer, and filling material of highway subgrade according to the different particle content. Zhong et al. [7] studied the reuse of excavated silty fine sand in shield synchronous grouting. He et al. [8] used shield muck sand (modulus at 1.88) from Shuangdun Station to Wangjiadun Center Station of Wuhan Metro Line 3 as fine aggregate to prepare environmentally friendly single-liquid synchronous grouting material with high performance through optimizing experimental mix proportion. In summary, although scholars have achieved fruitful results in the reuse of shield muck, there are few studies on the reuse of EPB shield muck as synchronous grouting material, and systematic research results have not been reported. Synchronous grouting is one of the key links in shield construction. Usually, grout is injected into the gaps between shield tail and segment to fill the gap and control strata deformation [9–12]. At the same time, the grouting protection layer is formed around the tunnel to improve the impermeability of shield tunnel [13–15]. The stratum that the shield passes through is usually soft soil stratum containing lots of clay particles and fine sand. If proper treatment is carried out, the discharged soil is one of the high quality raw materials for synchronous grouting.

In view of this, based on the construction example of a section crossing silty clay stratum and clayey silt stratum with sandy silt in Zhengzhou Metro, this paper discusses the feasibility about the reuse of shield muck, and the best synchronous grouting proportion with the reused shield muck is obtained. Moreover, the influences of grouting pressure and filling rate on surface settlement, soil plastic zone, and segment deformation are studied. The research results are expected to provide reference for similar projects.

#### 2. Project Overview

The main strata of a certain section of Zhengzhou Metro are clayey silt stratum with sandy silt, and the local strata are silty clay stratum and clayey silt stratum with sandy silt. According to geological prospecting, 56% of the strata traversed by the tunnel are silty clay, 38% are fine sand, and 6% are silty clay mixed with sandy silt, as shown in Figure 1.

Two Earth pressure balance shield machines with a diameter of 6250 mm were used for construction. Reinforced concrete segments are used as tunnel lining. The segment has an outer diameter of 6.2 m, a thickness of 0.35 m, and a width of 1.5 m. The theoretical muck output per ring is 49.4 m³. The theoretical grouting volume per ring is 4.1 m³, and the filling factor is 1.5–2.0.

#### 3. Feasibility of Shield Muck as Grouting Material

##### 3.1. Engineering Features Test of the Shield Muck

The detailed composition and phase analysis of shield muck, mainly including density, water content, specific gravity, limit moisture content, void ratio, saturation, and grading feature and other physical characteristics tests, are carried out. Shield muck sample is shown in Figure 2.

**(a)**

**(b)**

###### 3.1.1. Basic Physical Properties

Basic physical properties of shield muck can be obtained by physical property measurement, and the results are listed in Table 1.

According to the measurement results, the muck sample has moderate moisture content and a viscosity index of 0.5. The soil sample is a low-liquid limit silty clay in a plastic state.

###### 3.1.2. Grading Curve of the Muck

In the analysis of particle size composition, the particle size distribution of soil above 0.075 mm was determined by sieving method, and the particle size distribution of soil below 0.075 mm was measured by densimeter method considering the existence of fine particles in the muck of silty clay. The particle size composition of the muck is manifested in Table 2. The grading curve is portrayed in Figure 3.

According to Table 2 and Figure 3, the particle size distribution of silty clay samples is uneven, and the particle gradation is good. There are a lot of fine particles in the sample. The fine particles with particle size less than 0.075 mm account for 39.32%, and the clay particles with particle size less than 0.005 mm account for 8.45%.

Only the sieving method is used to determine the particle size distribution, and the fineness modulus and mud content are used to evaluate the quality of sand samples in the analysis of particle size composition of fine sand in muck. The particle size composition of the fine sand in muck is shown in Table 3, and the gradation curve is demonstrated in Figure 4.

According to Table 3 and Figure 4, the particle size distribution of fine sand samples is uniform and poor, mostly in the range of 0.6–0.15, and the composition is mostly composed of medium sand and fine sand.

###### 3.1.3. Mineral Composition Analysis of the Muck

Multipurpose X-ray diffractometer (XRD) is applied to analyze the phase of silty clay and fine sand. The mineral composition of silty clay in muck is mainly composed of quartz (21.35%), feldspar (16.38%), and montmorillonite (31.03%). The mineral composition of fine sand mainly consists of quartz (35.03%) and feldspar (27.35%).

##### 3.2. Feasibility Analysis of Shield Muck as Grouting Material

The clay mineral composition of the silty clay in muck is high with favorable water absorption and expansion properties. Besides, the particle size of soil sample is small, and the specific surface area of soil particles per unit mass is large, which can effectively improve the viscosity and stability of the grouting liquid. Therefore, it is feasible to replace the bentonite in the raw materials of synchronous grouting liquid based on the features like physical properties, mineral composition, and particle size composition of the muck in the section of Zhengzhou Metro. Moreover, the mineral components of the fine sand in muck are mainly quartz and feldspar with little mud content. It has the basic conditions for recycling and the utilization as engineering sand, so it is feasible to replace the fine sand in the raw materials of synchronous grouting liquid.

#### 4. Optimization of Ratio of the Reused Muck Grout

Based on the engineering situation of the shield section, cement, fly ash, silty clay in muck (substitute for bentonite), fine sand in muck (substitute for river sand), bentonite, and water are chosen as single-liquid grouting materials with active cement mortar.

By using the basic theory and method of uniform design, a large number of tests can be effectively avoided. Uniform design was initiated by Chinese statistician Professor Fang Kaitai and member of Chinese Academy of Sciences Wang Yuan, and it is the premier method to conduct multifactor and multilevel experimental design. It can complete the research and development of complex scientific research projects and new products with fewer experiments. The biggest feature is that the number of experiments can be equal to the maximum level number, rather than the square of the number of experimental factors. Moreover, the number of experiments is only related to the number of factors to be investigated. However, the number of trials should be about 3 times the number of experimental factors, which is conducive to modeling and optimization.

The synchronous grouting liquid is a kind of hard single liquid, and its influencing factors can be generally divided into four properties, namely, water-binder ratio, binder-sand ratio, bentonite-water ratio, and fly ash-cement ratio. The water-binder ratio is the ratio of water to cementitious materials. The bentonite-water ratio is the ratio of bentonite to water. The cement-sand ratio is the ratio of cementitious materials to the sand. The fly ash-cement ratio is the ratio of fly ash to cement. To analyze the influence of each influencing factor on the construction performance of the grout and design the optimal mix proportions, five levels of tests were set for each influencing factor, and the number of tests was three times the actual level of the test. Therefore, 15 groups of tests were arranged, namely, the uniform design of 4 factors and 15 levels. Referring to Appendix I of Fang Kaitai’s *uniform design and uniform design table* [16], the uniform design was carried out by choosing columns 1, 2, 4, and 6 of U15 ∗ (157). The material proportions of the uniform design are listed in Table 4. The items of the grouting material performance test mainly include specific gravity test, consistency test, fluidity test, bleeding rate test, setting time test, stone shrinkage test, and unconfined compression test. The test results are presented in Table 5.

Based on the results of uniform test, SPSS software was used for regression analysis with regression method, and quadratic multiple regression model was adopted to analyze the interaction among variables.where *x*_{i} and *x*_{k} are test factors, *b*_{0}, *b*_{k}, and *b*_{jk} are coregression coefficients, *x*_{1} is the water-binder ratio, *x*_{2} is the binder-sand ratio, *x*_{3} is the bentonite-water ratio, and *x*_{4} is the fly ash-cement ratio. The regression equation of each factor is obtained according to the test results.

Based on the results of regression analysis, each influencing factor of grouting material and each performance of grout are analyzed. The factor was changed within its value range, and the influences of this factor on the specific gravity of mortar, stone shrinkage rate, initial fluidity, 3 h bleeding rate, setting time, and compressive strength are studied. Simultaneously, the range of independent variables in the multiobjective programming problem is determined according to the level interval of uniform test factors, and the fmincon function in MATLAB optimization toolkit is used to solve the multiobjective programming problem with ideal point method. The optimal proportions of synchronous grouting material prepared from shield muck are listed in Table 6. The performances of the grouting material with the optimal proportion are listed in Table 7.

#### 5. Numerical Simulation Research on the Effect of Shield Muck Synchronous Grouting

##### 5.1. Calculation Model and Parameters

The 201 ring to 210 ring of the left line in the interval are used as the test section of muck reuse. The tunnel section at the 205 ring of the left line is chosen as the calculation section to carry out the sensitivity analysis of synchronous grouting construction parameters by using FLAC3D finite difference software. Based on the geological survey data of this project, the main tunneling strata are fine sand and silty clay. The specific stratum parameters are shown in Table 8. The buried depth of underground stable water level is about 16 m, and the tunnel roof is 18 m away from the surface. The distance between the two lines at the section is about 13.5 m. Only the influence of grouting construction in the left line is considered in the calculation considering that the right line has been excavated.

According to the parameters of the shield machine, it can be calculated that the theoretical overbreak gap (25 mm) is half of the difference between the cutter diameter (6480 mm) and the outer diameter of the shield tail. The theoretical gap (70 mm) in shield tail is half of the difference between the inner diameter of the shield tail and the outer diameter of the segment (6200 mm). The theoretical thickness of the ring filled by grout is half of the difference between the cutter diameter and the outer diameter of the segment (140 mm).

Due to the randomness and complexity in the physical and mechanical properties of geotechnical materials, the following calculation conditions and assumptions are used in the numerical simulation: assuming that the thickness of each stratum is uniform and there is no mutual intrusion, surrounding rock and supporting structure are homogeneous and isotropic materials; the strata below groundwater depth are permeable; the water level remains unchanged during tunnel excavation; the grout is incompressible, isotropic, and nondilutable and does not block.

According to the Saint-Venant principle, the displacement and stress changes caused by excavation can be ignored outside the excavation area with 3 to 5 times the tunnel span. Therefore, the left, right, upper, and lower dimensions of the tunnel excavation area are not less than 3–5 times the tunnel span in the actual numerical simulation process. The calculation model is shown in Figure 5. The boundary constraints are the horizontal constraints on the front, back, left, and right of the model, the vertical constraints at the bottom of the model, and the free boundary on the top of the model. For the fluid boundary in seepage, the fluid exchange in the left, right, and bottom boundaries is constrained, and no fluid constraints are imposed on the front and back boundaries.

The stratum soil is simulated with solid element, and it obeys the M-C failure criterion. The shield segment is simulated with solid element, and it conforms to the elastic constitutive model. The segment joint is considered by setting the stiffness reduction coefficient as 0.7. The steel plate of shield machine is simulated by shell element in line with the elastic constitutive model. The grouting layer is simulated by solid element in accordance with elastic constitutive model and is regarded as permeable material. Since the grout is used to fill the gaps between the excavation diameter and the shield body, the filling layer is simulated with solid element with the elastic constitutive model, and it is regarded as impermeable material.

##### 5.2. Simulation Steps of Construction

The calculation steps simulate the three stages of shield construction: (1) shield machine passing stage: the shield shell unit and the overexcavation filling layer in the model are activated; (2) segment out of the shield tail stage: the shield shell element and the overexcavated filling layer are passivated, and the shield tail grouting layer and the segment lining in the model are activated. The grouting pressure of the whole ring in shield tail is applied. Considering the shear movement and consolidation of the grout (the change rate of the grout parameters in the early stage of consolidation is huge), attention should be paid to change the hardening parameters of the grouting layer at the appropriate time and reduce the grouting pressure. (3) Grout hardening stage: time-dependent hardening parameters of grouting layer should also be applied at this stage.

To simulate the field construction more accurately, the three stages of shield construction are divided into six calculation steps (Table 9) according to the hardening process of grouting body considering the construction stage and the influence on surface deformation. In the calculation process, different stratum stress release coefficients are set for different construction steps.

##### 5.3. Calculation of Time-Dependent Parameters of the Grout

To obtain the time-dependent hardening parameters of the materials in this study, the relationships between the elastic modulus and Poisson’s ratio of the grout and time are fitted according to the measurement results of the grout performance parameters at different ages, and the relationships are as follows:where *E* is elastic modulus, MPa, *λ* is Poisson’s ratio, and *t* is the curing days of the grout for 0–28 d.

The hardening parameters of the grout are obtained (Table 10) according to the consolidation time point in construction simulation steps 2∼6.

##### 5.4. Construction Effect Analysis under Different Grouting Pressure

###### 5.4.1. Range and Vertical Distribution of Grouting Pressure

It is assumed that the grouting pressure and grout are evenly distributed in the whole ring with a diffusion way of filling method. The filling rate of the shield tail void is 100%. To reasonably and fully study the influence of grouting pressure on tunnel structure and surface deformation, the range of uniform grouting pressure in the whole ring is set to 0.10 MPa∼0.30 MPa, and the linear gradient is 0.04 MPa with a total of 6 working conditions.

Due to the circumferential and longitudinal shear movement of the grout when grouting, and the grout consolidating over time, it is considered that the grouting pressure is almost reduced to 0 within 1 d after the detachment of the ring (i.e., 9 m away from the shield tail) combined with the engineering practice and literature research. The reduction coefficients of longitudinal grouting pressure are listed in Table 11.

###### 5.4.2. Analysis of Numerical Simulation Results

*(1) Surface Settlement*. The surface settlement under different grouting pressures is shown in Figure 6. The lateral surface settlement curves under different grouting pressures are basically the same and nearly in line with the characteristics of Gaussian curve. The change of grouting pressure has a great influence on surface settlement. The maximum surface settlement gradually decreases with the increase of grouting pressure, but the width of settlement trough remains unchanged. When the grouting pressure is greater than 0.22 MPa, the uplift occurs at the surface about 25 m away from the axis on both sides of the shield, and the uplift increases with the increase of grouting pressure.

As shown in Figure 7, the segment deforms elliptically and sinks as a whole under the effect of water-soil pressure and grouting pressure. The segment deformation of the feature points changes greatly under different grouting pressure, and the overall settlement of the segment decreases with the increase of grouting pressure. The change values of segment convergence are manifested in Figure 8. With the increase of grouting pressure, the vertical convergence of segment gradually increases, while the horizontal convergence of segment shows a trend of first increasing and then decreasing. When the grouting pressure is set as 0.18 MPa, the horizontal convergence of segment reaches the maximum value of 4.16 mm.

*(2) Distribution of Plastic Zone*. Table 12 shows the distribution of the soil plastic zone under different grouting pressures. The plastic zone is mainly distributed on both sides of the tunnel, and the size decreases with the increase of grouting pressure. When the value of grouting pressure is less than 0.18 MPa, there are stretched failure units around the tunnel.

The control value of surface settlement in tunnel monitoring project is 10 mm according to the technical specifications of monitoring on urban rail transit engineering and the first-level requirements of engineering monitoring. Figure 6 shows that the minimum grouting pressure is 0.18 MPa. Figure 8 presents that, with the increase of grouting pressure, the convergence of segment clearance gradually increases. Excessive grouting pressure can easily lead to segment dislocation and cracking, so the uniformly distributed grouting pressure in the whole ring should be set to 0.18 MPa.

##### 5.5. Construction Effect Analysis with Different Grouting Quantity

###### 5.5.1. Filling Grouting

*(1) Simulation Method*. The convergent units are set up in the calculation model in order to simulate the effect that the grout incompletely fills the gaps of the shield tail. As demonstrated in Figure 9, when the segment is just separated from the shield tail, the soil tends to converge to the tunnel. It is assumed that when the grout fills 80% of the shield tail clearance, the remaining 20% is occupied by the convergent soil. To ensure the continuity of the grid in the numerical simulation, the convergent units are used to replace the empty units, and these units are given small strength parameters for approximate simulation.

**(a)**

**(b)**

**(c)**

*(2) Working Conditions*. Assuming the filling and grouting rate is 80%∼100%, the specific conditions are shown in Table 13.

*(3) Analysis of Numerical Results* ① Surface settlement. The surface settlement under different backfill grouting conditions is shown in Figure 10. The curves of horizontal surface settlement under different backfill grouting rates are basically the same, which is nearly in line with the feature of Gaussian. The change of backfill grouting rate has a great influence on surface settlement. With the increase of backfill grouting rate, the maximum value of surface settlement decreases gradually, and the variation is considerable. The width of the setting tank remains unchanged when the backfill grouting rate changes. ② Segment deformation. Figure 11 shows the deformation of segment under different backfill grouting conditions. The segment deforms elliptically and sinks as a whole under the effect of water-soil pressure and grouting pressure. The segment deformation of the feature points changes greatly under different grouting pressure, and the overall settlement of the segment decreases with the increase of grouting pressure. The relationship between backfill grouting rate and segment convergence is portrayed in Figure 12. With the increase of backfill grouting rate, the vertical convergence and horizontal convergence of the segment show an increasing trend; that is, the ellipticity of segment increases with the increase of backfill grouting rate. ③ Plastic zone distribution. Table 14 is the distribution of soil plastic zone without backfill grouting. The plastic zone is mainly distributed on the left and right sides of the tunnel, and the size decreases with the increase of filling grouting rate. The plastic deformation areas above the left and right of the tunnel are large when the backfill grouting rate is less than 100%, and it is easy to cause the instability of the soil above the tunnel.

The control value of surface settlement in tunnel monitoring project is 10 mm according to the technical specifications of monitoring on urban rail transit engineering and the first-level requirements of engineering monitoring. Figure 10 shows that when the backfill grouting rate is less than 100%, the maximum surface settlement is greater than 10 mm. Considering the surface settlement and plastic zone distribution, the backfill grouting rate should be 100%.

###### 5.5.2. Penetration Grouting

The influences of the change of grout injection rate on surface settlement, segment deformation, and plastic zone distribution are simulated simultaneously, and the grout injection rate of 100%–200% is set for calculation and analysis. It is assumed that when the grout injection rate exceeds 100%, penetration grouting happens, and the whole ring is uniformly penetrated. Meanwhile, based on the method in Lei et al. [17], the elastic modulus of grout-soil mixture can be roughly calculated according to (3), and Poisson’s ratio is adopted according to grout of corresponding calculation stage. The specific conditions are presented in Table 15.where *E*_{3} is the elastic modulus of grout-soil mixture, MPa, *E*_{1} is the elastic modulus of the stratum soil where the tunnel is located, MPa, *E*_{2} is the elastic modulus of the corresponding grout, MPa, and *n*_{f} is the effective porosity of the soil stratum where the tunnel is located, which is 25.44%.

The simulation results are shown in Figures 13-15. The simulation results indicate that the change of grout injection rate has little effect on surface settlement, segment deformation, and plastic zone distribution. With the increase of grout injection rate, the maximum surface settlement decreases. The ellipticity of the segment increases and then decreases with the increase of grout injection rate. When the grout injection rate is 120%, the vertical convergence and horizontal convergence of the segment reach the maximum.

#### 6. Conclusions

(1)The mineral composition of silty clay in the shield muck is mainly composed of quartz (16.38%), feldspar (21.35%), and montmorillonite (31.03%). The mineral composition of fine sand mainly consists of quartz (35.03%) and feldspar (27.35%). It is feasible to use the silty clay and fine sand in shield muck to replace the bentonite and fine sand in raw materials of synchronous grouting, respectively.(2)Based on the uniform test results, the SPSS multiple regression analysis and MATLAB are used to solve the nonlinear programming problem. The optimal proportions of shield muck synchronous grouting material are water-cement ratio of 0.74, binder-sand ratio of 0.84, bentonite-water ratio of 0.11, and fly ash-cement ratio of 2.75.(3)The finite difference method is used to analyze the influences of the grouting pressure and the filling rate by considering the combined action of the seepage field, the stress field, and the time-dependent properties of the grout. From the perspective of calculation results and economic benefits, when the grouting pressure is constraint at 0.18 MPa and the grouting rate is arranged at 120%∼150%, the surface settlement and the deformation of the segment are better controlled when muck grout is adopted.#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was sponsored by the State Key Program of National Natural Science of China (no. 51938005). The authors are grateful to this institution for the support.