張祝鴻，王保云，孔艷

（云南師范大學(xué)信息學(xué)院，昆明 650500）

0 引言

泥石流通常產(chǎn)生于溝谷或坡地上，呈現(xiàn)形式為飽含大量石塊、泥沙和巨礫的固液流體，是各種自然因素和人類活動的綜合產(chǎn)物，常具有突發(fā)性強、破壞力大的特點，對人類的生命財產(chǎn)安全造成極大的威脅[1]。中國擁有70%左右的山區(qū)面積，復(fù)雜的地形地貌以及顯著的季風(fēng)氣候，這為泥石流的形成發(fā)育提供了充足的物源條件、地質(zhì)條件和水源條件，因此我國多泥石流頻發(fā)地區(qū)。

截止到20世紀末期，由于泥石流災(zāi)害，我國每年死亡近千人，直接經(jīng)濟損失10億多人民幣[2]。并且，隨著人類社會經(jīng)濟活動的增強，泥石流發(fā)生頻率呈上升的發(fā)展趨勢。據(jù)統(tǒng)計資料顯示，我國有8500多條災(zāi)害性泥石流，其分布面積約占國土總面積的18.6%，其中有771個縣（市）有泥石流活動，140多座城鎮(zhèn)建設(shè)在泥石流堆積區(qū)上[3]。例如，2010年8月18日云南怒江傈僳族自治州貢山獨龍族怒族自治縣普拉底鄉(xiāng)突發(fā)泥石流災(zāi)害，造成37人死亡，53人失蹤，屬于特大型災(zāi)害事件[4]。2010年8月7日甘肅舟曲特大泥石流，沖毀房屋5500多間，造成1744人死亡和失蹤[5]。2017年6月，江西省中北部因超強降雨，導(dǎo)致泥石流、滑坡等災(zāi)害接連發(fā)生，共計328起，造成5人死亡，直接經(jīng)濟損失多達1700萬元。此外，強烈地震后會導(dǎo)致斜坡穩(wěn)定性降低，使山區(qū)地質(zhì)環(huán)境更加脆弱，一旦遭遇暴雨就會加劇泥石流災(zāi)害的發(fā)生[6]，如汶川地震后，暴雨誘發(fā)的北川縣發(fā)生72處泥石流災(zāi)害事件。

綜上所述，泥石流災(zāi)害對我國生產(chǎn)生活極具破壞力，如何對泥石流災(zāi)害的危險性進行準確的評價，對泥石流防災(zāi)減災(zāi)工作具有重要意義。目前，在泥石流危險性研究領(lǐng)域已經(jīng)提出了諸多泥石流危險性評價方法，并取得了較好的效果。本文在總結(jié)國內(nèi)外泥石流危險性評價研究進展的基礎(chǔ)上，對基于指標統(tǒng)計的泥石流危險性評價流程進行了概述，著重介紹了7種基于指標統(tǒng)計的泥石流危險性評價方法并總結(jié)其優(yōu)缺點。最后對目前存在的問題提出一些看法，為今后的研究奠定基礎(chǔ)。

1 國內(nèi)外泥石流危險性評價研究概況

泥石流的危險性評價是指在特定的地域和時間內(nèi)，泥石流發(fā)生的強度和可能性，并用危險度來定量地表示，危險度通常以區(qū)間范圍為[0,1]概率值表示[7-8]。自十九世紀以來，國內(nèi)外研究學(xué)者紛紛展開了相關(guān)的研究，并取得了卓而有效的成果。國外泥石流性評價較早研究較早，研究成果相對較多[9]。19世紀后半期俄國Б.И.斯塔科特夫邦基最初提到了泥石相關(guān)成因和危險度的問題[10]。1976年聯(lián)合國指派國際工程地質(zhì)聯(lián)合會負責(zé)滑坡泥石流危險度評價的基礎(chǔ)研究[11]，其研究成果為各個國家開展滑坡泥石流研究提供了一種規(guī)范性方法。此后的20世紀70～90年代，各個國家的研究者們陸續(xù)提出各種開創(chuàng)性的泥石流危險性評價方法[12-16]，這些成果為泥石流的后續(xù)研究奠定了基礎(chǔ)。90年代以后，各種適用于非線性計算的數(shù)學(xué)方法被引入到了泥石流的危險性評價并發(fā)揮了其作用，隨著計算機科學(xué)與技術(shù)的發(fā)展，紅外、遠紅外、無人機、地理信息系統(tǒng)GIS、遙感技術(shù)RS和全球定位系統(tǒng)GPS等也被運用到了泥石流的危險性評價中，泥石流空間數(shù)據(jù)的精確性和分析效率得到有效的提高[17-18]。泥石流危險性的研究逐步進入了一個階段，相關(guān)的研究成果很豐碩[19-23]，所取得的研究成果對人們的生產(chǎn)生活意義重大。

相比國外學(xué)者的研究，國內(nèi)起步較晚。20世紀90年代以前，國內(nèi)對泥石流仍處于一個認知的階段，研究的相關(guān)工作主要依托地質(zhì)水文等實地勘測，得到的結(jié)果也主要是對泥石流災(zāi)害的定性評論。1982王禮先年首次在對泥石流的定義進行了闡述[24]。1986年譚炳炎首次對泥石流的判斷及其嚴重程度的判定條件進行了研究，1991年發(fā)表的綜合評判結(jié)果在鐵路部門的得到了推廣應(yīng)用[25-26]。20世紀80年代后期，越來越多的學(xué)者投身到了泥石流危險性研究領(lǐng)域，也取得了一定成果[27-28]。20世紀90年代以后，越來越多的學(xué)者將數(shù)學(xué)模型引入到了國內(nèi)泥石流的研究中，這使得泥石流危險度的劃分更加合理客觀，隨著研究的深入，我國的研究成果不斷涌現(xiàn)，泥石流的規(guī)范性研究也在全國范圍內(nèi)全面展開[11]，如全國滑坡、泥石流災(zāi)害分布與災(zāi)害區(qū)劃圖[29]、中國云南小江流域泥石流危險度評估[30-31]和黃河上游崩塌滑坡泥石流危險度評估[32]等。隨著新技術(shù)的出現(xiàn)發(fā)展，泥石流危險性評價與其他學(xué)科、新技術(shù)、新方法的交叉研究逐漸引領(lǐng)了新的研究潮流。例如，2016年尹超等[33]利用遺傳算法和云模型改進AHP算法進行了公路沿線泥石流災(zāi)害危險性空間分析，并基于ArcGIS，以公路沿線泥石流災(zāi)害危險度為主要依據(jù)，完成了中國鐵沿線泥石流災(zāi)害危險性區(qū)劃方案。2017年劉希林等[34]比較了ICG模型和CF&LR組合模型應(yīng)用于我國大尺度區(qū)域滑坡泥石流災(zāi)害危險性評價結(jié)果的合理性和準確性，其研究結(jié)論為今后泥石流危險性評價研究提供了一定的參考。

2 基于指標統(tǒng)計的泥石流危險性評價的研究基礎(chǔ)和方法

在泥石流危險性評價方法中，基于指標統(tǒng)計的評價方法已成為較為成熟的研究方法。這類方法在應(yīng)用時主要分為單溝泥石流評價和區(qū)域性泥石流評價兩類。使用這類方法對泥石流進行危險性評價時需要先做一些基礎(chǔ)準備工作，即評價因子的選取、分級和評價單元的劃分。然后，在此基礎(chǔ)上構(gòu)建評價模型對泥石流的危險性做定量的評估；最終得到泥石流的危險性評價。圖1為基于指標的泥石流危險性評價方法的流程圖。

圖1 基于指標統(tǒng)計的泥石流危險性評價方法流程圖

2. 1 研究基礎(chǔ)

（1）評價因子的選取及分級

據(jù)相關(guān)資料統(tǒng)計，對泥石流的發(fā)生、發(fā)展、運動和堆積過程有影響的因子累計達70多種[35]。國內(nèi)外對泥石流危險評價因子的選取問題，一般從地形條件、工程地質(zhì)、水文條件三個方面來考慮。在選取泥石流危險性評價指標時，通常需要考慮評價指標的科學(xué)性、正確性、全面性、代表性、簡便性和實用性[9]。在實際研究中，一般會根據(jù)前人的經(jīng)驗或歷史文獻資料來選取評價因子，但要特別注意的是：一次泥石流的發(fā)生與當(dāng)?shù)乇旧憝h(huán)境背景有關(guān)系，在選取評價因子時要考慮當(dāng)?shù)貙嶋H狀況。例如：如果研究對象為遠離人群的荒僻地域時，認為影響因素不必要成為評價因子；當(dāng)研究對象的流域面積較小，各評價單元內(nèi)的降水差別不大時，降水這一因素可不作為評價因子。

在泥石流危險性評價中，為保證泥石流危險度評價的客觀性與科學(xué)性，需要根據(jù)評價因子對泥石流危險性影響的嚴重程度進行范圍的劃分，就是給每個評價因子分級。劃分的準則一般基于地質(zhì)學(xué)家對現(xiàn)場的整體認識或基于巖土工程，在這里不對劃分規(guī)則作具體闡述，可參考文獻[36-39]等。

（2）評價單元的劃分

式中：分別為兩級最小差及兩級最大差；ρ∈[0，1]為分辨系數(shù)，一般ρ取0.5，ρ越大，分辨率越大；ρ越小，分辨率越小。

泥石流危險性評價可分為單溝泥石流危險性評價和區(qū)域泥石流危險性評價。對于區(qū)域泥石流的評價，需要預(yù)先進行評價單元的劃分。Guzzetti將評價單元劃分為柵格、流域、地貌、地形以及均一單元5種[40]。在區(qū)域泥石流評價中，根據(jù)研究地域的實際特點科學(xué)客觀選擇不同的評價單元。

2. 2 評價方法

（1）層次分析法

在泥石流危險性評價中，層次分析法是一種較為主觀的賦權(quán)方法，主要用來確定因子的權(quán)重。岳溪柳[41]利用層次分析法結(jié)合GIS技術(shù)對貴州畢節(jié)地區(qū)進行了泥石流危險性評價，但存在氣象數(shù)據(jù)空間化，評價指標空間匹配的不確定性、評價因子選取存在一定的主觀性的問題。楊小鳳等[42]利用層次分析法對四川七盤溝進行了泥石流危險性評價，但層次分析法的實驗結(jié)果與因子分析法、灰色關(guān)聯(lián)法的評價結(jié)果相差較大。

（2）熵值法

熵是信息無序度的度量，信息熵越大，信息的效用值越小[43]。熵值法就是利用信息熵來確定各因子對泥石流的貢獻程度，某評價因子對應(yīng)的熵值越小，評價因子對應(yīng)的權(quán)值越大。夏晨皓利用熵值法對龍溪河流域泥石流進行了危險性評價，得到的結(jié)果與劉希林的公式進行匹配驗證表明了該評價方法對龍池地區(qū)的適用性[44]。在楊小鳳[42]對四川七盤溝地區(qū)的研究中，熵值法與因子分析法、層次分析法相比體現(xiàn)出更高的精度。

（3）灰色關(guān)聯(lián)法

（4）因子分析法

因子分析法用少數(shù)幾項不可觀測的公共的評價因子來描述評價指標的關(guān)系，目的是用盡可能最少個數(shù)的公共因子的線性函數(shù)與特殊因子之和來描述原來觀測的每一個評價因子，從而簡化分析過程。李雄峰等[48]利用因子分析法得到影響泥石流危險性的主因子，并進行了主次影響影響因子的重要性等級排序，得到了泥石流主要影響因子對泥石流危險性的貢獻度。牛岑岑用因子分析法確定了評價指標的客觀權(quán)重[9]。

（5）人工神經(jīng)網(wǎng)絡(luò)法

在泥石流危險性評價中，涉及的因素眾多，這些因素彼此之間相互作用，具有一定的復(fù)雜性和不確定性。而前向型人工神經(jīng)網(wǎng)絡(luò)在不知道數(shù)據(jù)的分布形式和變量之間關(guān)系的情況下可實現(xiàn)高度的非線性映射[49]。趙源等[50]基于前人的研究基礎(chǔ)得出整個BP神經(jīng)網(wǎng)絡(luò)的建立過程，對隱含層的層數(shù)的確定和加快收斂速度兩個方面均給出了比較明確的方法。曹祿來等[51]結(jié)合人工神經(jīng)網(wǎng)絡(luò)與T-S模糊系統(tǒng)理論建立的模糊神經(jīng)網(wǎng)絡(luò)模型，作者將實驗結(jié)果與劉希林的灰色聚類法和可拓物元方法進行縱向?qū)Ρ?，與線性內(nèi)插產(chǎn)生訓(xùn)練數(shù)據(jù)方法進行橫向?qū)Ρ?，結(jié)果證明了該模糊神經(jīng)網(wǎng)絡(luò)應(yīng)用于泥石流危險性評價的有效性和可行性。

（6）基于可拓理論的泥石流危險性評價方法

可拓理論運用在泥石流危險性評價中主要是用可拓集合的關(guān)聯(lián)函數(shù)值和關(guān)聯(lián)度的大小確定指標參數(shù)間的從屬關(guān)系，從而把屬于或不屬于某個危險性的定性問題轉(zhuǎn)化為定量描述。劉金鵬等[52]用可拓層次分析法對區(qū)域泥石流進行了研究。該方法在解決比較方案兩者之間的重要性賦值問題時，避免了只考慮兩種極端狀況的限制問題，并且使判斷矩陣更加合理。劉金鵬等[53]運用熵值法對評價因子進行權(quán)重賦值，然后在此基礎(chǔ)上運用改進的可拓方法——最大貼近度法建立了泥石流危險度評價模型，該方法在傳統(tǒng)可拓方法的基礎(chǔ)上增加了反映評價對象偏離評價等級程度的作用。

（7）組合評價的方法

組合評價的方法就是利用兩種或兩種以上主觀和客觀分析法綜合對泥石流進行危險性評價。如何將主觀分析法和客觀分析法得到的權(quán)值進行組合得到組合權(quán)值是關(guān)鍵。牛岑岑引入距離函數(shù)來反映兩種權(quán)重值的差異程度與對應(yīng)分配系數(shù)的差異程度一致這一關(guān)系，并采用線性加權(quán)方式，求解相應(yīng)的分配系數(shù)[9]。楊忠梅等[54]、孟凡奇等[55]運用最優(yōu)化理論與方法建立兩種權(quán)值之間的函數(shù)關(guān)系式，得到組合權(quán)重值。杜國梁等[56]在對岷縣板達溝進行泥石流危險性評價時，組合權(quán)重的值是通過AHP法得到的權(quán)重去修正熵值法得到的權(quán)重得到的。路遙等[57]結(jié)合灰色系統(tǒng)理論和和二元比較法，提出了一種關(guān)聯(lián)度-二元比較分析法，同時基于博弈論集結(jié)模型理論對關(guān)聯(lián)度-二元比較法確定的主觀權(quán)重與熵值法確定的客觀權(quán)重進行組合賦權(quán)。

2. 3 泥石流評價方法的優(yōu)缺點

下面對上述幾種基于指標統(tǒng)計的泥石流危險性評價方法的優(yōu)缺點進行闡述，見表1。旨在為改進這些方法提供一定的信息參考。

表1 評價方法優(yōu)缺點比較

3 結(jié)語

泥石流危險性評價是防災(zāi)減災(zāi)系統(tǒng)的重要組成部分，它的重要意義在于對泥石流的災(zāi)前評估，進而有利于減輕泥石流災(zāi)害的經(jīng)濟損失，制定防護、救援和避難措施，救災(zāi)甚至防災(zāi)都具有重要意義。泥石流的危險性評價方法應(yīng)盡可能提高其評價精度方能體現(xiàn)其存在意義。為了致力于提高泥石流危險性評價結(jié)果的精度，筆者綜合前人的研究提出以下問題并對研究的發(fā)展趨勢進行展望。

（1）文中詳述的一類基于評價因子統(tǒng)計分析的方法，發(fā)展較為成熟，使用廣泛，操作簡單，但評價因子的選取對其評價結(jié)果有直接影響。國內(nèi)外雖然都形成了一些選取評價因子的準則，但如何針對不同情況下的泥石流合理選擇評價因子仍是一個問題。

（2）泥石流的發(fā)生本身就涉及到多種因素，且這些因素與因素之間，與泥石流的主體之間的關(guān)系具有復(fù)雜性和不確定性。而基于評價因子統(tǒng)計的方法，本就是為了迎合這個復(fù)雜不確定的系統(tǒng)才引入的，它致力于描述系統(tǒng)中的這些關(guān)系，卻忽視了泥石流的動態(tài)活動機理。從某種程度上來說，對問題考慮地不夠全面導(dǎo)致得出的危險度值也不夠精確。

因此，為了提高評價的精度，就要使方法盡可能可靠。一是選擇評價因子時要盡可能減少主觀因素的影響，評價因子要向更能代表泥石流的破壞能力和未來的發(fā)展趨勢的方向靠近。二是選取的數(shù)學(xué)理論模型要盡可能較全面地描述這個復(fù)雜系統(tǒng)內(nèi)的關(guān)系。最后不僅僅要考慮基于評價因子統(tǒng)計的靜態(tài)評價，更應(yīng)該致力于對泥石流動態(tài)活動規(guī)律的把握，從泥石流的內(nèi)在活動機理入手對泥石流進行危險性評價。這樣泥石流的危險性評價將邁入另外一個層面，逐漸補足基于指標評價中一直以來缺失的那個部分。目前，從泥石流內(nèi)在動態(tài)機理入手的主要的評價方法是基于計算機數(shù)值模擬的泥石流危險性評價方法?；谠u價因子統(tǒng)計的評價方法和基于數(shù)值模擬的評價方法分別從兩個方面對泥石流進行了危險度評價，能否將二者聯(lián)合起來綜合二者優(yōu)點對泥石流進行危險性評價，這也是這個領(lǐng)域可以研究的一個方向。

此外，將GIS、RS、無人機等其他計算機技術(shù)不斷引入到泥石流的危險性評價中來將成為一種趨勢。一方面，這些技術(shù)將有利于大大縮短研究周期，從而更加有利于實時評價的實現(xiàn)并節(jié)省大量的人力物力。另一方面，這些技術(shù)將有利于探索到以前研究中對泥石流認識不夠清楚的地方，從而提高研究的精度。最后，泥石流的發(fā)生是一個復(fù)雜的過程，其中涉及到各種各樣的其他領(lǐng)域的知識，基于多學(xué)科的泥石流危險性評價方法對提高評價結(jié)果的精度有著不可忽視的作用。