velocity参数重新赋值_二维模型糙率参数自动率定——以海南昌化江入海口段为计算实例...

原标题：

基于BP神经网络的HEC-RAS二维模型糙率参数自动率定

摘 要：

水面宽阔河段、洪泛区和潮汐河口段的水文及水动力学计算模拟多采用二维模型,然而因模型复杂的非线性映射关系常导致参数率定效率和精度不高,甚至出现“异参同效”问题。为此,提出一种应用BP神经网络自动率定二维模型糙率参数的方法。以海南昌化江入海口段为计算实例进行数值检验,结果表明:编写的Python程序能够准确率定HEC-RAS二维模型实例各糙率分区的糙率参数,各糙率分区的率定糙率值均在参数取值范围内,且均在观测糙率值±0.011内,未出现“异参同效”现象;水位率定计算误差满足设定允许误差要求,设定允许误差0.20 m、0.15 m和0.10 m的平均率定计算误差小于0.10 m,设定允许误差0.05 m的平均率定计算误差接近0 m;设定允许误差为0.15 m时,率定程序的运算效率最高。选用海南201409号台风“威马逊”时期的昌化江实测资料进行验证,证明该方法有效可靠,能够实现HEC-RAS软件二维模型糙率率定功能,可推广应用于其他模型参数率定。

关键词：

糙率; 自动率定; BP神经网络; HEC-RAS; Python; 水位; 水力特性; 数值计算;

作者简介：

夏铭辉(1987—),男,工程师,学士,主要从事水利工程设计及技术咨询工作。E-mail:xia_mh@aliyun.com;

*秦景(1984—),男,教授级高级工程师,博士,主要从事水利工程设计、防洪减灾及水文水资源技术咨询和研究工作。E-mail:qinjing@iwhr.com;

基金：

国家重点研发计划(2018YFC1508203);

引用：

夏铭辉，秦景，牛文龙，等． 基于 BP 神经网络的 HEC-ＲAS 二维模型糙率参数自动率定［J］. 水利水电技术，2020，51( 5) : 38-46.

XIA Minghui，Qin Jing，NIU Wenlong，et al． BP neural network-based automatic calibration of roughness parameter for HEC-ＲAS 2-D model［J］. Water Ｒesources and Hydropower Engineering，2020，51( 5) : 38-46.

---

0 引 言

算法调参是数学模型参数率定问题的主要解决方式,各类算法的应用研究较多,特别是在新安江模型的算法调参应用上。

随着数值模拟技术的快速发展,集成多种数学模型的封装数值模型软件开始应用于水文、水动力和水质的计算模拟。其中,丹麦水力学研究所的MIKE系列软件、美国国家环境保护局的SWMM和WASP软件、美国陆军工程兵团水文工程中心HEC系列软件及中国水利水电科学研究院的IFMS/Urban应用较多。然而,这类模型软件因其模型更复杂、参数更多,常导致调参效率和精度不高,甚至出现“异参同效”问题。

为了避免“异参同效”现象和追求更高的率定效率及精度,因封装数值模型软件不易修改源代码,部分学者针对软件和其率定参数的特性选择算法,尝试从外部编程环境下进行参数率定,并取得了较好的效果。GOODELL在Excel中应用VBA编程实现HEC-RAS一维模型的调用、运算和率定。张质明等利用Matlab环境下Simulink工具实现应用GLUE法在WASP软件中进行多指标水质模型参数率定,并成功应用于北运河水质模型的参数率定。郭仪等采用BP神经网络对MIKE SHE模型进行模型参数自动率定,以丹麦Karup流域为例进行数值试验,结果表明该方法比MIKE SHE自带的参数率定程序具有更高的准确性。但目前封装数值模型软件参数率定研究主要集中在一维模型上,二维模型的算法调参相关研究甚少。

HEC-RAS为美国陆军工程兵团水文工程中心编写的封装模型软件,具有出色的计算功能和解算效率,被国内外广泛应用于水文及水动力学分析,特别是洪水模拟计算中。该软件自版本5.0.1起有二维模型计算程序,但缺少糙率率定功能。糙率是二维模型计算中的灵敏参数,实现糙率率定功能才能保证模型模拟的精度。BP神经网络是一种常用的神经网络算法,多用于径流、洪水和水位等水文预报研究。该算法具有较强的非线性映射能力,在二维模型参数率定上有较好的应用前景。

本文提出一种应用BP神经网络自动率定二维模型糙率参数方法,利用BP神经网络算法较强的非线性映射能力提高参数率定的效率和精度,避免出现“异参同效”问题。通过编写Python率定程序实现了该方法在HEC-RAS二维模型中的应用,经数值试验和案例验证,该方法具备一定的率定效率和效果,可为相关二维模型软件的算法调参提供借鉴和指导。

1 方法原理

1.1 HEC-RAS 二维模型基本方程

本文计算采用HEC-RAS二维模型构建水动力模型,模型网格为对偶网格。采用二维圣维南方程组通过半隐式有限体积法计算求解