本发明涉及自然灾害防治领域，公开了一种滑坡‑堵江‑堰塞坝‑溃决洪水灾害链模拟方法，本发明充分考虑了滑坡‑碎屑流在沿沟道内运动的高速流动性及滑坡碎屑流入江时与水流的相互作用，通过第一流体力学计算模型对滑坡碎屑流堵江流体力学过程进行数值计算模拟，得到滑坡碎屑流时空演化的具体分布形态；在堵江形成堰塞坝之后，考虑了坝体物质的非均匀性及溃决过程中水土物质的相互作用，通过第二流体力学计算模型对堰塞坝溃决流体力学过程进行数值计算模拟，得到堰塞坝溃决流体时空演化的具体分布形态；因此，本发明预测结果计算精度高；同时有利于为滑坡‑堵江‑堰塞坝‑溃决洪水灾害链的风险预测提供定量评价和科学支撑。