本發(fā)明涉及到地質(zhì)災害防治工程技術領域，尤其涉及一種基于降雨的地質(zhì)災害發(fā)生頻率的計算方法。

背景技術：

基于降雨頻率的地質(zhì)災害發(fā)生頻率是強震區(qū)地質(zhì)災害較為有效的非工程減災措施之一。然而，目前基于降雨頻率的地質(zhì)災害發(fā)生頻率的計算還沒有得到完善。國內(nèi)外雖已有大量關于地質(zhì)災害發(fā)生頻率的研究成果，但前人對地質(zhì)災害發(fā)生頻率的評價大多只提出了空間概率，并未對潛在受災民眾指出地質(zhì)災害在什么時候可能會發(fā)生，即缺少時間概率。此缺陷致使地質(zhì)災害發(fā)生頻率的評價方法未能滿足地質(zhì)災害發(fā)生的“天地耦合”條件。

如公開號為cn104318103a，公開日為2015年01月28日的中國專利文獻公開了一種滑坡災害監(jiān)測預警降雨閾值判定方法，其特征在于，其包括以下步驟：(1)已有勘查和降雨資料統(tǒng)計分析：主要包括滑坡破壞時的降雨量數(shù)據(jù)，滑坡變形破壞時地形地貌形態(tài)、破壞堆積范圍、崩滑過程、巖土體物理力學參數(shù)等方面資料；(2)建立滑坡破壞模型：結合典型降雨滑坡現(xiàn)場監(jiān)測數(shù)據(jù)，建立臨界降雨量與地下水位、孔隙水壓力、體積含水量之間的滲流水位模型，分析變形特點及形成機理，建立不同降雨滑坡破壞模型；(3)室內(nèi)模擬試驗：針對野外監(jiān)測數(shù)據(jù)的不足，須利用模型試驗的方法，對不同降雨情況下的典型降雨滑坡的變形破壞機制進行模擬試驗研究，并結合野外監(jiān)測數(shù)據(jù)初步建立降雨型滑坡的臨界降雨閾值；(4)試驗觀測資料整理分析及補充實驗：將上述各種資料包括調(diào)查、模型實驗及原型觀測等資料進行整理，分析降雨型滑坡變形破壞過程、分析破壞時的臨界雨量、地表位移及深部位移等數(shù)據(jù)，為理論分析和模型的建立奠定基礎；資料整理分析后如發(fā)現(xiàn)不足，需進行補充實驗及觀測；(5)理論分析臨界值確定：綜合運用工程地質(zhì)學、土力學、滲流力學等理論方法對降雨滑坡變形破壞過程及機理進行綜合分析；結合模型試驗數(shù)據(jù)和現(xiàn)場監(jiān)測數(shù)據(jù)，建立典型降雨滑坡變形破壞過程的降雨臨界值。

又如公開號為cn104036153a，公開日為2014年09月10日的中國專利文獻公開了一種泥石流發(fā)生的危險性定量檢測評價方法，其特征在于，包括以下步驟：a.根據(jù)統(tǒng)計的降雨量與泥石流發(fā)生可能性大小，建立降雨量大小與泥石流發(fā)生可能性相關關系，并按一定權重賦予降雨量相應分值，降雨量分值大小代表降雨量因素引起泥石流的可能性大?。籦.根據(jù)統(tǒng)計的泥石流溝縱坡坡降與泥石流發(fā)生可能性大小，建立泥石流縱坡坡降與泥石流發(fā)生可能性相關關系，并按一定權重賦予泥石流溝縱坡坡降相應分值，縱坡坡降分值大小代表縱坡坡降因素引起泥石流的可能性大小；c.根據(jù)統(tǒng)計的泥石流溝內(nèi)不穩(wěn)定物源量與泥石流發(fā)生可能性大小，建立不穩(wěn)定物源量與泥石流發(fā)生可能性相關關系，并按一定權重賦予不穩(wěn)定物源量相應分值，不穩(wěn)定物源量分值大小代表不穩(wěn)定物源量因素引起泥石流的可能性大??；d.綜合降雨量、泥石流溝縱坡坡降和泥石流溝內(nèi)不穩(wěn)定物源量這三個因素對泥石流發(fā)生的可能性大小的影響，采用相關分析方法，建立三個因素得分值之和與三個因素單項值大小之間的相關關系式；并根據(jù)總分值大小劃出不同危險性大小等級。

以上述專利文獻為代表的現(xiàn)有技術，對地質(zhì)災害發(fā)生頻率的評價大多只提出了空間概率，并未對潛在受災民眾指出地質(zhì)災害在什么時候可能會發(fā)生，即缺少時間概率，致使地質(zhì)災害預測精度較低，降低了應對災害的主動性。

技術實現(xiàn)要素：

本發(fā)明為了克服上述現(xiàn)有技術的缺陷，提供一種基于降雨的地質(zhì)災害發(fā)生頻率的計算方法，本發(fā)明明確了地質(zhì)災害發(fā)生的空間概率和時間概率，極大的提高了地質(zhì)災害的預測精度和應對災害的主動性。

本發(fā)明通過下述技術方案實現(xiàn)：

一種基于降雨的地質(zhì)災害發(fā)生頻率的計算方法，其特征在于，包括以下步驟：

a、通過中國氣象局確定研究區(qū)境內(nèi)的雨量監(jiān)測站點，通過中小流域暴雨洪水計算手冊中最大24小時暴雨均值等值線圖確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的點雨量均值變差系數(shù)cv，通過cs＝3.5cv適線確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的離均系數(shù)kp，通過式1確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p，通過畫圖軟件和克里金插值確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量等值線圖；

b、通過地理信息系統(tǒng)空間分析功能，提取滑坡的內(nèi)部影響因子，通過滑坡地質(zhì)災害相對密度確定滑坡的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定滑坡的各個內(nèi)部影響因子的邏輯回歸系數(shù)；通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)滑坡空間發(fā)生概率；

其中，p為某一評價因子相對于某事件的發(fā)生概率，取值范圍為(0，1)；x1、x2、x3、……、xn分別為各評價因子的取值；β1、β2、β3、……、βn為多元回歸系數(shù)，表示評價因子對事件發(fā)生p概率的說明程度；

c、通過地理信息系統(tǒng)空間分析功能，提取泥石流的內(nèi)部影響因子，確定泥石流的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定泥石流的各個內(nèi)部影響因子的邏輯回歸系數(shù)，通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)泥石流空間發(fā)生概率；

d、最后通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率；

rp＝s×h24p式3；

其中，rp為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下地質(zhì)災害發(fā)生頻率，s為研究區(qū)內(nèi)地質(zhì)災害的敏感性，h24p為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量。

所述地理信息系統(tǒng)空間分析功能為arcgis中的arctoolbox-3danalyst-柵格表面功能。

所述步驟b中滑坡的內(nèi)部影響因子為距斷層距離、地形坡度、地層巖性、距水系的距離、海拔高度和地震烈度。

所述步驟b中滑坡的各個內(nèi)部影響因子的無量綱值，通過式4計算；

其中，s為研究區(qū)面積，km²，si為某評價因子i等級所占的面積，km²，n為研究區(qū)內(nèi)滑坡的總數(shù)，ni為某評價因子i等級中滑坡的個數(shù)。

所述步驟c中泥石流的內(nèi)部影響因子為流域坡度、溝壑密度、物源條件和流域高差。

所述步驟c中泥石流的各個內(nèi)部影響因子的無量綱值，通過式5計算；

其中，xi(i,j)為j評價因子中第i個評價單元歸一化后的數(shù)值，x(i,j)為j評價因子中第i個評價單元的數(shù)值，min(x(i,j))為評價因子j中存在的最小值，max(x(i,j))為評價因子j中存在的最大值。

所述步驟d中通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率，是指通過arcgis中的arctoolbox-3danalyst-柵格計算功能編輯式3，計算出p頻率降雨下地質(zhì)災害的發(fā)生頻率。

本發(fā)明，適用于現(xiàn)有的水文站點密度程度及掌握的降雨相關數(shù)據(jù)資料不足以建立足夠精度的水文模型的區(qū)域，且研究區(qū)面積較大的強震區(qū)域的地質(zhì)災害發(fā)生頻率預測。

本發(fā)明，適用于不同降雨頻率的地質(zhì)災害發(fā)生頻率預測。

本發(fā)明的有益效果主要表現(xiàn)在以下方面：

一、本發(fā)明，采用“a-d”步驟形成的一個完整技術方案，回避了水文站點密度程度和雨量數(shù)據(jù)監(jiān)測無法有效覆蓋山區(qū)復雜環(huán)境的問題，通過研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p、滑坡內(nèi)部影響因子和泥石流內(nèi)部影響因子，提供了一種基于降雨的地質(zhì)災害發(fā)生頻率計算方法；明確了地質(zhì)災害的發(fā)生空間概率和時間概率，針對國內(nèi)外地質(zhì)災害發(fā)生頻率中時間尺度的不足，能夠向潛在受災人群揭示地質(zhì)災害“哪里可能發(fā)生”，即空間性；以及“什么時候可能發(fā)生”，即時間性；極大的提高了地質(zhì)災害的預測精度和應對災害的主動性。

二、本發(fā)明，研究區(qū)境內(nèi)的雨量監(jiān)測站點通過中國氣象局查詢獲??；研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的點雨量均值變差系數(shù)cv通過中小流域暴雨洪水計算手冊獲得；研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p通過式1獲得，具有參數(shù)獲取容易、計算過程簡便的效果。

三、本發(fā)明，步驟b中，地理信息系統(tǒng)是一種專門用于采集、存儲、管理、分析和表達空間數(shù)據(jù)的信息系統(tǒng)，是一種基于計算機的工具，它可以對空間信息進行分析和處理。

四、本發(fā)明，步驟b中，通過式4計算滑坡的各個內(nèi)部影響因子的無量綱值，具有計算目標明確、計算結果獲取簡便的效果。

五、本發(fā)明，步驟c中，通過式5計算泥石流的各個內(nèi)部影響因子的無量綱值，不僅計算目標明確，而且其計算結果獲取相當簡便。

六、本發(fā)明，步驟d中通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率，是指通過arcgis中的arctoolbox-3danalyst-柵格計算功能編輯式3，計算出p頻率降雨下地質(zhì)災害的發(fā)生頻率，整個計算過程通過程序化處理，其計算結論客觀性更強。

七、本發(fā)明，適用于現(xiàn)有的水文站點密度程度及掌握的降雨相關數(shù)據(jù)資料不足以建立足夠精度的水文模型的區(qū)域，且研究區(qū)面積較大的強震區(qū)域的地質(zhì)災害發(fā)生頻率預測，具有數(shù)據(jù)預處理流程簡潔、計算過程簡便、預測值直觀、降雨數(shù)據(jù)獲取方法易掌握的特點。

八、本發(fā)明，適用于不同降雨頻率的地質(zhì)災害發(fā)生頻率預測，預測時間指標落腳于不同的降雨頻率，空間指標落腳于現(xiàn)有的地質(zhì)災害庫存，具有預測結果直觀、預測發(fā)生災害的時間和地點易掌握的特點。

附圖說明

下面將結合說明書附圖和具體實施方式對本發(fā)明作進一步的具體說明，其中：

圖1為本發(fā)明推導年最大日降雨量的示意圖；

圖2為本發(fā)明推導不同降雨頻率的地質(zhì)災害發(fā)生頻率示意圖。

具體實施方式

實施例1

一種基于降雨的地質(zhì)災害發(fā)生頻率的計算方法，包括以下步驟：

a、通過中國氣象局確定研究區(qū)境內(nèi)的雨量監(jiān)測站點，通過中小流域暴雨洪水計算手冊中最大24小時暴雨均值等值線圖確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的點雨量均值變差系數(shù)cv，通過cs＝3.5cv適線確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的離均系數(shù)kp，通過式1確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p，通過畫圖軟件和克里金插值確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量等值線圖；

b、通過地理信息系統(tǒng)空間分析功能，提取滑坡的內(nèi)部影響因子，通過滑坡地質(zhì)災害相對密度確定滑坡的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定滑坡的各個內(nèi)部影響因子的邏輯回歸系數(shù)；通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)滑坡空間發(fā)生概率；

其中，p為某一評價因子相對于某事件的發(fā)生概率，取值范圍為(0，1)；x1、x2、x3、……、xn分別為各評價因子的取值；β1、β2、β3、……、βn為多元回歸系數(shù)，表示評價因子對事件發(fā)生p概率的說明程度；

c、通過地理信息系統(tǒng)空間分析功能，提取泥石流的內(nèi)部影響因子，確定泥石流的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定泥石流的各個內(nèi)部影響因子的邏輯回歸系數(shù)，通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)泥石流空間發(fā)生概率；

d、最后通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率；

rp＝s×h24p式3；

其中，rp為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下地質(zhì)災害發(fā)生頻率，s為研究區(qū)內(nèi)地質(zhì)災害的敏感性，h24p為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量。

采用“a-d”步驟形成的一個完整技術方案，回避了水文站點密度程度和雨量數(shù)據(jù)監(jiān)測無法有效覆蓋山區(qū)復雜環(huán)境的問題，通過研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p、滑坡內(nèi)部影響因子和泥石流內(nèi)部影響因子，提供了一種基于降雨的地質(zhì)災害發(fā)生頻率計算方法；明確了地質(zhì)災害的發(fā)生空間概率和時間概率，針對國內(nèi)外地質(zhì)災害發(fā)生頻率中時間尺度的不足，能夠向潛在受災人群揭示地質(zhì)災害“哪里可能發(fā)生”，即空間性；以及“什么時候可能發(fā)生”，即時間性；極大的提高了地質(zhì)災害的預測精度和應對災害的主動性。

實施例2

一種基于降雨的地質(zhì)災害發(fā)生頻率的計算方法，包括以下步驟：

a、通過中國氣象局確定研究區(qū)境內(nèi)的雨量監(jiān)測站點，通過中小流域暴雨洪水計算手冊中最大24小時暴雨均值等值線圖確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的點雨量均值變差系數(shù)cv，通過cs＝3.5cv適線確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的離均系數(shù)kp，通過式1確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p，通過畫圖軟件和克里金插值確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量等值線圖；

b、通過地理信息系統(tǒng)空間分析功能，提取滑坡的內(nèi)部影響因子，通過滑坡地質(zhì)災害相對密度確定滑坡的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定滑坡的各個內(nèi)部影響因子的邏輯回歸系數(shù)；通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)滑坡空間發(fā)生概率；

其中，p為某一評價因子相對于某事件的發(fā)生概率，取值范圍為(0，1)；x1、x2、x3、……、xn分別為各評價因子的取值；β1、β2、β3、……、βn為多元回歸系數(shù)，表示評價因子對事件發(fā)生p概率的說明程度；

c、通過地理信息系統(tǒng)空間分析功能，提取泥石流的內(nèi)部影響因子，確定泥石流的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定泥石流的各個內(nèi)部影響因子的邏輯回歸系數(shù)，通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)泥石流空間發(fā)生概率；

d、最后通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率；

rp＝s×h24p式3；

其中，rp為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下地質(zhì)災害發(fā)生頻率，s為研究區(qū)內(nèi)地質(zhì)災害的敏感性，h24p為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量。

所述地理信息系統(tǒng)空間分析功能為arcgis中的arctoolbox-3danalyst-柵格表面功能。

所述步驟b中滑坡的內(nèi)部影響因子為距斷層距離、地形坡度、地層巖性、距水系的距離、海拔高度和地震烈度。

所述步驟b中滑坡的各個內(nèi)部影響因子的無量綱值，通過式4計算；

其中，s為研究區(qū)面積，km²，si為某評價因子i等級所占的面積，km²，n為研究區(qū)內(nèi)滑坡的總數(shù)，ni為某評價因子i等級中滑坡的個數(shù)。

所述步驟c中泥石流的內(nèi)部影響因子為流域坡度、溝壑密度、物源條件和流域高差。

研究區(qū)境內(nèi)的雨量監(jiān)測站點通過中國氣象局查詢獲?。谎芯繀^(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的點雨量均值變差系數(shù)cv通過中小流域暴雨洪水計算手冊獲得；研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p通過式1獲得，具有參數(shù)獲取容易、計算過程簡便的效果。

步驟b中，地理信息系統(tǒng)是一種專門用于采集、存儲、管理、分析和表達空間數(shù)據(jù)的信息系統(tǒng)，是一種基于計算機的工具，它可以對空間信息進行分析和處理。

步驟b中，通過式4計算滑坡的各個內(nèi)部影響因子的無量綱值，具有計算目標明確、計算結果獲取簡便的效果。

實施例3

一種基于降雨的地質(zhì)災害發(fā)生頻率的計算方法，包括以下步驟：

a、通過中國氣象局確定研究區(qū)境內(nèi)的雨量監(jiān)測站點，通過中小流域暴雨洪水計算手冊中最大24小時暴雨均值等值線圖確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的點雨量均值變差系數(shù)cv，通過cs＝3.5cv適線確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的離均系數(shù)kp，通過式1確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量h24p，通過畫圖軟件和克里金插值確定研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量等值線圖；

b、通過地理信息系統(tǒng)空間分析功能，提取滑坡的內(nèi)部影響因子，通過滑坡地質(zhì)災害相對密度確定滑坡的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定滑坡的各個內(nèi)部影響因子的邏輯回歸系數(shù)；通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)滑坡空間發(fā)生概率；

其中，p為某一評價因子相對于某事件的發(fā)生概率，取值范圍為(0，1)；x1、x2、x3、……、xn分別為各評價因子的取值；β1、β2、β3、……、βn為多元回歸系數(shù)，表示評價因子對事件發(fā)生p概率的說明程度；

c、通過地理信息系統(tǒng)空間分析功能，提取泥石流的內(nèi)部影響因子，確定泥石流的各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定泥石流的各個內(nèi)部影響因子的邏輯回歸系數(shù)，通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)泥石流空間發(fā)生概率；

d、最后通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率；

rp＝s×h24p式3；

其中，rp為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下地質(zhì)災害發(fā)生頻率，s為研究區(qū)內(nèi)地質(zhì)災害的敏感性，h24p為研究區(qū)境內(nèi)雨量監(jiān)測點p頻率降雨下的年最大日降雨量。

所述地理信息系統(tǒng)空間分析功能為arcgis中的arctoolbox-3danalyst-柵格表面功能。

所述步驟b中滑坡的內(nèi)部影響因子為距斷層距離、地形坡度、地層巖性、距水系的距離、海拔高度和地震烈度。

所述步驟b中滑坡的各個內(nèi)部影響因子的無量綱值，通過式4計算；

其中，s為研究區(qū)面積，km²，si為某評價因子i等級所占的面積，km²，n為研究區(qū)內(nèi)滑坡的總數(shù)，ni為某評價因子i等級中滑坡的個數(shù)。

所述步驟c中泥石流的內(nèi)部影響因子為流域坡度、溝壑密度、物源條件和流域高差。

所述步驟c中泥石流的各個內(nèi)部影響因子的無量綱值，通過式5計算；

其中，xi(i,j)為j評價因子中第i個評價單元歸一化后的數(shù)值，x(i,j)為j評價因子中第i個評價單元的數(shù)值，min(x(i,j))為評價因子j中存在的最小值，max(x(i,j))為評價因子j中存在的最大值。

所述步驟d中通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率，是指通過arcgis中的arctoolbox-3danalyst-柵格計算功能編輯式3，計算出p頻率降雨下地質(zhì)災害的發(fā)生頻率。

步驟c中，通過式5計算泥石流的各個內(nèi)部影響因子的無量綱值，不僅計算目標明確，而且其計算結果獲取相當簡便。

步驟d中通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率，是指通過arcgis中的arctoolbox-3danalyst-柵格計算功能編輯式3，計算出p頻率降雨下地質(zhì)災害的發(fā)生頻率，整個計算過程通過程序化處理，其計算結論客觀性更強。

采用本發(fā)明對汶川進行分析，計算不同降雨頻率下滑坡、泥石流的發(fā)生頻率：

汶川縣位于四川盆地阿壩自治州東南部，經(jīng)度和緯度分別為北緯30°45′-31°43′、東經(jīng)102°51′-103°44′，總面積約4084平方千米。區(qū)域內(nèi)龍門山脈、邛崍山脈發(fā)育，西北高東南低，呈傾斜狀，地勢起伏，高差較大，屬于典型的高山峽谷地形。區(qū)域內(nèi)有岷江干流、岷江支流雜谷腦河、漁子溪等流域，發(fā)育茂汶斷層、北川——映秀斷層。

a、通過中國氣象局阿壩藏族羌族自治州氣象局，確定汶川縣境內(nèi)的32個雨量監(jiān)測站點；通過《四川省中小流域暴雨洪水計算手冊》中最大24小時暴雨均值等值線圖，確定汶川縣境內(nèi)32個雨量監(jiān)測點5％、2％、1％頻率降雨下的點雨量均值變差系數(shù)cv；通過cs＝3.5cv適線，確定5％、2％、1％頻率降雨下的離均系數(shù)kp；通過式1計算，確定32個雨量監(jiān)測點5％、2％、1％頻率降雨下的年最大日降雨量h24p；通過畫圖軟件和克里金插值確定汶川縣5％、2％、1％頻率降雨下年最大日降雨量等值線圖；

汶川縣32個雨量監(jiān)測站點最大日降雨量kp及h24p數(shù)值

b、通過地理信息系統(tǒng)空間分析功能，提取滑坡的距斷層距離、地形坡度、地層巖性、距水系的距離、海拔高度、地震烈度；通過滑坡地質(zhì)災害相對密度計算，確定各個內(nèi)部影響因子的無量綱值，通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)實現(xiàn)邏輯回歸分析過程，確定滑坡的各個內(nèi)部影響因子的邏輯回歸系數(shù)；通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)滑坡空間發(fā)生概率；

c、通過地理信息系統(tǒng)空間分析功能，提取泥石流的流域坡度、溝壑密度、物源條件、流域高差，確定泥石流的各個內(nèi)部影響因子的無量綱值；通過統(tǒng)計產(chǎn)品與服務解決方案軟件實現(xiàn)邏輯回歸分析過程，確定泥石流的各個內(nèi)部影響因子的邏輯回歸系數(shù)，通過邏輯回歸模型公式式2計算研究區(qū)境內(nèi)泥石流空間發(fā)生概率；

d、最后通過式3和地理信息系統(tǒng)空間疊置運算確定基于降雨頻率的地質(zhì)災害發(fā)生頻率。

通過室內(nèi)遙感解譯和野外復核，確定汶川縣境內(nèi)共發(fā)育6354處地質(zhì)災害點，其中滑坡6194處，泥石流溝160條。

表明：本發(fā)明將預測時間指標落腳于不同的降雨頻率，空間指標落腳于現(xiàn)有的地質(zhì)災害庫存，明確了地質(zhì)災害發(fā)生的空間概率和時間概率，具有預測結果直觀、預測發(fā)生災害的時間和地點易掌握的特點，能夠有效提高地質(zhì)災害的預測精度和應對災害的主動性。