Page 111 - 2024年第55卷第10期

形。综上所述，库水位下降至 １６５ｍ时为滑坡阶跃变形的启动阈值。
表 ２ ６次阶跃变形对应库水位高度情况

变形次序 变形启动时间 变形结束时间 开始变形时库水位高度?ｍ 累积位移?ｍｍ
１ ２０１７ － ３ － ２２ ２０１７ － ９ － １３ １６５．２２ １５５．６
２ ２０１８ － ３ － ２６ ２０１８ － ９ － １２ １６３．５８ １１４．９
３ ２０１９ － ４ － ２ ２０１９ － ８ － ２５ １６５．２０ ５３．１

４ ２０２０ － ３ － １２ ２０２０ － ９ － １８ １６６．５２ １２２．３
５ ２０２１ － ３ － １６ ２０２１ － ９ － ２５ １６５．５４ ４４１．６
６ ２０２２ － ４ － ２０ ２０２２ － ９ － ２０ １６５．３５ １４５．８

图 ６ 滑坡阶跃变形启动时对应库水位高度

４ 滑坡涌浪分析

４．１ 数值模型构建 前文分析指出，黑石板滑坡运动模式为牵引式，前缘失稳后，后缘失去阻滑力随
之失稳。１６５ｍ库水位是滑坡体下滑的启动阈值。本文考虑最不利情况，即滑坡整体失稳且库水位高
度在 １６５ｍ的情况下激发的涌浪特征和致灾范围，以得到相对保守的计算结果。采用基于计算流体动
力学的 ＦＬＯＷ－ ３Ｄ软件开展滑坡失稳入水的流固耦合计算。黑石板滑坡为土质滑坡，失稳运动过程中
必然会破碎解体，呈现出类流体的流动特征，因此选用颗粒流模型来模拟滑坡体。颗粒流模型采用密
度大于水的类牛顿流体来模拟滑坡体，可以较为准确表达滑坡运动过程中的解体、类流体运动和激浪
特征，该模型已被用于多个滑坡运动冲击过程，取得了良好的应用效果 ［２３，２９］ 。采用 ＲＮＧ湍流模型来
模拟水体，相比于标准 ｋ － ε模型，ＲＮＧ湍流模型 能够更 好 地模 拟 高应 变 率 和流 向 弯 曲程 度较大的
流动 ［３０］ 。
基于滑坡区域的 ＤＥＭ数据，建立黑石板滑坡在 １６５ｍ库水位条件下的数值计算模型（图 ７）。计算
域选取长和宽分别为 ５０８０和 ２７０９ｍ的部分山谷，网格大小设置为 ８ｍ。在涌浪生成和传播过程中设
置 Ａ—Ｄ共 ５个监测点，其中监测点 Ａ位于滑坡与水体的交界处，用以记录滑坡失稳入水后的激浪特
征，监测点 Ｂ位于坡体对岸。剖面Ⅰ － Ⅰ沿滑坡主滑方向设置，并经过监测点 Ａ、Ｂ。监测点 Ｃ和 Ｄ沿
２
河道设置，分别位于良斗河的上下游。整个系统采用垂直重力，重力加速度为 ９．８ｍ?ｓ，滑坡初速度
设置为 ０ｍ?ｓ，计算时间设置为 ６６０ｓ。
４．２ 参数与边界条件 滑坡的计算参数基于勘查数据和室内试验成果确定。流体设置为纯水，密度为
－ １
３
１０００ｋｇ?ｍ ，水的黏度为 ０．００１ｋｇ?（ｍ·ｓ ），其他参数如表 ３所示。

２
— １ ４ １ —