王小亞，馬 諾，李?；?，毛子卿，張?jiān)苹?/p>

（1.伊吾縣氣象局，新疆 伊吾 839300；2.新疆氣象臺(tái)，新疆 烏魯木齊 830002；3.喀什地區(qū)氣象局，新疆 喀什 844000）

暴雨洪澇災(zāi)害每年汛期均有發(fā)生，是我國(guó)最主要的氣象災(zāi)害之一[1-2]。暴雨災(zāi)害風(fēng)險(xiǎn)評(píng)估始于20世紀(jì)30年代，由美國(guó)田納西河流域管理局提出洪水災(zāi)害評(píng)價(jià)理論方法[3]。王艷君等[4]利用1984—2012年中國(guó)暴雨洪澇災(zāi)害和社會(huì)經(jīng)濟(jì)數(shù)據(jù)，分析了災(zāi)害暴露度范圍，暴露度總體呈顯著增加趨勢(shì)，高暴露區(qū)域主要在沿海?。ㄊ校粸?zāi)害脆弱性高值區(qū)主要分布在長(zhǎng)江中游沿岸。謝五三等[5]利用氣象、水文、地理信息、社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)及歷史災(zāi)情數(shù)據(jù)，運(yùn)用TOPMODEL水文模型進(jìn)行洪水淹沒(méi)模擬得到不同重現(xiàn)期下洪水淹沒(méi)圖，再疊加不同承載體信息，得出不同的風(fēng)險(xiǎn)區(qū)劃圖。研究人員[6-13]基于風(fēng)險(xiǎn)評(píng)估原理，從致災(zāi)因子危險(xiǎn)性、孕災(zāi)環(huán)境敏感性、承災(zāi)體易損性和防災(zāi)減災(zāi)能力等方面，利用氣象資料、自然地理和社會(huì)經(jīng)濟(jì)等數(shù)據(jù)，采用加權(quán)綜合評(píng)價(jià)法和層次分析法，得出本區(qū)域暴雨災(zāi)害綜合風(fēng)險(xiǎn)區(qū)劃圖，對(duì)防災(zāi)減災(zāi)起到了一定的指導(dǎo)作用。

新疆雖然是干旱半干旱區(qū)，近36 a有18 a暴雨洪澇災(zāi)害的直接經(jīng)濟(jì)損失居所有氣象災(zāi)害損失的第一位[14]，具有暴雨洪澇種類多、直接經(jīng)濟(jì)損失明顯上升的規(guī)律[15-18]。而南疆地區(qū)由于干旱少雨且暴雨次數(shù)雖不多，但一次強(qiáng)降水往往能沖壞公路、鐵路、橋梁、阻斷交通、淹沒(méi)農(nóng)田，造成巨大損失。南疆洪水多以暴雨洪澇為主，隨著氣候變化和區(qū)域氣候變暖增濕，洪災(zāi)發(fā)生的不確定性和突發(fā)性也有所增加[19-20]。降水量可作為新疆易災(zāi)暴雨的風(fēng)險(xiǎn)區(qū)劃指標(biāo)[21]；暴雨洪澇災(zāi)害損失與年降水量、GDP指標(biāo)均有顯著正相關(guān)[22]；1961—2019年新疆暴雨山洪災(zāi)害強(qiáng)度呈上升趨勢(shì)，年災(zāi)損指數(shù)也呈線性增長(zhǎng)趨勢(shì)[23]。魏艷英等[24]利用降水?dāng)?shù)據(jù)和災(zāi)情資料等對(duì)南疆暴雨洪澇災(zāi)害進(jìn)行分析，得出暴雨日數(shù)大值區(qū)主要在克州及山區(qū)；喀什地區(qū)暴雨洪澇災(zāi)情發(fā)生頻次最多；阿克蘇地區(qū)受災(zāi)人口、農(nóng)作物受災(zāi)面積、直接經(jīng)濟(jì)損失均最多。上述研究從不同角度對(duì)新疆暴雨災(zāi)害進(jìn)行分析，為新疆暴雨洪澇災(zāi)害分析評(píng)估提供了參考，但從致災(zāi)因子和孕災(zāi)環(huán)境方面對(duì)南疆暴雨洪澇災(zāi)害風(fēng)險(xiǎn)區(qū)劃分析較少。為此，本文充分利用2010—2020年南疆氣象觀測(cè)站逐小時(shí)降水、南疆各縣（市）暴雨洪澇災(zāi)情數(shù)據(jù)，結(jié)合孕災(zāi)環(huán)境，基于GIS技術(shù)對(duì)南疆暴雨洪澇災(zāi)害風(fēng)險(xiǎn)進(jìn)行研究，以期提升對(duì)南疆暴雨洪澇災(zāi)害的客觀認(rèn)識(shí)及危險(xiǎn)性評(píng)價(jià)水平，為合理制定防災(zāi)減災(zāi)措施以及防洪避險(xiǎn)規(guī)劃提供科學(xué)參考和決策依據(jù)。

1 研究區(qū)概況

南疆由巴音郭楞蒙古自治州（以下簡(jiǎn)稱“巴州”）、阿克蘇地區(qū)、喀什地區(qū)、克孜勒蘇柯爾克孜自治州（以下簡(jiǎn)稱“克州”）、和田地區(qū)5個(gè)行政區(qū)組成（圖1）?？偯娣e為106.34×106hm2，約占新疆總面積的64%。南疆是典型的大陸性干旱氣候，塔里木盆地及周邊平原區(qū)是暖溫帶大陸性干旱氣候和綠洲氣候，山間盆地呈溫帶、中溫帶大陸性氣候，而周邊的高山區(qū)因氣候干燥寒冷，與盆地氣候相反。天山的高山帶是半干旱半濕潤(rùn)高寒山地氣候，帕米爾高原、昆侖山、阿爾金山的高山帶是干旱半干旱高寒山地氣候。南疆各地降水差異大，總的分布特點(diǎn)為北部多，南部少，山區(qū)多，平原少，盆地沙漠區(qū)更少，山區(qū)迎風(fēng)坡多，背風(fēng)坡少。

圖1 南疆5地（州）分布

2 數(shù)據(jù)與方法

2.1 資料

利用2010—2020年南疆5地（州）44個(gè)國(guó)家站和623個(gè)區(qū)域站降水資料；暴雨洪澇災(zāi)害來(lái)自中國(guó)氣象局災(zāi)情直報(bào)系統(tǒng)中的災(zāi)情數(shù)據(jù)，考慮到數(shù)據(jù)的連續(xù)性和完整性，對(duì)收集到的災(zāi)情數(shù)據(jù)進(jìn)行逐條整理和核對(duì)，剔除有誤的數(shù)據(jù)，有終報(bào)的以終報(bào)為準(zhǔn)，無(wú)終報(bào)數(shù)據(jù)以初報(bào)或續(xù)報(bào)為主（1條災(zāi)情統(tǒng)計(jì)1次）。

2.2 數(shù)據(jù)處理方法

（1）相關(guān)分析方法是Pearson相關(guān)，軟件是SPSS。

Pearson計(jì)算公式為：

式中，xi為變量x的第i個(gè)值，yi為變量y的第i個(gè)值為變量x的樣本均值，為變量y的樣本均值，n 為樣本容量。在給定顯著性水平下，對(duì)計(jì)算出的相關(guān)系數(shù)根據(jù)相關(guān)系數(shù)檢驗(yàn)表進(jìn)行顯著性檢驗(yàn)[22]。

（2）歸一化處理方法。

式中，Di是歸一化數(shù)值，Xi是樣本數(shù)據(jù)，Xmin和Xmax分別是樣本數(shù)據(jù)中的最小值和最大值。

（3）百分位數(shù)法，又稱為百分位數(shù)，是數(shù)據(jù)統(tǒng)計(jì)中一種常用方法。把一組統(tǒng)計(jì)數(shù)據(jù)按其數(shù)值從小到大順序排列，并按數(shù)據(jù)個(gè)數(shù)100等分。在第P個(gè)分界點(diǎn)（稱為百分位點(diǎn)）上的數(shù)值，稱為第P個(gè)百分位數(shù)（P=1，2，…，99）。在第P個(gè)分界點(diǎn)到第P+1個(gè)分界點(diǎn)之間的數(shù)據(jù)，稱為處于第P個(gè)百分位數(shù)。百分位數(shù)計(jì)算公式如下：

式中，Pm為第m個(gè)百分位數(shù)，N為總頻次，L為Pm所在組的下限，U為Pm所在組的上限，f為Pm所在組的次數(shù)，F(xiàn)h為小于L的累積次數(shù)，F(xiàn)n為大于U的累積次數(shù)，i為組距。

（4）自然斷點(diǎn)法，是一種地圖分級(jí)算法[25]。

2.3 結(jié)果分析

暴雨洪澇災(zāi)害風(fēng)險(xiǎn)是致災(zāi)因子易發(fā)性（主要指氣象因子）、孕災(zāi)環(huán)境敏感性（包括地形、水系、下墊面類型等方面）、承災(zāi)體易損性（人口密度和經(jīng)濟(jì)價(jià)值總量等方面）和防災(zāi)抗災(zāi)能力（資金投入的多少及交通的發(fā)達(dá)程度等方面）綜合作用的結(jié)果，其中致災(zāi)因子易發(fā)性和孕災(zāi)環(huán)境敏感性為自然屬性，共同構(gòu)成暴雨洪澇災(zāi)害的危險(xiǎn)性，他們決定了暴雨洪澇災(zāi)害發(fā)生的可能性和自然災(zāi)害發(fā)生的嚴(yán)重程度[26]。

2.3.1 致災(zāi)因子識(shí)別

選擇表達(dá)單站暴雨過(guò)程強(qiáng)度的指標(biāo)組成致災(zāi)因子庫(kù)，包括1 h最大降水量R1h、3 h最大降水量R3h、6 h最大降水量R6h、12 h最大降水量R12h、24 h最大降水量R24h、暴雨持續(xù)天數(shù)Iday。

由于暴雨洪澇災(zāi)害事件的災(zāi)情要素較多，利用數(shù)據(jù)相對(duì)完整的受災(zāi)人口、直接經(jīng)濟(jì)損失、農(nóng)作物受災(zāi)面積3項(xiàng)災(zāi)損指標(biāo)（一共有122個(gè)樣本，每個(gè)樣本均有受災(zāi)人口、直接損失、農(nóng)作物受災(zāi)面積和對(duì)應(yīng)逐1、3、6、12、24 h降水?dāng)?shù)據(jù)），使用Pearson相關(guān)系數(shù)計(jì)算方法，分別計(jì)算歸一化處理的災(zāi)損指標(biāo)與歸一化處理后的暴雨過(guò)程各個(gè)特征量的相關(guān)系數(shù)，選取通過(guò)顯著性檢驗(yàn)（α=0.05）的因子作為暴雨評(píng)估指標(biāo)。通過(guò)災(zāi)損指標(biāo)和致災(zāi)因子庫(kù)中各因子的關(guān)系分析，確定致災(zāi)因子（表1）。

表1 災(zāi)損指標(biāo)與暴雨評(píng)估指標(biāo)相關(guān)系數(shù)

由表1可知，暴雨評(píng)估指標(biāo)R6h、R12h、R24h均與各災(zāi)損指標(biāo)相關(guān)性較好，因此選取R6h、R12h、R24h作為南疆暴雨洪澇的致災(zāi)因子。

2.3.2 災(zāi)損指標(biāo)等級(jí)劃分

災(zāi)損指標(biāo)進(jìn)行歸一化處理后再根據(jù)百分位法，分別取50%、75%、90%，得到災(zāi)損指標(biāo)4級(jí)對(duì)應(yīng)閾值（表2），根據(jù)等級(jí)劃分標(biāo)準(zhǔn)得到南疆各縣（市）災(zāi)損指標(biāo)等級(jí)圖（圖2）。

表2 各災(zāi)損指標(biāo)等級(jí)劃分標(biāo)準(zhǔn)

從圖2可知，受災(zāi)人口特重的區(qū)域均在和田地區(qū)（洛浦縣、墨玉縣和于田縣）；受災(zāi)嚴(yán)重區(qū)域前三位是皮山縣、葉城縣、和田市；受災(zāi)較重區(qū)域前三位是岳普湖縣、策勒縣、烏什縣。直接經(jīng)濟(jì)損失特重區(qū)域?yàn)楹挽o縣、沙雅縣、烏什縣；嚴(yán)重區(qū)域前三位為洛浦縣、莎車縣、溫宿縣；較重區(qū)域前三位為阿克蘇市、阿克陶縣、拜城縣。農(nóng)作物受災(zāi)特重的區(qū)域?yàn)榘⒖颂K地區(qū)的沙雅縣、喀什地區(qū)的英吉沙縣和岳普湖縣；嚴(yán)重區(qū)域前三位為阿克陶縣、麥蓋提縣、皮山縣；較重區(qū)域前三位為拜城縣、和靜縣、和碩縣。近10 a和田地區(qū)的災(zāi)損指標(biāo)最嚴(yán)重。

圖2 南疆各縣暴雨洪澇災(zāi)損指標(biāo)等級(jí)

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

致災(zāi)危險(xiǎn)性評(píng)估主要考慮暴雨事件和孕災(zāi)環(huán)境，由年雨澇指數(shù)和暴雨孕災(zāi)環(huán)境影響兩部分組成。

（1）雨澇指數(shù)計(jì)算。

單站暴雨過(guò)程強(qiáng)度指數(shù)IR的計(jì)算：

式中，R6h、R12h、R24h均做歸一化處理，A、B、C為權(quán)重系數(shù)。結(jié)合區(qū)域自動(dòng)站資料，選取1 568個(gè)樣本數(shù)據(jù)，利用信息熵權(quán)法，得到權(quán)重系數(shù)分別為0.5、0.3、0.2，計(jì)算各站暴雨過(guò)程強(qiáng)度指數(shù)。累加當(dāng)年逐場(chǎng)暴雨過(guò)程強(qiáng)度值，得到年雨澇指數(shù)。該指數(shù)用來(lái)計(jì)算致災(zāi)因子危險(xiǎn)性指數(shù)。

（2）暴雨孕災(zāi)環(huán)境影響系數(shù)。

暴雨孕災(zāi)環(huán)境指暴雨影響下，對(duì)形成洪澇、泥石流、滑坡、城市內(nèi)澇等次生災(zāi)害起作用的自然環(huán)境。暴雨孕災(zāi)環(huán)境對(duì)暴雨成災(zāi)危險(xiǎn)性起擴(kuò)大或縮小作用。暴雨孕災(zāi)環(huán)境主要考慮地形、河網(wǎng)水系易發(fā)條件等。

對(duì)地形、水系易發(fā)條件進(jìn)行加權(quán)綜合來(lái)計(jì)算暴雨孕災(zāi)環(huán)境影響系數(shù)Ie（圖3a）。

式中，Icmin為區(qū)域內(nèi)的最小暴雨孕災(zāi)環(huán)境綜合指數(shù)，Icmax為區(qū)域內(nèi)的最大暴雨孕災(zāi)環(huán)境綜合指數(shù)，C為常數(shù)0.2。

暴雨孕災(zāi)環(huán)境綜合指數(shù)Ic（圖3b）。

圖3 Ie值（a）和Ic值（b）分布

式中，Wh、Wr分別為地形影響系數(shù)Ph、水系影響系數(shù)Pr的系數(shù)，按照南疆地形及水系特點(diǎn)，Wh和Wr分別定為0.5。

水系影響系數(shù)Pr主要考慮評(píng)估點(diǎn)距離水體的遠(yuǎn)近和所在區(qū)域的水網(wǎng)密度，分別用水體距離系數(shù)Prd和水網(wǎng)密度系數(shù)Prm表示。一般來(lái)說(shuō)，距離河流、湖泊、大型水庫(kù)等越近的地方，河網(wǎng)越密集的區(qū)域遭受洪澇災(zāi)害的風(fēng)險(xiǎn)越大。水系影響系數(shù)Pr為水體距離系數(shù)Prd和水網(wǎng)密度系數(shù)Prm的平均值。水體距離系數(shù)根據(jù)表3按距離水體（河流、湖泊、水庫(kù)）的遠(yuǎn)近取相應(yīng)的數(shù)值，水網(wǎng)密度系數(shù)為單位面積內(nèi)自然和人工河道的總長(zhǎng)度，按表4進(jìn)行賦值。

表3 水體面積、距離水體距離、水體距離系數(shù)

表4 水網(wǎng)密度與對(duì)應(yīng)的水網(wǎng)密度系數(shù)

2.3.4 暴雨洪澇災(zāi)害危險(xiǎn)性區(qū)劃

致災(zāi)因子危險(xiǎn)性指數(shù)=（1+暴雨孕災(zāi)環(huán)境影響系數(shù)）×年雨澇指數(shù)。

在收集整理資料的同時(shí)，綜合考慮影響南疆暴雨洪澇危險(xiǎn)性的致災(zāi)因子、暴雨孕災(zāi)環(huán)境的地形因子、水系因子各指標(biāo)，利用ArcGIS系統(tǒng)的提取、標(biāo)識(shí)等功能，將各個(gè)因子的值對(duì)應(yīng)到每個(gè)柵格，利用自然斷點(diǎn)法，得到每個(gè)致災(zāi)因子危險(xiǎn)性指數(shù)，將南疆暴雨洪澇災(zāi)害危險(xiǎn)性指數(shù)劃為4個(gè)等級(jí)：低危險(xiǎn)區(qū)、中低危險(xiǎn)區(qū)、中高危險(xiǎn)區(qū)、高危險(xiǎn)區(qū)（圖4）。

圖4 南疆暴雨洪澇災(zāi)害致災(zāi)因子風(fēng)險(xiǎn)區(qū)劃

南疆暴雨洪澇災(zāi)害風(fēng)險(xiǎn)整體分布呈現(xiàn)北部高于南部，西部高于東部，山區(qū)高于平原。暴雨洪澇災(zāi)害風(fēng)險(xiǎn)高的前三位是阿克蘇地區(qū)烏什縣英阿瓦提北村、和田地區(qū)于田縣蘭干鄉(xiāng)昆侖渠首村，奧依托格拉克鄉(xiāng)山洪村和沙漠玫瑰基地，皮什蓋村和阿克蘇地區(qū)拜城縣的大橋鄉(xiāng)牧場(chǎng)3隊(duì)（災(zāi)損指標(biāo)等級(jí)對(duì)應(yīng)）；中高風(fēng)險(xiǎn)的前三位是阿克蘇地區(qū)拜城縣布隆鄉(xiāng)8大隊(duì)6小隊(duì)，賽里木鎮(zhèn)，蘇杭村、溫巴什鄉(xiāng)3大隊(duì)、和田地區(qū)策勒縣努爾水庫(kù)；中低風(fēng)險(xiǎn)的區(qū)域前三位是和田地區(qū)洛浦縣阿其克公路到石灰礦岔路口、墨玉縣英阿瓦提村和阿克蘇地區(qū)溫宿縣古勒阿瓦提鄉(xiāng)四大隊(duì)；暴雨洪澇災(zāi)害風(fēng)險(xiǎn)低的地區(qū)為巴州盆地（輪臺(tái)縣、尉犁縣、和碩縣）、和田地區(qū)中西部、喀什地區(qū)南部山區(qū)。

2.3.5 2021年災(zāi)害評(píng)估結(jié)果

基于上述風(fēng)險(xiǎn)區(qū)劃成果，統(tǒng)計(jì)2021年126次暴雨洪澇災(zāi)情數(shù)據(jù)（表5），按災(zāi)損指標(biāo)閾值進(jìn)行等級(jí)劃分（表6），對(duì)2021年南疆暴雨洪澇災(zāi)害進(jìn)行了評(píng)估（圖5）。

圖5 2021年南疆暴雨洪澇災(zāi)害評(píng)估結(jié)果

表5 2021年南疆地區(qū)出現(xiàn)暴雨洪澇災(zāi)害時(shí)間、地點(diǎn)和災(zāi)情

和田地區(qū)洛浦縣局地暴雨洪澇災(zāi)害等級(jí)為高危險(xiǎn)區(qū)，與表6相對(duì)應(yīng)，其直接經(jīng)濟(jì)損失、受災(zāi)人口、農(nóng)作物受災(zāi)面積均為最嚴(yán)重，進(jìn)一步驗(yàn)證了前述風(fēng)險(xiǎn)區(qū)劃的結(jié)果。于田縣局地暴雨洪澇災(zāi)害等級(jí)為高危險(xiǎn)區(qū)，但災(zāi)損指標(biāo)等級(jí)不高，說(shuō)明于田縣孕災(zāi)環(huán)境較高。和田地區(qū)均處于中高風(fēng)險(xiǎn)區(qū)，2021年夏季和田地區(qū)大部出現(xiàn)了極端降水。洛浦、墨玉、和田3站日降水量均為百年一遇，打破測(cè)站日降水量歷史極值并超過(guò)其年平均降水量。暴雨造成城市積澇并引發(fā)山洪、泥石流地質(zhì)災(zāi)害，對(duì)交通、水利設(shè)施、農(nóng)業(yè)生產(chǎn)和居民生活等造成不利影響，直接損失達(dá)1億元。巴州和靜縣的孕災(zāi)環(huán)境風(fēng)險(xiǎn)等級(jí)為中高危險(xiǎn)區(qū)，實(shí)際業(yè)務(wù)中應(yīng)引起關(guān)注。阿克蘇地區(qū)溫宿縣局地為中高危險(xiǎn)區(qū)，農(nóng)作物受災(zāi)等級(jí)較重，沙雅縣雖然處于低危險(xiǎn)區(qū)，但一出現(xiàn)暴雨洪澇災(zāi)害，其經(jīng)濟(jì)損失和農(nóng)作物受災(zāi)面積的等級(jí)較重，在實(shí)際業(yè)務(wù)中需進(jìn)一步關(guān)注。

表6 2021年南疆地區(qū)暴雨洪澇災(zāi)損指標(biāo)等級(jí)

3 結(jié)論與討論

利用近10 a南疆地區(qū)暴雨洪澇災(zāi)情數(shù)據(jù)，根據(jù)百分位法得到各縣（市）的災(zāi)損指標(biāo)、根據(jù)相關(guān)系數(shù)確定了致災(zāi)因子、信息熵權(quán)法確定了權(quán)重系數(shù)?；贏rcGIS技術(shù)，采用重采樣、鄰域分析和重分類等方法結(jié)合地形、水系等環(huán)境因子得到暴雨洪澇的區(qū)劃圖。主要結(jié)論如下：

（1）南疆暴雨洪澇災(zāi)害的致災(zāi)因子為R6h、R12h、R24h；南疆各縣（市）受災(zāi)人口特重區(qū)域在和田地區(qū)洛浦縣、墨玉縣和于田縣；嚴(yán)重區(qū)域在皮山縣、葉城縣、和田市。直接經(jīng)濟(jì)損失特重區(qū)域在和靜縣、沙雅縣、烏什縣；嚴(yán)重區(qū)域在洛浦、莎車縣、溫宿縣。農(nóng)作物受災(zāi)特重區(qū)域在阿克蘇地區(qū)在沙雅縣、喀什地區(qū)的英吉沙縣和岳普湖縣；嚴(yán)重區(qū)域在阿克陶縣、麥蓋提縣、皮山縣。

（2）南疆暴雨災(zāi)害風(fēng)險(xiǎn)區(qū)劃分為低風(fēng)險(xiǎn)、中低風(fēng)險(xiǎn)、中高風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)4個(gè)等級(jí)。風(fēng)險(xiǎn)高的區(qū)域主要集中在和田地區(qū)于田縣南部山區(qū)、阿克蘇地區(qū)西部北部山區(qū)；風(fēng)險(xiǎn)低的地區(qū)巴州盆地東部，和田地區(qū)中西部，喀什地區(qū)南部山區(qū)。

（3）通過(guò)對(duì)2021年南疆暴雨洪澇災(zāi)害評(píng)估，進(jìn)一步驗(yàn)證了本文風(fēng)險(xiǎn)區(qū)劃成果，也為今后開展新疆暴雨洪澇災(zāi)害風(fēng)險(xiǎn)預(yù)警業(yè)務(wù)提供了參考。