一种水库库岸滑坡易发性评价方法技术

本发明专利技术公开了一种水库库岸滑坡易发性评价方法，包括以下步骤：S1、采集地质物理参数；S2、根据地质物理参数，通过静态地质评价模型，计算水库蓄水后的库岸稳定性；S3、根据水库蓄水后的SLC影像，得到地理坐标系下的平均形变速率；S4、根据地理坐标系下的平均形变速率以及水库蓄水后的库岸稳定性，得到水库库岸滑坡易发程度；本发明专利技术解决了现有滑坡易发性评价方法，均采用单一模型对滑坡易发性进行评价，存在对滑坡易发性评价准确度不高的问题。在对滑坡易发性评价准确度不高的问题。在对滑坡易发性评价准确度不高的问题。

【技术实现步骤摘要】
一种水库库岸滑坡易发性评价方法


[0001]本专利技术涉及滑坡评估
，具体涉及一种水库库岸滑坡易发性评价方法。

技术介绍

[0002]库水位变动导致斜坡岩土体灾害突出，尤其滑坡灾害成为库区广泛存在的地质问题是库区经济和社会发展的最大威胁和阻碍。水库蓄水、放水可诱发大规模库岸滑坡，易发性评价可利用多种类评价模型结合各种评价因子掌控滑坡可能发生位置，因此及时识别库岸潜在滑坡隐患位置对水电站安全运行具有重要意义。
[0003]目前，关于滑坡易发性已提出很多评价方法，但是其均采用单一模型对滑坡易发性进行评价，存在对滑坡易发性评价准确度不高的问题。

技术实现思路

[0004]针对现有技术中的上述不足，本专利技术提供的一种水库库岸滑坡易发性评价方法解决了现有滑坡易发性评价方法，均采用单一模型对滑坡易发性进行评价，存在对滑坡易发性评价准确度不高的问题。
[0005]为了达到上述专利技术目的，本专利技术采用的技术方案为：一种水库库岸滑坡易发性评价方法，包括以下步骤：S1、采集地质物理参数；S2、根据地质物理参数，通过静态地质评价模型，计算水库蓄水后的库岸稳定性；S3、根据水库蓄水后的SLC影像，得到地理坐标系下的平均形变速率；S4、根据地理坐标系下的平均形变速率，以及水库蓄水后的库岸稳定性，得到水库库岸滑坡易发程度。
[0006]进一步地，所述步骤S1中地质物理参数包括：水库蓄水后，水库库岸附近岩土体的黏聚力、内摩擦角、容重和孔隙水压力。
[0007]进一步地，所述步骤S2包括以下分步骤：S21、在水库地形上方建立检索矩阵；S22、在检索矩阵内生成若干个球面检索点；S23、根据球面检索点生成球面检索点的球心；S24、将潜在滑坡体网格化为若干个三维柱体；S25、根据球面检索点的球心和三维柱体，以及地质物理参数，构建静态地质评价模型；S26、根据静态地质评价模型，计算水库蓄水后的库岸稳定性。
[0008]进一步地，所述步骤S25中静态地质评价模型的公式为：
其中，为水库蓄水后的库岸稳定性，为球心到潜在滑坡面的距离，为库岸潜在滑坡面表面积，为库岸潜在滑坡面的视倾角，为库岸潜在滑坡面的倾角，为球心到三维柱体质心的距离，为岩土体的容重，为岩土体的孔隙水压力，为水平振动荷载，为岩土体的黏聚力，为岩土体的内摩擦角，为比例常数。
[0009]进一步地，所述步骤S3包括以下分步骤：S31、根据水库蓄水后的SLC影像，生成M对干涉对，计算LOS向形变相位；S32、根据LOS向形变相位，进行两次相位解缠和滤波，得到去误差干涉对；S33、根据去误差干涉对，反演出干涉对的平均形变速率；S34、对干涉对的平均形变速率进行地理编码，得到地理坐标系下的平均形变速率。
[0010]进一步地，所述步骤S31包括以下分步骤：S311、采集水库蓄水后多张SLC影像；S312、选择一张SLC影像作为主影像，采用主影像与剩余SLC影像进行配准，生成M对干涉对；S313、对每个干涉对进行差分处理，得到M幅多视差分干涉图；S314、根据M幅多视差分干涉图的干涉相位，计算LOS向形变相位。
[0011]进一步地，所述步骤S312中M对干涉对的范围为：其中，为SLC影像的数量。
[0012]上述进一步地方案的有益效果为：将水库蓄水后多张SLC影像进行干涉处理，获取水库蓄水后的干涉相位图以便于查看较为明显的形变点。
[0013]进一步地，所述步骤S314中计算LOS向形变相位的公式为：其中，为LOS向形变相位，为干涉图的干涉相位，为地形相位，为大气相位，为噪声相位。
[0014]进一步地，所述步骤S32包括以下分步骤：S321、采用最小费用流法对LOS向形变相位进行相位解缠，得到干涉解缠相位集；S322、修正干涉解缠相位集的轨道参数，并去除干涉解缠相位集中的残余相位和相位跃变，得到解缠修正干涉对；S323、对解缠修正干涉对中未解缠相位进行再次解缠，得到二次解缠干涉对；S324、对二次解缠干涉对进行在空间上的高通滤波和在时间上的低通滤波，得到去误差干涉对。
[0015]上述进一步地方案的有益效果为：将干涉处理后获取的干涉相位经过最小费用流法、滤波等步骤去除地形相位、大气相位等误差相位。
[0016]综上，本专利技术的有益效果为：本专利技术一方面通过地质物理参数和静态地质评价模
型，计算水库蓄水后的库岸稳定性，另一方面通过SLC影像，生成M对干涉对，并通过两次相位解缠和滤波，实现对干涉对去除残余的相位斜坡或地形相位，从而精确得到地理坐标系下的平均形变速率，相比于现有单一的模型，本专利技术综合了地理坐标系下的平均形变速率，以及水库蓄水后的库岸稳定性，得到水库库岸滑坡易发程度，更能准确反应水库库岸滑坡的易发性情况。
附图说明
[0017]图1为一种水库库岸滑坡易发性评价方法的流程图。
具体实施方式
[0018]下面对本专利技术的具体实施方式进行描述，以便于本
的技术人员理解本专利技术，但应该清楚，本专利技术不限于具体实施方式的范围，对本
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本专利技术的精神和范围内，这些变化是显而易见的，一切利用本专利技术构思的专利技术创造均在保护之列。
[0019]如图1所示，一种水库库岸滑坡易发性评价方法，包括以下步骤：S1、采集地质物理参数；所述步骤S1中地质物理参数包括：水库蓄水后，水库库岸附近岩土体的黏聚力、内摩擦角、容重和孔隙水压力。
[0020]S2、根据地质物理参数，通过静态地质评价模型，计算水库蓄水后的库岸稳定性；所述步骤S2包括以下分步骤：S21、在水库地形上方建立检索矩阵；S22、在检索矩阵内生成若干个球面检索点；S23、根据球面检索点生成球面检索点的球心；S24、将潜在滑坡体网格化为若干个三维柱体；S25、根据球面检索点的球心和三维柱体，以及地质物理参数，构建静态地质评价模型；所述步骤S25中静态地质评价模型的公式为：其中，为水库蓄水后的库岸稳定性，为球心到潜在滑坡面的距离，为库岸潜在滑坡面表面积，为库岸潜在滑坡面的视倾角，为库岸潜在滑坡面的倾角，为球心到三维柱体质心的距离，为岩土体的容重，为岩土体的孔隙水压力，为水平振动荷载，为岩土体的黏聚力，为岩土体的内摩擦角，为比例常数。
[0021]S26、根据静态地质评价模型，计算水库蓄水后的库岸稳定性。
[0022]S3、根据水库蓄水后的SLC影像，得到地理坐标系下的平均形变速率；所述步骤S3包括以下分步骤：S31、根据水库蓄水后的SLC影像，生成M对干涉对，计算LOS向形变相位；SLC为单视复数图像（Single Look Complex），LOS为视距向（Line of Sight）。
[0023]所述步骤S31包括以下分步骤：S311、采集水库蓄水后多张SLC影像；S312、选择一张SLC影像作为主影像，采用主影像与剩余SLC影像进行配准，生成M对干涉对；所述步骤S312中M对干涉对的范围为：其中，为SLC影像的数量。
[0024]本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种水库库岸滑坡易发性评价方法，其特征在于，包括以下步骤：S1、采集地质物理参数；S2、根据地质物理参数，通过静态地质评价模型，计算水库蓄水后的库岸稳定性；S3、根据水库蓄水后的SLC影像，得到地理坐标系下的平均形变速率；S4、根据地理坐标系下的平均形变速率，以及水库蓄水后的库岸稳定性，得到水库库岸滑坡易发程度。2.根据权利要求1所述的水库库岸滑坡易发性评价方法，其特征在于，所述步骤S1中地质物理参数包括：水库蓄水后，水库库岸附近岩土体的黏聚力、内摩擦角、容重和孔隙水压力。3.根据权利要求2所述的水库库岸滑坡易发性评价方法，其特征在于，所述步骤S2包括以下分步骤：S21、在水库地形上方建立检索矩阵；S22、在检索矩阵内生成若干个球面检索点；S23、根据球面检索点生成球面检索点的球心；S24、将潜在滑坡体网格化为若干个三维柱体；S25、根据球面检索点的球心和三维柱体，以及地质物理参数，构建静态地质评价模型；S26、根据静态地质评价模型，计算水库蓄水后的库岸稳定性。4.根据权利要求3所述的水库库岸滑坡易发性评价方法，其特征在于，所述步骤S25中静态地质评价模型的公式为：其中，为水库蓄水后的库岸稳定性，为球心到潜在滑坡面的距离，为库岸潜在滑坡面表面积，为库岸潜在滑坡面的视倾角，为库岸潜在滑坡面的倾角，为球心到三维柱体质心的距离，为岩土体的容重，为岩土体的孔隙水压力，为水平振动荷载，为岩土体的黏聚力，为岩土体的内摩擦角，为比例常数。5.根据权利要求1所述的水库库岸滑坡易发性评价方法，其特征在于，所述步骤S3包括以下分步骤：S31、根据水库蓄水后的SLC影像，...

【专利技术属性】
技术研发人员：史先琳，戴可人，陈晨，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：