The Simulation of Grouting Behavior in the Pea Gravel Filling Layer Behind a Double-Shield TBM Based on the Level Set Method

In double-shield TBM tunnel construction, grouting plays a vital role in consolidating the gravel backfill and maintaining the integrity of the segmental lining. To investigate the permeation behavior of grout within the pea gravel layer, a fluid dynamics model was developed in this study. The model directly simulates the flow of grout through the porous medium by solving the Navier–Stokes equations and employs the level set method to track the evolving interface between the grout and air phases. Unlike conventional continuum approaches, this model incorporates particle-scale heterogeneity, allowing for a more realistic analysis of grout infiltration through the non-uniform pore structures formed by gravel packing. Three different grouting port positions and two boundary conditions are considered in the simulation. The results indicate that under pressure boundary conditions, the grout flow rate increases rapidly in the initial stage, and then decreases and stabilizes, with the flow rate peak increasing as the grout port moves upward. Under velocity boundary conditions, the injection pressure grows slowly in the early stage but accelerates with time. Additionally, the rate of pressure change is faster when the grout port is located lower in the backfilling layer. Through theoretical analysis, the existing analytical formula was extended by introducing a gravitational correction term. When the grouting port is near the upper part of the tunnel, the analytical solution aligns well with the numerical simulation results, but as the grout port moves downward, the discrepancy between the two increases.