一种暴雨诱发的浅层滑坡灾害预警阈值确定方法与流程

本发明属于滑坡灾害预警领域，具体地说，涉及一种暴雨诱发的浅层滑坡灾害预警阈值确定方法。

背景技术

滑坡预警是防灾减灾工作的主要手段和关键环节，国内外学者针对降雨诱发滑坡之间的关系做了大量研究，降雨诱发滑坡阈值(临界降雨量)被广泛应用于滑坡预警工作中，确定发生滑坡灾害的临界降雨量值是滑坡预警工作的关键。目前常用滑坡灾害模型可划分为统计学模型和物理机制模型。利用统计模型进行预警的前提是必须拥有长期、比较完整的降雨资料，具有布置密度比较合理的降雨遥控监测网络和先进的数据传输系统。但在华南山区，研究基础薄弱，历史记录少，往往缺乏良好的滑坡发生时间及其对应的累积雨量等资料，统计样本很难形成具有统计意义的数据体系，利用统计方法开展预警工作其误差较大。物理模型法是从降雨入渗、水文水力条件和边坡岩土体失稳机理上着手推导公式，分析得出基于滑坡过程模型的物理性降雨阈值，或者通过研究地质灾害的运动机理和径流过程，确定滑坡临界水文量。但因暴雨诱发的滑坡的机理有别于普通的降雨型滑坡，因此，临界降雨阈值确定方法也与一般降水滑坡不同。一般降雨型滑坡成灾机制认为暴雨主要是通过降水转化为地下水对斜坡稳定性产生影响，反映在增大坡体质量、抬高坡体内地下水位、影响岩土体的抗剪强度、改变水压等几个方面。但是，按照现有的产汇流理论，入渗率主要取决于岩土体的导水系数，而土体饱和后，渗透能力会明显减弱，这从理论上说明，在暴雨期间地下水对降水并不敏感，因此暴雨期间滑坡的产汇流机制区别于一般的降雨型滑坡。因此现行的基于物理机制模型的临界雨量确定方法未考虑暴雨于一般降水在滑坡水文过程中的差异，严重影响了预警结果的准确性。

技术实现要素：

有鉴于此，本发明提供了一种暴雨诱发的浅层滑坡灾害预警阈值确定方法，将斜坡单元和暴雨影响系数引入暴雨诱发的浅层滑坡体测量建模中，提高滑坡预警的工作效率，改进滑坡预警的空间的精度，为滑坡防治研究提供基础数据。

为了解决上述技术问题，本发明公开了一种暴雨诱发的浅层滑坡灾害预警阈值确定方法，包括以下步骤：

步骤一：收集研究区地形数据，建立高精度的数字高程模型dem，并基于dem进行斜坡单元的划分；

步骤二：基于斜坡对暴雨的响应机制，引入大孔隙系数，利用稳态水文模型的地形湿度指数与无限边坡模型相融合的方法，通过假设斜坡处于特定预警等级状态，将地形湿度指数与斜坡稳定性指数计算公式联立，建立基于斜坡单元和临界入渗量的稳定性分析预警模型；

步骤三：对应国家应急救援分级响应体系，设置五个预警等级，分别为ⅰ、ⅱ、ⅲ、ⅳ、ⅴ级，相应的预警信号分布为红色预警、橙色预警、黄色预警、蓝色预警和绿色预警，根据fs指数表征的含义，反推五个等级对应的斜坡稳定状态，计算每个五种预警等级的临界入渗量阈值；

步骤四：通过斜坡体对暴雨的响应机制分析，建立基于斜坡单元的滑坡降雨-入渗模型，计算特定降雨条件下每个斜坡单元的入渗量与降雨量的关系词，通过临界入渗量反推每个斜坡单元的临界降雨量。

进一步地，步骤一中，建立高精度的数字高程模型dem的构建方法，使用地形测绘仪器进行实地野外测绘获得密集高程点数据，在此基础上建立dem；

或，利用航空摄影测量手段，通过航空影像立体相对提取dem；或，利用现有的纸质地形图，通过软件数字化地形图矢量方式dem。

进一步地，步骤一中，基于dem进行斜坡单元的划分方法，采用源头切割法，以山脊线及沟谷线为边界，在gis平台上，经执行斜坡单元划分程式实现斜坡单元划分，划分流程如下：第一步，进行dem水流方向矩阵的计算，查找和充填洼地，生成无洼地dem数据；第二步，从无洼地dem数据求得各像元的流向和汇流量；第三步，通过gis的水文分析功能提取分水线；第四步，通过设定流经某一像元的最小汇流量求得栅格水系，通过栅格-矢量转换生成矢量水系；第五步，确定最小级别集水区域的出水口位置；第六步，由矢量水系分割各沟谷单元生成斜坡单元。

进一步地，步骤二中，斜坡的安全系数为抗滑力和下滑力之比，根据极限平衡原理在考虑滑体重力、静水压力、地震力、及动水压力(渗透压力)时，滑坡稳定性指数表示为：

其中：fs—安全系数；c—内聚力(kpa)；γ—滑坡土体容重(kn/m³)；γw—水容重(kn/m³)；z—滑坡土体垂直厚度(m)；zw—地下水位距滑面的高度(m)；θ—滑面倾角(o)；φ—内摩擦角(o)；fe为水平地震力(n)，fw＝zwlγwcosθsinθ为渗透压力(动水压力)；

当忽略地震力作用，并将分子分母同时除以zlγs：

令表示滑坡体水与土的密度比值；代表滑坡体的饱和状态，为饱和因子，将上公式进一步简化为：

根据暴雨入渗机制分析，对于暴雨型滑坡，应以形成地表径流为节点，考虑两种情况，在没形成地表径流时，地下水入渗量和稳定性系数受降雨量控制，在形成地表径流后，地下水入渗量取决于大孔隙容量；在暴雨形成径流后，土体饱和，降水入渗量主要受饱和土体导水系数影响，在这种情况下，假设为滑坡体为非渗透性，而滑面以下为临时承压水含水层，地下水沿承压水层进行深层流动循环；在暴雨作用下，大孔隙流产生的作用于滑坡面的浮托力(fb)作用于饱和的岩土体，浮托力是暴雨型滑坡不可忽略的因子。

fb＝γwzwlcosθ

设大孔隙系数为λ(％)，

则大孔隙产生的静水压力为：

ft＝λγwzlcos²θ

因此在考虑临时地下水位形成后滑体重力、静水压、渗透压力和浮托力综合作用下，则稳定性系数为：

将以上公式进一步简化为：

综合以上两种情况可得：

其中：

公式中，q为入渗量，a为上坡集水面积；b为集水区排泄口宽度；θ为边坡倾角，t为导水系数。

当fs＝1时：

则fs＝1时：

进一步地，步骤三中，根据fs指数表征的含义，定义稳定指数对应的四个预警等级；斜坡稳定性指数是判断斜坡稳定性的指标，当fs>1.0时，斜坡基本处于稳定状态，1.0>fs>0.5时，斜坡失稳的可能性低于50％，0.5>fs>0时，斜坡失稳可能性高于50％，fs<＝0时，斜坡处于极不稳定状态。因此，将概率函数引入滑坡灾害预警过程中，一方面可以解决偏远地区的资料缺乏的问题，另一方面也通过概率区间的运算减少预警过程的数据量，得到合理性的预警结果。定义滑坡稳定性指标si定义为稳定性系数fs的概率函数

si＝prob(fs＞1)

令q/t＝x,tanφ＝t,当内聚力与摩擦角的最小值及降雨参数x的最大值导致了斜坡失稳的最有利条件，即fsmin，如果在这种情况下，稳定性系数值仍然大于1，那么该斜坡为无条件稳定斜坡，无论降雨量多大，都为稳定状态，即si＝prob(fs＞1)＝1。

内聚力与摩擦角的最大值及降雨参数x的最小值是保证斜坡稳定的最有利条件，即fsmax，如果在这种情况下，稳定性系数值仍然小于1，那么该斜坡为无条件不稳定斜坡，即为极不稳定斜坡，即si＝prob(fs＞1)＝0。

基于以上判别标准，来界定各稳定状态下的临界入渗量值。因此假设滑坡处于各临界状态时，即fs＝1，fs＝0.5，时，将以上的降雨入渗量关系和入渗量与斜坡稳定性关系进行关联，即可获取引发滑坡的临界降雨入渗量值。

当fs＝1时，将式上式联立得到：

则fs＝1时：

以此类推，可得到fs＝0.5临界雨量值：

从以上公式可知，当φ<θ时，q<0，也就是说在天然状态下，无论有无降水，斜坡已经处于不稳定状态；实际上，在θ<50时，地势平缓，斜坡处于稳定状态；在5<θ<φ时，斜坡处于概率不稳定状态，可将实际降雨量与以上临界降雨量做对比，判断斜坡的稳定状态。

进一步地，步骤四中，建立的基于斜坡单元的滑坡降雨-入渗模型，主要由以下方法：

降雨下渗容量是岩土体导水能力与有效降雨的函数，即当有效降雨大于导水能力时，按下渗容量入渗，当有效降雨小于导水时，按有效降雨量入渗。在暴雨期间，有效降雨大于斜坡土体导水能力，但在暴雨冲刷作用下，斜坡大空隙通到打开，下渗容量取决于大空隙体积大小。

大气降水降落到地面主要分为三个部：植被截流、下渗补充地下水和形成地表径流。在持续暴雨条件下，空气湿度较大，蒸发量较小，树冠蓄水饱和，暴雨期间，可忽略树冠截流和蒸散发作用。当降雨量小于初损时，降雨量全部转化为入渗量，当降雨量大于初损时，入渗量与降雨量的关系取决于大孔隙容量，设斜坡大孔隙系数为λ，则入渗量与降雨量的关系为：

其中：

上式中，q—入渗量，mm；p—次降雨总量，mm；λ—大孔隙系数，a—斜坡单元面积；s—流域当时的最大可能滞留量，也称为保持系数，它是后损的上限，mm。cn为反映降雨前流域特征的参数，考虑了下垫面因素，并与流域前期土壤湿度、坡度、植被、土壤类型和土壤利用现状等有关，反映流域前期土壤湿润程度、坡度、土壤类型和土地利用现状的综合特性，可以较好地反映下垫面条件对产汇流过程的影响。

与现有技术相比，本发明可以获得包括以下技术效果：

1)本发明选用斜坡单元作为滑坡灾害预警的基本空间单元，斜坡单元一方面能够保证成灾过程的完整性，一方面又可以提高灾害预警的空间精度。

2)本发明是基于灾害的形成机理分析，引入暴雨激发的大孔隙系数，选用将气象水文模型和水文地质模型耦合的方法计算滑坡预警临界降雨阈值，可提高灾害预警精度。

3)本发明提供的基于斜坡单元的预警阈值确定方法可以应用于山区滑坡灾害预警工作中，为山区防灾减灾工作提供科学依据，在防灾决策领域具有广阔的应用前景。

当然，实施本发明的任一产品并不一定需要同时达到以上所述的所有技术效果。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是本发明暴雨诱发的浅层滑坡灾害预警阈值确定方法流程图；

图2是本发明实施例中研究区dem；

图3是本发明实施例的斜坡单元划分结果；

图4是本发明实施例的斜坡单元平均坡度图；

图5是本发明实施例的斜坡单元的比集水面积图；

图6是本发明实施例的斜坡单元土层平均厚度分布图；

图7是本发明实施例的斜坡单元的导水系数分布图；

图8是本发明实施例的斜坡单元有效粘聚力分布图；

图9是本发明实施例的斜坡单元暴雨影响系数图；

图10是本发明实施例的各斜坡单元预警级临界降雨量。

具体实施方式

以下将配合实施例来详细说明本发明的实施方式，藉此对本发明如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。

本发明公开了一种暴雨诱发的浅层滑坡灾害预警阈值确定方法，如图1所示，包括以下步骤：

步骤一：收集研究区地形数据，建立高精度的数字高程模型dem，并基于dem进行斜坡单元的划分；

建立高精度的数字高程模型dem的构建方法，使用地形测绘仪器进行实地野外测绘获得密集高程点数据，在此基础上建立dem；

或，利用航空摄影测量手段，通过航空影像立体相对提取dem；或，利用现有的纸质地形图，通过软件数字化地形图矢量方式dem。

基于dem进行斜坡单元的划分方法：采用源头切割法，以山脊线及沟谷线为边界，在gis平台上，经执行斜坡单元划分程式实现斜坡单元划分，划分流程如下：第一步，进行dem水流方向矩阵的计算，查找和充填洼地，生成无洼地dem数据；第二步，从无洼地dem数据求得各像元的流向和汇流量；第三步，通过gis的水文分析功能提取分水线；第四步，通过设定流经某一像元的最小汇流量求得栅格水系，通过栅格-矢量转换生成矢量水系；第五步，确定最小级别集水区域的出水口位置；第六步，由矢量水系分割各沟谷单元生成斜坡单元。

步骤二：基于斜坡对暴雨的响应机制，引入大孔隙系数，利用稳态水文模型的地形湿度指数与无限边坡模型相融合的方法，通过假设斜坡处于特定预警等级状态，将地形湿度指数与斜坡稳定性指数计算公式联立，建立基于斜坡单元和临界入渗量的稳定性分析预警模型；

斜坡的安全系数为抗滑力和下滑力之比。根据极限平衡原理在考虑滑体重力、静水压力、地震力、及动水压力(渗透压力)时，滑坡稳定性指数表示为：

其中：fs—安全系数；c—内聚力(kpa)；γ—滑坡土体容重(kn/m³)；γw—水容重(kn/m³)；z—滑坡土体垂直厚度(m)；zw—地下水位距滑面的高度(m)；θ—滑面倾角(o)；φ—内摩擦角(o)；fe为水平地震力(n)，fw＝zwlγwcosθsinθ为渗透压力(动水压力)。

当忽略地震力作用，并将分子分母同时除以zlγs：

令表示滑坡体水与土的密度比值；代表滑坡体的饱和状态，为饱和因子。将上公式进一步简化为：

根据暴雨入渗机制分析，对于暴雨型滑坡，应以形成地表径流为节点，考虑两种情况，在没形成地表径流时，地下水入渗量和稳定性系数受降雨量控制，在形成地表径流后，地下水入渗量取决于大孔隙容量。在暴雨形成径流后，土体饱和，降水入渗量主要受饱和土体导水系数影响，在这种情况下，可假设为滑坡体为非渗透性，而滑面以下为临时承压水含水层，地下水沿承压水层进行深层流动循环。在暴雨作用下，大孔隙流产生的作用于滑坡面的浮托力(fb)作用于饱和的岩土体，浮托力是暴雨型滑坡不可忽略的因子。

fb＝γwzwlcosθ

设大孔隙系数为λ(％)，

则大孔隙产生的静水压力为：

ft＝λγwzlcos²θ

因此在考虑临时地下水位形成后滑体重力、静水压、渗透压力和浮托力综合作用下，则稳定性指数为：

将以上公式进一步简化为：

综合以上两种情况可得：

其中：

公式中，q可定为入渗量，a为上坡集水面积；b为集水区排泄口宽度；θ为边坡倾角，t为导水系数。

当fs＝1时：

则fs＝1时：

步骤三：对应国家应急救援分级响应体系，设置五个预警等级，分别为ⅰ、ⅱ、ⅲ、ⅳ、ⅴ级，相应的预警信号分布为红色预警、橙色预警、黄色预警、蓝色预警和绿色预警(无预警)，根据fs指数表征的含义，反推五个等级对应的斜坡稳定状态，计算每个五种预警等级的临界入渗量阈值；

根据fs指数表征的含义，定义稳定指数对应的五个预警等级；斜坡稳定性指数是判断斜坡稳定性的指标，当fs>1.5时，斜坡处于及稳定状态，当fs>1.0时，斜坡基本处于稳定状态，1.0>fs>0.5时，斜坡失稳的可能性低于50％，0.5>fs>0时，斜坡失稳可能性高于50％，fs<＝0时，斜坡处于极不稳定状态。因此，将概率函数引入滑坡灾害预警过程中，一方面可以解决偏远地区的资料缺乏的问题，另一方面也通过概率区间的运算减少预警过程的数据量，得到合理性的预警结果。将滑坡稳定性指标si定义为稳定性系数fs的概率函数：

si＝prob(fs＞1)

令q/t＝x，tanφ＝t，当内聚力与摩擦角的最小值及降雨参数x的最大值导致了斜坡失稳的最有利条件，即fsmin，如果在这种情况下，稳定性系数值仍然大于1，那么该斜坡为无条件稳定斜坡，无论降雨量多大，都为稳定状态，即si＝prob(fs＞1)＝1。

内聚力与摩擦角的最大值及降雨参数x的最小值是保证斜坡稳定的最有利条件，即fsmax，如果在这种情况下，稳定性系数值仍然小于1，那么该斜坡为无条件不稳定斜坡，即为极不稳定斜坡，即si＝prob(fs＞1)＝0。将si分为几个区间分别代表滑坡的不同稳定性，分级标准如表(1)。

表1斜坡稳定性特征与预警等级关系表

基于以上判别标准，来界定各稳定状态下的临界入渗量值。因此假设滑坡处于各临界状态时，即fs＝1，fs＝0.5时，将以上的降雨入渗量关系和入渗量与斜坡稳定性关系进行关联，即可获取引发滑坡的临界降雨入渗量值。

当fs＝1时，将上式联立得到：

则fs＝1时：

以此类推，可得到fs＝0.5临界雨量值：

从以上公式可知，当φ<θ时，q<0，也就是说在天然状态下，无论有无降水，斜坡已经处于不稳定状态；实际上，在θ<50时，地势平缓，斜坡处于稳定状态；在5<θ<φ时，斜坡处于概率不稳定状态，可将实际降雨量与以上临界降雨量做对比，判断斜坡的稳定状态。其他状态可以此类推。

步骤四：通过斜坡体对暴雨的响应机制分析，建立基于斜坡单元的滑坡降雨-入渗模型，计算特定降雨条件下每个斜坡单元的入渗量与降雨量的关系词，通过临界入渗量反推每个斜坡单元的临界降雨量。

建立基于斜坡单元的滑坡降雨-入渗模型，主要由以下方法：

降雨下渗容量是岩土体导水能力与有效降雨的函数，即当有效降雨大于导水能力时，按下渗容量入渗，当有效降雨小于导水时，按有效降雨量入渗。在暴雨期间，有效降雨大于斜坡土体导水能力，但在暴雨冲刷作用下，斜坡大空隙通到打开，下渗容量取决于大空隙体积大小。

大气降水降落到地面主要分为三个部：植被截流、下渗补充地下水和形成地表径流。在持续暴雨条件下，空气湿度较大，蒸发量较小，树冠蓄水饱和，暴雨期间，可忽略树冠截流和蒸散发作用。当降雨量小于初损时，降雨量全部转化为入渗量，当降雨量大于初损时，入渗量与降雨量的关系取决于大孔隙容量，设斜坡大孔隙系数为λ，则入渗量与降雨量的关系为：

其中：

上式中，q—入渗量，mm；p—次降雨总量，mm；λ—为大孔隙系数，a—为斜坡单元面积；s—流域当时的最大可能滞留量，也称为保持系数，它是后损的上限，mm。cn—反映降雨前流域特征的参数，考虑了下垫面因素，并与流域前期土壤湿度、坡度、植被、土壤类型和土壤利用现状等有关，反映流域前期土壤湿润程度、坡度、土壤类型和土地利用现状的综合特性，可以较好地反映下垫面条件对产汇流过程的影响。

上述各参数的获取方法如下：p为次降雨总量可通过实测降雨量或者预测降雨量获取。a斜坡单元面积是通过gis空间分析力的面积统计工具，统计每个斜坡单元的面积。λ大空隙系数可通过野外实验如示踪法获取，或者按照地表覆盖和土壤类型模拟计算获取。cn值可以由历史自然降水径流水文资料统计得到，也根据所在流域的土地利用方式、处理情况、水文条件以及土壤类别，利用cn经验值，其值介于0-100。

在实际运算中，需获取包含斜坡的水文、气象及下垫面等参数，模型中的地形参数坡度(θ)和比集水面积(a)，利用dem，通过gis空间分析功能运算获得。因地质参数的大小主要取决于岩性类别，因此可按照岩性类型将研究区进行分区，在此基础上进行地质参数的反算工作。根据查阅历史文献和现场调查实验，确定了各地质区的平均地质参数。其中土壤容重(γ)和内摩擦角(φ)在每个地质区内基本稳定，可直接采纳各个分区的平均值。导水系数(t)根据渗透系数和坡度值推算、粘聚力(c)根据岩土体的粘聚力和植物粘聚力综合运算能获取、前期土壤湿润程度(amc)、斜坡特征参数cn，最大滞留量s，初始入渗量ia根据斜坡岩土体水文特征和scs模型附表。将每个斜坡单元的14个参数分别提取，构成预警模型数据库，然后运用gis的空间代数功能，按照前述的计算过程，分别计算临界有效降雨量(q)和临界降雨量(p)。

实施例

(1)根据地形图资料构建目标区域dem

小流域dem的数据运算是在arcgis平台上进行。dem的建立主要方法：1)矢量化1：1万地形图，建立空间拓扑关系；2)采用三角化不规则网(tin)方法，将矢量线划图转换成tin数据结构，进行空间插值；3)tin数据结构转换成格网结构；4)对格网数据进行重采样形成栅格dem(grid格式)，生成小流域dem；5)将dem数据再内插生成等高线与等高线扫描影像分色叠合，对生成的dem要进行粗差检查，检查有无超限偏移。研究区dem如附图2。

(2)斜坡单元的划分

以1:1万地形图生成的dem为斜坡单元划分的基础，采用源头切割法，在gis平台上，分别提取最低级别的水文单元的以山脊线及沟谷线，并以此为边界，构成斜坡单元。图3为斜坡单元划分结果。

(3)模型参数设定

模型中的地形参数坡度(θ)(图4为斜坡单元平均坡度图)和比集水面积(a)(图5为斜坡单元的比集水面积)，利用dem，通过gis空间分析功能运算获得。因地质参数的大小主要取决于岩性类别，根据查阅历史文献和现场调查实验，确定了各地质区的平均地质参数，图6为斜坡单元土层平均厚度分布图。土壤容重(γ)和内摩擦角(φ)在每个地质区内基本稳定，因此，直接采纳各个分区的平均值。导水系数(t)(图7为斜坡单元的导水系数分布图)、粘聚力(c)(图8为斜坡单元有效粘聚力分布图)值通过gis和相关的研究结果进行计算和反推。导水系数t＝kh＝kzcosθ，k通过实验测序而来，因土层厚度和坡度具有较高的相关性，为求得松岗河小流域的土层厚度分布情况，进行了土壤深度的野外量测工作，主要的量测项目为土层厚度和地表坡度，根据调查数据对松岗河小流域地区坡度和土层厚度之间的关系进行回归模拟得到y＝-2.864ln(x)+10.719。根据单元平均坡度值可反推土层厚度分布，进而获得导水系数分布。土壤的有效粘聚力是土壤粘聚力与植物根系粘聚力的组合，因此本研究根据植被覆盖的情况，结合几个地质分区的土壤类型及粘聚力数值分布情况，运算土壤的有效粘聚力。以植被指数推算有效粘聚力的方式确定有效粘聚力的运算：

其中c＇为有效粘聚力，cmax为最大有效粘聚力，ndvi为植被指数。式中cmax通过查阅文献资料和野外调查获得表1，ndvi通过遥感图像反演计算获得。最后经过gis的地图代数功能运算每个斜坡单元的土壤有效粘聚力。

大孔隙系数(λ)对斜坡的作用主要表现为大孔隙流的影响，如图9。根据文献调查和根据研究区的岩土特征，假设斜坡单元的影响系数为0-0.04。在研究区中部的相对平缓地区，大孔隙的影响作用不明显。

根据小流域的土地利用方式和土壤、水文情况，小流域内土地利用方式有道路和村落、耕地、林地、耕地几种土地处理情况，根据土壤渗透速率及现场调查确定土壤类型为d类；因计算的是极端气候条件下的滑坡的预警范围，因此土壤湿润程度为湿，经过综合考虑，根据所在流域的土地利用方式、处理情况、水文条件以及土壤类别，确定和调整流域的cn值，并推算出s值和ia值。

(4)计算临界降雨量值

新建一个存储临界降雨量值的q字段，运用属性表的字段运算(fieldcalculator)功能，按照公式在上述参数设定模式下，每个斜坡单元的临界降雨量值。结合气象预报对区域未来一个时段(通常为24h)降水预报，结合该区域地质灾害的潜在变化，应用临界降雨量值，对该时段内地质灾害发生的危害等级做出的预报。运算过程是将预报降雨量值域临界降雨量值相比较，依据松岗河小流域的灾害发生情况，运算不同等级条件下的临界降雨量q值，如图10。按照输入降雨量对单元进行预警。也就是以斜坡为单元，将气象预报的降雨量p与q1，q2，q3，q4，q5分别做对比，预警等级及显示标准,当p>q1为警报级，当q1>p>q2为预报级，当q2>p>q3为预警级，当q3>p>q4为预测级，当q4>p>q5为常规级。

利用arcgis的栅格计算功能，将上述参数代入如前所述的预警模型中，计算每个斜坡单元的临界降雨量，研究区临界降雨量值范围为0-517mm。临界降雨量较低的区域，也是滑坡灾害的高易发区。

上述说明示出并描述了发明的若干优选实施例，但如前所述，应当理解发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离发明的精神和范围，则都应在发明所附权利要求的保护范围内。