基于MIKE FLOOD模型的湟水河上游洪水風險評估

韓 嶺，蓋永崗，劉 楊，張瑞海

(黃河勘測規(guī)劃設計有限公司，鄭州 450003)

0 引 言

我國是世界上洪水災害多發(fā)的國家之一，近年來，受全球氣候變化的影響，暴雨極端天氣頻發(fā)，對于山區(qū)性河流，這種突發(fā)暴雨所引起的洪水對其影響更為顯著，對河流兩岸居民的生產生活造成了巨大影響，加之部分河流的防洪工程建設滯后，防洪應急管理能力不足，出現(xiàn)了嚴重的洪水災害。隨著防洪理念由控制洪水向洪水管理的轉變，洪水數(shù)值模擬已成為洪水風險圖繪制、洪水風險分析、洪水損失評估的一個重要手段，是洪水管理的科學依據[1-3]。

目前比較流行的水動力數(shù)值模擬軟件有荷蘭的Delft 3D模型[4,5]、英國的InfoWorks RS模型[6]、丹麥的MIKE模型[7,8]及其他自主開發(fā)的軟件模型[9-11]，雖然這些軟件在功能實現(xiàn)、計算算法和建模方法上各有不同，但都得到了廣泛的使用和工程驗證，具有很高的可信性。其中MIKE模型因其強大的前、后處理功能，是目前世界范圍內應用最為廣泛的水力學模型之一[12-18]。針對湟水河上游段的河道特征及地形地貌，本文擬采用MIKE FLOOD 模型建立該河段的一二維耦合水動力學模型，對不同量級的洪水方案進行模擬，獲得不同水文要素的時空分布及洪水淹沒信息。

1 研究區(qū)概況

湟水流域位于青藏高原與黃土高原的過渡地帶，地勢西北高、東南低，境內高山、丘陵交錯分布，高差懸殊，區(qū)域內地形多變，有高山、中山、黃土覆蓋的低山丘陵和河谷盆地，古老基地局部隆起形成峽谷，分隔了中生代斷陷盆地，各盆地呈串珠展布。湟水流域夏秋多暴雨，沿途溝道發(fā)育，山區(qū)地面坡降陡，匯流時間短，易在短時間內形成洪峰，引起河道水位暴漲，使人民生命財產遭受嚴重損失。根據湟水干流的地形地貌特征，湟水干流分為上、中、下游三段。西寧以上為上游，河段長174 km左右，面積9 022 km2，落差1 976 m，平均比降1.14%；西寧至民和為中游，民和以下為下游。

本次研究區(qū)域為湟水上游的東大灘水庫溢洪道壩下至湟中縣多巴鎮(zhèn)黑嘴橋河段，河道長約56.2 km，研究范圍40.8 km2，河道細長彎曲屬于典型的山區(qū)性河道。湟水河上游河道干流主要流經湟源縣和湟中縣，其中湟源縣境內河道長約41.9 km，已建堤防(護岸)長度34.68km，其中左岸14.4 km，右岸20.28 km，防洪標準為50年一遇；湟中縣境內河道長約14.26 km，已建堤防長度7.7 km，其中左岸4.03 km，右岸3.67 km，防洪標準為30年一遇，防洪工程建設情況見表1。

表1 湟水河上游段防洪工程統(tǒng)計表Tab.1 The upper reaches of Huangshui River flood control projects statistical

2 模型原理

2.1 一維模型原理

河道一維模型采用MIKE11水動力模塊，模型基本方程為圣維南方程組，包括質量守恒方程(連續(xù)方程)和動量方程：

(2)

式中：q為旁側流量；Q為總流量；s為距離坐標；V為斷面平均流速；h為水深；A為過水斷面面積；i為渠底坡降；C為謝才系數(shù)。

Mike11模型采用6點Abbott～Ionescu有限差分格式對圣維南方程組求解，該模型基于反映有關物理定律的微分方程組，并對微分方程組進行線性化的有限差分格式，以及求解線性方程組的算法。

2.2 二維模型原理

河道外二維區(qū)域采用MIKE21 FM (非結構化網格)水動力模塊建模，其原理基于二維不可壓縮流體雷諾平均應力方程，服從布辛涅斯克假設和靜水壓力假設：

h=η+d

(3)

連續(xù)方程：

(4)

動量方程：

(5)

(6)

(7)

2.3 模型耦合

河道一維模型(MIKE11)和河道外二維模型(MIKE21)通過MIKE FLOOD模塊連接在一起，成為動態(tài)耦合地模型系統(tǒng)，使模型既繼承了一維和二維模型的優(yōu)點，又避免了采用單一模型時遇到的網格精度和準確性方面的問題。本文通過側向連接方式將一維模型連接到二維模型中，在計算時，一二維模型同時計算，在連接處進行動量傳遞，可以較好地模擬洪水匯入河道的匯流過程和洪水從河道漫溢到堤防外的淹沒過程。

3 模型建立與驗證

3.1 一維模型構建

根據研究區(qū)域的基礎地理信息數(shù)據提取水系圖層文件，構建湟水河上游段河道一維模型(MIKE11)的河網文件。為保證建立模型的合理性，在一維水動力學模型中，堤防(護岸)在河道斷面文件中以堤岸的形式體現(xiàn)在模型中；引橋路基則根據距離該處最近的河道上下游斷面形態(tài)，通過適當修正以河道斷面的形式加入模型中；攔河閘壩通過添加河道建筑物中的溢流堰加入模型中。

根據河道斷面設計資料以及相關實測資料，湟水河上游河道平均300 m設置一個斷面，共設置河道斷面144個，其中河道斷面116處，引橋路基斷面28處；概化攔河閘壩13座，一維水動力學河網概化見圖1。

圖1 一維水動力學模型河網概化圖Fig.1 River network generation of 1D model

3.2 二維模型構建

湟水河上游段河道外采用二維水動力學模型(MIKE21)進行搭建。湟水河上游二維水動力學模型計算區(qū)域面積為40.8 km2。將計算范圍邊界線及網格剖分的分區(qū)線圖層文件導入MIKE21 FM模型，為滿足模型計算精度，分區(qū)設定網格剖分最大面積、最大角度等參數(shù)，然后生成網格。二維模型網格節(jié)點數(shù)37 533個，網格數(shù)為60 435個，平均網格面積1 587 m2，單網格最大面積為10 000 m2，最小面積46 m2，網格剖分結果見圖2。

圖2 二維水動力學模型網格剖分圖Fig.2 Mesh generation of 2D model

對各計算分區(qū)進行網格剖分后，采用1∶2 000比例尺的DEM高程數(shù)據，進行網格高程插值，得到二維模型的地形云圖。

在湟水河上游二維水動力學模型建立過程中，高速、國道等阻水線狀地物均以堤防建筑物的形式加入模型中，不考慮風、浪、潮汐、溫度、鹽度等因素的影響。

3.3 邊界條件

研究區(qū)域內的MIKE11一維水動力學模型上游開邊界為東大灘水庫尾水斷面的設計流量過程線，下游開邊界為湟水河黑嘴橋河道大斷面的水位流量關系，見圖4。

圖3 黑嘴橋斷面水位流量關系Fig.3 Relationship between water level and water level in the section of HEI Zui bridge

湟水干流兩岸支流眾多，上下游干支流的洪水常有遭遇的情況，設計洪峰沿程逐漸增加。為保證湟水河上游干流河段洪水分析的連續(xù)性，且各斷面均能達到設計洪峰流量數(shù)值，在洪水分析時，干支流交匯處均假設以支流匯合口上下游斷面同頻，各支流(區(qū)間)相應。根據湟水河干流上段洪水設計成果資料，分別在胡旦橋、大路橋、石崖莊、西納川入匯口處將區(qū)間設計洪水過程作為點源加入，使干流上下游洪水與設計洪水保持一致。東大灘尾水斷面和各點源洪峰流量，見表2。

表2 上邊界及點源洪峰流量 m3/s

研究區(qū)域內的MIKE21二維水動力學模型中除開邊界外，其余邊界為固定邊界。開邊界通常指計算區(qū)的出流和入流邊界。研究區(qū)域內河道外二維水動力模型的入流設置通過MIKE FLOOD耦合模型一二維側向連接實現(xiàn)，二維模型本身不再設置出入流邊界。

3.4 耦合參數(shù)

由于一維模型中的構筑物高程和二維模型中的網格地形存在偏差，且本次研究區(qū)域內的地形采用較為精細的1∶2 000的地形圖，為了更為準確的模擬洪水漫堤的情況，耦合模型中堰的底高程采用二維模型的單元高程。

3.5 模型參數(shù)率定

湟水河上游段河道一維水動力學模型率定選用2010年石崖莊站(湟源站)實測的水位進行率定。采用2010年7月7日場次洪水流量過程作為模型上游邊界條件，進行洪水演進模型計算，提取計算成果中石崖莊水文站觀測斷面所對應模型河道斷面的水位與流量結果文件，分析整理后與石崖莊水文站實測資料進行對比分析，分析結果見表3。

表3 石崖莊站斷面水位流量關系實測與計算對比表Tab.3 Relationship between water level and water level in the section of SHI Yazhuang

從表3中可以看出：模型計算水位與實測水位相差較小，基本一致。考慮到模型本身等誤差，可以認為湟水河上游段一維水動力學模型模擬效果良好，模型參數(shù)設置合理。

二維水動力學模型建模中，需要著重分析論證的是二維模型地形糙率。根據《洪水風險圖編制技術細則》[19]、《水力計算手冊》[20]等資料，對不同下墊面賦予不同的糙率值。為保證模型計算精度，根據研究區(qū)域內的土地利用類型為每個網格設定不同的糙率值，通過插值得到計算區(qū)域糙率值云圖，見圖4。研究區(qū)域內的主要土地利用類型包括耕地、居民地、灘地、林地、草地等，其分布比例見表4。

圖4 二維模型地面糙率值云圖Fig.4 Surface roughness map of two dimensional model

表4 土里利用類型分布情況表 %

3.6 模型驗證

根據洪水資料的收集情況，本次模型驗證采用相關設計報告中并已正在實施的湟水河八敗溝至石板溝河段的設計水面線，對河道模型進行參數(shù)驗證。

選取湟水河八敗溝至石板溝河段9個斷面進行模型驗證，通過對比分析，設計報告水位和本次計算水位相差比較小，可以認為模型搭建合理。水面線對比見表5。

表5 湟水河上游段一維河道模型水面線成果對比表Tab.5 Contrast of water line results of the one-dimensional channel model of the upper reaches of Huangshui River

4 計算結果與風險分析

根據洪水分析計算方案，研究區(qū)域內發(fā)生10年一遇洪水時，主要在下胡旦橋至大路橋河段及果米灘河段、石板溝村河段、扎麻隆村以下河段漫溢；發(fā)生50年一遇洪水時，主要在下胡旦橋至大路橋河段及石崖莊至多巴鎮(zhèn)河段漫溢；發(fā)生100年一遇洪水時，主要在下胡旦橋至大路橋河段及湟源縣城內的立達橋至東風橋段、石崖莊至黑嘴橋全段漫溢，且漫溢最大范圍都發(fā)生在洪峰到達河道斷面的時刻。

研究區(qū)域內10年、50年和100年一遇淹沒水深面積分別為0.47、2.26和2.46 km2，隨著洪水量級的增加，淹沒總面積逐漸增大。不同頻率洪水淹沒主要集中在0.05～2.5 m之間，淹沒面積占到總淹沒面積的95%以上，研究區(qū)域內的洪水淹沒風險情況見圖5。

圖5 湟水河上游段洪水淹沒風險圖Fig.5 Inundation risk map of the upper reaches of Huangshui River flood

結合目前研究區(qū)域內防洪工程的現(xiàn)狀，湟水河大路橋—八敗溝段堤防防洪標準為50年一遇，八敗溝—石板溝段堤防防洪標準為10年一遇，多巴鎮(zhèn)—黑嘴橋堤防防洪標準為30年一遇，由圖5可知，當發(fā)生10年一遇及50年一遇洪水時，湟源縣大路橋至湟源縣城段洪水不出槽；在遭遇超標準洪水(100年一遇洪水)時，在湟源縣城段局部洪水漫溢出河槽進入城區(qū)；在遭遇10年、50年、100年一遇量級洪水時，均在局部河段存在漫溢情況，并且隨量級增大其洪水要素也逐漸增大。

5 結 語

(1)當發(fā)生10年、50年一遇洪水時，湟水河上游段防洪工程基本上能防御相應級別的洪水，能較好地發(fā)揮防洪作用。特別是湟源縣城段河道防洪工程的修建對區(qū)域防洪安全起到了非常重要的作用。

(2)當發(fā)生100年一遇超標準洪水時，由于主槽的過流能力有限，在湟源縣城段局部洪水漫溢出河槽進入城區(qū)，但淹沒范圍相對很小且平均水深不足0.5 m；在申中鄉(xiāng)、多巴鎮(zhèn)河段洪水大面積漫溢，淹沒河道兩岸耕地和村莊，防洪能力略顯不足。

(3)建議根據洪水分析的結果，完善各地區(qū)防洪搶險預案，合理安排工程建設尤其湟水河上游段的申中鄉(xiāng)、多巴鎮(zhèn)河段的部分河道防洪工程建設。

(4)根據對計算結果的分析可知，MIKE FLOOD耦合模型合理的模擬了湟水河上游段遭遇不同頻率洪水時的洪水匯入河道的匯流過程和洪水從河道漫溢到堤防外的淹沒過程，為該地區(qū)的防洪救災提供了技術支撐，表明MIKE FLOOD耦合模型在山區(qū)性河流的風險評估中具有良好的適用性。

□

[1] 李 娜, 向立云, 程曉陶. 國外洪水風險圖制作比較及對我國洪水風險圖制作的建議[J].水利發(fā)展研究, 2005,12,(6):28-32.

[2] 馬建明, 許 靜, 朱云楓, 等. 國外洪水風險圖編制綜述[J].中國水利, 2005,10(17):29-31.

[3] 葉愛民, 劉曙光, 韓 超, 等. MIKE FLOOD耦合模型在杭嘉湖流域嘉興地區(qū)洪水風險圖編制工作中的應用[J].中國防汛抗旱, 2016,26(2):56-60.

[4] 左書華. Delft3D 在鰲江口外平陽咀海域流暢模擬中的應用[J].水文, 2007,27(6):55-58.

[5] 劉 成, 何 耘, 李行偉, 等. 上海市污水排放口污染物運動軌跡模擬[J].水利學報, 2003,13(4):114-118.

[6] 肖 魁. 基于InfoWorks RS 的鄱陽湖洲灘水情變化分析[D].武漢:長江科學院, 2011.

[7] 紀忠華, 胡勐乾, 王 璐, 等. 基于MIKE21模型的濱海核電廠廠址防洪評價研究[J].水利水電技術, 2016,47(1):132-137.

[8] 王崇浩,曹文洪,張世奇. 黃河口潮流與泥沙輸移過程的數(shù)據研究[J].水利學報, 2008,39(10):1 256-1 263.

[9] 宋利祥,周建中,王光謙, 等. 潰壩水流數(shù)值計算的非結構有限體積模型[J].水科學進展, 2011,22(3):378-381.

[10] 張秋霞, 王義成. 二維非恒定流洪水演進模擬模型開發(fā)及應用[J].水利水電技術, 2009,40(3):62-64,73.

[11] 魏文禮,沈永明,孫廣才,等. 二維潰壩洪水波演進的數(shù)值模擬[J].水利學報,2003,17(9):43-47.

[12] 欒 慕, 袁文秀, 劉 俊, 等. 基于SWMM-MIKE11耦合模型的桐廬縣內澇風險評估[J].水資源保護, 2016,32,(2): 57-61.

[13] 郭鳳清, 屈寒飛, 曾 輝, 等. 基于MIKE21FM模型的蓄洪區(qū)洪水演進數(shù)值模擬[J].水電能源科學,2013,31(5):34-37.

[14] 楊 洵, 梁國華, 周惠成. 基于MIKE11的太子河觀—葠河段水文水動力模型研究[J].水電能源科學, 2010,28(11):84-87.

[15] 袁雄燕, 徐德龍. 丹麥MIKE21 模型在橋渡壅水計算中的應用研究[J].人民長江, 2006,37(4):31-32,52.

[16] 修海峰, 吳聯(lián)志. 基于MIKE 11的平原河網洪水演進水動力研究[J].南水北調與水利科技, 2012,10(12):151-153.

[17] 李毓湘, 逄 勇. 珠江三角洲地區(qū)河網水動力學模型研究[J].水動力學研究與進展, 2011,16(2):143-155.

[18] 郭鵬程, 蔡 明, 閆大鵬. 基于MIKE21 模型的人工生態(tài)湖優(yōu)化設計[J].人民黃河, 2014,36(4):56-58.

[19] 水利部. 洪水風險圖編制技術細則(試行)[R]. 2013.

[20] 李 煒. 水力計算手冊[M]. 北京: 中國水利水電出版社, 2006.