水庫(kù)標(biāo)準(zhǔn)內(nèi)及超標(biāo)準(zhǔn)洪水下泄下游淹沒分析

鄧 歡，周 迎，高 波

（山東省水利勘測(cè)設(shè)計(jì)院,山東 濟(jì)南 250013）

水庫(kù)作為重要的水利工程,在興利除害、防汛減災(zāi)中發(fā)揮重要作用。水庫(kù)按照不同標(biāo)準(zhǔn)洪水進(jìn)行泄洪時(shí),下游可能會(huì)出現(xiàn)潰堤或漫堤等險(xiǎn)情,造成區(qū)域內(nèi)出現(xiàn)洪水淹沒。為科學(xué)的指導(dǎo)防汛搶險(xiǎn)及下游群眾避險(xiǎn)轉(zhuǎn)移,為水庫(kù)防汛調(diào)度提供決策依據(jù),采用科學(xué)的分析方法對(duì)水庫(kù)下游的可能淹沒范圍進(jìn)行洪水淹沒分析是有必要的[1]。本次以山東省日照市莒縣青峰嶺水庫(kù)為例,采用MIKE軟件構(gòu)建了水庫(kù)下游三維數(shù)學(xué)模型,對(duì)水庫(kù)下游淹沒情況進(jìn)行分析。

1 工程概況

青峰嶺水庫(kù)位于日照市莒縣縣城西北30 km處,主壩位于洛河、棋山二鎮(zhèn)交界處的東盧家岔河村。水庫(kù)于1960 年建成,是一座防洪、灌溉為主,發(fā)電、供水等結(jié)合的大（二）型水庫(kù)[2]。

青峰嶺水庫(kù)下游為沭河干流,沭河保護(hù)人口54.5 萬人,保護(hù)耕地48.4 萬畝。沂沭泗河洪水東調(diào)南下續(xù)建工程實(shí)施完成后,沭河的潯河口～湯河口河段防洪標(biāo)準(zhǔn)已達(dá)到20 年一遇,湯河口以下至新沂河河段防洪標(biāo)準(zhǔn)可達(dá)到50 年一遇。沭河中下游防洪保護(hù)區(qū)防洪條件得到較大改善。但沭河上游（潯河口以上段）未進(jìn)行過系統(tǒng)治理,目前堤防斷面不足,險(xiǎn)工險(xiǎn)段較多,防洪標(biāo)準(zhǔn)偏低[3]。

2 基礎(chǔ)資料

2.1 地形資料

分析范圍（見圖2）內(nèi)的1∶10000 地形數(shù)據(jù)、沭河條帶圖等資料。

2.2 水文成果

青峰嶺水庫(kù)100 年一遇,5000 年一遇泄洪過程線,見圖1。

圖1 青峰嶺水庫(kù)泄洪過程線

2.3 防洪設(shè)施

本次淹沒分析范圍內(nèi)防洪設(shè)施主要包括現(xiàn)狀的橋梁、壩、閘（考慮其尺寸及調(diào)度運(yùn)用方式）,沭河及相關(guān)支流的河道堤防（主要考慮其堤頂高程）。

3 模型建立

3.1 模擬的范圍

根據(jù)青峰嶺水庫(kù)下游的地形、沭河及相關(guān)支流的分布情況、堤防相應(yīng)高程、阻水建筑物的尺寸及分布,本次模擬區(qū)域北側(cè)以水庫(kù)大壩為界；南側(cè)以日照市境界為界；東西兩側(cè)計(jì)算范圍取沭河兩側(cè)大致5 km為界并結(jié)合等高線進(jìn)行優(yōu)化。該區(qū)域南北向最寬處為52.3 km,東西向最寬處為15.7 km,總面積603.5 km2,共涉及日照市莒縣、五蓮縣,臨沂市莒南縣等多個(gè)區(qū)縣,模擬范圍見圖2。

圖2 模擬范圍示意圖

3.2 網(wǎng)格的剖分

采用Triangular網(wǎng)格對(duì)模擬區(qū)域進(jìn)行有限元?jiǎng)澐?在地勢(shì)平坦、高程變化不大的區(qū)域, 為縮短模擬時(shí)間適當(dāng)加大了網(wǎng)格間距;在重要的地物附近,對(duì)網(wǎng)格進(jìn)行加密以提高模擬精度。本次模擬采用的網(wǎng)格平均面積為5000 m2,共劃分網(wǎng)格28.8萬個(gè)。

3.3 模擬中地物處理的體現(xiàn)

模型充分考慮模擬范圍內(nèi)高出地面1.0 m的線狀地物,如堤防、鐵路、公路等,并在網(wǎng)格模型中體現(xiàn),以體現(xiàn)其阻水和導(dǎo)流作用。

3.4 模擬中地形插值處理

沭河干流河道內(nèi)高程由實(shí)測(cè)地形圖處理后獲得,周邊區(qū)域高程點(diǎn)采用山東省國(guó)土測(cè)繪院近期調(diào)繪的1/10000 地理信息數(shù)據(jù)。插值后的地形見圖3及圖4。

圖3 插值后的地形三維視圖

圖4 插值后的局部網(wǎng)格

3.5 糙率選取

為反映下墊面類型對(duì)洪水演進(jìn)的影響,對(duì)模擬區(qū)域內(nèi)的河槽、居民地、道路、耕地等各類下墊面分別設(shè)置不同的糙率數(shù),本地區(qū)缺少糙率實(shí)測(cè)資料,故采用《水力計(jì)算手冊(cè)》中的糙率建議值[4],見表1。

表1 模擬區(qū)域內(nèi)糙率取值

3.6 邊界條件

定義青峰嶺水庫(kù)南側(cè)溢洪道出口為流量入口邊界,泄洪流量過程見圖1。給定沭河支流潯河處6月~9月多年平均水位為出口邊界。

3.7 方案匯總

各工況組合見表2。

表2 計(jì)算方案匯總

4 模擬結(jié)果

依據(jù)模型計(jì)算結(jié)果進(jìn)行對(duì)比分析,青峰嶺水庫(kù)按照不同標(biāo)準(zhǔn)進(jìn)行洪水下泄時(shí),洪水演進(jìn)的路徑基本一致,但不同標(biāo)準(zhǔn)下泄的淹沒范圍、淹沒水深、淹沒歷時(shí)等有較大的差異。

4.1 洪水演進(jìn)過程

以青峰嶺水庫(kù)校核洪水標(biāo)準(zhǔn)5000 年一遇洪水標(biāo)準(zhǔn)進(jìn)行泄洪為例,對(duì)下游洪水演進(jìn)歷程及淹沒情況分時(shí)段進(jìn)行分析。各時(shí)段淹沒范圍及淹沒水深見圖5。

圖5 水庫(kù)5000年一遇泄洪洪水演進(jìn)過程

第10 h：青峰嶺水庫(kù)溢洪道下泄流量541 m3/s。洪水在沭河主槽中演進(jìn),基本限制在主河道內(nèi)。

第20 h：青峰嶺水庫(kù)溢洪道下泄流量4320 m3/s。洪水部分漫溢,主要發(fā)生在曲坊社區(qū)至薛家城子村段,從南至北淹沒后站家岔河村、水由村、西地村、章莊村、曲坊社區(qū)、岔河村等。

第30 h：青峰嶺水庫(kù)溢洪道下泄流量3860 m3/s,洪水部分漫溢,上階段所述曲坊社區(qū)至薛家城子村段東側(cè)淹沒范圍變化不大,西側(cè)淹沒大范圍擴(kuò)大、此處涉及村莊多、面積大；何家莊子至前山后村發(fā)生漫溢,但此處淹沒村莊較少。

第40 h：青峰嶺水庫(kù)溢洪道下泄流量2320 m3/s,淹沒范圍達(dá)到最大。曲坊社區(qū)至薛家城子村段往南演進(jìn)至王家課莊,淹沒范圍進(jìn)一步擴(kuò)大；何家莊子至前山后村段淹沒范圍變化不大。

第50 h：青峰嶺水庫(kù)溢洪道下泄流量967 m3/s,淹沒面積逐漸減小,淹沒水深逐漸降低。

第60 h：青峰嶺水庫(kù)溢洪道下泄流量909 m3/s,淹沒面積繼續(xù)減小,淹沒水深繼續(xù)降低。

第70 h：青峰嶺水庫(kù)溢洪道下泄流量186 m3/s,泄洪過程趨于結(jié)束,淹沒面積繼續(xù)減小,淹沒水深繼續(xù)降低。

4.2 最大流速分布

青峰嶺水庫(kù)按5000 年一遇設(shè)計(jì)標(biāo)準(zhǔn)洪水泄洪時(shí),下游各區(qū)域最大流速分布見圖6（a）。

從圖中看出,該標(biāo)準(zhǔn)洪水泄洪時(shí),流速大于0.6 m/s的區(qū)域多聚集在溢洪道至岳家村段,最大流速可達(dá)3.0 m/s。

4.3 淹沒歷時(shí)

青峰嶺水庫(kù)按5000 年一遇標(biāo)準(zhǔn)洪水泄洪時(shí),下游淹沒歷時(shí)見圖6（b）。

圖6 最大流速及淹沒歷時(shí)分布（5000a一遇）

4.4 不同標(biāo)準(zhǔn)洪水下泄淹沒情況對(duì)比

由圖7可以看出,青峰嶺水庫(kù)按不同標(biāo)準(zhǔn)泄洪時(shí),5000年一遇洪水相應(yīng)的淹沒面積與100 年一遇洪水相比,影響范圍在嶠山鎮(zhèn)至小店鎮(zhèn)區(qū)域內(nèi)有所增大,其它區(qū)域影響范圍大致相同；同時(shí)淹沒水深較大的區(qū)域,在高標(biāo)準(zhǔn)洪水較百年一遇標(biāo)準(zhǔn)下泄時(shí)顯著增加。

圖7 青峰嶺水庫(kù)不同標(biāo)準(zhǔn)泄洪下游淹沒范圍圖

5 高風(fēng)險(xiǎn)區(qū)域評(píng)估

通過對(duì)青峰嶺水庫(kù)泄洪時(shí)下游區(qū)域的淹沒水深、行進(jìn)流速、淹沒歷時(shí)等數(shù)據(jù)進(jìn)行分析,識(shí)別出青峰嶺水庫(kù)下游高風(fēng)險(xiǎn)區(qū)域有以下3 處：其中一處位于沭河章村-柳河村段兩岸,該部分地勢(shì)平坦、無有效擋水建筑物,東西跨度大,淹沒范圍廣,涉及村莊40 余個(gè),淹沒深度集中于0.5 m~1.5 m,局地達(dá)到3 m,災(zāi)情嚴(yán)重；一處位于西旺疃村-金墩村,該處范圍小,涉及村莊4個(gè)；一處位于盧家孟堰村-公書社區(qū),該處范圍小,涉及村莊5個(gè)。

6 結(jié)語

本次構(gòu)建了日照市莒縣青峰嶺水庫(kù)下游3 D數(shù)學(xué)模型,通過比較水庫(kù)按照不同標(biāo)準(zhǔn)進(jìn)行洪水下泄時(shí),下游區(qū)域洪水的演進(jìn)過程、淹沒的范圍、淹沒的水深、行進(jìn)流速、總計(jì)淹沒歷時(shí)等數(shù)據(jù),得出以下結(jié)論：青峰嶺水庫(kù)水庫(kù)按照100 年一遇及5000 年一遇標(biāo)準(zhǔn)進(jìn)行洪水下泄時(shí),下游的淹沒區(qū)域均比較大,水深均較大。當(dāng)水庫(kù)需要按設(shè)計(jì)洪水標(biāo)準(zhǔn)以上的洪水標(biāo)準(zhǔn)進(jìn)行泄洪時(shí),應(yīng)充分做好淹沒范圍內(nèi)群眾的轉(zhuǎn)移避險(xiǎn)工作,盡可能的減少人民群眾的生命和財(cái)產(chǎn)損失。

本數(shù)字模擬方法可為水庫(kù)的調(diào)度運(yùn)用及防汛搶險(xiǎn)、轉(zhuǎn)移避險(xiǎn)等提供技術(shù)支撐。