孫 濤

(大連金普新區(qū)農(nóng)業(yè)農(nóng)村發(fā)展服務(wù)中心，遼寧 大連 116199)

0 引 言

在下墊面和降雨因素共同作用下引發(fā)的山丘區(qū)溪溝間暴漲暴落的地表徑流即為山洪，通常存在沖刷破壞力強(qiáng)、流速大、覆蓋范圍廣等典型特征，極易造成泥石流、滑坡等災(zāi)害。由于預(yù)防和預(yù)測(cè)較為困難使得山洪災(zāi)害易產(chǎn)生人員傷亡，對(duì)人身健康和國(guó)民經(jīng)濟(jì)造成重大威脅。實(shí)踐表明，山丘區(qū)的下墊面環(huán)境和降雨條件為驅(qū)動(dòng)山洪災(zāi)害的關(guān)鍵要素，一般易引發(fā)山洪的小流域其土壤質(zhì)地不利于水分吸收且比降較大[1-6]。采取科學(xué)有效的減災(zāi)措施和管理政策重要基礎(chǔ)是對(duì)風(fēng)險(xiǎn)的準(zhǔn)確評(píng)估，通過(guò)對(duì)山洪災(zāi)害驅(qū)動(dòng)因子的分析開(kāi)展系統(tǒng)性、前瞻性的風(fēng)險(xiǎn)評(píng)價(jià)可大大提升人類(lèi)抵御自然災(zāi)害的能力，為構(gòu)建綜合風(fēng)險(xiǎn)管理體系及促進(jìn)區(qū)域可持續(xù)發(fā)展提供重要依據(jù)，例如史培軍等從人為、環(huán)境和自然3個(gè)方面分解了區(qū)域?yàn)?zāi)害系統(tǒng)[7]。

近年來(lái)，將山洪災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)與GIS地理信息系統(tǒng)相結(jié)合的研究越來(lái)越多，如羅日洪等通過(guò)構(gòu)建曹江上游小流域山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估模型獲取了區(qū)域風(fēng)險(xiǎn)區(qū)劃圖；靳夢(mèng)等從山洪災(zāi)害易損度和危險(xiǎn)度的角度提取評(píng)價(jià)因子，利用GIS軟件疊置處理功能獲取山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃圖[8-9]。一般情況下山丘區(qū)為山洪災(zāi)害易發(fā)區(qū)，所以山洪災(zāi)害重點(diǎn)防治區(qū)為風(fēng)險(xiǎn)評(píng)估范圍。另外，現(xiàn)有研究還未形成較為成熟、完善的山洪災(zāi)害評(píng)估體系，且涉及山洪災(zāi)害形成原因的分析不夠系統(tǒng)全面。文章選取大連市山洪災(zāi)害防治區(qū)為研究對(duì)象，從人為、環(huán)境和自然災(zāi)害系統(tǒng)選擇土地利用類(lèi)型、河網(wǎng)密度、地形指數(shù)、年最大6h降雨量等10項(xiàng)因子構(gòu)建風(fēng)險(xiǎn)評(píng)價(jià)體系，采用Arc GIS軟件和層次分析法疊加、分配各項(xiàng)因子及其權(quán)重，從而獲取區(qū)域山洪災(zāi)害風(fēng)險(xiǎn)等級(jí)。

1 區(qū)域概況

大連市地處黃渤海之濱，介于E120°58′-123°31′、N38°43′-40°10′之間，大連市總面積12573.85km2，地勢(shì)呈南窄北寬、南低北高的特征，地貌形態(tài)復(fù)雜且大部分為低緩丘陵和山地，一些山間谷底及河流入海處零星分布有平原低地，海蝕地貌和喀斯特地貌比較發(fā)育且?guī)r溶地形隨處可見(jiàn)。

研究區(qū)具有溫暖帶大陸性季風(fēng)氣候，近鄰海洋海岸線較長(zhǎng)且氣候適宜，年均氣溫10.5℃，年均降水量550-950mm，夏季多暴雨且集中全年60%-70%的降水，日照時(shí)數(shù)2500-2800h。黃海和渤海為該地區(qū)的主要水系，贊子河、大沙河、清水河、莊河、英那河及碧流河等為注入黃海了的主要河流，且以碧流河最大；三十里堡河、李官村河及復(fù)州河等為注入渤海的較大河流，其它小河近200余條。統(tǒng)計(jì)資料顯示，大連市淡水、地下水和地表水分別為37.86億m3、8.84億m3、34.2億m3，其中重復(fù)水為5.8億m3。近年來(lái)，特殊的地形條件、復(fù)雜的地質(zhì)構(gòu)造及降水環(huán)境使得大連山洪災(zāi)害頻發(fā)，對(duì)當(dāng)?shù)亟煌?、通訊和生產(chǎn)生活產(chǎn)生嚴(yán)重的影響，為進(jìn)一步提升區(qū)域山洪災(zāi)害科學(xué)管理水平，有必要采取科學(xué)的方法評(píng)估山洪災(zāi)害風(fēng)險(xiǎn)狀況。

2 研究方法

2.1 山洪災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)體系

山洪屬于一種自然災(zāi)害對(duì)山丘地區(qū)的經(jīng)濟(jì)發(fā)展和生命安全產(chǎn)生巨大威脅，區(qū)域?yàn)?zāi)害系統(tǒng)理論認(rèn)為環(huán)境、自然、人為災(zāi)害風(fēng)險(xiǎn)的疊加即為山洪災(zāi)害風(fēng)險(xiǎn)，因此考慮從此方面構(gòu)造評(píng)估體系。

致災(zāi)因子的風(fēng)險(xiǎn)性即為自然災(zāi)害的風(fēng)險(xiǎn)性，由于大連市具有明顯的山丘區(qū)域差異，高強(qiáng)度、短歷時(shí)暴雨為造成山洪災(zāi)害的關(guān)鍵因子，選取年最大24h、年最大6h、年最大10min和多年平均降雨量4個(gè)指標(biāo)突出暴雨致災(zāi)因子的空間差異性；區(qū)域下墊面環(huán)境的穩(wěn)定性即為環(huán)境致災(zāi)風(fēng)險(xiǎn)性，由于降雨匯水。地形與山洪災(zāi)害之間的聯(lián)系較為密切，選擇的環(huán)境災(zāi)害系統(tǒng)指標(biāo)包括河網(wǎng)密度、地形指數(shù)和地形坡度；承載體的脆弱性即為人為災(zāi)害風(fēng)險(xiǎn)性，經(jīng)濟(jì)發(fā)展水平較高和人口密集區(qū)的脆弱性較低，選取人為災(zāi)害系統(tǒng)指標(biāo)包括土地利用類(lèi)型、人口分布密度和人均GDP。

根據(jù)選取的10項(xiàng)指標(biāo)構(gòu)建包括指標(biāo)層、準(zhǔn)則層、目標(biāo)層的山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估體系如表1所示，且山洪災(zāi)害風(fēng)險(xiǎn)性高低與各指標(biāo)因子的數(shù)值大小呈正相關(guān)性。依據(jù)相關(guān)研究資料和自然間斷法，將山洪災(zāi)害風(fēng)險(xiǎn)等級(jí)劃分為1-5級(jí)，風(fēng)險(xiǎn)程度隨等級(jí)的增大而增加，各風(fēng)險(xiǎn)因子分界值如表2所示。

表1 大連市山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估體系

表2 山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估等級(jí)

2.2 數(shù)據(jù)的處理

根據(jù)遼寧省中小流域暴雨圖集中的線矢量插值獲取年最大24h、6h、10min和多年平均降雨量數(shù)據(jù)，由地理空間數(shù)據(jù)云平臺(tái)下載精度為90m的數(shù)字高程模型DEM數(shù)據(jù)；全國(guó)山洪災(zāi)害項(xiàng)目組下發(fā)獲取土地利用格柵數(shù)據(jù)和中小河流線矢量，精度為30m；中科院資源環(huán)境科學(xué)數(shù)據(jù)中心提供人口密度和人均GDP數(shù)據(jù)，通過(guò)相關(guān)公式計(jì)算獲取河網(wǎng)密度、地形指數(shù)和地形起伏度。其中，地形起伏度V利用Arc GIS格柵計(jì)算器和DEM數(shù)據(jù)源獲取，用于描述單位相鄰格柵高程的起伏程度，數(shù)學(xué)表達(dá)式為：

V=Hmax-Hmin

(1)

式中：Hmin、Hmax為鄰域內(nèi)高程最小、最大值；

地形指數(shù)TI體現(xiàn)了在空間分布上的流域飽和缺水量，用于描述地表產(chǎn)流受地形的影響作用，利用Arc GIS格柵處理工具和DEM數(shù)據(jù)源獲取，其計(jì)算式為：

TI=lnA/S

(2)

式中：S、A為任一點(diǎn)i的坡度及其在徑流坡面上的單寬集水面積。

單位相鄰面積a上的平均河長(zhǎng)l即為河網(wǎng)密度D，若劃分為n個(gè)河段或河流的流域面積為A、總長(zhǎng)度為L(zhǎng)，則l=L/n、a=A/n。運(yùn)用Arc GIS線密度分析工具和以中小河流矢量數(shù)據(jù)源獲取河網(wǎng)密度D，其表達(dá)式為：

(3)

確保各項(xiàng)參評(píng)因子的分辨率和投影坐標(biāo)系的統(tǒng)一為疊置分析的基本條件，將所有數(shù)據(jù)利用Arc GIS中的重采樣工具和投影格柵投影值分辨率為1km的WGS_1984_UTM_ Zone_50N。依據(jù)以上處理方法和數(shù)據(jù)來(lái)源，提取大連市各區(qū)縣山洪災(zāi)害風(fēng)險(xiǎn)各因子初始值，如表3所示。

表3 研究區(qū)山洪災(zāi)害風(fēng)險(xiǎn)評(píng)判因子值

2.3 確定指標(biāo)權(quán)重

根據(jù)AHP層次分析法獲取各參評(píng)因子權(quán)重，依據(jù)區(qū)域山洪災(zāi)害實(shí)際情況選擇評(píng)價(jià)因子構(gòu)建多要素、多目標(biāo)層次結(jié)構(gòu)；然后采用專(zhuān)家打分法對(duì)兩兩要素的重要度評(píng)分，經(jīng)適當(dāng)調(diào)整構(gòu)造判斷矩陣，經(jīng)一致性檢驗(yàn)確定各評(píng)價(jià)指標(biāo)權(quán)重見(jiàn)表1，其中CI值為0.0185滿足<0.1的要求。強(qiáng)度大、短歷時(shí)暴雨為驅(qū)動(dòng)山洪災(zāi)害的主要因子，且較大地面起伏區(qū)很容易產(chǎn)生洪澇，對(duì)居民生命安全和經(jīng)濟(jì)發(fā)展產(chǎn)生嚴(yán)重威脅，所以權(quán)重較大的指標(biāo)有人口密度、地形起伏度和年最大6h降雨量。

2.4 風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果

將10項(xiàng)指標(biāo)數(shù)據(jù)利用Arc GIS軟件提取，經(jīng)柵格化處理后分類(lèi)各個(gè)評(píng)價(jià)指標(biāo)，然后對(duì)山洪災(zāi)害風(fēng)險(xiǎn)等級(jí)運(yùn)用自然間斷法分為1-5級(jí)，各區(qū)域風(fēng)險(xiǎn)等級(jí)如表4所示，最后，經(jīng)疊加計(jì)算最終確定風(fēng)險(xiǎn)評(píng)價(jià)值。

表4 自然、環(huán)境和人工災(zāi)害系統(tǒng)風(fēng)險(xiǎn)評(píng)價(jià)

續(xù)表4 自然、環(huán)境和人工災(zāi)害系統(tǒng)風(fēng)險(xiǎn)評(píng)價(jià)

采用正態(tài)分布取值的方法獲取山洪災(zāi)害風(fēng)險(xiǎn)級(jí)別的閾值，將風(fēng)險(xiǎn)級(jí)別劃分為很高、較高、中等、較低、很低5個(gè)等級(jí)，所對(duì)應(yīng)的風(fēng)險(xiǎn)評(píng)價(jià)最終值為100%-85%、85%-60%、60%-40%、40%-15%、15%-0%，最終的評(píng)價(jià)結(jié)果見(jiàn)表5。結(jié)果顯示，大連市山洪災(zāi)害風(fēng)險(xiǎn)以中等-較高-很高為主，由于大連市防洪災(zāi)害防治的財(cái)力、物力和人力有限，為進(jìn)一步提升山洪災(zāi)害管理效率應(yīng)優(yōu)先考慮較高和很高風(fēng)險(xiǎn)區(qū)，同時(shí)兼顧中等和較低風(fēng)險(xiǎn)區(qū)域。

表5 大連市山洪災(zāi)害風(fēng)險(xiǎn)分級(jí)及面積占比

2.5 驗(yàn)證分析

驗(yàn)證分析數(shù)據(jù)來(lái)源于大連市歷史山洪災(zāi)害點(diǎn)相關(guān)數(shù)據(jù)，在風(fēng)險(xiǎn)等級(jí)和點(diǎn)工具上利用Arc GIS軟件、Arc toolbox工具箱標(biāo)注提取，歷史山洪災(zāi)害點(diǎn)統(tǒng)計(jì)占比與數(shù)量如表6所示。根據(jù)表6可知，歷史上發(fā)生很高-較高-中等山洪災(zāi)害風(fēng)險(xiǎn)級(jí)別的區(qū)域占比達(dá)到91.26%，由此表明評(píng)價(jià)結(jié)果具有較高的準(zhǔn)確率，分析結(jié)果與區(qū)域歷史統(tǒng)計(jì)數(shù)據(jù)保持良好的一致性，可為未來(lái)開(kāi)展區(qū)域山洪災(zāi)害防治提供科學(xué)的依據(jù)和指導(dǎo)。

表6 各個(gè)風(fēng)險(xiǎn)級(jí)別歷史災(zāi)害點(diǎn)占比

3 結(jié) 論

1)從人為、自然和環(huán)境3大子系統(tǒng)選擇山洪災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)因子，結(jié)合災(zāi)害系統(tǒng)理論構(gòu)建綜合評(píng)價(jià)體系，并對(duì)各指標(biāo)權(quán)重利用AHP法求解，選取指標(biāo)信息提取空間為Arc GIS系統(tǒng)軟件，從而獲取大連市各區(qū)域山洪風(fēng)險(xiǎn)等級(jí)。結(jié)果發(fā)現(xiàn)，達(dá)到很高-較高-中等山洪風(fēng)險(xiǎn)等級(jí)的區(qū)域占比超過(guò)90%，在空間分布上揭示了各區(qū)域山洪風(fēng)險(xiǎn)的差異性特征。

2)根據(jù)歷史統(tǒng)計(jì)資料驗(yàn)證評(píng)價(jià)結(jié)果，結(jié)果顯示評(píng)價(jià)結(jié)果具有較高的準(zhǔn)確度，風(fēng)險(xiǎn)分果與區(qū)域歷史統(tǒng)計(jì)數(shù)據(jù)保持良好的一致性，較為客觀、準(zhǔn)確的揭示區(qū)域山洪災(zāi)害分布規(guī)律，可為未來(lái)開(kāi)展區(qū)域山洪災(zāi)害防治提供科學(xué)的依據(jù)和指導(dǎo)。

3)山洪災(zāi)害的發(fā)生受多個(gè)驅(qū)動(dòng)因素的影響，且各因子影響程度和作用角度存在較大差異，較為準(zhǔn)確、系統(tǒng)的分析山洪災(zāi)害風(fēng)險(xiǎn)還存在相應(yīng)的難度，為提高評(píng)價(jià)結(jié)果的參考價(jià)值和針對(duì)性，未來(lái)應(yīng)繼續(xù)完善評(píng)價(jià)模型和優(yōu)化評(píng)價(jià)體系。