基于集合经验模态分解的自适应调度期末水位预测方法技术

本发明专利技术公开了一种基于自适应集合经验模态分解的调度期末水位预测方法，涉及水库运行调度及水利信息化技术领域。该方法充分考虑了调度期末水位序列的非稳态性，通过使用集合经验模态分解的方法将调度期末水位转化为多组稳态序列，实现了水文序列稳态化，为常规的预报方法提供了最基础的数据条件。另外，利用本发明专利技术提供的自适应预报模型进行预报时，是一种滚动预报作业，使得模型实现了实时校正，保证了模型的适应性，为精准预报保证了模型基础。所以，采用本发明专利技术提供的方法，建立的预测模型用于调度期末水位预测时，具有较高的精准性。

【技术实现步骤摘要】
基于集合经验模态分解的自适应调度期末水位预测方法
本专利技术涉及水库运行调度及水利信息化
，尤其涉及一种基于自适应集合经验模态分解的调度期末水位预测方法。
技术介绍
调度期末水位是水库调度的重要组成部分，调度期末水位预测对水库调度的余留效益、风险评估等都有重要的意义。水库调度期末水位具有趋势性、周期性和随机性，呈非稳态序列，同时该水位受气候变量、来水量多少、下垫面情况、调度决策者的调度方式等多因素影响，隐藏多种确定性的和不确定性的多层信息，故要准确的预报该水位，需要对多种影响因素进行深入挖掘。当前，对调度期末水位预测方法较少，主要是通过预报径流量，调度规则推求调度期末的水位，或者通过不同的智能优化算法如人工神经网络、支持向量机、线性回归等进行自回归模型预测水位。然而这些智能优化算法是一个黑箱模型，对水位的非稳态识别比较困难，另外，由于这些智能优化算法忽略了水位过程中的水文物理机理，所以，对水位进行自相关分析，会导致预报过拟合，出现失真现象。
技术实现思路
本专利技术的目的在于提供一种基于自适应集合经验模态分解的调度期末水位预测方法，从而解决现有技术中存在的前述问题。为了实现上述目的，本专利技术采用的技术方案如下：一种基于集合经验模态分解的自适应调度期末水位预测方法，包括如下步骤：S1，数据预处理：提取一列长度为L连续多调度周期的调度期末水位序列为自变量LS，利用调度初期的初水位L0、t时刻水位Lt和t时刻余留期间的来水量W组成因变量集合X＝[L0,Lt,W]；S2，自变量稳态化分解：将自变量LS序列，采用集合经验模态分解的方法分解，生成多组稳态本征模态函数IMF(n)和一组残余序列Res，共同组成自变量组合Y＝[IMF(n),Res]；S3，确定预报因子：利用因变量集合X＝[L0,Lt,W]与自变量组合Y的相关关系确定预报决策因子，按照如下步骤进行实施：S301，构建备选预报因子集合：将因变量集合X＝[L0,Lt,W]中的三个影响因子序列L0、Lt和W，分别以单独，两两、共同的方式组合，共构成七组备选的预报因子集合，分别为F1＝[L0]、F2＝[Lt]、F3＝[W]、F4＝[L0,Lt]、F5＝[L0,W]、F6＝[Lt,W]、F7＝[L0,Lt,W]；S302，分别对七组备选的预报因子集合与自变量组合Y＝[IMF(n),Res]中每组数据进行相关性分析，得到相关因子；S303，根据S302得到的相关因子，选取与Y＝[IMF(n),Res]中每组序列相关性最大的备选预报因子作为每组自变量的最终预报因子F，完成预报因子的识别工作；S4，根据每组自变量的最终预报因子构建自适应调度期末水位预报模型；S5，根据所述预报模型对调度期末水位进行预报，得到自适应调度期末水位预报结果。优选地，S303中，通过T检验，选取与Y＝[IMF(n),Res]中每组序列相关性最大的备选预报因子作为每组自变量的最终预报因子F。优选地，S4包括如下步骤：S401，建立训练样本集：将预报根据样本数据量的大小确定模型训练期长度M和验证期的长度(L-M)，由于模型为滚动预报，故预见期的步长为1个时间步长；S402，建立三层BP人工神经Y＝[IMF(n),Res]预测模型：将调度期末期水位序列自变量LS，基于经验模态方法分解出的IMF(n),Res稳态序列，分别与各自的预报因子F建立n+1个预报模型，根据这些预报模型分别对调度期末期水位序列自变量LS进行预报，得到每组的预报值，分别为IMFf(n)、Resf；S403，根据预报出来的每层预报值，按照如下公式计算M+1时刻预测的调度周期末的水库水位L(M+1)f：其中，IMFf(j)为第j层本征模态函数分量的预测值，Resf为预测的残余量；S404，将预报出来L(M+1)f添加到训练样本中，跳至S401，进行下一轮的水位预报，直到完成L时刻的预测，得到L时刻预测的调度周期末的水库水位LLf；S405，对比预报出的L(M+1)f...LLf与(M+1)到L时刻的调度期末水库水位LS，对其预报效果进行评价；S406，若评价效果符合设定的阈值，则自适应模型建立完成，该模型用于对未来时段调度期末水库水位进行预报，否则，重新调整BP神经网络的模型参数，重新建模直至模型完成建立。优选地，S405中，所述对其预报效果进行评价，评价指标包括纳什效率系数、相对误差和合格率；所述纳什效率系数Nash按照如下公式进行计算：其中，Lf为预测的调度期末水位，为实际调度期末水位的均值，Nash越接近1，预报越精准；所述相对误差MARE按照如下公式进行计算：其中，N为的序列长度，MARE越接近0，说明实测与预报值越接近，预报效果越精准，常认为MARE＜20％时，预报效果较好；所述合格率QR按照如下公式进行计算：其中，n为合格预报次数；m为预报总次数，当QR>80％时，认为预报效果较好。本专利技术的有益效果是：本专利技术实施例提供的基于自适应集合经验模态分解的调度期末水位预测方法，充分考虑了调度期末水位序列的非稳态性，通过使用集合经验模态分解的方法将调度期末水位转化为多组稳态序列，实现了水文序列稳态化，为常规的预报方法提供了最基础的数据条件。另外，利用本专利技术提供的自适应预报模型进行预报时，是一种滚动预报作业，使得模型实现了实时校正，保证了模型的适应性，为精准预报保证了模型基础。所以，采用本专利技术提供的方法，建立的预测模型用于调度期末水位预测时，具有较高的精准性。附图说明图1为本专利技术提供的基于自适应集合经验模态分解的调度期末水位预测方法流程示意图；图2为本专利技术采用的集合经验模态分解方法流程示意图；图3为基于集合经验模态分解方法对水位序列进行分解后的序列示意图；图4为自适应调度期末水位预测效果示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不用于限定本专利技术。为了解决水位为非稳态序列限制常规预报方法的预报，来水不确定性无法量化的问题。本专利技术提供了一种自适应集合经验模态分解的调度期末水位预测方法，有效地将调度期末水位稳态化，挖掘了响应其水位的影响因素，显性化和量化了来水不确定性对水库调度期末水位的影响，同时自适应的预报模型，使预报模型保持不断地更新，适应调度期末的水位变化过程，保持模型的稳定性。如图1所示，本专利技术实施例提供了一种基于集合经验模态分解的自适应调度期末水位预测方法，包括如下步骤：S1，数据预处理：提取一列长度为L连续多调度周期的调度期末水位序列为自变量LS，利用调度初期的初水位L0、t时刻水位Lt和t时刻余留期间的来水量W组成因变量集合X＝[L0,Lt,W]；S2，自变量稳态化分解：将自变量LS序列，采用集合经验模态分解的方法分解，生成多组稳态本征模态函数IMF(n)和一组残余序列Res，共同组成自变量组合Y＝[IMF(n),Res]；S3，确定预报因子：利用因变量集合X＝[L0,Lt,W]与自变量组合Y的相关关系确定预报决策因子，按照如下步骤进行实施：S301，构建备选预报因子集合：将因变量集合X＝[L0,Lt,W]中的三个影响因子序列L0、Lt和W，分别以单独，两两、共同的方式组合，共构成七组备选的本文档来自技高网...

【技术保护点】
1.一种基于集合经验模态分解的自适应调度期末水位预测方法，其特征在于，包括如下步骤：S1，数据预处理：提取一列长度为L连续多调度周期的调度期末水位序列为自变量LS，利用调度初期的初水位L0、t时刻水位Lt和t时刻余留期间的来水量W组成因变量集合X＝[L0,Lt,W]；S2，自变量稳态化分解：将自变量LS序列，采用集合经验模态分解的方法分解，生成多组稳态本征模态函数IMF(n)和一组残余序列Res，共同组成自变量组合Y＝[IMF(n),Res]；S3，确定预报因子：利用因变量集合X＝[L0,Lt,W]与自变量组合Y的相关关系确定预报决策因子，按照如下步骤进行实施：S301，构建备选预报因子集合：将因变量集合X＝[L0,Lt,W]中的三个影响因子序列L0、Lt和W，分别以单独，两两、共同的方式组合，共构成七组备选的预报因子集合，分别为F1＝[L0]、F2＝[Lt]、F3＝[W]、F4＝[L0,Lt]、F5＝[L0,W]、F6＝[Lt,W]、F7＝[L0,Lt,W]；S302，分别对七组备选的预报因子集合与自变量组合Y＝[IMF(n),Res]中每组数据进行相关性分析，得到相关因子；S303，根据S302得到的相关因子，选取与Y＝[IMF(n),Res]中每组序列相关性最大的备选预报因子作为每组自变量的最终预报因子F，完成预报因子的识别工作；S4，根据每组自变量的最终预报因子构建自适应调度期末水位预报模型；S5，根据所述预报模型对调度期末水位进行预报，得到自适应调度期末水位预报结果。...

【技术特征摘要】
1.一种基于集合经验模态分解的自适应调度期末水位预测方法，其特征在于，包括如下步骤：S1，数据预处理：提取一列长度为L连续多调度周期的调度期末水位序列为自变量LS，利用调度初期的初水位L0、t时刻水位Lt和t时刻余留期间的来水量W组成因变量集合X＝[L0,Lt,W]；S2，自变量稳态化分解：将自变量LS序列，采用集合经验模态分解的方法分解，生成多组稳态本征模态函数IMF(n)和一组残余序列Res，共同组成自变量组合Y＝[IMF(n),Res]；S3，确定预报因子：利用因变量集合X＝[L0,Lt,W]与自变量组合Y的相关关系确定预报决策因子，按照如下步骤进行实施：S301，构建备选预报因子集合：将因变量集合X＝[L0,Lt,W]中的三个影响因子序列L0、Lt和W，分别以单独，两两、共同的方式组合，共构成七组备选的预报因子集合，分别为F1＝[L0]、F2＝[Lt]、F3＝[W]、F4＝[L0,Lt]、F5＝[L0,W]、F6＝[Lt,W]、F7＝[L0,Lt,W]；S302，分别对七组备选的预报因子集合与自变量组合Y＝[IMF(n),Res]中每组数据进行相关性分析，得到相关因子；S303，根据S302得到的相关因子，选取与Y＝[IMF(n),Res]中每组序列相关性最大的备选预报因子作为每组自变量的最终预报因子F，完成预报因子的识别工作；S4，根据每组自变量的最终预报因子构建自适应调度期末水位预报模型；S5，根据所述预报模型对调度期末水位进行预报，得到自适应调度期末水位预报结果。2.根据权利要求1所述的基于自适应集合经验模态分解的调度期末水位预测方法，其特征在于，S303中，通过T检验，选取与Y＝[IMF(n),Res]中每组序列相关性最大的备选预报因子作为每组自变量的最终预报因子F。3.根据权利要求1所述的基于自适应集合经验模态分解的调度期末水位预测方法，其特征在于，S4包括如下步骤：S401，建立训练样本...

【专利技术属性】
技术研发人员：王旭，雷晓辉，王佳，王超，权锦，田雨，王浩，蔡思雨，张城，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京,11