劉躍猛

(海城市河務(wù)管理處，遼寧 海城 114200)

基于MIKE 11法對現(xiàn)狀堤防防洪能力研究

劉躍猛

(海城市河務(wù)管理處，遼寧 海城 114200)

MIKE11以求解圣維南方程組作為理論基礎(chǔ)，包含了眾多的水文模擬方法，是動態(tài)模擬河流和水道水力的世界級標準，具有無限的河流模擬能力，它是國際多個國家模擬河流使用驗證的標準工具，文章對現(xiàn)狀河堤防防洪能力分析的方式、方法進行論述，研究結(jié)果將為MIKE11在河道中現(xiàn)狀防洪等問題提供參考借鑒。

MIKE11；防洪能力；分析；河流模擬；工具

1 研究情況

Mike一維洪水演進主要是應(yīng)用海城河河道斷面，建立完整的河網(wǎng)模型，輸入水文資料，確定邊界條件，進行海城河水面線的推求。模型結(jié)構(gòu)見圖1。以海城河東三道橋至高鐵橋段，以此為例進行研究討論。

2 模型建立

2.1 河網(wǎng)文件

生成河網(wǎng)文件有兩種方法，①利用DEM數(shù)據(jù)；②利用BMP地圖進行數(shù)字化。海城河河網(wǎng)文件如圖2。

圖1 基本資料與MIKE 11文件關(guān)系圖

2.2 斷面文件

將某年實測斷面數(shù)據(jù)(河道各斷面x，z數(shù)據(jù))轉(zhuǎn)換為MIKE 11需要的文檔格式，導入斷面文件，見圖3。

圖3 海城河斷面文件

其中各斷面里程從CAD圖中計算得出，起始點東三道明水橋，樁號為22+527。海城河斷面多為復(fù)式斷面，呈階梯型形，主槽多位于河道的中心[2]。本次計算中，糙率根據(jù)實地勘測河道地形及組成確定。自東三道至東外環(huán)橋糙率系數(shù)n主槽選定為0.03，灘地為0.04，東外環(huán)橋至鐵路橋之間，主槽糙率選定為0.025，灘地糙率定為0.065；鐵路橋至高鐵橋，主槽糙率選定為0.03，灘地糙率確定為0.065。

2.3 邊界條件和時間序列文件

本次模擬區(qū)域下游高鐵橋處，距離河口較遠，該處水位隨著流量的變化而變化，且無實時水文監(jiān)測資料[3]，下游邊界不宜設(shè)為Z(t)、Q(t)。

2.4 參數(shù)文件

本次計算中，糙率根據(jù)高分辨率衛(wèi)星圖片和多光譜實地勘測河道地形確定[4]。自東三道至東外環(huán)橋主槽糙率系數(shù)n選定為0.035，邊灘糙率定為0.04；東外環(huán)至鐵路橋之間河道經(jīng)過幾次整修已形成較為規(guī)整的堤防，糙率選定為0.025，灘地糙率定為0.065；鐵路橋至高鐵橋，該段內(nèi)邊灘多為樹木和農(nóng)作物，主槽糙率選定為0.03，邊灘糙率系數(shù)選定為0.065[5]。

2.5 模擬文件

模擬文件共有5個菜單，下面就逐一介紹各個菜單的設(shè)定。

1)模型菜單。模型類型選擇水動力模型，模擬方式選擇非恒定流。

2)輸入菜單。在這里可以引入剛才生成的所有文件：河網(wǎng)文件、斷面文件、邊界文件、HD參數(shù)文件，并且可以隨時編輯這些文件。

3)模擬菜單。時間步長的確定經(jīng)常要通過反復(fù)試算調(diào)整，與河床地形與邊界條件密切相關(guān)，我們選擇固定時間步長為30sec。

模擬的起始時刻和結(jié)束時刻要在邊界條件設(shè)置的時間之內(nèi)，不然模型無法起算。接下來選擇初始條件設(shè)定。本次計算中。糙率在斷面文件和參數(shù)文件中均進行了設(shè)立[6]。穩(wěn)態(tài)和參數(shù)文件的選擇并不能改變模型模擬的結(jié)果。計算中，我們采用穩(wěn)態(tài)(steady state)進行模擬。

4)計算結(jié)果菜單。定義輸出結(jié)果文件名和保存頻率。在模擬菜單中計算時間步長為30sec，我們不需要這么密的計算值，0.5h一個結(jié)果已經(jīng)足夠，那么就可以定義保存頻率為6個時間步長，即計算6個時間步長保存一次結(jié)果。這樣可以減小結(jié)果文件大小。

5)開始菜單。進入Start菜單，準備開始計算。在驗證狀態(tài)框內(nèi)都是綠燈時，就可以按Start鍵開始計算了；如有紅燈，那么在下面的驗證信息框內(nèi)就會出現(xiàn)相應(yīng)的出錯信息，按照提示進行修改，直到出現(xiàn)綠燈。

計算運行時，有時會出現(xiàn)警告信息，但這些問題還不至于影響到模型的運行。若是出現(xiàn)錯誤信息，則會停止運行。按照錯誤信息的提示修改各文件，或是調(diào)節(jié)某個參數(shù)，直到模型成功運行到100%。這樣模型就建立成功，只需要進行參數(shù)的率定驗證了[7]。

2.6 模型驗證

采用2002年海城水文站實測水文資料對模型進行檢驗。對比海城水文站，模型計算水位值與實測水位值見表1及圖4。

表1 實測水位與MIKE11模擬水位結(jié)果表

圖4 模擬水位與實測水位對比

分析表1及圖4發(fā)現(xiàn)，模型模擬結(jié)果與實測結(jié)果吻合較好，兩者最大差值為0.41m發(fā)生在洪水漲水期，在洪峰期，模擬值與實測值基本一致。計算實測-模擬相關(guān)系數(shù)，達到0.994，證明模型具有一定的準確性[8]。

3 模型推求洪水水面線

采用2002年作為典型洪水過程線，以水位-流量作為下游邊界條件根據(jù)，利用mike11建立海城站一維洪水演進模型。

3.1 計算水面線

采用2011年實測斷面，2002年典型年設(shè)計洪水為計算條件，模擬計算50a遇洪水水面線。在不靠考慮風浪影響的情況下，各頻率設(shè)計洪水，模型計算水面線與左右岸堤防關(guān)系，如圖5、6所示，海城河典型地段水面線計算成果見表2。

圖5 50a一遇洪水水面線與左岸堤防對比

圖6 50a一遇洪水水面線與右岸堤防對比

樁號位置模擬洪水水面線現(xiàn)狀堤頂高程P=2%左岸右岸最深河底高程橋梁底部高程22+527東三道橋374385441316914+728東外環(huán)橋28428729233145+583二臺子大橋2022522252310+000高速鐵路橋16318751918858

3.2 漫溢險情分析

根據(jù)所繪制的河道行洪水面線，結(jié)合各個斷面的左右岸堤防，可以確定河道的堤防淹沒情況。遭遇50a一遇洪水時，15+605斷面至11+611斷面，除11+864斷面右岸沒有被淹沒外，其余斷面左右岸均有不同程度的漫堤；22+369、21+279、18+832、16+476、10+113、5+677和4+123斷面的右岸也出現(xiàn)漫堤情況。遭遇20a一遇洪水時，14+663斷面至11+611斷面右岸漫堤，其中14+568和13+539斷面輕微淹沒；左岸漫堤的斷面有18+832、12+539、14+663和11+968斷面，13+539、5+677和4+123斷面輕微淹沒。發(fā)生10a一遇洪水時，14+663斷面和12+539至11+611斷面右岸有漫堤情況；左岸淹沒的斷面有18+832、12+539和11+968斷面。

建立洪水演進模型對規(guī)劃堤防進行防洪能力評價分析。在考慮風浪情況下，城區(qū)段分別采用100a一遇和50a一遇設(shè)計洪水水面線，與堤防進行對比。對比結(jié)果顯示原規(guī)劃堤防共有12處斷面是由于左岸堤防過低發(fā)生漫溢情況；有1處斷面是由于右岸堤防過低，發(fā)生漫溢險情、樁號為12+347和11+968斷面兩岸均發(fā)生了漫堤。

由于項目段主要位于在市區(qū)內(nèi)，兩岸用地都較為緊張，河道拓寬成本較高，難以實現(xiàn)。在河道重新規(guī)劃過程中，主要采用加筑堤防，拓寬、清淤主槽增大河道過流能力等方式達到防洪要求。

綜上所述，12+539至11+611斷面右岸為主要漫堤線，應(yīng)采取加高堤防或拓寬河道等規(guī)劃措施。

3.3 橋梁的壅水分析

50a一遇橋梁壅水分析見圖7。

圖7 50a一遇洪水橋梁壅水分析

圖7中藍線為在有橋梁情況下，50a一遇設(shè)計洪水演進水面線；粉線為為有橋梁情況下，50a一遇設(shè)計洪水演進水面線。由于橋梁的壅水作用，整體上水位模擬結(jié)果包含橋梁的模型略高于未包含橋梁的模型。

3.4 典型河段流場分析

河道走勢的變化，會導致水流方向的改變。而對于河流走勢發(fā)生較大的堤段，由于水流的沖刷作用，往往會導致漫堤、甚至潰堤的可能。利用MIKE Zero中的Result Viewer工具，可以將模擬時間段河道內(nèi)的流場變化呈現(xiàn)出來。以下示意圖分別是遭遇50年一遇的洪水時，典型斷面及河道轉(zhuǎn)彎處的流場分布情況[9]。15+605斷面流場分布見圖8。

圖8 15+605斷面流場分布

4 結(jié) 論

本文總結(jié)了采用一維數(shù)值模型對河道內(nèi)洪水演進進行數(shù)值模擬，并取得了較好的效果。研究思路初步涵蓋各種防洪影響情況，具有普遍借鑒意義，但一維模型難以考慮河道的局部特征，并且不能模擬河道流場變化規(guī)律，因此，往往需要二維模型進行對比計算[10]。

在對河道一維洪水數(shù)值模擬中，也充分證明mike11模型在該區(qū)域洪水演進模擬中具有切實的可行性，為規(guī)劃斷面的一維模擬提供了依據(jù)?？梢钥闯觯痪S模型建模簡單，可以隨時根據(jù)河道斷面更新數(shù)據(jù)對模型進行修改，同時模擬時間較短[11]。

[1]胡朝陽, 劉建衛(wèi), 朱林.糙率改變對城市河道水位及堤防防洪能力的影響[J].水電能源科學, 2012(10):47-50.

[2]姜彪.基于洪水數(shù)值模擬的堤防安全評價與對策研究[D].大連：大連理工大學, 2010.

[3]蔣書偉, 武永新.基于MIKE11與HEC-RAS的南渡江防洪能力對比分析[J].中國農(nóng)村水利水電, 2014(02):46-49.

[4]馮金鵬.大洋河岫巖城市段河道防洪能力評估研究[J].東北水利水電, 2016, 34(06)：26-28.

[5]王小丹.定向鉆穿越河道防洪評價的要點[J].水利技術(shù)監(jiān)督,2015,23(04):26-28.

[6]趙希嶺.防洪評價報告編制中有關(guān)問題的探討[J].水利技術(shù)監(jiān)督,2011,19(02):34-35.

[7]才慶欣.南票區(qū)女兒河項目區(qū)水土保持工程治理效益探析[J].水利規(guī)劃與設(shè)計,2015(05):108-110.

[8]李滿剛.城市小匯水區(qū)域設(shè)計洪水計算方法應(yīng)用研究[J].水利規(guī)劃與設(shè)計,2012(02):11-14.

[9]包宇.大凌河流域朝陽區(qū)段水環(huán)境恢復(fù)治理對策[J].水土保持應(yīng)用技術(shù)，2011(06)：39-40.

[10]袁崢.渭河下游干流堤防防洪能力分析[J].地下水，2005(09)：32-36

[11]梁在潮，李泰來.江河堤防防洪能力的風險分析[J].長江科學院院報，2001(04)：67-68.

1007-7596(2016)11-0025-04

2016-10-20

劉躍猛(1981-)，男，遼寧海城人，助理工程師，研究方向為河道治理與管理、水利工程施工、水利工程管理等。

TV871

B