王翔弘紳,余建華,胡桂勝,楊志全

(1.中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,四川 成都 610041;2.昆明理工大學(xué)公共安全與應(yīng)急管理學(xué)院,云南 昆明 650093;3.四川省公路規(guī)劃勘察設(shè)計研究院有限公司,四川 成都 610041;4.高原科學(xué)與可持續(xù)發(fā)展研究院,青海 西寧 810016)

多數(shù)地震能夠削弱山體結(jié)構(gòu)強度,并激發(fā)滑坡災(zāi)害,為泥石流形成提供物源,并在集中強降雨條件作用下,滑坡可能進一步發(fā)展為泥石流災(zāi)害,因此在強震后的地震災(zāi)區(qū)具有常形成滑坡-泥石流災(zāi)害鏈的特征[1-2]。以2008年5月12日汶川地震和2013年4月20日蘆山地震為例,地震作用下激活了山體間部分已存在的滑坡體,產(chǎn)生并激發(fā)了其他新的地質(zhì)災(zāi)害[3-5]。而在汶川地震后的幾年里,由于山體坍滑破壞形成了大量的松散堆積物源,在局部降雨的作用下形成了許多災(zāi)難性的泥石流災(zāi)害[6-7],如汶川縣映秀鎮(zhèn)2010年8月14日、2013年7月10日和2019年8月20日泥石流。盡管震后滑坡-泥石流災(zāi)害鏈已經(jīng)得到初步的研究和揭示,但是風(fēng)險評價工作依然缺乏,因此對滑坡-泥石流災(zāi)害鏈特征與風(fēng)險評價的研究具有重要的意義。

我國對多災(zāi)害災(zāi)害鏈的研究尚處于起步階段,近期的研究都集中在我國西南山區(qū),特別是以四川省汶川縣為中心,如尹云鶴等[2]建立了西南地區(qū)地震-滑坡-泥石流災(zāi)害鏈風(fēng)險防范措施框架,陸研等[8]和雨德聰?shù)萚9]運用信息流理論將降雨指標(biāo)與滑坡-泥石流災(zāi)害進行了因果分析。以往的研究更注重對災(zāi)害數(shù)據(jù)的統(tǒng)計分析,但是需要投入大量的時間與精力,且效率和準(zhǔn)確率都偏低。隨著GIS技術(shù)的發(fā)展,復(fù)雜山地災(zāi)害的研究進程不斷推進,已有研究主要利用災(zāi)害耦合預(yù)測對災(zāi)害鏈進行預(yù)測分析[10-11]。但是目前國外學(xué)者在這方面的研究較少,國內(nèi)學(xué)者針對不同災(zāi)害鏈提出了相應(yīng)的理論模型與災(zāi)害鏈防治措施框架,如:胡凱衡等[12]建立了滑坡-泥石流災(zāi)害鏈風(fēng)險評估的理論框架和概念體系;馬金珠[13]通過遙感解譯和科學(xué)實驗,在白龍江流域建立了多尺度滑坡-泥石流綜合風(fēng)險評估體系與技術(shù)示范點;喬建平等[14]以工程效益四要素(風(fēng)險區(qū)劃、風(fēng)險概率﹑風(fēng)險受損率、風(fēng)險防御)為基礎(chǔ)建立了小流域滑坡-泥石流災(zāi)害鏈風(fēng)險綜合評估模型,為更多災(zāi)害鏈的研究提供了參考。然而,對于相同災(zāi)害類型和不同災(zāi)害類型個體之間或不同災(zāi)害類型之間的相互作用過程和形成機理方面的研究仍然相對薄弱,針對地質(zhì)災(zāi)害鏈體系,其評價方法、研究方法和技術(shù)體系還不完善和系統(tǒng)化[15],在這一領(lǐng)域仍然有待深入研究。

針對上述滑坡-泥石流災(zāi)害鏈研究和風(fēng)險分析的不足,本文以四川省石棉縣回隆村馬頸子溝低頻泥石流(2013年7月4日18時30分左右,馬頸子溝發(fā)生泥石流,其沖出固體物質(zhì)總量達30萬m3)為研究對象,分析了該滑坡-泥石流災(zāi)害鏈特征,并基于多災(zāi)害理論對馬頸子溝滑坡-泥石流災(zāi)害鏈進行了綜合風(fēng)險評價,以為類似研究提供理論依據(jù)和經(jīng)驗。

1 研究區(qū)概況

馬頸子溝流域地處青藏高原橫斷山脈東部石棉縣(雅安市西南部,北緯28°51′～29°32′,東經(jīng)101°55′～102°34′)回隆村(石棉縣縣城東南向6.8 km),南椏河右岸,溝口地理坐標(biāo)為北緯29°10′52.55"、東經(jīng)102°23′26.14",該泥石流溝位置如圖1所示。馬頸子溝流域內(nèi)地形陡峭,屬深高中山地區(qū),總體上以“V”型谷為主,流域面積為8.4 km2,相對高差1 957 m,溝床平均縱坡降為310‰。該流域趨勢呈上陡下緩,上游溝段(清水區(qū))溝道相對較陡,地形起伏較大,溝道縱坡降為487‰;中游溝段(形成區(qū))溝道縱坡降為250‰左右;中下游溝段為流通區(qū),縱坡降約為128‰,下游至溝口為堆積區(qū),縱坡降約為107‰,最終匯入南椏河[16]。馬頸子溝流域內(nèi)發(fā)育支溝4條,分布于上游溝段清水區(qū),如圖2所示。

圖1 馬頸子溝流域位置圖Fig.1 Location map of Majingzi Gully watershed

圖2 馬頸子溝流域特征Fig.2 Characteristics of Majingzi Gully watershed

石棉縣多年氣象統(tǒng)計數(shù)據(jù)及其線性分布,如表1和圖3所示。石棉縣多年平均年降水量為770 mm,日最大降雨量為105.6 mm,當(dāng)?shù)亟邓恐饕杏?至8月份內(nèi),占全年降水量的70%以上[16-17],經(jīng)常以陣雨或暴雨的形式出現(xiàn),降雨多集中在山地區(qū)域,河谷地帶次之。

表1 石棉縣多年氣象統(tǒng)計數(shù)據(jù)表

圖3 石棉縣多年氣象數(shù)據(jù)統(tǒng)計分布圖Fig.3 Statistical distribution of annual meteorological data in Shimian County

3) 第四系崩積層(Qcol):主要分布于溝谷兩岸坡腳地帶,由不穩(wěn)定斜坡和崩塌等堆積而成,物質(zhì)松散。

結(jié)合現(xiàn)場調(diào)查得知,區(qū)域土石比約為54.76∶45.24,綜合評定區(qū)域巖性為軟硬相間。

2 馬頸子溝滑坡-泥石流災(zāi)害鏈特征

2.1 地震和滑坡對泥石流災(zāi)害的影響分析

2.1.1 地震對泥石流災(zāi)害的影響分析

文廣超等[19]對汶川地震前后四川省主要地質(zhì)災(zāi)害時空發(fā)育規(guī)律進行了統(tǒng)計分析,結(jié)果表明:“5·12”汶川地震之前,川西區(qū)域構(gòu)造整體穩(wěn)定,滑坡、泥石流災(zāi)害數(shù)量較少,災(zāi)害規(guī)模較小;震后以汶川為中心,周邊區(qū)域出現(xiàn)了大量滑坡體、泥石流災(zāi)害等。與震前相比,震后四川省(近汶川區(qū)域為中心)地質(zhì)災(zāi)害強度在不斷增加,活躍性在不斷增強。

地震會進一步激發(fā)不良地質(zhì)體的破壞,土石下移至溝道,改變溝道微地貌,填埋甚至阻斷河床,嚴(yán)重情況下會形成堰塞壩[20]。研究區(qū)主要地震活動包括2008年5月12日發(fā)生的汶川地震,同年7月底在該縣發(fā)生的一次4.9級地震,以及隨后在2013年4月20日發(fā)生的蘆山地震,地震影響烈度為5～7度,尤其是2013年蘆山地震,馬頸子溝內(nèi)出現(xiàn)了輕微變形破壞跡象,地震導(dǎo)致山體松動,利于降雨入滲,成為汛期出現(xiàn)坡體變形破壞和泥石流的重要因素[16-17]。汶川地震和蘆山地震后,馬頸子溝流域內(nèi)未發(fā)現(xiàn)堰塞湖,僅溝床和兩側(cè)地形發(fā)生輕微變化。

地震的影響還會在極短時間內(nèi)以一種極端的方式加速區(qū)域內(nèi)巖土體轉(zhuǎn)化為松散物質(zhì)的過程。地震對馬頸子溝泥石流物源特征的影響顯著,其不光表現(xiàn)在物源總量上的增加,還對物源的分布特征和物質(zhì)結(jié)構(gòu)特征有著明顯的改變[20],馬頸子溝泥石流物源點分布位置和物源量統(tǒng)計,如圖4和表2所示。這也是繼“4·20”蘆山地震后發(fā)生“7·4”泥石流沖出固體物質(zhì)總量達30×104m3的重要因素,其總量約占可啟動物源量的42.4%[16]。另外,研究區(qū)內(nèi)人類活動(房屋修建、墾植開荒)活躍,對坡面物源和溝道物源有顯著的促進作用。

表2 馬頸子溝泥石流物源量統(tǒng)計表

圖4 馬頸子溝泥石流物源點分布位置圖Fig.4 Distribution of source points of debris flow in Majingzi Gully

2.1.2 滑坡對泥石流災(zāi)害的影響分析

在地震作用下,馬頸子溝兩側(cè)出現(xiàn)了多處崩滑體(圖5),這些崩滑體上部巖土體隨滑動面逐漸滑移至溝床,由于滑體后緣基巖裸露,表面光滑,導(dǎo)致雨水滲流和匯流的條件發(fā)生改變,在局部暴雨情況下,水流不斷沖刷下方松散堆積土體,水體下滲,崩滑體下緣土體流動性增強,易被溝道水流攜帶,形成泥石流災(zāi)害,沿途摻雜溝道旁側(cè)崩滑物源,擴大泥石流的規(guī)模,增加泥石流的危險性[20]。

圖5 2013年后馬頸子溝流域內(nèi)增加的典型滑坡區(qū)域分布圖Fig.5 Distribution of typical landslide area in Majingzi Gully watershed after 2013

通過現(xiàn)場調(diào)查,在馬頸子溝道中下游周邊出現(xiàn)了許多失穩(wěn)滑坡體。自“4·20”蘆山地震以來,馬頸子溝流域內(nèi)斜坡上不良地質(zhì)體在局部降雨形成的地表徑流(水動力)與風(fēng)蝕的作用下逐年增加,流域內(nèi)先后發(fā)生了7次滑坡(2.8×104～5.5×104m3),溝道內(nèi)的松散固體物源也隨之增多,泥石流易發(fā)程度提高,危險性增加。根據(jù)現(xiàn)場調(diào)查和統(tǒng)計,馬頸子溝流域內(nèi)2013年以后增加的典型滑坡區(qū)域分布,如圖5所示。

2.2 滑坡-泥石流災(zāi)害鏈演化過程分析

自“7·4”泥石流以來,馬頸子溝流域內(nèi)多為洪水沖刷、局部小型泥石流和水石流。通過調(diào)查統(tǒng)計該泥石流溝域內(nèi)可能參與泥石流活動的物源量為70.7×104m3(表2),這些物源并非同時參與一次泥石流活動,且一次參與泥石流活動的松散固體物質(zhì)也并非都會全部沖出泥石流溝進入主河,參照近期有史以來規(guī)模最大的泥石流活動(2013年7月4日泥石流)情況,“7·4”泥石流主要淤積段為G05、G07和溝口段(圖4)。研究發(fā)現(xiàn)這種情況主要受到以下幾個方面因素的影響:

1) 馬頸子溝降雨分布的不均一特性也決定了物源啟動參與泥石流活動的不均一性(當(dāng)泥石流物源分布集中區(qū),尤其馬頸子溝道中游-中下游區(qū)域,出現(xiàn)集中降雨或暴雨洪水等條件時,崩滑體或坡積物才可能啟動參與泥石流活動),馬頸子溝流域內(nèi)崩滑體及其他物源等要分多次參與泥石流活動。

2) “7·4”泥石流災(zāi)害中溝床堆積物參與泥石流活動的方式主要為揭底沖刷,在泥石流運動過程中,隨著溝道寬度和縱坡降的變化,泥石流在G05、G07和溝口段發(fā)生水沙分離,未沖出泥石流溝而是部分固體物質(zhì)沿溝道堆積。調(diào)查點G05和G07現(xiàn)場泥石流淤積特征如圖6所示。

圖6 馬頸子溝道特征點G05和G07泥石流淤積特征Fig.6 Debris flow deposition characteristics of G05 and G07 in Majingzi Gully

3) 該泥石流發(fā)生當(dāng)日,在特大暴雨情況下,溝道兩側(cè)崩滑體侵蝕剝落土石進入主溝道,并攜卷溝床堆積物,同時形成的泥石流又對兩側(cè)殘坡積物進行沖刷和侵蝕,從而形成滑坡-泥石流災(zāi)害鏈。

按水源成因和物源成因,馬頸子溝泥石流屬典型的暴雨溝谷型泥石流,通過總結(jié)近期泥石流情況可知,該泥石流規(guī)模主要與溝域內(nèi)松散固體物源的累積和動態(tài)變化情況以及與引發(fā)泥石流的暴雨情況相關(guān),當(dāng)溝域內(nèi)松散固體物源累積較多,且遇到集中暴雨時,往往發(fā)生局部巖土體崩滑,進而演變?yōu)槟嗍鳛?zāi)害。馬頸子溝滑坡-泥石流災(zāi)害鏈演化過程,如圖7所示。

圖7 馬頸子溝滑坡-泥石流災(zāi)害鏈演化過程Fig.7 Evolution process of landslide-debris flow disaster chain in Majingzi Gully

目前,馬頸子溝泥石流的危害對象主要為回隆村(溝口)和石龍村(中下游)靠近溝道兩側(cè)的住戶及公路、G108國道等(圖8),據(jù)統(tǒng)計,受威脅戶數(shù)77戶,受威脅人數(shù)251人。

圖8 石龍村“7·4”泥石流淤埋房屋和道路Fig.8 Roads and houses buried by “7·4” debris flow in Shilong Village

3 馬頸子溝滑坡-泥石流災(zāi)害鏈綜合風(fēng)險評價

3.1 評價方法與標(biāo)準(zhǔn)

災(zāi)害風(fēng)險科學(xué)研究的對象是由災(zāi)害、地理環(huán)境和暴露單元組成的“災(zāi)害系統(tǒng)”,具有區(qū)域性、互聯(lián)性、耦合性和復(fù)雜性的特點。區(qū)域多災(zāi)害、災(zāi)害鏈和災(zāi)害復(fù)合可能會導(dǎo)致復(fù)雜的影響,擴大或減弱災(zāi)害強度,并改變受影響地區(qū)的范圍[21]。

自聯(lián)合國發(fā)起“國際減災(zāi)十年(IDNDR)”和致力于減少災(zāi)害風(fēng)險的全球聯(lián)合努力以來,“降低災(zāi)害風(fēng)險(DRR)”科學(xué)和技術(shù)有了顯著的發(fā)展,科學(xué)家們針對區(qū)域多災(zāi)害、災(zāi)害鏈和災(zāi)害復(fù)合的情況,提出了定性災(zāi)害鏈的多災(zāi)害風(fēng)險理論公式[21]。災(zāi)害鏈中災(zāi)害之間存在某些鏈?zhǔn)疥P(guān)系,災(zāi)害的危險度也隨著災(zāi)害鏈長度的增加而增加,針對地區(qū)滑坡-泥石流災(zāi)害鏈(自汶川地震后,成為川西主要地質(zhì)災(zāi)害鏈和研究對象)的風(fēng)險評價,提出該類型災(zāi)害鏈風(fēng)險度計算公式[22]如下:

Rdc=Req+Rls+Rdf

=HeqVeqEeq+HlsVlsEls+HdfVdfEdf

=HeqVeqEeq+HlsVls(Eeq-Req)+
HdfVdf(Els-Rls)

(1)

式中:Rdc為震后泥石流災(zāi)害鏈造成的綜合損失風(fēng)險;Req、Rls、Rdf分別為地震、滑坡、泥石流災(zāi)害鏈中的風(fēng)險度;Veq、Vls、Vdf分別為地震、滑坡、泥石流造成損失的脆弱性;Eeq、Els、Edf分別為地震、滑坡、泥石流區(qū)域的人口密度(人/km2);Heq為地震危險度;Hls為滑坡危險度,主要通過分析滑坡發(fā)生的概率和運動得到;Hdf為泥石流危險度,由構(gòu)建的評價方法計算得到。

危險/風(fēng)險度分級標(biāo)準(zhǔn)如表3所示[23]。

表3 地震-滑坡-泥石流災(zāi)害鏈中危險/風(fēng)險度分級標(biāo)準(zhǔn)

災(zāi)害鏈風(fēng)險的表現(xiàn)為主要災(zāi)害風(fēng)險和次要災(zāi)害風(fēng)險的累積。據(jù)現(xiàn)場調(diào)查,地震對于石棉縣的影響僅限于促進當(dāng)?shù)乇阑w結(jié)構(gòu)的穩(wěn)定性破壞(間接增加了流域物源量),并沒有造成房屋倒塌和人員、財產(chǎn)的直接損失,不排除將來回隆鎮(zhèn)發(fā)生大型地震的可能性。由于該區(qū)域地震的相關(guān)數(shù)據(jù)缺失和計算存在較大誤差,區(qū)域災(zāi)害鏈風(fēng)險分析對于Req(地震風(fēng)險)級別定為低,針對研究區(qū)特征,上述公式(1)中Req取級別低對應(yīng)數(shù)(0～0.35],綜合取值Req=0,即本次針對馬頸子溝滑坡-泥石流災(zāi)害鏈的風(fēng)險分析不加入地震風(fēng)險分析。

另外,公式(1)中計算滑坡和泥石流影響下人口密度取決于地震威脅范圍內(nèi)人口的密度,該模式主要針對于區(qū)域性災(zāi)害鏈的計算,對于單溝并不受用。通過對往年馬頸子溝流域已發(fā)生的滑坡與崩塌的位置進行分析,滑坡與崩塌并沒有對當(dāng)?shù)鼐幼〉木用褡》炕蚓用竦娜松戆踩拓敭a(chǎn)造成直接的影響和危害,只是作為泥石流形成的物源條件,因此計算中Eeq、Els、Edf的值存在誤差,在此省略該值。綜合分析下,最后嘗試采用如下公式對研究區(qū)滑坡-泥石流災(zāi)害鏈的風(fēng)險度進行計算:

Rdc=Rls+Rdf=HlsVls+HdfVdf

(2)

該公式(2)分別對滑坡的危險度和泥石流的危險度進行計算,最后耦合為滑坡-泥石流災(zāi)害鏈的風(fēng)險度。其中,采用目前工程中評價滑坡危險度的手段,通過構(gòu)建層次分析法(AHP)結(jié)構(gòu)模型對馬頸子溝滑坡危險度進行計算[24];針對泥石流危險度的計算則采用泥石流危險指標(biāo)進行估算。對于滑坡和泥石流危險度的計算和分析采用上述兩種方式得出的結(jié)果均符合表3的分級標(biāo)準(zhǔn),可以實行。

3.2 滑坡風(fēng)險度(Rls)計算

3.2.1 構(gòu)建AHP結(jié)構(gòu)模型

從影響滑坡穩(wěn)定性的內(nèi)因和外因出發(fā),并考慮各因子量化的可操作性,以泥石流溝兩側(cè)山體地形地貌特征(B1)、工程地質(zhì)特征(B2)、其他特征(B3)為二級指標(biāo),以坡形(C1)、坡度(C2)、坡高(C3)、巖土類型(C4)、結(jié)構(gòu)面發(fā)育程度(C5)、巖土體內(nèi)摩擦角(C6)、巖土體黏聚力(C7)、年平均降雨量(C8)、植被覆蓋率(C9)、巖土風(fēng)化程度(C10)、地震烈度(C11)11個評價因子作為三級指標(biāo),建立AHP結(jié)構(gòu)模型,見圖9[24]。

圖9 AHP結(jié)構(gòu)模型Fig.9 AHP structure analysis model

3.2.2 計算指標(biāo)權(quán)重

采用1～9標(biāo)度法[24-25]構(gòu)造判斷矩陣,并根據(jù)上述建立的AHP結(jié)構(gòu)模型計算出各判斷矩陣中各評價因子(指標(biāo))的權(quán)重值(表4至表7),最后由下而上將C級目標(biāo)與B級目標(biāo)權(quán)重值相乘,對11個評價因子權(quán)重由大到小進行排序,見表8。

表4 A-B判斷矩陣和評價因子權(quán)重計算值

表5 B1-C判斷矩陣和評價因子權(quán)重計算值

表6 B2-C判斷矩陣和評價因子權(quán)重計算值

表7 B3-C判斷矩陣和評價因子權(quán)重計算值

表8 評價因子權(quán)重總排序

3.2.3 建立滑坡危險度綜合評價數(shù)學(xué)模型

運用德爾斐法[26]以及結(jié)合實際情況,邀請專家對評價因子進行打分得到其分級指數(shù)Xi,并通過滑坡危險度綜合評價數(shù)學(xué)模型將Xi與權(quán)重Wi組合,計算得到滑坡綜合危險分級指數(shù)Hls。滑坡危險度綜合評價數(shù)學(xué)模型如下:

(3)

式中:Hls為滑坡危險度;Wi為各評價因子的權(quán)重;Xi為各評價因子的分級指數(shù),分別取值為0.2、0.4、0.6、0.8、1.0五個等級[21]。

將各評價因子的Wi(表8)和分級指數(shù)Xi取值分別代入公式(1),經(jīng)計算可得:

Hls=0.232×0.8+0.203×0.6+0.201×0.6+0.129×0.6+0.058×0.4+0.053×0.4+0.040×0.4+0.033×0.4+0.027×0.2+0.012×0.2+0.012×0.2=0.589

由表3可知,綜合評定研究區(qū)滑坡災(zāi)害(目前調(diào)查有6處)屬于中度危險,與現(xiàn)場調(diào)查結(jié)果一致。

根據(jù)區(qū)域滑坡特征,將滑坡災(zāi)害造成損失的脆弱性(Vls)分為4個易損性等級,即:低易損性(Vls取值0～0.3)、中易損性(Vls取值0.3<～0.4)、較高易損性(Vls取值0.4<～0.6)、高易損性(Vls取值0.6<～1)[27],Vls計算采用紀(jì)虹等[28]構(gòu)建的滑坡災(zāi)害易損性等級計算模型進行計算,其計算公式如下:

(4)

通過計算可得xi=3,按照滑坡災(zāi)害易損性等級劃分標(biāo)準(zhǔn)為等級三,即滑坡災(zāi)害點易損性為較高易損性區(qū)域,取Vls=0.6。據(jù)調(diào)查(圖5)顯示,研究區(qū)典型滑坡多集中在形成區(qū)和流通區(qū),對于溝口居民區(qū)的影響較小,因此研究區(qū)滑坡風(fēng)險度定性趨于低風(fēng)險度,目前滑坡處于穩(wěn)定狀態(tài)。

3.3 泥石流風(fēng)險度(Rdf)計算

對馬頸子溝泥石流易損度的確定,根據(jù)文獻[22-27,29]所統(tǒng)計的2006年、2009年和2015年岷江上游84個鄉(xiāng)鎮(zhèn)泥石流災(zāi)害造成損失的易損性(代表岷江上游段不同區(qū)域泥石流災(zāi)害造成損失的易損性數(shù)值),發(fā)現(xiàn)易損性數(shù)值大小排序為2009年>2006年>2015年。馬頸子溝大規(guī)模泥石流繼2013年“4·20”蘆山地震幾個月后發(fā)生,與映秀鎮(zhèn)牛圈溝2008年發(fā)生的大規(guī)模泥石流類似,均屬于典型的震后泥石流,且映秀鎮(zhèn)和回隆鎮(zhèn)地形地貌相仿,主要分布花崗巖地層。結(jié)合兩地多種條件相似的情況,且目前距2013年馬頸子溝大規(guī)模泥石流已有9 a,該區(qū)域易損性在相應(yīng)的防治措施建設(shè)后已有所改善,因此依據(jù)岷江上游鄉(xiāng)鎮(zhèn)泥石流災(zāi)害易損性調(diào)查統(tǒng)計數(shù)據(jù),選擇2015年映秀鎮(zhèn)易損性數(shù)值[29]作為馬頸子溝泥石流易損性取值,即取

Vdf=0.230 6

(5)

這里應(yīng)用泥石流危險區(qū)的指標(biāo)來估算泥石流的危險度[22],則有:

Hdf=y1w1+y2w2+…+ynwn

(6)

其中,根據(jù)現(xiàn)場調(diào)查和計算數(shù)據(jù),將泥石流峰值流量(y1)、泥石流流速(y2)、相對高差(y3)、地形起伏度(y4)、巖性(y5)、植被覆蓋率(y6)、可啟動物源量(y7)和人類活動(y8)作為泥石流災(zāi)害評價的因子,根據(jù)各評價因子不同的危險閾值將其分為從一級至五級5個評價等級,分別賦值為0.2、0.4、0.6、0.8和1.0[21],相應(yīng)的危險閾值取值如表9所示;wi為各評價因子權(quán)重。泥石流災(zāi)害各評價因子yi的權(quán)重wi(采用AHP[30-31]法確定)取值如下表10所示[16-17,32]。

表9 各泥石流災(zāi)害危險閾值等級劃分和取值

表10 各評價因子權(quán)重計算匯總

將表9和表10數(shù)據(jù)代入公式(6),計算得出泥石流的危險度Hdf=0.757,即馬頸子溝泥石流的危險度為高度危險區(qū),與現(xiàn)場調(diào)查結(jié)果相符。

3.4 災(zāi)害鏈風(fēng)險分析

匯總上述計算結(jié)果所得Hls=0.589、Rls=0.35,Hdf=0.757,Rdf=0.17,研究區(qū)滑坡災(zāi)害為中危險度、中風(fēng)險度,泥石流災(zāi)害為高危險度、低風(fēng)險度,研究區(qū)滑坡災(zāi)害風(fēng)險度大于泥石流災(zāi)害風(fēng)險度的原因主要受區(qū)域災(zāi)害脆弱性取值的影響,綜合分析區(qū)域滑坡災(zāi)害屬于易發(fā),馬頸子溝為低頻泥石流溝(20～100 a一遇)[33],泥石流災(zāi)害不易發(fā),Vdf取值參考岷江上游數(shù)據(jù)統(tǒng)計結(jié)果,綜合判定結(jié)果可取。

綜合上述,馬頸子溝滑坡和泥石流危險度和區(qū)域災(zāi)害脆弱性相關(guān)計算數(shù)據(jù)代入公式(2),求得Rdc=0.46,馬頸子溝滑坡-泥石流災(zāi)害鏈為中風(fēng)險度。

目前,當(dāng)出現(xiàn)極端氣候,激發(fā)馬頸子溝道旁側(cè)滑坡體或演變新的不穩(wěn)定滑坡體,該區(qū)域依然具有形成滑坡-泥石流災(zāi)害鏈的可能性,若再次發(fā)生大規(guī)模的泥石流災(zāi)害,有堵塞馬頸子溝流域中下游和溝口處的可能,在沒有防治整改措施的情況下,對G108國道的暢通和溝口處居民的人身和財產(chǎn)安全存在較大的威脅。馬頸子溝口堆積扇淤埋房屋和G108國道現(xiàn)狀,如圖10所示。

圖10 馬頸子溝口堆積扇淤埋房屋和G108國道現(xiàn)狀Fig.10 Present situation of silted houses with Majingzi Gully accumulation fan and national highway G108

4 結(jié) 論

本文基于多災(zāi)害風(fēng)險理論,采用災(zāi)害鏈風(fēng)險評價方法結(jié)合多種計算模型分析了馬頸子溝滑坡-泥石流災(zāi)害鏈危險性,并通過資料收集和計算結(jié)果分析了馬頸子溝滑坡-泥石流災(zāi)害鏈風(fēng)險度。得到主要結(jié)論如下:

1) 自2013年4月20日蘆山地震以來,馬頸子溝流域內(nèi)先后發(fā)生了7次滑坡(2.8×104～5.5×104m3),溝道內(nèi)的松散固體物源增多,泥石流易發(fā)程度提高,危險性增加。隨后2013年7月4日爆發(fā)的大規(guī)模泥石流災(zāi)害,一次固體物質(zhì)沖出量約為30萬m3,約占可啟動物源量的42.4%。

2) 構(gòu)建AHP模型計算得到馬頸子溝滑坡災(zāi)害的危險度Hls=0.589,評價等級為中危險度。

3) 應(yīng)用泥石流危險區(qū)的指標(biāo)來估算泥石流的危險度,通過泥石流災(zāi)害評價因子危險閾值和權(quán)重統(tǒng)計,計算得出馬頸子溝泥石流災(zāi)害的危險度Hdf=0.757,評價等級為高危險度。

4) 基于多災(zāi)害風(fēng)險理論的災(zāi)害鏈風(fēng)險評價計算模型得到Rdc=0.46,表明馬頸子溝滑坡-泥石流災(zāi)害鏈為中風(fēng)險度。

5) 基于多災(zāi)害風(fēng)險理論的風(fēng)險評價是將不同災(zāi)種的風(fēng)險度疊加耦合而得出區(qū)域災(zāi)害鏈的風(fēng)險度,并根據(jù)實際情況,針對特定的災(zāi)害鏈進行分析,找出鏈?zhǔn)綖?zāi)害中的主要災(zāi)害和次要災(zāi)害進行綜合分析。目前僅應(yīng)用于地震-滑坡、崩塌-泥石流相關(guān)的多災(zāi)害鏈中,其他路徑還有待創(chuàng)新。本文在滑坡-泥石流災(zāi)害鏈中進行了相關(guān)嘗試,其評價結(jié)果與現(xiàn)場調(diào)查結(jié)果相符,相關(guān)方法后期還需進一步改進。