工程科学与技术2025，Vol.57Issue(4)：39-51,13.DOI:10.12454/j.jsuese.202300948
非均匀结构堰塞坝溃决过程的3维数值模拟
Three-Dimensional Numerical Simulation of the Breaching Process of Landslide Dams with Heterogeneous Structures
摘要
Abstract
Objective Landslide dams are natural blockages formed by the rapid accumulation of slope failure materials,typically resulting from mass move-ments such as rockslides or debris avalanches.Due to gravitational sorting during downslope transport and complex topographic constraints in narrow valleys,these dams generally exhibit inherently three-dimensional(3D),spatially heterogeneous internal structures.However,most exist-ing breach simulation models assume material homogeneity for computational convenience,neglecting the influence of real-world structural non-uniformity.This simplification results in significant deviations in breach morphology and peak discharge predictions,ultimately impairing the reli-ability of hazard assessments and early warning systems for downstream communities.This study clarifies the fundamental influence of 3D mate-rial heterogeneity on the breach process of landslide dams,enhancing the predictive capabilities of numerical models in geohazard risk ma-nagement. Methods A novel 3D landslide dam breach model was developed by coupling Large Eddy Simulation(LES)with sediment mass conservation equations that incorporated phase transitions between solid and suspended states,accurately capturing the erosion dynamics in heterogeneous dams.The model considered the non-directional transport and deposition behavior of sediment under complex flow conditions.A Volume-of-Fluid(VOF)method was employed to simulate the evolution of the free surface and interface,enabling the precise tracking of water-sediment in-teractions.The model was validated through physical model experiments simulating three structural types:homogeneous,vertically heteroge-neous,and laterally heterogeneous dams.Key breach parameters,including incision rate,breach geometry,and outflow hydrographs,were com-pared to experimental data to verify the model's reliability and accuracy. Results For homogeneous dams,breach behavior varied significantly with the material grain size.Fine-grained dams exhibited rapid,layered ero-sion,forming triangular longitudinal profiles,with breach durations under 60 seconds and peak discharges reaching 3.32 L/s.In contrast,coarse-grained dams underwent multistage headcut erosion,requiring higher flow shear stress for sediment entrainment.This delayed breach develop-ment resulted in longer times to peak(up to 120 s)and reduced peak discharge(2.1 L/s).Reverse vortices formed at the headcut bases enhanced local scouring but exerted limited influence on overall erosion rates.Medium-grained dams exhibited intermediate characteristics,characterized by single-stage headcuts and corresponding breach metrics that fell between those of fine-and coarse-grained dams.As the median grain size in-creased,breach morphology transitioned from uniform scouring to progressively complex headcut erosion patterns,with delayed peak times and reduced peak flows.For vertically heterogeneous dams,the breach process was susceptible to the configuration of layered materials.The upper layer(V1)influenced breach initiation:in tests where it consisted of fine particles,breach formation was accelerated;in contrast,coarse-grained V1 layers delayed erosion onset and increased upstream impoundment volumes by up to 40％.The middle layer(V2)governed vertical incision rates:fine-grained layers accelerated downward erosion,whereas coarse layers formed headcuts that restricted further deepening.The lower layer(V3)controlled basal stability:coarse-grained foundations enhanced dam resistance to scouring,while fine-grained ones raised undercutting and subsequent collapse of the overlying mass.In addition,an interactive effect was observed between the various layer combinations.For example,"inverse grading"structures(coarse-over-fine)tended to form armor layers on the surface,delaying breach initiation and leading to sudden failure modes with elevated peak discharges.In contrast,"normal grading"(fine-over-coarse)favored headcut erosion.These structural patterns critically influenced both the erosion mechanisms and the hydraulic response.In laterally heterogeneous dams,spatial variability in material properties across dam zones significantly altered breach morphology and flow dynamics.The breach core zone(C1)determined incision and discharge char-acteristics,while adjacent zones(C2 and D1)influenced breach asymmetry and propagation direction.The Dl zone played a major role in control-ling breach widening,whereas D2 exhibited minimal influence.Comparative simulations(Tests 7～9)revealed differences:Test 7 yielded a wide breach,Test 8 produced a narrow and deep breach,and Test 9 exhibited slow incision but extensive lateral expansion.In addition,the deposition of coarse particles within the breach resulted in a reduction of the slope angle and flow velocity,which suppressed the initiation of further coarse particle entrainment and impeded upstream headward erosion.These feedback mechanisms resulted in a lower peak discharge and a more gradual evolution of the breach. Conclusion This study demonstrates that the spatial heterogeneity of dam materials is a primary factor governing the evolution of breaches in landslide dams.Vertically inverse grading structures(coarse-over-fine)contribute to delayed breach initiation and energy storage,often resulting in abrupt failure events with peak discharges exceeding those of homogeneous dams by more than 30％.In laterally heterogeneous dams,the sort-ing of materials on the overtopping side directly influences the breach depth-to-width ratio(ranging from 0.8 to 2.5),controlling the discharge ca-pacity.A positive feedback mechanism associated with coarse sediment deposition is identified:the accumulation of sediment reduces the breach slope,which in turn decreases flow velocity and further enhances deposition,transforming the breach mode from rapid incision to a slow-release erosion regime.The developed 3D VOF-LES-based hydro-sediment coupled model overcomes the limitations of the traditional homogeneous as-sumption and,for the first time,enables high-resolution simulation of breach morphology evolution under realistic heterogeneous conditions.With prediction errors maintained within 10％under complex experimental conditions,this model provides a robust tool for enhancing risk assess-ments and emergency planning in regions prone to landslide dam breaches.关键词
堰塞坝/非均质结构/3维数值模拟/溃决侵蚀/溃决流量Key words
landslide dam/heterogeneous structure/three-dimensional numerical simulation/breaching erosion/outflow discharge分类
天文与地球科学引用本文复制引用
胡宪睿,彭铭,付小莉,杨鸽,朱艳,石振明,张公鼎..非均匀结构堰塞坝溃决过程的3维数值模拟[J].工程科学与技术,2025,57(4):39-51,13.基金项目
国家自然科学基金项目(42061160480 （）
42071010 （）
U23A2044) （）
