CN113887787A - 一种基于长短时记忆网络和nsga-ii算法的洪水预报模型参数多目标优化方法 - Google Patents

Abstract

一种基于长短时记忆网络和NSGA‑II算法的洪水预报模型多目标优化方法，首先，确定长短时记忆LSTM网络洪水预报模型的结构，将收集与整理的研究流域场次洪水资料划分为训练集、测试集场次，截取每次洪水的前期降雨序列使每个样本的输入降雨序列长度相等，得到训练集、测试集样本的输入降雨序列。其次，根据防洪保护对象的实际需求设计多个目标函数，构建LSTM洪水预报模型参数多目标优化框架，将LSTM洪水预报模型的前向计算嵌套至多目标优化框架内，循环更新与评价模型参数组，得到最优参数组。最后，将最优参数组输入LSTM洪水预报模型，并评估分析LSTM洪水预报模型的模拟、预报效果。本发明能够满足不同场景下LSTM洪水预报模型参数多目标优化的需求，为山洪灾害预报预警工作提供了新的技术支撑。

Description

一种基于长短时记忆网络和NSGA-II算法的洪水预报模型参 数多目标优化方法

技术领域

本发明属于洪水预报技术领域，涉及一种基于长短时记忆网络和NSGA-Ⅱ算法的洪水预报模型参数多目标优化方法。

背景技术

随着进化算法的发展和高性能计算革命的兴起，机器学习以其高效的数据处理能力和强大的泛化能力引起了广大水文工作者的注意，机器学习模型与洪水预报的交叉研究已成为当下国内外水文领域的研究热点之一。长短时记忆(LSTM)网络作为一种特殊的循环神经网络(RNN)，借助其内部特殊控制门和细胞单元状态能够学习输入序列数据的长期依赖关系，有效避免了梯度消失的问题，目前已在语音识别、机器翻译、计算机视觉等领域取得成功应用，也逐渐被应用于洪水预报领域。但受限于TensorFlow、Keras和Pytorch等主流深度学习框架下机器学习模型构建、训练和测试的计算模式，LSTM洪水预报模型的参数优化过程只能设置单一目标(损失)函数，无法满足LSTM洪水预报模型多目标参数优化的实际需求，从而难以反映流域水文系统的不同动态行为特征，且无法兼顾流域内不同防洪保护对象对洪水预报模型输出值的需求。因此，要实现LSTM洪水预报模型参数的多目标优化必须摆脱现有深度学习框架的限制，重新构建LSTM洪水预报模型参数多目标优化框架。

带精英策略的快速非支配排序遗传算法(NSGA-II)是传统水文模型参数优化中应用较为广泛的多目标优化算法，该算法最大的特点是适用于低维多目标优化问题，降低了非劣排序遗传算法的复杂性，通过引入精英策略提高优化结果的精度，具有运行速度快，解集收敛性好的优点。为此，本发明立足于提高山丘区中小流域洪峰值的预报精度，提出一种新的模型参数多目标优化方法。以NSGA-II作为洪水预报模型的多目标优化算法，构建LSTM洪水预报模型参数的多目标优化框架，将LSTM洪水预报模型计算的多个目标函数值反馈给优化算法来调整与更新洪水预报模型参数组，实现基于长短时记忆网络和NSGA-II算法的洪水预报模型参数多目标优化，以提高洪水预报模型的峰值模拟和预报性能。

发明内容

本发明针对现有技术存在的问题，提供了一种基于长短时记忆网络和NSGA-II算法的洪水预报模型参数多目标优化方法。

为了达到上述目的，本发明采用的技术方案如下：

一种基于长短时记忆网络和NSGA-II算法的洪水预报模型参数多目标优化方法，包括以下步骤：

第一步，确定长短时记忆(LSTM)网络洪水预报模型的结构。

1.1)设计模型输入输出结构：模型输入为流域内各雨量站长序列实测降雨，输入长度为n个时段，每次洪水输入的长度n均相同。模型输出为与场次洪水历时相等的流量值序列，输出长度为l个时段，l≤n；

1.2)给定长短时记忆(LSTM)网络洪水预报模型的隐藏层层数和隐藏层神经元节点数量。

第二步，生成样本。

收集与整理研究流域场次洪水资料，并将所有场次洪水划分为训练集、测试集场次。其中，训练集场次洪水用于优化LSTM洪水预报模型的参数；测试集场次洪水用于检验训练后模型的外推能力。根据第一步设计的模型输入输出结构确定样本长度，每次洪水对应一个样本，截取每次洪水的前期降雨序列使每个样本的输入降雨序列长度等于n，获得输入降雨序列P＝[P₁,P₂,…,P_t,…,P_n]。其中，n为LSTM洪水预报模型输入降雨序列长度(h)，每一时间步的输入降雨P_t包含流域内各雨量站第t个时段的降雨量。

第三步，设计多目标函数。

不同的防洪保护对象在洪水预报中具有不同的关注重点。山丘区中小流域的洪水预报重点关注流域出口或河道控制断面的峰值是否达到预警流量，相应的洪水预报模型应侧重考虑洪峰相对误差，尤其大洪水的洪峰误差；入库洪水预报中入库洪量和洪水过程是水库管理人员关注的重点，洪水预报模型既要考虑产流量(次洪量)预报精度，又要关注洪峰流量大小和场次洪水过程；若以河道防洪为主，则应侧重考虑河道控制断面的洪峰流量与峰现时间预报结果。因此，洪水预报模型需根据防洪保护对象的实际需求设置多个目标函数进行参数优化。流域洪水预报常用的单目标函数包括：场次洪水峰值合格率、峰现时间合格率、纳什效率系数合格率、均方误差和加权均方误差等。

第四步，LSTM洪水预报模型的前向计算。主要步骤如下：

4.1)输入训练集样本数据，给定一组LSTM洪水预报模型的控制门、细胞单元状态的权重矩阵和偏置向量参数W_f、U_f、b_f、W_c、U_c、b_c、W_i、U_i、b_i、W_o、U_o、b_o以及输出层权重矩阵和参数W_out、b_out。其中，W_f、U_f、b_f分别为遗忘门输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_c、U_c、b_c分别为细胞单元更新状态中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_i、U_i、b_i分别为输入门中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_o、U_o、b_o分别为输出门中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_out、b_out分别为输出层的权重矩阵和偏置向量；

4.2)初始化模型内部状态，令输入门、遗忘门、输出门、细胞单元状态和隐藏状态的初始状态i₀、f₀、o₀、h₀、C₀均等于零向量，向量长度等于神经元节点数量；

4.3)逐时段循环计算，第t时段的计算涉及遗忘门的门控状态f_t、细胞单元潜在更新状态

输入门的门控状态i_t、输出门的门控状态o_t、细胞单元状态C_t和隐藏状态h_t，其中t＝1,2,…,n；

4.4)提取LSTM洪水预报模型的输出隐藏状态序列h＝[h₁,h₂,h₃,…,h_n]；

4.5)计算LSTM洪水预报模型的输出流量值序列Q＝[Q₁,Q₂,Q₃,…,Q_n]；

4.6)提取每个样本对应的模型输出流量值序列Q＝[Q_n-l+1,Q_n-l+2,Q_n-l+3,…,Q_n]，l为样本的输出长度，其值大小等于场次洪水历时。

第五步，构建LSTM洪水预报模型参数多目标优化框架，求解多目标优化问题。

基于Platypus库调用NSGA-II算法，通过设置算法的种群规模、评价次数、二进制交叉算子(SBX)及多项式变异概率(PM)等参数取值，搭建NSGA-II多目标优化框架。将第四步中LSTM洪水预报模型的前向计算嵌套至多目标优化框架内，以第三步设计的多个目标函数作为LSTM洪水预报模型参数优化的目标函数，更新一组LSTM洪水预报模型参数(权重矩阵和偏置向量)，并提取LSTM洪水预报模型前向计算中每个样本对应的模型输出流量值序列，在此基础上计算出多个目标函数值来评价当前这组模型参数的优劣。以NSGA-Ⅱ算法的评价次数作为循环次数，重复上述参数更新与评价过程，直至多个目标函数结果收敛，使模型总体的拟合性能达到最优，从而得到最优参数组。

第六步，将第五步优选得到的参数输入LSTM洪水预报模型，利用训练集、测试集样本数据评估分析多目标函数优化的LSTM洪水预报模型的模拟、预报效果。

本发明提出了一种基于长短时记忆网络和NSGA-II算法的洪水预报模型参数多目标优化方法，将LSTM洪水预报模型的前向计算嵌套至NSGA-II多目标优化框架内，实现LSTM洪水预报模型参数多目标优化，提升机器学习洪水预报模型的性能。

上述基于长短时记忆网络和NSGA-II算法的多目标优化方法应用于流域洪水预报。

本发明的有益效果为：基于主流的深度学习框架下构建长短时记忆(LSTM)网络洪水预报模型无法实现模型参数多目标优化，本发明将LSTM洪水预报模型的前向计算嵌套至多目标优化框架内，利用NSGA-II算法求解多目标优化问题，构建了LSTM洪水预报模型参数多目标优化框架，能够满足不同场景下模型参数多目标优化的实际需求；同时针对山丘区中小流域洪水预报侧重于洪峰预报的实际需求，开展了LSTM洪水预报模型参数多目标优化示范应用，有效提高了山丘区中小流域场次洪水峰值模拟、预报精度，为山洪灾害预报预警工作提供了新的技术支撑。

附图说明

图1是本发明实例应用所采用的安和流域图；

图2是本发明生成样本示意图；

图3是本发明LSTM洪水预报模型参数多目标优化框架图；

图4是本发明LSTM洪水预报模型参数优选计算过程示意图；

图5是本发明多目标函数方案下LSTM洪水预报模型场次洪水洪峰流量模拟与预报结果图，其中图5(a)为训练集场次洪水洪峰流量模拟与预报结果图，图5(b)为测试集场次洪水洪峰流量模拟与预报结果图；

图6是本发明多目标函数方案下LSTM洪水预报模型典型场次洪水的模拟、预报流量过程图，其中图(a)为19970620场次洪水模拟结果，图(b)为20040705次洪水预报结果，图(c)为20010609场次洪水模拟结果，图(d)为20050619场次洪水预报结果，图(e)为20010801场次洪水模拟结果；图(f)为20060505场次洪水预报结果；图(g)为20020728场次洪水模拟结果；图(h)为20060725场次洪水预报结果。

具体实施方式

以下结合具体实施例对本发明作进一步说明。

本发明提出了一种基于长短时记忆网络和NSGA-II算法的洪水预报模型参数多目标优化方法。研究流域实例应用中训练集、测试集场次洪水的结果分别代表了LSTM洪水预报模型的模拟、预报性能。下面通过实施案例，并结合附图，对本发明做进一步说明。

安和流域位于江西省赣州市上犹县，地处东经114°-114°40′，北纬25°42′-26°01′，流域面积为251km²。该流域内植被发育良好，属于亚热带丘陵山区湿润季风气候区，雨量充沛，年均降雨量1497mm，平均气温18.8℃，安和流域地形及水文测站分布见图1。选取该山丘区中小流域作为研究实例进行洪水预报，实现基于长短时记忆网络和NSGA-II算法的洪水预报模型多目标优化。主要步骤如下：

第一步，确定长短时记忆(LSTM)网络洪水预报模型的结构。

1.1)将流域内各雨量站长序列实测降雨作为模型输入，输入长度为486h；输出流量值序列长度等于场次洪水历时；

1.2)LSTM洪水预报模型的隐藏层为单层，隐藏层神经元节点数量为5个。

第二步，生成样本。

收集与整理得到安和流域86场实测洪水过程。其中，61场为训练集场次，用于优化LSTM洪水预报模型的参数；25场为测试集场次，用于检验训练后模型的外推能力。截取每次洪水的前期降雨序列使每个样本的输入降雨序列长度等于486h，获得输入降雨序列P＝[P₁,P₂,…,P_t,…,P₄₈₆]。其中，每一时间步的输入降雨P_t包含流域内8个雨量站第t个时段的降雨量。以场次洪水19940613为例，选取的计算样本如图2所示。

第三步，设计多目标函数。

针对山丘区中小流域洪水预报侧重于场次洪水峰值预报的实际需求，设置组合加权均方误差(WMSE)作为优化目标函数之一，将每次洪水过程中的流量值划分为高、中、低量级流量，划分阈值按照场次洪水实测洪峰值的2/3、1/3进行设置。其中，大于洪峰值2/3的为高流量，小于洪峰值1/3的为低流量，其余为中流量。MSE公式见式(1)，WMSE计算公式见式(2)。

WMSE＝k₁HMSE+k₂MMSE+k₃LMSE (2)

式中：Q_ij,sim、Q_ij,obs分别为第i场洪水第j时段的预报流量值和实测流量值，单位为m³/s；HMSE、MMSE和LMSE分别为高流量、中流量和低流量对应的均方误差MSE计算值；k₁、k₂和k₃为权重系数，三者取值大小需符合k₁＞k₂＞k₃的约束。采用试错法选取权重系数的若干组合，最终根据模型在训练集场次洪水中的表现确定取值，其中k₁＝10.0，k₂＝1.0，k₃＝0.1。

设置两种多目标函数优化方案，并与TensorFlow深度学习框架下单目标函数MSE方案进行对比：

1)QRP-QRT-NSE：以场次洪水的峰值合格率QRP、峰现时间合格率QRT和纳什效率系数NSE合格率作为目标函数，三者均越大越优；

2)WMSE-QRP：以加权均方误差WMSE为目标函数，同时兼顾场次洪水峰值合格率QRP。加权均方误差WMSE越小越优，场次洪水峰值合格率QRP越大越优；

3)TensorFlow-MSE：深度学习框架TensorFlow下以均方误差MSE为目标函数。

其中，洪峰流量相对误差≤20％视为合格；峰现时间绝对误差≤2h视为合格；场次洪水模拟、实测过程纳什效率系数NSE大于0.7视为合格。峰值合格率QRP计算见式(3)，峰现时间合格率QRT计算见式(4)，纳什效率系数NSE计算见式(5)。

式中：N_P表示场次洪水洪峰流量合格数量，N_T表示场次洪水峰现时间合格数量，N表示场次洪水总数；Q_ij,sim、Q_ij,obs分别为第i场洪水第j时段的预报流量值和实测流量值，单位为m³/s；

为场次洪水平均实测流量，单位为m³/s。

第四步，LSTM洪水预报模型的前向计算。主要步骤如下：

4.1)输入训练集样本数据，给定一组LSTM洪水预报模型的控制门、细胞单元状态的权重矩阵和偏置向量参数W_f、U_f、b_f、W_c、U_c、b_c、W_i、U_i、b_i、W_o、U_o、b_o以及输出层权重矩阵和参数W_out、b_out。其中，W_f、U_f、b_f分别为遗忘门输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_c、U_c、b_c分别为细胞单元更新状态中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_i、U_i、b_i分别为输入门中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_o、U_o、b_o分别为输出门中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_out、b_out分别为输出层的权重矩阵和偏置向量；

4.2)初始化模型内部状态，令输入门、遗忘门、输出门、细胞单元状态和隐藏状态的初始状态i₀、f₀、o₀、h₀、C₀均等于零向量，向量长度等于神经元节点数量；

4.3)逐时段循环计算，第t时段的计算涉及遗忘门的门控状态f_t、细胞单元潜在更新状态

输入门的门控状态i_t、输出门的门控状态o_t、细胞单元状态C_t和隐藏状态h_t，其中t＝1,2,…,486；

4.4)提取LSTM洪水预报模型的输出隐藏状态序列h＝[h₁,h₂,h₃,…,h₄₈₆]；

4.5)计算LSTM洪水预报模型的输出流量值序列Q＝[Q₁,Q₂,Q₃,…,Q₄₈₆]；

4.6)提取每个样本对应的模型输出流量值序列Q＝[Q_486-l+1,Q_486-l+2,Q_486-l+3,…,Q₄₈₆]，l为样本的输出长度，其值大小等于场次洪水历时。

第五步，构建LSTM洪水预报模型参数多目标优化框架，求解多目标优化问题。

基于Platypus库调用NSGA-II算法，将算法中种群规模设置为100，评价次数为5万，二进制交叉算子(SBX)取值为(1.0，15.0)，多项式变异概率(PM)取值为(0.125，20.0)，搭建NSGA-II多目标优化框架，见图3。将第四步中LSTM洪水预报模型的前向计算嵌套至多目标优化框架内，以第三步设计的多个目标函数作为LSTM洪水预报模型参数优化的目标函数，更新一组LSTM洪水预报模型参数(权重矩阵和偏置向量)，并提取LSTM洪水预报模型前向计算中每个样本对应的模型输出流量值序列，在此基础上计算出多个目标函数值来评价这一组模型参数的优劣。以NSGA-Ⅱ算法的评价次数作为循环次数，重复上述参数更新与评价过程，直至多个目标函数结果收敛，使模型总体的拟合性能达到最优，得到最优参数组。多目标优化框架下LSTM洪水预报模型参数优选计算过程见图4。

第六步，将第五步优选的参数输入LSTM洪水预报模型，对比分析单目标及两种多目标函数方案下的LSTM洪水预报模型的场次洪水模拟、预报结果。

评价指标采用场次洪水的峰值合格率QRP、峰现时间合格率QRT、纳什效率系数NSE、平均绝对误差MAE以及平均均方根误差RMSE。MAE计算见式(6)，RMSE计算见式(7)。

式中：Q_ij,sim、Q_ij,obs分别为第i场洪水第j时段的预报流量值和实测流量值，单位为m³/s。

安和流域多目标函数方案LSTM洪水预报模型场次洪水模拟、预报统计结果见表1。就山丘区中小流域重点关注的场次洪水洪峰流量而言，多目标函数方案优化得到的LSTM洪水预报模型对场次洪水的峰值模拟、预报能力整体上要优于基于TensorFlow框架训练的LSTM洪水预报模型，尤其场次洪水的峰值合格率QRP、峰现时间合格率QRT明显优于TensorFlow深度学习框架下MSE损失函数取得的结果。此外，以场次洪水峰值合格率QRP、加权均方差WMSE组成的多目标函数方案优化的LSTM洪水预报模型的洪峰合格率最高，安和流域场次洪水峰值模拟、预报合格率分别可达到0.902和0.88，说明加大高流量的权重能显著提升模型对洪峰值的拟合与预报能力。

表1多目标函数方案下LSTM洪水预报模型洪水预报结果

图5为上述两种多目标函数方案下训练集、测试集场次洪水洪峰流量模拟、预报值与实测值对比图，并给出了计算值与实测值的相关性系数。由图5明显可以看出，以峰值合格率QRP、加权均方差WMSE组成的多目标函数方案优化的LSTM洪水预报模型对场次洪水的峰值模拟、预报效果均更好，训练集、测试集场次洪水模拟、预报洪峰值与实测值相关性系数分别0.946、0.931，且该方案下LSTM洪水预报模型对较大洪峰流量值(大于200m³/s)的模拟、预报结果均更优。

表2为多目标函数方案优化的LSTM洪水预报模型的典型场次洪水的模拟、预报统计结果。图6为多目标函数方案优化的LSTM洪水预报模型的典型场次洪水模拟、预报洪水过程图。由表2和图6可明显看出，以峰值合格率QRP、加权均方差WMSE组成的多目标函数(QRP-WMSE)方案优化的LSTM洪水预报模型的洪峰相对误差较其余两种目标函数方案更小，训练集、测试集场次洪水的模拟、预报洪峰值更接近于实测洪峰值，且低流量阶段的场次洪水模拟、预报流量过程线稍高于其余两种目标函数方案，说明通过多目标优化方法得到的LSTM洪水预报模型对山丘区中小流域场次洪水高流量值的预报能力更强。

表2不同目标函数下LSTM洪水预报模型典型场次洪水模拟、预报结果

以上结果表明本发明提出的基于长短时记忆网络和NSGA-II算法的洪水预报模型多目标优化方法能够有效提高LSTM洪水预报模型在山丘区中小流域的场次洪水峰值和峰现时间模拟、预报精度。

以上所述实施例仅表达本发明的实施方式，但并不能因此而理解为对本发明专利的范围的限制，应当指出，对于本领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些均属于本发明的保护范围。

Claims (1)

1.一种基于长短时记忆网络和NSGA-II算法的洪水预报模型参数多目标优化方法，其特征在于，所述洪水预报模型参数多目标优化方法能够应用于流域洪水预报，包括以下步骤：

第一步，确定长短时记忆LSTM网络洪水预报模型的结构；

1.1)设计模型输入输出结构：模型输入为流域内各雨量站长序列实测降雨，输入长度为n个时段，每次洪水输入的长度n均相同；模型输出为与场次洪水历时相等的流量值序列，输出长度为l个时段，l≤n；

1.2)给定长短时记忆LSTM网络洪水预报模型的隐藏层层数和隐藏层神经元节点数量；

第二步，生成样本；

收集与整理研究流域场次洪水资料，并将所有场次洪水划分为训练集、测试集场次；其中，训练集场次洪水用于优化LSTM洪水预报模型的参数；测试集场次洪水用于检验训练后模型的外推能力；根据第一步设计的模型输入输出结构确定样本长度，每次洪水对应一个样本，截取每次洪水的前期降雨序列使每个样本的输入降雨序列长度等于n，获得输入降雨序列P＝[P₁,P₂,…,P_t,…,P_n]；其中，n为LSTM洪水预报模型输入降雨序列长度(h)，每一时间步的输入降雨P_t包含流域内各雨量站第t个时段的降雨量；

第三步，洪水预报模型根据防洪保护对象的实际需求设置多个目标函数进行参数优化；

第四步，LSTM洪水预报模型的前向计算；主要步骤如下：

4.1)输入训练集样本数据，给定一组LSTM洪水预报模型的控制门、细胞单元状态的权重矩阵和偏置向量参数W_f、U_f、b_f、W_c、U_c、b_c、W_i、U_i、b_i、W_o、U_o、b_o以及输出层权重矩阵和参数W_out、b_out；其中，W_f、U_f、b_f分别为遗忘门输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_c、U_c、b_c分别为细胞单元更新状态中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_i、U_i、b_i分别为输入门中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_o、U_o、b_o分别为输出门中输入特征的权重矩阵、隐藏状态的权重矩阵和偏置向量，W_out、b_out分别为输出层的权重矩阵和偏置向量；

4.2)初始化模型内部状态，令输入门、遗忘门、输出门、细胞单元状态和隐藏状态的初始状态i₀、f₀、o₀、h₀、C₀均等于零向量，向量长度等于神经元节点数量；

4.3)逐时段循环计算，第t时段的计算涉及遗忘门的门控状态f_t、细胞单元潜在更新状态

输入门的门控状态i_t、输出门的门控状态o_t、细胞单元状态C_t和隐藏状态h_t，其中t＝1,2,…,n；

4.4)提取LSTM洪水预报模型的输出隐藏状态序列h＝[h₁,h₂,h₃,…,h_n]；

4.5)计算LSTM洪水预报模型的输出流量值序列Q＝[Q₁,Q₂,Q₃,…,Q_n]；

4.6)提取每个样本对应的模型输出流量值序列Q＝[Q_n-l+1,Q_n-l+2,Q_n-l+3,…,Q_n]，l为样本的输出长度，其值大小等于场次洪水历时；

第五步，构建LSTM洪水预报模型参数多目标优化框架，求解多目标优化问题；

基于Platypus库调用NSGA-II算法，通过设置算法的种群规模、评价次数、二进制交叉算子SBX及多项式变异概率PM参数取值，搭建NSGA-II多目标优化框架；

将第四步中LSTM洪水预报模型的前向计算嵌套至多目标优化框架内，以第三步设计的多个目标函数作为LSTM洪水预报模型参数优化的目标函数，更新一组LSTM洪水预报模型的权重矩阵和偏置向量，并提取LSTM洪水预报模型前向计算中每个样本对应的模型输出流量值序列，在此基础上计算出多个目标函数值来评价当前这组模型参数的优劣；

以NSGA-Ⅱ算法的评价次数作为循环次数，重复上述参数更新与评价过程，直至多个目标函数结果收敛，使模型总体的拟合性能达到最优，从而得到最优参数组；

第六步，将第五步优选得到的参数输入LSTM洪水预报模型，利用训练集、测试集样本数据评估分析多目标函数优化的LSTM洪水预报模型的模拟、预报效果。

Images