《工程科学学报》：降雨和地震条件下浅层黄土滑坡三维稳定性评价

工程科学学报 Chinese Journal of Engineering 降雨和地晨条件下浅层黄土滑坡三维稳定性评价 李泊良张帆宇 Three-dimensional stability evaluation of shallow loess landslides under rainfall and earthquake conditions LI Bo-liang.ZHANG Fan-yu 引用本文： 李泊良，张帆宇.降雨和地震条件下浅层黄土滑坡三维稳定性评价.工程科学学报，2022,44(3)：440-450.doi: 10.13374j.issn2095-9389.2020.10.08.002 LI Bo-liang,ZHANG Fan-yu.Three-dimensional stability evaluation of shallow loess landslides under rainfall and earthquake conditions[J].Chinese Journal of Engineering,2022,44(3):440-450.doi:10.13374/j.issn2095-9389.2020.10.08.002 在线阅读View online::htps:/ldoi.org10.13374.issn2095-9389.2020.10.08.002 您可能感兴趣的其他文章 Articles you may be interested in 三维确定性模型在浅层黄土滑坡稳定性预测中的应用 Application of a 3D deterministic model for predicting shallow loess landslide stability 工程科学学报.2018,40(4：397 https:1doi.org10.13374j.issn2095-9389.2018.04.002 基于MSR300雷达监测的凹山采场降雨条件下的边坡变形及滑坡 Investigation of deformation and failure in washan slope considering rainfall conditions based on MSR300 radar monitoring 工程科学学报.2018.40(4：407htps:/doi.org10.13374.issn2095-9389.2018.04.003 强降雨条件下孔隙气压作用的高台阶排土场渗流与稳定性 Seepage and stability analysis of pore air pressure on a high-bench dump under heavy rainfall 工程科学学报.2021,43(3)：365 https::/1doi.0rg/10.13374j.issn2095-9389.2020.09.01.005 渗透剪切作用下黄土的力学特征 Mechanical characteristics of loess under seepage shear 工程科学学报.2018,40(5：639 https:1doi.org10.13374.issn2095-9389.2018.05.015 三维软硬互层边坡的破坏模式与稳定性研究 Numerical analysis of the failure modes and stability of 3D slopes with interbreeding of soft and hard rocks 工程科学学报.2017,392：182 https::/1doi.org/10.13374.issn2095-9389.2017.02.003 变渣皮厚度条件下铜冷却壁应力分布规律及挂渣稳定性 Stress distribution law and adherent dross stability of the copper cooling stave with variable slag coating thickness 工程科学学报.2017,393：389htps:/doi.org10.13374issn2095-9389.2017.03.011

降雨和地震条件下浅层黄土滑坡三维稳定性评价 李泊良 张帆宇 Three-dimensional stability evaluation of shallow loess landslides under rainfall and earthquake conditions LI Bo-liang, ZHANG Fan-yu 引用本文: 李泊良, 张帆宇. 降雨和地震条件下浅层黄土滑坡三维稳定性评价[J]. 工程科学学报, 2022, 44(3): 440-450. doi: 10.13374/j.issn2095-9389.2020.10.08.002 LI Bo-liang, ZHANG Fan-yu. Three-dimensional stability evaluation of shallow loess landslides under rainfall and earthquake conditions[J]. Chinese Journal of Engineering, 2022, 44(3): 440-450. doi: 10.13374/j.issn2095-9389.2020.10.08.002 在线阅读 View online: https://doi.org/10.13374/j.issn2095-9389.2020.10.08.002 您可能感兴趣的其他文章 Articles you may be interested in 三维确定性模型在浅层黄土滑坡稳定性预测中的应用 Application of a 3D deterministic model for predicting shallow loess landslide stability 工程科学学报. 2018, 40(4): 397 https://doi.org/10.13374/j.issn2095-9389.2018.04.002 基于MSR300雷达监测的凹山采场降雨条件下的边坡变形及滑坡 Investigation of deformation and failure in washan slope considering rainfall conditions based on MSR300 radar monitoring 工程科学学报. 2018, 40(4): 407 https://doi.org/10.13374/j.issn2095-9389.2018.04.003 强降雨条件下孔隙气压作用的高台阶排土场渗流与稳定性 Seepage and stability analysis of pore air pressure on a high-bench dump under heavy rainfall 工程科学学报. 2021, 43(3): 365 https://doi.org/10.13374/j.issn2095-9389.2020.09.01.005 渗透剪切作用下黄土的力学特征 Mechanical characteristics of loess under seepage shear 工程科学学报. 2018, 40(5): 639 https://doi.org/10.13374/j.issn2095-9389.2018.05.015 三维软硬互层边坡的破坏模式与稳定性研究 Numerical analysis of the failure modes and stability of 3D slopes with interbreeding of soft and hard rocks 工程科学学报. 2017, 39(2): 182 https://doi.org/10.13374/j.issn2095-9389.2017.02.003 变渣皮厚度条件下铜冷却壁应力分布规律及挂渣稳定性 Stress distribution law and adherent dross stability of the copper cooling stave with variable slag coating thickness 工程科学学报. 2017, 39(3): 389 https://doi.org/10.13374/j.issn2095-9389.2017.03.011

工程科学学报.第44卷，第3期：440-450,2022年3月 Chinese Journal of Engineering,Vol.44,No.3:440-450,March 2022 https://doi.org/10.13374/j.issn2095-9389.2020.10.08.002;http://cje.ustb.edu.cn 降雨和地震条件下浅层黄土滑坡三维稳定性评价 李泊良，张帆宇四 兰州大学土木工程与力学学院，兰州730000 通信作者，E-mal:Zhangfy@lzu.edu.cn 摘要以兰州市大沙沟流域浅层黄土滑坡为研究对象，整合降雨渗透模型到黄土斜坡三维确定性模型，评价不同降雨和地 震耦合效应下浅层黄土滑坡稳定性，并用混淆矩阵法和受试者工作特征曲线法(ROC)评价稳定性评价预测的结果.研究结 果发现，耦合降雨入渗和地震的三维确定性模型，对流域尺度浅层黄土滑坡稳定性评价具有较好效果，能作为降雨和地震诱 发黄土滑坡灾害评价和早期预警的工具，对加强不同尺度极端事件下黄土滑坡灾害时空灾害评价和预测具有重要参考价值. 关键词浅层黄土滑坡：地震：降雨：渗透模拟：稳定性评价 分类号P642.22 Three-dimensional stability evaluation of shallow loess landslides under rainfall and earthquake conditions LI Bo-liang,ZHANG Fan-yu College of Civil Engineering and Mechanics,Lanzhou University,Lanzhou 730000,China Corresponding author,E-mail:Zhangfy@lzu.edu.cn ABSTRACT Loess is widely distributed in the Northwest Plateau of China.One-third of the landslides in China occur in the loess area. Shallow loess landslides are especially widespread and frequent geological disasters,causing serious casualties and huge property damage.Under rainfall and loading,loess is prone to structural collapse and strength reduction.Therefore,shallow loess landslides distribute widely and occur frequently.Usually,rainfall and earthquakes are the frequent and active triggers for loess landslides.In recent years,a large number of loess landslides have been induced by the coupling of rainfall and earthquakes on the Loess Plateau. Although the coupling effect of earthquake and rainfall will seriously aggravate the instability probability and disaster risk of shallow loess landslides,there is still a lack of quantitative disaster evaluation research on such landslide events.This study chose the shallow loess landslide as the research object in the Dashagou catchment of Lanzhou city.The rainfall penetration model was integrated into a three-dimensional deterministic model of the loess slope,and the stability of the shallow loess landslide was evaluated in the study area with different rainfall and seismic coupling effects.The confusion matrix and the receiver operating characteristic(ROC)curve were used to evaluate the results of the stability evaluation prediction.Results of this study reveal that the integration of a three-dimensional deterministic model of rainfall infiltration and earthquake effects has a good impact on the stability evaluation of shallow loess landslides at the watershed scale.Moreover,this model can be used as a tool for the assessment and early warning of rainfall and earthquake- induced loess landslides.The employment of the three-dimensional deterministic model considering a complicated slope and rainfall situation has great significance in the acquisition of results that are more accordant with the actual situation.It is of great reference value to strengthen the spatiotemporal disaster assessment and prediction of loess landslide disasters under different scales of extreme events. KEY WORDS shallow loess landslide;earthquake;rainfall;infiltration simulation;stability evaluation 收稿日期：2021-09-24 基金项目：国家自然科学基金资助项目(42090053.41977212)：兰州大学中央高校基本科研业务费专项资金资助项目(zujbky-2021-ct04)

降雨和地震条件下浅层黄土滑坡三维稳定性评价 李泊良，张帆宇苣 兰州大学土木工程与力学学院, 兰州 730000 苣通信作者， E-mail: Zhangfy@lzu.edu.cn 摘    要    以兰州市大沙沟流域浅层黄土滑坡为研究对象，整合降雨渗透模型到黄土斜坡三维确定性模型，评价不同降雨和地 震耦合效应下浅层黄土滑坡稳定性，并用混淆矩阵法和受试者工作特征曲线法（ROC）评价稳定性评价预测的结果. 研究结 果发现，耦合降雨入渗和地震的三维确定性模型，对流域尺度浅层黄土滑坡稳定性评价具有较好效果，能作为降雨和地震诱 发黄土滑坡灾害评价和早期预警的工具，对加强不同尺度极端事件下黄土滑坡灾害时空灾害评价和预测具有重要参考价值. 关键词    浅层黄土滑坡；地震；降雨；渗透模拟；稳定性评价 分类号    P642.22 Three-dimensional stability evaluation of shallow loess landslides under rainfall and earthquake conditions LI Bo-liang，ZHANG Fan-yu苣 College of Civil Engineering and Mechanics, Lanzhou University, Lanzhou 730000, China 苣 Corresponding author, E-mail: Zhangfy@lzu.edu.cn ABSTRACT    Loess is widely distributed in the Northwest Plateau of China. One-third of the landslides in China occur in the loess area. Shallow loess landslides are especially widespread and frequent geological disasters, causing serious casualties and huge property damage. Under rainfall and loading, loess is prone to structural collapse and strength reduction. Therefore, shallow loess landslides distribute widely and occur frequently. Usually, rainfall and earthquakes are the frequent and active triggers for loess landslides. In recent years, a large number of loess landslides have been induced by the coupling of rainfall and earthquakes on the Loess Plateau. Although the coupling effect of earthquake and rainfall will seriously aggravate the instability probability and disaster risk of shallow loess landslides, there is still a lack of quantitative disaster evaluation research on such landslide events. This study chose the shallow loess landslide as the research object in the Dashagou catchment of Lanzhou city. The rainfall penetration model was integrated into a three-dimensional deterministic model of the loess slope, and the stability of the shallow loess landslide was evaluated in the study area with different rainfall and seismic coupling effects. The confusion matrix and the receiver operating characteristic (ROC) curve were used to evaluate the results of the stability evaluation prediction. Results of this study reveal that the integration of a three-dimensional deterministic model of rainfall infiltration and earthquake effects has a good impact on the stability evaluation of shallow loess landslides at the watershed scale. Moreover, this model can be used as a tool for the assessment and early warning of rainfall and earthquake￾induced loess landslides. The employment of the three-dimensional deterministic model considering a complicated slope and rainfall situation has great significance in the acquisition of results that are more accordant with the actual situation. It is of great reference value to strengthen the spatiotemporal disaster assessment and prediction of loess landslide disasters under different scales of extreme events. KEY WORDS    shallow loess landslide；earthquake；rainfall；infiltration simulation；stability evaluation 收稿日期: 2021−09−24 基金项目: 国家自然科学基金资助项目（42090053，41977212）；兰州大学中央高校基本科研业务费专项资金资助项目（lzujbky-2021-ct04） 工程科学学报，第 44 卷，第 3 期：440−450，2022 年 3 月 Chinese Journal of Engineering, Vol. 44, No. 3: 440−450, March 2022 https://doi.org/10.13374/j.issn2095-9389.2020.10.08.002; http://cje.ustb.edu.cn

李泊良等：降雨和地震条件下浅层黄土滑坡三维稳定性评价 441. 黄土遇水和加载下，极易发生结构性崩塌和 滑动体沿坡体方向移动，这与实际情况不完全相 强度降低，导致黄土滑坡广泛分布且频发.统计数 同1为了弥补这些二维确定性模型存在的不足， 据显示，我国有1/3的滑坡灾害发生在黄土地区. 有学者指出三维模型得到的预测结果能更符合实 已经造成大量的人员伤亡和财产损失-降雨和 际情况.例如Bromhead等P指出二维模型很难准 地震是浅层黄土滑坡最常见的诱因.降雨在黄土 确地表示岩土性质、地下水条件的空间变化和滑 高原诱发了大量浅层黄土滑坡，诸如兰州、天水、 动面的形状，三维模型则能更准确的表示这些条 宝鸡等地主要是降雨型滑坡B-6然而，这些地区 件.此外，Xie等2指出三维模型得出的结果要比 还是历史地震的高发区和强影响区.诸如1718年 二维方法更可靠.因此，三维确定性模型当前成为 甘肃通渭地震、1920年宁夏海原地震、1995年甘 滑坡稳定性评价中更为关注的方法 肃永登地震都引起了大量的黄土滑坡⑧-0作为滑 Scoops3D是美国地质勘探局开发的三维确定 坡的两大重要的诱因，两者耦合作用下能诱发更 性模型，用来评价和预测研究区内每个栅格的稳 大范围、更广分布、更强破坏的黄土滑坡事件.近 定性2.目前，Scoops.3D已经成功预测了沿海地区 年来，这类灾害链效应在黄土高原已有发生.2013年 火山灰地层的滑坡稳定性2而黄土地层和火山 天水持续降雨诱发大量浅层黄土滑坡后，该地区 灰地层具有一定的相似性，所以考虑利用该模型 又突遇7.22彰县岷县6.6级地震.震后调查显示诱 评价黄土滑坡的稳定性.此外，该模型还可将其他 发以浅层黄土滑坡为主的斜坡灾害2000多处-2 水文模拟软件输出结果作为直接参数输人，具有 已有研究表明，震前的极端降雨事件会严重加剧 将区域复杂水文条件耦合到稳定性评级模型中的 同震滑坡的发生概率和灾害程度因此，加强 优势I)Tran等2应用TRIGRS计算降雨过程中 降雨和地震耦合效应下滑坡灾害的评价，对于不 斜坡动态水文条件并输入到Scoops.3D,对韩国 同尺度滑坡的防灾、减灾、应急都具有重要的指 Umyeon山区进行降雨条件下的稳定性评价，实现 导作用.然而，这类研究在黄土高原滑坡灾害评价 了降雨诱发滑坡的实时早期预警.Scoops.3D模型 中仍不多见 还引入了水平地震加速度系数，用来计算地震条 滑坡灾害危险评价通常包括定性和定量两类 件下的滑坡稳定性，因此，采用Scoops3D三维确 方法.定性方法主要是基于研究区地质条件和专 定模型耦合降雨和地震条件，进而评价黄土地区 家的经验对滑坡进行评级，而定量方法是依据物 的滑坡稳定性， 理模型和数学模型对滑坡进行分析)]定性方法 本文首先介绍了TRIGRS模型和Scoops3D模 依赖专家的主观性，缺少重复性.尽管定量方法可 型的基本原理，并分析了模型各类参数在稳定性 不同程度消除主观性的影响，且具有更好的适用 评价中的敏感性.然后，选择兰州市大沙沟流域作 性，但是需要大量数据的支持和输入，尽管如此， 为研究对象，结合TRIGRS和Scoops3D评价降 定量方法已成为滑坡灾害评价中主要的方法剧 雨、地震和两者耦合条件下区内浅层黄土滑坡的 在各类定量评价的技术中，基于数字高程模型并 稳定性.最后，结合滑坡点分布图，用混淆矩阵法 结合岩土力学参数的确定性模型，已成为流域尺 和受试者工作特征(ROC)曲线法对预测结果进行 度滑坡灾害评价的重要方法 评价，检验了降雨与地震触因及其耦合条件下对 物理确定性模型是通过无限斜坡稳定性模 黄土高原流域尺度黄土滑坡稳定性评价结果的适 型，计算不同水文条件下的坡体稳定性，通常能获 用性 得准确客观的结果9根据计算斜坡的维度，可 1 模型方法 分为二维和三维确定性模型.目前大多数确定性 模型都是基于二维无限斜坡模型评估斜坡稳定 1.1降雨渗透场模拟 性，常见的有基于区域栅格的瞬时降雨人渗的 TRIGRS模型是一种基于栅格数据计算降雨 TRIGRS与考虑稳态水文分布的SHALSTAB202四 入渗的斜坡稳定性模型，该程序能计算降雨入渗 这些模型均可模拟降雨过程中土体中水的渗流 引起的孔隙水压力变化以及随之变化的安全系 场，评估降雨条件下斜坡的稳定性.尽管这些二维 数.该模型结合降雨入渗和地表径流的方法计 确定性模型能一定程度考虑地质条件、水文条件 算饱和与非饱和条件下的一维垂向渗流过程，能 和气候条件，还可较为准确预测滑坡稳定性，但是 够获得降雨过程中土体压力水头和含水率的空间 在计算过程中他们不考虑滑动面的方向，只假设 分布规律.本文采用了降雨入渗的非饱和无限斜

黄土遇水和加载下，极易发生结构性崩塌和 强度降低，导致黄土滑坡广泛分布且频发. 统计数 据显示，我国有 1/3 的滑坡灾害发生在黄土地区， 已经造成大量的人员伤亡和财产损失[1−2] . 降雨和 地震是浅层黄土滑坡最常见的诱因. 降雨在黄土 高原诱发了大量浅层黄土滑坡，诸如兰州、天水、 宝鸡等地主要是降雨型滑坡 [3−6] . 然而，这些地区 还是历史地震的高发区和强影响区[7] . 诸如 1718 年 甘肃通渭地震、1920 年宁夏海原地震、1995 年甘 肃永登地震都引起了大量的黄土滑坡[8−10] . 作为滑 坡的两大重要的诱因，两者耦合作用下能诱发更 大范围、更广分布、更强破坏的黄土滑坡事件. 近 年来，这类灾害链效应在黄土高原已有发生. 2013 年 天水持续降雨诱发大量浅层黄土滑坡后，该地区 又突遇 7.22 彰县岷县 6.6 级地震. 震后调查显示诱 发以浅层黄土滑坡为主的斜坡灾害 2000 多处[11−12] . 已有研究表明，震前的极端降雨事件会严重加剧 同震滑坡的发生概率和灾害程度[13−16] . 因此，加强 降雨和地震耦合效应下滑坡灾害的评价，对于不 同尺度滑坡的防灾、减灾、应急都具有重要的指 导作用. 然而，这类研究在黄土高原滑坡灾害评价 中仍不多见. 滑坡灾害危险评价通常包括定性和定量两类 方法. 定性方法主要是基于研究区地质条件和专 家的经验对滑坡进行评级，而定量方法是依据物 理模型和数学模型对滑坡进行分析[17] . 定性方法 依赖专家的主观性，缺少重复性. 尽管定量方法可 不同程度消除主观性的影响，且具有更好的适用 性，但是需要大量数据的支持和输入. 尽管如此， 定量方法已成为滑坡灾害评价中主要的方法[18] . 在各类定量评价的技术中，基于数字高程模型并 结合岩土力学参数的确定性模型，已成为流域尺 度滑坡灾害评价的重要方法. 物理确定性模型是通过无限斜坡稳定性模 型，计算不同水文条件下的坡体稳定性，通常能获 得准确客观的结果[17−19] . 根据计算斜坡的维度，可 分为二维和三维确定性模型. 目前大多数确定性 模型都是基于二维无限斜坡模型评估斜坡稳定 性 ，常见的有基于区域栅格的瞬时降雨入渗的 TRIGRS 与考虑稳态水文分布的 SHALSTAB[20−22] . 这些模型均可模拟降雨过程中土体中水的渗流 场，评估降雨条件下斜坡的稳定性. 尽管这些二维 确定性模型能一定程度考虑地质条件、水文条件 和气候条件，还可较为准确预测滑坡稳定性，但是 在计算过程中他们不考虑滑动面的方向，只假设 滑动体沿坡体方向移动，这与实际情况不完全相 同[23] . 为了弥补这些二维确定性模型存在的不足， 有学者指出三维模型得到的预测结果能更符合实 际情况. 例如 Bromhead 等[24] 指出二维模型很难准 确地表示岩土性质、地下水条件的空间变化和滑 动面的形状，三维模型则能更准确的表示这些条 件. 此外，Xie 等[25] 指出三维模型得出的结果要比 二维方法更可靠. 因此，三维确定性模型当前成为 滑坡稳定性评价中更为关注的方法. Scoops3D 是美国地质勘探局开发的三维确定 性模型，用来评价和预测研究区内每个栅格的稳 定性[26] . 目前，Scoops3D 已经成功预测了沿海地区 火山灰地层的滑坡稳定性[26] . 而黄土地层和火山 灰地层具有一定的相似性，所以考虑利用该模型 评价黄土滑坡的稳定性. 此外，该模型还可将其他 水文模拟软件输出结果作为直接参数输入，具有 将区域复杂水文条件耦合到稳定性评级模型中的 优势[27] . Tran 等[28] 应用 TRIGRS 计算降雨过程中 斜坡动态水文条件并输入 到 Scoops3D，对韩 国 Umyeon 山区进行降雨条件下的稳定性评价，实现 了降雨诱发滑坡的实时早期预警. Scoops3D 模型 还引入了水平地震加速度系数，用来计算地震条 件下的滑坡稳定性. 因此，采用 Scoops3D 三维确 定模型耦合降雨和地震条件，进而评价黄土地区 的滑坡稳定性. 本文首先介绍了 TRIGRS 模型和 Scoops3D 模 型的基本原理，并分析了模型各类参数在稳定性 评价中的敏感性. 然后，选择兰州市大沙沟流域作 为研究对象 ，结 合 TRIGRS 和 Scoops3D 评价降 雨、地震和两者耦合条件下区内浅层黄土滑坡的 稳定性. 最后，结合滑坡点分布图，用混淆矩阵法 和受试者工作特征（ROC）曲线法对预测结果进行 评价，检验了降雨与地震触因及其耦合条件下对 黄土高原流域尺度黄土滑坡稳定性评价结果的适 用性. 1    模型方法 1.1    降雨渗透场模拟 TRIGRS 模型是一种基于栅格数据计算降雨 入渗的斜坡稳定性模型，该程序能计算降雨入渗 引起的孔隙水压力变化以及随之变化的安全系 数[20] . 该模型结合降雨入渗和地表径流的方法计 算饱和与非饱和条件下的一维垂向渗流过程，能 够获得降雨过程中土体压力水头和含水率的空间 分布规律. 本文采用了降雨入渗的非饱和无限斜 李泊良等： 降雨和地震条件下浅层黄土滑坡三维稳定性评价 · 441 ·

442 工程科学学报，第44卷，第3期 坡模型，在计算过程中将基底面以上的土体分为 滑动面上的栅格安全系数.该模型最初是为了 非饱和区和毛细饱和区，并通过Richards方程计算 评估火山灰地层中的滑坡稳定性，在尼加拉瓜的 降雨过程中压力水头空间分布.Alvioli等29开发 卡西塔和圣克里斯托瓦尔火山等多个地点进行了 的TRIGRS2.1提供了输出空间分布的孔隙水压力 成功测试阿随后还应用到华盛顿等区域的滑坡 和体积含水量的功能，并可直接输入到Scoops3.D 稳定性评价中，结果显示也具有较好效果2叨 模型使用 Scoops3D通过圆弧检索法定义滑动面，模型 1.2 Scoops.3D模型 根据检索参数定义若干检索点，以每个检索点为 Scoops3:D模型由美国地质勘探局开发，是基 球心绘制半径递增的球体，球体与数字高程模型 于数字高程模型，通过三维简化毕肖普法计算各 相交，相交面定义为滑动面，如图1所示 ZA Horizontal search increment ▣Gid Dem range 。Search point Slip surface Radius increment 家 图1 Scoops3D检素原理图 Fig.I Retrieval principle of Scoops3D 在该模型的弯矩平衡法中，安全系数可以表 示为抗剪强度s与剪应力τ之比，安全系数越大则 Slide azimuth 滑块越稳定，安全系数计算如下式 Center of search sphere (1) Center of geometric element ■Slide surface 土体的抗剪强度s通过摩尔库伦准则计算： →Slide direction Slide surface direction s=c+(On-u)tand (2) W:Slide mass gravity 式中，c为土的黏聚力，中为土的内摩擦角，0m为滑 N:Normal force 动体受到的正应力，！为切面的孔隙水压力. S:Shear force Scoops3.D稳定性模型对地震荷载是通过定义 Ko W:Horizontal seismic load R:Search radius 水平地震加速度系数Kg实现的，被视为作用在滑 动体上WK的水平力.根据单元栅格内滑动体受 Horizontal plane 力矩平衡（图2），重力爪、地震荷载WKa、切向力 图2单元栅格示意图 S=tA对于球心产生的力矩和为0，e为球心到滑体 Fig.2 Schematic diagram of cell grid 重心的垂直距离，其余参数已在示意图中标出： 带入式(2)和式(1)整理后得到三维稳定性系数计 Ti.jAi.jRi.j-Wi.jRi.jsinQi.j-Wi.jKeqei.j=0 (3) 算公式： 根据垂直方向的力平衡，重力W、法向力N=oA Fs= 的垂向分力、切向力S的垂向分力的合力为0： Ri.jci.jAi.jcossi.j+Wi.j-ui.jAi.jcossi.j)tani.j/maij i.jAi.jcossi.j+Ti.jAi.jsinai.j-Wi.j=0 (4) ∑Wi.jRi.jsin@i.j+Keqei 通过式(3)和式(4)求解剪应力x和正应力o, (5)

坡模型，在计算过程中将基底面以上的土体分为 非饱和区和毛细饱和区，并通过 Richards 方程计算 降雨过程中压力水头空间分布. Alvioli 等[29] 开发 的 TRIGRS 2.1 提供了输出空间分布的孔隙水压力 和体积含水量的功能，并可直接输入到 Scoops3D 模型使用. 1.2    Scoops3D 模型 Scoops3D 模型由美国地质勘探局开发，是基 于数字高程模型，通过三维简化毕肖普法计算各 滑动面上的栅格安全系数[26] . 该模型最初是为了 评估火山灰地层中的滑坡稳定性，在尼加拉瓜的 卡西塔和圣克里斯托瓦尔火山等多个地点进行了 成功测试[26] . 随后还应用到华盛顿等区域的滑坡 稳定性评价中，结果显示也具有较好效果[27] . Scoops3D 通过圆弧检索法定义滑动面，模型 根据检索参数定义若干检索点，以每个检索点为 球心绘制半径递增的球体，球体与数字高程模型 相交，相交面定义为滑动面,如图 1 所示. Y X Z X Horizontal search increment Radius increment Vertical search increment Grid Dem range Search point Slip surface 图 1 Scoops3D 检索原理图 Fig.1 Retrieval principle of Scoops3D 在该模型的弯矩平衡法中，安全系数可以表 示为抗剪强度 s 与剪应力 τ 之比，安全系数越大则 滑块越稳定，安全系数计算如下式 Fs = s τ （1） 土体的抗剪强度 s 通过摩尔库伦准则计算： s = c+(σn −u)tanϕ （2） 式中，c 为土的黏聚力，ϕ 为土的内摩擦角，σn 为滑 动体受到的正应力，u 为切面的孔隙水压力. Scoops3D 稳定性模型对地震荷载是通过定义 水平地震加速度系数 Keq 实现的，被视为作用在滑 动体上 WKeq 的水平力. 根据单元栅格内滑动体受 力矩平衡（图 2），重力 W、地震荷载 WKeq、切向力 S=τA 对于球心产生的力矩和为 0，e 为球心到滑体 重心的垂直距离，其余参数已在示意图中标出： τi, jAi, jRi, j − Wi, jRi, jsinαi, j − Wi, jKeqei, j = 0 （3） 根据垂直方向的力平衡，重力 W、法向力 N=σnA 的垂向分力、切向力 S 的垂向分力的合力为 0： σi, jAi, j cosεi, j +τi, jAi, jsinαi, j − Wi, j = 0 （4） 通过式（3）和式（4）求解剪应力 τ 和正应力 σn， 带入式（2）和式（1）整理后得到三维稳定性系数计 算公式： Fs = ∑ Ri, j [ ci, jAi, j cosεi, j+ ( Wi, j −ui, jAi, j cosεi, j ) tanϕi, j ] /mai, j ∑ Wi, j [ Ri, jsinαi, j+Keqei, j ] （5） Center of search sphere Center of geometric element Slide surface Slide direction Slide surface direction W: Slide mass gravity N: Normal force S: Shear force KeqW: Horizontal seismic load R: Search radius Horizontal plane Slide azimuth x y z N KeqW W R S α ε 图 2 单元栅格示意图 Fig.2 Schematic diagram of cell grid · 442 · 工程科学学报，第 44 卷，第 3 期

李泊良等：降雨和地震条件下浅层黄土滑坡三维稳定性评价 443. maij =cossi.j+sinai.jtandi.j/Fs (6) 293.3mm,降雨主要集中在7~9月.研究区出露 模型输人参数敏感性分析可直接地反映各 地层主要为第四系马兰黄土，研究区内新构造运 参数对模型影响程度的大小，从而确定每个参数 动较为强烈，以垂直升降运动为主，具有明显的继 的重要性0.本文基于稳定性系数计算公式，开展 承性、差异性特点.基于1：5000地形图，通过详 参数的敏感性分析，目的是分析各输入参数对 细的现场调查和GPS现场测量，绘制了61处浅层 Scoops3D模型影响程度程度，更重要的是指导在 滑坡点，如图4.根据王恭先滑坡滑动面分类标准 试验测试基础上确定输入参数的取值.根据研究 B别，这些滑坡厚度均小于10m,属于浅层黄土滑 区的实际地层分布和试验测试的结果，确定了输 坡.据历史地震资料记载，兰州市历史上曾多次遭 入参数的标准值和取值范围（表1），保持其他参数 到地震破坏，地震基本烈度为Ⅷ度，地震动峰值加 为标准值，在取值范围内改变待分析参数，得到系 速度系数为0.2B2 数变化率，得到F对输入变量的相对敏感性变化 规律（图3）. 表1敏感性分析参数变化范围 Table 1 Variation range of sensitivity analysis parameters Item R/m C/kPa WIN /Kea ukPa Range3-1510-5010-2.0×10100-2.0×1010-300-0.2-50-100 Standard 9 30 1.0x×104 1.1×105 200.1 25 100 60 Pore water pr 40 20 Slope,a 0 -20 Landslide points -40 2023m -60 -80 1575m 100-80-60 -40-20020406080.100 Parameter change rate/% 图4研究区的数字高程模型和滑坡分布图 图3参数敏感性分析结果 Fig.4 Study area digital elevation model and landslide distribution map Fig.3 Sensitivity analysis results of parameter 3材料和方法 安全系数随黏聚力、内摩擦角、滑动面面积、 球体半径、坡度的增大而增大，随地震加速度、孔 3.1数据和参数 隙水压力、滑块重量的增大而减小，这表明地形地 Scoops3D模型和TRIGRS需要研究区的强度 貌和土体性质共同控制着斜坡稳定性（图3）.在土 参数，本文参考了wen和YanB对马兰黄土强度 体性质和地形条件基本确定的情况下，孔隙水压 的研究，两个模型所需的黄土参数设置如下表2 力和地震加速度在更大程度的控制斜坡的稳定 性.而降雨和地震也是触发滑坡最常见的要素.因 表2土参数取值表 此，下文重点分析了研究区降雨、地震及两者耦合 Table 2 Soil parameter value table Permeability coefficient,Hydraulic diffusivity, 条件下斜坡的稳定性 kPa (kN'm)( k/(m's) ha /(m3s 1) 30 15 29 2.4×10t 2.4×104 2研究区概况 选择甘肃省兰州市大沙沟流域作为研究区 降雨增加土体含水量，从而降低土体强度，斜 它处于黄土高原，地貌单元复杂，坡体坡度在0°到 坡稳定性随之降低.土的抗剪强度在斜坡中起着 61°之间，海拔在1575m到2023m之间，区域面积 重要作用，而非饱和黄土的基质吸力对抗剪强度 100km2.研究区属于半干旱气候，平均年降水为 有重要影响，Scoops.3D中通过含水率与土水特征

mai j = cosεi, j +sinαi, j tanϕi, j/Fs （6） 模型输入参数敏感性分析可直接地反映各 参数对模型影响程度的大小，从而确定每个参数 的重要性[30] . 本文基于稳定性系数计算公式，开展 参数的敏感性分析 ，目的是分析各输入参数对 Scoops3D 模型影响程度程度，更重要的是指导在 试验测试基础上确定输入参数的取值. 根据研究 区的实际地层分布和试验测试的结果，确定了输 入参数的标准值和取值范围（表 1），保持其他参数 为标准值，在取值范围内改变待分析参数，得到系 数变化率，得到 Fs 对输入变量的相对敏感性变化 规律（图 3）. 表 1 敏感性分析参数变化范围 Table 1 Variation range of sensitivity analysis parameters Item R /m C/kPa A/m2 W/N ϕ /(°) Keq u/kPa Range 3–15 10–50 10–2.0×104 100–2.0×106 10–30 0–0.2 −50–100 Standard 9 30 1.0×104 1.1×106 20 0.1 25 100 80 60 40 20 0 −20 −40 −60 −80 −100 −80 −60 −40 −20 0 Parameter change rate/% 20 40 60 80 100 Safety factor change rate/ % Cohesion, C Slide surface area, A Internal friction angle, ϕ Sphere radius, R Pore water pressure, U Earthquake acceleration, Keq Slide mass gravity, W Slope, α 图 3 参数敏感性分析结果 Fig.3 Sensitivity analysis results of parameter 安全系数随黏聚力、内摩擦角、滑动面面积、 球体半径、坡度的增大而增大，随地震加速度、孔 隙水压力、滑块重量的增大而减小，这表明地形地 貌和土体性质共同控制着斜坡稳定性（图 3）. 在土 体性质和地形条件基本确定的情况下，孔隙水压 力和地震加速度在更大程度的控制斜坡的稳定 性. 而降雨和地震也是触发滑坡最常见的要素. 因 此，下文重点分析了研究区降雨、地震及两者耦合 条件下斜坡的稳定性. 2    研究区概况 选择甘肃省兰州市大沙沟流域作为研究区. 它处于黄土高原，地貌单元复杂，坡体坡度在 0°到 61°之间，海拔在 1575 m 到 2023 m 之间，区域面积 100 km2 . 研究区属于半干旱气候，平均年降水为 293.3 mm，降雨主要集中在 7～9 月. 研究区出露 地层主要为第四系马兰黄土. 研究区内新构造运 动较为强烈，以垂直升降运动为主，具有明显的继 承性、差异性特点. 基于 1∶5000 地形图，通过详 细的现场调查和 GPS 现场测量，绘制了 61 处浅层 滑坡点，如图 4. 根据王恭先滑坡滑动面分类标准 [31] ，这些滑坡厚度均小于 10 m，属于浅层黄土滑 坡. 据历史地震资料记载，兰州市历史上曾多次遭 到地震破坏，地震基本烈度为Ⅷ度，地震动峰值加 速度系数为 0.2 [32] . 2023 m Landslide points 1575 m 0 1 2 km N 图 4 研究区的数字高程模型和滑坡分布图 Fig.4 Study area digital elevation model and landslide distribution map 3    材料和方法 3.1    数据和参数 Scoops3D 模型和 TRIGRS 需要研究区的强度 参数，本文参考了 Wen 和 Yan[33] 对马兰黄土强度 的研究，两个模型所需的黄土参数设置如下表 2. 表 2 土参数取值表 Table 2 Soil parameter value table C / kPa γ / (kN·m−3) ϕ/ (°) Permeability coefficient, k /(m·s−1) Hydraulic diffusivity, hd /(m3 ·s−1) 30 15 29 2.4×10−6 2.4×10−4 降雨增加土体含水量，从而降低土体强度，斜 坡稳定性随之降低. 土的抗剪强度在斜坡中起着 重要作用，而非饱和黄土的基质吸力对抗剪强度 有重要影响，Scoops3D 中通过含水率与土水特征 李泊良等： 降雨和地震条件下浅层黄土滑坡三维稳定性评价 · 443 ·