(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210146901.6 (22)申请日 2022.02.17 (71)申请人 成都理工大 学 地址 610059 四川省成 都市成华区二仙桥 东三路一号 (72)发明人 曾鹏　王宇豪　李天斌　孙小平　 张琳　 (74)专利代理 机构 北京润泽恒知识产权代理有 限公司 1 1319 专利代理师 苟冬梅 (51)Int.Cl. G06F 30/23(2020.01) G06F 30/27(2020.01) G06N 3/12(2006.01) G06T 13/60(2011.01)G06T 17/05(2011.01) G06T 17/20(2006.01) G06V 10/762(2022.01) G06K 9/62(2022.01) G06F 111/04(2020.01) G06F 111/06(2020.01) G06F 111/10(2020.01) G06F 119/14(2020.01) (54)发明名称 一种黄土滑坡稳定性预测方法、 电子 设备以 及存储介质 (57)摘要 本申请提供一种黄土滑坡 稳定性预测方法、 电子设备以及存储介质， 属于稳定性预测的领 域。 所述方法包括： 获取目标地区的抗剪强度参 数， 并构建抗剪强度参数均值误差函数； 根据目 标滑坡第一次失稳的滑面观测信息， 以及通过有 限差分数值模拟计算确定的临界滑动面， 构建滑 面位置误差函数； 通过有限差分强度折减法计算 得到滑坡 稳定性系数， 构建滑坡 稳定性系数误差 函数； 结合构建的三个误差函数， 基于多组随机 抽样样本利用遗传算法进行多次迭代， 得到最优 参数组； 根据目标滑坡的其他参数与最优参数 组， 构建模拟滑动面， 在所述一致性表征正确的 情况下， 预测目标滑坡的稳定性预测结果。 本申 请旨在获得更可靠的反分析结果， 并更准确地预 测滑坡的稳定性。 权利要求书2页 说明书11页 附图6页 CN 114692441 A 2022.07.01 CN 114692441 A 1.一种黄 土滑坡稳定性预测方法， 其特 征在于， 所述方法包括： 获取目标地区的抗剪强度参数， 所述抗剪强度参数包括多组天然黄土的有效黏聚力、 多组天然黄土的有效内摩擦角、 多组饱和黄土的有效黏聚力以及多组饱和黄土的有效内摩 擦角的数据， 根据所述 抗剪强度参数构建抗剪强度参数均值 误差函数； 根据目标滑坡的第 一次失稳的滑面观测信 息， 以及通过有限差分数值模拟计算确定的 临界滑动面， 构建滑面 位置误差函数； 对所述抗剪强度参数进行抽样， 根据抽样得到的黏聚力数据与内摩擦角数据， 通过有 限差分强度折减法计算得到滑坡稳定性系数， 构建滑坡稳定性系数误差函数； 结合所述抗剪强度参数均值误差函数、 所述滑面位置误差函数以及所述滑坡稳定性系 数误差函数， 基于多组随机抽样样本利用遗传算法进行多次迭代， 得到最优参数组， 其中， 最优参数组包括最优的天然黄土的有效黏聚力、 天然黄土的有效内摩擦角、 饱和黄土的有 效黏聚力以及饱和黄 土的有效内摩擦角； 根据所述目标滑坡的其他参数与所述最优参数组， 在Flac模型中构建模拟滑动面， 对 比所述模拟滑动面与所述目标滑坡的第一失稳的滑面观测信息的一 致性； 在所述一致性表征正确的情况下， 通过Flac模型得到所述目标滑坡的稳定性预测结 果， 所述预测结果包括后续失稳的稳定性系数与滑面信息 。 2.根据权利要求1所述的方法， 其特征在于， 根据 所述抗剪强度参数构建抗剪强度参数 均值误差函数， 包括： 分别计算多组天然黄土的有效黏聚力、 多组天然黄土的有效内摩擦角、 多组饱和黄土 的有效黏聚力以及多组饱和黄土的有效内摩擦角各自的均值与标准差， 并假设均值与标准 差服从正态分布； 根据计算得到的均值与标准差， 构建抗剪强度参数均值 误差函数。 3.根据权利要求2所述的方法， 其特征在于， 构建的所述抗剪强度参数均值误差函数 为： 式中， x为抗剪强度参数向量， μx为各参数的均值， 表示各参数的协方差矩阵的逆矩 阵， T表示矩阵转置 。 4.根据权利要求1所述的方法， 其特征在于， 根据目标滑坡第一次失稳的滑面观测信 息， 以及通过有限差分数值模拟计算确定的临界滑动面， 构建滑面 位置误差函数， 包括： 根据目标滑坡第 一次失稳的滑面观测信 息， 选取第 一次失稳的滑面上的多个特征坐标 点； 通过多次有限差分数值模拟计算确定多个临界滑动 面， 分别在每个所述临界滑面上获 取多个特 征坐标点； 根据所述第一次失稳的滑面上的多个特征坐标点与每个所述临界滑面上的多个特征 坐标点， 构建滑面 位置误差函数。 5.根据权利要求1所述的方法， 其特征在于， 根据目标滑坡第一次失稳的滑面观测信 息， 以及通过有限差分数值模拟计算确定的临界滑动面， 构建滑面 位置误差函数， 包括： 根据目标滑坡第一次失稳的滑面观测信息， 选取第一次失稳的滑面上的5个特征坐标权　利　要　求　书 1/2 页 2 CN 114692441 A 2点， 记为(ai， bi)， 其中， i ＝1， 2， 3， 4， 5； 通过多次有限差分数值模拟计算确定多个临界滑动面， 获取每个临界滑面上的5个特 征坐标点， 记为(xi， yi)， 其中， i ＝1， 2， 3， 4， 5； 根据第一次失稳的滑面上的5个特征坐标点与临界滑面上的5个特征坐标点， 构建的滑 面位置误差函数为： 6.根据权利要求1所述的方法， 其特征在于， 对所述抗剪强度参数进行抽样， 根据抽样 得到的黏聚力数据与内摩擦角数据， 通过有限差分强度折减法计算得到滑坡稳定性系 数， 构建滑坡稳定性系数误差函数， 包括： 对所述抗剪强度参数进行抽样， 根据抽样得到的黏聚力数据与内摩擦角数据， 通过有 限差分强度折减法计算得到滑坡稳定性系数； 假设所述目标滑坡失稳时的稳定性系数为1； 构建的滑坡稳定性系数误差函数为： f3(x)＝(FS(x) ‑1) 式中， FS为有限差分强度折减法计算得到滑坡稳定性系数。 7.一种电子设备， 包括存储器、 处理器及存储在所述存储器上并可在所述处理器上运 行的计算机程序， 其特征在于， 所述处理器执行所述计算机程序时实现权利要求1至6中任 一项所述的方法。 8.一种计算机可读存储介质， 其特征在于， 所述计算机可读存储介质上存储计算机程 序， 所述计算机程序被处 理器执行时实现如权利要求1至 6中任一项所述的方法。权　利　要　求　书 2/2 页 3 CN 114692441 A 3