摘要: 应用HEC-RAS和MIKE11模型,对成都市锦江流域的南河和磨底河进行了模拟分析。结果表明:同一模型不同计算方法时,恒定流和非恒定流计算结果差值较小;不同模型同一计算方法时,计算差值较大;在同一小流量下,HEC-RAS计算的磨底河水位比MIKE11的高;在同一大流量下,HEC-RAS计算的南河水位比MIKE11的低。非恒定流在模拟水流运动状态时,计算结果相对更加精准;建议锦江干流如府河、南河选用MIKE11计算水位,锦江支流如磨底河、浣花溪、沙河、桃花江等选用HEC-RAS计算水位。
Abstract:
In this paper, HEC-RAS and MIKE11 models are used to calculate the water surface line of the Nanhe and Modihe Rivers in Jinjiang basin. The results show as follows: Using the same model with different methods, the difference between constant flow and unsteady flow is small; Using different model with the same method, the difference is large; At the same low flow rate, the water surface line of Modihe River calculated by HEC-RAS is higher than MIKE11; At the same high flow rate, the water surface line of Nanhe River calculated by HEC-RAS is lower than MIKE11; The results of unsteady flow are more accurate simulating the flow movement. Calculating the water surface line of Jinjiang main stream, MIKE11 model is recommended; Calculating the water surface line of other tributaries like Modihe, Huanhuaxi, Shahe, and Taohuajiang Rivers, HEC-RAS model is recommended.
1. 引言
水面线推求是河道整治的基础工作,常用的方法如恒定流法和非恒定流法,常用的计算软件如天然河道水面线系统、HEC-RAS软件、SOBEK软件等 [1] 。近几年,MIKE11软件、InfoWorks RS也较常应用于水面线的推求 [2] 。如林杰等 [3] 以大富水应城至南垸河段为例,应用HEC-RAS软件,分析得出恒定流计算结果总体大于非恒定流计算结果,且流量越大差值呈现越大的趋势;蒋书武等 [4] 通过对HEC-RAS和MIKE11的对比分析,得出MIKE11具有较高的精度和稳定性;代兴勇等 [5] 通过MIKE11和HEC-RAS模拟水力半径对过流能力的影响,得出MIKE11软件计算得到的水面线与河工模型试验研究成果更为接近;蒋楠等 [6] 通过MIKE11和HEC-RAS模拟锦江水面线,得出MIKE11模型计算的水面线精度更高;安莉莉 [7] 对HEC-RAS和MIKE11的水位数值成果进行比较分析,得出当流量、糙率、纵坡为单一变量时,在小流量下,HEC-RAS比MIKE11的水位高,大流量时,MIKE11比HEC-RAS的水位高。
在以上研究的基础上,本文以成都锦江流域的干支流为例,研究HEC-RAS和MIKE11软件计算结果的不同并加以分析,推荐适合的软件计算成都锦江水面线。
2. 材料与方法
2.1. 试验区概况
锦江是岷江流经成都市区的两条主要河流——府河、南河的合称,其水源均来自都江堰。此外,流经成都市中心城区还有磨底河、浣花溪、桃花江、西郊河、饮马河等,皆称为锦江支流,成都市中心城区水系图见图1。
锦江城区堤防于1988年建成至今,年代较久,目前锦江正在开展流域治理工作,复核防洪水位成了其中重要的环节。本文主要以成都锦江流域的南河和磨底河为例,对HEC-RAS模型和MIKE11模型的计算结果进行对比分析。
2.2. 试验方法
2.2.1. 使用软件
HEC-RAS由美国陆军工程兵团水文工程中心开发,主要用于天然河道或明渠一维恒定流和非恒定流的水面线计算 [8] 。恒定流基于一维能量方程,通过迭代法逐段向上(下)游推算;非恒定流基于连续方程和动量方程,采用Newton-Raphson迭代法计算 [6] 。
![](//html.hanspub.org/file/7-2411203x8_hanspub.png?20230314084447800)
Figure 1. River system map of downtown Chengdu metropolis
图1. 成都市中心城区水系图
MIKE模型系列是由丹麦水资源及水环境研究所开发的软件,由MIKE11、MIKE21、MIKE SHE等模块,其中,MIKE11是一维河网模型,多运用于一维水面线推求。基于连续方程和动量方程,采用Abbott-Ionescu六点隐式有限差分格式进行求解 [6] 。
2.2.2. 计算对象和相关参数
本文选取锦江干流的南河和锦江支流的磨底河为分析对象。模型计算相关基本条件见表1。
![](Images/Table_Tmp.jpg)
Table 1. Relevant conditions for model calculation
表1. 模型计算相关条件
3. 结果与分析
3.1. 计算结果
为分析比较同一模型恒定流和非恒定流之间的差异及同一计算方法不同模型计算结果之间的差异,本文选取锦江流域的南河和磨底河进行分析。
1) 同一模型不同计算方法之间的对比
采用HEC-RAS和MIKE11的恒定流和非恒定流分别计算,对比同一模型不同计算方法的成果及对比,见表2、表3的序列 ⑤ 和序列 ⑥ 以及图2、图3。
![](Images/Table_Tmp.jpg)
Table 2. Comparison of water level of the same model with different methods in the Nanhe River
表2. 同一模型不同计算方法的南河水位对比
![](Images/Table_Tmp.jpg)
Table 3. Comparison table of water level of the same model with different methods in the Modihe River
表3. 同一模型不同计算方法的磨底河水位对比
结合表2和表3的序列 ⑤ 和序列 ⑥ 可以看出,同一模型恒定流与非恒定流的差异较小,MIKE11的恒定流计算结果较非恒定流高,最大差值为0.02 m;而HEC-RAS的恒定流计算结果较非恒定流低,最大差值为0.06 m。
2) 同一计算方法不同模型之间的对比
采用HEC-RAS和MIKE11的恒定流和非恒定流计算,对比不同模型同一计算方法的成果,见表2、表3的序列 ⑦ 和序列 ⑧ 以及图2、图3。
![](//html.hanspub.org/file/7-2411203x9_hanspub.png?20230314084447800)
Figure 2. Comparison water level of Nanhe River
图2. 南河计算水位对比图
![](//html.hanspub.org/file/7-2411203x10_hanspub.png?20230314084447800)
Figure 3. Comparison water level of Modihe River
图3. 磨底河计算水位对比图
结合以上图表可以看出,不同模型相同计算方法之间差值较大,南河用HEC-RAS的计算结果较MIKE11低,最大差值为0.71 m;而磨底河用HEC-RAS的整体计算结果较MIKE11高,最大差值为0.41 m,仅下游起推段出现低于MIKE11的情况。
3.2. 成因分析
1) 同一模型恒定流和非恒定流计算结果差异较小。其中,MIKE11恒定流计算比非恒定流高,而HEC-RAS恒定流计算结果比非恒定流低。主要原因在于MIKE11恒定流和非恒定流的解法均为六点隐式差分法,而非恒定流考虑了一定的河道调蓄,故恒定流的水位较高;而HEC-RAS恒定流与非恒定流的解法不同,是造成差异的主要成因。
2) 同一计算方法MIKE11与HEC-RAS的计算结果存在较大差异。在同一小流量下,HEC-RAS计算的磨底河水位比MIKE11的高;在同一大流量下,HEC-RAS计算的南河水位比MIKE11的低,与参考文献 [7] 的结果相符。主要原因在于:两种模型定义的水力半径不同,HEC-RAS采用水力半径,MIKE11采用阻力半径。根据参考文献 [5] ,阻力半径适用于断面形状突变的河道,能够避免洪水刚上滩时流量模数减少的情况,而水力半径适用于断面形状狭窄的河道;两种模型的解法不同,HEC-RAS是基于四点隐式差分格式进行计算,MIKE11则采用六点隐式差分法求解。
3) HEC-RAS计算的磨底河水位在下游起推段比MIKE11的水位低,其主要原因是MIKE11受下游起推水位的影响更加敏感。
4. 结论
结合对成都锦江其他河道的水面线推算和对比,总结如下:
1) 恒定流是假设每一计算时间步长内流量是恒定均匀的,不考虑河道调蓄、回水影响等。反之非恒定流考虑了河道调蓄和回水影响,且运动要素随时间变化,能够模拟现实中不同状况下水流的运动状态。建议在工程实践中对采用非恒定流进行计算,计算结果相对更加精准。
2) 锦江干流河道尺度较大,河床存在较多滩地的情况,故MIKE11的阻力半径更加适用;而锦江支流尺度较小,断面深窄,HEC-RAS的水力半径更加适用。
3) 在同一小流量下,HEC-RAS的水位比MIKE11的水位高;在同一大流量下,MIKE11的水位比HEC-RAS的水位高。建议锦江干流如府河、南河选用MIKE11计算水位;锦江支流如磨底河、浣花溪、沙河、桃花江等选用HEC-RAS计算水位。
参考文献