遼寧暴雨誘發(fā)山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃研究

耿世波 聶哲 王瀛 方緣 肖光梁 李鴻強(qiáng)紀(jì)永明 朱憲龍 周曉宇

（1.中國(guó)氣象局沈陽(yáng)大氣環(huán)境研究所，遼寧 沈陽(yáng) 110166；2.遼寧省氣象災(zāi)害監(jiān)測(cè)預(yù)警中心，遼寧 沈陽(yáng) 110166；3.營(yíng)口市氣象局，遼寧 營(yíng)口 115000；4.中國(guó)氣象局氣象干部培訓(xùn)學(xué)院遼寧分院，遼寧 沈陽(yáng) 110166；5.大連市金州區(qū)氣象局，遼寧 大連 116600；6.沈陽(yáng)區(qū)域氣候中心，遼寧 沈陽(yáng) 110166）

引言

山洪災(zāi)害是一種洪水災(zāi)害［1］。山洪由于爆發(fā)迅猛、影響因素眾多以及精細(xì)山洪資料缺乏等原因，對(duì)山洪的預(yù)報(bào)預(yù)警難度很大［2-3］。受全球氣候變暖影響，近年來(lái)遼寧省強(qiáng)降水的極端天氣呈現(xiàn)多發(fā)頻發(fā)趨勢(shì)，如2017年8月4日鞍山市岫巖縣石灰窯鎮(zhèn)5 h降水達(dá)到481.7 mm，超過(guò)了遼寧24 h過(guò)程降水歷史極值；2020年8月31日夜間，遼寧省莊河市突降特大暴雨，1 h雨強(qiáng)為133 mm，達(dá)到歷史極值。局地的短時(shí)強(qiáng)降水極易造成山洪和地質(zhì)災(zāi)害，嚴(yán)重威脅社會(huì)經(jīng)濟(jì)發(fā)展和人民財(cái)產(chǎn)安全。也做了一些工作。李吉順等［7］根據(jù)歷史災(zāi)害資料，編制了暴雨洪澇災(zāi)害區(qū)劃指南；國(guó)內(nèi)的一些學(xué)者也先后編寫(xiě)了自然災(zāi)害風(fēng)險(xiǎn)評(píng)估技術(shù)方法和氣象風(fēng)險(xiǎn)預(yù)警業(yè)務(wù)技術(shù)指南［8-10］，同時(shí)也有學(xué)者開(kāi)展了洪水淹沒(méi)模型的模擬災(zāi)損評(píng)估［11-13］。

總的來(lái)說(shuō)，山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃是對(duì)山洪危險(xiǎn)性、承災(zāi)體暴露度和易損性等進(jìn)行綜合評(píng)價(jià)［4］。國(guó)外有關(guān)學(xué)者對(duì)暴雨洪澇災(zāi)害評(píng)估做了大量研究。Ahmad和Simonovic［5］、Kim和Marcouiller［6］分析暴雨洪澇災(zāi)害的形成是承災(zāi)體脆弱性、致災(zāi)因子和暴露度等多方面因素綜合作用的結(jié)果。

關(guān)于山洪災(zāi)害風(fēng)險(xiǎn)性區(qū)劃，中國(guó)學(xué)者在近些年

近年來(lái)，針對(duì)遼寧地區(qū)暴雨災(zāi)害風(fēng)險(xiǎn)評(píng)估，韓秀君等［14］利用歷史暴雨災(zāi)情資料，構(gòu)建了暴雨致災(zāi)指標(biāo)與災(zāi)害影響預(yù)評(píng)估的關(guān)系。劉元?jiǎng)伲?5］利用降水?dāng)?shù)據(jù)，建立了遼寧省山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃；陳殿強(qiáng)等［16］利用遼寧省的地質(zhì)災(zāi)害數(shù)據(jù)，詳細(xì)分析了遼寧省山洪地質(zhì)災(zāi)害特點(diǎn)及空間分布規(guī)律；曹永強(qiáng)等［17］利用GIS技術(shù)，繪制了遼寧山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃圖。以上研究多關(guān)于遼寧山洪風(fēng)險(xiǎn)區(qū)劃，僅針對(duì)降水或?yàn)?zāi)情空間變化進(jìn)行統(tǒng)計(jì)，對(duì)于遼寧綜合暴雨致災(zāi)風(fēng)險(xiǎn)性、承災(zāi)體暴露度和脆弱性的山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估與區(qū)劃研究較少。因此，本文利用多源災(zāi)害調(diào)查數(shù)據(jù)、高精度的降水?dāng)?shù)據(jù)等綜合分析，嘗試構(gòu)建遼寧暴雨誘發(fā)山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃模型，為遼寧山洪災(zāi)害精確防治提供參考。

1 資料與方法

1.1 資料來(lái)源

所使用的數(shù)據(jù)包括氣象數(shù)據(jù)、地理信息數(shù)據(jù)、山洪溝數(shù)據(jù)和風(fēng)險(xiǎn)普查數(shù)據(jù)。氣象數(shù)據(jù)為遼寧地區(qū)2005—2019年1639個(gè)區(qū)域自動(dòng)氣象站逐時(shí)降水?dāng)?shù)據(jù)（圖1），來(lái)源于遼寧省氣象信息中心。地理信息數(shù)據(jù)為30 m空間分辨率的DEM地形數(shù)據(jù)，來(lái)源于美國(guó)國(guó)家航空航天局（NASA）。山洪溝數(shù)據(jù)為來(lái)源于遼寧水利部門。風(fēng)險(xiǎn)普查數(shù)據(jù)為遼寧地區(qū)1951—2013年山洪災(zāi)情資料，來(lái)源于沈陽(yáng)區(qū)域氣候中心暴雨洪澇災(zāi)害精細(xì)化風(fēng)險(xiǎn)普查資料。

圖1 遼寧地形及選取的自動(dòng)觀測(cè)氣象站分布Fig.1 Topography of Liaoning province and distribution of selected automatic observation meteorological stations

1.2 暴雨洪澇災(zāi)害風(fēng)險(xiǎn)評(píng)估

基于風(fēng)險(xiǎn)評(píng)估理論［8-9］，以遼寧省山洪溝作為基本研究單元，先后從暴雨致災(zāi)因子危險(xiǎn)性、山洪溝本身孕災(zāi)環(huán)境敏感性、山洪溝附近承災(zāi)體暴露度和承災(zāi)體災(zāi)損敏感性來(lái)選取指標(biāo)，通過(guò)構(gòu)建遼寧省暴雨誘發(fā)山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃模型得到綜合風(fēng)險(xiǎn)指數(shù)，進(jìn)而表征山洪氣象風(fēng)險(xiǎn)等級(jí)。

（1）致災(zāi)因子。選取6 h暴雨強(qiáng)度和頻率的綜合指標(biāo)進(jìn)行暴雨致災(zāi)危險(xiǎn)性評(píng)估，選取6 h暴雨作為評(píng)估致災(zāi)危險(xiǎn)性的依據(jù)將在致災(zāi)危險(xiǎn)性計(jì)算中說(shuō)明。

（2）孕災(zāi)環(huán)境敏感性。主要選擇山洪溝溝口位置海拔高度、溝床比降、河網(wǎng)密度等指標(biāo)。

（3）承災(zāi)體暴露度。由于承災(zāi)體年變化較大且準(zhǔn)確評(píng)估其客觀價(jià)值很困難，為了簡(jiǎn)便評(píng)估，選取了易評(píng)估且年變化比較小的耕地比例、人口密度兩個(gè)物理暴露度指標(biāo)。

2 結(jié)果分析

2.1 評(píng)估指標(biāo)計(jì)算

2.1.1致災(zāi)危險(xiǎn)性

為分析遼寧省山洪與暴雨的關(guān)系，以歷史山洪災(zāi)害為樣本進(jìn)行統(tǒng)計(jì)。由于2010—2013年山洪災(zāi)情資料較為完整，因此在這4 a內(nèi)隨機(jī)選取撫順縣、清原縣、新賓縣、鳳城市、蓋州市5個(gè)縣市共106次發(fā)生山洪災(zāi)害樣本和20次無(wú)災(zāi)害樣本，分別用山洪災(zāi)害發(fā)生前3 d的降水量（3d-rain）、前1 d的降水量（1 d-rain）、當(dāng)日最大12 h的降水量（12 h-rain），當(dāng)日最大6 h的降水量（6 h-rain）、當(dāng)日最大3 h的降水量（3 h-rain）、當(dāng)日最大1 h的降水量（1 hrain）分別與樣本序列計(jì)算相關(guān)系數(shù)。

山洪災(zāi)害與不同時(shí)效的降水相關(guān)性見(jiàn)表1。由表1可以看出，發(fā)生山洪災(zāi)害的樣本與當(dāng)日最大6 h、12 h的降水量相關(guān)性較好，相關(guān)系數(shù)分別達(dá)到0.268和0.222，通過(guò)了置信度0.01的顯著性檢驗(yàn)。當(dāng)日最大6 h暴雨量與發(fā)生山洪相關(guān)系數(shù)最大，這與文明章等［18］統(tǒng)計(jì)福建暴雨誘發(fā)山洪的強(qiáng)降水持續(xù)時(shí)間一般為4—6 h大體一致，因此本文選取6 h作為山洪致災(zāi)臨界雨量的累計(jì)時(shí)間。

統(tǒng)計(jì)2005—2019年遼寧省自動(dòng)氣象站連續(xù)6 h降水量大于50 mm的資料，并將所有連續(xù)6 h降水量大于50 mm自動(dòng)氣象站排序，按照第98、第95、第90、第80、第60百分位數(shù)分別確定其對(duì)應(yīng)的雨量閾值［19］，這樣劃分連續(xù)6 h暴雨量的分級(jí)強(qiáng)度見(jiàn)表2，1—5級(jí) 分 別 對(duì) 應(yīng)6 h雨 強(qiáng)R≥147.4 mm、20.1 mm≤R＜147.4 mm、101.6 mm≤R＜120.1 mm、84.4 mm≤R＜101.6 mm、68.4 mm≤R＜84.4 mm，這樣得到6 h雨量臨界閾值要比劉鳴彥等［20］對(duì)清原“0816”洪水模擬2 h閾值明顯偏大，與臨界閾值的時(shí)效差別有關(guān)。采用加權(quán)綜合評(píng)價(jià)法，應(yīng)用乘積模型對(duì)暴雨致災(zāi)危險(xiǎn)性指數(shù)（H）計(jì)算

式（1）中，ωi為分級(jí)權(quán)重系數(shù)，結(jié)合高展等［19］的研究，1—5級(jí)分別取5/15、4/15、3/15、2/15、1/15為權(quán)重系數(shù)；RtDi為標(biāo)準(zhǔn)化分級(jí)6 h暴雨，Rni為分級(jí)6 h暴雨頻次，其中RtDi標(biāo)準(zhǔn)化計(jì)算公式如下

式（2）中，Rti、Rtmaxi、Rtmini分別為6 h平均暴雨量實(shí)際值、最大值和最小值，Rnmaxi、Rnmini分別為年均6 h分級(jí)暴雨頻次的最大值、最小值。

表2 暴雨強(qiáng)度（R）的劃分Table 2 Division of rainstorm intensity

利用GIS反距離插值方法，得到2005—2019年遼寧省6 h分級(jí)暴雨頻次與暴雨量空間分布（圖2）。由圖2可看出，1—5級(jí)暴雨強(qiáng)度高頻次均在遼寧省東南部山區(qū)（大連、丹東地區(qū)、岫巖縣、蓋州市），其次為遼寧省西南部山區(qū)（葫蘆島地區(qū)），再次零星分布在遼河流域和東部山區(qū)，遼寧省北部地區(qū)處于暴雨強(qiáng)度低頻區(qū)，這種空間分布與陳傳雷等［21］統(tǒng)計(jì)分析定義遼寧長(zhǎng)歷時(shí)降水空間分布大體一致，區(qū)別是通過(guò)此次方法分類，發(fā)現(xiàn)遼寧東南部的暴雨高頻區(qū)更加突出。從分級(jí)來(lái)看，強(qiáng)度最大（1級(jí)）的暴雨位于遼寧南部和東部的莊河市、瓦房店市、普蘭店市、岫巖縣、蓋州市、丹東市、東港市和葫蘆島市（圖2b），1級(jí)暴雨強(qiáng)度最大頻次為0.188—0.378次·a-1，2—4級(jí)暴雨強(qiáng)度發(fā)生頻次分布與總體分布基本一致，5級(jí)暴雨強(qiáng)度最大頻次出現(xiàn)在遼寧東部寬甸為1.334—2.746次·a-1。總體來(lái)看，遼寧最強(qiáng)暴雨基本發(fā)生在東南山區(qū)面臨黃海的迎風(fēng)坡上，隨海拔的升高而增加。

圖2 2005—2019年遼寧省6 h暴雨強(qiáng)度分級(jí)（a）、1級(jí)暴雨頻次（b）、2級(jí)暴雨頻次（c）、3級(jí)暴雨頻次（d）、4級(jí)暴雨頻次（e）、5級(jí)暴雨頻次（d）空間分布Fig.2 6-hour rainstorm intensity classification（a），1st grade rainstorm frequency（b），2nd grade rainstorm frequency（c），3rd grade rainstorm frequency（d），4th grade rainstorm frequency（e），5th grade rainstorm frequency frequency（f）in Liaoning province from 2005 to 2019

2.1.2孕災(zāi)環(huán)境敏感性

依據(jù)水沿著斜坡最陡方向流動(dòng)的原理［22］，將遼寧省流域面積≤200 km2劃分為519個(gè)山洪溝。山洪的爆發(fā)與山洪溝的地形地貌有很大關(guān)系，地形地貌特征影響著區(qū)域內(nèi)匯水與排水條件。山洪溝自然風(fēng)險(xiǎn)屬性特征主要由溝口的高程、溝床比降、河網(wǎng)密度3個(gè)指標(biāo)來(lái)表示匯水與排水的強(qiáng)弱。統(tǒng)計(jì)結(jié)果表明，遼寧省山洪溝溝口高程范圍在18—742 m，溝床比降為0—115，如圖3所示。參考文獻(xiàn)［23］給出了溝口的高程和溝床比降地形因子賦值（表3），這樣可得出遼寧省山洪溝地形地貌特征的風(fēng)險(xiǎn)區(qū)劃。

根據(jù)許海軍和陳守煜［24］定義，河網(wǎng)密度為單位面積內(nèi)河流數(shù)量，如果在一個(gè)區(qū)域內(nèi)的河流越多，匯水的概率增大，出現(xiàn)山洪災(zāi)害的風(fēng)險(xiǎn)也更大。依據(jù)張小峰等［26］河網(wǎng)密度與河道總長(zhǎng)度關(guān)系，將流域面積記為A，在面積A內(nèi)有河流或河段條數(shù)記為n，總長(zhǎng)度記為L(zhǎng)；則有

式（3）中，l為河流長(zhǎng)度；a為河流面積。

那么該區(qū)域的河網(wǎng)密度方程可表示為

圖3 遼寧省山洪溝溝口高程（a）、溝床比降（b）和河網(wǎng)密度（c）空間分布Fig.3 Spatial distribution ofmountain torrent gully mouth elevation（a），gully bed gradient（b），river network density（c）in Liaoning province

根據(jù)式（3）—式（4），計(jì)算出遼寧省山洪溝區(qū)域內(nèi)的河網(wǎng)密度為0—6，最大河網(wǎng)密度的山洪溝為北崴子河，河網(wǎng)密度為6。按照自然斷點(diǎn)法，將河網(wǎng)密度分別賦予相應(yīng)等級(jí)指數(shù)，得到遼寧省山洪溝河網(wǎng)密度指數(shù)的風(fēng)險(xiǎn)區(qū)劃，應(yīng)用加權(quán)求和模型進(jìn)行山洪溝孕災(zāi)環(huán)境敏感指數(shù)（S）估算

式（5）中，Te和Ri分別為地形因子和河網(wǎng)密度；ω1和ω2為其對(duì)應(yīng)的權(quán)重系數(shù)。歸一化計(jì)算公式如下［23］

式（6）中，D為經(jīng)標(biāo)準(zhǔn)化后某個(gè)評(píng)價(jià)指標(biāo)值；A、Amin、Amax分別為某評(píng)價(jià)指標(biāo)在標(biāo)準(zhǔn)化前實(shí)際值、最小值和最大值［26］。文獻(xiàn)［27］和本地專家打分確定各因子權(quán)重系數(shù)，得到地形因子、河網(wǎng)密度的權(quán)重ω1、ω2依次為0.6和0.4，因此將孕災(zāi)環(huán)境劃分為5個(gè)等級(jí)敏感指數(shù)，依次為低敏感區(qū)、次低敏感區(qū)、中等敏感區(qū)、次高敏感區(qū)和高敏感區(qū)。

2.1.3承災(zāi)體暴露度

承災(zāi)體易損性指數(shù)V按照如下方法計(jì)算

式（7）—式（8）中，Ve為承災(zāi)體暴露度風(fēng)險(xiǎn)指數(shù)；Vd為承災(zāi)體災(zāi)損敏感系數(shù)；P、C分別為標(biāo)準(zhǔn)化后人口密度和耕地面積比例系數(shù)；ω1、ω2分別為相應(yīng)指標(biāo)的系數(shù)權(quán)重，對(duì)各因子標(biāo)準(zhǔn)化方法依據(jù)式（6）。而權(quán)重確定依據(jù)專家打分，確定人口密度和耕地比例系數(shù)的權(quán)重ω1、ω2分別為0.7和0.3。越容易發(fā)生山洪災(zāi)害的地方，越容易造成人員傷亡、農(nóng)田沖毀和漬澇，根據(jù)統(tǒng)計(jì)的1951—2013年山洪災(zāi)害發(fā)生頻率，構(gòu)建災(zāi)損敏感系數(shù)Vd（表4）。最后根據(jù)式（7），建立遼寧省山洪災(zāi)害承災(zāi)體分析指標(biāo)，同樣采用自然斷點(diǎn)的方法依次將山洪溝劃分低風(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ū)5個(gè)等級(jí)。

2.1.4山洪災(zāi)害風(fēng)險(xiǎn)綜合區(qū)劃

根據(jù)文獻(xiàn)［26］暴雨洪澇災(zāi)害風(fēng)險(xiǎn)指數(shù)（D）模型

式（9）中，H、S、V依次為標(biāo)準(zhǔn)化后的致災(zāi)危險(xiǎn)性指數(shù)、孕災(zāi)環(huán)境敏感性指數(shù)和易損性指數(shù)；ω3、ω4、ω5依次為相應(yīng)指標(biāo)的系數(shù)權(quán)重。系數(shù)權(quán)重確定方法如下：根據(jù)本地專家打分，確定各項(xiàng)指標(biāo)系數(shù)權(quán)重分別為0.6、0.2和0.2，從而得出了綜合風(fēng)險(xiǎn)區(qū)劃指數(shù)D。按照綜合風(fēng)險(xiǎn)區(qū)劃指數(shù)D計(jì)算結(jié)果，將遼寧山洪災(zāi)害風(fēng)險(xiǎn)劃依次劃分為低風(fēng)險(xiǎn)、次低風(fēng)險(xiǎn)、中等風(fēng)險(xiǎn)、次高風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)5個(gè)等級(jí)。

表4 遼寧省暴雨誘發(fā)山洪承災(zāi)體災(zāi)損敏感系數(shù)Table 4 Assignment of sensitivity coefficient of disaster damage ofmountain floods induced by torrential rain

2.2 致災(zāi)、孕災(zāi)、易損性評(píng)估

2.2.1致災(zāi)危險(xiǎn)性評(píng)估

根據(jù)遼寧暴雨致災(zāi)危險(xiǎn)性區(qū)劃（圖5a），高危險(xiǎn)和次高危險(xiǎn)區(qū)主要位于遼寧南部地區(qū)的大連（莊河市、瓦房店市），遼寧中西部地區(qū)的葫蘆島市與朝陽(yáng)縣交界、錦州（北鎮(zhèn)市、黑山縣、義縣）、營(yíng)口市、岫巖縣，遼寧北部和東部地區(qū)的鐵嶺大部（鐵嶺市、鐵嶺縣、開(kāi)原市、西豐縣）、桓仁縣、寬甸縣。遼寧東部地區(qū)（丹東市、東港市）、遼寧中部地區(qū)（遼陽(yáng)地區(qū)）、遼寧西部地區(qū)（阜新與朝陽(yáng)交界）致災(zāi)危險(xiǎn)性也較高；而遼寧東部地區(qū)（撫順東部、本溪中西部）、遼寧西部地區(qū)（阜新、朝陽(yáng)西北部、葫蘆島市南部）、遼寧南部地區(qū)的大連為暴雨致災(zāi)低危險(xiǎn)性區(qū)，表明遼寧地區(qū)暴雨致災(zāi)危險(xiǎn)性較為分散。

圖5 遼寧省致災(zāi)暴雨危險(xiǎn)性區(qū)劃（a）、山洪溝孕災(zāi)環(huán)境敏感性區(qū)劃（b）、承災(zāi)體易損性區(qū)劃（c）Fig.5 Disaster-causing torrential rain hazard zoning（a），mountain torrent gully disaster-pregnant environment sensitivity zoning（b），and hazard-bearing body vulnerability zoning（c）in Liaoning province

2.2.2孕災(zāi)環(huán)境評(píng)估

由遼寧暴雨致洪孕災(zāi)環(huán)境因子溝口高程、溝床比降和河網(wǎng)密度可以看出（圖3），低值區(qū)地勢(shì)遼河流域（≤50 m）主要位于遼寧中部的沈陽(yáng)、遼陽(yáng)、盤錦等地。由于遼寧西、南、東北三面環(huán)山，山區(qū)的地勢(shì)起伏大，在溝谷內(nèi)山洪溝低洼地帶比較容易出現(xiàn)山洪災(zāi)害，從山洪溝的溝口位置高程來(lái)看，遼寧東南部、西南部地區(qū)的山洪溝溝口位置海拔高度較低，遼寧東北部山洪溝溝床比降較高，具有發(fā)生山洪的較好孕災(zāi)環(huán)境。河網(wǎng)密集區(qū)（圖3c）通常也多是山洪溝的匯水區(qū)，因此出現(xiàn)山洪的可能性較大［24］，遼寧東部的撫順東部和營(yíng)口蓋州和鞍山的岫巖縣河網(wǎng)密度較大，其中河網(wǎng)密度高值區(qū)（≥0.8 km·km-2）是遼寧重要河流如碧流河、英那河、大洋河、哨子河、大遼河、柴河、凡河、蘇子河等匯集區(qū)。從遼寧省山洪溝的孕災(zāi)敏感性區(qū)劃（圖5b）來(lái)看，山洪溝孕災(zāi)環(huán)境高敏感區(qū)域主要位于遼寧東南部和錦州的黑山縣一帶，遼寧西部的葫蘆島、遼寧北部的鐵嶺山區(qū)處于次高敏感區(qū)，遼寧西部的朝陽(yáng)地區(qū)、遼寧東部的撫順地區(qū)東南部新賓縣的山洪溝孕災(zāi)環(huán)境敏感性較低，其他地區(qū)山洪溝位于處于中等敏感區(qū)或次敏感區(qū)。

圖6 遼寧省暴雨誘發(fā)山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃（a）和1951—2013年山洪災(zāi)害次數(shù)分布（b）Fig.6 Risk division ofmountain torrent disasters induced by heavy rain（a）and statistical distribution map ofmountain torrent disasters from 1951 to 2013（b）in Liaoning province

2.2.3易損性評(píng)估

根據(jù)遼寧省承災(zāi)體風(fēng)險(xiǎn)暴露度分布可知（圖4a），遼寧東南部（鞍山市、岫巖縣、莊河市）、遼寧西部（葫蘆島、錦州北部）、遼寧北部（鐵嶺）部分地區(qū)的山洪溝有比較高的人口密度，鐵嶺大部分地區(qū)、阜新縣、岫巖縣的山洪溝有比較高耕地比例。由遼寧省承災(zāi)體的易損性風(fēng)險(xiǎn)區(qū)劃（圖5c）可知，高、中、低風(fēng)險(xiǎn)區(qū)非常分散，每個(gè)具有山洪溝的縣、市，均有高、中、低不同等級(jí)易損性風(fēng)險(xiǎn)區(qū)域。

2.3 暴雨誘發(fā)山洪災(zāi)害綜合風(fēng)險(xiǎn)區(qū)劃

由遼寧暴雨誘發(fā)山洪災(zāi)害綜合風(fēng)險(xiǎn)區(qū)劃可以看出（圖6a），遼寧東南部（大連北部、營(yíng)口東部、鞍山南部、丹東中西部、本溪東部）和遼寧西部（錦州北部）為高風(fēng)險(xiǎn)區(qū)，遼寧西部（葫蘆島）、遼寧北部（鐵嶺）地區(qū)處于次高風(fēng)險(xiǎn)區(qū)，朝陽(yáng)大部、撫順東南部等地為低風(fēng)險(xiǎn)區(qū)，遼寧其他地區(qū)的山洪溝分別處于次低風(fēng)險(xiǎn)或中等風(fēng)險(xiǎn)區(qū)。

以山洪溝為單位計(jì)量的風(fēng)險(xiǎn)區(qū)劃與1951—2013年發(fā)生的山洪災(zāi)害空間區(qū)域比較來(lái)看（圖6b），1951—2013年遼寧省發(fā)生山洪災(zāi)害高頻區(qū)大多位于遼寧東南部山區(qū)，這也是本次暴雨誘發(fā)洪災(zāi)害高風(fēng)險(xiǎn)或次高風(fēng)險(xiǎn)區(qū)，錦州北部、葫蘆島、鐵嶺等地山洪溝在歷史上發(fā)生過(guò)1—2次山洪災(zāi)害，為本次風(fēng)險(xiǎn)區(qū)劃次高風(fēng)險(xiǎn)或中等風(fēng)險(xiǎn)區(qū)。值得一提的是，本溪東部桓仁縣發(fā)生過(guò)山洪災(zāi)害的山洪溝與本次評(píng)估為次高風(fēng)險(xiǎn)或中風(fēng)險(xiǎn)區(qū)域基本一致。因此，通過(guò)對(duì)比結(jié)果來(lái)看，本次風(fēng)險(xiǎn)區(qū)劃為高風(fēng)險(xiǎn)區(qū)確實(shí)發(fā)生過(guò)山洪災(zāi)害災(zāi)害并且出現(xiàn)的頻率相對(duì)較高，有的山洪溝出現(xiàn)山洪災(zāi)害次數(shù)較少或者未出現(xiàn)山洪災(zāi)害，但風(fēng)險(xiǎn)區(qū)劃仍評(píng)估具有中風(fēng)險(xiǎn)或者次高風(fēng)險(xiǎn)的潛在風(fēng)險(xiǎn)，表明本次暴雨誘發(fā)山洪災(zāi)害綜合風(fēng)險(xiǎn)的結(jié)果較為可信。

3 結(jié)論與討論

（1）針對(duì)誘發(fā)山洪災(zāi)害從致災(zāi)因子危險(xiǎn)性、孕災(zāi)環(huán)境敏感性和承災(zāi)體易損性狀況，形成了一套基于GIS的遼寧省山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃評(píng)估方法。通過(guò)統(tǒng)計(jì)分析，以6 h暴雨作為遼寧省山洪致災(zāi)因子更為合適，進(jìn)一步構(gòu)建了綜合6 h分級(jí)暴雨強(qiáng)度的暴雨頻次作為暴雨致災(zāi)危險(xiǎn)性評(píng)價(jià)指標(biāo)；通過(guò)暴雨致災(zāi)危險(xiǎn)性指標(biāo)與孕災(zāi)環(huán)境、承載體暴露度等有關(guān)的空間指標(biāo)，進(jìn)行遼寧省山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估與區(qū)劃。與已有遼寧山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃研究成果相比，本次風(fēng)險(xiǎn)區(qū)劃綜合利用了暴雨致災(zāi)危險(xiǎn)性、孕災(zāi)環(huán)境的敏感性和承災(zāi)體的易損性，且以山洪溝為單元，區(qū)劃結(jié)果較為精細(xì)。

（2）遼寧省暴雨誘發(fā)山洪災(zāi)害綜合風(fēng)險(xiǎn)的高風(fēng)險(xiǎn)區(qū)域主要位于3個(gè)地區(qū)，遼寧東南部（大連北部、營(yíng)口東部、鞍山南部、丹東中西部、本溪東部）、遼寧西部（錦州北部，葫蘆島）和遼寧北部（鐵嶺）地區(qū)，與1951—2013年發(fā)生的山洪災(zāi)害高發(fā)區(qū)基本一致，本次風(fēng)險(xiǎn)區(qū)劃反映了遼寧省暴雨誘發(fā)山洪災(zāi)害的潛在風(fēng)險(xiǎn)，區(qū)劃結(jié)果可為遼寧省山洪災(zāi)害預(yù)報(bào)、預(yù)警及防治提供參考。

（3）本文應(yīng)用山洪風(fēng)險(xiǎn)普查資料，已經(jīng)對(duì)數(shù)據(jù)質(zhì)量進(jìn)行質(zhì)控，但仍存在個(gè)別數(shù)據(jù)不完整現(xiàn)象，可能會(huì)對(duì)暴雨誘發(fā)山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃的研究結(jié)果有影響。此外，山洪災(zāi)害風(fēng)險(xiǎn)影響因素較多［28-29］，風(fēng)險(xiǎn)區(qū)劃模型權(quán)重選擇存在一定的主觀性，也可能影響山洪災(zāi)害危險(xiǎn)性指標(biāo)計(jì)算和評(píng)估準(zhǔn)確度。以后將收集更加完善的承災(zāi)體暴露度資料和更合理的權(quán)重參數(shù)，進(jìn)一步提高風(fēng)險(xiǎn)區(qū)劃的準(zhǔn)確率。