(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202211036619.9 (22)申请日 2022.08.29 (65)同一申请的已公布的文献号 申请公布号 CN 115114807 A (43)申请公布日 2022.09.27 (73)专利权人 成都理工大 学 地址 610059 四川省成 都市成华区二仙桥 东三路1号 (72)发明人 史先琳　戴可人　陈晨　 (74)专利代理 机构 成都正德明志知识产权代理 有限公司 513 60 专利代理师 陈瑶 (51)Int.Cl. G06F 30/20(2020.01) G06F 30/23(2020.01) G06Q 10/06(2012.01) G06T 7/30(2017.01) G06T 17/20(2006.01) G06F 119/14(2020.01)(56)对比文件 CN 113866764 A,2021.12.31 CN 102608584 A,2012.07.25 WO 2021217 763 A1,2021.1 1.04 张路等.基 于时间序列I nSAR分析的西部山 区滑坡灾害隐患早期识别― ―以四川丹巴为例. 《武汉大 学学报(信息科 学版)》 .2018,第43卷(第 12期),20 39-2049. 李萌等.基 于地表变形 数据的潜在滑坡识别 研究. 《科技创新与应用》 .2020,(第0 6期),9-13. 史先琳等.高山峡谷区滑坡灾害隐患I nSAR 早期识别— —以雅砻江中段为例. 《雷达学报》 .2020,第9卷(第3期),5 54-568. Keren Dai 等.Interpretati on and sensitivity analysis of the I nSAR line of sight displacements i n landsl ide measurements. 《GIScience & Remote Sensing》 .2022,第59卷(第1期),12 26-1242. 审查员 田娜 (54)发明名称 一种水库 库岸滑坡易发性评价方法 (57)摘要 本发明公开了一种水库库岸滑坡易发性评 价方法， 包括以下步骤： S1、 采集地质物理参数； S2、 根据地质物理参数， 通过静态地质评价模型， 计算水库蓄水后的库岸稳定性； S3、 根据水库蓄 水后的SLC影像， 得到地理坐标系下的平均形变 速率； S4、 根据地理坐标系下的平均形变速率以 及水库蓄水后的库岸稳定性， 得到水库库岸滑坡 易发程度； 本发 明解决了 现有滑坡易发性评价方 法， 均采用单一模型对滑坡易发性进行评价， 存 在对滑坡易发性评价 准确度不高的问题。 权利要求书2页 说明书5页 附图1页 CN 115114807 B 2022.12.20 CN 115114807 B 1.一种水库 库岸滑坡易发性评价方法， 其特 征在于， 包括以下步骤： S1、 采集地质物理参数； S2、 根据地质物理参数， 通过静态地质评价模型， 计算水库蓄水后的库岸稳定性； 步骤S2包括以下分步骤： S21、 在水库地形 上方建立检索矩阵； S22、 在检索矩阵内生成若干个 球面检索点； S23、 根据球面检索点 生成球面检索点的球心； S24、 将潜在滑坡体网格化 为若干个三维柱体； S25、 根据球面检索点的球心和三维柱体， 以及地质物理参数， 构建静态地质评价模型； 步骤S25中静态地质评价模型的公式为： 其中， Fs为水库蓄水后的库岸稳定性， R为球心到潜在滑坡面的距离， A为库岸潜在滑坡 面表面积， α 为库岸潜在滑坡面的视倾角， β 为库岸潜在滑坡面的倾角， e为球心到三维柱体 质心的距离， W为岩土体的容重， u为岩土体的孔隙水压力， k为水平振动荷载， c为岩土体的 黏聚力， 为岩土体的内摩擦角， T为比例常数； S26、 根据静态地质评价模型， 计算水库蓄水后的库岸稳定性； S3、 根据水库蓄水后的SLC影 像， 得到地理坐标系下的平均形变速率； 步骤S3包括以下分步骤： S31、 根据水库蓄水后的SLC影 像， 生成M对干涉对， 计算 LOS向形变相位； S32、 根据LOS向形变相位， 进行两次相位 解缠和滤波， 得到去误差 干涉对； S33、 根据去误差 干涉对， 反演出干涉对的平均形变速率； S34、 对干涉对的平均形变速率进行地理编码， 得到地理坐标系下的平均形变速率； S4、 根据地理坐标系下的平均形变速率， 以及水库蓄水后的库岸稳定性， 得到水库库岸 滑坡易发程度； 所述步骤S4中计算水库 库岸滑坡易发程度的公式为： Wi＝2.5(i‑3)2+3 Wj＝2.5(j‑3)2+3 其中， X为水库库岸滑坡易发程度， i,j为大于0的整数， Vi为平均形变速率所在的第i个 等级， Wi为平均形变速率所在的第i个等级Vi的权重， Sj为水库蓄水后的库岸稳定性所在的 第j个等级， Wj为水库蓄水后的库岸稳定性所在的第j个等级Sj的权重。 2.根据权利要求1所述的水库库岸滑坡易发性评价方法， 其特征在于， 所述步骤S1中地 质物理参数包括： 水库蓄水后， 水库库岸附近岩土体的黏聚力、 内摩擦角、 容重和孔隙水压 力。 3.根据权利要求1所述的水库库岸滑坡易发性评价方法， 其特征在于， 所述步骤S31包 括以下分步骤：权　利　要　求　书 1/2 页 2 CN 115114807 B 2S311、 采集水库蓄水后多张SLC影 像； S312、 选择一张SLC影像作为主影像， 采用主影像与剩余SLC影像进行配准， 生成M对干 涉对； S313、 对每 个干涉对进行差分处 理， 得到M幅多视差分干涉图； S314、 根据M幅多视 差分干涉图的干涉相位， 计算 LOS向形变相位。 4.根据权利要求3所述的水库库岸滑坡易发性评价方法， 其特征在于， 所述步骤S312中 M对干涉对的范围为： N/2≤M≤(N‑1)N/2 其中， N为SLC影像的数量。 5.根据权利要求3所述的水库库岸滑坡易发性评价方法， 其特征在于， 所述步骤S314中 计算LOS向形变相位的公式为： 其中， 为LOS向形变相位， 为干涉图的干涉相位， 为地形相位， 为大 气相位， 为噪声相位。 6.根据权利要求1所述的水库库岸滑坡易发性评价方法， 其特征在于， 所述步骤S32包 括以下分步骤： S321、 采用最小费用流法对LOS向形变相位进行相位 解缠， 得到 干涉解缠相位 集； S322、 修正干涉解缠相位集的轨道参数， 并去除干涉解缠相位集中的残余相位和相位 跃变， 得到解缠修 正干涉对； S323、 对解缠修 正干涉对中未解缠相位进行 再次解缠， 得到二次解缠干涉对； S324、 对二次解缠干涉对进行在空间上的高通滤波和在时间上的低通滤波， 得到去误 差干涉对。权　利　要　求　书 2/2 页 3 CN 115114807 B 3