Numerical Simulation of Rainfall-Induced Erosion on Infiltration and Slope Stability

Abstract

In slopes where a mixture of coarse and fine particles is present, the infiltration of rainfall can cause the migration of fine particles. This migration alters the hydraulic properties of the soil and has implications for slope stability. In this study, the slope under investigation is a tailings dam composed of loosely consolidated soil with a wide particle size distribution. Due to rainfall infiltration, fine particles tend to migrate within the voids of the coarse particle framework, leading to changes in hydraulic properties and inducing slope instability. The classical internal erosion constitutive model, known as the Cividini and Gioda erosion criterion, is commonly used to predict the behavior and effects of fine particle erosion in geotechnical engineering. However, certain parameters in this erosion criterion equation, such as long-term density, are challenging to obtain through experiments. To investigate the coupled evolution of seepage and erosion within landfill slopes under the influence of rainfall infiltration and to understand the mechanisms of slope instability, this research assumes the erosion of fine particle suspension and adopts the Worman and Olafsdottir erosion criterion to establish a coupled model of unsaturated seepage and internal erosion. The developed model simulates the coupled response of seepage and erosion in unsaturated landfill slopes under three different rainfall intensities. It is then combined with the infinite slope model to quantitatively analyze the impact of fine particle migration on soil permeability and slope stability. The numerical simulations provide the following findings: The Worman and Olafsdottir erosion criterion, unlike the Cividini and Gioda erosion criterion, only requires the determination of the soil’s gradation curve to estimate the erosion rate. Internal erosion primarily occurs within the leading edge of moisture penetration, accelerating the advancement of the wetting front and reducing slope stability. When the rainfall intensity is lower than the saturated permeability coefficient, the influence of internal erosion can be disregarded. However, under rainfall intensities equal to or greater than the saturated permeability coefficient, considering internal erosion results in a difference in the depth of the wetting front of up to 34.2 cm after 6 h in the R2 scenario. The safety factor without considering internal erosion is 1.12, whereas considering internal erosion yields safety factors between 1.08 and 1.09. In the R3 scenario, the difference in the depth of the wetting front reaches 53.8 cm after 6 h, with a safety factor of 1.12 without considering internal erosion and safety factors between 1.06 and 1.07 when considering internal erosion.