基于MIKE 11模型的鵝頭分汊型河道水面線及分流比計算

顧 盼，張洪霞，耿文杰，連雷雷，左 建

（長江水利委員會水文局漢江水文水資源勘測局，湖北 襄陽 441000）

1 研究背景

丹江口水庫樞紐蓄水運用后，漢江中下游水沙條件及河勢發生了根本變化，河道由亞分汊淤積型向穩定分汊侵蝕型演變，河床沖刷，主支易位，魚梁洲灘尾下延，遇大水年份，深泓線變遷，頂沖點遷移，但由于人為活動的干預，修筑河堤，實施河道護岸工程、建造河勢控制工程，致使河道外形已不大可能發生大的變化，唯江心河灘及汊道進行沖淤的交替。另外，建庫后襄崔河段來水長期較為平穩，水流的流向、流路變化較小，左汊曾經經歷過長期清水沖刷，河床沖刷達到一定程度后很難繼續發展，且2009年10月崔家營航電樞紐蓄水運用前，對魚梁洲洲體進行了整治，在洲灘上修建了防洪子堤，開挖了南北汊河道及魚梁洲洲灘，子堤的高程為64.3 m。因此，漢江魚梁洲河段未來河勢不會發生大的變化。

漢江魚梁洲河段屬典型的鵝頭分汊道型河道，在丹江口水庫和崔家營水庫調度的雙重影響下，河段水流狀態較復雜。本文采用2015年最新實測固定斷面資料，對不同來水條件下魚梁洲河段水面線及左、右汊分流比進行了計算和研究。研究結果對于此河段分汊河道設計修建涉水工程可提供一定的參考。

2 河段概況

襄陽～崔家營河段，河道穿越城區將襄陽市一分為二，右岸為襄城區，左岸為樊城區。北河段為彎曲的河道，形狀為開口向西南方向的“U”形。漢江河床在由西向東流經襄陽市區中部呈Y字型分流，主河道開始轉向北偏東30°左右的流向，在漢江鐵路橋以下，河流由東折向東南，河床寬1000 m～6000 m。在其分流處下游約3 km小清河從自西北方向匯入，下游約6 km與唐白河匯合后轉向南偏東流向而后逐漸轉向南，在兩河口下游12 km處與分流河道匯合。分流河道流向東偏南方向后逐漸轉為南偏東、最終向南與主河道匯合。由于小清河、唐白河自西北方向匯入，形成回流，經歷次洪水的沖刷、淤積、切割，演變成現在的魚梁洲。魚梁洲洲灘高程一般在62.0 m～66.5 m，枯水期其東西寬約5.3 km，南北長約10.6 km，全洲面積約為27 km2，洲頭距襄陽水文站約2.5 km。

20世紀60年代以前左汊為主汊，20世紀60年代末丹江口水庫蓄水后清水下泄，河床下切，主流歸槽，70年代以后，右汊逐漸發展為主汊。襄陽水文站斷面1974年至1992年實測河床平均沖深達1.79 m。目前該河段因水流沖走床面泥沙中的較細顆粒，使床面的泥沙組成粗化，沖刷基本平穩，河勢相對穩定。魚梁洲是河段內最大江心洲，崔家營水電樞紐建成蓄水后，魚梁洲的可開發面積約13 km2。現階段，魚梁洲的開發利用應堅持“不礙洪、穩河勢、保民生、促發展”的原則，在崔家營航電樞紐正常蓄水位62.73 m（黃海）以上的部分洲灘進行適度利用。

3 MIKE 11模型

3.1 模型原理

MIKE 11模型由丹麥水力研究所（DHI）開發，包括HD模塊、RR模塊、DA模塊、FF模塊、ST模塊、AD模塊和SO模塊。MIKE 11HD作為MIKE 11模型建立過程中的核心模塊，可以對一維地表明渠流的水動力進行動態模擬，適用于山溪性、平原性河流的水力計算，應用范圍包括單一河道河流、汊狀河網及環狀河網。

模型的基本方程為圣維南方程組，其連續方程、動量守恒方程分別為：

式中：A為河道過水面積；Q為流量；u為側向流在河道方向的流速；t為時間；x為沿水流方向的水平坐標；q為河道的側向流量；α為動量修正系數；g為重力加速度；y為水位；Sf為糙率系數。

MIKE 11模型采用Abbott-Ionesco[3]六點中心隱式差分格式對圣維南方程進行數值離散化，并應用“追趕法”求解差分方程。由于Abbott六點中心隱式差分算法為無條件穩定方法，具有較高的計算穩定性，誤差率＜2%，且計算速度較快，一般比常規四點差分格式快5～6倍。

3.2 模型的構建

MIKE 11軟件主要用于河口、河流、灌溉系統和其他內陸水域的水文、水力、水質和泥沙傳輸模擬，在防汛、洪水預報、水資源水量、水質管理、水利工程規劃、設計、論證等領域得到廣泛應用。

基于MIKE 11系統的模型構建主要包括河網概化、斷面文件處理、邊界與初始條件設置、模型參數率定和驗證等。MIKE 11 HD包含河網文件、斷面數據、邊界條件和模型參數四個文件，前三個為建模基礎數據資料，模型參數文件為需率定的參數，根據模型原理可知，率定的參數為糙率系數。MIKE 11 HD的模型結構，見圖1。

圖1 MIKE 11 HD模型結構

3.2.1 模型輸入

（1）河網及實測斷面

本次計算以襄崔河段1∶5000實測河道地形圖為基礎資料，生成計算區域的河網文件，建立左、右汊河道的拓撲結構，進行水系概化；采用河段2015年實測的11個固定斷面為基本控制斷面，示意圖如圖2。

圖2 襄陽～崔家營河段河道及實測斷面布置及洪痕位置示意圖

（2）邊界與初始條件

根據河道分汊情況，確定模型的邊界條件，選取襄陽水文站上游約1 km處HZ49-1為上游邊界輸入斷面，崔家營電廠大壩壩址處HZ55-1為下游邊界出口斷面；上邊界輸入流量，下邊界輸入對應控制水位。

3.2.2 模型率定

以2017年10月在河段5個位置處的實測水位以及左、右汊分流比為控制目標，對模型中反映河床糙率的糙率系數進行率定。經過反復率定，漢江右汊糙率為0.025，左汊糙率為0.021。模擬水位與實測斷面水位相差不超過0.05 m，見表1；模型率定水面線與實測洪水位見圖3；模擬的左汊流量1785 m3/s、右汊流量6225 m3/s，占總流量比例為22.3%和77.7%，實測左汊流量1800 m3/s、右汊流量6210 m3/s，占總流量比例為22.5%和77.5%，見表2。

表1 模型率定水位與實測水位

圖3 模型率定水面線與實測洪水位

表2 模型率定分流比與實測分流比

3.2.3 模型的驗證

模型的驗證采用2011年9月漢江洪水，襄陽站流量13000 m3/s，崔家營水位62.86 m，用實測水位與分流比對模型進行驗證，其計算成果與實測成果的比較，見表3、表4和圖4。結果表明，模擬水位與實測水位偏差較小，相差不超過0.03 m，模擬的左汊流量、右汊流量占總流量比例為22.8%、77.2%，實測左汊流量、右汊流量站總流量比例為23.0%：77.0%。可見率定的MIKE 11模型具有較高精度，可用于本河段中高型典型洪水水面線及分流比的推求運用。

表3 模型率定水位與實測水位

圖4 模型驗證水面線與實測洪水位

表4 模型驗證分流比與實測分流比

4 典型洪水水面線及分流比推算

新集水電站位于本河段上游約15 km，區間無其它支流，因此根據新集水電站設計洪水流量成果，采用率定的糙率系數計算襄陽～崔家營河段5年、20年、50年、100年和漢江1964年洪水的水面及漢江左、右汊分流比。各流量條件下設計洪水水面線，見表 5～6，分流比見表 7。

表5 不同來水條件下漢江右汊水面線（干流+右汊）

表6 不同來水條件下漢江左汊水面線（干流+左汊）

表7 不同來水條件下漢江魚梁洲左汊分流比

5 結論

（1）各級不同來水條件下，漢江魚梁洲左、右汊分流比比較穩定，左汊分流量占總流量的23%左右，右汊分流量占77%。因此，當漢江魚梁洲左、右汊河道設計建設涉水工程時，設計洪水可采用此分流比推求兩汊相應的流量。

（2）當漢江干流總來水流量增大時，漢江魚梁洲左汊的分流比在一定程度上將隨之增大，但增加幅度有限，5年～100年一遇設計洪水條件下，左汊分流比從22.6%增加至23.3%，僅增加0.7%。