唐克銀，趙 芳，趙 釗，張 晨

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

自1968年美國(guó)國(guó)會(huì)通過(guò)全國(guó)洪水保險(xiǎn)法推動(dòng)洪水風(fēng)險(xiǎn)圖以來(lái)，世界各國(guó)逐漸開始展開洪水風(fēng)險(xiǎn)圖的編制與應(yīng)用工作，形成了歷史洪水調(diào)查法、地貌學(xué)法和非恒定流數(shù)據(jù)模擬法等一系列洪水風(fēng)險(xiǎn)圖編制方法。2006年，歐洲24個(gè)國(guó)家或組織成立了歐洲洪水風(fēng)險(xiǎn)圖交流圈(EXCIMAP)，旨在匯集歐洲所有經(jīng)驗(yàn)和專有知識(shí)并改善洪水風(fēng)險(xiǎn)圖繪制實(shí)踐。20世紀(jì)80年代中期，我國(guó)開始了洪水風(fēng)險(xiǎn)圖的編制研究。1984年，中國(guó)水利水電科學(xué)研究院與海河水利委員會(huì)合作對(duì)永定河洪泛區(qū)洪水演進(jìn)特征進(jìn)行分析，繪制了我國(guó)第一張洪水風(fēng)險(xiǎn)圖。2005年1月，國(guó)家防辦在七大流域部署了36個(gè)試點(diǎn)項(xiàng)目。2013年，國(guó)家防辦在全國(guó)組織實(shí)施重點(diǎn)地區(qū)洪水風(fēng)險(xiǎn)圖編制工作，至今已形成大量的成果，但多為專題研究成果，分析計(jì)算條件單一。文章以洙萬(wàn)片防洪保護(hù)區(qū)洙趙新河右堤曾劉涵洞險(xiǎn)工段發(fā)生50年一遇標(biāo)準(zhǔn)洪水為例，對(duì)復(fù)雜條件下平原保護(hù)區(qū)內(nèi)洪水淹沒(méi)特征值進(jìn)行分析與研究，供洪水風(fēng)險(xiǎn)圖編制工作人員借鑒參考。

1 洙萬(wàn)片防洪保護(hù)區(qū)基本情況

洙萬(wàn)片防洪保護(hù)區(qū)位于魯西南地區(qū)，黃河下游，區(qū)內(nèi)主要為黃河沖積平原，地勢(shì)平坦，土層深厚，屬華北平原新沉降盆地的一部分，自西南向東北呈簸萁形逐漸降低，海拔高度為37～68m，平均地面坡降為1/4700。保護(hù)區(qū)是以洙趙新河右堤、萬(wàn)福河左堤、東魚河北支左堤、黃河右堤及南四湖堤防為界，總面積為2950km2，涉及菏澤市牡丹區(qū)、東明縣、定陶縣、鄆城縣、巨野縣、成武縣，濟(jì)寧市任城區(qū)、嘉祥縣、金鄉(xiāng)縣九個(gè)縣(區(qū))行政區(qū)域的66個(gè)鄉(xiāng)鎮(zhèn)。

2 MIKE軟件計(jì)算方法

2.1 河道一維洪水演進(jìn)計(jì)算方法

MIKE11是一維河道水動(dòng)力數(shù)學(xué)模型，采用6點(diǎn)Abbott-Ionescu有限差分格式對(duì)圣維南方程求解，其連續(xù)方程、動(dòng)量方程分別為

(1)

(2)

式中，A—河道過(guò)水面積；Q—流量；u—側(cè)向流在河道方向的流速；t—時(shí)間；x—沿水流方向的水平坐標(biāo)；q—河道側(cè)向流量；α—?jiǎng)恿啃拚禂?shù)；g—重力加速度；y—水位；Sf—摩阻坡降。

2.2 保護(hù)區(qū)二維洪水演進(jìn)計(jì)算方法

MIKE21是平面二維自由表面流模型，采用隱式交替方向(ADI)逐行法對(duì)連續(xù)方程及動(dòng)量方程分別進(jìn)行時(shí)空上的積分，每個(gè)方向及每個(gè)單獨(dú)的網(wǎng)格線產(chǎn)生的方程矩陣采用追趕法求解，其連續(xù)方程、動(dòng)量方程分別為

(3)

(4)

(5)

式中，H—水深；B—地面高程；Z—水位，Z=H+B；M、N—x和y方向的單寬流量；u、v—x和y方向上的流速分量；n—糙率；g—重力加速度；q—源匯項(xiàng)。

2.3 一、二維耦合計(jì)算方法

通過(guò)MIKE FL00D側(cè)向連接功能，實(shí)現(xiàn)MIKE11的特定河段與MIKE21邊緣網(wǎng)格連接，實(shí)現(xiàn)一維與二維區(qū)域之間自由水體交換，從而模擬河道洪水漫溢后發(fā)生潰決過(guò)程，并在保護(hù)區(qū)內(nèi)模擬洪水演進(jìn)過(guò)程。

3 模型建立與驗(yàn)證

3.1 洙趙新河一維水動(dòng)力模型

洙趙新河一維水動(dòng)力模型采用MIKE11模擬河道行洪演進(jìn)過(guò)程。根據(jù)河道斷面、糙率、上下邊界條件以及河道主要控制建筑物等要素，同時(shí)考慮河道一維與二維的耦合計(jì)算校驗(yàn)潰口參數(shù)和分洪流量等。

(1)河道工況條件

根據(jù)洙趙新河近幾年治理情況以及近期治理規(guī)劃，模型中河道徐河口以下斷面采用設(shè)計(jì)斷面，以上斷面采用治理后的實(shí)測(cè)斷面。

(2)糙率擬定

模型計(jì)算對(duì)計(jì)算精度影響較大的參數(shù)是河道糙率值。山東省水利勘測(cè)設(shè)計(jì)院等有關(guān)單位對(duì)洙趙新河的糙率進(jìn)行了分析驗(yàn)證，其中，主槽糙率為0.0225，灘地糙率為0.03。洙趙新河河道較為規(guī)整順直，計(jì)算時(shí)采用整體糙率0.027。

(3)模型驗(yàn)證

模型驗(yàn)證選取洙趙新河梁山閘2005年實(shí)測(cè)洪水進(jìn)行了參數(shù)率定，并用2006年7月2—8日實(shí)測(cè)洪水對(duì)模型進(jìn)行了驗(yàn)證。根據(jù)擬定的糙率參數(shù)、建筑物布置、實(shí)測(cè)洪水進(jìn)行水位模擬，其中閘壩按暢泄控制，梁山閘實(shí)測(cè)與模擬水位對(duì)照如圖1所示。

圖1 梁山閘下2006年大洪水實(shí)測(cè)、模擬水位對(duì)照

3.2 保護(hù)區(qū)二維水動(dòng)力模型

(1)網(wǎng)格劃分

綜合考慮保護(hù)區(qū)面積、模擬精度、計(jì)算時(shí)間及軟件性能等因素，洪水演進(jìn)分析采用三角形不規(guī)則網(wǎng)格。模型中以流域面積在100km2以上的支流水面線、堤防線和其他河流及高速公路、鐵道等線狀地物的中心線作為網(wǎng)格劃分控制線，并以地形和控制線間距確定劃分精度。對(duì)于河道、堤防、道路及其他地形變化劇烈的區(qū)域，網(wǎng)格適當(dāng)加密，最小網(wǎng)格為0.00007km2。

(2)地形插值

洪水演進(jìn)分析結(jié)果與地形精度密切相關(guān)，模型中地形插值高程點(diǎn)由1∶1萬(wàn)DEM提取。根據(jù)地物復(fù)雜程度提取不同密度的高程散點(diǎn)，村莊散點(diǎn)間距為20m，鐵路、公路、河流、堤防散點(diǎn)間距為5m，其他地物散點(diǎn)間距為50m。

(3)地物概化

洪水演進(jìn)分析需考慮線狀地物對(duì)洪水演進(jìn)的阻水和導(dǎo)流作用，以及沿程缺口、橋涵的過(guò)水作用。模型中地物概化主要考慮高于地面0.5m以上的線狀地物，包括省道、高速、鐵路及堤防，線狀地物沿程缺口及橋涵。

(4)糙率設(shè)置

根據(jù)《洪水風(fēng)險(xiǎn)圖編制細(xì)則》，糙率取值一般應(yīng)利用實(shí)測(cè)洪水資料進(jìn)行率定，無(wú)實(shí)測(cè)資料的地區(qū)可根據(jù)《水力學(xué)手冊(cè)》確定，或參考采用相似條件地區(qū)的糙率。本地區(qū)缺少潰堤洪水演進(jìn)實(shí)測(cè)資料，故采用《水力學(xué)手冊(cè)》中的建議值。對(duì)保護(hù)區(qū)內(nèi)的村莊、道路、耕地、河流等地物設(shè)置不同的糙率，以反映保護(hù)區(qū)下墊面對(duì)洪水演進(jìn)的影響，各地物糙率取值見表1。

表1 保護(hù)區(qū)內(nèi)各地物糙率參照表

(5)排水設(shè)置

根據(jù)當(dāng)?shù)貙?shí)際內(nèi)澇風(fēng)險(xiǎn)特點(diǎn)，模型在洪水退水階段考慮閘門及泵站的外排作用。依據(jù)全國(guó)第一次水利大普查及現(xiàn)場(chǎng)調(diào)查成果，部分泵站位于保護(hù)區(qū)內(nèi)部，遠(yuǎn)離支流及邊界，另有部分泵站盡管靠近支流，但無(wú)法發(fā)揮外排功能，模型計(jì)算考慮規(guī)模以上泵站50余座。

(6)蒸發(fā)與下滲

綜合本區(qū)域下墊面情況、氣候及實(shí)測(cè)資料，將蒸發(fā)與下滲一并考慮。模型計(jì)算采用的蒸發(fā)與下滲值取為8mm/d。

4 洪水演進(jìn)結(jié)果分析

4.1 洙趙新河一維洪水分析成果

通過(guò)MIKE FLOOD側(cè)向連接方式進(jìn)行一維、二維模型耦合計(jì)算，獲取一維河道分洪流量作為二維模型潰口入流流量。洙趙新河20年一遇和50年一遇標(biāo)準(zhǔn)洪水各潰口分洪流量過(guò)程如圖2所示。

圖2 洙趙新河標(biāo)準(zhǔn)洪水各潰口分洪流量過(guò)程圖

4.2 保護(hù)區(qū)二維洪水分析

(1)內(nèi)澇處理

將MIKE SHE計(jì)算得到的內(nèi)澇成果折算為均勻分布的降雨，在MIKE21中重新演算。澇水在平面上流動(dòng)趨向東，受線狀地物的阻攔及分隔，主要積聚在低洼處形成較大的水面。

(2)潰堤洪水演進(jìn)

潰堤起始時(shí)刻，潰口進(jìn)洪量達(dá)225.7m3/s，最大水深為1.49m，洪水在嘉祥縣境內(nèi)沿省道S252向南行進(jìn)，最遠(yuǎn)到蔡河，部分洪水匯入蔡河，沿河道流向下游；潰堤12h，洪水沿蔡河向下游行進(jìn)，在金鄉(xiāng)縣滿硐鄉(xiāng)東部、胡集鎮(zhèn)北部匯集，最大水深達(dá)1.2m；潰堤24h，洪水匯集于滿硐鄉(xiāng)東部、胡集鎮(zhèn)北部，形成較大積水區(qū)；潰堤2d，洙趙新河干流洪峰消退，潰堤過(guò)程基本結(jié)束；潰堤6d，洙趙新河干流洪峰消退，總進(jìn)洪量達(dá)6592.32萬(wàn)m3，受地形西高東低影響，洪水沿蔡河緩慢東流，匯集于蔡河下游，水深超過(guò)2m，洪水演進(jìn)過(guò)程如圖3所示。

圖3 50年一遇洪水演進(jìn)過(guò)程

表2 50年一遇潰堤洪水各時(shí)刻進(jìn)洪流量

(3)洪水特征值

洪水特征值主要包括洪水水位、流速、流量等。通過(guò)在洪水演進(jìn)路線上設(shè)置水位監(jiān)測(cè)點(diǎn)(如圖4所示)，對(duì)局部區(qū)域水位過(guò)程進(jìn)行分析。潰堤后，各監(jiān)測(cè)點(diǎn)水位均經(jīng)歷兩個(gè)上漲和一個(gè)消退過(guò)程，即因降雨和潰堤洪水造成的上漲和支流排水、蒸發(fā)下滲及洪水演進(jìn)的消退。水位上漲的程度和速度大于水位下降的程度和速度，符合洪水傳播機(jī)制。

圖4 水位監(jiān)測(cè)點(diǎn)設(shè)置分布圖

因預(yù)設(shè)水位監(jiān)測(cè)點(diǎn)與潰口遠(yuǎn)近關(guān)系，水位依次上漲，其中潰口處監(jiān)測(cè)點(diǎn)t1洪峰最早出現(xiàn)，水位在潰堤時(shí)刻急劇上升2m以上，但是消退較快，歷時(shí)較短；監(jiān)測(cè)點(diǎn)t2、t3、t4為洼地或者演進(jìn)路線點(diǎn)，距離潰口較遠(yuǎn)，洪水流速較慢，消退時(shí)間長(zhǎng)；t5、t6是洪水最終匯集地，距離潰口二三十千米，尤其是監(jiān)測(cè)點(diǎn)t6水位上漲較慢，且受集中排水作用，有明顯的水位消退過(guò)程。

潰口最大進(jìn)洪流量出現(xiàn)在潰決后第23h，為478.4m3/s。受兩山形成的狹窄通道影響，最大流速出現(xiàn)在潰口下游250m范圍內(nèi)，為1.065m/s；潰口外1.2km的范圍內(nèi)，流速迅速降至0.5m/s；潰口下游2.0km，流速降到0.3m/s以下。

(4)淹沒(méi)面積

根據(jù)淹沒(méi)結(jié)果統(tǒng)計(jì)，洪水影響范圍內(nèi)，水深大于0.5m的最大淹沒(méi)面積為132.55km2，占保護(hù)區(qū)總面積的8.38%。洪水退水以后，保護(hù)區(qū)內(nèi)尚有部分洪水無(wú)法排出，洪水淹沒(méi)范圍及歷時(shí)如圖5所示。

圖5 保護(hù)區(qū)洪水淹沒(méi)范圍及歷時(shí)圖

對(duì)保護(hù)區(qū)內(nèi)模擬發(fā)生50年一遇暴雨內(nèi)澇洪水的淹沒(méi)面積與1957年歷史記錄大洪水淹沒(méi)范圍圖進(jìn)行對(duì)比分析，發(fā)現(xiàn)其主要淹沒(méi)區(qū)域、范圍、水深等特征基本一致。

5 結(jié)論

(1)基于MIKE11建立的洙趙新河一維水動(dòng)力模型，采用2005年實(shí)測(cè)洪水進(jìn)行了參數(shù)率定，用2006年實(shí)測(cè)洪水對(duì)模型進(jìn)行了驗(yàn)證，河道實(shí)測(cè)水位和模擬水位存在一定誤差，模擬結(jié)果與實(shí)測(cè)或調(diào)查數(shù)據(jù)基本吻合。

(2)基于MIKE21建立的保護(hù)區(qū)二維水動(dòng)力模型，利用1∶10000地形圖和全國(guó)第一次水利大普查數(shù)據(jù)，對(duì)保護(hù)區(qū)內(nèi)構(gòu)筑物進(jìn)行概化建模，計(jì)算結(jié)果與1957年歷史大洪水淹沒(méi)范圍圖進(jìn)行對(duì)比分析，模擬結(jié)果與實(shí)測(cè)或調(diào)查數(shù)據(jù)基本吻合。

(3)通過(guò)洪水淹沒(méi)特征值分析與對(duì)比，模型對(duì)保護(hù)區(qū)的概化合理，能從整體上反映保護(hù)區(qū)發(fā)生暴雨或河道洪水引起的淹沒(méi)范圍分布。率定和驗(yàn)證后的模型可以用于其他方案的洪水分析計(jì)算，為該區(qū)域洪水影響分析和洪水風(fēng)險(xiǎn)圖編制工作提供依據(jù)。