Dashiling Gully, located in Miyun District, Beijing, exhibits a high susceptibility to debris flow due to its unique geological and topographical characteristics. The area is characterized by well-developed rock joints and fissures, intense weathering, a steep gradient, and a constricted gully morphology. These factors contribute to the accumulation of surface water and loose sediment, significantly increasing the risk of debris flow events. Following a comprehensive field geological investigation of Dashiling Gully, key parameters for simulation were obtained, including fluid weight, volume concentration, and rainfall. The formation and development conditions of potential mudslides were analyzed, and numerical simulations were conducted using FLO-2D software (version 2009) to assess scenarios with rainfall probabilities of 1 in 30, 50, and 100 years. The simulations accurately reconstructed the movement velocity, deposition depth, and other critical movement characteristics of mudslides under each rainfall scenario. Using ArcGIS, pre- and post-treatment hazard zoning maps were generated for Dashiling Gully. Furthermore, the efficacy of implementing a retaining wall as a mitigation measure was evaluated through additional numerical simulations. The results indicated that mudslide velocities ranged from 0 to 3 m/s, with deposition depths primarily between 0 and 3 m. The maximum recorded velocity reached 3.5 m/s, corresponding to a peak deposition depth of 4.31 m. Following the implementation of the retaining wall, the maximum deposition depth significantly decreased to 1.9 m, and high-risk zones were eliminated, demonstrating the intervention’s effectiveness. This study provides a rigorous evaluation of mudslide movement characteristics and the impact of mitigation measures within Dashiling Gully. The findings offer valuable insights and serve as a reference for forecasting and mitigating similar mudslide events triggered by heavy rainfall in gully mudslides.