長江黃石段堤防決口經濟損失估算

【摘要】我國是洪水災害嚴重的國家，近幾年隨著社會經濟的快速發展，流域發生大洪水導致堤防決口情況下產生的洪災損失也越來越嚴重。依托全國重點地區洪水風險圖編制工作，以黃石地區為研究對象，對區域進行分析計算，采用MIKE FLOOD洪水演進一二維耦合數學模型模擬河道潰堤洪水在區域內的演進情況，分析洪水風險，在洪水風險信息計算的基礎上，對各量級洪水導致的居民財產、農林牧漁、工商企業、交通運輸、水利設施等方面的直接損失估算分析，為黃石地區洪水風險管理及防洪調度決策提供技術參考。

【關鍵詞】數學模型，洪水風險，損失評估

1、引言

防洪保護區洪水風險圖是標示區內各處受洪水災害危險程度的一種防洪非工程措施，是抵御洪澇災害、保證防洪安全的重要參考。通過對可能發生洪水堤防潰決后，洪水演進路線、到達時間、淹沒水深、淹沒范圍、流速大小及避洪轉移范圍、路線、安置場所等要素進行展示，為各級防汛指揮部門的抗洪搶險救災提供決策依據，指導洪水威脅地區的建設，在發生洪水的緊急情況下選擇正確路線和地點疏散群眾，并為洪澇災情評估提供基礎依據。

隨著人們水患意識的加強，減災工作受到重視，不斷有新方法和新技術應用到評估分析中。基于GIS的洪災損失評估方法則是運用GIS技術研究洪水災害，該方法不僅能準確快速地提供研究所需要的數據，而且能對數據進行進一步的處理。

2、研究區域簡介

黃石市堤防洪保護區處于長江中游南岸湖北省黃石市境內，屬亞熱帶大陸型氣候區。區域年降水量豐沛，歷年平均降雨量為1413mm，最大值為2258mm（1954年）。年降水時空分布不均，多集中在每年4至8月份（占全年平均雨量的65%），雨季時間長，且暴雨頻繁。保護區內有黃石市及鄂州市，還有武九鐵路、山南鐵路和大廣高速、滬渝高速、G45等重要的交通設施，一旦發生洪澇災害，損失較大，歷來是湖北省防御洪澇災害的重點地區。

長江位于保護區東面，汛期長江黃石港水位一般為22.5～23.5m（1954年為24.54m），而黃石市區地面高程為16.5～18.1m，低于長江水位5～6m；因此，每年汛期，整個黃石無論是市區還是郊區均處于長江干流洪水威脅之下，因此有必要計算長江干流發生大洪水黃石市堤堤防潰決情況下，防洪保護區內的洪水演進和潰堤后的經濟損失情況，為黃石地區洪水風險管理及防洪調度決策提供技術參考。

3、洪水風險計算

3.1 洪水風險計算模型

洪水風險分析方法主要有水文學方法、實際水災法以及水力學法，水力學法是較為經典的計算方法，在洪水演進模擬方面得到了廣泛的運用，是較為成熟的分析方法。長江干堤黃石段防洪保護區堤防決口洪水風險計算采用為水力學法，通過洪水分析獲得防洪保護區洪水淹沒過程。

洪水演進計算采用丹麥水動力學研究所開發的MIKE FLOOD軟件，建立整體一、二維嵌套的非恒定流洪水演進數學模型，在這個嵌套模型中，流域內水流運動概化為河道、防洪保護區及其之間的聯系三部分。河道水流采用一維非恒定方法模擬，建立河網模型；防洪保護區概化為二維區域進行模擬；對于河系中控制水流運動的堰、閘、泵、潰口等概化為聯系，聯系的過流能力使用水力學上的相應計算公式。在耦合計算時，一維與二維模型同時計算，在連接處會出現水量交換，可以很好地模擬洪水匯入河道的匯流過程和潰堤洪水從河道泛濫到堤防外的淹沒過程，從而滿足洪水風險圖繪制過程中的功能需求。

（1）一維模型

MIKE 11是用于一維河道內洪水演進過程的模擬，計算范圍為漢口至湖口（八里江）段長江干流，約264km江段，沿江有府澴河、灄水河、倒水河、舉水河、巴水河、浠水河、蘄水、富水、鄱陽湖等支流匯入，支流匯入均按點源方式加入一維模型中。

（2）二維模型

MIKE21是二維平面區域內的水力學計算模型，模擬長江黃石段右堤決口后基于保護區下墊面條件的洪水演進過程，得出保護區的洪水淹沒過程及流場、流速、洪水到達時間、淹沒歷時、洪峰到達時間等數據。計算范圍涉及到黃石市、大冶市、陽新縣和鄂州市的鄂城區，總面積約1224.59m2。

（3）設計洪水

對于長江干流1954年實際洪水，采用長江中下游大湖演進模型，并采用長江中下游歷年來的實測資料進行率定和驗證。上邊界考慮三峽水庫采用優化調度方案的調洪過程，陸水、府澴河、灄水、倒水、舉水、巴水、浠水、蘄水等支流來水采用相應年份的實測過程按照洪量縮放倍比進行縮放后匯入，沿程各河段按照《長江洪水調度方案》（國汛[2011]22號）的洪水調度安排進行調度，得到中下游各主要控制站不同計算工況的水位、流量過程。

（4）堤防潰決條件

綜合考慮河道、防洪保護區地形、地質狀況、工程狀況、行政區劃、歷史出險、洪水風險情況等因素，確定長江干流黃石段堤防決口位于青山湖處，結合長江干流歷史堤防潰決情況，堤防決口寬度定為300m。

由于現狀長江干堤堤防建設標準較高，長江發生1954年實際洪水時，黃石市堤處長江水位沒有達到設計水位，與長江堤頂還有較大距離，約2.7m左右。因此堤防決口時機確定為為到達最高水位的時刻。

3.2 潰堤洪水風險計算

通過洪水演進計算，長江干流發生1954年實際洪水時，黃石段在青山湖處決口，決口寬度300m，決口時間為1954年7月31日3：00，最大進洪流量分別為4997m3/s，堤外保護區內總進洪水量分別為10.7億m3。

堤防決口后，潰堤洪水進入黃石市區，瞬時峰量達4997m3/s，且流速較大，潰堤后6小時內，洪水主要淹沒黃石市區，慈湖、花馬湖湖區均受淹，由于淹沒水深較大，達到4m左右，城區內民房、建筑物、交通設施等財產受到嚴重威脅。之后，受黃石城區周邊地形及黃石市區建筑物及保護區內洪水位的頂托影響，潰堤洪水主要向鄂州方向演進，潰堤后42個小時，洪水沿江邊地勢較低處進入鄂州城區。隨著潰堤洪水的持續演進，洪水大部分進入鄂州市，至青山附近受地形影響不再向北影響，潰堤48小時后，保護區內（含鄂州）受淹面積增加緩慢，主要在鄂州城區以南方向及花馬湖周邊擴展。

4、堤防決口經濟損失估算

4.1 洪災影響分析

耕地、公路、鐵路等基礎圖層、經濟數據圖層以及洪水風險圖層等基于GIS平臺進行空間疊加分析，得到長江干流1954年實際洪水的黃石段決口堤防內防洪保護區淹沒面積、淹沒農田面積、受影響GDP、受影響人口、受影響公路和受影響鐵路等的統計結果。

堤外防洪保護區受淹主要行政區域為團城山街道、花湖街道、黃石港街道，受淹行政區11個，淹沒面積34.66km2，淹沒居民地面積263.85萬m2，淹沒耕地面積17hm2，受影響公路長度101.84km，受影響鐵路長度20.6km，受影響的重點單位及設施585個。

受淹行政區總人口29.44萬人，受影響人口15.16萬人。受影響人口最大的黃石港區勝陽港街道，其中1.0～2.0m淹沒水深下受影響人口為0，2.0～3.0m淹沒水深下受影響人口為0.08萬人，3m以上淹沒水深下受影響人口為4.11萬人。受淹行政區總GDP為285.01億元，受影響GDP為118.14億元。

4.2 財產損失率

洪災損失評估的關鍵是估算不同淹沒水深（歷時）條件下，各類財產洪災損失率，建立淹沒水深（歷時）與各類財產洪災損失率關系。洪災損失率估算利用前人調查研究的資料，并借鑒相似地區資料，及防洪保護區的歷史洪災損失資料予以確定。水深分級按0.05～1.0m、1.0～2.0m、2.0～3.0m、3.0m以上四級劃分，淹沒歷時按小時計。采用經驗數據，如表1所示。

4.3 洪災損失計算結果

在確定了各類承災體受淹程度、災前價值之后，根據洪災損失率關系，即可進行分類洪災直接經濟損失估算，本次洪災損失分別對城鄉居民住房財產損失、城鄉工礦及商業企業損失、農林牧漁業損失、交通等進行估算。

青山湖潰口1954年實際洪水洪災損失為155.58億元，其中居民房屋損失為10.26億元，家庭財產損失為12.9億元，農業損失為0.07億元，工業資產損失為30.57億元，工業產值損失為56.43億元，商貿業資產損失為25.71億元，商貿業主營收入損失為17.79億元，道路損失為1.47億元，鐵路損失為0.4億元。

洪災損失合計最大的行政區為下陸區團城山街道，共損失34.76億元。居民房屋損失最大的行政區為黃石港區黃石港街道，損失2.8億元；家庭財產損失最大的行政區為黃石港區黃石港街道，損失3.84億元；農業損失最大的行政區為下陸區團城山街道，損失0.06億元；工業資產損失最大的行政區為黃石港區黃石港街道，損失8.4億元；工業產值損失最大的行政區為下陸區團城山街道，損失18.36億元；商貿業資產損失最大的行政區為下陸區團城山街道，損失8.29億元；商貿業主營收入損失最大的行政區為下陸區團城山街道，損失5.79億元；道路損失最大的行政區為黃石港區黃石港街道，損失0.44億元；鐵路損失最大的行政區為西塞山區澄月街道，損失0.19億元。

5、結語

地理信息系統已經廣泛的應用于水利領域的各個角落，通過洪水的情景演進數學模擬推求各類洪水風險信息基礎數據圖層，利用地理信息系統的空間分析功能，可以迅速實現各種洪水災情一定程度上的監測、災害面積的確定、災情的預估等，從而為防洪減災提供及時、準確的信息，為地方洪水風險管理及防洪調度決策提供可靠的技術參考。

參考文獻：

[1] 李娜，向立云，程曉陶. 國外洪水風險圖制作比較及對我國洪水風險圖制作的建議[J]。水利發展研究，2005（6）：28-32。

[2] 曹永強，黃林顯，苗迪，等. 我國洪水風險圖繪制方法研究[J].人民黃河，2008（8）：6-7。

[3] Zhang W， Zhou J， Liu Y， et al. Emergency evacuation planning against dike-break flood： a GIS-based DSS for flood detention basin of Jingjiang in central China[J]. Natural Hazards， 1-19.