第47卷第2期2023年6月doi:10.3969/j.issn.1674-3636.2023.02.012雅鲁藏布江色东普沟冰崩一堵江一溃决灾害数值模拟田杨杨，姜亮²，郭江²（1.西藏地勘局区域地质调查大队，西藏拉萨851414；2.四川省地质环境调查研究中心，四川成都610036)摘要：为了揭示雅鲁藏布江色东普沟2018年10月17日冰崩一堵江一溃决灾害链的动力演化过程，基于Massflow数值模拟仿真平台，使用Fortran编程语言，根据研究区域地质条件特征对程序进行二次开发以优化Voellmy模型，模拟冰崩一泥石流动力过程；将模拟泥石流得到的堰塞坝体嵌入地形中，运用ArcGIS计算堰塞湖范围及体积，通过Manning模型模拟堰塞湖溃决洪水动力过程。采用分段模拟法再现冰崩一泥石流一堵江一堰塞湖一溃坝的完整动力过程，对泥石流运动过程中的流速、流深，坝体高度，溃决洪水的流深、流速等参数进行定量化研究，为色东普流域的防灾减灾工作提供有效支撑。关键词：冰崩；泥石流；堰塞湖；灾害链;Massflow数值模拟；色东普沟;雅鲁藏布江中图分类号：P642.230引言堵江事件是指高海拔斜坡上的冰川在重力作用下断裂解体后迅速滑塌、铲刮沟道并裹挟碎屑物堵塞河道形成堰塞坝，在坝体溃决后又引发洪水的灾害,属于典型的冰崩一泥石流一堵江一溃决链式灾害。为揭示该链式灾害的动力过程，对灾害链进行模拟反演。基于动力过程的数值模型主要有连续介质力学模型、离散介质模型及混合介质模型（王钟文，2020）。Savage等（1989）最早提出了深度积分连续介质力学的方法概念，并在干性沙土的流动试验中得到了验证。连续介质运动方向特征尺度远大于垂直方向特征尺度，使得复杂的三维地表过程可通过二维计算求解。该方法提出以来,崩滑碎屑流、洪水、泥石流等地表过程动力学数值模拟研究均取得了突破性进展。深度积分连续介质力学方法在泥石流等山地灾害动力学过程分析中具有独特优势：①可输出多种泥石流特征参数，不仅能为定量化评价提供更多的参考依据,而且可以与已有的野外监测数据、大规模试验数据进行对比分析,提高模型的准确性(Iverson收稿日期：2022-09-30；修回日期：2022-11-09；编辑：陈露，朱明君基金项目：西藏自治区自然资源项目“监测示范工程"（SCXX-XZCG-20098）作者简介：田杨杨（1981一），男，工程师，主要从事地质矿产、地质灾害动力演化过程研究，E-mail：35155447@qq.com通信作者：姜亮（1994—），男，助理工程师，硕士，主要从事地质灾害动力演化特征研究，E-mail：1130183255@qq.com地质学...
