融合InSAR形变量与降雨量的滑坡危险性评价方法

本发明涉及一种滑坡危险性评价方法，具体涉及一种融合insar形变量与降雨量的滑坡危险性评价方法。

背景技术：

1、滑坡危险性评价是滑坡灾害风险防控的前提和基础，国内外对滑坡危险性评价开展了大量研究，在评价指标因子选取上，多为地形地貌、地层岩性、地质构造、人类工程活动，对于地质结构与地下水在地球物理场的响应特征、岩层与坡向组合关系和矿山开采情况未见纳入评价指标。在易发性评价方法上，多以指标权重计算，概括起来有定性方法、数理统计法、机器学习法，其常用的为层次分析法、专家打分法、确定性系数、信息量法、证据权法、神经网络、逻辑回归、随机森林等，其中定性方法评价精度较低，且个人主观性较高，不同的专家评价结果不一。在危险性评价方法上，对于滑坡灾害多以静态因子评价，比如降雨量以年均降雨量为评价因子，将insar形变量作为危险性评价因素，已有研究虽然有考虑insar形变量，但是也是将其考虑为静态评价因子，比如形变速率。综上所述，已有研究在评价指标因子选取上未考虑地质结构与地下水在地球物理场的响应特征、岩层与坡向组合关系和矿山开采情况，危险性评价以静态评价为主。

技术实现思路

1、为了解决现有技术中指标因子不全面、评价方法单一的问题，本发明提供了一种融合insar形变量与降雨量的滑坡危险性评价方法。该方法分别考虑了insar的形变量以及有效降雨量，基于滑坡易发性，采用矩阵分析法将insar地表形变量的滑坡危险性评价结果与降雨量阈值的滑坡危险性评价结果叠加运算进行等级划分，使得滑坡危险性评价更加合理，实现随降雨量变化、地表变形的滑坡危险性动态评价。

2、本发明的目的是通过以下技术方案实现的：

3、一种融合insar形变量与降雨量的滑坡危险性评价方法，包括如下步骤：

4、步骤s1：多源异构数据获取：

5、获取多源异构数据，作为滑坡危险性评价的评价指标，其中：多源异构数据包括：(1)sar影像，(2)地形地貌，(3)地球物理场特征，(4)基础地质，(5)水文环境，(6)人类工程活动，其中：地球物理场特征选择电阻率数据来表征；地形地貌数据包括高程、坡度、坡向、起伏度、地形曲率、剖面曲率；基础地质包括地层岩性、断层破碎带、岩层与坡向组合关系(逆向坡、顺向坡、切向坡)；水文环境包括水系河流、归一化植被指数、土地利用类型；人类工程活动包括道路切坡、矿山开采密度；

6、步骤s2：滑坡易发性评价：构建滑坡易发性评价模型，基于步骤s1获取的评价指标与历史滑坡数据进行滑坡易发性评价；

7、步骤s3：insar技术地表形变量计算与分级：基于步骤s1获取的sar影像数据，运用insar技术计算形变量，计算得到每个历史滑坡发生前后形变量，进行数理统计分析，按滑坡概率在某insar形变量区间内的概率进行等级划分；

8、步骤s4：有效降雨量计算与分级：计算研究区有效降雨量，结合历史滑坡数据发生时间，统计分析有效降雨量达到阈值时，对应滑坡发生概率，按概率值分布划定有效降雨量区间；

9、步骤s5：基于insar地表形变量的滑坡危险性评价：采用矩阵分析法将步骤s2的滑坡易发性评价结果与步骤s3计算的insar形变量叠加运算进行等级划分，得到基于insar地表形变量的滑坡危险性评价结果；

10、步骤s6：基于降雨量阈值的滑坡危险性评价：采用矩阵分析法将步骤s2的滑坡易发性评价结果与步骤s4的降雨量阈值计算结果叠加运算进行等级划分，得到基于降雨量阈值的滑坡危险性评价结果；

11、步骤s7：滑坡危险性综合评价：采用矩阵分析法将步骤s5基于insar地表形变量的滑坡危险性评价结果与步骤s6基于降雨量阈值的滑坡危险性评价结果叠加运算进行等级划分，得到滑坡危险性综合评价结果。

12、相比于现有技术，本发明具有如下优点：

13、1、将岩层与坡向组合关系(逆向坡、顺向坡、切向坡)作为评价因子更加全面的考虑了滑坡的孕灾环境；

14、2、斜坡结构以及地下水是滑坡孕灾环境的重要因素，然而在已有方法中难以在区域性滑坡危险性评价中量化评价，而斜坡结构以及地下水的变化在地球物理响应特征上存在差异性，将区域电阻率平面等值线图作为评价因子也是更加全面的考虑了滑坡的孕灾环境；

15、3、矿山开采存在大量削坡切坡、地下空洞等，导致边坡失稳、采空区坍塌易引发滑坡，矿山周围存在大量人类工程活动，尤其是重型车托运，对周围斜坡增加了动荷载，因而将矿山开采密度作为评价因子，使评价因子更加全面；

16、4、本发明的方法分别考虑了insar的形变量以及有效降雨量，基于滑坡易发性，采用矩阵分析法将insar地表形变量的滑坡危险性评价结果与有效降雨量的滑坡危险性评价结果叠加运算进行等级划分，使得滑坡危险性评价更加合理。

技术特征：

1.一种融合insar形变量与降雨量的滑坡危险性评价方法，其特征在于所述方法包括如下步骤：

2.根据权利要求1所述的融合insar形变量与降雨量的滑坡危险性评价方法，其特征在于所述步骤s1中，地球物理场特征选择电阻率数据来表征；地形地貌数据包括高程、坡度、坡向、起伏度、地形曲率、剖面曲率；基础地质包括地层岩性、断层破碎带、岩层与坡向组合关系；水文环境包括水系河流、归一化植被指数、土地利用类型；人类工程活动包括道路切坡、矿山开采密度。

3.根据权利要求1所述的融合insar形变量与降雨量的滑坡危险性评价方法，其特征在于所述s1的具体步骤如下：

4.根据权利要求1所述的融合insar形变量与降雨量的滑坡危险性评价方法，其特征在于所述步骤s2的具体步骤如下：

5.根据权利要求1所述的融合insar形变量与降雨量的滑坡危险性评价方法，其特征在于所述步骤s3的具体步骤如下：

6.根据权利要求1所述的融合insar形变量与降雨量的滑坡危险性评价方法，其特征在于所述步骤s4的具体步骤如下：

技术总结
本发明公开了一种融合InSAR形变量与降雨量的滑坡危险性评价方法，所述方法包括如下步骤：步骤S1：多源异构数据获取；步骤S2：滑坡易发性评价；步骤S3：InSAR技术地表形变量计算与分级；步骤S4：有效降雨量计算与分级；步骤S5：基于InSAR地表形变量的滑坡危险性评价；步骤S6：基于降雨量阈值的滑坡危险性评价；步骤S7：滑坡危险性综合评价。本发明解决了现有技术滑坡灾害危险性评价方法的不足，同时本发明实施简单方便，结果呈现在GIS里面以图像形式呈现，根据降雨量变化，形变量变化，可进一步开发自动计算软件，评估结果能为当地应急管理部门提供决策指导意义，为防灾减灾相关部门提供防治科学依据。

技术研发人员：丛晟亦,凌贤长,蔡德钩,闫宏业,张钟远,唐亮,毛小刚,孙辉,赵香萍,田爽,孔祥勋
受保护的技术使用者：哈尔滨工业大学
技术研发日：
技术公布日：2024/1/15