基于GIS的水庫洪水風(fēng)險(xiǎn)圖編制

王 軍,梁忠民,施 曄

(河海大學(xué)水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098)

基于GIS的水庫洪水風(fēng)險(xiǎn)圖編制

王 軍,梁忠民,施 曄

(河海大學(xué)水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098)

利用水文學(xué)方法、水庫潰壩洪水計(jì)算模型及二維水力學(xué)數(shù)學(xué)模型分別對水庫庫區(qū)洪水風(fēng)險(xiǎn)和水庫潰壩洪水風(fēng)險(xiǎn)進(jìn)行分析,利用Microsoft Visual Basic 6.0對ArcGIS進(jìn)行二次開發(fā),將風(fēng)險(xiǎn)信息按照一定的方式自動(dòng)顯示在地圖上,從而實(shí)現(xiàn)了水庫洪水風(fēng)險(xiǎn)圖的自動(dòng)編制,并將該方法應(yīng)用于江蘇省石梁河水庫洪水風(fēng)險(xiǎn)圖的編制.應(yīng)用結(jié)果表明,通過該方法編制的水庫洪水風(fēng)險(xiǎn)圖可以為當(dāng)?shù)胤姥床块T制訂科學(xué)合理的防洪調(diào)度方案和抗洪救災(zāi)措施提供科學(xué)依據(jù)和技術(shù)支撐.

水庫;潰壩洪水;GIS;洪水風(fēng)險(xiǎn)圖

洪水風(fēng)險(xiǎn)圖是一種融合地理、社會(huì)經(jīng)濟(jì)信息、洪水特征信息,通過資料調(diào)查、洪水計(jì)算和成果整理,以地圖形式直觀反映某一地區(qū)發(fā)生某一頻率的洪水后可能淹沒的范圍和水深,從而可以預(yù)知和分析不同量級洪水可能造成的風(fēng)險(xiǎn)和災(zāi)害的專題地圖,它可以為相關(guān)部門的防洪減災(zāi)決策提供指導(dǎo)和決策依據(jù)[1-2].洪水風(fēng)險(xiǎn)圖的編制在發(fā)達(dá)國家已發(fā)展較為成熟[3-4],而在我國的發(fā)展卻剛剛起步.隨著我國治水理念的改變,編制洪水風(fēng)險(xiǎn)圖將成為我國防洪減災(zāi)迫切需要開展的重點(diǎn)工作之一[5].洪水風(fēng)險(xiǎn)圖按照編制對象的不同可分為江河湖泊洪水風(fēng)險(xiǎn)圖、蓄滯洪區(qū)洪水風(fēng)險(xiǎn)圖、水庫洪水風(fēng)險(xiǎn)圖[1].其中,水庫洪水風(fēng)險(xiǎn)圖是反映水庫庫區(qū)及水庫潰壩情況下洪水風(fēng)險(xiǎn)信息的專題地圖.水庫建成后,存在的洪水風(fēng)險(xiǎn)主要有2種:(a)庫區(qū)的淹沒風(fēng)險(xiǎn).因?yàn)樗畮焖灰话憔S持在防洪高水位以下,在防洪高水位以上的庫區(qū)范圍會(huì)有農(nóng)作物等的分布,當(dāng)發(fā)生較大洪水時(shí),農(nóng)作物等將有可能被淹.(b)水庫潰壩洪水對下游區(qū)域造成的淹沒風(fēng)險(xiǎn).因此,水庫洪水風(fēng)險(xiǎn)圖可以分為水庫庫區(qū)洪水風(fēng)險(xiǎn)圖和潰壩洪水風(fēng)險(xiǎn)圖.

地理信息系統(tǒng)(GIS)是利用計(jì)算機(jī)強(qiáng)大的計(jì)算能力實(shí)現(xiàn)自動(dòng)或半自動(dòng)存儲(chǔ)、處理、分析和管理空間數(shù)據(jù)的一種技術(shù)[6-8].它融合了數(shù)據(jù)庫、軟件工程、人工智能、網(wǎng)絡(luò)技術(shù)等[6]主流技術(shù),憑借其強(qiáng)大的空間處理、圖形展現(xiàn)、地圖處理等眾多優(yōu)點(diǎn),已經(jīng)成為很多重要領(lǐng)域的有力工具,并且在洪水風(fēng)險(xiǎn)的自動(dòng)編制上,GIS也大有作為,利用GIS技術(shù)可高效、準(zhǔn)確地半自動(dòng)化自動(dòng)化編制洪水風(fēng)險(xiǎn)圖.同時(shí)基于GIS進(jìn)行二次開發(fā),還可以開發(fā)出符合用戶特殊需求的洪水風(fēng)險(xiǎn)圖管理系統(tǒng).

我國水庫眾多[9-11],現(xiàn)有大壩總數(shù)超過8.5萬座,居世界首位.由于我國大壩大都修建于20世紀(jì)50～70年代,設(shè)計(jì)標(biāo)準(zhǔn)低,工程質(zhì)量差[12],而且很多大壩年久失修,帶病運(yùn)行,每逢遭遇大水,必定存在潰壩隱患,并可能導(dǎo)致庫區(qū)及水庫下游地區(qū)重大經(jīng)濟(jì)損失和生命財(cái)產(chǎn)損失.因此,水庫洪水風(fēng)險(xiǎn)圖的編制對防汛部門了解水庫發(fā)生重大洪水災(zāi)害時(shí)可能產(chǎn)生的后果,并制訂和采取有效的避險(xiǎn)措施具有重要的意義.

隨著GIS技術(shù)的興起與成熟,基于GIS的水庫洪水風(fēng)險(xiǎn)圖的編制漸漸成為人們研究的熱點(diǎn).本文利用水文學(xué)方法、水庫潰壩洪水計(jì)算模型及二維水力學(xué)數(shù)學(xué)模型分別對水庫庫區(qū)洪水風(fēng)險(xiǎn)和水庫潰壩洪水風(fēng)險(xiǎn)進(jìn)行分析,利用Microsoft Visual Basic 6.0對ArcGIS進(jìn)行二次開發(fā),將分析所得風(fēng)險(xiǎn)信息按照一定的方式自動(dòng)顯示在地圖上,從而實(shí)現(xiàn)了水庫洪水風(fēng)險(xiǎn)圖的自動(dòng)編制.

1 水庫洪水風(fēng)險(xiǎn)分析方法

洪水風(fēng)險(xiǎn)是指未來可能引起災(zāi)害性后果的洪水發(fā)生的概率或頻率.洪水風(fēng)險(xiǎn)分析的主要內(nèi)容為:(a)引起災(zāi)害性后果的不同頻率洪水的特征值或洪水過程;(b)不同頻率洪水引起的災(zāi)害性后果,如淹沒范圍、淹沒水深、淹沒歷時(shí)、洪災(zāi)損失等.根據(jù)水庫洪水風(fēng)險(xiǎn)可能發(fā)生的實(shí)際情況,可將水庫洪水風(fēng)險(xiǎn)分為庫區(qū)洪水風(fēng)險(xiǎn)和水庫潰壩洪水風(fēng)險(xiǎn).

1.1 庫區(qū)洪水風(fēng)險(xiǎn)

庫區(qū)洪水風(fēng)險(xiǎn)信息主要是指不同頻率洪水對應(yīng)的淹沒范圍及淹沒水深,庫區(qū)洪水風(fēng)險(xiǎn)分析的步驟如下:(a)通過水庫調(diào)洪演算,確定不同頻率洪水的庫區(qū)水位過程及淹沒水位;利用庫區(qū)水量平衡方程和庫容曲線,聯(lián)合求解不同頻率洪水庫區(qū)水位過程.(b)由水位數(shù)據(jù)和地面高程數(shù)據(jù)推求庫區(qū)不同位置的淹沒水深.(c)結(jié)合數(shù)字地面高程模型確定不同頻率洪水時(shí)水庫庫區(qū)淹沒范圍及淹沒水深分布.

1.2 水庫潰壩洪水風(fēng)險(xiǎn)

水庫潰壩洪水風(fēng)險(xiǎn)信息包括潰壩條件下水庫下游區(qū)域的洪水最大淹沒范圍、最大淹沒水深分布、淹沒歷時(shí)分布、流速分布,以及不同時(shí)刻洪水演進(jìn)過程的不同水文要素信息等.水庫一旦潰壩,潰壩洪水會(huì)沖出河道,四處奔襲,這種下墊面地物復(fù)雜、水流路徑不定的情況需要用二維水力學(xué)數(shù)學(xué)模型進(jìn)行模擬和分析.而二維水力學(xué)數(shù)學(xué)模型所模擬的潰壩處的入流過程就是潰口發(fā)展相應(yīng)的洪水過程,因此要運(yùn)用二維水力學(xué)數(shù)學(xué)模型,首先要分析潰壩潰口洪水過程,才能推求出潰壩洪水引起的下游淹沒信息.實(shí)際上,潰壩有瞬間潰壩和漸進(jìn)潰壩2種形式,這里僅研究漸進(jìn)潰壩的洪水風(fēng)險(xiǎn)問題.

1.2.1 潰壩洪水過程計(jì)算

1.2.1.1 大壩潰口寬度確定

如有相關(guān)調(diào)查資料和統(tǒng)計(jì)資料,則根據(jù)資料確定大壩最終潰口寬度;如缺乏資料,根據(jù)《洪水風(fēng)險(xiǎn)圖編制導(dǎo)則(試行)》[1]中的規(guī)定,潰口形狀按梯形斷面考慮,并參考陸吉康等[13]提出的經(jīng)驗(yàn)公式進(jìn)行計(jì)算.

最終潰口平均寬度計(jì)算的經(jīng)驗(yàn)公式為

式中:Hb——潰決有效高度(水庫大壩潰決時(shí)刻水位減壩址斷面平均底高程);Vr——水庫有效下泄庫容;Bm——最終潰口的平均寬度,即梯形潰口上下寬度的平均值;K——修正系數(shù),漫頂造成的潰決K=1,管涌造成的潰決K=1.4.

t時(shí)潰口平均寬度Bt的計(jì)算公式為

式中:B0——初始寬度;Tf——潰口發(fā)展歷時(shí);t——潰口時(shí)間.

1.2.1.2 大壩潰口底高程確定

在潰口發(fā)展歷時(shí)Tf內(nèi),大壩潰口底高程是不斷變化的,假定在潰口發(fā)展歷時(shí)Tf內(nèi),大壩潰口縱向按線性比例擴(kuò)大,則潰口底高程

式中:hb——潰口底高程;hd——壩頂高程;hbm——最終潰口底高程;ρ0——反映潰口非線性程度的參數(shù),取值范圍為1≤ρ0≤4.

1.2.1.3 潰壩泄流流量計(jì)算

假定潰壩斷面水流流態(tài)與寬頂堰流流態(tài)類似,則潰壩泄流流量可按堰流公式進(jìn)行計(jì)算:

式中:m——流量系數(shù),近似取0.385;Ht——t時(shí)潰口底部高程以上水頭,即t時(shí)庫水位與hb的差值.

1.2.2 潰壩洪水沿程演進(jìn)計(jì)算

在求出水庫潰壩潰口洪水過程的基礎(chǔ)上,采用二維水力學(xué)數(shù)學(xué)模型進(jìn)行水庫潰壩洪水的沿程演進(jìn)計(jì)算,以推求潰壩洪水在下游區(qū)域形成的淹沒水深、淹沒范圍、流速流量分布、淹沒歷時(shí)等風(fēng)險(xiǎn)信息.

二維水力學(xué)數(shù)學(xué)模型由連續(xù)方程和水流沿x和y方向的動(dòng)量方程組成.連續(xù)方程為

二維水力學(xué)數(shù)學(xué)模型的求解步驟如圖1所示.

圖1 二維水力學(xué)數(shù)學(xué)模型的求解步驟Fig.1 Generic steps for solving two-dimensional hydraulic model

2 GIS支持下水庫洪水風(fēng)險(xiǎn)圖編制方法

由于水庫風(fēng)險(xiǎn)分為庫區(qū)洪水風(fēng)險(xiǎn)和水庫潰壩洪水風(fēng)險(xiǎn),因此分析時(shí)二者要分別進(jìn)行考慮.如果要使庫區(qū)洪水風(fēng)險(xiǎn)圖能夠反映洪水風(fēng)險(xiǎn)信息,則需要繪制不同頻率設(shè)計(jì)洪水的庫區(qū)淹沒范圍圖和淹沒水深分布圖.由于洪水的到達(dá)時(shí)間、淹沒歷時(shí)、淹沒水深變化過程等信息是防汛部門組織民眾避險(xiǎn)的重要依據(jù),所以水庫潰壩洪水風(fēng)險(xiǎn)圖需包括最大淹沒范圍圖、最大水深分布圖、最大流速分布圖、到達(dá)時(shí)間分布圖、淹沒歷時(shí)分布圖以及潰壩洪水演進(jìn)過程中不同時(shí)刻的淹沒范圍圖、淹沒水深分布圖、流速分布圖、流向分布圖等.以上各種洪水風(fēng)險(xiǎn)圖需采用一定的顏色標(biāo)注才可以清晰地反映洪水風(fēng)險(xiǎn)信息.根據(jù)《洪水風(fēng)險(xiǎn)圖編制導(dǎo)則(試行)》[1]的規(guī)定 ,采用灰(RGB:128 128 128)、藍(lán)(RGB:0 0 255)、綠(RGB:0 255 0)、黃(RGB:255 255 0)、橙(RGB:255 165 0)、紅(RGB:255 0 0)6種顏色分別標(biāo)注5年一遇、10年一遇、20年一遇、50年一遇、100年一遇、1000年一遇洪水風(fēng)險(xiǎn)區(qū)域及最大淹沒范圍.其他洪水風(fēng)險(xiǎn)圖用相應(yīng)淹沒范圍的顏色分級標(biāo)注.同時(shí),各洪水風(fēng)險(xiǎn)圖層只有與防洪工程分布圖、水系圖層、地形圖層、行政區(qū)劃圖層等疊加顯示,才能為防汛管理提供依據(jù).

各種GIS軟件都具有圖層渲染功能,為了實(shí)現(xiàn)洪水風(fēng)險(xiǎn)圖的自動(dòng)編制,筆者采用世界GIS領(lǐng)袖ESRI公司的ArcGIS桌面軟件處理原始空間數(shù)據(jù),利用ArcObject(AO)提供的高層開發(fā)控件和組件以及可視化開發(fā)語言VB進(jìn)行應(yīng)用程序開發(fā).洪水風(fēng)險(xiǎn)圖編制的技術(shù)路線為:

a.利用AO提供的接口(Ilayer)新建一個(gè)用來存儲(chǔ)各頻率設(shè)計(jì)洪水或潰壩洪水淹沒信息的洪水風(fēng)險(xiǎn)圖層.

b.對新建圖層進(jìn)行網(wǎng)格剖分.需要指出的是,剖分網(wǎng)格時(shí)要充分考慮研究區(qū)域內(nèi)的阻水建筑物、擋水建筑物、過水建筑物等復(fù)雜邊界條件.計(jì)算網(wǎng)格中的地形信息即高程值是計(jì)算的基礎(chǔ),因此剖分網(wǎng)格時(shí)要融合對應(yīng)網(wǎng)格的高程信息.可以根據(jù)資料實(shí)際擁有情況,通過對DEM、高程點(diǎn)、等高線的處理來獲得計(jì)算網(wǎng)格的高程信息.完整洪水風(fēng)險(xiǎn)原始圖層,可以是矢量格式的,也可以是柵格格式的.

c.將洪水風(fēng)險(xiǎn)分析計(jì)算得到的各個(gè)網(wǎng)格的洪水風(fēng)險(xiǎn)信息,包括最大水深、最大流速、到達(dá)時(shí)間、淹沒歷時(shí)、洪水演進(jìn)過程中不同時(shí)刻的淹沒水深、流速等,分別作為洪水風(fēng)險(xiǎn)圖層的一個(gè)屬性信息寫入該圖層屬性表的相應(yīng)字段中.

d.針對不同的洪水風(fēng)險(xiǎn)圖,利用AO的圖層描述接口(IRenderer)、顏色方案接口(IcolorRamp),并按照設(shè)定的色彩方案,以屬性表中相應(yīng)字段的取值進(jìn)行分類即可對該圖層進(jìn)行渲染,從而實(shí)現(xiàn)洪水風(fēng)險(xiǎn)圖層的自動(dòng)編制.

e.將渲染完畢的洪水風(fēng)險(xiǎn)圖層與其他基礎(chǔ)圖層進(jìn)行疊加,即可得到相應(yīng)洪水要素的洪水風(fēng)險(xiǎn)圖.洪水風(fēng)險(xiǎn)圖可以通過開發(fā)的應(yīng)用程序查看,也可以通過ArcGIS軟件查看.

3 水庫洪水風(fēng)險(xiǎn)圖編制實(shí)例

江蘇省石梁河水庫位于新沭河中游,蘇魯兩省的贛榆、東海、臨沭縣三縣交界處,是一座具有綜合功能的大(2)型均質(zhì)土壩水庫.水庫按100年一遇設(shè)計(jì),2000年一遇校核.石梁河水庫一旦潰壩,對下游的贛榆、東海兩縣產(chǎn)生極大的威脅,并可能使兩縣遭受巨大的社會(huì)經(jīng)濟(jì)、生命財(cái)產(chǎn)、生態(tài)環(huán)境等損失,同時(shí)可能導(dǎo)致下游新沭河堤防潰決并威脅下游的連云港市.因此,石梁河水庫洪水風(fēng)險(xiǎn)圖對相關(guān)部門指導(dǎo)當(dāng)?shù)厝嗣穹罏?zāi)減災(zāi)具有重要意義.

3.1 庫區(qū)洪水風(fēng)險(xiǎn)圖

采用上述庫區(qū)洪水風(fēng)險(xiǎn)分析方法,結(jié)合石梁河水庫實(shí)際水文數(shù)據(jù),分析得到的不同頻率洪水庫區(qū)水位如下:50年一遇洪水庫區(qū)最高水位為25.99m;100年一遇洪水庫區(qū)最高水位為26.72m;1000年一遇洪水庫區(qū)最高水位為27.53m;2000年一遇洪水庫區(qū)最高水位為27.86m.

結(jié)合石梁河水庫洪水頻率分析結(jié)果,用不同頻率最高水位減去庫區(qū)DEM,即可得到庫區(qū)洪水風(fēng)險(xiǎn)信息,按重現(xiàn)期從小到大分別采用不同的色彩方案對相應(yīng)頻率洪水風(fēng)險(xiǎn)信息進(jìn)行渲染,從而得到不同頻率洪水庫區(qū)洪水風(fēng)險(xiǎn)圖.庫區(qū)不同頻率疊加淹沒范圍和2000年一遇洪水淹沒水深分布分別如圖2和圖3所示.

圖2 庫區(qū)不同頻率疊加淹沒范圍Fig.2 Superposed submerged area in reservoir region with floods of different frequencies

圖3 庫區(qū)2000年一遇洪水淹沒水深分布Fig.3 Distribution of water depth in reservoir region under flood with return period of 2000 years

3.2 潰壩洪水風(fēng)險(xiǎn)圖

根據(jù)上述二維水力學(xué)數(shù)學(xué)模型的求解步驟,結(jié)合石梁河水庫實(shí)際情況,劃定研究區(qū)域(圖4),并剖分網(wǎng)格.網(wǎng)格剖分采用研究區(qū)域內(nèi)的新沭河堤防、海堤、汾灌高速公路等確定計(jì)算邊界,并從地圖上提取該邊界.網(wǎng)格剖分采用專業(yè)軟件GAMBIT,根據(jù)地形數(shù)據(jù)的實(shí)際情況,將網(wǎng)格大小設(shè)計(jì)為500m×500m.石梁河水庫有主壩和老副壩,以老副壩為例,其潰壩洪水計(jì)算網(wǎng)格如圖4所示.

用式(5)～(7)推求潰壩洪水過程,并將此洪水過程作為二維水力學(xué)數(shù)學(xué)模型的入流條件,根據(jù)下墊面地物分布情況設(shè)定糙率,運(yùn)行模型,可以求出每一個(gè)計(jì)算網(wǎng)格的淹沒要素值,進(jìn)而可根據(jù)潰壩洪水風(fēng)險(xiǎn)分析結(jié)果,繪制最大水深分布圖、最大淹沒范圍圖、最大流速分布圖、到達(dá)時(shí)間分布圖以及洪水演進(jìn)過程任意時(shí)刻水文要素的分布圖.老副壩潰壩洪水最大水深分布和到達(dá)時(shí)間分別如圖5和圖6所示.

圖4 研究區(qū)域及老副壩潰壩計(jì)算網(wǎng)格Fig.4 General situation of study area and calculation grids of dam-break

圖5 老副壩潰壩洪水最大水深分布Fig.5 Distribution of maximumwater depth of dam-break flood

圖6 老副壩潰壩洪水到達(dá)時(shí)間Fig.6 Distribution of arrival time of dam-break flood

從圖5和圖6可以看出洪水淹沒水深分布和到達(dá)時(shí)間分布.由于研究區(qū)域歷史上沒有發(fā)生過潰壩事件,不可能通過實(shí)測數(shù)據(jù)來定量地驗(yàn)證模型計(jì)算結(jié)果,故定性分析如下:(a)如圖5所示,洪水最大水深分布符合研究區(qū)域內(nèi)水庫附近為丘陵、下游為平原、周邊是堤防、從西北向東南方向高程遞減的特點(diǎn),因此該計(jì)算結(jié)果是合理的.(b)如圖6所示,由于汾灌高速公路的阻水作用,汾灌高速公路左側(cè)和右側(cè)的洪水到達(dá)時(shí)間有較為明顯的差別,該公路右側(cè)附近有一片顏色較深的區(qū)域,從地形信息可知,該處高程較高,所以洪水到達(dá)時(shí)間較長.因此,該計(jì)算結(jié)果也是合理的.

4 結(jié) 語

本研究基于GIS組件式技術(shù),利用Microsoft Visual Basic 6.0對ArcGIS進(jìn)行二次開發(fā),并將水文學(xué)方法與水力學(xué)洪水分析方法相耦合,從而實(shí)現(xiàn)了水庫洪水風(fēng)險(xiǎn)圖的自動(dòng)編制.應(yīng)用結(jié)果表明,該水庫洪水風(fēng)險(xiǎn)圖編制方法合理可行,采用此方法編制的石梁河水庫洪水風(fēng)險(xiǎn)圖可以為連云港市制訂科學(xué)合理的防洪調(diào)度方案和抗洪救災(zāi)措施提供科學(xué)依據(jù)和技術(shù)支撐.但是,從石梁河水庫洪水風(fēng)險(xiǎn)圖編制的研究中也發(fā)現(xiàn)了如下一些問題:(a)實(shí)際上,水庫庫容包括靜庫容和動(dòng)庫容,而動(dòng)庫容的確定歷來都是庫容計(jì)算的難點(diǎn),因此水庫尤其是大型水庫庫區(qū)洪水風(fēng)險(xiǎn)圖的編制需考慮動(dòng)庫容并對動(dòng)庫容進(jìn)行分析,才能使洪水風(fēng)險(xiǎn)圖信息更為準(zhǔn)確.(b)水庫潰壩有瞬間潰和漸進(jìn)潰2種形式,這2種潰壩形式在物理機(jī)制上有很大差異,形成的泄流過程截然不同,因而洪水演進(jìn)過程會(huì)產(chǎn)生很大差異.由于正確的洪水分析方法是洪水風(fēng)險(xiǎn)圖信息準(zhǔn)確的基礎(chǔ),是保證洪水風(fēng)險(xiǎn)圖具有科學(xué)性、合理性、實(shí)用性的基礎(chǔ),因此這2種潰壩形式的發(fā)生條件及泄流公式值得進(jìn)一步研究.

[1]國家防汛抗旱總指揮部辦公室.洪水風(fēng)險(xiǎn)圖編制導(dǎo)則(試行)[K].北京:國家防汛抗旱總指揮部辦公室,2005.

[2]曹東,金東春.洪水風(fēng)險(xiǎn)圖及其作用[J].東北水利水電,1998(8):8-10.(CAO Dong,JIN Dong-chun.Flood risk map and its functions[J].Water Resourse&Hydropower of Northeast China,1998(8):8-10.(in Chinese))

[3]馬建明,許靜,朱云楓,等.國外洪水風(fēng)險(xiǎn)圖編制綜述[J].中國水利,2005(17):29-31.(MA Jian-ming,XU Jing,ZHU Yun-feng,et al.General introduction of flood hazard mapping in other countries[J].China Water Resources,2005(17):29-31.(in Chinese))

[4]李娜,向立云,程曉陶.國外洪水風(fēng)險(xiǎn)圖制作比較及對我國洪水風(fēng)險(xiǎn)圖制作的建議[J].水利發(fā)展研究,2005(6):28-32.(LI Na,XIANG Li-yun,CHENG Xiao-tao.Comparsion of foreign drawing flood riskmap and suggestionsfor how to draw flood risk map in our country[J].Water Resources Development Research,2005(6):28-32.(in Chinese))

[5]張志彤.落實(shí)“兩個(gè)轉(zhuǎn)變” 大力推進(jìn)洪水風(fēng)險(xiǎn)圖編制工作[J].中國水利,2005(17):5-6.(ZHANG Zhi-tong.Implementing“two transfers”and making great efforts on flood hazard mapping work[J].China Water Resources,2005(17):5-6.(in Chinese))

[6]李娜.基于GIS的洪災(zāi)風(fēng)險(xiǎn)管理系統(tǒng)[D].北京:中國水利水電科學(xué)研究院,2002.

[7]王靜,王軍,梁忠民.基于GIS的連云港市區(qū)城市洪水風(fēng)險(xiǎn)圖管理系統(tǒng)研究[C]//高丹盈,左其亭.中國水論壇第四屆學(xué)術(shù)研討會(huì)論文集.北京:中國水利水電出版社,2006:1112-1116.

[8]黨安榮,賈海峰,易善楨,等.地理信息系統(tǒng)應(yīng)用指南[M].北京:清華大學(xué)出版社,2002.

[9]楊旭.水庫洪水風(fēng)險(xiǎn)圖的制作[J].江西水利科技,1999(增刊1):82-87.(YANG Xu.Work out the flood risk figures of reservoirs[J].Jiangxi Hydraulic Sciences&Technology,1999(S1):82-87.(in Chinese))

[10]唐永祥,夏雨豐.二維潰壩洪水計(jì)算方法在洪水風(fēng)險(xiǎn)圖編制中的應(yīng)用[J].浙江水利科技,2004(3):32-33.(TANG Yongxiang,XIA Yu-feng.Application of two-dimensional calculation method for dam breach flood to making flood risk diagram[J].Zhejiang Hydrotechnics,2004(3):32-33.(in Chinese))

[11]姜樹海,范子武,吳時(shí)強(qiáng).洪災(zāi)風(fēng)險(xiǎn)評估和防洪安全決策[M].北京:中國水利水電出版社,2005.

[12]李雷,王仁鐘,盛金保,等.大壩風(fēng)險(xiǎn)評價(jià)與風(fēng)險(xiǎn)管理[M].北京:中國水利水電出版社,2006.

[13]陸吉康,林秉南,龔振氵云,等.棱柱體河道的潰壩波[J].江西水利科技,1985(1):57-60.(LU Ji-kang,LIN Bing-nan,GONG Zhen-yun,et al.Dam-break flood spread of prism channel[J].Jiangxi Hydraulic Sciences&Technology,1985(1):57-60.(in Chinese))

Mapping of flood risk of reservoirs using GIS technology

WANG Jun,LIANG Zhong-min,SHI Ye
(StateKey Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University,Nanjing210098,China)

Hydrological methods,dam-break flood models and two-dimensional hydraulic models were employed to analyze the risk of flooding in reservoir regions and flooding caused by dam breakage.The secondary development of ArcGIS was performed using Microsoft Visual Basic 6.0,and risk information was automatically displayed on a map in a certain way so as to realize the automatic mapping of flooding risk of reservoirs.The proposed method was applied to the Shilianghe Reservoir in Jiangsu Province.The results show that the flood risk map obtained with the proposedmethod can provide a scientific basis and technical support for local departments to compile schemes for flood prevention and operation and to take measures for flood control and disaster mitigation.

reservoir;dam-break flood;GIS;flood risk map

TV122.4

A

1000-1980(2010)01-0020-06

10.3876/j.issn.1000-1980.2010.01.005

2008-11-10

江蘇省創(chuàng)新基金(CX08B_105Z);國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)(2007CB714104);國家自然科學(xué)基金(50779013);教育部博士點(diǎn)基金(20070294018);高等學(xué)校學(xué)科創(chuàng)新引智計(jì)劃(B08048)

王軍(1982—),男,安徽亳州人,博士研究生,主要從事防災(zāi)減災(zāi)及水文不確定性分析研究.E-mail:wfmingyu@hhu.edu.cn