为克服传统的河道洪水演进模型中存在的人为误差和模拟精度低的问题,引入MIKE11对河道洪水过程进行数值模拟,分析区间洪水传播时间和流量演进规律。以汉江上游河段为例,分别建立MIKE11HD和马斯京根的洪水演进模型。结果对比表明:MIKE11不仅可应用于汉江流域的洪水模拟,且模拟结果达到了甲级预报精度,计算结果验证了MIKE11模型的可靠性、准确性。将两种模型结果和实测数据进行对比,绘制了安康电站–蜀河电站的流量与传播时间的函数关系,揭示了洪水演进规律。研究成果为实现下游河道防洪、旬阳水库建设、梯级水库调度运行的安全提供科学依据。 In order to overcome the shortcomings of artificial error and low precision in the traditional flood routing models, the MIKE11 was used to simulate the flood process of watercourse and analyze the flood propagation time and discharge evolution rule. Taking the upstream of Hanjiang River as an example, the MIKE11 HD flood evolution model and Muskingum method were established, respectively. The results show that the MIKE11 HD not only can be applied to flood simulation in the upstream of Hanjiang River, but also reaches the first-order accuracy. The calculated results confirm the reliability and accuracy of the MIKE11 HD model. By comparing the results of the two models with the measured data, the function relationship between the flow and the propagation time of Ankang-Shuhe hydropower stations are plotted, and the law of flood evolution is revealed. These results provide scientific basis for the safety of downstream flood control, construction of Xunyang reservoir, and operation of cascade reservoirs.
高超,白涛*,杨旺旺,畅建霞,张飒
西安理工大学,西北旱区生态水利工程国家重点实验室,陕西 西安
收稿日期:2017年3月29日;录用日期:2017年4月15日;发布日期:2017年4月18日
为克服传统的河道洪水演进模型中存在的人为误差和模拟精度低的问题,引入MIKE11对河道洪水过程进行数值模拟,分析区间洪水传播时间和流量演进规律。以汉江上游河段为例,分别建立MIKE11HD和马斯京根的洪水演进模型。结果对比表明:MIKE11不仅可应用于汉江流域的洪水模拟,且模拟结果达到了甲级预报精度,计算结果验证了MIKE11模型的可靠性、准确性。将两种模型结果和实测数据进行对比,绘制了安康电站–蜀河电站的流量与传播时间的函数关系,揭示了洪水演进规律。研究成果为实现下游河道防洪、旬阳水库建设、梯级水库调度运行的安全提供科学依据。
关键词 :MIKE11模型,马斯京根,糙率率定,流量演进,传播时间
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
河道洪水演进是根据断面的入流条件,以及洪水波的传播、坦化、变形规律,采用各种技术手段求得下断面的出流过程。科学、准确地预报河道洪水传播时间、模拟洪水涨退过程、揭示洪水演进规律,能够尽早地发布洪水警报,使洪水受淹区的居民及时撤离,最大程度地降低洪水灾害。
目前,分析河道流量演进常用马斯京根法。该方法以其简单的模型结构,可靠的模拟效果,成为使用广泛的河道洪水演算的水文学方法。随着科学技术的发展,MIKE11水动力学模型开始进入应用领域。其中MIKE11 HD模块根据水力学原理,借助计算机工具数值求解圣维南方程组,在国内外科学研究和工程应用的广泛领域都获得了令人满意的效果,已成为多个国家的标准工具。近年来,国内用马斯京根法和MIKE11模型在不同流域进行了有关河道流量演进规律的研究:严伏朝等 [
汉江是长江最大支流,上游正在建设实施引汉济渭跨流域调水工程,中下游的南水北调中线工程通水两年,累计供水超过60亿m3。加上汉江上游在建的旬阳水电站,汉江上游已成为我国重点水利建设和跨流域调水的核心区域。但针对汉江上游洪水演进规律的研究成果较少,特别是汛期不同量级洪水预报、演进规律的研究少之又少。因此,亟需开展汉江上游河段洪水演进的研究。本文针对当前安康市防洪安全、旬阳水电站建设期的防洪安全、梯级水库群洪水资源化以及汉江上游流域防洪风险面临的实际问题,拟采用MIKE11 HD水动力学模型、马斯京根法模拟河道流量演进过程,以获得安康水库不同下泄流量到达下游安康城区、旬阳、蜀河的传播时间;阐明不同遇见期洪水到达下游各断面的传播时间函数关系,揭示区间洪水的传播规律。研究成果对于蜀河电站的防洪预警、合理调度和安康城区防洪撤离、旬阳电站预停施工具有重要的参考价值,可为梯级水库群洪水资源化以及梯级水库的安全调度运行提供科学依据。
MIKE11系列软件是由丹麦水动力研究所开发研制,属分布式水文模型,软件的功能涉及范围从降雨、产流、河流、城市、河口、近海、深海等水文循环的整个过程,目前,MIKE11软件系列在国内很多水利工程中成功应用,包括长江流域、松花江流域、黄河流域等 [
洪水波的演进与变形可用一元非恒定流的基本方程组描述,即圣维南方程组,其连续性和运动方程具如下:
式中: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 [
在连续方程中,Q仅对x求偏导,所以方程可以写成以h为中心的形式,见图2;动量方程则以Q为中心,见图3。
MIKE11软件内部通过河网文件(.nwk11)、断面文件(.xns11)、边界条件(.bnd11)、模型参数文件(.hd11)和时间序列文件(.dfs0)进行组织和实现文件输入。具体资料主要包括:1) 河网形状;2) 河道和滩区断面数据;3) 模拟边界处水文实测数据,本文采用水位过程和流量过程。具体的建模步骤如图4所示。
马斯京根流量演算法是以运动波理论为基础,主要使用连续方程和简化或是近似处理的动力方程联解。通过线性假定建立槽蓄方程,联合水量平衡方程可得如下流量演算方程 [
其中:
图1. 计算河道上节点布置
图2. 连续方程Abbott六点隐式差分格式
图3. 动量方程Abbott六点隐式差分格式
图4. MIKE11 HD水动力模型建模流程图
式中:Q1,Q2,I1,I2分别为研究河段时段始末出流断面流量和时段始末入流断面流量,K为槽蓄曲线坡度,x为流量因子。
安康水库位于汉江上游安康市境内,是一座以发电为主,兼顾防洪、航运等综合利用的大型水利枢纽,也是汉江上游梯级水库的龙头水库。坝址位于安康市城西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。
本文以汉江流域DEM数据为基础资料,应用ARCGIS10.2软件建立河网文件底图。通过查阅相关水文手册和水电站实测资料,得到详细的河道断面资料和边界处的小时尺度流量和水位资料。数据资料包括安康电站出库径流资料、安康电站–蜀河电站区间支流径流资料、旬阳电站坝址径流资料和蜀河电站坝址径流资料。研究河段为汉江上游安康电站至蜀河电站,区间共计四条支流汇入,分别为月河、黄洋河、坝河和旬河。支流水文站依次为长枪铺、县河口、桂花园和向家坪,安康电站下游18 km处设有安康水文站。
图5. 汉江安康水电站–蜀河水电站河段示意图
蜀河电站2010年建成,运行时间较短,收集到的电站运行资料时间长度为2010~2014年。依据河道洪水演算传播时间规律,同样选取上游安康电站2010~2014年的运行资料。
MIKE11作为一维河网水力模型,在进行河道洪水演算时,需要对反映河床表面阻水作用的糙率n值进行参数率定 [
确定糙率后,对模型精度进行验证。设置安康水文站断面为对比断面,选取20120907场次洪水模拟安康水文站断面流量过程,结果如图7所示。
通过分析图7得到:MIKE11 HD模型对安康电站-安康水文站河段的汛期洪水过程的模拟结果和实测较为接近,模拟洪水过程线与实测洪水拟合较好,率定的参数能够很好的模拟洪水过程。采用MIKE11软件对安康电站-安康水文站河段的汛期洪水过程的模拟和实测较为接近,模拟洪水流量过程线与实测洪水拟合较好,能够
图6. 不同糙率对应确定性系数
图7. 20120907次洪水模拟过程
很好的模拟洪水过程,且洪水退水段模拟精度高于涨水段模拟精度。
由上述MIKE11 HD模型的建立和验证可知,MIKE11 HD模型在河道流量的数值仿真模拟上合理可靠而且精度较高。因此,本文采用上述模型,选取安康电站–安康城区河段8场不同洪水量级的洪水进行模拟分析,结果见表1。
通过分析MIKE11洪水模拟结果得到:
1) 模拟洪水总量最大误差为7.26%,最小误差仅为−0.72%,模拟洪峰与实测值最大相对误差为8.8%,最小值为1.55%,洪水总量和洪峰值模拟结果与实际比较接近。
2) 因为模型选取的洪水资料为小时尺度,洪峰传达时差均小于2 h,多数洪水模拟的时差小于1 h,所以洪峰传达历时结果相对准确。
综合比较得到结论:MIKE11模型可以很好的应用于研究区河道洪水演进问题中。
受资料长度的影响,本文采用分段马斯京根数学模型对安康–安康城区河段洪水进行模拟,釆用洪峰流量、洪量、洪峰时差、确定性系数4个指标作为评定预报结果质量的标准分析,结果见表2。
场次 评价项 | 洪峰流量/(m3/s) | 洪量/(亿m3) | 洪峰时差/(h) | 确定性系数 | ||||
---|---|---|---|---|---|---|---|---|
模拟 | 实测 | 误差(%) | 模拟 | 实测 | 误差(%) | |||
20110917 | 19819 | 20700 | −4.26 | 23.76 | 25.63 | −7.30 | <1 | 0.95 |
20120901 | 7234 | 7480 | −3.29 | 7.68 | 8.23 | −6.88 | 2 | 0.94 |
20120707 | 14361 | 13200 | 8.8 | 29.67 | 31.39 | −5.48 | 1 | 0.95 |
20120909 | 8355 | 8080 | 3.40 | 10.99 | 11.07 | −0.72 | 1 | 0.97 |
20140928 | 6947 | 6670 | 4.15 | 5.90 | 5.82 | 1.25 | <1 | 0.98 |
20140914 | 3649 | 3770 | −3.22 | 1.48 | 1.59 | −6.94 | <1 | 0.95 |
20110707 | 2764 | 2650 | 4.30 | 1.13 | 1.06 | 7.26 | <1 | 0.93 |
20130722 | 7799 | 7680 | 1.55 | 9.09 | 9.32 | −2.44 | 2 | 0.95 |
表1. MIKE11洪水模拟结果
场次 评价项 | 洪峰流量/(m3/s) | 洪量/(亿m3) | 洪峰时差/(h) | 确定性系数 | ||||
---|---|---|---|---|---|---|---|---|
模拟 | 实测 | 误差(%) | 模拟 | 实测 | 误差(%) | |||
20110917 | 13997 | 12556 | 11.5 | 28.0 | 28.5 | −1.8 | 1 | 0.92 |
20120901 | 7493 | 7931 | −5.5 | 8.13 | 8.51 | −4.5 | <1 | 0.93 |
20120707 | 8223 | 7301 | 12.6 | 10.1 | 9.1 | 11.0 | 1 | 0.87 |
20120909 | 8475 | 8450 | 0.3 | 17.4 | 17.8 | −2.2 | 1 | 0.87 |
20140928 | 2887 | 3127 | −7.7 | 2.56 | 2.92 | −12.3 | 2 | 0.74 |
20140914 | 7061 | 7214 | −2.1 | 12.9 | 13.7 | −5.8 | 3 | 0.81 |
20110707 | 5469 | 6201 | −11.8 | 4 | 4.23 | −5.4 | 2 | 0.81 |
20130722 | 7581 | 6903 | 9.8 | 11.4 | 10.5 | 8.6 | <1 | 0.86 |
表2. 马斯京根模拟结果分析表
将上述两种模型计算的传播时间和实测数据进行对比,如图8所示。绘制区间流量–时间函数关系图,如图9所示,得到了不同量级洪水从安康水库到蜀河水库的传播时间方程。考虑两种方法所选洪水场次受资料的尺度和方法限制,故同等量级的洪水场次选取的不完全相同。尽管同等量级选取不同场次洪水进行模拟,但对于洪水传播时间的规律分析的影响可忽略不计。通过比较流量和传播时间关系得到:
通过比较流量和传播时间关系得到:
1) MIKE11 HD模型演算结果相比马斯京根模型演算结果偏大,实测结果比MIKE11 HD模型结果偏大,对于汛期的河道洪水过程,MIKE11 HD的结果和实测的相关系数约为0.98,略高于马斯京根法的结果,其计算的洪水传播时间比马斯京根模型更接近实测值。
2) 由于马斯京根模型需要同时知道上下游实测流量过程才能模拟,对资料依赖程度高,而实际中间断面无实测资料,导致马斯京根模型受限。而MIKE11模型只需要上下游边界的资料就能计算中间各断面的流量和水位过程,所以MIKE11模型的应用更为便捷。
3) 通过MIKE11 HD模型计算结果得到了安康–蜀河区间流量与传播时间对应的函数关系,通过这个函数
图8. 马斯京根和MIKE11结果比较
图9. 区间流量–时间函数关系
关系可以计算出不同下泄流量在区间内的传播时间。
本文分别采用传统的马斯京根模型和MIKE11模型对不同量级的安康出库洪水进行模拟,得到以下结论:
1) MIKE11模型对场次洪水进行模拟,模拟得到洪水过程到达安康断面的流量过程,与实际流量过程进行对比得出该河段MIKE11模型的构建合理可靠,所率定的参数能够反映流域的真实情况,模拟效果良好。
2) 通过对马斯京根法和MIKE11模型模拟的效果进行对比,结果表明:两种方法都能够用于预报模拟河段的洪水过程;但MIKE11模型具有明确的物理意义,准确地反映模拟河段的河道特征,可以加入溢流堰、涵洞、测站等水工建筑物,对河道的刻画更加真实、灵活,与马斯京根模拟对比,MIKE11模拟的平均确定性系数为0.95,优于马斯京根模拟的0.85,验证了MIKE11 HD模型的优越性,模拟效果优于马斯京根法,能够更好地适用于山区河流的河道洪水预报。
3) 得到了洪水在不同断面传播时间和流量的函数关系、各个断面的流量、水位变化的过程,揭示了洪水在不同断面的演进规律。研究成果为河道洪水预报、汉江上游城市防洪预警、下游水电站的安全运行以及施工提供了技术参考和科学依据。
国家自然科学基金(51409210,91325201);水利部公益性行业科研专项(201501058);陕西省水利科技计划项目(2016slkj-8,2017slkj-16),西安理工大学水利水电学院科研项目(2016ZZKT-15)。
高 超,白 涛,杨旺旺,畅建霞,张 飒. 基于MIKE11的汉江上游洪水演进规律研究 Study on Flood Propagation in the Upstream of Hanjiang River Using MIKE11 Model[J]. 水资源研究, 2017, 06(02): 156-165. http://dx.doi.org/10.12677/JWRR.2017.62019