1、第 54 卷 第 2 期2 0 2 3 年 2 月人民长江YangtzeiverVol 54，No 2Feb，2023收稿日期:2022 04 20基金项目:国家自然科学基金面上项目(51779078);江苏省高校优势学科建设工程项目(PAPD)作者简介:易梓杨，男，硕士研究生，主要从事水环境数值模拟研究。E mail:997994238 qq com通信作者:钱进，男，教授，博士，主要从事水环境系统规划理论与方法、河湖水污染治理与生态修复等方面的研究。E mail:hhuqj hhu edu cn文章编号:1001 4179(2023)02 0141 06引用本文:易梓杨，郦建锋，钱进，等
2、基于 BP 与 PSO 的河口双边界条件逆向推求方法 J 人民长江，2023，54(2):141 146基于 BP 与 PSO 的河口双边界条件逆向推求方法易 梓 杨1，2，郦 建 锋3，4，钱进1，2，陆 卞 和1，2，张 语 航1，2，李 丰 铎1，2(1 河海大学 环境学院，江苏 南京 210098;2 河海大学 浅水湖泊综合治理与资源开发教育部重点实验室，江苏 南京 210098;3 中国电建集团华东勘测设计研究院有限公司，浙江 杭州 311122;4 浙江省华东生态环境工程研究院，浙江 杭州 311122)摘要:在入海河道的防洪调度计算过程中，经常会遇到在确定河口计算区间某一断面水位
3、条件下逆向推求双边界(上游流量、下游潮位)条件的情形。针对基于 MIKE 11 模型的传统试算法逆向推求河口双边界条件时计算效率与搜索效率低的问题，提出采用 BP 神经网络和 PSO 算法来改进传统试算法的优化方法。该方法首先通过MIKE 11 模型建立河口双边界与某一断面水位的离线数据库，然后利用 BP 神经网络建立河口双边界与该断面水位的高精度非线性映射关系，最后以该断面的确定水位为优化目标，采用 PSO 算法逆向推求进而确定双边界条件的映射关系。晋江河口的实例研究表明:该优化方法与传统试算法相比，在保障计算精度的前提下，计算时长约减少至原来的十分之一，大大提高了计算效率。研究成果可为河口
4、地区防洪调度提供参考依据。关键词:河口双边界;逆向推求;MIKE 11 模型;BP 神经网络;PSO 算法;传统试算法;晋江中图法分类号:TP183;P333文献标志码:ADOI:10 16232/j cnki 1001 4179 2023 02 0210引 言台风期的强降雨与天文潮会给中国东部沿海河口地区带来严重的洪涝灾害，造成大量经济损失1 2。为保证河口地区的防洪安全，水利部门经常会采用数学模型来预测河口地区的水文过程。MIKE 11 模型以其简单干净的数据处理界面、可靠的边界条件输入、模拟板块的多样化等特点在水文模拟中被广泛应用3 5。在河口地区水文模拟的过程中，一般是以上游流量和下游
5、潮位为出发点，探究双边界条件对河口地区的影响6 9，而以限制河口计算区间某一断面水位为出发点，逆向推求双边界(上游流量、下游潮位)条件的研究较少。然而在逆向推求河口双边界条件时，基于 MIKE 11 软件的传统试算法是通过手动调节双边界条件来完成该任务，搜索随机性、效率、时间成本成为传统试算法亟待解决的问题。反馈式神经网络模型，即 BP(Back Propagation)神经网络模型是最成熟及应用范围最广的人工神经网络模型，适合解决模拟软件计算效率低的问题10。粒子群优化(Particle Swarm Optimizer，PSO)算法是一种新型的全局随机优化算法，其优点是结构简单、全局搜索能力
6、强、鲁棒性好、易于实现和优化效率高11。BP 神经网络和 PSO 算法以其网络结构优良和搜索效率高在各种复杂的优化问题中得到广泛应用。例如，赵文刚等12 基于 BP 神经网络与 PSO 算法通过相关性分析和因子贡献率分析，建立西洞庭湖南咀水文站月平人民长江2023 年均径流预报模型，实现了在保障拟合精度的前提下确定影响径流预报因子最少个数，相比于传统试算法，提高了计算效率;何胜男等13 以调蓄池截污总量最大为优化目标，以污水处理厂处理规模和可用征地为约束条件，耦合 BP 神经网络与 PSO 算法推求各调蓄池的优化截污率，实现了满足占地面积和污水处理厂处理规模下的截污效率最大化，极大节约了计算时
7、间。本文基于 BP 神经网络与 PSO 算法，提出一种在河口某一断面水位确定条件下逆向推求双边界(上游流量和下游潮位)条件的优化方法，并以晋江河口为例，将该方法的计算结果与传统试算法的计算结果进行对比分析。1传统的河口双边界条件逆向推求过程MIKE 11 水动力模型由于特有的数据输入格式与数据集成平台的限制，只能根据输入边界得到输出结果，为得到河道水位首先需要将上、下游边界的数据修改成MIKE 11 可识别的dfs0 时间序列文件，然后将时间序列文件导入 MIEK 11 的边界文件，最后运行模型得到河口沿程水面线计算结果文件 14 16，计算效率较低。因而，以河口某一断面限制水位为约束条件，采
8、用传统试算法逆向推求双边界条件，首先需确定河口下游高潮位，接着根据下游高潮位对上游最大流量不断地进行调整，并且在 MIKE 11 模型中一一计算，直到求得满足该断面水位条件的上游最大流量，最后输出上游流量、下游潮位双边界条件，具体计算流程如图 1所示。由于下游高潮位是一个变化的区间，因此试算法需要根据下游高潮位的变化反复调整上游最大流量，搜索效率低，灵活性较差。图 1传统试算法逆向推求河口双边界条件流程Fig 1Flow chart of traditional trial algorithm for inversederivation of estuarine double boundary
9、 conditions2基于 BP 神经网络和粒子群优化算法的河口双边界条件逆向推求过程MIKE 11 水动力模型从输入双边界条件到输出某一断面水位结果是一个繁琐的过程，因此采用 BP 神经网络建立输入双边界条件与某一断面水位的非线性映射关系。非线性映射关系建立后，从边界条件输入到水位结果输出的时间大大减少。同时为了改进试算法手动修改边界、手动计算不同工况等重复性工作，引入粒子群优化算法，以某一断面水位为优化目标，根据下游高潮位的变化自动搜索上游最大流量进而提高搜索效率。21BP 神经网络建立输入与输出的非线性映射关系BP 神经网络是人工神经网络的一种，主要由输入层、隐含层和输出层组成，具有多
10、层前馈网络和误差反向传播的学习机制17 18，图 2 为建立的河口双边界与某一断面水位的 BP 神经网络结构图。BP 神经网络具有高度的非线性映射能力和较强的泛化能力，其基于梯度下降法原理，预测精度高19 20。只要隐含层和隐节点足够多，可以以任意精度逼近一个非线性映射关系21。图 2BP 神经网络模型结构Fig 2Diagram structure of BP neural network modelBP 神经网络结构的拟合精度受样本数量和样本代表性的影响。样本数量根据实际情况确定。为提高样本数据的代表性，对双边界条件采用拉丁超立方抽样的方法，防止抽取数据过于集中22;接着将抽取的样本用 M
11、IKE batch 模块并行计算;最后利用 MATLAB提取指定断面的水位结果，进而得到双边界条件和指定断面水位的离线数据库。BP 神经网络的训练函数选取的是 L M(Leven-berg Marquardt)算法的 training 函数，设置 BP 神经网络学习率，训练误差为。将离线数据库分为训练集和验证集。BP 神经网络结构的建立分为两部分:BP 神经网络的训练。通过 BP 神经网络各层之间的权值与阈值，不断调整训练集中预测输出与期望输出的偏差，直到误差小于预设值，训练结束。BP 神经网络的验证。将验证集的双边界条件带入 BP 神经网络得到预测输出，对比预测输出与期望输出的差异以评价网络
12、结构的优良性。采取 MSE(均方误差)、2(相关系数)来评判 BP 神经网络的拟合效果。241第 2 期易梓杨，等:基于 BP 与 PSO 的河口双边界条件逆向推求方法2 2粒子群优化算法实现河口双边界条件的逆向推求根据 BP 神经网络建立的河口双边界条件与指定断面水位的映射关系，引入粒子群优化算法，实现对边界条件的自动搜寻，提高搜索效率23 24。PSO 是科学家研究鸟类捕食行为而得出来的一种算法结构，其原理为:PSO 初始化为一组随机粒子(随机解)，在每次搜寻过程中粒子通过两个极值(Pbest，Gbest)来更新自己的位置25，各个粒子根据对比自身最优解与群体最优解来调整自己的位置和速度，
13、直至到达预设的精度或者预设的算法迭代次数26。通过对 PSO 变量中初始种群数、学习因子、迭代次数以及对河口某一断面水位(优化目标)的设定，构建 PSO 算法结构。设置河口某一断面水位为优化目标，初始化的粒子随机分布在样本空间中，粒子不断调整自身空间位置以搜索满足优化目标的最优解。2 3优化方法与传统试算法对比图 3 为优化方法与传统试算法的计算流程对比图。优化方法相比于传统试算法，主要从两个方面进行了改进:基于 BP 神经网络构建河口双边界条件与某一断面水位的非线性映射关系，从边界条件的输入到断面水位结果的输出只需要几秒，提高了传统试算法的计算效率;基于粒子群优化算法，通过优化目标的设定，实
14、现双边界条件的自动搜索，提高了传统试算法的搜索效率。图 3优化法与传统试算法计算流程对比Fig 3Comparison of the calculation processes of the optimizationmethod and the traditional trial method3实例应用3 1研究区域概况晋江位于福建省东南沿海，是中国台风登陆的主要路线之一，台风活动极为频繁，多发生在 7 9 月。如遭遇天文大潮，其危害性更大。晋江作为主要的泄洪通道，在雨水频繁期，承载着巨大的泄洪压力。查阅晋江水文资料，金鸡水闸 100 a 一遇的洪峰流量为11 100 m3/s;前浦站100
15、a 一遇的最高潮位为5 16 m，多年平均最高潮位为 4 00 m。本次研究范围从金鸡水闸到前浦站，共 19 30 km，具体如图 4 所示。为了保证晋江河口两岸人民的安全，需要限制防洪安全空间不足地区的最高水位。根据福建省水利水电勘测设计研究院2020 年编制的 晋江下游防洪标准复核报告 可知，晋江河口地区民俗主题文化公园处的防洪措施最为薄弱。在 100 a 一遇洪水和 100 a 一遇潮位的情形下，民俗主题文化公园处计算水位为7 832 m。为了给民俗主题文化公园预留足够的防洪空间，本文以民俗主题文化公园水位 6 832 m 为优化目标，对双边界条件进行逆向推求。图 4研究区域概况Fig
16、4Overview of the research area3 2水动力模型构建以 MIKE 11 为计算平台，对研究区域建立水动力模型，概化河长为 19 30 km，从金鸡水闸到前浦站，共概化断面 55 个。为保证模型计算时的稳定性，根据河道水位的实测值，设置初始水位为 2 00 m。泉州大桥(桥下)水位站为国家基本水文站，有逐日水位观测数据。晋江下游沿程设有水尺，自动化记录水位数据。2016 年 9 月 15 日为晋江下游近年来最大洪水，2019 年 6 月 3 日为当年最大洪水，根据水文站及水尺的记录数据，采用 2016 年 9 月 14 16 日测量水位数据对建立的水动力模型进行率定，选择 2019341人民长江2023 年年 6 月 2 5 日测量水位数据对率定的水动力模型进行验证。模型率定的结果表明:河道主槽糙率在0 036 0 040 之间，滩地糙率在 0 042 0 050 之间。模型验证的结果表明:河道沿程水位的实测值与计算值拟合较为良好，水位绝对误差小于5 cm(见表1);泉州大桥(桥下)水位站实测水位与计算水位的绝对误差小于 8 cm(见图 5)。水利部颁布的 S