MIKE系列模型在蓄滯洪區(qū)洪水模擬中的應用研究

劉曉琴，劉國龍，王 振

(河南省水利勘測設計研究有限公司，河南 鄭州 450016)

0 引 言

蓄滯洪區(qū)作為大江大河大洪水時的分滯洪區(qū)域，是江河防洪體系中重要的組成部分，在防汛應急中發(fā)揮了一定的削峰滯洪作用，是防洪調度的重要手段和有效措施[1,2]。在我國洪水災害頻發(fā)[3]的情況下，蓄滯洪區(qū)遭受洪水淹沒頻率較高，準確合理的蓄滯洪區(qū)洪水模擬可為防洪調度的科學決策提供支撐；為洪水資源化利用提供基礎；為流域防洪規(guī)劃，以及交通、石油、通訊、電力等工程穿越蓄滯洪區(qū)的影響評價和工程設計等[4-6]提供參考；同時對指導蓄滯洪區(qū)安全建設、避險轉移和洪水保險服務等具有重要意義。

丹麥DHI公司研發(fā)的MIKE系列模型是國際上應用較為廣泛的一款商業(yè)水力學計算軟件，曾在丹麥、埃及以及中國等地區(qū)得到了廣泛使用和成功驗證[7-9]。通常采用MIKE11構建行洪河道一維模型，MIKE21構建蓄滯洪區(qū)洪水演進二維模型，MIKE FLOOD構建一二維動態(tài)耦合模型。本文針對多條河流自然漫溢進洪、干流河道上控制閘控制泄洪的蓄滯洪區(qū)，以楊莊蓄滯洪區(qū)為例，基于實際DEM數(shù)據(jù)，通過局部加密網(wǎng)格將河溝概化到MIKE21模型中，模擬洪水進洪及演進過程；基于實際調度運用方式，通過MIKE11模型中可控建筑物的設置，模擬退洪過程；采用MIKE FLOOD進行一二維動態(tài)耦合，針對關鍵技術問題進行研究，以期為蓄滯洪區(qū)模型的構建和MIKE系列模型的應用提供參考。

1 MIKE系列模型原理與方法

1.1 MIKE11一維河道水力學模型

MIKE11水動力學模型假定水流是不可壓縮、均質流體，流動過程基本是一維流態(tài)，河道坡降小、縱向斷面變化幅度小，服從靜水壓力分布。采用一維明渠非恒定流微分方程，即圣維南(Saint-Venant)方程組。

連續(xù)方程：

(1)

動量方程:

(2)

式中：t為時間坐標；x為河道沿程坐標；Q為流量；Z為水位；A為過水斷面面積；B為水面寬度；K為流量模數(shù)；g為重力加速度；q為旁側入流流量。

(3)

式中：n為河道糙率；R為水力半徑。

MIKE11采用六點Abbott-Ionescu有限差分格式對圣維南方程組進行離散。離散格式在每一個網(wǎng)格點并不同時計算水位和流量，而是按順序交替計算水位或流量，分別稱為h點和Q點。通過其上可控建筑物的設置，模擬各類可控建筑物調度工況。

1.2 MIKE21二維洪水演進模型

MIKE21水動力學模型依據(jù)水流運動的二維非恒定流方程組，水流連續(xù)性方程、水流沿x方向的動量方程及水流沿y方向的動量方程：

(4)

以上方程組利用迭代法求解即可得到每一時刻在(x，y)處的水位z、水深h以及x，y方向的流速u，v。模型中考慮了各線狀阻水如堤防、公路、鐵路以及區(qū)內涵閘影響，洪水過涵洞采用同一維流相同的公式計算，洪水漫堤采用寬頂堰公式計算。

1.3 MIKE FLOOD 一二維耦合模型

MIKE FLOOD是把一二維模型連接在一起，進行動態(tài)耦合的模型系統(tǒng)[10-12]，耦合模型本身不進行水力計算，只是通過耦合方式、位置及相關參數(shù)設置，使一二維模型實時交換數(shù)據(jù)、協(xié)同工作、靈活并穩(wěn)定計算。MIKE11和MIKE21基于MIKE FLOOD可建立5種不同的耦合連接形式，即標準連接、側向連接、建筑物連接、零流動連接和側向建筑物連接，蓄滯洪區(qū)一二維耦合模型常采用標準連接方式。

對于非結構化網(wǎng)格，標準連接是將連接線映射到一個或多個二維網(wǎng)格的邊上(稱為耦合線)，MIKE11模型為MIKE21模型提供流量值Q作為二維模型的邊界條件，將Q值分布到MIKE21計算單元的各節(jié)點上，然后MIKE21將耦合單元格邊上的平均水位返回到MIKE11模型中，較好地模型洪水匯入河道的匯流過程和洪水從河道漫溢到蓄滯洪區(qū)的淹沒過程。

2 楊莊蓄滯洪區(qū)洪水模擬

2.1 區(qū)域概況及計算方案

楊莊蓄滯洪區(qū)位于駐馬店西平縣城西20 km的小洪河干流上，由東部和南部大壩、西部及北部自然高地合圍而成，總面積82 km2，設計蓄洪水位71.54 m、蓄洪量2.03 億m3。區(qū)內主要河道有小洪河、洪溪河和吉斗河分別長15.86、7.6和7.19 km，地形西高東低，地面高程67～72.15 m。

根據(jù)洪水來源、量級、組合方式以及進退洪口情況，擬定3個計算方案，即發(fā)生20年一遇、50年一遇和100年一遇洪水時，小洪河干流及其主要支流洪溪河和吉斗河根據(jù)河道實際過流能力自然漫溢進洪，楊莊水文站洪峰流量超過650 m3/s時，楊莊泄洪閘控泄650 m3/s，水位達到50年一遇水位71.54 m時，泄洪閘敞泄(最大1 500 m3/s)，多余洪水滯于蓄滯洪區(qū)內，計算方案見表1，框架示意見圖1。

表1 楊莊蓄滯洪區(qū)計算方案表Tab.1 Calculating schemes table for YangZhuang-detention basin

圖1 楊莊蓄滯洪區(qū)計算方案框架示意圖Fig.1 The computational schemes for YangZhuang

2.2 模型構建及參數(shù)選取和驗證

2.2.1 一維模型

小洪河干流橫貫楊莊蓄滯洪區(qū)，其上楊莊泄洪閘控制退洪，構建小洪河MIKE11模型，模擬退洪過程。根據(jù)模型需要及河道走勢，取小洪河楊莊泄洪閘至其下游2.4 km河道作為一維退水河道模型的計算范圍，以流量控制方式設置可控建筑物模擬泄洪閘調度工況；斷面數(shù)據(jù)采用小洪河近期規(guī)劃治理并結合實地查勘整理；模型上邊界為固定水位邊界(暫定邊界，耦合后自動轉化為一二維模型耦合處即時水位)，下邊界為河道末端的水位流量關系，邊界條件見表2；模擬時間步長取20 s。

2.2.2 二維模型

蓄滯洪區(qū)面積較小、無連續(xù)堤防分割，計算區(qū)按整體計算；

表2 小洪河退水一維模型邊界條件Tab.2 Boundary condition of Xiao-Hong-he MIKE11 model for the flood fall

區(qū)內DEM密度不大于50 m，單一網(wǎng)格面積不大于1 200 m2，邊長不大于30 m，利用第三方軟件(Bentley公司的GEOPAK Civil)對區(qū)內主要河溝進行三維開挖，將開挖數(shù)據(jù)導入模型插值，河溝處網(wǎng)格加密到10 m左右，采用三角形網(wǎng)格剖分，網(wǎng)格數(shù)11.5萬個。

上游開邊界為小洪河干流、洪溪河和吉斗河入流洪水過程，時間跨度為18 d；下游開邊界為小洪河干流下游出流，為固定水位(暫定邊界，耦合后自動轉化為一二維模型耦合處即時水位)；其余為系統(tǒng)默認的陸地邊界，計算區(qū)外圍均為閉邊界；阻水地物考慮可能對洪水行進阻滯的0.5 m以上線狀地物，區(qū)內僅有合水村圍堤，堤頂高程72.0 m，圍堤長2.66 km，堤頂寬3 m，邊坡1∶2.5，二維模型信息見表3。

表3 二維模型信息統(tǒng)計表Tab.3 Information statistics table of Mike21 model

二維模型計算時間根據(jù)來水和退洪全過程確定為17～18 d，計算步長根據(jù)精度要求和模型收斂需要采用0.01～2 s可變步長、輸入邊界與一維模型銜接采用20 s，輸出時間步長根據(jù)結果要求設定為1 h；干濕邊界綜合考慮模型計算精度和穩(wěn)定性確定，干水深、浸沒水深和濕水深分別取0.002、0.05和0.1 m。

2.2.3 一二維耦合模型

結合調度運用方式，選取標準連接方式，在二維模型泄洪閘附近建立一個耦合點與一維模型起點進行耦合，局部地形及耦合情況見圖2。

圖2 楊莊蓄滯洪區(qū)MIKE FLOOD模型局部圖Fig.2 The section view of MIKE FLOOD for YangZhuang

主要參數(shù)有動量因子、水深調整和指數(shù)平滑因子。在不影響一二維流量和水位交互情況下，為了模型穩(wěn)定，忽略能量傳遞，動量因子取0；在天然地形或渠道上設置的連接，流過連接的水流會基于MIKE21單元水深分配水流到一維單元中，水深調整為激活；綜合模型穩(wěn)定性和水位傳遞精度要求：指數(shù)平滑因子選中間值0.5，耦合點處一二維模型水位相差小于0.05 m，滿足精度要求。

2.2.4 模型參數(shù)驗證

模型涉及的主要水力參數(shù)即一二維模型的糙率值，河道糙率選用其近期治理設計值；二維模型糙率依據(jù)《水力學計算手冊》并結合土地利用取值，如表4。

表4 模型糙率值表Tab.4 The roughness value of model

蓄滯洪區(qū)歷史啟用均無實測進洪過程，且隨著社會經(jīng)濟發(fā)展，區(qū)域內建筑物、構筑物等變化較大，模型相關參數(shù)難以用實測數(shù)據(jù)率定，故結合理論分析與實際經(jīng)驗類比選取，與歷史運用總體情況對比驗證。

楊莊滯洪工程在1998年建成后，2000年汛期連續(xù)兩次滯洪，區(qū)內經(jīng)濟損失巨大，后期與老王坡滯洪區(qū)聯(lián)合運用后未再進洪，故驗證洪水取2000年的2場洪水。6月24日至28日大到暴雨首次滯洪，簡稱“6.26”洪水；7月13日至15日再次遭遇大到暴雨，簡稱“7.13”洪水，經(jīng)比對分析洪水重現(xiàn)期分別相當于20年和50年一遇，故選取方案1和2模擬成果與實際運用情況對比，結果見表5。相同滯洪水位下滯洪量和淹沒面積相對誤差分別為-5.8%、10.5%和-9.7%、4.7%，存在誤差主要原因是本次模擬和2000年實測區(qū)域內前期影響雨量和來水過程不完全相同，致使進洪過程不完全一樣；同時由于經(jīng)濟社會的發(fā)展，區(qū)域內的構筑物以及下墊面已經(jīng)發(fā)生了一定的變化，致使洪水演進過程亦不完全一樣。但總的洪水模擬情況基本一致，且相對誤差較小，因此參數(shù)選用較為合理。

表5 與2000年運用情況對比表Tab.5 Comparative Table for results between simulation and 2000 observation

表6 楊莊蓄滯洪區(qū)洪水模擬計算結果Tab.6 The results of flood simulation for YangZhuang-detention basin

2.3 模擬結果

參數(shù)驗證后，模擬楊莊蓄滯洪區(qū)擬定的3個計算方案，各方案洪水演進的最大淹沒面積、水深、歷時、滯洪水位以及滯蓄洪量等結果見表6；每天零點滯洪區(qū)退洪流量見表7，退洪過程見圖3。

表7 楊莊泄洪閘下泄即滯洪區(qū)的退洪過程Tab.7 The flood discharge of YangZhuang-sluice

圖3 3種方案楊莊泄洪閘下泄流量過程線Fig.3 The hydrograph of YangZhuang-sluice for three schemes

3 成果合理性分析

3.1 淹沒演進合理性分析

從地形走勢與洪水演進趨勢匹配度、同一方案洪水演進過程和不同方案淹沒情況對比分析成果合理性。

地形云圖的高度走勢與洪水淹沒演進趨勢對比見圖4，二者完全匹配；同一量級洪水不同時刻洪水演進對比見圖5，變化趨勢符合實際洪水行進過程；不同量級洪水同一時刻演進對比見圖6，隨著洪水量級增大，淹沒范圍增大、水深加深，遵循一般規(guī)律。因此，楊莊蓄滯洪區(qū)洪水模擬模型洪水淹沒演進合理。

3.2 退洪合理性分析

(1)下泄過程合理性分析。對比模型計算成果(圖3和表7)和泄洪閘調度原則，泄洪閘下泄流量過程符合調度原則，退洪合理。

圖4 地形走勢與洪水演進趨勢對比圖Fig.4 Comparison diagram for topographic trend and flood routing trend

圖5 同一量級洪水不同時刻淹沒對比圖Fig.5 Comparison diagram for flood routing at different time in same recurrence period

圖6 不同量級洪水同一時刻淹沒對比圖Fig.6 Comparison diagram for flood routing in the same time for different recurrence period

(2)水量平衡合理性分析。各方案計算得入流總量、出流總量及滯留量成果見表8。考慮到一、二維模型本身無水量損失，累計各時刻入流總量為出流總量和滯留量之和，即滿足總體水量平衡的規(guī)律。其中滯留量是由于局部洼地在退洪時地勢較低積水無法排出導致，還有部分滯留量為河道占用水量。因此，從水量平衡角度分析該計算成果合理。

3.3 與原設計方案比較分析

模擬成果與楊莊滯洪區(qū)樞紐主體工程初步設計中確定的調度原則及設計成果對比，見表9。

表9 本次計算成果與原設計成果對比表Tab.9 Comparative table for simulated results and original design results

由表9可知，計算成果與原設計成果有一定差別。本次模擬區(qū)內水位動態(tài)變化，最大滯洪水位在不同位置略有差別，選取壩前水位略有偏低，同時受區(qū)內地物以及下墊面變化等影響，模擬值與設計值必然存在一定的差別，但水位絕對差小于0.2 m，在允許誤差范圍內，成果較為合理。

4 結 論

(1)MIKE系列模型構建蓄滯洪區(qū)洪水模擬模型，須結合蓄滯洪區(qū)實際情況，構建合適的一二維模型并聯(lián)合使用[13,14]。對于多條河流自然漫溢進洪、控制閘控制泄洪的蓄滯洪區(qū)，宜通過局部加密網(wǎng)格將河溝概化到MIKE21模型中，模擬洪水進洪及演進過程；實際調度運用方式通過MIKE11模型可控建筑物的設置，模擬退洪過程。MIKE FLOOD一二維耦合模型中采用標準連接方式可使一二維模型實時交換數(shù)據(jù)、穩(wěn)定計算，并能夠較準確、靈活地模擬蓄滯洪區(qū)的調度運用方式。

(2)蓄滯洪區(qū)面積相對較小，二維水動力學模型構建中網(wǎng)格剖分應結合區(qū)內現(xiàn)有溝渠的尺寸確定網(wǎng)格尺寸的大小，在溝渠位置適當加密才能夠保證洪水模擬的精度。

(3)為了更好地發(fā)揮數(shù)學模型的決策支持作用，需要進一步研究MIKE SHE模塊建立降雨徑流模型[15,16]，并與現(xiàn)有模型進行耦合，從而直接根據(jù)降雨推測洪水，實現(xiàn)在線實時洪水預報與洪水演進同時模擬計算；同時需進一步研究蓄滯洪區(qū)分蓄洪水的資源化利用[17-20]，區(qū)內排水順暢保證人民群眾生命財產安全的前提下，適時適地研究論證分蓄洪水的蓄水工程建設，用于枯水期灌溉、生態(tài)等用水需求，推動水資源高效利用。

□