李 航，丁 曼，徐貴新

(吉林省水利水電勘測(cè)設(shè)計(jì)研究院，長(zhǎng)春 130021)

飲馬河是吉林省第二松花江左岸最大支流，發(fā)源于四平市伊通縣地局子鄉(xiāng)老爺嶺東南，流經(jīng)伊通縣、磐石市、永吉市、雙陽(yáng)區(qū)、九臺(tái)區(qū)、德惠市、農(nóng)安縣和長(zhǎng)春市二道區(qū)等8個(gè)市縣(區(qū))，于農(nóng)安縣靠山鎮(zhèn)紅石壘屯匯入第二松花江，集水面積18 247 km2，其中石頭口門水庫(kù)以上集水面積4 944 km2。飲馬河河道長(zhǎng)度387.5 km，流域平均高程300 m左右。伊通河流域其上游屬丘陵崗地，下游地形波狀起伏，兩岸河谷較開闊，地勢(shì)較高部分大都開墾為旱田，河谷平地大部分開墾為水田。流域水系發(fā)育、支流較多，較大支流依次有伊通河、雙陽(yáng)河、霧開河、岔路河、三道溝河和小南河等[1]。

分析范圍為飲馬河石頭口門水庫(kù)至飲馬河河口，包括沿岸的防洪保護(hù)區(qū)。其中，一維河道模型的構(gòu)建范圍為石頭口門水庫(kù)到飲馬河河口，范圍內(nèi)的主干河道長(zhǎng)度為110 km；二維水流模型的構(gòu)建范圍為飲馬河沿岸的防洪保護(hù)區(qū)。

圖1 飲馬河(石頭口門水庫(kù)以下)水系示意圖

1 洪水分析方法及模型構(gòu)建

1.1 洪水分析方法

飲馬河石頭口門水庫(kù)放流與區(qū)間小南河、霧開河與伊通河等支流組成編制區(qū)域的洪水來(lái)源，干流來(lái)水控制站為石頭口門水庫(kù)站和德惠站，支流來(lái)水控制站為九臺(tái)站、農(nóng)安站等。

采用MIKE ZERO系列洪水模擬水動(dòng)力軟件進(jìn)行洪水分析，包含河道洪水演進(jìn)、潰堤洪水以及淹沒(méi)區(qū)洪水演進(jìn)，分別建立MIKE11一維非恒定流模型模擬飲馬河河道洪水演進(jìn)過(guò)程，DAMBRK方法計(jì)算堤防潰堤洪水，MIKE21二維非恒定流模型模擬飲馬河防洪保護(hù)區(qū)的洪水演進(jìn)過(guò)程[2]。由于二維水流模擬的上邊界為堤防潰口處的流量過(guò)程，而潰口流量又取決于潰口上下游(河道與淹沒(méi)區(qū))的水位差，因此將上述3種模型相互耦合、聯(lián)立求解，系統(tǒng)模擬潰堤后對(duì)防洪保護(hù)區(qū)的影響。

1.2 模型原理

1.2.1 一維水動(dòng)力學(xué)模型原理

圣維南(Saint-Venant)方程組為MIKE11一維非恒定流模型中的河道洪水運(yùn)動(dòng)的基本方程：

動(dòng)量方程：

式中：Q為斷面流量，m3/s；A為過(guò)水?dāng)嗝婷娣e，m2；q為源匯的單寬流量，m2/s；x為距離坐標(biāo)，m；t為時(shí)間坐標(biāo)，s；h為水位，m；C為謝才系數(shù)；R為水力半徑，m；g為重力加速度，m/s2；α為動(dòng)量校正系數(shù)。

1.2.2MIKE21二維水動(dòng)力模型

二維非恒定流方程如下，為平面二維水流模型，由連續(xù)方程及動(dòng)量方程組成：

式中：h為水深，m；Z為水位，m；u為X方向沿垂線平均的水平流速分量，m/s；v為Y方向沿垂線平均的水平流速分量，m/s；g為重力加速度，m/s2；n為糙率；q為源匯項(xiàng)，m/s。

1.3 模型構(gòu)建

1.3.1 網(wǎng)格剖分

收集編制范圍1∶10 000矢量地形圖及數(shù)字高程，確定模型中防洪保護(hù)區(qū)分析范圍，并設(shè)置開邊界和閉邊界，開邊界為水流進(jìn)出邊界，閉邊界為阻水邊界?？紤]無(wú)結(jié)構(gòu)網(wǎng)格適應(yīng)能力強(qiáng)并且靈活的特點(diǎn)，為提高邊界處模擬精度，采用無(wú)結(jié)構(gòu)三角形網(wǎng)格剖分[3]，控制三角形網(wǎng)格最大面積不超過(guò)0.01 km2，且網(wǎng)格大小均勻，盡量剖分成等邊三角形。見圖2、圖3。

圖2 網(wǎng)格劃分示意圖

圖3 局部網(wǎng)格劃分示意圖

1.3.2 計(jì)算時(shí)間和步長(zhǎng)

河道斷面間距、連接位置以及參數(shù)驗(yàn)證的水文站點(diǎn)位置等因素與一維河道模型空間步長(zhǎng)設(shè)定息息相關(guān)，飲馬河斷面間距最小為50 m，作為MIKE模型中水位計(jì)算點(diǎn)間距，為了保證模型穩(wěn)定準(zhǔn)確運(yùn)行，此次時(shí)間步長(zhǎng)設(shè)定為30 s。

1.3.3 邊界條件設(shè)置

根據(jù)編制范圍河網(wǎng)分布，上邊界位于石頭口門壩址，根據(jù)石頭口門水庫(kù)調(diào)度規(guī)則，調(diào)節(jié)天然入庫(kù)洪水過(guò)程線得到水庫(kù)泄流過(guò)程，下邊界位于飲馬河河口，根據(jù)河口處實(shí)測(cè)斷面分析計(jì)算得到水位-流量關(guān)系。區(qū)間小南河、霧開河與伊通河等支流作為點(diǎn)源以洪水過(guò)程線形式加載。

1.3.4 模型參數(shù)選取與驗(yàn)證

根據(jù)飲馬河各斷面河道河床的實(shí)際組成情況，并通過(guò)石頭口門站、德惠站洪痕調(diào)查以及H-Q關(guān)系曲線反推糙率驗(yàn)證而得到的糙率成果，主槽值選用0.03，邊灘選用0.05。

將防洪保護(hù)區(qū)按照土地利用圖分為6種地表類型，各類下墊面糙率參考《洪水風(fēng)險(xiǎn)圖編制導(dǎo)則》有關(guān)數(shù)值，并根據(jù)典型區(qū)域?qū)崪y(cè)大洪水淹沒(méi)資料進(jìn)行分析率定。見表1。

表1 飲馬河防洪保護(hù)區(qū)糙率表

1.3.5 模型驗(yàn)證

以飲馬河石頭口門水庫(kù)2010年7月21日實(shí)測(cè)泄流為上邊界，結(jié)合水文年鑒中洪水要素摘要實(shí)測(cè)流量，模擬大汛期洪水，進(jìn)行模型參數(shù)驗(yàn)證。通過(guò)調(diào)整飲馬河各斷面的河床主槽及邊灘的糙率參數(shù)，使干流區(qū)間內(nèi)的德惠水文站處經(jīng)模型模擬得到的流量和水位與實(shí)測(cè)值接近吻合。當(dāng)洪峰流量Qm和水位H的誤差分別控制在10%和20 cm以內(nèi)[4]，說(shuō)明河道參數(shù)取值較為合理。見圖4。

利用已經(jīng)參數(shù)驗(yàn)證后的一維河道模型，采用飲馬河干流及各支流50年一遇設(shè)計(jì)洪水(以1953年作為典型年，干流洪峰流量Qm為700 m3/s)作為邊界條件，分析計(jì)算飲馬河石頭口門水庫(kù)壩下至河口段河道水面線，并與已通過(guò)上級(jí)水利主管部門審查的50年一遇設(shè)計(jì)水面線成果進(jìn)行對(duì)比(圖5)，模擬水面線成果與堤防設(shè)計(jì)成果平均誤差為10 cm，最大誤差為21.4 cm，驗(yàn)證結(jié)果符合要求。

圖4 德惠站模擬值與實(shí)測(cè)值對(duì)比圖(左圖為流量，右圖為水位)

圖5 模擬與設(shè)計(jì)水面線對(duì)比圖(P=2%)

2 洪水分析成果

2.1洪水分析計(jì)算

以陳家崴子潰口20年一遇洪水方案為例，飲馬河干流按20年一遇洪水控制，各支流相應(yīng)洪水組合，洪水過(guò)程2015年8月2日0時(shí)至9月20日6時(shí)。8月28日8時(shí)，德惠市農(nóng)村段陳家崴子險(xiǎn)工處河道達(dá)到特征水位，堤防開始發(fā)生潰決，最大潰口寬度50 m，潰口最大流量Qm=297 m3/s。8月24日12時(shí)，飲馬河下游里程75 606 m處有小部分水漫出，其余堤防均未發(fā)生漫堤現(xiàn)象。至洪水過(guò)程結(jié)束，共進(jìn)入飲馬河防洪保護(hù)區(qū)總水量4 451×104m3，最大淹沒(méi)范圍61.22 km2，最大淹沒(méi)水深達(dá)4.42 m，最大洪水流速達(dá)1.29 m/s。見圖6、圖7。

2.2 水量平衡分析

2.2.1 淹沒(méi)區(qū)水量平衡分析

計(jì)算模型中的進(jìn)水量、出流量以及防洪保護(hù)區(qū)內(nèi)淹沒(méi)的水量，并將入流量減去出流量后得到的差值，與防洪保護(hù)區(qū)內(nèi)淹沒(méi)的水量對(duì)比，若誤差在1%以內(nèi)，說(shuō)明洪水模擬結(jié)果合理。經(jīng)計(jì)算，陳家崴子潰口方案進(jìn)水量11 209.78×104m3，出水量6 759.26×104m3，最終水量4 451.41×104m3，相對(duì)誤差0.02%。

圖6 飲馬河陳家崴子潰口流量過(guò)程圖

圖7 飲馬河陳家崴子潰口最大淹沒(méi)水深示意圖

2.2.2 耦合模型水量平衡分析

在一、二維耦合模型計(jì)算中，潰堤模型作為橋梁，將一維水動(dòng)力模型水量傳遞給二維水動(dòng)力模型。由于是三者耦合計(jì)算，如模型設(shè)置不合理，會(huì)造成水量缺失現(xiàn)象，影響結(jié)果的可信度。通過(guò)對(duì)比一維模型輸入二維模型水量與二維模型接收一維模型水量情況檢驗(yàn)計(jì)算合理性，陳家崴子潰口發(fā)生20年一遇洪水時(shí)，一維水動(dòng)力模型河道出流水量11 209.78×104m3，二維水動(dòng)力模型邊界入流水量11 209.78×104m3，水量相差為0，說(shuō)明模型計(jì)算合理[5]。

3 結(jié) 語(yǔ)

通過(guò)MIKE系列軟件構(gòu)建一、二維水動(dòng)力模型，可以直觀展示洪水在發(fā)生潰堤后產(chǎn)生的最大淹沒(méi)范圍、最大淹沒(méi)水深及最大洪水流速等信息，同時(shí)也可展示洪水演進(jìn)動(dòng)態(tài)過(guò)程，能夠提取洪水到達(dá)時(shí)間、洪水淹沒(méi)歷時(shí)等洪水風(fēng)險(xiǎn)要素。本文以飲馬河陳家崴子潰口20年一遇洪水方案為例，介紹了MIKE11和MIKE21模型原理，模型構(gòu)建中的網(wǎng)格剖分、設(shè)置邊界條件、模型參數(shù)驗(yàn)證等，建立具有較高模擬精度的水動(dòng)力耦合模型，模擬成果經(jīng)水量平衡分析后合理可靠，成果可為后續(xù)的洪災(zāi)損失評(píng)估和防洪減災(zāi)提供技術(shù)支撐。