一种考虑动库容影响的水库超短期水位预测方法、系统及存储介质

技术领域

本发明属于水库调度技术领域，具体涉及一种考虑动库容影响的水库超短期水位预测方法、系统及存储介质。

背景技术

在水库实际运行过程中，考虑水库来水信息和发电计划的水库超短期水位预测对于水电站发电计划制定、水库安全运行具有重要意义，基于水位预测结果可在符合调度规程的框架下对水库出力计划进行适当的修正，生产满足调度需求且发电效益最大的出力计划。

对于一般水库而言，常通过基于静库容原理的水量平衡法对水库水位进行预测，在实际应用中可以取得良好的水位预测结果。在该方法中，基于“入库水量-出库水量＝水库蓄量变化”的水量平衡关系计算水库蓄水量的变化，从而推算坝前水位变化，其中入库水量采用的是静库容入库水量，且计算中只考虑水库静库容的变化，并假定水库水面为水平面，其涨落为水平涨落。但是，对于河道型水库而言，水库动库容(最远回水断面到坝前的总蓄水量)的存在使水库水位预测面临着两方面问题：

(1)河道型水库的水面具有一定的纵向比降，水库的实际水面与坝前水位的水平线之间会形成楔形水库容积，水库蓄水量变化不仅会反映到坝前水位的变化中，也会影响库区水面线，且库区水面线变化的水动力学特征复杂，因此根据水库蓄水量变化难以直接推算水库坝前水位变化；

(2)水库的动库容入库水量仅能反映水库回水断面处的水量，其需通过库区水面线变化反映到坝前水位的变化中，因此动库容入库水量对于坝前水位的影响具有一定的滞后性。

总体而言，河道型水库动库容的存在使水库蓄量变化难以准确计算，动库容入库水量对水库坝前水位影响的滞后性也造成了水位预测的困难。

发明内容

针对现有技术存在的缺陷，本发明提供了一种考虑动库容影响的水库超短期水位预测方法，此方法能够在水库超短期水位预测计算中考虑动库容的影响，实现快速、高精度的水库超短期水位预测，从而为水电站发电计划制定及水库安全运行提供更好的支撑。

为了实现上述的技术特征，本发明的目的是这样实现的：一种考虑动库容影响的水库超短期水位预测方法，包括以下步骤：

步骤S1：数据分析及处理：

收集相应时间尺度的水库运行数据，根据设定参数和收集的水库运行数据构建数据集(X

步骤S2：建立预报模型：

基于步骤S1构建的数据集(X

步骤S3：产生水位预测结果：

基于步骤S2建立的预报模型，产生水位预测结果。

所述步骤S1具体包括：

步骤S101：根据实际业务需求设定时间步长step、预见期flt；

步骤S102：确定库区水面线水位对坝前水位影响的最大滞后期lags和入库流量对坝前水位影响的滞后时间QLags；

步骤S103：根据设定的时间步长step收集相应时间尺度的水库运行数据；

所述水库运行数据包括电站出力数据N、库区水面线水位数据WL、坝前水位数据WLUp、下游水位数据WLDown、动库容入库流量数据Q以及弃水流量数据QOut，并截取同时段的数据，截取后数据长度记为num

步骤S104：根据设定参数和收集数据构建数据集，数据集中包含num＝num

Y

Y

X

X

X

X

所述步骤S2具体包括：

步骤S201：数据预处理：

采用归一化或者标准化方法将构建的数据集(X

步骤S202：数据划分：

将标准数据集按照一定的比例随机划分为训练集(TrainSX

步骤S203：模型结构及参数设计：

根据所构建的数据集的特征，利用深度学习中的长短期记忆层、卷积层、全连接层设计能够提取其中库区水面线水位数据、流量数据、出力数据变化特征的模型结构，设定模型的优化目标，并设定初步的模型超参数集；

步骤S204：模型训练：

在给定的模型结构和超参数集条件下，基于训练集(TrainSX

步骤S205：模型评估：

基于训练好的模型和优化目标，评估在验证集(TestSX

步骤S206：确定预报模型：

确定模型评估结果是否满足预测精度需求，若否，则调整模型结构、优化目标和超参数，重复执行步骤S204、步骤S205，直至获取符合预测精度需求的预报模型。

所述步骤S3具体包括：

步骤S301：在进行预报的时刻T，按照步骤S103、步骤S104收集相应数据，构建该时刻T对应的样本点的输入X

步骤S302：采用步骤S201中的数据预处理方法，将输入X

步骤S303：采用步骤S206确定的预报模型，预测得到预测结果SYHat

步骤S304：根据步骤S201中采用的数据预处理方法，将预测结果SYHat

一种考虑动库容影响的水库超短期水位预测系统，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现所述的考虑动库容影响的水库超短期水位预测方法。

一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的考虑动库容影响的水库超短期水位预测方法。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，本发明提供的一种考虑动库容影响的水库超短期水位预测方法，主要具有以下有益效果：

1、针对河道型水库动库容影响下的水位预测问题，本发明在预报输入中充分考虑动库容入库流量、库区水面线对坝前水位影响的滞后性以及前期库区水面线水位特征反映的蓄水量，从而更好地反映坝前水位变化的影响因素，提升预测的精度和可靠性。

2、针对库区水面线变化情况下水库蓄量难以求解的问题，本发明通过构建深度学习模块以提取历史运行过程中不同条件下的库区水面线变化特征，从而实现在复杂水动力学条件下的坝前水位预测计算。

附图说明

下面结合附图和实施例对本发明作进一步说明。

图1为本发明提供的考虑动库容影响的水库超短期水位预测方法的流程示意图。

图2为实施例中的优化后的模型结构图。

图3为实施例中某一次预报的预报水位和实测水位对比图。

图4为实施例中全年预报的整体误差随预见期的变化图。

图5为实施例中全年预报的整体误差分布图。

具体实施方式

下面结合附图对本发明的实施方式做进一步的说明。

实施例1：

如图1所示，本发明一方面提供一种考虑动库容影响的水库超短期水位预测方法，其通过对库区水面线时空信息的提取以反映水库动库容的影响，通过对入库径流滞后以反映动库容入库径流对坝前水位变化影响的滞后性。所述方法包括如下步骤：

步骤S1：数据分析及处理：收集相应时间尺度的水库运行数据，根据设定参数和收集的水库运行数据构建数据集(X

步骤S101：根据实际业务需求设定计算的时间步长step、预见期flt(指预测后续flt步)。

在本实施例中，设置时间步长为1小时，预见期为48小时，即step＝1，flt＝48。

步骤S102：通过机理分析或数据相关性分析，确定库区水面线水位对坝前水位影响的最大滞后期lags和入库流量对坝前水位影响的滞后时间QLags。

在本实施例中，根据动库容水动力学特性，确定库区水面线水位对坝前水位影响的最大滞后期lags为24小时，入库流量对坝前水位影响的滞后时间QLags为8小时，即lags＝24，QLags＝8。

步骤S103：根据设定的时间步长step收集相应时间尺度的水库运行数据，其中，所述水库运行数据包括电站出力数据N、库区水面线水位数据WL、坝前水位数据WLUp、下游水位数据WLDown、动库容入库流量数据Q以及弃水流量数据QOut，并截取同时段的数据，截取后数据长度记为num

在本实施例中，收集了某水电站(记为A电站)2015-2019年运行期间的电站出力数据N、库区水面线的9个点的水位数据WL、坝前水位数据WLUp、下游水位数据WLDown、动库容入库流量(可采用预报入库流量)数据Q以及弃水流量数据QOut，由于不同数据的时间精度不同，因此，通过数据插补或降尺度方法将其处理为1小时尺度数据，处理后的数据长度为num

步骤S104：根据设定参数和收集数据构建数据集。数据集中包含num＝num

Y

Y

X

X

X

X

在本实施例中，在2015-2019年期间选择一时间点，以2019年3月15日0时为例，其待预报数据为2019年3月15日1时至2019年3月17日0时期间的坝前水位数据(48×1矩阵)，其输入包含前期水位数据和同期出力、流量数据，其中前期水位数据为2019年3月14日1时至2019年3月15日0时期间的水面线9个点的水位以及坝前水位、下游水位，为24×11的矩阵，其同期的出力、流量数据为2019年3月15日1时至2019年3月17日0时期间的出力数据、弃水流量数据以及2019年3月14日17时至2019年3月16日16时期间的入库流量数据，为48×3的矩阵。通过这种构建方式形成的数据集中共包含43825-24-48+1＝43754个样本点，将所有样本结合起来构成的数据集包含输出Y和输入X＝(X

步骤S2：基于构建的数据集(X

步骤S201：数据预处理。在本实施例中，为了避免数据差异过大导致的模型训练问题，采用归一化方法将数据处理到0和1之间，处理方式如下：

式中，v为14个变量(11个点的水位以及出力、入库流量、弃水流量)中的任意一个，v

经归一化处理后，数据集(X

步骤S202：数据划分：将数据按照一定的比例随机划分为训练集(TrainSX

步骤S203：模型结构及参数设计。在本实施例中，初始模型结构采用最简单的LSTM模型结构，首先，为了保证输入数据的结构一致，将SX

式中，n为数据个数，Y

步骤S204：模型训练：在给定的模型结构和超参数集条件下，基于模型设定训练集(TrainSX

步骤S205：模型评估：

基于训练好的模型和优化目标，评估在验证集(TestSX

步骤S206：确定模型评估结果是否满足预测精度需求，若否，则调整模型结构、优化目标和超参数，重复执行步骤S204、步骤S205，直至获取符合预测精度需求的预报模型。在本实施例中，根据实际效果，不断优化后确定模型结构如图2所示。其中Input 1和步骤S203步骤中的输入一致，Input 2为未来的出力、流量数据(48×3矩阵)，Outpt 1的结构为48×1的向量，和Input 2的输入拼接成48×4的矩阵作为后续的输入，最终通过TimeDistributed层形成输出。经优化调整后，训练集和验证集的均方误差分别为1.54×10

步骤S3：基于步骤S2步骤中所建立的预报模型，产生水位预测结果。具体包括如下步骤：

步骤S301：在进行预报的时刻T，按照步骤S103、步骤S104步骤中所述方法收集相应数据，构建该时刻对应的样本点的输入X

步骤S302：采用步骤S201中同样的数据预处理方法，将输入X

步骤S303：采用步骤S206步骤中确定的模型，预测得到SYHat

步骤S304：根据步骤S201中采用的数据预处理方法，将SYHat

YHat

其中，v

最终该实施例的预报结果如图3所示。

为了进一步说明本发明所提供的考虑动库容影响的水库超短期水位预测方法的实际应用效果，进一步分析该方法在A电站2020年全年的预测效果。在误差分析时以绝对误差和平均绝对误差为分析指标，平均绝对误差的计算方式如下：

其中，

图4中显示了该方法在2020年全年的整体预测误差随预见期的变化情况，从图中可以看出，整体预测的平均绝对误差随预见期的增长而逐渐增长，在预见期为48小时时，预报的平均绝对误差为0.063m。图5展示了2020年全年预测的绝对误差分布情况，从图中可知，随着预见期的增长，误差较大的点的数量逐渐增多，但均有95％的样本点的绝对误差在0.2m以内。无论是平均绝对误差还是绝对误差的分布，均表明预测精度相对较小，基本能够满足精度要求。

实施例2：

本发明另一方面提供了一种考虑动库容影响的水库超短期水位预测系统，包括：计算机可读存储介质和处理器；

进一步的，所述计算机可读存储介质用于存储可执行指令；

进一步的，所述处理器用于读取所述计算机可读存储介质中存储的可执行指令，执行第一方面所述的考虑动库容影响的水库超短期水位预测方法。

实施例3：

本发明另一方面提供了一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现第一方面所述的考虑动库容影响的水库超短期水位预测方法。