[发明专利]一种震后泥石流沟敏感性判别方法及其应用

摘要

本发明公开了一种震后泥石流沟敏感性判别方法及其应用。该方法先确定研究区内所有崩滑体在6个评价因子的任意分级条件下的概率综合判别值P；根据P值大小对原本任意的分级条件进行整合，得到6个评价因子的敏感性高低分级特定区间，并建立相应打分标准；然后通过单沟在6个评价因子的分布情况，参照打分标准得到该沟在6个评价因子上的分值；最后结合灰色关联法确定评价因子权重，建立单沟泥石流敏感性判别模型，进行泥石流预测预报。与现有技术相比，本发明充分考虑了单沟所在研究区的区域特征，避免了人为划分单沟泥石流的评价因子敏感性分级区间的误差，能快速准确地对评价因子进行特定等级划分，从而实现对单沟泥石流敏感性的有效判别。

主权项

一种震后泥石流沟敏感性判别方法，所述震后泥石流沟的泥石流暴发是由地震引起的崩滑体起动形成，其特征在于：所述震后泥石流沟的敏感性判别要考虑所在研究区的区域特征，具体判别方法步骤如下：A.通过1:5万地形图矢量化得到整个研究区的高程因子、沟床纵比降因子、坡度因子、坡向因子和沟壑密度因子；通过1:20万地形图矢量化得到整个研究区的地层岩性因子；利用ARCGIS软件对得到的6个评价因子分别进行任意的区间分级，得到6个评价因子的不同分级条件；B.通过1:5万地形图得到整个研究区的总面积，通过遥感影像解译出整个研究区内所有崩滑体面积；针对步骤A中得到的6个评价因子的不同分级条件，通过遥感影像分别解译出整个研究区内、6个评价因子不同分级条件下的崩滑体面积，并通过1:5万地形图分别得到6个评价因子不同分级条件下的研究区面积；C.通过以下公式分别计算得到6个评价因子不同分级条件下的概率综合判别值P： P = p 1 - p 2 p 1 * ( 1 - p 2 ) p 1 ≥ p 2 p 1 - p 2 p 2 * ( 1 - p 1 ) p 1 < p 2 式中，P—6个评价因子不同分级条件下的概率综合判别值；P1—整个研究区内、6个评价因子不同分级条件下的崩滑体面积占该级别研究区面积的比例；整个研究区内、6个评价因子不同分级条件下的崩滑体面积，6个评价因子不同分级条件下的研究区面积，均由步骤B确定；P2—整个研究区内所有崩滑体面积占整个研究区总面积的比例；整个 研究区内所有崩滑体面积，整个研究区的总面积，均由步骤B确定；D.针对步骤C中得到的6个评价因子不同分级条件下的概率综合判别值P，逐一选取其中1个评价因子不同分级条件下的概率综合判别值P进行如下a）‑d）步骤，划分敏感性高低分级的特定区间，并建立打分标准：a)确定特定区间的个数，将特定区间个数减1，得到特定区间个数的间隔；将不同分级条件下的概率综合判别值P按大小排序；b)将P值最小值对应的分级条件作为最低敏感性区间，将P值最大值对应的分级条件作为最高敏感性区间；c)将P值次大值与P值次小值进行差值运算，再除以特定区间个数的间隔，得到分级区间递增值；从P值次小值开始，根据分级区间递增值依次确定特定区间的组合范围，得到除最低敏感性区间和最高敏感性区间外的其余分级特定区间；d)从最低敏感性区间开始，根据P值递增顺序，依次对得到的特定区间进行打分，分值由小到大，建立泥石流在所选评价因子的敏感性高低分级特定区间上的打分标准；最终得到6个评价因子的敏感性高低分级特定区间，及相应打分标准；E.针对研究区中待进行敏感性判别的震后泥石流沟，通过1:5万地形图矢量化得到该条泥石流沟的高程因子、沟床纵比降因子、坡度因子、坡向因子和沟壑密度因子，通过1:20万地形图矢量化得到该条泥石流沟的地层岩性因子；F.根据步骤E中得到的该沟在6个评价因子上的分布情况，参照步骤D中得到的6个评价因子的敏感性高低分级特定区间及相应打分标准，得到该沟在6个评价因子上的分值；G.根据灰色关联法基本公式，计算得到6个评价因子的权重值；H.通过以下公式确定待进行敏感性判别的震后泥石流沟的敏感性大小： R = Σ i = 1 6 x i ( k ) w i 式中，xi(k)为待进行敏感性判别的震后泥石流沟在6个评价因子上的分值，由步骤F确定；wi为6个评价因子的权重值，由步骤G确定；R为敏感性大小判别值；当R小于等于2.0，则待进行敏感性判别的震后泥石流沟为低敏感性泥石流流域；当R大于2.0同时小于2.6，则待进行敏感性判别的震后泥石流沟为中敏感性泥石流流域；当R大于等于2.6，则待进行敏感性判别的震后泥石流沟为高敏感性泥石流流域。