薛志航,鄧　創(chuàng)，樊曉一,喬建平,王　萌,吳彩燕

(1. 國(guó)網(wǎng)四川省電力公司電力科學(xué)研究院，四川 成都 610072；2. 國(guó)網(wǎng)四川省電力公司電力應(yīng)急中心，四川 成都 610000；3.西南科技大學(xué) 土木工程與建筑學(xué)院，四川 綿陽(yáng) 621010；4.中國(guó)科學(xué)院地表過(guò)程與山地災(zāi)害重點(diǎn)實(shí)驗(yàn)室,中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所, 四川 成都610041)

?

汶川地震擾動(dòng)區(qū)小流域地質(zhì)災(zāi)害風(fēng)險(xiǎn)概率問(wèn)題*

薛志航1,鄧創(chuàng)2，樊曉一3,喬建平4,王萌4,吳彩燕3

(1. 國(guó)網(wǎng)四川省電力公司電力科學(xué)研究院，四川 成都 610072；2. 國(guó)網(wǎng)四川省電力公司電力應(yīng)急中心，四川 成都 610000；3.西南科技大學(xué) 土木工程與建筑學(xué)院，四川 綿陽(yáng) 621010；4.中國(guó)科學(xué)院地表過(guò)程與山地災(zāi)害重點(diǎn)實(shí)驗(yàn)室,中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所, 四川 成都610041)

摘要:汶川地震擾動(dòng)區(qū)地質(zhì)災(zāi)害的特點(diǎn)是滑坡崩塌極為發(fā)育，泥石流物源豐富，災(zāi)害發(fā)生頻率高。以都江堰白沙河小流域?yàn)槔?，在完成小流域風(fēng)險(xiǎn)區(qū)劃的基礎(chǔ)上，采用Logistic 回歸統(tǒng)計(jì)模型分別對(duì)震后小流域中滑坡泥石流發(fā)生的臨界降雨量及時(shí)間概率，不同等級(jí)風(fēng)險(xiǎn)區(qū)(高、較高、中等、較低，共4級(jí)風(fēng)險(xiǎn)區(qū))滑坡泥石流發(fā)生的臨界雨量及空間概率進(jìn)行分析統(tǒng)計(jì)。研究表明當(dāng)日降雨量是觸發(fā)滑坡泥石流的時(shí)間條件，發(fā)生概率為0.7～0.9。大暴雨是觸發(fā)高、較高風(fēng)險(xiǎn)區(qū)滑坡泥石流的空間條件，發(fā)生概率可達(dá)到0.9。發(fā)生風(fēng)險(xiǎn)損失的臨界雨量為大暴雨。

關(guān)鍵詞:汶川地震；擾動(dòng)區(qū)；小流域；地質(zhì)災(zāi)害；臨界雨量；風(fēng)險(xiǎn)概率

汶川8.0級(jí)大地震發(fā)生后，地震擾動(dòng)區(qū)在強(qiáng)震作用下滑坡泥石流的發(fā)育狀況超乎尋常。頻率高、規(guī)模大、分布廣成為基本特點(diǎn)。這些次生地質(zhì)災(zāi)害基本都與降雨相關(guān)，亦即降雨量控制了災(zāi)害的發(fā)生頻率。因?yàn)轱L(fēng)險(xiǎn)評(píng)估中，必須涉及到滑坡泥石流概率發(fā)生的問(wèn)題，否則風(fēng)險(xiǎn)管理將陷入盲目。但到目前為止，對(duì)地震擾動(dòng)區(qū)小流域在降雨條件下地質(zhì)災(zāi)害發(fā)生的概率問(wèn)題[1-20]進(jìn)行的探討比較少。強(qiáng)震作用對(duì)斜坡穩(wěn)定性具有長(zhǎng)期的影響作用，特別是連續(xù)的強(qiáng)降雨使震后的泥石流接連發(fā)生。但是，對(duì)于強(qiáng)震后的滑坡泥石流活動(dòng)性能夠持續(xù)到何時(shí)才能趨于減弱，至今還沒(méi)有實(shí)際而理想的趨勢(shì)預(yù)測(cè)模型。由于缺少有效監(jiān)測(cè)數(shù)據(jù)來(lái)建立隨時(shí)間的后續(xù)降雨量和滑坡泥石流之間的關(guān)系，震后滑坡泥石流的發(fā)育規(guī)律和活動(dòng)特征還不能被人們完全掌握。因此，通過(guò)研究地震擾動(dòng)區(qū)小流域滑坡泥石流的發(fā)生概率，預(yù)測(cè)滑坡泥石流發(fā)育趨勢(shì)具有重要的價(jià)值。筆者已對(duì)小流域滑坡泥石流進(jìn)行了風(fēng)險(xiǎn)區(qū)劃研究，分析了汶川地震擾動(dòng)區(qū)滑坡泥石流發(fā)育的特點(diǎn)，闡述了研究小流域滑坡泥石流的意義，劃分了5級(jí)風(fēng)險(xiǎn)區(qū)[21]。在此基礎(chǔ)上，本文將統(tǒng)計(jì)分析小流域滑坡泥石流風(fēng)險(xiǎn)的時(shí)間概率和空間概率，為下一階段的風(fēng)險(xiǎn)評(píng)估提供依據(jù)[22]。

1回歸統(tǒng)計(jì)模型

對(duì)于滑坡泥石流而言，從統(tǒng)計(jì)學(xué)角度出發(fā)，對(duì)其產(chǎn)生影響的各因子數(shù)據(jù)均可作為自變量，而災(zāi)害發(fā)生與不發(fā)生可以作為分類因變量。由于均不是連續(xù)變量，因此，線性回歸并不適用于推導(dǎo)此類自變量和因變量之間的關(guān)系。所以，通常采用對(duì)數(shù)線性模型(log-linear model)，而Logistic 回歸模型即對(duì)數(shù)線性模型的一種特殊形式。設(shè)p為某事件發(fā)生的概率，取值范圍為[0，1]，(1-p)即為該事件不發(fā)生的概率，將其兩者的比值取自然對(duì)數(shù)ln[p/(1-p)]，以p為因變量，xm為自變量，建立線性Logistic回歸方程：

(1)

式中：α為常數(shù)，βi(i =1，2，…，m)為邏輯回歸系數(shù)。Logistic回歸模型是普通多元線性回歸模型的推廣，但它的誤差項(xiàng)服從二項(xiàng)分布而非正態(tài)分布。

由式(1)可得

(2)

利用式(2)可以對(duì)滑坡泥石流發(fā)生概率進(jìn)行預(yù)測(cè)。但是，對(duì)于滑坡泥石流而言，其難點(diǎn)主要在于自變量的選擇及雨量數(shù)據(jù)的獲取上。

2雨量數(shù)據(jù)選擇

滑坡泥石流分析的雨量數(shù)據(jù)一般是當(dāng)天及前幾天(每天)的雨量記錄，有些地區(qū)也選用小時(shí)甚至分鐘雨量來(lái)進(jìn)行分析，但是考慮到滑坡泥石流災(zāi)害發(fā)生特點(diǎn)及多數(shù)地區(qū)實(shí)際監(jiān)測(cè)情況，當(dāng)天及前幾天的雨量則成為最重要、最通用的分析數(shù)據(jù)。但是地表徑流的產(chǎn)生、水分的蒸發(fā)等過(guò)程的存在，會(huì)使進(jìn)入巖土體的雨量小于實(shí)際記錄雨量，即記錄到的雨量特別是前期降雨不能全部對(duì)滑坡泥石流的發(fā)生產(chǎn)生影響。故本文采用前期有效降雨量的概念。所謂“前期有效降雨量”，是指災(zāi)害發(fā)生

前的降雨過(guò)程中對(duì)斜坡巖土體穩(wěn)定性產(chǎn)生作用的雨量。國(guó)外學(xué)者對(duì)此已作過(guò)相應(yīng)的研究，并提出了計(jì)算進(jìn)入巖土體雨量的經(jīng)驗(yàn)公式[18]：

ra0=kr1+k2r2+…+knrn。

(3)

式中：ra0為災(zāi)害發(fā)生時(shí)的前期有效降雨量，k為有效降雨量系數(shù)(k≤1)，rn為前第n天的降雨量。有效降雨量系數(shù)k表示巖土體對(duì)雨水滯留能力的大小，由區(qū)域內(nèi)巖土體總體性質(zhì)決定，一般取0.84。將前期有效降雨量和當(dāng)日雨量作為自變量，利用雨量和滑坡泥石流的歷史記錄數(shù)據(jù)就可以進(jìn)行滑坡泥石流臨界雨量的分析，得到不同降雨情況下滑坡泥石流的發(fā)生概率。

3時(shí)間概率問(wèn)題

小流域內(nèi)滑坡泥石流誘發(fā)因素主要為降雨條件。臨界降雨量是時(shí)間預(yù)警的關(guān)鍵。降雨量一般可以由前期雨量和當(dāng)日雨量組成。根據(jù)筆者統(tǒng)計(jì)，汶川地震擾動(dòng)區(qū)的滑坡泥石流主要都受到當(dāng)日雨量控制。當(dāng)日降雨誘發(fā)的滑坡泥石流，基本對(duì)應(yīng)了事件發(fā)生的時(shí)間。在此時(shí)間內(nèi)，多大降雨量觸發(fā)滑坡泥石流可能性有多大，這種可能性便稱為風(fēng)險(xiǎn)發(fā)生的時(shí)間概率。

3.1降雨量及災(zāi)害事件統(tǒng)計(jì)

根據(jù)野外實(shí)地調(diào)查和氣象局的資料，表1給出了白沙河小流域滑坡泥石流和當(dāng)日、前期降雨量數(shù)據(jù)。

表1　白沙河小流域流域2009-2013年發(fā)生滑坡泥石流的時(shí)間及雨量

3.2當(dāng)日降雨量誘發(fā)災(zāi)害的概率

依據(jù)災(zāi)害發(fā)生當(dāng)日的降雨量等值線圖進(jìn)行雨量站點(diǎn)數(shù)據(jù)和災(zāi)害發(fā)生地雨量的對(duì)應(yīng)，選取災(zāi)害發(fā)生當(dāng)日降雨量作為災(zāi)害雨量，災(zāi)害發(fā)生前的降雨數(shù)據(jù)作為未發(fā)生災(zāi)害的樣本。利用Logistic回歸模型進(jìn)行分析，由式(2)得到災(zāi)害發(fā)生時(shí)間概率p與當(dāng)日降雨量R0的概率關(guān)系(表2)，即：

(4)

由表2得到當(dāng)日降雨量條件下，滑坡泥石流發(fā)生概率分布圖(圖1)。圖中表明，當(dāng)日降雨量達(dá)到30～40 mm時(shí)，災(zāi)害發(fā)生的概率約為0.5，當(dāng)日降雨量達(dá)到50 mm時(shí)，災(zāi)害發(fā)生的概率為0.7(表3)。

圖1　白沙河小流域當(dāng)日降雨量誘發(fā)災(zāi)害時(shí)間概率分布圖

雨量等級(jí)當(dāng)日降雨量/mm災(zāi)害發(fā)生概率小雨<10<0.1中雨10～24.90.1～0.3大雨25～49.90.3～0.5暴雨50～99.90.5～0.7大暴雨100～249.90.7～0.9特大暴雨>250>0.9

4空間概率問(wèn)題

因?yàn)樾×饔騼?nèi)不同的地形條件、承載體類型各異，風(fēng)險(xiǎn)的程度不同，發(fā)生災(zāi)害的概率也有所差異。所以，統(tǒng)計(jì)小流域滑坡泥石流發(fā)生時(shí)間概率后，還應(yīng)該根據(jù)風(fēng)險(xiǎn)區(qū)劃等級(jí)劃分結(jié)果[]，如：高風(fēng)險(xiǎn)區(qū)、較高風(fēng)險(xiǎn)區(qū)、中等風(fēng)險(xiǎn)區(qū)、較低風(fēng)險(xiǎn)區(qū)和低風(fēng)險(xiǎn)區(qū)(一般不發(fā)生滑坡泥石流災(zāi)害，可不統(tǒng)計(jì))，按不同的風(fēng)險(xiǎn)等級(jí)區(qū)域空間統(tǒng)計(jì)滑坡泥石流的發(fā)生概率，分析不同空間區(qū)域的概率特點(diǎn)。這種概率就稱為風(fēng)險(xiǎn)區(qū)空間概率。

4.1高風(fēng)險(xiǎn)區(qū)概率

以高風(fēng)險(xiǎn)區(qū)的銀洞子溝泥石流和紅色干溝泥石流兩處災(zāi)害點(diǎn)為典型代表，根據(jù)這兩處災(zāi)害發(fā)生的時(shí)間，獲得白沙河小流域高風(fēng)險(xiǎn)區(qū)降雨量分布數(shù)據(jù)(表4)。

根據(jù)高風(fēng)險(xiǎn)區(qū)滑坡泥石流當(dāng)日降雨量數(shù)據(jù)，利用式(2)建立當(dāng)日降雨量條件下發(fā)生滑坡泥石流概率模型，即：

(5)

由此獲得到白沙河小流域高風(fēng)險(xiǎn)區(qū)滑坡泥石流當(dāng)日降雨量與災(zāi)害發(fā)生概率的分布圖(圖2)。結(jié)果表明，當(dāng)日降雨量達(dá)到40～50mm時(shí)，高風(fēng)險(xiǎn)區(qū)災(zāi)害發(fā)生概率為0.5，當(dāng)日降雨量達(dá)到60mm時(shí)，災(zāi)害發(fā)生概率大于0.75(表5)。

表4　高風(fēng)險(xiǎn)區(qū)發(fā)生滑坡泥石流降雨量

圖2　高風(fēng)險(xiǎn)區(qū)當(dāng)日降雨量誘發(fā)災(zāi)害概率分布圖

4.2較高風(fēng)險(xiǎn)區(qū)概率

較高風(fēng)險(xiǎn)包括上坪西村泥石流、觀鳳溝泥石流、高原村4組老干溝泥石流、大干溝泥石流和深溪溝泥石流。依據(jù)這5處災(zāi)害點(diǎn)發(fā)生的時(shí)間，獲得白沙河小流域較高風(fēng)險(xiǎn)區(qū)滑坡泥石流降雨量分布數(shù)據(jù)(表6)。

表5　高風(fēng)險(xiǎn)區(qū)當(dāng)日降雨誘發(fā)災(zāi)害概率統(tǒng)計(jì)表

利用式(2)建立較高風(fēng)險(xiǎn)區(qū)滑坡泥石流與當(dāng)日降雨量概率統(tǒng)計(jì)模型，即：

(6)

由此獲得較高風(fēng)險(xiǎn)區(qū)滑坡泥石流當(dāng)日降雨量與災(zāi)害發(fā)生概率的分布圖(圖3)。圖中表明，當(dāng)日降雨量達(dá)到50mm時(shí)，災(zāi)害發(fā)生的概率約為0.5，當(dāng)日降雨量達(dá)到75mm時(shí)，災(zāi)害發(fā)生的概率大于0.75(表7)。

表6　白沙河流域較高風(fēng)險(xiǎn)區(qū)發(fā)生滑坡泥石流降雨量

圖3　較高風(fēng)險(xiǎn)區(qū)當(dāng)日降雨量誘發(fā)災(zāi)害概率分布圖

4.3中等風(fēng)險(xiǎn)區(qū)概率

中等風(fēng)險(xiǎn)區(qū)僅有鍋圈巖泥石流屬于中等風(fēng)險(xiǎn)區(qū)。依據(jù)鍋圈巖泥石流災(zāi)害點(diǎn)發(fā)生的時(shí)間，獲得白沙河流域中等風(fēng)險(xiǎn)區(qū)滑坡泥石流降雨量分布數(shù)據(jù)(表8)。

表7　較高風(fēng)險(xiǎn)區(qū)當(dāng)日降雨誘發(fā)災(zāi)害概率統(tǒng)計(jì)表

利用式(2)建立中等風(fēng)險(xiǎn)區(qū)滑坡泥石流與當(dāng)日降雨量概率統(tǒng)計(jì)模型，即：

(7)

由此得到白沙河流域中等風(fēng)險(xiǎn)區(qū)滑坡泥石流當(dāng)日降雨量與災(zāi)害發(fā)生概率的分布圖(圖4)。圖中表明，當(dāng)日降雨量達(dá)到70 mm時(shí)，災(zāi)害發(fā)生的概率約為0.5，當(dāng)日降雨量達(dá)到130 mm時(shí)，災(zāi)害發(fā)生的概率大于0.75(表9)。

圖4　中等風(fēng)險(xiǎn)區(qū)當(dāng)日降雨量誘發(fā)災(zāi)害概率分布圖

4.4較低風(fēng)險(xiǎn)區(qū)概率

較低風(fēng)險(xiǎn)區(qū)的滴水巖滑坡崩塌、峰子坨滑坡崩塌、丁木樹(shù)溝泥石流、泡桐槽泥石流和甜竹坪滑坡屬于較低風(fēng)險(xiǎn)區(qū)。依據(jù)災(zāi)害點(diǎn)發(fā)生的時(shí)間，獲得白沙河流域較低風(fēng)險(xiǎn)區(qū)滑坡、泥石流降雨量分布數(shù)據(jù)(表10)。

表9　中等風(fēng)險(xiǎn)區(qū)當(dāng)日降雨誘發(fā)災(zāi)害概率統(tǒng)計(jì)表

表10　白沙河流域較低風(fēng)險(xiǎn)區(qū)發(fā)生滑坡泥石流降雨量

利用式(2)建立較低風(fēng)險(xiǎn)區(qū)滑坡泥石流與當(dāng)日降雨量概率統(tǒng)計(jì)模型，即：

(8)

由此得到白沙河流域較低風(fēng)險(xiǎn)區(qū)滑坡泥石流當(dāng)日降雨量與災(zāi)害發(fā)生概率的分布圖(圖5)。圖中表明，當(dāng)日降雨量達(dá)到100mm時(shí)，災(zāi)害發(fā)生的概率約為0.5(表11)。

圖5　較低風(fēng)險(xiǎn)區(qū)當(dāng)日降雨量誘發(fā)災(zāi)害概率分布圖

雨量等級(jí)當(dāng)日降雨量/mm發(fā)生概率小雨<10<0.03中雨10～24.90.03～0.10大雨25～49.90.10～0.26暴雨50～99.90.26～0.47大暴雨100～249.90.47～0.67特大暴雨>2500.67～0.72

4.5全流域空間概率統(tǒng)計(jì)分析

參照?qǐng)D2～圖5所示的結(jié)果，建立全流域高風(fēng)險(xiǎn)區(qū)、較高風(fēng)險(xiǎn)區(qū)、中等風(fēng)險(xiǎn)區(qū)和較低風(fēng)險(xiǎn)區(qū)4個(gè)區(qū)域(空間)與當(dāng)日(時(shí)間)降雨量誘發(fā)災(zāi)害發(fā)生的概率關(guān)系(圖6)，可見(jiàn)高風(fēng)險(xiǎn)區(qū)和較高風(fēng)險(xiǎn)區(qū)的滑坡泥石流發(fā)生概率明顯大于其他風(fēng)險(xiǎn)區(qū)，并且隨風(fēng)險(xiǎn)區(qū)的等級(jí)減小，滑坡泥石流的發(fā)生概率降低。高風(fēng)險(xiǎn)區(qū)的滑坡泥石流發(fā)生的概率與較高風(fēng)險(xiǎn)區(qū)較為接近。

流域內(nèi)在小雨條件下，只有高風(fēng)險(xiǎn)區(qū)滑坡泥石流發(fā)生概率大于0.1，其他風(fēng)險(xiǎn)區(qū)滑坡泥石流發(fā)生概率都小于0.1。在中雨條件下，高風(fēng)險(xiǎn)區(qū)滑坡泥石流發(fā)生概率為0.11～0.24，大于其他風(fēng)險(xiǎn)等級(jí)區(qū)的滑坡泥石流概率。在大雨條件下，高風(fēng)險(xiǎn)區(qū)、較高風(fēng)險(xiǎn)區(qū)、中等風(fēng)險(xiǎn)區(qū)與較低風(fēng)險(xiǎn)區(qū)滑坡泥石流發(fā)生最大概率分別為0.61、0.52、0.38和0.26。在暴雨條件下，高風(fēng)險(xiǎn)區(qū)、較高風(fēng)險(xiǎn)區(qū)滑坡泥石流發(fā)生最大概率大于0.80。在特大暴雨條件下，高風(fēng)險(xiǎn)區(qū)、較高風(fēng)險(xiǎn)區(qū)和較高風(fēng)險(xiǎn)區(qū)滑坡泥石流發(fā)生最大概率大于0.80(表12)。

圖6　當(dāng)日降雨條件下白沙河全流域不同風(fēng)險(xiǎn)區(qū)滑坡泥石流概率分布圖

5結(jié)語(yǔ)

(1)地震擾動(dòng)區(qū)小流域滑坡泥石流風(fēng)險(xiǎn)與其他區(qū)域風(fēng)險(xiǎn)相同，都應(yīng)該包括風(fēng)險(xiǎn)的概率問(wèn)題。發(fā)生滑坡泥石流概率越大，意味著風(fēng)險(xiǎn)的存在可能性越高，反之亦然。因此進(jìn)行風(fēng)險(xiǎn)評(píng)估時(shí)，必須回答評(píng)估區(qū)風(fēng)險(xiǎn)概率的類型和大小。

(2)白沙河小流域滑坡泥石流風(fēng)險(xiǎn)概率包括時(shí)間概率和空間概率兩種類型。其中時(shí)間概率主要與當(dāng)日降雨量相關(guān)，當(dāng)日降雨量達(dá)到暴雨時(shí)存在發(fā)生滑坡泥石流風(fēng)險(xiǎn)概率。當(dāng)日降雨量達(dá)到或超過(guò)大暴雨和特大暴雨時(shí)，滑坡泥石流風(fēng)險(xiǎn)概率迅速提高?？臻g概率因風(fēng)險(xiǎn)等級(jí)區(qū)域而異，高、較高風(fēng)險(xiǎn)區(qū)滑坡泥石流受臨界降雨量控制，發(fā)生風(fēng)險(xiǎn)損失的概率最高。

(3)白沙河小流域可能出現(xiàn)滑坡泥石流風(fēng)險(xiǎn)損失的臨界降雨量為大暴雨。降雨量達(dá)到大暴雨時(shí)，風(fēng)險(xiǎn)概率曲線基本穩(wěn)定在0.9水平。建議采用大暴雨作為風(fēng)險(xiǎn)預(yù)警臨界指標(biāo)。

(4)研究區(qū)的地震滑坡泥石流應(yīng)具有衰減性。隨時(shí)間推移，發(fā)生頻率和規(guī)模都將逐漸衰減，因此風(fēng)險(xiǎn)概率也應(yīng)該出現(xiàn)變化。此類問(wèn)題還應(yīng)該通過(guò)長(zhǎng)期觀測(cè)，積累有價(jià)值的資料數(shù)據(jù)進(jìn)一步研究，得出科學(xué)結(jié)論。

表12　當(dāng)日降雨條件下白沙河全流域不同風(fēng)險(xiǎn)區(qū)誘發(fā)災(zāi)害概率匯總表

參考文獻(xiàn):

[1]唐川.汶川地震區(qū)暴雨滑坡泥石流活動(dòng)趨勢(shì)預(yù)測(cè)[J].山地學(xué)報(bào)，2010，28(3)：341-349.

[2]喬建平,王萌,石莉莉.區(qū)域滑坡風(fēng)險(xiǎn)評(píng)估中的風(fēng)險(xiǎn)區(qū)劃與概率分析[J].自然災(zāi)害學(xué)報(bào),2012,21(2): 51-56.

[3]侯圣山,李昂,周平根.四川雅安市雨城區(qū)地質(zhì)災(zāi)害預(yù)警系統(tǒng)研究[J].地學(xué)前緣,2007,14(6):160-165.

[4]彭貴芬,劉盈曦.基于模糊信息概率區(qū)間數(shù)的突發(fā)地質(zhì)災(zāi)害降水因子評(píng)價(jià)研究[J].災(zāi)害學(xué),2015,30(1)：1-4.

[5]萬(wàn)石云,趙寧坤,段瑋.云南省滑坡泥石流與強(qiáng)降水的相關(guān)性和多時(shí)間尺度特征研究[J].災(zāi)害學(xué),2015,30(1)：45-50.

[6]雷發(fā)洪,胡凱衡,胡云華,等.汶川地震災(zāi)區(qū)震后泥石流激發(fā)雨量研究[J].災(zāi)害學(xué),2014,29(2)：199-203.

[7]Xu W B, Jing S C, Yu W J, et al. A comparison between Bayes discriminant analysis and logisticregression for prediction of debris flow in southwest Sichuan, China[J]. Geomorphology, 2013, 201: 45-51.

[8]Nikolopoulos I E, Crema S A B, Marchi L, et al. Impact of uncertainty in rainfall estimation on the identification of rainfall thresholds for debris flow occurrence[J]. Geomorphology, 2014, 221:286-297.

[9]Chang Ming, Tang Chuan, Zhang Dandan,et al. Debris flow susceptibility assessment using a probabilistic approach: A case study in the longchi area, Sichuan province, China[J]. Journal of Mountain Science, 2014, 11(4): 1001-1014.

[10]喬建平,楊宗佶,田宏嶺. 降雨滑坡預(yù)警的概率分析方法[J].工程地質(zhì)學(xué)報(bào),2009,17(3):43-348.

[12]Glade T. Applying probability determination to refine landslide-triggering rainfall thresholds using empirical antecedent daily rainfall model [J]. Pure and Applied Geophysics, 2000,157:1059-1079.

[13]Glade T. Modelling landslide-triggering rainfalls in different regions of new zealand the soil water status model[J]. Zeitschrift fur Geomorphologie N E, 2000, 122: 63-84.

[14]Guidicini G. Tentative correlation between rainfall and landslides in a humid tropical environment[J]. Bulletin of IAEC,1997,16:13-20.

[15]Saito H. Relationship Between the initiation of a shallow landslide and rainfall intensity duration thresholds in Japan[J]. Geomorphology, 2010, 118: 167-175.

[16]Finlay P J. The relationship between the probability of landslide occurrence and rainfall[J]. Canadian geotechnical,1997,34(6):811-824.

[17]謝守益,張年學(xué),許兵.長(zhǎng)江三峽庫(kù)區(qū)典型滑坡降雨誘發(fā)的概率分析[J].工程地質(zhì)學(xué)報(bào),1995,3(2):60-69.

[18]陳劍,楊志法,李曉.三峽庫(kù)區(qū)滑坡發(fā)生概率與降水條件的關(guān)系[J].巖石力學(xué)與工程學(xué)報(bào),2005,24(17)：3052-3056.

[19]陳靜靜,姚蓉,文強(qiáng),等.湖南省降雨型地質(zhì)災(zāi)害致災(zāi)雨量閾值分析[J].災(zāi)害學(xué),2014,29(2)：42-47.

[20]陳麗霞,殷坤龍,劉長(zhǎng)春.降雨重現(xiàn)期及其用于滑坡概率分析的探討[J].工程地質(zhì)學(xué)報(bào), 2012, 20(5):745-750

[21]喬建平,王萌,吳彩燕.汶川地震擾動(dòng)區(qū)小流域滑坡泥石流風(fēng)險(xiǎn)區(qū)劃[J].災(zāi)害學(xué)，2016，31(2)：1-5.

[22]喬建平,王萌.滑坡風(fēng)險(xiǎn)分類及層次鏈[J].工程地質(zhì)學(xué)報(bào),2010,18(1)：84-90.

*收稿日期：2015-12-21修回日期：2016-03-16

基金項(xiàng)目：國(guó)家電網(wǎng)公司科技項(xiàng)目(521999150031)；國(guó)家自然科學(xué)基金項(xiàng)目(41272297)；科技部重點(diǎn)國(guó)際合作項(xiàng)目(2013DFA21720)；國(guó)家青年科學(xué)基金(41301592); 中國(guó)科學(xué)院山地災(zāi)害與地表過(guò)程重點(diǎn)實(shí)驗(yàn)室自主支持基金項(xiàng)目資助(2013年度)；長(zhǎng)江科學(xué)院開(kāi)放研究基金資助項(xiàng)目(CKWV2015228/KY)

第一作者簡(jiǎn)介：薛志航(1987-)，男，四川德陽(yáng)人，碩士，工程師，主要從事電網(wǎng)防災(zāi)減災(zāi).E-mail：xuezhihang2@163.com 通訊作者：樊曉一(1974-)，男，四川成都人，博士，教授，主要從事巖土工程及地質(zhì)災(zāi)害方面的教學(xué)與研究工作. E-mail：xyfan1003@126.com

中圖分類號(hào)：X43；P642

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-811X(2016)03-0039-07

doi:10.3969/j.issn.1000-811X.2016.03.007

Study on Geological Hazards of Small Watershed Risk Probability in Wenchuan Eathquake Disturbance Area

XUE Zhihang1, DENG Chuang2, FAN Xiaoyi3, QIAO Jianping4, WANG Meng4, WU Caiyan3

(1.StateGridSichuanElectricPowerResearchInstitute,Chengdu610072,China;2.PowerEmergencyCenterofStateGridPowerCompanyofSichuanProvince,Chengdu610000,China;3.SchoolofCivilEngineeringandArchitecture,SouthwestUniversityofScienceandTechnology,Mianyang62010,China; 4.KeyLaboratoryofMountainHazardsandEarthSurfaceProcesses,InstituteofMountainHazardsandEnvironment,CAS,Chengdu610041,China)

Abstract:The characteristics of geological hazards in Wenchuan earthquake disturbance area are numerous landslides and rock falls, rich debris flow materials, and high occurrence frequency of hazards. We selected Dujiangyan City as the study area. Based on small watershed risk zonation, logistic regression model was adopted to make a statistics analysis on critical rainfall and time probability aiming at landslides and debris flows of small watershed after earthquake and them in different grade risk zones. The study results showed that intraday rainfall was the time condition to trigger landslide and debris flow with probability of 0.7 to 0.9. Heavy rainfall was spatial contidion to trigger landslide and debris flow in higher and high risk areas with occurrence probability of 0.9. Critical rainfall to make risk loss was heavy rainfall.

Key words:Wenchuan earthquake; small watershed; hazard; critical rainfall; probability

薛志航,鄧創(chuàng)，樊曉一，等. 汶川地震擾動(dòng)區(qū)小流域地質(zhì)災(zāi)害風(fēng)險(xiǎn)概率問(wèn)題[J]. 災(zāi)害學(xué)，2016，31(3)：39-45. [XUE Zhihang,DENG Chuang, FAN Xiaoyi, et al. Geological Hazards of Small Watershed Risk Probability in Wenchuan Eathquake Disturbance Area[J].Journal of Catastrophology，2016，31(3)：39-45.]