## 1 Introduction  
With the development of society and the economy, it is inevitable that a high fill slope will be formed through engineering construction. Compared to general slopes, high fill slopes are more susceptible to engineering disasters such as debris flows, landslides, and collapses due to external factors. The instability mechanism of these slopes is also more complex. Therefore, it is essential to conduct a detailed study and analysis of the stability of high fill slopes based on the specific conditions of the project. Currently, numerous scholars have conducted extensive research on the determination of fill slope parameters and stability analyses. The diversity in analyzing slope geotechnical parameters underscores the fundamental role of determining these parameters for stability analyses [1; 2; 3]. The main methods for slope stability analyses include the limit equilibrium method, based on the rigid body theory of equilibrium, and the numerical analysis method, based on finite element (EM) and finite difference (FDM) [4; 5]. In a paper by [PERSON] et al. [6], the deformation process of high-loess-filled slopes under different rainfall disaster mitigation measures was studied using the Flume Test (which includes the following three parts: an experimental slope system, a rainfall system, and a data acquisition system). It was found that the deformation process of the slope is highly correlated with the variation of the pore water pressure. When the slope fails, the pore water pressure at each position in the slope will experience a sudden fluctuation in response to the instantaneous excitation. [PERSON] et al. [7], using three types of rainfall combined with actual rainfall events and data, analyzed the stability of high embankment fill bodies on slopes, considering the effects of water deterioration. According to the simulation analyses, the negative pore water pressure (matrix suction) exhibits a greater area of dissipation under the multi-peak type of rainfall with the assistance of the antecedent rainfall, and the evident decrease in safety factors with the smallest value is also found with the upward movement of the potential sliding surface. When considering the deterioration effect of water, the stability of embankment slopes is more easily affected by the infiltration of rainwater, and the position of the potential sliding surface is more consistent with the dispersal of matrix suction. [PERSON] et al. [8], combining indoor tests and simulation calculations, analyzed the embankment slope of a certain area in Baoshan, Yunnan Province, which contains a soft weak layer, focusing on the influence of hygroscopic swelling on the embankment slope of expansive soil. A proposal was made to improve the stability of the slope by directly reinforcing the soft weak layer with rigid piles. [PERSON] et al. [9] studied the stability of a vertical slope on a high embankment at a certain site using the Geo-Studio, Liren software, and Ansys software, and conducted stability analyses using the M-P method, Bishop method, and strength reduction method. [PERSON] et al. [10] used the unbalanced thrust method, simplified Bishop strip method, and finite element strength reduction method to compare the stability of a typical section of the Wuhan Shennvfeng Airport high embankment under four conditions of natural, rainstorm, cleared natural, and cleared rainstorm, and the results showed that the engineering measures of clearing the embankment surface could significantly enhance the stability of the high embankment. [PERSON] et al. [11] conducted centrifugal model tests on a high embankment using an improved, generalized proportionality method to study its dynamic response. The study found that the large dynamic displacements near the slope surface caused the development of lateral cracks.