1. 引言
河道洪水演进是根据断面的入流条件,以及洪水波的传播、坦化、变形规律,采用各种技术手段求得下断面的出流过程。科学、准确地预报河道洪水传播时间、模拟洪水涨退过程、揭示洪水演进规律,能够尽早地发布洪水警报,使洪水受淹区的居民及时撤离,最大程度地降低洪水灾害。
目前,分析河道流量演进常用马斯京根法。该方法以其简单的模型结构,可靠的模拟效果,成为使用广泛的河道洪水演算的水文学方法。随着科学技术的发展,MIKE11水动力学模型开始进入应用领域。其中MIKE11 HD模块根据水力学原理,借助计算机工具数值求解圣维南方程组,在国内外科学研究和工程应用的广泛领域都获得了令人满意的效果,已成为多个国家的标准工具。近年来,国内用马斯京根法和MIKE11模型在不同流域进行了有关河道流量演进规律的研究:严伏朝等 [1] 利用马斯京根法的正向演算和基于经验的断面流量相关法的反向演算研究渭河下游小流量演进规律;白涛等 [2] 在分析防凌期流量演进特点的基础上,考虑封开河日期、槽蓄水增量、气温等凌情影响因子提出了分期分河段的混合算法以改进传统演进方法;Rahman等 [3] 在孟加拉国提斯塔河应用频率分析软件和MIKE11模型确定了不同重现期的设计洪水流量和相应的洪水演进过程;吴天蛟等 [4] 应用MIKE11软件建立有区间入流的三峡库区水力学洪水演进模型,得到对应的出库流量、坝前水位的模拟值,预报精度达到甲级标准;林波等 [5] 将MIKE11水动力模型和降雨径流模型(NAM)耦合,应用到三江平原挠力河流域,结果表明MIKE11/NAM模型可以对挠力河流域降雨径流过程进行较好的模拟。
汉江是长江最大支流,上游正在建设实施引汉济渭跨流域调水工程,中下游的南水北调中线工程通水两年,累计供水超过60亿m3。加上汉江上游在建的旬阳水电站,汉江上游已成为我国重点水利建设和跨流域调水的核心区域。但针对汉江上游洪水演进规律的研究成果较少,特别是汛期不同量级洪水预报、演进规律的研究少之又少。因此,亟需开展汉江上游河段洪水演进的研究。本文针对当前安康市防洪安全、旬阳水电站建设期的防洪安全、梯级水库群洪水资源化以及汉江上游流域防洪风险面临的实际问题,拟采用MIKE11 HD水动力学模型、马斯京根法模拟河道流量演进过程,以获得安康水库不同下泄流量到达下游安康城区、旬阳、蜀河的传播时间;阐明不同遇见期洪水到达下游各断面的传播时间函数关系,揭示区间洪水的传播规律。研究成果对于蜀河电站的防洪预警、合理调度和安康城区防洪撤离、旬阳电站预停施工具有重要的参考价值,可为梯级水库群洪水资源化以及梯级水库的安全调度运行提供科学依据。
2. 模型介绍
2.1. MIKE11 HD水动力模型
MIKE11系列软件是由丹麦水动力研究所开发研制,属分布式水文模型,软件的功能涉及范围从降雨、产流、河流、城市、河口、近海、深海等水文循环的整个过程,目前,MIKE11软件系列在国内很多水利工程中成功应用,包括长江流域、松花江流域、黄河流域等 [6] [7] [8] [9] [10] 。本文河道洪水演算主要应用MIKE11 HD水动力模块。
2.1.1. 模型原理和求解方法
洪水波的演进与变形可用一元非恒定流的基本方程组描述,即圣维南方程组,其连续性和运动方程具如下:
(1)
式中:A为过水断面,m2;Q为流量,m3/s;t为时间,s;x为距离,m;h为断面平均水深,m;C为谢才系数,m1/2/s;g为重力加速度,m2/s;q为侧向流入流量,m3/s;R为水力半径,m。
MIKE11对于一维河道非恒定流的基本方程组式(1)应用Abbot六点隐式差分进行求解。河道上横断面(节点)按照水位–流量–水位的顺序交替布置,如图1所示,即h和q不在同一个计算节点同时出现,q点总是布置在相邻的h点之间,距离可以不相同。然后,在每个时间步长内,利用隐式格式的有限分法交替计算q和h [11] 。
在连续方程中,Q仅对x求偏导,所以方程可以写成以h为中心的形式,见图2;动量方程则以Q为中心,见图3。
2.1.2. MIKE11建模流程图
MIKE11软件内部通过河网文件(.nwk11)、断面文件(.xns11)、边界条件(.bnd11)、模型参数文件(.hd11)和时间序列文件(.dfs0)进行组织和实现文件输入。具体资料主要包括:1) 河网形状;2) 河道和滩区断面数据;3) 模拟边界处水文实测数据,本文采用水位过程和流量过程。具体的建模步骤如图4所示。
2.2. 马斯京根流量演进模型
马斯京根流量演算法是以运动波理论为基础,主要使用连续方程和简化或是近似处理的动力方程联解。通过线性假定建立槽蓄方程,联合水量平衡方程可得如下流量演算方程 [12]
(2)
其中:
![](//html.hanspub.org/file/7-2410485x12_hanspub.png)
![](//html.hanspub.org/file/7-2410485x13_hanspub.png)
Figure 1. Node arrangement on river reach
图1. 计算河道上节点布置
![](//html.hanspub.org/file/7-2410485x14_hanspub.png)
Figure 2. Continuous equation Abbott six-point implicit difference scheme
图2. 连续方程Abbott六点隐式差分格式
![](//html.hanspub.org/file/7-2410485x15_hanspub.png)
Figure 3. Momentum equation Abbott six-point implicit difference scheme
图3. 动量方程Abbott六点隐式差分格式
![](//html.hanspub.org/file/7-2410485x16_hanspub.png)
Figure 4. Flow chart for the MIKE11 HD hydrodynamic model
图4. MIKE11 HD水动力模型建模流程图
式中:Q1,Q2,I1,I2分别为研究河段时段始末出流断面流量和时段始末入流断面流量,K为槽蓄曲线坡度,x为流量因子。
3. 实例应用
3.1. 区域概况
安康水库位于汉江上游安康市境内,是一座以发电为主,兼顾防洪、航运等综合利用的大型水利枢纽,也是汉江上游梯级水库的龙头水库。坝址位于安康市城西18 km处,控制流域集水面积35,700 km2。目前,安康水库至蜀河水库区间河段旬阳处正在兴建汉江上游规划的第五级水电站——旬阳水电站,其位于陕西省旬阳县城南约2 km,上距已建的安康水电站65 km,下距蜀河水电站55 km。本文研究河段为以安康水电站为入流控制断面,以蜀河水电站为出流控制断面的安康_蜀河河段。安康水库至蜀河水库全长109.36 km,占汉江干流全河段的7%,落差61 m,河床平均比降为0.06%。
安康至蜀河河段包括月河、黄洋河、坝河、旬河四条支流,支流均设有水文站。已建的水电站有安康电站和蜀河电站,在建的水电站有旬阳电站。由于流域内山高坡陡,岩层透水性差,加之河槽又缺乏良好的调节作用,洪水过程消涨极快,峰形尖瘦、一次洪水过程一般历时5~7天。洪水年际变化较大,安康1983年实测最大洪峰流量为31,000 m3/s,1966年最大洪峰流量仅3500 m3/s,前者是后者的8.7倍,本流域洪水变差系数较大。洪水季节性特点表现为:夏季强度大、历时短,秋季洪量大、历时长,河段示意如图5。
3.2. MIKE11模型及其计算结果
3.2.1. 数据处理
本文以汉江流域DEM数据为基础资料,应用ARCGIS10.2软件建立河网文件底图。通过查阅相关水文手册和水电站实测资料,得到详细的河道断面资料和边界处的小时尺度流量和水位资料。数据资料包括安康电站出库径流资料、安康电站–蜀河电站区间支流径流资料、旬阳电站坝址径流资料和蜀河电站坝址径流资料。研究河段为汉江上游安康电站至蜀河电站,区间共计四条支流汇入,分别为月河、黄洋河、坝河和旬河。支流水文站依次为长枪铺、县河口、桂花园和向家坪,安康电站下游18 km处设有安康水文站。
![](//html.hanspub.org/file/7-2410485x17_hanspub.png)
Figure 5. Schematic diagram of the reach of Ankang-Shuhe Hydropower Stations
图5. 汉江安康水电站–蜀河水电站河段示意图
蜀河电站2010年建成,运行时间较短,收集到的电站运行资料时间长度为2010~2014年。依据河道洪水演算传播时间规律,同样选取上游安康电站2010~2014年的运行资料。
3.2.2. 参数率定
MIKE11作为一维河网水力模型,在进行河道洪水演算时,需要对反映河床表面阻水作用的糙率n值进行参数率定 [13] 。为确定模型对糙率的敏感程度,n值设定初值为0.042,选取步长值0.002依次进行模拟,选取20120707、20120907、20100822三场洪水进行河道洪水演算,对模型参数n进行率定,通过对比安康断面的实测水位和模拟水位,根据GB/T22482-2008《水文预报情报规范》 [14] 计算确定性系数Dc,对模拟结果进行评价,结果如图6所示。随着糙率n值从0.042增加到0.048,安康断面模拟水位过程不断变化,在n = 0.045时,模拟的水位过程和实际水位过程拟合相对较好,综合确定0.045为模型参数糙率的最优值。
确定糙率后,对模型精度进行验证。设置安康水文站断面为对比断面,选取20120907场次洪水模拟安康水文站断面流量过程,结果如图7所示。
通过分析图7得到:MIKE11 HD模型对安康电站-安康水文站河段的汛期洪水过程的模拟结果和实测较为接近,模拟洪水过程线与实测洪水拟合较好,率定的参数能够很好的模拟洪水过程。采用MIKE11软件对安康电站-安康水文站河段的汛期洪水过程的模拟和实测较为接近,模拟洪水流量过程线与实测洪水拟合较好,能够
![](//html.hanspub.org/file/7-2410485x18_hanspub.png)
Figure 6. Deterministic coefficient of different roughness
图6. 不同糙率对应确定性系数
![](//html.hanspub.org/file/7-2410485x19_hanspub.png)
Figure 7. Simulated flood process of 20120907
图7. 20120907次洪水模拟过程
很好的模拟洪水过程,且洪水退水段模拟精度高于涨水段模拟精度。
3.2.3. 计算结果
由上述MIKE11 HD模型的建立和验证可知,MIKE11 HD模型在河道流量的数值仿真模拟上合理可靠而且精度较高。因此,本文采用上述模型,选取安康电站–安康城区河段8场不同洪水量级的洪水进行模拟分析,结果见表1。
通过分析MIKE11洪水模拟结果得到:
1) 模拟洪水总量最大误差为7.26%,最小误差仅为−0.72%,模拟洪峰与实测值最大相对误差为8.8%,最小值为1.55%,洪水总量和洪峰值模拟结果与实际比较接近。
2) 因为模型选取的洪水资料为小时尺度,洪峰传达时差均小于2 h,多数洪水模拟的时差小于1 h,所以洪峰传达历时结果相对准确。
综合比较得到结论:MIKE11模型可以很好的应用于研究区河道洪水演进问题中。
3.3. 马斯京根模型计算结果
受资料长度的影响,本文采用分段马斯京根数学模型对安康–安康城区河段洪水进行模拟,釆用洪峰流量、洪量、洪峰时差、确定性系数4个指标作为评定预报结果质量的标准分析,结果见表2。
![](Images/Table_Tmp.jpg)
Table 1. MIKE11 HD model flood simulation results
表1. MIKE11洪水模拟结果
![](Images/Table_Tmp.jpg)
Table 2. Simulated results using Muskingum method
表2. 马斯京根模拟结果分析表
4. 流量传播规律
将上述两种模型计算的传播时间和实测数据进行对比,如图8所示。绘制区间流量–时间函数关系图,如图9所示,得到了不同量级洪水从安康水库到蜀河水库的传播时间方程。考虑两种方法所选洪水场次受资料的尺度和方法限制,故同等量级的洪水场次选取的不完全相同。尽管同等量级选取不同场次洪水进行模拟,但对于洪水传播时间的规律分析的影响可忽略不计。通过比较流量和传播时间关系得到:
通过比较流量和传播时间关系得到:
1) MIKE11 HD模型演算结果相比马斯京根模型演算结果偏大,实测结果比MIKE11 HD模型结果偏大,对于汛期的河道洪水过程,MIKE11 HD的结果和实测的相关系数约为0.98,略高于马斯京根法的结果,其计算的洪水传播时间比马斯京根模型更接近实测值。
2) 由于马斯京根模型需要同时知道上下游实测流量过程才能模拟,对资料依赖程度高,而实际中间断面无实测资料,导致马斯京根模型受限。而MIKE11模型只需要上下游边界的资料就能计算中间各断面的流量和水位过程,所以MIKE11模型的应用更为便捷。
3) 通过MIKE11 HD模型计算结果得到了安康–蜀河区间流量与传播时间对应的函数关系,通过这个函数
![](//html.hanspub.org/file/7-2410485x20_hanspub.png)
Figure 8. Comparison of Muskingum and MIKE11 HD results
图8. 马斯京根和MIKE11结果比较
![](//html.hanspub.org/file/7-2410485x21_hanspub.png)
Figure 9. Relationship of interval discharge-time function
图9. 区间流量–时间函数关系
关系可以计算出不同下泄流量在区间内的传播时间。
5. 结论
本文分别采用传统的马斯京根模型和MIKE11模型对不同量级的安康出库洪水进行模拟,得到以下结论:
1) MIKE11模型对场次洪水进行模拟,模拟得到洪水过程到达安康断面的流量过程,与实际流量过程进行对比得出该河段MIKE11模型的构建合理可靠,所率定的参数能够反映流域的真实情况,模拟效果良好。
2) 通过对马斯京根法和MIKE11模型模拟的效果进行对比,结果表明:两种方法都能够用于预报模拟河段的洪水过程;但MIKE11模型具有明确的物理意义,准确地反映模拟河段的河道特征,可以加入溢流堰、涵洞、测站等水工建筑物,对河道的刻画更加真实、灵活,与马斯京根模拟对比,MIKE11模拟的平均确定性系数为0.95,优于马斯京根模拟的0.85,验证了MIKE11 HD模型的优越性,模拟效果优于马斯京根法,能够更好地适用于山区河流的河道洪水预报。
3) 得到了洪水在不同断面传播时间和流量的函数关系、各个断面的流量、水位变化的过程,揭示了洪水在不同断面的演进规律。研究成果为河道洪水预报、汉江上游城市防洪预警、下游水电站的安全运行以及施工提供了技术参考和科学依据。
基金项目
国家自然科学基金(51409210,91325201);水利部公益性行业科研专项(201501058);陕西省水利科技计划项目(2016slkj-8,2017slkj-16),西安理工大学水利水电学院科研项目(2016ZZKT-15)。