VMD-PSO-LSTM模型的日径流多步预测.pdf

DOI:10.12170/20220208001王秀杰，王玲，滕振敏，等.VMD-PSO-LSTM 模型的日径流多步预测 J.水利水运工程学报，2023（4）：81-90.（WANG Xiujie,WANG Ling,TENG Zhenmin,et al.Research on multi-step forecast of daily runoff based on VMD-PSO-LSTM modelJ.Hydro-Science and Engineering,2023（4）:81-90.（in Chinese）VMD-PSO-LSTM 模型的日径流多步预测王秀杰1，王 玲1，滕振敏2，田福昌1，袁佩贤3，苑希民1(1.天津大学 水利工程仿真与安全国家重点实验室，天津 300072；2.崇左市左江治旱工程管理中心，广西 崇左532200；3.北京城建设计发展集团股份有限公司，北京 100032)摘要:为了弱化径流时间序列的非线性和非平稳性，提高不同预见期的日径流预测精度，提出了一种新的VMD-PSO-LSTM 多步预测组合模型。首先采用变分模态分解（VMD）方法将原始日径流序列分解为子序列，通过粒子群优化算法（PSO）对长短期记忆（LSTM）模型参数进行优化，对各子序列建立 PSO-LSTM 模型，各分量的预测值重构集成预测结果。将 VMD-PSO-LSTM 模型应用于黄河下游花园口和利津站的日径流多步预测，采用 Nash sutcliffe 效率系数（ENS）、相关系数（R）和均方根误差（ERMS）3 个定量评价指标对模型预测结果进行评价。结果表明：在预见期为 1、2、3 d 的情况下，两个测站的 Nash sutcliffe 效率系数和相关系数均在 0.90 以上。与 CEEMD-PSO-LSTM 和 PSO-LSTM 模型的预测结果对比表明，该模型能够有效提高日径流多步预测精度，是一种高效稳定的径流预报模型。关键词：非平稳序列；日径流多步预测；长短期记忆模型；变分模态分解方法中图分类号：TV124 文献标志码：A 文章编号：1009-640X(2023)04-0081-10 准确可靠的径流预测对合理利用水资源、提升防洪减灾能力及充分发挥水库的综合效益具有重要意义。但受气候变化和人类活动等因素的综合影响，径流变化具有较强的非线性和不确定性，为水文预报带来困难。因此，提高水文预报精度一直是水文水资源领域的研究重点和难点。目前，常用的基于数据驱动的径流预测方法有线性预测模型和机器学习模型。线性预测模型基于线性回归理论，结构简单，难以应对复杂非平稳的径流信息1。机器学习的发展为复杂径流的预测提供了新思路，其中传统浅层的机器学习方法虽然能够成功进行复杂非平稳的径流预测，但其精度仍需进一步改进2。随着人工智能技术的不断发展，深度学习被引入预测领域。其中长短期记忆（LSTM）模型具有较强的非线性预测能力、较快的收敛速度和捕捉时间序列长期相关性的能力，是一种具有特殊门控设计的循环神经网络（RNN），可以长时间保留历史有用信息3，已有研究4-5证明了其在径流预测中的优势。虽然机器学习模型在水文预报领域已经取得了丰硕成果，但传统单一的神经网络模型学习过程易产生局部最优和过拟合问题6。为解决此问题，一些学者采用自适应调节人工蚁群（ARACS）、灰狼优化算法（GWO）、PSO 等算法优化神经网络参数，如白继中等7提出的自适应调节人工蚁群-径向基函数神经网络混合预测模型(ARACS-RBF)，径流预测效果明显优于RBF 模型；王立辉等 8提出灰狼优化的长短期记忆模型（GWO-LSTM），将其应用于径流预测，结果表明泛化性能较好；郑恺原等9提出的粒子群算法-加权随机森林模型（PSO-WRF），与 WRF 模型相比，PSO-收稿日期：2022-02-08基金项目：国家重点研发计划资助项目（2018YFC1508403）；科技部重点领域创新团队项目（2014RA031）；国家自然科学基金委创新团队项目（51621092）作者简介：王秀杰（1973），女，辽宁葫芦岛人，副教授，博士，主要从事水文系统预测与模拟研究。E-mail： 通信作者：田福昌（E-mail：）第 4 期水利水运工程学报No.42023 年 8 月HYDRO-SCIENCE AND ENGINEERINGAug.2023WRF 模型优化效果显著，具有良好的预测精度和泛化能力。由于径流序列高度的非线性与非稳定性，使得运用单一模型对径流序列进行预测，难以完全捕捉径流序列中的非线性因素。利用适当的时频信号分解技术分解原始水文时间序列建立集成模型，可有效提取水文时间序列中的有效信息，提高时间序列的平稳性和可预报性，从而提高模型的预报精度。许多学者对“分解-预测-重构”的组合预测模型进行了大量研究。周婷等10提出了小波分解-支持向量机-PSO 算法（WD-SVM-PSO）混合预测模型，基于淮河流域响洪甸水库 19592014 年径流过程进行年径流预测，验证了时间序列分解预处理的有效性；孙娜等11提出了一种将 EEMD、样本熵（SE）和正则化极限学习机（RELM）相结合的非平稳日径流预测方法（ES-RELM），基于金沙江下游控制站屏山站 20002010 年的日径流实测数据进行日径流预测，表明分解后的模型预测精度更高；王文川等12采用具有自适应序列特征的时变滤波经验模态分解（TVF-EMD）与 LSTM 耦合，构成 TVF-EMD-LSTM 预测模型进行洛河流域长水水文站月径流预测；黄景光等13提出一种小波分析-支持向量机（WD-SVM）特征分类组合预测模型，基于宜昌站 20132017 年日径流序列进行预测，结果表明相对于单一模型，新模型预测精度和稳定性明显提高；He 等14建立了变分模态分解(VMD)-深度神经网络（DNN）混合预测模型预报水库径流，预测性能优于 EMD-DNN 和EEMD-DNN 等模型。小波分解需要根据经验预先设定母波函数及分解层数，受人为影响大；EMD、EEMD则可根据时间序列的局部特征自适应将序列分解为一系列固有模态分量，与 WT 相比更具有普适性，但是容易存在伪分量和模态混叠问题，制约了预测精度的提高；VMD 是一种自适应维纳滤波器组，相比 EMD 的递归筛选模式，VMD 是一种完全非递归的变分模态分解，能减少伪分量和模态混叠，并具有比 EMD 和EEMD 更好的噪声鲁棒性15，是分解复杂非线性非平稳时间序列的一种有效手段。从以上研究可以看出，目前径流预测研究多是单步预测，利用“分解-预测-重构”方法对日径流多步预测的研究并不多见，而考虑不同预见期的径流预测模型能够对流域的防洪减灾起到更好的预防警示作用。常用的多步预测方法包括递推多步预测和直接多步预测。递归多步预测通过单步预测模型多次迭代来实现，将存在误差的单步预测输出作为下一步预测的输入，因此，随着多步预测的推进，预测误差不断放大，且在滞后较大的动态非线性系统中误差更为显著。而直接多步预测，根据模型直接得到多步预测值，不存在误差传播的问题，预测精度较高。如 Atiya 等16对比了递推多步预测和直接多步预测的预测效果，结果表明直接多步预测比递推多步预测的预测效果更好、精度更高。综合以上各研究方法的优势，本文将 VMD 分解方法与 LSTM 结合，并基于 PSO 算法优化 LSTM 模型参数，进行不同预见期下的日径流直接多步预测研究。1 研究方法 1.11.1变分模态分解（VMD）VMD 是一种自适应、非递归的变分模态信号分解方法，其自适应性体现在可以根据输入信号的特性确定所给序列的模态分解个数 k，并匹配每种模态的最佳中心频率和有限带宽，以便平衡处理序列各部分存在的噪声17。VMD 的目标是将一个原始输入信号分解成一组离散的具有特殊稀疏性的子信号（模态），在保证分解序列各模态的带宽之和最小的前提下，所有模态之和与原始输入信号相等，对模态函数进行希尔伯特（Hilbert）变换并进行变分约束，表达式如下：minuk,wkk?t(t)+jt)uk(t)ejwkt?22使kuk=f(1)(t)式中：uk和wk分别为原信号分解得到的所有模态及其中心频率的集合；t 为时间变量；为狄拉克函数；82水 利 水 运 工 程 学 报2023 年 8 月uk(t)f(t)+jt)uk(t)uk(t)为分解序列模态调幅调频（AM-FM）信号；为径流序列；是的希尔伯特变换结果。对式（1）进行重构，同时引入拉格朗日系数，转化为非约束变分问题，扩展的拉格朗日表达式为：L(uk,wk,)=k?t(t)+jt)uk(t)ejwkt?2+?f(t)kuk(t)?22+(t),f(t)kuk(t)(2)ukwk式中：为拉格朗日乘子；为惩罚因子。将上述拉格朗日函数进行时-频域变换，并进行相应的极值求解，即可得到模态分量及中心频率的频域表达式：un+1k(w)=f(w)i,k ui(w)+(w)21+2(wwk)2(3)wn+1k=w0w|uk(w)|2dww0|uk(w)|2dw(4)un+1k(w)f(w)i,k ui(w)wn+1k uk(w)uk(t)ukwk式中：相当于当前剩余量的维纳滤波；为当前模态函数功率谱的重心，对进行傅里叶逆变换，实部则为。最后，在迭代子优化序列中运用乘子交替方向法(ADMM)，不断更新和得到最优解，进而求解分解层数 k 与各模态分量。1.21.2长短期记忆网络（LSTM）LSTM 神经网络是一种特殊的 RNN 网络，LSTM 通过其独特的门控设计，克服了传统 RNN 学习长期依赖关系的问题，适合于处理和预测时间序列中间隔和延迟非常长的重要事件18。在原始 RNN 中，只有一种内部状态 ht，对于短期序列的输入比较敏感，LSTM 的循环单元更为复杂，与 RNN 简单的 tanh 层相比，其拥有 4 个相互交替的网络层，如图 1 所示。遗忘门AA1ht1ht+1Xt+1Xt1ht1Ct1+tanhtanhXtOtitCtfthtCtht23输入门输出门图 1 LSTM 网络结构Fig.1 LSTM network structure diagram LSTM 每一个循环单元都有一个新增的单元状态 C，将前期输入信息长期保留，这是 LSTM 更适用于长期序列的原因。在 t 时刻，网络的输入包括 t-1 时刻的输出值 ht-1、t 时刻的输入值及 t-1 时刻的单元状态Ct-1；网络的输出包括 t 时刻的单元状态 Ct和输出值 ht。LSTM 的关键是对长期状态 C 的控制，LSTM 的神经网络模块使用了一种“门”（gate）结构，包括：输入门、遗忘门和输出门。输入门可将前序单元的信息输入到长期状态 C 中长期保留；遗忘门控制前一状态Ct-1中被遗忘的元素；输出门控制长期状态 C 在 t 时刻 LSTM 模块中的输出。第 4 期王秀杰，等：VMD-PSO-LSTM 模型的日径流多步预测83ft 1.2.1遗 忘 门遗忘门 通过读取上一时刻输出值 ht1及当前时刻输入值 xt，输出 0 到 1 之间的数值用于决定上一时刻单元状态 Ct1中信息保留比重。1 表示完全保留，0 表示完全舍弃。计算式如下：ft=(Wfht1xt+bf)(5)ftWfht1xt式中：为遗忘门的输出；为 sigmoid 函数；为遗忘门的权重；为上一时刻输出值；为当前时刻输入值；bf为遗忘门的偏置项。CtCt 1.2.2输 入 门输入门 it确定了当前输入的单元状态中用于更新当前时间步的最终单元状态 Ct的信息。其中，是由上一时刻输出 ht1和当前时刻输入 xt计算得出的当前输入下的单元状态信息，Ct取值范围为（1，1），it取值范围为（0，1），计算式如下：it=(Wiht1xt+bi)(6)Ct=tanh(W cht1xt+b c)(7)Ct=ftCt1+itCt(8)itWiW cCtb cCtCt1式中：为输入门的输出；为输入门的权重矩阵；bi为输入门的偏置项；为单元状态的权重矩阵；为单元状态的偏置项；Ct为当前时刻的单元状态；为上一时刻的单元状态。1.2.3输 出 门输出门 ot决定了 Ct中用于作为当前时刻输出 ht的信息。Ct经过 tanh 层处理，得到（1，1）之间的值。ot和 ht的计算式如下：ot=(Woht1xt+bo)(9)ht=ottanh(Ct)(10)otWoht式中：为输出门的输出；为输出门的权重矩阵；bo为输出门的偏置项；为当前时刻输出值。另外，本文采用 PSO 算法优化 LSTM 模型中隐含层节点数和学习率两个参数。PSO 是一种群体协作的随机搜索算法，通过不断迭代计算，协调粒子自身和群体的最优值对粒子当前运动方向和运动速度的影响，全局极值不断更新来找出问题的最优解。具体原理见文献 19，这里不再赘述。1.31.3VMD-PSO-LSTM 径流预测模型VMD-PSO-LSTM 径流预测模型综合了 VMD 分解方法自适应、完全非递归的特性及 PSO-LSTM 神经网络参数自动寻优、长期数据记忆的优势。模型的建立包含以下 4 个步骤：步骤 1：VMD 序列分解。使用 VMD 算法将原始径流序列分解成 k 个离散的固有模态分量（SIMF），使非平稳原序列分解得到多个较稳定的子序列。VMD 分解效率和预测精度主要受分解层数 k 和惩罚因子等参数的影响，VMD 分解虽具有较好的分解效果，但存在确定分解层数 k 的不确定性问题，因此本文采用排列熵优化算法（permutation entroy optimization，PEO）实现分解层数 k 的自适应确定。该算法具体步骤参考相关文献20。步骤 2：PSO 参数优化。利用 PSO 算法对 LSTM 网络的参数进行优化，使得网络隐含层节点数（NHN）和学习率（）达到最优。步骤 3：LSTM 网络训练。在参数优化的基础上建立预测模型，使用步骤 2 得到的最优参数建立LSTM 网络，对每个子序列进行 LSTM 网络训练。在训练之前需将分解后的各 SIMF进行归一化处理（由于径流序列具有高度非线性和非平稳性，因此必须对输入数据进行归一化处理，以避免训练过程中数据出现大的波动）。其中，直接多步预测在模型中以滑动窗口机制实现，假定一段径流时间序列为 X=x1,x2,xn，模型输入步长为 L，即滑动窗口的当前窗口长度为 L；滑动窗口的移动步长为 1；当预测预见期为 a（本文预见期设84水 利 水 运 工 程 学 报2023 年 8 月为 1、2、3 d）时，当前 t 时刻向后预测 a 天的直接多步预测模型输入输出结构如图 2 所示。步骤 4：序列预测。利用 PACF 偏自相关分析法确定的输入变量，在训练好的 LSTM 网络的基础上对每个子序列进行预测，各 SIMF预测值叠加即为最终预测结果。2 实例分析本文采用黄河下游花园口和利津站 20032012 年的逐日径流资料作为研究数据（数据均来源于黄河水利委员会），以 20032009 年的日径流序列为训练数据，20102012 年的日径流序列为验证数据。花园口和利津站的日径流序列过程见图 3。VMD 分解方法具有自适应性和强大的降噪性能，针对径流序列的非平稳性，将频率不同且成分复杂的径流时间序列分解成具有固有模态函数的平稳性序列分量 SIMF，采用 VMD 法分别对花园口和利津站的日径流序列进行分解，经过排列熵算法，自适应得到最优 k 值，从而获取 k 个序列分量。花园口和利津站的日径流序列经分解得到的子序列个数分别为 6 和 8 个。参考相关文献15,21，本文模型惩罚因子取 2 000。以花园口站为例，研究数据的分解结果如图 4 所示，可以看出，通过 VMD 分解，将原始日径流序列分为了由高频到低频的周期性模态分量与表示趋势的残差，揭示出日径流数据内的隐藏信息，即周期性震荡变化和趋势。这样不仅使模型能更好地理解周期性的信息，同时也有利于提高模型的训练和预测精度。对于上述每个 SIMF分别建立 PSO-LSTM 预测模型，各子序列预测结果的叠加即原径流序列的预测结果。其中模型的输入步长利用偏自相关函数（PACF）来确定，可突出时间滞后对当前时段径流的影响。假设输出变量为 xi，当滞后 L 的偏自相关系数 fPACF超出 95%置信区间时，前 L 个变量即为输入变量22。以花园口站为例，原始日径流序列经过 VMD 分解得到 6 个子序列，由 PACF 计算得到各 SIMF的输入步长分别为 1、6、4、4、8、11，如图 5。具体的输入步长及输入变量如表 1 所示。LSTM 网络隐含层节点数（NHN）和学习率（）利用 PSO 算法确定，隐含层数为 2 层。以花园口站为例，每个 SIMF在 13 d 预见期对应的最优NHN 和 值见表 2。模型的输入层节点数等于输入变量数，输出层节点数为 1。将当前时段 t 之前的数据作为模型的输入，分别将 t 时刻、t+1 时刻、t+2 时刻的预测数据作为输出构建预测模型，实现预见期 13 d 的多步预测。输入输出t=Kt=K+1t=K+2xK LxK 2xK 1xK 1xKxKxK+1xK L+1xK L+1xK L+2xK+a 1xK+axK+a+1xK L+2xK L+3 图 2 直接多步预测示意Fig.2 Schematic diagram of direct multi-step prediction 0500 1 000 1 500时间/d利津花园口2 000 2 500 3 000 3 500 4 0001 0002 0003 000日径流/(m3s1)4 0005 000训练数据验证数据 图 3 花园口和利津站的日径流序列Fig.3 The daily runoff sequence at Huayuankouand Lijin Stations 第 4 期王秀杰，等：VMD-PSO-LSTM 模型的日径流多步预测85 表 1 花园口和利津站各SIMF的输入步长及输入变量Tab.1 Input steps and input variables of each SIMF at Huayuankou and Lijin Stations站名SIMF输入步长/d输入变量花园口SIMF11x1(t-1)SIMF26x2(t-1)(t-2)(t-3)(t-4)(t-5)(t-6)SIMF34x3(t-1)(t-2)(t-3)(t-4)SIMF44x4(t-1)(t-2)(t-3)(t-4)SIMF58x5(t-1)(t-2)(t-3)(t-4)(t-5)(t-6)(t-7)(t-8)SIMF611x6(t-1)(t-2)(t-3)(t-4)(t-5)(t-6)(t-7)(t-8)(t-9)(t-10)(t-11)利津SIMF11x1(t-1)SIMF21x2(t-1)SIMF33x3(t-1)(t-2)(t-3)SIMF44x4(t-1)(t-2)(t-3)(t-4)SIMF53x5(t-1)(t-2)(t-3)SIMF64x6(t-1)(t-2)(t-3)(t-4)SIMF78x7(t-1)(t-2)(t-3)(t-4)(t-5)(t-6)(t-7)(t-8)SIMF88x8(t-1)(t-2)(t-3)(t-4)(t-5)(t-6)(t-7)(t-8)表 2 花园口站各SIMF对应隐含层节点数（NHN）和学习率（）Tab.2 The number of hidden layer nodes(NHN)and learning rate()corresponding to SIMF at Huayuankou StationSIMF1 d2 d3 dNHNNHNNHNSIMF1157,1910.003 9199,1530.006 6190,900.004 9SIMF2164,1810.003 8194,930.005 1157,1910.003 9SIMF3137,1820.004 6190,880.005 0157,1910.003 9SIMF4140,1870.004 6157,1910.003 9190,900.004 9SIMF5189,1120.005 1191,530.005 3192,1110.005 8SIMF6191,450.008 4190,900.004 9195,1490.005 0注：因采用的LSTM模型隐含层个数为2，逗号前后分别为其对应的最优隐含层节点数。SIMF41 00001 0000500 1 0001 5002 0002 5003 0003 5004 000SIMF21 50001 5000500 1 0001 5002 0002 5003 0003 5004 000SIMF31 50001 5000500 1 0001 5002 0002 5003 0003 5004 000SIMF580008000500 1 0001 5002 0002 5003 0003 5004 000800SIMF608000500 1 0001 5002 0002 5003 0003 5004 000SIMF1时间/d0500 1 0001 5002 0002 5003 0003 5004 0001 5003 000 图 4 花园口站 VMD 分解各分量结果Fig.4 VMD decomposition results of each component atHuayuankou Station 110fPACF1110fPACF2110fPACF3110fPACF4110fPACF51051015时滞/d20253005101520253005101520253005101520253005101520253005101520253010fPACF6 图 5 花园口站各子序列的 fPACFFig.5 PACF method input step recognition of sub-sequence atHuayuankou Station 86水 利 水 运 工 程 学 报2023 年 8 月 3 结果分析基于 VMD 的分解结果，根据输入步长，将 t 时刻之前的数据作为输入变量，构建每个 SIMF的 PL 模型，将 t 时刻、t+1 时刻、t+2 时刻的数据作为输出结果，分别将各子序列的预测结果进行加和，实现 13 d 的VPL 多步预测。为了进一步验证 VPL 模型在径流多步预测中的优越性，建立 PSO-LSTM（PL）、CEEMD-PSO-LSTM(CPL)多步预测模型进行对比。建模所用的训练数据和验证数据同 VPL 模型。其中，日径流经CEEMD 分解，花园口得到13 个子序列，利津12 个子序列，比VMD 分解的子序列个数分别多了7 和4 个。VPL模型、CPL 模型和 PL 模型的预测结果评价指标见表 3。表 3 不同模型的预测结果对比Tab.3 Comparison of prediction results of different models测站预见期/dVPL模型CPL模型PL模型ENSERMS/(m3s1)RENSERMS/(m3s1)RENSERMS/(m3s1)R花园口10.954159.8560.9820.908226.3040.9560.907227.1090.95020.934192.4890.9700.879259.3210.9530.803331.6550.89630.902233.1600.9560.854285.6730.9540.671427.9490.833利津10.959150.7960.9810.946172.5330.9760.940182.5580.97320.914218.1970.9600.906228.7880.9590.849288.8760.92530.918212.9000.9700.888249.6730.9520.739380.3140.869 对比 VPL 模型、CPL 模型和 PL 模型的预测结果指标可知，总体上 3 种模型的预测精度：VPL 模型CPL 模型PL 模型。在验证期，VPL 模型的预测精度最高，在 13 d 预见期下，其 ENS、ERMS、R 基本上均优于其他两种模型。另外，经过分解后的 VPL、CPL 组合模型预测精度明显优于未经分解的 PL 模型，在13 d 预见期下，VPL 模型对 PL 模型各指标的改进效果比 CPL 模型更为显著。随着预见期增长，分解后模型各个指标的改进效果更显著。在 13 d 的预见期下，VPL 模型比 CPL 模型的 ENS提高了 5.1%、9.3%、2.6%，ERMS降低了 29.4%、32.6%、10.1%，R 值提高了 2.7%、1.8%、0.3%，可能的原因是 VMD 通过自适应 k 值，减少了子序列的个数，从而减少了子序列预测结果叠加过程中产生的累计误差。随着预见期的增长，VPL 和 CPL 模型的预测精度都出现了降低的趋势，以花园口站为例，VPL 模型在 13 d 的预见期下，ENS能够保持 0.9 以上，而 CPL 已经下降到 0.854。这表明 VPL 模型保证了预见期增长时预测精度的稳定性。以花园口站为例，VPL 模型、CPL 模型、PL 模型的预测径流过程和散点图如图 6 所示。从径流过程对比图可以看出，VPL 模型的预测径流过程最接近实测径流过程，且 VPL 模型的峰值预测精度要明显优于CPL、PL 模型，虽然随着预见期的增长，3 种模型的预测峰值在时间和大小上与实测值相比都有了不同程度的偏移，但是 VPL 模型的预测峰值的偏移程度相对微弱。由散点图可知 VPL 模型在 13 d 预见期下的R2值比 CPL、PL 模型更高，预测结果更加接近实测径流，尤其是中高流量，VPL 预测结果散点图相对 CPL、PL 模型更加集中于 45线，说明 VPL 模型对中高流量的预测精度较好。综合以上分析，与 CPL、PL 模型预测结果相比，在预见期为 13 d 的情况下，VPL 模型均表现出更稳定、更精确的优势，实现了高精度的日径流多步预测。第 4 期王秀杰，等：VMD-PSO-LSTM 模型的日径流多步预测87 4 结语针对日径流时间序列的强非线性和非平稳性，本文提出了 VMD-PSO-LSTM 混合多步预测模型。为了验证该混合模型在日径流多步预测中的适用性，对黄河下游花园口和利津站进行了预见期 13 d 的多步预测，与 PSO-LSTM、CEEMD-PSO-LSTM 模型预测结果对比，得出以下结论：（1）与未分解模型相比，“分解-预测-重构”模式的组合预测模型的预测精度明显提高，且随着预见期增长，分解后模型的优势更加显著。可见信号分解方法有助于预测模型更好地学习径流序列中的周期、趋势等本质特征，提高日径流的预测精度。（2）与其他分解方法相比，VMD 方法更能有效提取日径流序列的特征信息，减弱原始日径流序列的非 02010-01-012010-03-012010-05-012010-07-012010-09-012010-11-01日期(a1)径流过程(1 d 预见期)2010-12-311 0002 0003 000日径流/(m3s1)4 0005 000实测值VPLCPLPL日期(a2)径流过程(2 d 预见期)02010-01-012010-03-012010-05-012010-07-012010-09-012010-11-012010-12-311 0002 0003 000日径流/(m3s1)4 0005 000实测值VPLCPLPL5 000日期(a3)径流过程(3 d 预见期)02010-01-012010-03-012010-05-012010-07-012010-09-012010-11-012010-12-311 0002 0003 000日径流/(m3s1)4 000实测值VPLCPLPL01 0006 0002 0004 000R=0.956R=0.982R=0.9502 0003 0004 000预测值/(m3s1)5 0006 000实测值/(m3s1)(b1)散点图(1 d 预见期)VPLCPLPL01 0006 0002 0004 000R=0.953R=0.970R=0.8962 0003 0004 000预测值/(m3s1)5 0006 000实测值/(m3s1)(b2)散点图(2 d 预见期)VPLCPLPL01 0006 000实测值/(m3s1)(b3)散点图(3 d 预见期)2 0004 000R=0.950R=0.956R=0.8332 0003 0004 000预测值/(m3s1)5 0006 000VPLCPLPL 图 6 花园口站 2010 年不同预见期各模型预测结果与实测值的比较Fig.6 Comparison of the predicted results and measured values of Huayuankou Station in different forecast periods in 2010 88水 利 水 运 工 程 学 报2023 年 8 月平稳性，进而得到更优的预测结果。因此该分解方法的频率分离性能更佳，对噪声干扰具有较好的鲁棒性，可以有效减少伪分量。（3）VMD-PSO-LSTM 混合多步预测模型在预测精度、预测稳定性方面展示出独特优势，尤其对中、高值流量预报表现出较强的跟踪预报能力，为实现复杂非平稳日径流时间序列的高精度多步预测提供了一种新途径。参考文献：梁浩,黄生志,孟二浩,等.基于多种混合模型的径流预测研究J.水利学报，2020，51(1)：112-125.（LIANG Hao,HUANGShengzhi,MENG Erhao,et al.Runoff prediction based on multiple hybrid modelsJ.Journal of Hydraulic Engineering,2020,51(1):112-125.（in Chinese）1 崔东文.改进Elman神经网络在径流预测中的应用J.水利水运工程学报，2013(2)：71-77.（CUI Dongwen.An improvedElman neural network and its application to runoff forecastJ.Hydro-Science and Engineering,2013(2):71-77.（in Chinese）2 HOCHREITER S,SCHMIDHUBER J.Long short-term memoryJ.Neural Computation,1997,9(8):1735-1780.3 HU C H,WU Q,LI H,et al.Deep learning with a long short-term memory networks approach for rainfall-runoff simulationJ.Water,2018,10(11):1543.4 YUAN X H,CHEN C,LEI X H,et al.Monthly runoff forecasting based on LSTM-ALO modelJ.Stochastic EnvironmentalResearch and Risk Assessment,2018,32(8):2199-2212.5 阚光远,洪阳,梁珂.基于耦合机器学习模型的洪水预报研究J.中国农村水利水电，2018(10)：165-169,176.（KANGuangyuan,HONG Yang,LIANG Ke.Research on the flood forecasting based on coupled machine learning modelJ.ChinaRural Water and Hydropower,2018(10):165-169,176.（in Chinese）6 白继中,师彪,冯民权,等.基于自适应调整蚁群-RBF神经网络模型的中长期径流预测J.自然资源学报，2011，26(6)：1065-1074.（BAI Jizhong,SHI Biao,FENG Minquan,et al.Adaptive regulation ant colony system algorithm-radial basisfunction neural network model and application in mid-long term runoff predictionJ.Journal of Natural Resources,2011,26(6):1065-1074.（in Chinese）7 王立辉,杨辉斌,王银堂,等.基于GWO-LSTM的丹江口水库入库径流预测J.水利水运工程学报，2021(6)：51-59.（WANG Lihui,YANG Huibin,WANG Yintang,et al.Prediction of inflow to the Danjiangkou reservoir based on GWO-LSTMJ.Hydro-Science and Engineering,2021(6):51-59.（in Chinese）8 郑恺原,向小华,吴晓玲.基于PSO-WRF的渭河中下游年径流预测模型J.水利水电技术，2020，51(8)：39-44.（ZHENGKaiyuan,XIANG Xiaohua,WU Xiaoling.Prediction model of annual runoff in the middle and lower reaches of Wei Riverbased on PSO-WRFJ.Water Resources and Hydropower Engineering,2020,51(8):39-44.（in Chinese）9 周婷,金菊良,李荣波,等.基于小波支持向量机的径流预测性能优化分析J.水力发电学报，2017，36(10)：45-55.（ZHOUTing,JIN Juliang,LI Rongbo,et al.Performance optimization analysis for inflow prediction using wavelet Support VectorMachineJ.Journal of Hydroelectric Engineering,2017,36(10):45-55.（in Chinese）10 孙娜,周建中.基于正则极限学习机的非平稳径流组合预测J.水力发电学报，2018，37(8)：20-28.（SUN Na,ZHOUJianzhong.Hybrid forecasting model for non-stationary runoff based on regularized extreme learning machineJ.Journal ofHydroelectric Engineering,2018,37(8):20-28.（in Chinese）11 王文川,高畅,徐雷.基于TVF-EMD与LSTM神经网络耦合的月径流预测研究J.中国农村水利水电，2022(2)：76-81,89.（WANG Wenchuan,GAO Chang,XU Lei.Monthly runoff forecast based on TVF-EMD and LSTM neural networkcouplingJ.China Rural Water and Hydropower,2022(2):76-81,89.（in Chinese）12 黄景光,吴巍,程璐瑶,等.基于小波支持向量机特征分类的日径流组合预测：以宜昌三峡水库为例J.中国农村水利水电，2018(6)：33-39.（HUANG Jingguang,WU Wei,CHENG Luyao,et al.Daily runoff combination prediction based onwavelet support ve