哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)空間演變模式

陳鵬張繼權(quán)孫瀅悅張立峰劉家福

摘要：[JP+1]由于洪災(zāi)突發(fā)性強(qiáng)、破壞性大、影響范圍廣，給居民生命安全、財(cái)產(chǎn)等造成重大損失，為了更好的預(yù)防洪災(zāi)及了解未來(lái)洪災(zāi)風(fēng)險(xiǎn)演變模式，利用自然災(zāi)害風(fēng)險(xiǎn)形成四因子理論，選取哈爾濱市區(qū)氣象、水文、經(jīng)濟(jì)及社會(huì)等18個(gè)指標(biāo)，構(gòu)建哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)模型，對(duì)哈爾濱市區(qū)2000年-2009年洪災(zāi)風(fēng)險(xiǎn)區(qū)劃，并通過(guò)GIS技術(shù)對(duì)哈爾濱市區(qū)2000年-2009年洪災(zāi)風(fēng)險(xiǎn)區(qū)劃分析，歸類(lèi)風(fēng)險(xiǎn)空間演變模式。結(jié)果為：哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)空間轉(zhuǎn)移升高模式包括邊緣接觸模式、擺動(dòng)模式、內(nèi)部突增模式、跨越模式4種；洪災(zāi)風(fēng)險(xiǎn)轉(zhuǎn)移降低模式包括邊緣退縮模式、片狀降低模式、內(nèi)部突降模式3種。研究結(jié)果可為減少洪災(zāi)對(duì)城市居民帶來(lái)的損失及城市防洪規(guī)劃提供決策依據(jù)。

關(guān)鍵詞：風(fēng)險(xiǎn)分析；演變模式；格網(wǎng)GIS；洪災(zāi)風(fēng)險(xiǎn)；哈爾濱

中圖分類(lèi)號(hào)：X45文獻(xiàn)標(biāo)志碼：A文章編號(hào)：16721683（2016）06002706

洪災(zāi)是較為常見(jiàn)的一種自然災(zāi)害，其特點(diǎn)是高發(fā)性、強(qiáng)破壞性及影響范圍廣，給人們生活造成嚴(yán)重影響，嚴(yán)重制約了我國(guó)城市可持續(xù)發(fā)展。哈爾濱市是我國(guó)東北北部的政治、文化、經(jīng)濟(jì)中心，也是國(guó)內(nèi)省轄市中陸地管轄面積最大、管轄總?cè)丝谳^多的特大城市。因哈爾濱市特殊的地理環(huán)境影響，每年受洪水災(zāi)害影響程度較大，僅2013年哈爾濱市因洪澇災(zāi)害造成經(jīng)濟(jì)損失達(dá)66億元，因此，對(duì)哈爾濱市洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)與演變模式研究是非常必要的，若采用科學(xué)的方法及遵循自然規(guī)律正確的對(duì)研究區(qū)洪災(zāi)風(fēng)險(xiǎn)及演變模式研究，并以洪災(zāi)風(fēng)險(xiǎn)空間演變模式適時(shí)修訂防洪規(guī)劃，可減少城市居民的生命、財(cái)產(chǎn)損失，同時(shí)對(duì)城市經(jīng)濟(jì)及可持續(xù)發(fā)展起到保障作用。

目前，有不少學(xué)者對(duì)洪災(zāi)風(fēng)險(xiǎn)區(qū)劃及評(píng)價(jià)進(jìn)行了深入的研究[12]，一部分學(xué)者研究著重從三個(gè)方面進(jìn)行，即致災(zāi)因子的危險(xiǎn)性、承災(zāi)體的暴露性與脆弱性三個(gè)方面來(lái)選取指標(biāo)，并構(gòu)建洪災(zāi)評(píng)價(jià)模型，該種評(píng)價(jià)方法在指標(biāo)選取上忽略了人對(duì)災(zāi)害防御方面的內(nèi)容[34]；另一部分學(xué)者則是先建立致災(zāi)因子強(qiáng)度與承災(zāi)體的脆弱性對(duì)應(yīng)關(guān)系，后通過(guò)擬合風(fēng)險(xiǎn)曲線(xiàn)，建立綜合災(zāi)損模型并計(jì)算洪災(zāi)造成的損失，以此作為洪災(zāi)風(fēng)險(xiǎn)，該種評(píng)價(jià)方法中的風(fēng)險(xiǎn)曲線(xiàn)擬合需通過(guò)大量的災(zāi)害數(shù)據(jù)及調(diào)查數(shù)據(jù)結(jié)合，所以在研究中相對(duì)獲取數(shù)據(jù)較為困難[58]，但此種方法對(duì)于洪災(zāi)保險(xiǎn)方面研究具有重要意義[9]。而在近些年出現(xiàn)一種從綜合角度出發(fā)，對(duì)洪災(zāi)風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)方法。其中具有代表性的是從自然災(zāi)害風(fēng)險(xiǎn)形成機(jī)理出發(fā)，以自然災(zāi)害風(fēng)險(xiǎn)四因子為基礎(chǔ)，即從致災(zāi)因子的危險(xiǎn)性、承災(zāi)體的暴露性與脆弱性及防災(zāi)減災(zāi)能力四個(gè)方面選取評(píng)價(jià)指標(biāo)，綜合構(gòu)建洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)模型，實(shí)現(xiàn)了洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)[1011]。但已有研究中僅對(duì)洪災(zāi)風(fēng)險(xiǎn)進(jìn)行了評(píng)價(jià)，并未考慮各年間風(fēng)險(xiǎn)轉(zhuǎn)移規(guī)律。本研究在洪災(zāi)風(fēng)險(xiǎn)區(qū)劃基礎(chǔ)上，對(duì)哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)空間演變模式進(jìn)行了研究，研究結(jié)果可為城市防洪提供決策依據(jù)。

1研究區(qū)概況

哈爾濱市地處E125°42′-130°10′、N44°04′-46°40′，全市土地面積為531萬(wàn)km2，轄8區(qū)10縣（市）?？?cè)丝? 0635萬(wàn)人，其中市區(qū)人口位587.9萬(wàn)人。境內(nèi)的大小河流均屬于松花江水系和牡丹江水系，降雨多集中在6月-9月，易發(fā)洪水災(zāi)害。如在1932年曾發(fā)生過(guò)較大洪水災(zāi)害，哈爾濱整個(gè)城市被淹，造成幾十萬(wàn)人受災(zāi)。

2研究方法與數(shù)據(jù)來(lái)源

2.1研究方法

[BT4]2.1.1自然災(zāi)害風(fēng)險(xiǎn)指數(shù)法

自然災(zāi)害風(fēng)險(xiǎn)是指未來(lái)一段時(shí)間內(nèi)可能達(dá)到的災(zāi)害程度及其發(fā)生的可能性。公認(rèn)的自然災(zāi)害風(fēng)險(xiǎn)是由四因子綜合構(gòu)成，即：致災(zāi)因子危險(xiǎn)性、承災(zāi)體的暴露性及脆弱性、防災(zāi)減災(zāi)能力綜合作用的結(jié)果[12]，可用下面公式表達(dá)：

自然災(zāi)害風(fēng)險(xiǎn)指數(shù)=危險(xiǎn)性×暴露性×脆弱性×防災(zāi)減災(zāi)能力1[JY]（1）

2.1.2網(wǎng)格GIS技術(shù)

[JP+2]網(wǎng)格GIS技術(shù)是由現(xiàn)有的網(wǎng)格技術(shù)、空間信息技術(shù)、地理信息技術(shù)綜合形成的空間信息管理與處理環(huán)境，是提高研究精度的重要工具[13]。目前，常用的網(wǎng)格劃分方法為規(guī)則網(wǎng)格與不規(guī)則網(wǎng)格兩種，網(wǎng)格精度（網(wǎng)格大?。┦苎芯?jī)?nèi)容影響，樣式比較靈活，是現(xiàn)代小尺度研究中常用的方法。本文將研究區(qū)按規(guī)則網(wǎng)格進(jìn)行劃分，網(wǎng)格大小為5 km×5 km，網(wǎng)格劃分?jǐn)?shù)量為 2 260個(gè)。其中各研究指標(biāo)數(shù)據(jù)依據(jù)網(wǎng)格展布原理展布到研究區(qū)網(wǎng)格中。

2.2數(shù)據(jù)來(lái)源與處理

[JP+2]研究數(shù)據(jù)來(lái)源于哈爾濱市氣象局，主要為歷年降雨數(shù)據(jù)（2000年-2009年逐年降雨量），利用數(shù)學(xué)中頻率計(jì)算方法計(jì)算得到研究區(qū)降雨頻率；研究區(qū)地形數(shù)據(jù)空間分辨率90 m，來(lái)源于地球系統(tǒng)科學(xué)數(shù)據(jù)共享網(wǎng)STRM系統(tǒng)，包括DRM與遙感影像，其中下墊面基礎(chǔ)數(shù)據(jù)是利用GIS數(shù)字化技術(shù)轉(zhuǎn)化為矢量數(shù)據(jù)，作為研究區(qū)基礎(chǔ)地理信息數(shù)據(jù)；生成的規(guī)則網(wǎng)格數(shù)據(jù)來(lái)源于ARCGIS軟件自動(dòng)生成，生成網(wǎng)格為矢量數(shù)據(jù)，具有對(duì)應(yīng)的網(wǎng)格編號(hào)及其相應(yīng)屬性，利用ARCGIS軟件中的疊加、相交等功能實(shí)現(xiàn)基礎(chǔ)地理信息數(shù)據(jù)與研究區(qū)網(wǎng)格數(shù)據(jù)關(guān)聯(lián)，使網(wǎng)格數(shù)據(jù)獲取下墊面基礎(chǔ)地理信息，以備后續(xù)研究使用；其它社會(huì)、經(jīng)濟(jì)等數(shù)據(jù)主要來(lái)源于《黑龍江省統(tǒng)計(jì)年鑒》，主要包括居民數(shù)量、人口密度、工業(yè)總產(chǎn)值、農(nóng)林牧副漁總產(chǎn)量、小企業(yè)數(shù)量、受災(zāi)面積等，數(shù)據(jù)年份為2000年-2009年，按年份分別關(guān)聯(lián)到每個(gè)網(wǎng)格中，為自然災(zāi)害風(fēng)險(xiǎn)指數(shù)計(jì)算提供數(shù)據(jù)基礎(chǔ)。

3指標(biāo)選取、模型建立與洪災(zāi)風(fēng)險(xiǎn)區(qū)劃

3.1指標(biāo)選取與格網(wǎng)化

首先，從洪災(zāi)形成的四因子角度出發(fā)，選取洪災(zāi)的危險(xiǎn)性、脆弱性、暴露性及防災(zāi)減災(zāi)能力四個(gè)方面出發(fā)，進(jìn)行初步篩選 [14]；其次，在分析洪災(zāi)發(fā)生特點(diǎn)的基礎(chǔ)上選取了10個(gè)副因子及18個(gè)指標(biāo)數(shù)據(jù)；最后，利用指標(biāo)數(shù)據(jù)空間展布方法將指標(biāo)數(shù)據(jù)展布到網(wǎng)格中[15]。各指標(biāo)數(shù)據(jù)均采用歸一化方法進(jìn)行處理，以避免數(shù)據(jù)單位不統(tǒng)一帶來(lái)的計(jì)算不便，各因子權(quán)重利用層次分析法確定，指標(biāo)體系與權(quán)重計(jì)算結(jié)果見(jiàn)表1。

3.2洪災(zāi)風(fēng)險(xiǎn)模型建立

本研究中的洪災(zāi)風(fēng)險(xiǎn)模型構(gòu)建方法是采用四因子法與加權(quán)綜合法構(gòu)建，其中引入防災(zāi)減災(zāi)能力以體現(xiàn)洪災(zāi)風(fēng)險(xiǎn)綜合評(píng)估中的重要環(huán)節(jié)，這與以往研究中從致災(zāi)因子的危險(xiǎn)性、承災(zāi)體的暴露性與脆弱性三個(gè)方面構(gòu)建洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)模型方法有著本質(zhì)區(qū)別。以自然災(zāi)害風(fēng)險(xiǎn)數(shù)學(xué)表達(dá)式與洪災(zāi)發(fā)生特點(diǎn)為基礎(chǔ)，采用加權(quán)綜合法及層次分析法，綜合建立洪水災(zāi)害風(fēng)險(xiǎn)指數(shù)模型[1516]：

FFDR=（HWH）（EWE）（VWV）（0.1（1-a）R+a）[JY]（2）

H=WH1XH1+WH2XH2+WH3XH3+WH4XH4+WH5XH5[JY]（3）

E=WE1XE1+WE2XE2+WE3XE3+WE4XE4[JY]（4）

V=WV1XV1+WV2XV2+WV3XV3+WV4XV4[JY]（5）

R=WR1XR1+WR2XR2+WR3XR3+WR4XR4+WR5XR5[JY]（6）

式中：FFDR表示洪災(zāi)風(fēng)險(xiǎn)指數(shù)，其值大小表征風(fēng)險(xiǎn)大??；WHi、WEi、WVi、WRi分別表示洪災(zāi)風(fēng)險(xiǎn)的危險(xiǎn)性、暴露性、脆弱性及防災(zāi)減災(zāi)能力的權(quán)重值；H、E、V、R表示危險(xiǎn)性、暴露性、脆弱性及防災(zāi)減災(zāi)能力值的大小。XHi、XEi、XVi、XRi分別代表各指標(biāo)量化后的值。

3.3哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)區(qū)劃

[JP+2]利用建好的洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)模型，實(shí)現(xiàn)了研究區(qū)洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)，并將評(píng)價(jià)結(jié)果劃分為五個(gè)等級(jí)，即：低風(fēng)險(xiǎn)、較低風(fēng)險(xiǎn)、中等風(fēng)險(xiǎn)、較高風(fēng)險(xiǎn)及高等風(fēng)險(xiǎn)，最高值為5，最低值為1[15]。研究區(qū)利用ARCGIS軟件生成規(guī)則網(wǎng)格，利用疊置分析、關(guān)聯(lián)分析等方法關(guān)聯(lián)屬性數(shù)據(jù)，并在洪災(zāi)風(fēng)險(xiǎn)指數(shù)計(jì)算的基礎(chǔ)上，結(jié)合繪制的網(wǎng)格及專(zhuān)題地圖制作功能，綜合繪制哈爾濱市區(qū)2000年-2009年的洪災(zāi)風(fēng)險(xiǎn)指數(shù)及風(fēng)險(xiǎn)等級(jí)圖（圖1），以此分析哈爾濱市洪水災(zāi)害風(fēng)險(xiǎn)空間演變規(guī)律，為城市建設(shè)和可持續(xù)發(fā)展提供決策依據(jù)。

4哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)空間演變模式分析

以上述研究結(jié)果為基礎(chǔ)，利用ARCGIS中的空間分析功能，分析哈爾濱市2000年-2009年10年間的空間變化過(guò)程，包括哈爾濱市洪災(zāi)風(fēng)險(xiǎn)變化在任意兩年間的風(fēng)險(xiǎn)變化模式，并對(duì)分析結(jié)果進(jìn)行了總結(jié)，得出哈爾濱市洪災(zāi)風(fēng)險(xiǎn)在空間上的變化特征，并總結(jié)出洪災(zāi)風(fēng)險(xiǎn)空間演變模式分別為升高模式與降低模式，具體見(jiàn)表2。

風(fēng)險(xiǎn)升高模式中的邊緣接觸式較為明顯的為2008年-2009年，出現(xiàn)在道外區(qū)中部，主要原因?yàn)樵谠袨?zāi)環(huán)境不變的情況下，承災(zāi)體增多導(dǎo)致洪災(zāi)風(fēng)險(xiǎn)將越來(lái)越大；內(nèi)部突變演變模式出現(xiàn)在2000年-2001年，主要原因?yàn)榈劳鈪^(qū)在偏南部的孕災(zāi)環(huán)境與承災(zāi)體都發(fā)生了變化，另外由于下墊面的道路變?yōu)闉r青道路，一些綠地變?yōu)閺V場(chǎng)，導(dǎo)致下墊面滲透率減少、糙率變小及下墊面產(chǎn)匯流時(shí)間變短，加之承災(zāi)體數(shù)量增多，最終導(dǎo)致洪災(zāi)風(fēng)險(xiǎn)呈現(xiàn)內(nèi)部突變式增長(zhǎng)；而擺動(dòng)式與跨越式洪災(zāi)風(fēng)險(xiǎn)增長(zhǎng)模式主要發(fā)生在道里區(qū)和道外區(qū)，其原因是由于道里區(qū)和道外區(qū)承災(zāi)體增多，導(dǎo)致洪災(zāi)風(fēng)險(xiǎn)增大，但在下一年由于防災(zāi)減災(zāi)能力增加，反而導(dǎo)致洪災(zāi)風(fēng)險(xiǎn)變小。在洪災(zāi)風(fēng)險(xiǎn)降低模式[JP+1]原因，發(fā)現(xiàn)三種模式中導(dǎo)致洪災(zāi)風(fēng)險(xiǎn)降低的原因相似，都是發(fā)生在道里區(qū)和道外區(qū)，因道里區(qū)和道外區(qū)是哈爾濱老城區(qū)，其孕災(zāi)環(huán)境相對(duì)變化不大，承災(zāi)體變化也相對(duì)穩(wěn)定，但政府部門(mén)防災(zāi)減災(zāi)能力逐年變強(qiáng)，導(dǎo)致道里區(qū)和道外區(qū)洪災(zāi)風(fēng)險(xiǎn)變小。上述洪災(zāi)風(fēng)險(xiǎn)演變模式在不同年份間因各地的孕災(zāi)環(huán)境、承災(zāi)體及防災(zāi)減災(zāi)能力變化不同，造成了道里區(qū)和道外區(qū)洪災(zāi)風(fēng)險(xiǎn)演變模式不同。

綜上所述，洪災(zāi)風(fēng)險(xiǎn)大小變化能夠體現(xiàn)不同年份間研究區(qū)孕災(zāi)環(huán)境、承災(zāi)體及防災(zāi)減災(zāi)能力的變化情況，該研究結(jié)果可為應(yīng)急管理部門(mén)在不同年份間采取不同應(yīng)急措施提供決策依據(jù)，同時(shí)也為掌握洪災(zāi)風(fēng)險(xiǎn)轉(zhuǎn)移規(guī)律提供參考。

5結(jié)語(yǔ)

以5 km×5 km大小的網(wǎng)格作為洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)單元，采用自然災(zāi)害指數(shù)法、層次分析法并結(jié)合GIS技術(shù)，實(shí)現(xiàn)了哈爾濱市區(qū)洪災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)，并繪制了網(wǎng)格尺度下的哈爾濱市區(qū)2000年-2009年洪災(zāi)風(fēng)險(xiǎn)區(qū)劃圖，實(shí)現(xiàn)了哈爾濱市區(qū)洪水災(zāi)害風(fēng)險(xiǎn)演變模式歸類(lèi)。歸類(lèi)結(jié)果分為兩種演變模式，其中升高模式包括邊緣接觸模式、內(nèi)部突增模式、擺動(dòng)模式和跨越模式；降低模式包括邊緣退縮模式、內(nèi)部突降模式和片狀降低模式。由于研究中搜集的數(shù)據(jù)資料有限，且洪災(zāi)風(fēng)險(xiǎn)空間演變受人文、社會(huì)及經(jīng)濟(jì)等因素綜合影響較大，因此，在進(jìn)行指標(biāo)數(shù)據(jù)空間網(wǎng)格展布時(shí)存在誤差，需在今后的工作中尋求提高指標(biāo)網(wǎng)格化精度的方法，更加精確的表達(dá)洪災(zāi)空間演變模式，即提高時(shí)間和空間上各數(shù)據(jù)分辨率，使洪災(zāi)風(fēng)險(xiǎn)演變模式歸類(lèi)更為精確。研究結(jié)果可為城市防洪規(guī)劃、應(yīng)急等方面提供決策依據(jù)。

參考文獻(xiàn)（References）：

[1]黃崇福，史培軍.城市地震災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)的數(shù)學(xué)模型[J].自然災(zāi)害學(xué)報(bào)，1995，4（2）：3037.（HUANG Chongfu，SHI Peijun.Mathematical model risk assessment of earthquake disaster in city[J].Journal of Natural Disasters，1995，4（2）：3037.（ in Chinese））

[2]黃詩(shī)峰.洪水災(zāi)害風(fēng)險(xiǎn)分析的理論與方法研究[D].北京：中國(guó)科學(xué)院地理研究所，1999.（HUANG Shifeng.The theory and method of flood disaster risk analysis[D].Beijing：Institute of Geography，Chinese Academy of Sciences，1999.（ in Chinese））

[3]周成虎，萬(wàn)慶，黃詩(shī)峰.基于GIS的洪水災(zāi)害風(fēng)險(xiǎn)區(qū)劃研究[J].地理科學(xué)，2000，55（1）：1524.（ZHOU Chenghu，WAN Qing，HUANG Shifeng.A GISbased approach to flood risk zona tion[J].Scientia Geographica Sinica，2000，55（1）：1524.（ in Chinese））