基于MIKE 11法對現狀堤防防洪能力研究

劉躍猛

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

基于MIKE 11法對現狀堤防防洪能力研究

劉躍猛

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

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

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

1 研究情況

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

2 模型建立

2.1 河網文件

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

圖1 基本資料與MIKE 11文件關系圖

2.2 斷面文件

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

圖3 海城河斷面文件

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

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

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

2.4 參數文件

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

2.5 模擬文件

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

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

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

3)模擬菜單。時間步長的確定經常要通過反復試算調整，與河床地形與邊界條件密切相關，我們選擇固定時間步長為30sec。

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

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

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

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

2.6 模型驗證

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

表1 實測水位與MIKE11模擬水位結果表

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

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

3 模型推求洪水水面線

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

3.1 計算水面線

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

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

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

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

3.2 漫溢險情分析

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

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

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

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

3.3 橋梁的壅水分析

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

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

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

3.4 典型河段流場分析

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

圖8 15+605斷面流場分布

4 結 論

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

在對河道一維洪水數值模擬中，也充分證明mike11模型在該區域洪水演進模擬中具有切實的可行性，為規劃斷面的一維模擬提供了依據。可以看出，一維模型建模簡單，可以隨時根據河道斷面更新數據對模型進行修改，同時模擬時間較短[11]。

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

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

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

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

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

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

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

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

[9]包宇.大凌河流域朝陽區段水環境恢復治理對策[J].水土保持應用技術，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