1、. 基于概率數(shù)學(xué)方法與GIS的泥石流敏感性分析及評(píng)價(jià)以汶川縣為例 基于概率數(shù)學(xué)方法與GIS的泥石流敏感性分析及評(píng)價(jià)以汶川縣為例 夏晨皓，朱 靜，常 鳴， 宇 (理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 610059) 摘 要：以汶川縣為研究區(qū)，通過對(duì)震后的高分辨遙感影像的遙感解譯，獲得了145條泥石流溝及流域中的崩滑物源的分布圖。選擇了坡度、坡向、地層巖性、地震峰值加速度、距水系距離、距斷層距離作為評(píng)價(jià)因子，利用概率數(shù)學(xué)方法與GIS的空間分析功能，獲得145條泥石流溝在上述6個(gè)評(píng)價(jià)因子上的分值；再利用層次分析法計(jì)算得到各評(píng)價(jià)因子權(quán)重，建立泥石流敏感性評(píng)價(jià)模型，得到了研究區(qū)泥石流溝敏

2、感性分級(jí)圖。結(jié)果說明： 研究區(qū)145條泥石流溝中，有60條泥石流溝屬于高敏感性，43條泥石流溝為中敏感性，其余42條泥石流溝敏感性較低，研究成果對(duì)于認(rèn)識(shí)該地區(qū)泥石流災(zāi)害現(xiàn)狀及開展趨勢(shì)與防災(zāi)減災(zāi)有一定的參考意義。 關(guān)鍵詞：泥石流；層次分析法；概率數(shù)學(xué)方法;敏感性；汶川縣；GIS 1 研究背景 2008年5月12日在省汶川縣發(fā)生了8.0級(jí)特震， 汶川縣作為震中心， 震后次生山地災(zāi)害頻發(fā)， 有學(xué)者研究說明1-3， 汶川地震后， 受強(qiáng)震影響的重災(zāi)區(qū)的地質(zhì)災(zāi)害活潑期將持續(xù)20 a以上。 由于汶川縣既處于地質(zhì)環(huán)境脆弱的高地區(qū)， 又受512地震影響， 境泥石流災(zāi)害頻發(fā)， 它時(shí)刻危害著當(dāng)?shù)厝嗣竦纳a(chǎn)生活平安，

3、 因此,對(duì)汶川縣進(jìn)展整體的區(qū)域泥石流敏感性評(píng)價(jià)對(duì)認(rèn)識(shí)該地區(qū)災(zāi)害現(xiàn)狀及開展趨勢(shì)與防災(zāi)減災(zāi)有積極的作用， 對(duì)于保障該地區(qū)的震后社會(huì)經(jīng)濟(jì)可持續(xù)開展具有現(xiàn)實(shí)意義。 泥石流敏感性評(píng)價(jià)是對(duì)預(yù)測(cè)評(píng)價(jià)*區(qū)域未來*段時(shí)間導(dǎo)致泥石流爆發(fā)的各致災(zāi)因子等各種敏感性指標(biāo)的變化情況進(jìn)展分析，估算它們的概率分布，完成它們的敏感性評(píng)價(jià)圖。目前3S技術(shù)在地質(zhì)災(zāi)害領(lǐng)域中廣泛運(yùn)用，成為了泥石流(滑坡)的敏感性分析評(píng)價(jià)中的有力工具，國(guó)外專家們也在此根底上運(yùn)用了許多敏感性分析模型。 Ayalew等42004年在日本Kakuda-Yahiko山區(qū)利用邏輯回歸模型開展了基于GIS的滑坡敏感性分析研究；唐川52005年在怒江流域利用條件概率

4、模型與GIS結(jié)合繪制了當(dāng)?shù)啬嗍髅舾行詫ｎ}地圖；Chang等62007年在臺(tái)北地區(qū)利用人工神經(jīng)網(wǎng)絡(luò)法為當(dāng)?shù)?71條泥石流溝作了定量危險(xiǎn)性評(píng)價(jià)研究，評(píng)價(jià)結(jié)果也根本準(zhǔn)確。洋等7在2013年利用信息量模型與層次分析法對(duì)省龍池地區(qū)48條泥石流溝進(jìn)展敏感性評(píng)價(jià),從而確定了該地區(qū)需要重點(diǎn)監(jiān)測(cè)防治的泥石流溝道。向靈芝等8在2015年將流域演化理論運(yùn)用到泥石流敏感性研究中，對(duì)汶川縣都汶公路和省道303沿線60條泥石流溝進(jìn)展敏感性評(píng)價(jià)，結(jié)果根本反映了地貌演化與泥石流開展的關(guān)系。 基于上述研究，考慮到研究區(qū)汶川縣是512地震強(qiáng)震區(qū)，震后誘發(fā)的大量的崩塌滑坡在強(qiáng)降雨情況下就可能很快轉(zhuǎn)變成泥石流或參與泥石流運(yùn)動(dòng)，筆者

5、認(rèn)為汶川縣泥石流的敏感性與汶川縣的同震滑坡可能有較強(qiáng)聯(lián)系。因此,本文以汶川縣為研究區(qū)，通過對(duì)其震后的高分辨遙感影像的遙感解譯，獲得145條泥石流溝及流域中的崩滑物源的分布圖。結(jié)合研究區(qū)特點(diǎn)，選擇適宜評(píng)價(jià)因子，采用概率數(shù)學(xué)方法與GIS空間分析功能得到每條泥石流溝各評(píng)價(jià)因子的分值;又利用層次分析法確定各評(píng)價(jià)因子權(quán)重，最終建立泥石流敏感性評(píng)價(jià)模型，使評(píng)價(jià)結(jié)果更加客觀可靠。 2 研究區(qū)概況 研究區(qū)汶川縣位于盆地西北部、阿壩藏族羌族自治州東南部，地理位置為3045N3143N,10251 E103 44 E之間,全縣面積4 084 km2，常住人口為101 533人(2013年末)。汶川縣以高地形為主，

6、地形坡度大，境溝壑縱橫，切割強(qiáng)烈，水系興旺，河流眾多，有大小支流達(dá)百余條。汶川縣降雨較多且集中、干濕季明顯，雨量528.71 332.2 mm9，全縣處于九頂山新華夏構(gòu)造帶,巖體地質(zhì)空間特征變化復(fù)雜，巖性差異變化較大，地層發(fā)育較為完整,只是除了奧系、志留系大局部缺失。斷層、褶皺較為發(fā)育,茂縣汶川斷裂帶和北川映秀斷裂帶2條主要大斷裂帶呈北東南西方向穿過汶川縣。其中映秀北川斷裂是512地震的發(fā)震斷裂。 汶川縣在2008年汶川地震以前就是地質(zhì)災(zāi)害高易發(fā)區(qū)域，特別是在岷江、漁子溪、雜古腦河沿岸斜坡地質(zhì)災(zāi)害嚴(yán)重，其中泥石流又是汶川縣最主要的地質(zhì)災(zāi)害類型10。根據(jù)唐邦興等11的野外調(diào)查，境有泥石流活動(dòng)歷史

7、的溝多達(dá)98條。由于512汶川地震為泥石流提供了豐富的松散固體物質(zhì)，同時(shí)根據(jù)周偉等12研究，地震后震區(qū)泥石流起動(dòng)的臨界雨量也較地震前有所降低，使得在強(qiáng)降雨情況下泥石流爆發(fā)的可能性大增。在512震后，汶川縣泥石流災(zāi)害頻發(fā)。2008年9月2226日，汶川縣進(jìn)入震后第1個(gè)雨季里就有大圍的泥石流活動(dòng)，2010年 8月13日的強(qiáng)降雨讓汶川縣多個(gè)鄉(xiāng)鎮(zhèn)泥石流群發(fā)。2013年7月10日汶川縣遭受了自地震以來最大規(guī)模的一次群發(fā)性泥石流災(zāi)害，全縣90%以上的村寨不同程度受災(zāi)。因此,可以看出汶川縣泥石流溝不僅分布密集，危害嚴(yán)重,且進(jìn)入到了一個(gè)新的活潑期。 3 數(shù)據(jù)獲取 研究區(qū)等高線與水系的矢量數(shù)據(jù)主要來源于省測(cè)繪地

8、理信息局的150 000地形圖，斷裂帶及巖性的矢量數(shù)據(jù)主要來源于中國(guó)地質(zhì)調(diào)查局的1100 000地質(zhì)圖，地震峰值加速度值來源于省地震局，研究區(qū)DEM數(shù)據(jù)通過150 000地形圖中的等高線在ArcGIS軟件中插值生成，網(wǎng)格大小為10 m10 m。 研究區(qū)崩滑體空間屬性數(shù)據(jù)庫(kù)主要利用工程組對(duì)汶川縣城、映秀鎮(zhèn)、岷江及漁子溪沿岸等重點(diǎn)區(qū)域的遙感解譯成果和戴福初等132011年對(duì)汶川縣全縣震后次生災(zāi)害遙感解譯成果建立的。其中采用的遙感影像數(shù)據(jù)主要有SPOT5衛(wèi)星影像(分辨率2.5 m)、WorldView-2衛(wèi)星影像(分辨率0.5 m)和美國(guó) IKNOS 衛(wèi)星影像(分辨率1 m)。 4 研究方法 4.1

9、 評(píng)價(jià)因子選取 控制和影響泥石流敏感性程度的因子有多種。因?yàn)殂氪h是地震強(qiáng)震區(qū)，地震對(duì)汶川縣地表及山坡穩(wěn)定性破壞強(qiáng)烈，導(dǎo)致崩塌、滑坡等地質(zhì)災(zāi)害的發(fā)育，為泥石流提供了豐富的松散固體物質(zhì)，增強(qiáng)了泥石流的形成條件，所以泥石流溝流域的崩滑體敏感性是影響研究區(qū)泥石流敏感性的重要因子，因此,選擇了泥石流溝流域及其崩滑體都適用的評(píng)價(jià)因子:坡度、坡向、地層巖性、地震峰值加速度、距水系距離、距斷層距離這6個(gè)因子。 (1) 當(dāng)斜坡坡度大于自然休止角并且沒有足夠大的抗滑力時(shí)，斜坡就會(huì)發(fā)生變形破壞14。通常來說，坡度越大，土體自重力沿滑動(dòng)面向下的分力也將會(huì)越大，這將使滑坡越不穩(wěn)定，但是坡度也不宜過大，因?yàn)檫^大的坡度也

10、不利于松散固體物質(zhì)的儲(chǔ)存堆積。因此,坡度一方面將會(huì)影響松散固體物質(zhì)的分布和聚集;另一方面還會(huì)影響斜坡應(yīng)力的分布，并對(duì)坡外表的地表水徑流和坡體部地下水的供應(yīng)和排泄起控制作用，同時(shí)較大的坡度為斜坡兩側(cè)崩滑體的能量轉(zhuǎn)換提供客觀條件。因此,坡體的坡度陡緩將間接影響泥石流的規(guī)模大小。 (2) 坡向?qū)τ诨聻?zāi)害以及泥石流災(zāi)害來說是一個(gè)較為重要的地形特征，坡向主要影響局部的小氣候，而小氣候又通過太陽(yáng)輻射、溫度、水分蒸發(fā)等渠道影響坡面上自然地理諸要素。研究區(qū)海拔較高，接收日照時(shí)間更長(zhǎng)，并且研究區(qū)地形起伏較大，相對(duì)高差較大，坡度也較大，導(dǎo)致了各坡向斜坡的自然地理諸要素的規(guī)律性差異更為巨大。同時(shí)，不同坡向造成巖石

11、風(fēng)化程度也有所不同，影響巖層的破碎程度，進(jìn)而可能在滑坡、崩塌過程中產(chǎn)生大量源物質(zhì)。 (3) 地層巖性是產(chǎn)生滑坡地質(zhì)災(zāi)害的物質(zhì)根底，并且在一定的程度上影響泥石流的性質(zhì)。泥石流通常在千枚巖、板巖、砂巖和花崗巖風(fēng)化殼地區(qū)分布較為集中，比方研究區(qū)里志留系的板巖、千枚巖因?yàn)閿嗔褞Ъ暗卣鹩绊懀瑤r石破碎、風(fēng)化強(qiáng)烈,為泥石流的爆發(fā)提供了豐富的松散物源,也為泥石流裹挾巨石長(zhǎng)距離運(yùn)移提供必要條件 15 。由各類松散土體及固結(jié)較差的軟質(zhì)巖石組成的斜坡，在強(qiáng)降雨后隨著抗剪強(qiáng)度的下降穩(wěn)定性也會(huì)顯著降低。研究區(qū)圍共有7個(gè)地質(zhì)年代的地層出露，分布最廣的3種地層是由花崗巖和閃長(zhǎng)巖組成的二疊紀(jì)、砂巖和千枚巖組成的三疊紀(jì)以及板巖

12、、千枚巖組成的志留系，它們構(gòu)造破碎、風(fēng)化嚴(yán)重，且受強(qiáng)震影響產(chǎn)生大量松散物源。 (4) 河流水系的沖蝕影響滑坡災(zāi)害的發(fā)生和泥石流災(zāi)害的規(guī)模，因?yàn)楹恿魉禌_蝕切割時(shí)會(huì)產(chǎn)生大量的臨空面，造成了許多的滑移控制面外露，降低了土體斜坡的穩(wěn)定性，進(jìn)而導(dǎo)致了滑坡，而這些滑坡又將為泥石流的發(fā)生提供豐富的松散物源。 (5) 在地質(zhì)構(gòu)造運(yùn)動(dòng)影響下斷裂帶附近巖層破碎嚴(yán)重，特別是在有差異性較大的升降運(yùn)動(dòng)時(shí)，巖層經(jīng)過擠壓破壞而穩(wěn)定性會(huì)受損，更易發(fā)生崩滑災(zāi)害，為泥石流發(fā)生提供豐富的松散物源;另外根據(jù)以前對(duì)集集地震同震滑坡的研究說明，地震崩滑體的空間分布與它的發(fā)震斷層關(guān)系密切16，汶川地震和1999年地區(qū)集集地震一樣，屬于逆

13、斷層型地震，汶川地震的發(fā)震斷層為映秀北川斷裂帶。 (6)因?yàn)榈卣鸹顒?dòng)導(dǎo)致泥石流災(zāi)害加劇甚至直接使泥石流爆發(fā)的事件較多，前人也已有研究17，地震活動(dòng)一方面直接使山體劇烈破壞，產(chǎn)生大量同震崩滑體，提供豐富的松散物源;另一方面，地震的巨大動(dòng)力還可以導(dǎo)致巖體強(qiáng)度降低，土體強(qiáng)度大量下降，最終導(dǎo)致斜坡失穩(wěn)參與泥石流活動(dòng)。因此地震峰值加速度可以選為泥石流敏感性因子。 4.2 各評(píng)價(jià)因子概率取值 確定評(píng)價(jià)因子后需要將這些因子進(jìn)展量化，以保證泥石流敏感性評(píng)價(jià)的科學(xué)性。采用概率數(shù)學(xué)方法來計(jì)算各個(gè)評(píng)價(jià)因子在不同的等級(jí)區(qū)間圍影響泥石流發(fā)生的概率，相對(duì)提高了敏感性評(píng)價(jià)的科學(xué)性18。具體采用式(1)進(jìn)展計(jì)算。 (1) 式

14、中：p為評(píng)價(jià)因子在不同等級(jí)區(qū)間圍概率綜合值；p1為評(píng)價(jià)因子不同等級(jí)區(qū)間圍崩滑體面積占該區(qū)間圍研究區(qū)總面積(泥石流溝流域面積)的比例，是條件概率；p2為整個(gè)研究區(qū)所有崩滑體面積占整個(gè)研究區(qū)流域總面積比例，是先驗(yàn)概率。 根據(jù)式(1):當(dāng)p1p2時(shí)，說明該評(píng)價(jià)因子在這個(gè)分級(jí)條件區(qū)間的條件概率大于先驗(yàn)概率，對(duì)泥石流的發(fā)生有明顯的推動(dòng)作用;當(dāng)p1p2時(shí)，說明該評(píng)價(jià)因子在這個(gè)分級(jí)條件區(qū)間的條件概率小于先驗(yàn)概率，對(duì)泥石流的發(fā)生推動(dòng)作用不明顯。 利用研究區(qū)DEM數(shù)據(jù)，使用ArcGIS軟件提取研究區(qū)地形坡度、坡向特征，并根據(jù)前人對(duì)汶川地震同震滑坡統(tǒng)計(jì)編錄及敏感性研究19 ，將距水系距離按200 m為間隔，距斷層

15、距離按15 km為間隔，分為4類，建立緩沖區(qū)，將各評(píng)價(jià)因子按DEM精度10 m10 m柵格化處理并進(jìn)展空間疊加分析參與泥石流敏感性評(píng)價(jià)分析，得到研究區(qū)各因子分級(jí)見圖1。 圖1 研究區(qū)各因子分級(jí)Fig.1 Zoning of the study area according to rating of factors 利用公式(1)劃分出各評(píng)價(jià)因子下的分級(jí)區(qū)間,并計(jì)算評(píng)價(jià)因子在不同等級(jí)區(qū)間圍的概率綜合值P,見表1。 4.3 敏感性評(píng)價(jià)因子綜合賦值及統(tǒng)計(jì) 根據(jù)上述對(duì)6個(gè)評(píng)價(jià)因子在不同等級(jí)區(qū)間圍里的概率綜合值統(tǒng)計(jì)分析整合，將每個(gè)敏感性評(píng)價(jià)因子劃分為不同分級(jí)圍區(qū)間并分別進(jìn)展權(quán)重賦值。其中，數(shù)值的大小表示

16、該因子在研究區(qū)區(qū)域的影響程度,具體見表2。參照表2里的各因子分級(jí)圍及綜合賦值情況,結(jié)合研究區(qū)145條泥石流溝單溝在各評(píng)價(jià)因子上的空間面積分布情況，在ArcGIS軟件中進(jìn)展統(tǒng)計(jì)量化計(jì)算，得到每條泥石流溝在6個(gè)評(píng)價(jià)因子上的敏感性量化分值。 表1 研究區(qū)各因子分級(jí)及泥石流概率綜合值統(tǒng)計(jì)表Table 1 Rating of factors and statistics of prehensive values of probability of debris flow in the study area因子分級(jí)滑坡體面積/km2流域面積/km2P1P2P坡度坡向地震峰值加速度地層巖性距斷層距離距水系距

17、離0,12)3.7054.210.0680.095-0.30212,24)9.86168.840.0580.095-0.40924,32)17.08295.190.0580.095-0.41532,40)37.76460.990.0820.095-0.1540,48)44.56374.700.1190.0950.22248,86)36.38219.680.1660.0950.471N(022.5,337.5360)16.01176.140.0910.095-0.047NE(22.567.5)15.72166.930.0940.095-0.01E(67.5112.5)19.95201.910.0

18、990.0950.043SE(112.5157.5)28.25231.790.1220.0950.245S(157.5202.5)23.42214.340.1090.0950.144SW(202.5247.5)17.04190.030.090.095-0.062W(247.5292.5)14.35200.570.0720.095-0.266NW(292.5337.5)14.60191.920.0760.095-0.2160.25g52.13571.820.0910.095-0.044侏羅紀(jì)砂巖1.3611.110.1220.0950.245三疊紀(jì)砂巖和千枚巖4.31176.390.0240.0

19、95-0.432二疊紀(jì)閃長(zhǎng)巖122.58997.390.1230.0950.251石炭系灰?guī)r1.3914.590.0960.0950.006泥盆紀(jì)石灰?guī)r6.45165.700.0390.095-0.614志留紀(jì)板巖9.94179.470.0550.095-0.441震旦砂巖和白云巖3.3128.950.1140.0950.1860,15)km68.46630.000.110.0950.13915,30)km51.96601.620.090.095-0.09930,45)km24.91292.880.090.095-0.11445km4.0149.110.080.095-0.1530,200)m

20、59.75494.110.1210.0950.237200,400)m41.11400.710.1030.0950.082400,600)m23.31285.300.0820.095-0.152600m25.18393.500.0640.095-0.349 注：滑坡體總面積為149.34 km2,流域總面積為1 573.61 km2 表2 研究區(qū)泥石流敏感性評(píng)價(jià)因子綜合賦值Table 2 prehensive scores of debris flow susceptibility影響因子因子分級(jí)圍坡度032324040484886坡向NW(292.5337.5),SW(202.5247.5)

21、,NE(22.567.5),SW(112.5157.5),W(247.5292.5)N(022.5,337.5360)E(67.5112.5)S(157.5202.5)地震峰值加速0.25g(0.20,0.25g0.15,0.20g距斷層距離45km(30,45km(15,30km0,15km地層巖性泥盆紀(jì)石灰?guī)r三疊紀(jì)砂巖和千枚巖,志留紀(jì)板巖石炭系灰?guī)r,震旦砂巖和白云巖侏羅紀(jì)砂巖,二疊紀(jì)閃長(zhǎng)巖距水系距離600m(400,600m(200,400m0,200m賦值1234 4.4 敏感性因子權(quán)重賦值 選用層次分析法來構(gòu)建判別矩陣，同樣考慮坡度、坡向、地層巖性、地震峰值加速度、距水系距離、距斷層距

22、離共6個(gè)評(píng)價(jià)指標(biāo)，通過查閱文獻(xiàn)，讓專家打分，并結(jié)合實(shí)地野外調(diào)查等，最終建立判別矩陣(表3)，并且求得CR=0.04 表3 泥石流敏感性評(píng)價(jià)因子層次分析法判別系數(shù)矩陣Table 3 parative matri* and weight of evaluation factors affecting debris flow susceptibility因子B1B2B3B4B5B6WB113121/21/20.16B21/311/21/21/31/30.07B312121/21/30.14B41/221/211/31/20.10B5232311/20.24B62332210.30 注：B1為距水系距

23、離; B2為坡向;B3為地層巖性; B4為距斷層距離; B5為地震峰值加速度; B6為坡度；W為權(quán)重 根據(jù)判別矩陣可以看到，對(duì)泥石流敏感性影響最大的因素為坡度,其次為地震峰值加速度，距水系的距離也對(duì)泥石流敏感性有較為明顯的影響，這結(jié)果也跟我們實(shí)地野外調(diào)查觀察大致一樣，大量崩滑體多發(fā)生在坡度較大的坡頂上，沿著山脊線呈串珠狀分布，同時(shí)也有沿河呈條帶狀分布的特征。 4.5 敏感性評(píng)價(jià)模型建立及評(píng)價(jià)結(jié)果 根據(jù)前文層次分析法中計(jì)算出來的研究區(qū)泥石流6個(gè)敏感性評(píng)價(jià)因子權(quán)重值及研究區(qū)145條泥石流溝的各因子分?jǐn)?shù)值建立敏感性評(píng)價(jià)模型，公式為 (2) 式中:S為泥石流溝的綜合敏感度值；i=1,2,6;*i(k)

24、為研究區(qū)145條泥石流溝因子賦值；Wi為研究區(qū)6個(gè)泥石流敏感性評(píng)價(jià)因子權(quán)重。 根據(jù)研究區(qū)145條泥石流溝判別系數(shù)矩陣表里的值，在ArcGIS軟件里對(duì)泥石流溝屬性表進(jìn)展字段賦值，然后在ArcGIS里的字段計(jì)算器模塊中利用公式(2)，計(jì)算得到泥石流溝的綜合敏感度值。參考前人研究結(jié)果6，考慮到本文單溝的各評(píng)價(jià)因子都在ArcGIS軟件中進(jìn)展了定量化的處理并不是按整數(shù)14賦值，并結(jié)合實(shí)地對(duì)研究區(qū)泥石流災(zāi)害野外調(diào)查統(tǒng)計(jì)等因素將研究區(qū)145條泥石流溝的敏感性劃分為3類，得到泥石流溝敏感性評(píng)價(jià)圖，見圖2。 圖2 研究區(qū)泥石流敏感性分級(jí)Fig.2 Classification of the study area

25、 according to levels of debris flow susceptibility (1)S2.3，為低度敏感性泥石流流域，主要分布在汶川縣的西南部，未來可能有低頻小規(guī)模的泥石流發(fā)生，隨著時(shí)間變化，山體、植被恢復(fù)，流域可能會(huì)慢慢向高含沙水流為主的山洪溝演變。 (2)2.3S (3)S2.9，為高度敏感性泥石流流域，主要在岷江與漁子溪流域交匯處集中分布，未來可能爆發(fā)高頻大規(guī)模的泥石流，造成較為嚴(yán)重的生命財(cái)產(chǎn)損失，長(zhǎng)時(shí)間需要綜合泥石流防治措施治理。 5 結(jié) 論 根據(jù)概率綜合判別方法并運(yùn)用層次分析法，建立了研究區(qū)145條泥石流溝敏感性評(píng)價(jià)模型，得到研究區(qū)泥石流敏感性分級(jí)圖。根據(jù)這一

26、模型最終得到了60條泥石流屬于高敏感性，易發(fā)生泥石流，43條泥石流溝為中敏感性，需要定期監(jiān)測(cè)警覺突然爆發(fā)，其余42條泥石流溝敏感性較低，較為平安。 考慮到研究區(qū)汶川縣是地震強(qiáng)震區(qū)，因此選擇了坡度、坡向、地層巖性、地震峰值加速度、距水系距離、距斷層距離6個(gè)評(píng)價(jià)因子作為評(píng)價(jià)指標(biāo)，使用概率數(shù)學(xué)方法對(duì)各評(píng)價(jià)因子進(jìn)展數(shù)據(jù)統(tǒng)計(jì)分析，從大量原始數(shù)據(jù)中找到了一定的規(guī)律性，提高敏感性評(píng)價(jià)的客觀性，得到每個(gè)評(píng)價(jià)因子的因子分級(jí)圍及奉獻(xiàn)分值,使泥石流敏感性評(píng)價(jià)結(jié)果更加精準(zhǔn)可靠。 根據(jù)6個(gè)評(píng)價(jià)因子的判別區(qū)間權(quán)重賦值標(biāo)準(zhǔn)，結(jié)合研究區(qū)145條泥石流溝在上述評(píng)價(jià)因子的分布情況，得到研究區(qū)每條泥石流溝的敏感性量化分值，單溝的各

27、評(píng)價(jià)因子都在ArcGIS軟件中進(jìn)展了定量化的處理，使評(píng)價(jià)結(jié)果更加客觀可靠。 參考文獻(xiàn)： 1 黃潤(rùn)秋.汶川地震地質(zhì)災(zāi)害后效應(yīng)分析J 工程地質(zhì)學(xué)報(bào)，2011，19(2):145-151 2 鵬,韋方強(qiáng),何思明,等512汶川地震誘發(fā)的山地災(zāi)害及減災(zāi)措施J山地學(xué)報(bào)，2008，26(3)：280-282. 3 唐 川.汶川地震區(qū)暴雨滑坡泥石流活動(dòng)趨勢(shì)預(yù)測(cè)J山地學(xué)報(bào)，2010，28(3)：341-349. 4 AYALEW L,YAMAGISHI H.The Application of GIS-based Logistic Regression for Landslide Susceptibility

28、Mapping in the Kakuda-Yahiko Mountains，Central JapanJElsvier，2005，65(1):15-31. 5 唐 川.怒江流域泥石流敏感性空間分析J.地理研究,2005,24(2): 178-185. 6 CHANG T C.Risk Degree of Debris Flow Applying Neural NetworksJ.Natural Hazards,2007,42(1):209-224. 7 洋,唐 川,為樂,等.省都江堰市龍池地區(qū)群發(fā)性泥石流物源敏感性分析J.災(zāi)害學(xué),2013,28(2):107-113. 8 向靈芝, 泳,洪凱

29、,等.基于流域演化的泥石流敏感性分析J.長(zhǎng)江流域資源與環(huán)境,2015,24(11):1984-1992. 9 胡寶榮.基于遙感與GIS技術(shù)的汶川縣地震前后生態(tài)環(huán)境質(zhì)量評(píng)價(jià)D.:理工大學(xué),2009. 10 程 思,易加強(qiáng).省汶川縣地質(zhì)災(zāi)害的成因及防治對(duì)策J.地質(zhì)災(zāi)害與環(huán)境保護(hù),2007,18(4):1-4. 11 唐邦興，柳素清 .省阿壩藏族羌族自治州泥石流及其防治研究M : 科技大學(xué)，1993: 21-25. 12 周 偉,唐 川，周春花.汶川震區(qū)暴雨泥石流激發(fā)雨量特征J.水科學(xué)進(jìn)展,2012,23(5): 650-655. 13 DAI F C,*U C,YAO *,et al.Spatia

30、l Distribution of Landslides Triggered by the 2008 Ms 8.0 Wenchuan Earthquake,ChinaJ.Journal of Asian Earth Sciences,2012,40(4):883-895. 14 KAMP U， GROWLEY B J， KHATTAK G A，et al GIS-based Landslide Susceptibility Mapping for the 2005 Kashmir Earthquake RegionJGeomorphology，2008,101:631-642. 15 唐 川,

31、梁京濤.汶川震區(qū)北川924暴雨泥石流特征研究J.工程地質(zhì)學(xué)報(bào),2008,(6): 752-758. 16 WANG W N，NAKAMURA H，TSUCHIYA S，et alDistributions of Landslides Triggered by the Chi-Chi Earthquake in Central Taiwan on September 21，1999JJournal of the Japan Landslide Society，2002，38( 4): 318-326. 17 樹德,任秀生,岳升陽(yáng),等.地震與泥石流活動(dòng)J.水土保持研究,2001,8(2):26-2

32、8. 18 CHANG Ming,TANG Chuan,ZHANG Dan-dan,et al.Debris Flow Susceptibility Assessment Using a Probabilistic Approach: A Case Study in the Longchi Area,Sichuan Province,ChinaJJournal of Mountain Science,2014,11(4):1001-1014. 19 LI Wei-le,HUANG Run-qiu,TANG Chuan，et alCo-seismic Landslide Inventory an

33、d Susceptibility Mapping in the 2008 Wenchuan Earthquake Disaster Area,ChinaJ Journal of Mountain Science,2013,10(3):339-354. (編輯：衛(wèi)兵) Susceptibility Assessment of Debris Flow Using a Probabilistic andGIS Approach: A Case Study on the Wenchuan County *IA Chen-hao,ZHU Jing,CHANG Ming,YANG Yu (State Key Laboratory of Geo-hazard Prevention and Geo-environment Protection,Chengdu University of Technology,Chengdu 610059,China) Abstract:In this paper we mapped co-seismic landslides and