基于可拓灰色模型的臺(tái)風(fēng)暴雨型泥石流易發(fā)性評(píng)價(jià)

胡富杭　熊朝正　石豫川　吉鋒

摘要：2019年8月10日，受臺(tái)風(fēng)“利奇馬”登陸影響，安徽寧國(guó)抽水蓄能電站匯水區(qū)域內(nèi)6條臺(tái)風(fēng)暴雨型泥石流暴發(fā)，對(duì)工程建筑設(shè)施及居民生命財(cái)產(chǎn)造成了巨大損害。臺(tái)風(fēng)暴雨型泥石流的易發(fā)性研究目前尚處于起步階段，在對(duì)研究區(qū)內(nèi)地質(zhì)環(huán)境調(diào)查的基礎(chǔ)上，總結(jié)了臺(tái)風(fēng)暴雨型泥石流特征，將可拓理論與灰色關(guān)聯(lián)度法相結(jié)合，采用定量的方式對(duì)區(qū)內(nèi)泥石流進(jìn)行易發(fā)性評(píng)價(jià)，并將評(píng)價(jià)結(jié)果與規(guī)范推薦的評(píng)分法和三要素法進(jìn)行對(duì)比分析。研究結(jié)果表明：臺(tái)風(fēng)暴雨型泥石流具有超強(qiáng)短歷時(shí)降雨啟動(dòng)、低頻性、群發(fā)性的特征;在可拓灰色關(guān)聯(lián)法的易發(fā)性評(píng)價(jià)結(jié)果中僅虹龍河為中等易發(fā)，其余溝谷均為輕度易發(fā);對(duì)比發(fā)現(xiàn)，在小流域面積泥石流溝的評(píng)價(jià)結(jié)果中，3種方法一致性較高，當(dāng)流域面積增大時(shí)，規(guī)范推薦的方法出現(xiàn)偏高和偏低，而可拓灰色關(guān)聯(lián)法的易發(fā)性評(píng)價(jià)結(jié)果更加穩(wěn)定可靠，與實(shí)際調(diào)查情況更吻合，進(jìn)而為臺(tái)風(fēng)暴雨型泥石流的易發(fā)性研究提供了一種新的思路。

關(guān)鍵詞：臺(tái)風(fēng)暴雨型泥石流; 泥石流易發(fā)性評(píng)價(jià); 可拓理論; 灰色關(guān)聯(lián)度法; 安徽寧國(guó)抽水蓄能電站

中圖法分類(lèi)號(hào)： P642.23

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

DOI：10.16232/j.cnki.1001-4179.2021.10.005

0引 言

泥石流是一種暴發(fā)突然，運(yùn)動(dòng)迅速，影響范圍廣，破壞性強(qiáng)，具有多次致災(zāi)能力的山地地質(zhì)災(zāi)害。其常發(fā)育于地形變化大，物源儲(chǔ)量豐富的狹窄溝谷，在持續(xù)性降雨或突發(fā)型暴雨的作用下啟動(dòng)，沿途不斷侵蝕并裹挾溝道物質(zhì)，對(duì)溝道及溝口所在區(qū)域內(nèi)的工程建筑及居民生命財(cái)產(chǎn)安全造成巨大威脅。泥石流的易發(fā)性研究是泥石流危險(xiǎn)性評(píng)價(jià)的前提，也是區(qū)域泥石流預(yù)警預(yù)報(bào)和防治的基礎(chǔ)。20世紀(jì)70年代后，泥石流易發(fā)性研究逐步發(fā)展完善，從定性評(píng)價(jià)開(kāi)始向半定量評(píng)價(jià)、定量評(píng)價(jià)發(fā)展[1]。目前采用的方法有層次分析法、多因子綜合評(píng)判模型[2]、信息量權(quán)數(shù)法、機(jī)器學(xué)習(xí)算法、人工神經(jīng)網(wǎng)絡(luò)模型[3]、貝葉斯模型、灰色關(guān)聯(lián)分析法、可拓學(xué)理論等。田運(yùn)濤[1]等在對(duì)白龍江泥灣溝泥石流地質(zhì)調(diào)查及動(dòng)力學(xué)特征計(jì)算的基礎(chǔ)上選取了18個(gè)評(píng)價(jià)因子，應(yīng)用層次分析法評(píng)價(jià)了該泥石流為高易發(fā)性。王金亮等[4]應(yīng)用雅魯藏布江冰川泥石流實(shí)例介紹了主成分分析理論的概念和方法，并評(píng)價(jià)了該方法的實(shí)用性。Xu等[5]結(jié)合GIS技術(shù)和信息量方法，選取了7個(gè)環(huán)境因子作為信息量模型對(duì)四川省泥石流易發(fā)性的帶狀分布特征進(jìn)行了深入研究。王念秦等[6]以模糊數(shù)學(xué)為基礎(chǔ)，結(jié)合泥石流易發(fā)性特點(diǎn)，建立了多級(jí)模糊綜合評(píng)價(jià)模型，并應(yīng)用于西安市大水溝泥石流易發(fā)性評(píng)價(jià)中。Gómez等[7]以委內(nèi)瑞拉安地斯山脈雅蓬諾薩河流域?yàn)槔?，介紹了人工神經(jīng)網(wǎng)絡(luò)技術(shù)在地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)中的應(yīng)用，結(jié)果表明其評(píng)價(jià)精度高達(dá)90%。此外，許多學(xué)者在研究過(guò)程中還將多種方法結(jié)合，取得了不少研究成果[8-9]。葛濤濤[10]結(jié)合遙感數(shù)據(jù)，GIS和機(jī)器學(xué)習(xí)算法，對(duì)西藏日喀則地區(qū)泥石流易發(fā)性進(jìn)行了深入研究，并總結(jié)出其中起主要影響作用的環(huán)境因子。蔡健玲[10]采用貝葉斯模型對(duì)滇西北區(qū)域泥石流易發(fā)性空間分布進(jìn)行了評(píng)價(jià)，同時(shí)結(jié)合全因子及區(qū)域特色logistic回歸模型進(jìn)行對(duì)比，表明了貝葉斯模型在大尺度范圍泥石流評(píng)價(jià)中的重要作用。崔傳峰、金星等[12-13]在可拓學(xué)理論基礎(chǔ)上利用灰色關(guān)聯(lián)分析法獲取評(píng)價(jià)因子權(quán)重，建立了一種新的泥石流災(zāi)害評(píng)價(jià)模型，并應(yīng)用于多條泥石流溝災(zāi)害評(píng)價(jià)，結(jié)果表明該方法具有很高的適用性。

中國(guó)地勢(shì)總體上西高東低，據(jù)統(tǒng)計(jì)，傳統(tǒng)的泥石流災(zāi)害多分布在中國(guó)地形陡峭、植被貧瘠的西北、西南地區(qū)，經(jīng)過(guò)多年的科學(xué)研究和工程經(jīng)驗(yàn)積累，已形成了相當(dāng)成熟的易發(fā)性評(píng)價(jià)體系。近年來(lái)，隨著東部地區(qū)臺(tái)風(fēng)登陸愈加頻繁，由臺(tái)風(fēng)暴雨引發(fā)的臺(tái)風(fēng)暴雨型泥石流越來(lái)越成為東部地質(zhì)災(zāi)害防治的重點(diǎn)。王一鳴[14]等結(jié)合臺(tái)風(fēng)暴雨特點(diǎn)，以浙江省溫州市為主要研究區(qū)域，深入研究了臺(tái)風(fēng)暴雨型泥石流啟動(dòng)機(jī)制及單溝危險(xiǎn)度，并提出了此類(lèi)特殊泥石流的風(fēng)險(xiǎn)評(píng)價(jià)方法。Chen [15]、甘建軍[16]等通過(guò)野外調(diào)查和室內(nèi)實(shí)驗(yàn)，對(duì)臺(tái)風(fēng)暴雨型泥石流的形成條件和致災(zāi)模式進(jìn)行了深入研究，總結(jié)了該類(lèi)泥石流的災(zāi)害連鎖效應(yīng)、低頻性、群發(fā)性等特點(diǎn)。

但是，在東部地區(qū)，地形較為平緩，因其靠近大海，空氣常年濕潤(rùn)，利于植被的生長(zhǎng)發(fā)育，泥石流溝谷普遍具有較小的縱比降和較大的植被覆蓋率，按傳統(tǒng)泥石流的易發(fā)性評(píng)價(jià)方法及工程經(jīng)驗(yàn)常認(rèn)為其難以爆發(fā)泥石流，考察因素不夠全面，忽略了在臺(tái)風(fēng)暴雨這種極端天氣條件下爆發(fā)的可能性，難以做到有效的預(yù)防。因此，針對(duì)東部地區(qū)，為了準(zhǔn)確科學(xué)地進(jìn)行臺(tái)風(fēng)暴雨型泥石流易發(fā)性評(píng)價(jià)，需要建立一種考慮多方面因素，且要將各因素進(jìn)行量化的泥石流易發(fā)性評(píng)價(jià)方法。本文在總結(jié)臺(tái)風(fēng)暴雨型泥石流特點(diǎn)的基礎(chǔ)上，結(jié)合工程實(shí)例，建立了可拓灰色關(guān)聯(lián)模型進(jìn)行臺(tái)風(fēng)暴雨型泥石流易發(fā)性評(píng)價(jià)。

1研究區(qū)概況

研究區(qū)為影響擬建安徽寧國(guó)抽水蓄能電站的整個(gè)匯水區(qū)域（見(jiàn)圖1（a）），位于安徽省宣城市寧國(guó)市境內(nèi)。2019年9號(hào)超強(qiáng)臺(tái)風(fēng)“利奇馬”于8月10日在浙江溫嶺市登陸，受其影響，研究區(qū)當(dāng)日降雨量高達(dá)285～306 mm，共導(dǎo)致區(qū)內(nèi)多條溝谷暴發(fā)泥石流，不僅對(duì)擬建電站工程場(chǎng)地造成嚴(yán)重破壞，還對(duì)當(dāng)?shù)鼐用竦纳?cái)產(chǎn)造成了巨大損害。根據(jù)現(xiàn)場(chǎng)調(diào)查及遙感解譯，暴發(fā)的泥石流溝主要為虹龍河、青龍溝、羅漢塘溝、大石里溝、記陽(yáng)坑溝、大河壩溝，分布如圖1（a）所示。

研究區(qū)處于天目山山脈的中段西北坡，上水庫(kù)區(qū)正常蓄水位高程633 m，位于青龍山凹地，青龍溝是庫(kù)內(nèi)主溝，庫(kù)區(qū)平坦寬闊，庫(kù)周分水嶺封閉較好，山脊走向呈NE～NNE向展布，地形較為陡峭。下水庫(kù)區(qū)位于虹龍河圣龍寺村河段，正常蓄水位208 m，庫(kù)區(qū)河道蜿蜒狹窄，斷面呈“V”字型，植被茂密，庫(kù)周山坡坡度一般30°～45°，局部見(jiàn)陡崖、陡坎。下水庫(kù)主要受虹龍河、大石里溝、羅漢塘溝和記陽(yáng)坑溝影響。

區(qū)內(nèi)地層巖性主要為雷公塢組（Z1l）含礫凝灰?guī)r，藍(lán)田組（Z2l）薄層硅質(zhì)泥巖、白云質(zhì)灰?guī)r，皮園村組（Z2p）硅質(zhì)泥巖、炭質(zhì)硅質(zhì)巖。第四系堆積體在溝道兩側(cè)廣泛分布，結(jié)構(gòu)較為松散，是泥石流暴發(fā)的物質(zhì)基礎(chǔ)。研究區(qū)屬亞熱帶季風(fēng)濕潤(rùn)氣候區(qū)，年平均降水量為1 455.9 mm。區(qū)內(nèi)水流主要以降雨的形式補(bǔ)給，以地表水的形式向虹龍河排泄。

在這次9號(hào)臺(tái)風(fēng)“利奇馬”影響下，研究區(qū)所處霞西鎮(zhèn)日降雨量高達(dá)285 mm，如圖2所示。據(jù)目擊者描述，16：00左右（約臺(tái)風(fēng)登陸14 h）虹龍河水位暴漲，在半小時(shí)內(nèi)泥石流暴發(fā)，沖毀多處居民住房，掩埋工程施工設(shè)備。虹龍河和青龍溝是此次暴發(fā)最為嚴(yán)重的泥石流溝，虹龍河因其溝道長(zhǎng)，整體縱比降小，泥流堆積物呈現(xiàn)出多段堆積特征，主要集中在汪家園及塔里附近停留堆積（見(jiàn)圖1（b）和圖1（c））。青龍溝泥石流堆積體主要集中在較為寬闊的上水庫(kù)庫(kù)區(qū)，整體粒徑較小，僅少量塊石粒徑大于30 cm，據(jù)計(jì)算堆積物方量約2.49萬(wàn)m3（見(jiàn)圖1（d））。

對(duì)研究區(qū)內(nèi)6條臺(tái)風(fēng)暴雨型泥石流進(jìn)行分析后，總結(jié)出以下特征：

（1） 在超強(qiáng)短歷時(shí)降雨下啟動(dòng)。根據(jù)降雨數(shù)據(jù)（見(jiàn)圖2），研究區(qū)的小時(shí)降雨量在臺(tái)風(fēng)登陸后14.25 h達(dá)到極值46 mm，是前一記錄點(diǎn)的4.6倍，超強(qiáng)短歷時(shí)降雨特征十分突出。根據(jù)目擊者描述的泥石流啟動(dòng)暴發(fā)時(shí)段在16：00左右，與超強(qiáng)短歷時(shí)降雨出現(xiàn)的時(shí)段高度吻合。

（2） 暴發(fā)頻率低。一般受臺(tái)風(fēng)影響嚴(yán)重的東部沿海地區(qū)地形變化相對(duì)較小，而且常年濕潤(rùn)的氣候使得山區(qū)植被覆蓋率極高，在一般條件下難以暴發(fā)大規(guī)模泥石流。據(jù)訪問(wèn)，研究區(qū)內(nèi)溝谷雖然常年流水，但百年以來(lái)從未暴發(fā)具有破壞性的泥石流。

（3） 集群式暴發(fā)。據(jù)記載，2004年14號(hào)臺(tái)風(fēng)“云娜”造成浙江省樂(lè)清市北部山區(qū)41處泥石流集中暴發(fā)，“莫蘭蒂”造成泰順縣泗溪鎮(zhèn)22處泥石流集中暴發(fā)[14]。在本次臺(tái)風(fēng)“利奇馬”作用下，面積僅16.86 km2的研究區(qū)內(nèi)共暴發(fā)6條主要泥石流溝，同樣具有顯著的集群式暴發(fā)特點(diǎn)。

2可拓灰色關(guān)聯(lián)模型的構(gòu)建

1983年蔡文教授創(chuàng)立了可拓學(xué)，將事件、事件特征、特征量值通過(guò)建立物元的方式進(jìn)行有機(jī)結(jié)合，把事物的質(zhì)和量進(jìn)行綜合考慮，用關(guān)聯(lián)度的大小將特征性質(zhì)與研究對(duì)象的從屬關(guān)系進(jìn)行定量描述，針對(duì)問(wèn)題本身來(lái)解決矛盾，是一種新的人工智能工具。根據(jù)前人研究成果，可拓理論能很好地應(yīng)用于泥石流易發(fā)性評(píng)價(jià)。

3結(jié)果分析

3.1評(píng)價(jià)因子選擇及等級(jí)劃分

泥石流的暴發(fā)是多種影響因子共同耦合作用的結(jié)果，在不同地質(zhì)條件下影響因子的作用機(jī)理差異巨大。然而在各種客觀因素的限制下，難以將所有影響因子都納入泥石流易發(fā)性及危險(xiǎn)性評(píng)價(jià)中加以考慮。因此，需要根據(jù)研究區(qū)泥石流溝實(shí)際情況選擇具有代表性和明確物理意義的參評(píng)因子，這些因子不僅要易于獲得，能做數(shù)值化處理，而且要相互獨(dú)立。

在泥石流的發(fā)育中，山體坡度越大，邊坡物源越不穩(wěn)定，更易為泥石流提供固體物質(zhì)，同時(shí)也更利于匯流的快速產(chǎn)生。主溝長(zhǎng)度與補(bǔ)給段長(zhǎng)度占比表明了泥石流溝的物源儲(chǔ)量情況，即主溝越長(zhǎng)，補(bǔ)給段占比越大，則泥石流物源儲(chǔ)量越豐富。主溝長(zhǎng)度與主溝高差決定了溝床縱比降，縱比降的增加會(huì)促使固體物質(zhì)的動(dòng)勢(shì)能轉(zhuǎn)換，利于泥石流的啟動(dòng)發(fā)生。降雨量及溝谷流域面積控制著影響泥石流發(fā)生的水源條件，降雨量和流域面積的增加均會(huì)加大溝床徑流強(qiáng)度，使泥石流在強(qiáng)大的水動(dòng)力條件下更易啟動(dòng)，此外邊坡松散物源的穩(wěn)定性還會(huì)隨著降雨量的增加顯著降低，利于產(chǎn)生滑坡為泥石流提供固體物質(zhì)。溝谷兩岸生長(zhǎng)的植被能消減降雨沖擊和坡表徑流對(duì)表層土體的破壞作用，同時(shí)通過(guò)根系的生長(zhǎng)不斷加固邊坡表層土體，因此較高的植被覆蓋率能有效減小泥石流的易發(fā)性。

對(duì)研究區(qū)內(nèi)在臺(tái)風(fēng)期間暴發(fā)的6條泥石流溝的實(shí)測(cè)數(shù)據(jù)及遙感解譯數(shù)據(jù)進(jìn)行統(tǒng)計(jì)，如表1所列。其中，在這次“利奇馬”臺(tái)風(fēng)影響下，根據(jù)研究區(qū)降雨統(tǒng)計(jì)數(shù)據(jù)，24 h降雨量為285 mm，為評(píng)價(jià)區(qū)內(nèi)泥石流溝在臺(tái)風(fēng)暴雨影響下的易發(fā)性，對(duì)24 h最大降雨量c6均取值為285 mm。參照前人根據(jù)多年工程經(jīng)驗(yàn)建立的易發(fā)性分級(jí)標(biāo)準(zhǔn)[12，17，22]，同時(shí)結(jié)合研究區(qū)實(shí)際調(diào)查情況及研究特性，將7個(gè)評(píng)價(jià)因子及易發(fā)性按照《泥石流災(zāi)害防治工程勘查規(guī)范》中的定名分為不易發(fā)生（N1），輕度易發(fā)（N2），中度易發(fā)（N3），極度易發(fā)（N4）4個(gè)等級(jí)，如表2所列。

3.2灰色關(guān)聯(lián)度法確定評(píng)價(jià)因子的權(quán)重

在泥石流災(zāi)害的易發(fā)性評(píng)價(jià)中，選取的各評(píng)價(jià)因子對(duì)泥石流易發(fā)性的貢獻(xiàn)程度是不一樣的[13]，即各因子所占權(quán)重不等，單以關(guān)聯(lián)函數(shù)來(lái)衡量易發(fā)性不夠準(zhǔn)確，因此對(duì)參評(píng)因子權(quán)重的確定是關(guān)乎評(píng)價(jià)結(jié)果準(zhǔn)確與否的另一核心問(wèn)題。

目前對(duì)權(quán)重的計(jì)算多采用經(jīng)驗(yàn)法[12]、簡(jiǎn)單關(guān)聯(lián)函數(shù)法、專(zhuān)家評(píng)價(jià)法[19]、模糊綜合評(píng)價(jià)法[20-22]、灰色關(guān)聯(lián)度法等。為避免評(píng)價(jià)過(guò)程中主觀性的影響，通過(guò)比較，本文選用適用性強(qiáng)的灰色關(guān)聯(lián)度法計(jì)算各參評(píng)因子的權(quán)重。

3.3物元的構(gòu)造與計(jì)算

根據(jù)表1數(shù)據(jù)及表2易發(fā)性分級(jí)標(biāo)準(zhǔn)，構(gòu)建研究區(qū)臺(tái)風(fēng)暴雨型泥石流溝易發(fā)性評(píng)價(jià)經(jīng)典物元（Ri）與節(jié)域物元（R）如表3所示。

根據(jù)式（4）對(duì)指定泥石流溝S進(jìn)行關(guān)聯(lián)函數(shù)的計(jì)算，以虹龍河為例，計(jì)算結(jié)果如表4所示。根據(jù)該計(jì)算方法，分別對(duì)研究區(qū)其余5條泥石流溝進(jìn)行關(guān)聯(lián)函數(shù)計(jì)算。

3.4泥石流易發(fā)性計(jì)算與分析

前文已求得研究區(qū)6條泥石流溝的易發(fā)性關(guān)聯(lián)函數(shù)值及7個(gè)參評(píng)因子的權(quán)重值，通過(guò)建立的可拓灰色關(guān)聯(lián)模型式（5）計(jì)算泥石流溝關(guān)聯(lián)度ki（S），運(yùn)用式（6）對(duì)泥石流溝易發(fā)程度等級(jí)進(jìn)行劃分，結(jié)果見(jiàn)表5。

為了驗(yàn)證本文評(píng)價(jià)結(jié)果，筆者應(yīng)用規(guī)范中泥石流易發(fā)程度評(píng)價(jià)方法推薦的評(píng)分法和三要素評(píng)價(jià)法對(duì)研究區(qū)內(nèi)6條泥石流溝進(jìn)行了易發(fā)性判定，結(jié)果如表5所列。為了明顯對(duì)比3種方法的評(píng)價(jià)結(jié)果，將其以柱狀圖表示（見(jiàn)圖3）。整體上可明顯觀察到流域范圍較小的大河壩溝、記陽(yáng)坑溝、大石里溝、羅漢塘溝在3種評(píng)價(jià)方法中的判定等級(jí)均為輕度易發(fā)，與實(shí)際情況符合，但青龍溝和虹龍河的判定結(jié)果開(kāi)始出現(xiàn)差異。

對(duì)于青龍溝可拓法和三要素法都評(píng)價(jià)為輕度易發(fā)，評(píng)分法評(píng)價(jià)其為中等易發(fā)。根據(jù)青龍溝野外調(diào)查情況，物源主要發(fā)育于溝道中上部，青龍溝中部（擬建水庫(kù)區(qū)）地形平坦開(kāi)闊，當(dāng)溝床出現(xiàn)匯流時(shí)，大部分物質(zhì)將會(huì)停留堆積于此，使得青龍溝爆發(fā)泥石流的可能性大大降低，因此青龍溝整體上更加符合輕度易發(fā)等級(jí)，可拓法和三要素法的評(píng)價(jià)結(jié)果較為合適。

對(duì)于虹龍河，可拓法和評(píng)分法都評(píng)價(jià)為中等易發(fā)，而三要素法評(píng)價(jià)其為不易發(fā)?！袄骜R”登陸期間，虹龍河研究區(qū)內(nèi)是泥石流暴發(fā)最強(qiáng)烈、破壞最嚴(yán)重的泥石流溝。根據(jù)后期現(xiàn)場(chǎng)調(diào)查，虹龍河流域面積巨大達(dá)16.86 km2，且支溝發(fā)育，匯水條件良好，流域內(nèi)不穩(wěn)定物源總量達(dá)163.35萬(wàn)m3，其易發(fā)性更加符合中等易發(fā)等級(jí)。三要素法的評(píng)價(jià)結(jié)果和實(shí)際情況明顯不符，這是因?yàn)槠鋬H從小時(shí)雨強(qiáng)，不穩(wěn)定物源量及縱坡降3個(gè)方面進(jìn)行評(píng)價(jià)，評(píng)價(jià)因子較少，忽略了虹龍河巨大的流域面積，且由于虹龍河主溝長(zhǎng)度大，縱坡降極小，使得結(jié)果嚴(yán)重偏小。因此可拓法和評(píng)分法的評(píng)價(jià)結(jié)果較為合適。

經(jīng)過(guò)上述分析可知，可拓灰色關(guān)聯(lián)法對(duì)研究區(qū)內(nèi)泥石流的易發(fā)性評(píng)價(jià)都是合理的。而評(píng)分法和三要素評(píng)價(jià)法在應(yīng)用于流域面積較大的泥石流溝時(shí)會(huì)出現(xiàn)評(píng)價(jià)結(jié)果不準(zhǔn)確的情況，這是由于2種方法幾乎都是通過(guò)地質(zhì)調(diào)查人員的主觀判斷來(lái)對(duì)泥石流溝的各種影響因子進(jìn)行量化，在這過(guò)程中存在較大不確定性，而且地質(zhì)人員對(duì)泥石流溝的調(diào)查及認(rèn)識(shí)是有限的，當(dāng)泥石流流域面積足夠大時(shí)，地質(zhì)人員難以查清確切的地質(zhì)條件情況，在主觀評(píng)價(jià)時(shí)就會(huì)加大誤差，從而使得評(píng)價(jià)結(jié)果不夠準(zhǔn)確。可拓灰色關(guān)聯(lián)分析法在獲取評(píng)價(jià)因子時(shí)秉持代表性，明確物理意義，容易獲取，相互獨(dú)立和能數(shù)值化處理的原則，很大程度上避免了人的主觀影響，因此在評(píng)價(jià)結(jié)果上與實(shí)際調(diào)查情況更加符合。

4結(jié) 論

（1） 在9號(hào)臺(tái)風(fēng)“利奇馬”影響下，研究區(qū)日降雨量高達(dá)285 mm，通過(guò)對(duì)區(qū)內(nèi)爆發(fā)的泥石流溝進(jìn)行分析，總結(jié)出臺(tái)風(fēng)暴雨型泥石流具有以下特征：超強(qiáng)短歷時(shí)降雨啟動(dòng)、暴發(fā)頻率低、集群式暴發(fā)的特點(diǎn)。

（2） 計(jì)算表明選取的7個(gè)泥石流易發(fā)性評(píng)價(jià)因子cj在全部評(píng)價(jià)因子中所占權(quán)重λj分別為0.158 0，0.149 1，0.156 2，0.160 3，0.070 3，0.162 4，0.143 7。

（3） 運(yùn)用可拓灰色關(guān)聯(lián)法對(duì)研究區(qū)內(nèi)6條臺(tái)風(fēng)暴雨型泥石流進(jìn)行易發(fā)性評(píng)價(jià)，計(jì)算結(jié)果表明虹龍河關(guān)聯(lián)度為0.064 8，屬中等易發(fā);青龍溝、羅漢塘溝、大石里溝、記陽(yáng)坑溝、大河壩溝關(guān)聯(lián)度分別為0.057 3，-0.166 9，-0.084 7，-0.079 7，-0.048 3，均屬輕度易發(fā)。

（4） 將可拓灰色關(guān)聯(lián)法與規(guī)范推薦的評(píng)分法和三要素評(píng)價(jià)法進(jìn)行對(duì)比時(shí)發(fā)現(xiàn)，在評(píng)價(jià)流域面積較小的泥石流溝時(shí)3種方法評(píng)價(jià)結(jié)果一致，而當(dāng)流域面積過(guò)大時(shí)，評(píng)分法偏高1個(gè)等級(jí)，三要素法偏低2個(gè)等級(jí)，分析表明本文采用的可拓灰色關(guān)聯(lián)法更加客觀科學(xué)。

（5） 影響泥石流暴發(fā)的因素很多，可拓灰色關(guān)聯(lián)法在泥石流易發(fā)性評(píng)價(jià)中可避免定性指標(biāo)的選取，減少人為主觀影響，可拓理論中關(guān)聯(lián)函數(shù)計(jì)算與灰色關(guān)聯(lián)分析法中權(quán)重計(jì)算相結(jié)合，將定性問(wèn)題定量化，使評(píng)價(jià)結(jié)果更加真實(shí)可靠，該方法具有很好的適用性，為臺(tái)風(fēng)暴雨型泥石流的易發(fā)性研究提供了一種新思路。

參考文獻(xiàn)：

[1]田運(yùn)濤，葉振南，王高峰，等.白龍江泥灣溝泥石流特征及易發(fā)性評(píng)價(jià)研究[J].人民長(zhǎng)江，2019，50（增2）：59-64.

[2]劉希林.泥石流危險(xiǎn)度判定的研究[J].災(zāi)害學(xué)，1988（3）：10-15.

[3]ELKADIRI R，SULTAN M，YOUSSEF A，et al.A remote sensing-based approach for debris-flow susceptibility assessment using artificial neural networks and logistic regression modeling[J].IEEE Journal of Selected Topics in Applide Earth Observations and Remote Sensing，2014，7（12）：4818-4835.

[4]王金亮，翁秀紅.基于主成分分析法的泥石流危險(xiǎn)度綜合評(píng)價(jià)[J].地質(zhì)災(zāi)害與環(huán)境保護(hù)，2013，24（1）：88-92.

[5]XU W B，YU W J，JING S C，et al.Debris flow susceptibility assessment by GIS and information value model in a large-scale region，Sichuan Province（China）[J].Natural Hazards，2013，65（3）：1379-1392.

[6]王念秦，姚勇.基于模糊數(shù)學(xué)和權(quán)的最小平方法的泥石流易發(fā)性評(píng)價(jià)方法[J].災(zāi)害學(xué)，2008（2）：5-9.

[7]GMEZ H，KAVZOGLU T.Assessment of shallow landslide susceptibility using artificial neural networks in Jabonosa River Basin，Venezuela[J].Engineering Geology，2005，78：11-27.

[8]LINA H，JI Q Z，YI C Z，et al.Applying a Series and Parallel Model and a Bayesian Networks Model to Produce Disaster Chain Susceptibility Maps in the Changbai Mountain area，China[J].Water，2019，11（10）：2144.

[9]FANG C A，QIAN D A，Xu S A，et al.Debris flow susceptibility evaluation of Karakoram Highway based on extenics and entropy weight[J].Journal of Natural Disasters，2016，25（6）：18-26.

[10]葛濤濤.基于機(jī)器學(xué)習(xí)算法的日喀則地區(qū)泥石流易發(fā)性研究[D].南京：南京信息工程大學(xué)，2019.

[11]蔡健玲.GIS技術(shù)與貝葉斯層次模型支持下的滇西北泥石流易發(fā)性評(píng)價(jià)[D].昆明：昆明理工大學(xué)，2019.

[12]崔傳峰，王俊豪，崔志超，等.基于灰色可拓模型的洮河下游泥石流易發(fā)性評(píng)價(jià)[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2019，30（6）：40-48.

[13]金星，劉高，魏蒙恩，等.基于可拓灰色關(guān)聯(lián)分析法的泥石流危險(xiǎn)性評(píng)價(jià)[J].水資源與水工程學(xué)報(bào)，2014，25（4）：79-83.

[14]王一鳴.臺(tái)風(fēng)暴雨型泥石流災(zāi)害風(fēng)險(xiǎn)研究[D].武漢：中國(guó)地質(zhì)大學(xué)，2018.

[15]CHEN H Z，SU Z M，XU L R，et al.Characteristics and damage of debris flow events along Guangle expressway during Typhoon Utor，2013[J].Journal of Natural Disasters，2016，25（1）：65-73.

[16]甘建軍，儲(chǔ)小東.中低山臺(tái)風(fēng)暴雨型泥石流形成機(jī)制和動(dòng)力特征：以江西德安杜家溝泥石流為例[J].災(zāi)害學(xué)，2020，35（1）：150-155.

[17]劉希林，莫多聞.泥石流風(fēng)險(xiǎn)評(píng)價(jià)[M].成都/烏魯木齊：四川科技技術(shù)出版社/新疆科技衛(wèi)生出版社，2003：100-104.

[18]鄭立峰，唐見(jiàn)兵.基于灰色關(guān)聯(lián)分析的系統(tǒng)仿真因素權(quán)重確定[J].計(jì)算機(jī)仿真，2007（9）：76-78.

[19]林原，戰(zhàn)仁軍，吳虎勝.基于猶豫度和相似度的專(zhuān)家權(quán)重確定方法及其應(yīng)用[J/OL].控制與決策：1-8[2020-07-02].https：∥doi.org/10.13195/j.kzyjc.2019.1382.

[20]陳維，趙鑫，張海太，等.基于模糊綜合評(píng)判法的滇西北公路泥石流危險(xiǎn)性評(píng)價(jià)[J].公路，2020，65（2）：25-29.

[21]陳明，王運(yùn)生，梁瑞鋒，等.基于云模型的單溝泥石流危險(xiǎn)性模糊綜合評(píng)價(jià)[J].人民長(zhǎng)江，2018，49（10）：66-71.

[22]崔志超，王俊豪，崔傳峰，等.基于層次分析法和模糊數(shù)學(xué)相結(jié)合的甘肅東鄉(xiāng)八丹溝泥石流易發(fā)性評(píng)價(jià)[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2020，31（1）：44-50.

[23]呂擎峰，李軍鵬，趙本山，等.基于灰色關(guān)聯(lián)度法的可拓理論在霍爾古吐水電站壩址區(qū)泥石流危險(xiǎn)性評(píng)價(jià)中的應(yīng)用[J].工程地質(zhì)學(xué)報(bào)，2016，24（2）：211-219.

（編輯：劉 媛）

Abstract：On August 10，2019，affected by the landing of typhoon “Lekima”，six typhoon-triggered debris flows broke out in the catchment area of Ningguo Pumped Storage Power Station in Anhui Province，causing great damage to the engineering construction facilities and residents′ lives and properties.At present，the research on the susceptibility of typhoon-triggered debris flow is still in its infancy.Based on the investigation of the geological environment in the research area，this paper summarized the characteristics of typhoon-triggered debris flow.Combining the extension theory with the grey correlation method，we evaluated the susceptibility of debris flow in this area，and the evaluated results were compared with the standard recommended scoring method and the three factors method.The results showed that the typhoon-triggered debris flow had the characteristics of super short duration rainfall initiation，low frequency and massive occurrence;in the evaluation results of the extension grey correlation method，only the Honglong river was moderately prone，and the other gullies were mildly prone.The comparison shows that the three methods give highly consistent evaluation results for debris flow gullies with small valley area，however when the valley area increases，the results of standard recommended method is high or low，the extension grey correlation method give more stable and reliable results，which is more consistent with the actual investigation situation，thus providing a new idea for the research on susceptibility of typhoon-triggered debris flows.

Key words：typhoon-triggered debris flow;susceptibility evaluation on debris flow;extension theory;grey correlation degree method;Ningguo Pumped Storage Power Station in Anhui Province