[发明专利]一种基于降雨的地质灾害发生频率的计算方法

权利要求书

1.一种基于降雨的地质灾害发生频率的计算方法，其特征在于，包括以下步骤：

a、通过中国气象局确定研究区境内的雨量监测站点，通过中小流域暴雨洪水计算手册中最大24小时暴雨均值等值线图确定研究区境内雨量监测点P频率降雨下的点雨量均值变差系数C_v，通过C_s＝3.5C_v适线确定研究区境内雨量监测点P频率降雨下的离均系数K_p，通过式1确定研究区境内雨量监测点P频率降雨下的年最大日降雨量H_24P，通过画图软件和克里金插值确定研究区境内雨量监测点P频率降雨下的年最大日降雨量等值线图；

b、通过地理信息系统空间分析功能，提取滑坡的内部影响因子，通过滑坡地质灾害相对密度确定滑坡的各个内部影响因子的无量纲值，通过统计产品与服务解决方案软件实现逻辑回归分析过程，确定滑坡的各个内部影响因子的逻辑回归系数；通过逻辑回归模型公式式2计算研究区境内滑坡空间发生概率；

其中，P为某一评价因子相对于某事件的发生概率，取值范围为(0，1)；x₁、x₂、x₃、……、x_n分别为各评价因子的取值；β₁、β₂、β₃、……、β_n为多元回归系数，表示评价因子对事件发生P概率的说明程度；

c、通过地理信息系统空间分析功能，提取泥石流的内部影响因子，确定泥石流的各个内部影响因子的无量纲值，通过统计产品与服务解决方案软件实现逻辑回归分析过程，确定泥石流的各个内部影响因子的逻辑回归系数，通过逻辑回归模型公式式2计算研究区境内泥石流空间发生概率；

d、最后通过式3和地理信息系统空间叠置运算确定基于降雨频率的地质灾害发生频率；

R_p＝S×H_24p 式3；

其中，R_p为研究区境内雨量监测点P频率降雨下地质灾害发生频率，S为研究区内地质灾害的敏感性，H_24P为研究区境内雨量监测点P频率降雨下的年最大日降雨量。

2.根据权利要求1所述的一种基于降雨的地质灾害发生频率的计算方法，其特征在于：所述地理信息系统空间分析功能为ArcGIS中的ArcToolbox-3D Analyst-栅格表面功能。

3.根据权利要求1所述的一种基于降雨的地质灾害发生频率的计算方法，其特征在于：所述步骤b中滑坡的内部影响因子为距断层距离、地形坡度、地层岩性、距水系的距离、海拔高度和地震烈度。

4.根据权利要求1所述的一种基于降雨的地质灾害发生频率的计算方法，其特征在于：所述步骤b中滑坡的各个内部影响因子的无量纲值，通过式4计算；

其中，S为研究区面积，km²，S_i为某评价因子i等级所占的面积，km²，N为研究区内滑坡的总数，N_i为某评价因子i等级中滑坡的个数。

5.根据权利要求1所述的一种基于降雨的地质灾害发生频率的计算方法，其特征在于：所述步骤c中泥石流的内部影响因子为流域坡度、沟壑密度、物源条件和流域高差。

6.根据权利要求1所述的一种基于降雨的地质灾害发生频率的计算方法，其特征在于：所述步骤c中泥石流的各个内部影响因子的无量纲值，通过式5计算；

其中，Xi(i,j)为j评价因子中第i个评价单元归一化后的数值，x(i,j)为j评价因子中第i个评价单元的数值，min(x(i,j))为评价因子j中存在的最小值，max(x(i,j))为评价因子j中存在的最大值。

7.根据权利要求1所述的一种基于降雨的地质灾害发生频率的计算方法，其特征在于：所述步骤d中通过式3和地理信息系统空间叠置运算确定基于降雨频率的地质灾害发生频率，是指通过ArcGIS中的ArcToolbox-3D Analyst-栅格计算功能编辑式3，计算出P频率降雨下地质灾害的发生频率。