王紅旗，馮靖雲(yún)，范子武

（1.華南農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院，廣州 510642；2.南京水利科學(xué)研究院，南京 210029）

0 引 言

蓄滯洪區(qū)是江河流域防洪減災(zāi)體系中的重要組成部分，是犧牲局部保全局，最大程度減小洪水災(zāi)害損失的重要舉措。國(guó)內(nèi)外學(xué)者從建設(shè)與管理、生態(tài)修復(fù)、洪水演進(jìn)數(shù)值模擬與風(fēng)險(xiǎn)圖繪制、避險(xiǎn)轉(zhuǎn)移等多個(gè)方面開(kāi)展了蓄滯洪區(qū)的研究［1-4］。王小笑等［5］構(gòu)建了二維數(shù)學(xué)模型，利用GIS 反演技術(shù)對(duì)蓄滯洪區(qū)進(jìn)行危險(xiǎn)程度分區(qū)并編制洪水風(fēng)險(xiǎn)圖，為蓄滯洪區(qū)管理提供科學(xué)依據(jù)；Qiang Liu 等［6］提出有限體積法進(jìn)行數(shù)值求解河道與蓄滯洪區(qū)一二維耦合水動(dòng)力模型，更真實(shí)模擬洪水演進(jìn)過(guò)程；果鵬等［7］提出受淹對(duì)象失穩(wěn)機(jī)制的洪水風(fēng)險(xiǎn)分析方法，評(píng)估蓄滯洪區(qū)的洪水風(fēng)險(xiǎn)等級(jí)，具有更好適用性。目前，蓄滯洪區(qū)的洪水風(fēng)險(xiǎn)分析逐漸成為眾多學(xué)者的研究重點(diǎn)［8，9］。

飛來(lái)峽水庫(kù)作為北江防洪體系的關(guān)鍵性工程，通過(guò)與潖江蓄滯洪區(qū)聯(lián)合運(yùn)用，可將下游北江大堤防洪標(biāo)準(zhǔn)由100年一遇提高到300年一遇，確保廣州、佛山等城市的防洪安全。飛來(lái)峽水庫(kù)自1999年投入運(yùn)行以來(lái)，先后成功抗擊了“05·6”珠江流域超100年一遇特大洪水，“06·7”和“13·8”北江近50年一遇洪水。飛來(lái)峽庫(kù)區(qū)人口總數(shù)多、分布廣，涉及的社會(huì)經(jīng)濟(jì)總量大，當(dāng)啟用庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)調(diào)洪，將產(chǎn)生不同程度的淹沒(méi)損失。在調(diào)蓄300年一遇洪水的情況下，飛來(lái)峽庫(kù)區(qū)臨時(shí)淹沒(méi)影響人口約27.7萬(wàn)人，淹沒(méi)耕地7 633 hm2。本文以飛來(lái)峽水庫(kù)為主要研究對(duì)象，構(gòu)建水力學(xué)模型模擬洪水演進(jìn)過(guò)程，分析了沿程斷面最高水位、洪水淹沒(méi)面積與洪峰抵達(dá)時(shí)間，分析了庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)不同年型與頻率的洪水風(fēng)險(xiǎn)，有利于增強(qiáng)臨時(shí)蓄滯洪區(qū)居民及管理工作人員的風(fēng)險(xiǎn)意識(shí)。

1 飛來(lái)峽庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)概況

飛來(lái)峽水庫(kù)位于廣東省北江干流中游清遠(yuǎn)市境內(nèi)，壩址控制流域面積34 097 km2，占北江流域面積的73%，總庫(kù)容19.04億m3，防洪庫(kù)容13.36 億m3，正常庫(kù)容4.23 億m3。庫(kù)區(qū)臨時(shí)淹沒(méi)涉及英城、波羅坑、連江口和社崗4 處蓄滯洪區(qū)，各臨時(shí)蓄滯洪區(qū)概況如表1。

表1 飛來(lái)峽庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)概況Tab.1 Overview of the temporary flood storage and detention area in Feaixia Reservoir

2 模型方法

本文對(duì)河道采用一維水動(dòng)力建模，對(duì)蓄滯洪區(qū)進(jìn)行二維水動(dòng)力建模，一維河道和二維洪泛區(qū)在堤岸處進(jìn)行耦合模擬洪水演進(jìn)過(guò)程。

2.1 河道內(nèi)洪水演進(jìn)模擬

在一維單寬水流情況下，圣維南方程組描述為：

式中：B為河面總寬度，m；H為水位，m；Q為流量，m3/s；s0為旁側(cè)入流流量，m3/s；A 為斷面面積，m2；Vx為旁側(cè)入流在水流方向上的流速分量，m/s；K 為流量模數(shù)；α 為動(dòng)量修正系數(shù)。對(duì)上式采用HLL格式的黎曼近似解計(jì)算一維模型通量。

2.2 洪泛區(qū)內(nèi)洪水演進(jìn)模擬

洪水在洪泛區(qū)的水流運(yùn)動(dòng)具有明顯的二維特征，采用守恒型式的洪流方程作為控制方程：

2.3 一二維耦合洪水模擬

一維河道與二維洪泛區(qū)水流運(yùn)動(dòng)控制方程的統(tǒng)一形式表示為：

式中：D 為系數(shù)矩陣；?F 表示變量F 對(duì)空間變量的偏導(dǎo)數(shù)和；S為源項(xiàng)；Ω為控制體范圍；t為時(shí)間變量；Ts為持續(xù)時(shí)間。對(duì)上式采用中心形式的有限體積法求解一二維耦合模型［11］。

3 模型的構(gòu)建及驗(yàn)證

3.1 模型的構(gòu)建范圍

本次模型的構(gòu)建范圍為飛來(lái)峽水庫(kù)庫(kù)區(qū)，北起北江干流馬競(jìng)寮站，南至水庫(kù)壩址，東起連江支流高道站，西至滃江支流長(zhǎng)湖壩下，總面積606.6 km2。其中一維動(dòng)庫(kù)調(diào)洪模型范圍從馬競(jìng)寮至壩址河段總長(zhǎng)約95.8 km；二維調(diào)洪模型的范圍為飛來(lái)峽水庫(kù)壩址以上的防洪保護(hù)區(qū)，包括英德、波羅坑、連江口及社崗4處臨時(shí)蓄滯洪區(qū)。

3.2 邊界條件

模型上邊界條件為馬徑寮站、高道站、長(zhǎng)湖站的實(shí)測(cè)流量；下邊界條件為按照飛來(lái)峽水庫(kù)調(diào)度規(guī)則確定的水位流量關(guān)系曲線。

3.3 模型率定與驗(yàn)證

3.3.1 參數(shù)率定

采用非恒定流變步長(zhǎng)法進(jìn)行多次試算，并綜合考慮模型收斂及計(jì)算時(shí)間等因素，對(duì)模型的計(jì)算參數(shù)進(jìn)行設(shè)定。北江干支流不同河段經(jīng)過(guò)調(diào)整優(yōu)化后的糙率取值見(jiàn)表2。

表2 北江流域干支流不同河段糙率取值Tab.2 Roughness value of different river sections of main tributaries in Beijiang River Basin

利用“68·6”和“82·5”兩場(chǎng)典型洪水進(jìn)行河道糙率率定，率定過(guò)程及最終模擬成果如下：“68·6”型洪水洪峰流量實(shí)測(cè)值為15 000 m3/s，計(jì)算洪峰流量為14 787 m3/s，相對(duì)誤差為1.42%；“82.5”型洪水實(shí)測(cè)洪峰流量為18 000 m3/s，計(jì)算洪峰流量為17 307 m3/s，相對(duì)誤差為3.85%。兩場(chǎng)年型洪水橫石站計(jì)算流量過(guò)程與實(shí)測(cè)流量過(guò)程基本一致，見(jiàn)圖1。

圖1 橫石站實(shí)測(cè)和計(jì)算的流量過(guò)程線比較Fig.1 Comparison of flow process line measured and calculated at Hengshi station

3.3.2 模型驗(yàn)證

模型的驗(yàn)證利用了北江流域2006年7月和2008年6月發(fā)生的兩場(chǎng)洪水的實(shí)測(cè)資料，驗(yàn)證結(jié)果見(jiàn)表3。

表3 干流主要水位站點(diǎn)驗(yàn)證結(jié)果Tab.3 Verification results of main water level stations in the main stream

模型驗(yàn)證結(jié)果表明，干流主要站點(diǎn)模擬水位的相對(duì)誤差均小于3%。其中2008年6月場(chǎng)次洪水驗(yàn)證中，連江口絕對(duì)誤差達(dá)0.59 m，與水位測(cè)量存在誤差有關(guān)?？傮w而言，模擬成果具有較高精度，所構(gòu)建的模型及設(shè)定的計(jì)算參數(shù)合理可行［12，13］。

4 洪水風(fēng)險(xiǎn)分析

根據(jù)率定的北江流域干支流不同河段糙率參數(shù)值（見(jiàn)表2），選取北江流域發(fā)生的兩場(chǎng)典型洪水“68·6”型和“82·5”型進(jìn)行動(dòng)庫(kù)調(diào)洪演算。對(duì)該兩場(chǎng)典型洪水分別設(shè)計(jì)300、200、100和50年一遇共8 種工況，選取沿程英德大橋、盲仔峽、連江口以及壩址等典型斷面，統(tǒng)計(jì)其沿程斷面水位過(guò)程、洪水淹沒(méi)面積以及洪峰抵達(dá)時(shí)間。結(jié)合基于DOM 及DEM 疊加圖統(tǒng)計(jì)的庫(kù)區(qū)人口與房屋調(diào)查數(shù)據(jù)，分析4處臨時(shí)蓄滯洪區(qū)的淹沒(méi)風(fēng)險(xiǎn)。

4.1 人口與房屋淹沒(méi)分析

在“68·6”和“82·5”兩種年型，重現(xiàn)期為50年、100年、200年和300年等不同頻率洪水條件下，各臨時(shí)蓄滯洪區(qū)所在處典型斷面最高水位統(tǒng)計(jì)見(jiàn)表4。各臨時(shí)蓄滯洪區(qū)涉及鄉(xiāng)鎮(zhèn)人口及房屋淹沒(méi)統(tǒng)計(jì)見(jiàn)表5。

表4 不同工況下沿程斷面最高水位 mTab.4 The highest water level along the section under different working conditions

表5 不同水位高程庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)涉及鄉(xiāng)鎮(zhèn)人口及房屋淹沒(méi)調(diào)查統(tǒng)計(jì)表Tab.5 Statistical table of temporary flood storage and detention areas in different water level elevation reservoir areas involving township population and housing flooding

可以看出，不同年型、不同重現(xiàn)期洪水沿程最高水位均呈相同的變化趨勢(shì)：不同設(shè)計(jì)洪水依次經(jīng)過(guò)英德大橋、盲仔峽、連江口及壩址處，沿程水位逐漸降低；相同年型下，洪水頻率越小，沿程斷面最高水位越大。以“68·6”年型300年一遇為例，壩址處最高洪水位達(dá)32.6 m，比水庫(kù)防洪高水位31.17 m 超高1.43 m，距校核洪水位33.17 m 僅0.57 m，啟用臨時(shí)蓄滯洪區(qū)防洪十分迫切。結(jié)合調(diào)查統(tǒng)計(jì)數(shù)據(jù)，洪水重現(xiàn)期越大，各城鎮(zhèn)的淹沒(méi)房屋面積與受影響人口越多，洪災(zāi)損失越嚴(yán)重。其中，連江口鎮(zhèn)所在地高程較低，其所處的連江口臨時(shí)蓄滯洪區(qū)最先遭受洪水威脅；當(dāng)遭遇100年一遇以上洪水時(shí)，英德鎮(zhèn)淹沒(méi)房屋面積及受影響人口倍增，面臨最嚴(yán)重的洪災(zāi)損失。

4.2 淹沒(méi)范圍影響分析

通過(guò)模擬淹沒(méi)分別得到8種設(shè)計(jì)工況的洪水淹沒(méi)面積見(jiàn)表6。圖2 給出了“68·6”年型50年、100年、200年和300年一遇洪水時(shí)，飛來(lái)峽庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)的淹沒(méi)范圍示意圖。

圖2 “68·6”年型不同重現(xiàn)期洪水淹沒(méi)范圍示意圖Fig.2 Schematic diagram of the“68·6”flood

表6 模擬洪水淹沒(méi)面積統(tǒng)計(jì)表Tab.6 Statistical table of simulated flood inundation areas

從淹沒(méi)數(shù)據(jù)及范圍示意圖可看出，兩種年型洪水的淹沒(méi)面積都隨著洪水重現(xiàn)期的增大而增加。當(dāng)發(fā)生300年一遇洪水時(shí)，整個(gè)庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)的洪水淹沒(méi)面積最大，在“68·6”年型洪水下，洪水淹沒(méi)面積最大可達(dá)290 km2。

4.3 洪峰抵達(dá)時(shí)間分析

對(duì)8 種設(shè)計(jì)工況條件下，以河道內(nèi)洪峰到達(dá)白石窯處的時(shí)刻為起始點(diǎn)，分別統(tǒng)計(jì)各設(shè)計(jì)工況條件下洪水由白石窯至英德大橋、盲仔峽、連江口及壩址的用時(shí)，見(jiàn)表7。

表7 典型斷面洪峰抵達(dá)時(shí)間統(tǒng)計(jì)表 hTab.7 Statistics of flood peak arrival time of typical section

從結(jié)果可看出，“68·6”型洪水洪峰抵達(dá)壩址時(shí)間較短，而“82·5”型洪水洪峰抵達(dá)壩址時(shí)間略長(zhǎng)。同時(shí)年型洪峰抵達(dá)時(shí)間隨洪水頻率的減小而縮短，如300年一遇洪水要比50年一遇洪水快2～3 h，即大洪水傳播速度更快，臨時(shí)蓄滯洪區(qū)啟用時(shí)間越緊張。當(dāng)遭遇300年一遇洪水，洪峰抵達(dá)英城臨時(shí)蓄滯洪區(qū)只需3.0～4.5 h，抵達(dá)波羅坑臨時(shí)蓄滯洪區(qū)需要5.5～6.5 h，抵達(dá)連江口蓄滯洪區(qū)需要6.0～7.0 h，抵達(dá)社崗臨時(shí)蓄滯洪區(qū)僅需10.0～12.0 h，洪峰抵達(dá)時(shí)間之快使蓄滯洪區(qū)調(diào)度困難，也給人員及時(shí)轉(zhuǎn)移帶來(lái)極大的挑戰(zhàn)。

5 結(jié)論與建議

（1）基于飛來(lái)峽庫(kù)區(qū)地形地貌及北江干支流洪水特征，構(gòu)建了洪水演進(jìn)水動(dòng)力數(shù)學(xué)模型。采用“68·6”和“82·5”兩場(chǎng)典型洪水率定模型參數(shù)，通過(guò)流量過(guò)程線比較，洪峰流量誤差較小；并利用2006年7月和2008年6月兩場(chǎng)實(shí)測(cè)洪水資料進(jìn)行模型驗(yàn)證，通過(guò)與干流主要水位站點(diǎn)的實(shí)測(cè)水位比較，水位誤差較小。所構(gòu)建的洪水演進(jìn)模型參數(shù)設(shè)置合理，具有較高的精度，模擬結(jié)果可用于飛來(lái)峽庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)洪水淹沒(méi)風(fēng)險(xiǎn)分析。

（2）運(yùn)用一維、二維及水力學(xué)耦合模型，模擬飛來(lái)峽庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)洪水演進(jìn)模型，設(shè)計(jì)“68.6”和“82.5”年型重現(xiàn)期分別為300年、200年、100年及50年共8種不同工況。以“68.6”年型300年一遇為例，沿程各典型斷面均超過(guò)32 m，各臨時(shí)蓄滯洪區(qū)人口與房屋均受不同程度的淹沒(méi)影響，其中連江口臨時(shí)蓄滯洪區(qū)地勢(shì)較低洼，最先遭遇洪水風(fēng)險(xiǎn)；英城蓄滯洪區(qū)工農(nóng)業(yè)總產(chǎn)值65.23 億元，其中英城街道是英德市政治、經(jīng)濟(jì)、文化、教育與信息中心，一旦啟用英城臨時(shí)蓄滯洪區(qū)，轉(zhuǎn)移遷安難度大且成本高［14-16］；整個(gè)飛來(lái)峽庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)洪水淹沒(méi)面積高達(dá)290 km2，洪峰抵達(dá)飛來(lái)峽庫(kù)區(qū)各臨時(shí)蓄滯洪區(qū)的時(shí)間不超過(guò)12 h，應(yīng)及時(shí)啟動(dòng)英城、波羅坑、連江口及社崗臨時(shí)蓄滯洪區(qū)安全轉(zhuǎn)移預(yù)案。

（3）近年來(lái)異常天氣多發(fā)，極端暴雨洪水事件頻現(xiàn)，在北江流域遭遇大洪水時(shí)，為減輕下游防洪壓力及保護(hù)重要城市安全，考慮啟用庫(kù)區(qū)臨時(shí)蓄滯洪區(qū)滯納部分洪水，區(qū)域內(nèi)受淹沒(méi)風(fēng)險(xiǎn)在所難免。因此，在非工程措施方面需要定期開(kāi)展洪水情景模擬，調(diào)度時(shí)充分考慮臨時(shí)蓄滯洪區(qū)的洪災(zāi)損失程度，并構(gòu)建飛來(lái)峽水庫(kù)實(shí)時(shí)預(yù)報(bào)調(diào)度模型，提高洪水預(yù)見(jiàn)性和科學(xué)決策水平，對(duì)淹沒(méi)風(fēng)險(xiǎn)進(jìn)行預(yù)判以提高風(fēng)險(xiǎn)防范意識(shí)，合理設(shè)計(jì)轉(zhuǎn)移撤退路線，提前預(yù)報(bào)預(yù)警并做好人員轉(zhuǎn)移避險(xiǎn)預(yù)案的演練工作，強(qiáng)化“四預(yù)”措施以減少人員傷亡和財(cái)產(chǎn)損失。