本发明涉及光伏发电，具体涉及一种基于cnn-lstm的短期太阳辐射预测方法及装置。
背景技术：
目前经济的快速发展主要依赖石油、煤炭等化石燃料，造成了严重的环境污染和温室效应，故可再生能源的开发利用越来越受到世界范围内的广泛关注。在这些能源中，太阳能是最有前途的选择之一，特别是在光伏发电中。然而，太阳辐照度因受到云层快速变化的影响，具有随机性。当大规模光伏发电并网时，这种随机性会对电网安全运行造成严重影响。因此，建立准确的短期太阳辐照度预测模型对于保证电力系统的最优调度和管理是至关重要的，特别是在光伏发电领域。
现有的太阳辐照度预测模型根据采用的输入数据可分为四类：数值天气预报(nwp)、基于图像的模型、基于历史数据的统计模型和基于多种类型输入数据的混合模型。其中，数值天气预报已广泛应用于提前6-48小时的太阳辐照度预测；然而，这种方法涉及到的计算成本较高，通常不适合进行短期太阳辐照度预测。且nwp模型不能很好地捕捉太阳辐照度与气象因子之间的实时关系，进一步限制了其短期预测能力。与nwp不同，基于图像的模型提供了有效的短期太阳辐照度预测水平，范围从几分钟到几个小时不等。然而，图像捕获设备的高成本限制了该模型的广泛使用。与前两种模型相比，统计模型基于历史太阳辐照度数据，通常用于5分钟到6小时的预测时间尺度。相比而言，基于时间序列数据进行太阳辐照度预测时，统计模型往往是较好的选择，但统计模型只能考虑历史数据，其他影响太阳辐照度的相关天气条件不能包含在预测过程中。这个缺点可以通过机器学习方法在混合模型中集成多种不同类型的输入数据来克服。机器学习方法具有从多个输入中提取复杂非线性特征的优点。人工神经网络和支持向量机等多种机器学习模型在数据样本量较小的情况下具有较好的预测性能，在太阳辐射预测中得到了广泛的应用。然而，传统的机器学习模型只能提取较浅层的特征，且通常需要复杂的特征工程。
技术实现要素：
发明目的：本申请的目的在于提供一种基于cnn-lstm的短期太阳辐射预测方法及装置，利用cnn-lstm模型预测短期太阳辐射，解决现有技术短期预测能力差，成本高，预测精度较低的缺陷。
技术方案：本发明一方面提供了一种基于cnn-lstm的短期太阳辐射预测方法，包括：
(1)将目标站点和邻近站点作为样本站点，采集样本站点的气象参数和辐射数据；
(2)基于所有样本站点的气象参数组成二维矩阵，重构空间特征；
(3)基于所有样本站点的历史太阳辐射序列重构时间特征；
(4)将空间特征和时间特征分别作为cnn-lstm混合模型的cnn部分和lstm部分的输入，预测目标站点太阳总辐射。
进一步地，步骤(2)包括：
(21)将所有样本站点待预测时刻的气象参数构成二维矩阵，每列表示所有站点的一种气象参数，每行表示一个站点的所有气象参数；二维矩阵的元素由下式表示：
fi,j(i＝1,2,…,s；j＝1,2,…,f)
其中，s表示站点数量，f表示气象参数的数量，fi,j表示第i个站点的第j个气象参数；
(22)特定的气象参数在不同的站点间存在第一相关性，具体表现为同一列中相邻两个元素之间的相关性，第一相关性corr1定义如下：
corr1＝g1(fi,j,fi 1,j)
其中g1(·)表示站点相关性函数；
同时，同一站点的不同气象参数之间存在第二相关性，具体表现为同一行相邻两个元素之间的相关性，第二相关性corr2定义如下：
corr2＝g2(fi,j,fi,j 1)
其中g2(·)表示气象参数相关性函数；第一相关性和第二相关性共同构成空间相关性。
进一步地，步骤(3)包括：
(31)将每个时间步的输入设置为所有样本站点在该时间步的辐射值；
(32)将所有样本站点的辐射序列按照时间步依次输入到lstm部分中，以提取历史太阳辐射序列中的时间特征。
进一步地，cnn-lstm混合模型包括cnn部分和lstm部分；将第一输入输入到cnn部分的卷积层中，经过池化层、平坦层和全连接层得到第一特征；将第二输入输入到lstm部分的lstm层中，经过全连接层后得到第二特征；将第一特征和第二特征通过concatenate操作融合，经过全连接层后得到最终的输出。
进一步地，步骤(4)包括：
(41)将重构的空间特征作为cnn部分的输入，以提取空间相关性；
(42)将重构的时间特征作为lstm部分的输入，以提取时间相关性；
(43)将提取的空间相关性和时间相关性相融合，通过神经网络层来预测目标站点太阳总辐射。
另一方面，本申请还提供了一种基于cnn-lstm的短期太阳辐射预测装置，包括：
数据采集模块，用于将目标站点和邻近站点作为样本站点，采集样本站点的气象参数和辐射数据；
空间特征重构模块，用于基于所有样本站点的气象参数组成二维矩阵，重构空间特征；
时间特征重构模块，用于基于所有样本站点的历史太阳辐射序列重构时间特征；
预测模块，用于将空间特征和时间特征分别作为cnn-lstm混合模型的cnn部分和lstm部分的输入，预测目标站点太阳总辐射。
进一步地，空间特征重构模块包括：
二维矩阵构建单元，用于将所有样本站点待预测时刻的气象参数构成二维矩阵，每列表示所有站点的一种气象参数，每行表示一个站点的所有气象参数；二维矩阵的元素由下式表示：
fi,j(i＝1,2,…,s；j＝1,2,…,f)
其中，s表示站点数量，f表示气象参数的数量，fi,j表示第i个站点的第j个气象参数；
相关性定义单元，用于定义第一相关性和第二相关性：特定的气象参数在不同的站点间存在第一相关性，具体表现为同一列中相邻两个元素之间的相关性，第一相关性corr1定义如下：
corr1＝g1(fi,j,fi 1,j)
其中g1(·)表示站点相关性函数；
同时，同一站点的不同气象参数之间存在第二相关性，具体表现为同一行相邻两个元素之间的相关性，第二相关性corr2定义如下：
corr2＝g2(fi,j,fi,j 1)
其中g2(·)表示气象参数相关性函数；第一相关性和第二相关性共同构成空间相关性。
进一步地，时间特征重构模块包括：
输入设置单元，用于将每个时间步的输入设置为所有样本站点在该时间步的辐射值；
时间特征提取单元，用于将所有样本站点的辐射序列按照时间步依次输入到lstm部分中，以提取历史太阳辐射序列中的时间特征。
进一步地，cnn-lstm混合模型包括cnn部分和lstm部分；将第一输入输入到cnn部分的卷积层中，经过池化层、平坦层和全连接层得到第一特征；将第二输入输入到lstm部分的lstm层中，经过全连接层后得到第二特征；将第一特征和第二特征通过concatenate操作融合，经过全连接层后得到最终的输出。
进一步地，预测模块包括：
空间相关性提取单元，用于将重构的空间特征作为cnn部分的输入，以提取空间相关性；
时间相关性提取单元，用于将重构的时间特征作为lstm部分的输入，以提取时间相关性；
融合预测单元，用于将提取的空间相关性和时间相关性相融合，通过神经网络层来预测目标站点太阳总辐射。
有益效果：与现有技术相比，本申请的方法和装置基于目标站点和周边站点的气象参数和辐射数据，重构空间特征和时间特征，输入到cnn-lstm混合模型中，对目标站点的太阳辐射进行预测，提高了太阳辐射预测的精度，保障了电力系统的最优调度和安全稳定运行。
附图说明
图1为本发明的预测方法流程图；
图2为本发明的预测装置的结构框图；
图3为本发明的空间特征重构示意图；
图4为本发明采用的cnn-lstm混合模型结构示意图；
图5为本发明实施例中不同模型的mae误差指标对比；
图6为本发明实施例中不同模型的rmse误差指标对比。
具体实施方式
下面结合附图和实施例对本发明做进一步描述：
实施例一
本发明一方面提供了一种基于cnn-lstm的短期太阳辐射预测方法，如图1所示，包括：
s101将目标站点和邻近站点作为样本站点，采集样本站点的气象参数和辐射数据；
s102基于所有样本站点的气象参数组成二维矩阵，重构空间特征；包括以下步骤：
(21)将所有样本站点待预测时刻的气象参数按照图3形式构成二维矩阵，每列表示所有站点的一种气象参数，每行表示一个站点的所有气象参数；所述二维矩阵的元素由下式表示：
fi,j(i＝1,2,…,s；j＝1,2,…，s；j＝1，2，…,f)
其中，s表示站点数量，f表示气象参数的数量，fi,j表示第i个站点的第j个气象参数；
(22)特定的气象参数在不同的站点间存在第一相关性，具体表现为同一列中相邻两个元素之间的相关性，第一相关性corr1定义如下：
corr1＝g1(fi,j,fi 1,j)
其中g1(·)表示站点相关性函数；
同时，同一站点的不同气象参数之间存在第二相关性，具体表现为同一行相邻两个元素之间的相关性，第二相关性corr2定义如下：
corr2＝g2(fi,j,fi,j 1)
其中g2(·)表示气象参数相关性函数。第一相关性和第二相关性共同构成空间相关性。
s103基于所有样本站点的历史太阳辐射序列重构时间特征；包括以下步骤：
(31)将每个时间步的输入设置为所有样本站点在该时间步的辐射值；
(32)将所有样本站点的辐射序列按照时间步依次输入到lstm部分中，以提取历史太阳辐射序列中的时间特征。
s104将空间特征和时间特征分别作为cnn-lstm混合模型的cnn部分和lstm部分的输入，预测目标站点太阳总辐射；其中，cnn-lstm混合模型结构如图4所示：cnn-lstm混合模型包括cnn部分和lstm部分；将第一输入输入到cnn部分的卷积层中，经过池化层、平坦层和全连接层得到第一特征；将第二输入输入到lstm部分的lstm层中，经过全连接层后得到第二特征；将第一特征和第二特征通过concatenate操作融合，经过全连接层后得到最终的输出；包括以下步骤：
(41)将重构的空间特征作为cnn部分的输入，以提取空间相关性；
(42)将重构的时间特征作为lstm部分的输入，以提取时间相关性；
(43)将提取的空间相关性和时间相关性相融合，通过神经网络层来预测目标站点太阳总辐射。
下面以德克萨斯州为例介绍本发明的方法：
选取德克萨斯州三个气候带(cfa、bsh和bsk)进行测试，其中cfa气候带共有12个站点，选取两个目标站点dallas和sanjacinto，bsh气候带共有8个站点，选取zapata为目标站点，bsk气候带共有14个站点，选取两个目标站点moore和lamb。每个站点的数据包括太阳总辐射、露点温度、太阳天顶角、风速、风向、可降水量、相对湿度和温度。所使用的数据是从2006年1月1日到2012年12月31日的每小时一测的数据。同时，由于辐射数据在夜间降为0，所以只选取每天6点到18点的数据进行训练和测试，提前一小时预测太阳总辐射。训练集为2006年1月1日到2010年12月31日，测试集为2011年1月1日到2012年12月31日，所采用的误差指标为mae、nmae、rmse和nrmse。
此外，除了所提出的cnn-lstm混合模型，选取了svm、ann、cnn、lstm、cnn-ann和ann-lstm六个基准模型进行对比，其中cnn-ann和ann-lstm是将cnn-lstm混合模型中的lstm部分和cnn部分分别换成了ann，测试的结果如表1、图5和图6所示：
表1三个气候带下不同模型的预测结果
从表1可以看出，不论在何种气候带下，不论是哪种误差指标，综合来说，本发明所提出的技术手段cnn-lstm混合模型的预测误差都较小，验证了基于cnn-lstm的短期太阳辐射预测方法的有效性。
图5和图6分别对比了不同模型的mae和rmse指标情况，不难发现，cnn-lstm的区域面积基本上小于其他对比模型的，进一步说明了cnn-lstm的预测效果较好。
实施例二
另一方面，本发明还提供了一种基于cnn-lstm的短期太阳辐射预测装置，如图2所示，包括数据采集模块201、空间特征重构模块202、时间特征重构模块203、预测模块204：
数据采集模块201用于将目标站点和邻近站点作为样本站点，采集样本站点的气象参数和辐射数据；
空间特征重构模块202用于基于所有样本站点的气象参数组成二维矩阵，重构空间特征；空间特征重构模块202包括二维矩阵构建单元和相关性定义单元；其中：
二维矩阵构建单元，用于将所有所述样本站点待预测时刻的气象参数构成二维矩阵，每列表示所有站点的一种气象参数，每行表示一个站点的所有气象参数；所述二维矩阵的元素由下式表示：
fi,j(i＝1,2,…,s；j＝1,2,…,f)
其中，s表示站点数量，f表示气象参数的数量，fi,j表示第i个站点的第j个气象参数；
相关性定义单元，用于定义第一相关性和第二相关性：特定的气象参数在不同的站点间存在第一相关性，具体表现为同一列中相邻两个元素之间的相关性，第一相关性corr1定义如下：
corr1＝g1(fi,j,fi 1,j)
其中g1(·)表示站点相关性函数；
同时，同一站点的不同气象参数之间存在第二相关性，具体表现为同一行相邻两个元素之间的相关性，第二相关性corr2定义如下：
corr2＝g2(fi,j,fi,j 1)
其中g2(·)表示气象参数相关性函数。第一相关性和第二相关性共同构成空间相关性。
时间特征重构模块203用于基于所有样本站点的历史太阳辐射序列重构时间特征；时间特征重构模块203包括输入设置单元和时间特征提取单元，其中：输入设置单元，用于将每个时间步的输入设置为所有样本站点在该时间步的辐射值；时间特征提取单元，用于将所有样本站点的辐射序列按照时间步依次输入到lstm部分中，以提取历史太阳辐射序列中的时间特征。
预测模块204用于将空间特征和时间特征分别作为cnn-lstm混合模型的cnn部分和lstm部分的输入，预测目标站点太阳总辐射。cnn-lstm混合模型结构如图4所示，包括cnn部分和lstm部分；将第一输入输入到cnn部分的卷积层中，经过池化层、平坦层和全连接层得到第一特征；将第二输入输入到lstm部分的lstm层中，经过全连接层后得到第二特征；将第一特征和第二特征通过concatenate操作融合，经过全连接层后得到最终的输出。
预测模块204包括：
空间相关性提取单元，用于将重构的空间特征作为cnn部分的输入，以提取空间相关性；
时间相关性提取单元，用于将重构的时间特征作为lstm部分的输入，以提取时间相关性；
融合预测单元，用于将提取的空间相关性和时间相关性相融合，通过神经网络层来预测目标站点太阳总辐射。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
1.一种基于cnn-lstm的短期太阳辐射预测方法，其特征在于，包括：
(1)将目标站点和邻近站点作为样本站点，采集所述样本站点的气象参数和辐射数据；
(2)基于所有所述样本站点的气象参数组成二维矩阵，重构空间特征；
(3)基于所有所述样本站点的历史太阳辐射序列重构时间特征；
(4)将所述空间特征和所述时间特征分别作为cnn-lstm混合模型的cnn部分和lstm部分的输入，预测目标站点太阳总辐射。
2.根据权利要求1所述的方法，其特征在于，步骤(2)包括：
(21)将所有样本站点待预测时刻的气象参数构成二维矩阵，每列表示所有站点的一种气象参数，每行表示一个站点的所有气象参数；所述二维矩阵的元素由下式表示：
fi,j(i＝1,2,…,s；j＝1,2,…,f)
其中，s表示站点数量，f表示气象参数的数量，fi,j表示第i个站点的第j个气象参数；
(22)特定的气象参数在不同的站点间存在第一相关性，具体表现为同一列中相邻两个元素之间的相关性，第一相关性corr1定义如下：
corr1＝g1(fi,j,fi 1,j)
其中g1(·)表示站点相关性函数；
同时，同一站点的不同气象参数之间存在第二相关性，具体表现为同一行相邻两个元素之间的相关性，第二相关性corr2定义如下：
corr2＝g2(fi,j,fi,j 1)
其中g2(·)表示气象参数相关性函数；第一相关性和第二相关性共同构成空间相关性。
3.根据权利要求2所述的方法，其特征在于，步骤(3)包括：
(31)将每个时间步的输入设置为所有样本站点在该时间步的辐射值；
(32)将所有样本站点的辐射序列按照时间步依次输入到lstm部分中，以提取历史太阳辐射序列中的时间特征。
4.根据权利要求3所述的方法，其特征在于，所述cnn-lstm混合模型包括cnn部分和lstm部分；将第一输入输入到cnn部分的卷积层中，经过池化层、平坦层和全连接层得到第一特征；将第二输入输入到lstm部分的lstm层中，经过全连接层后得到第二特征；将所述第一特征和所述第二特征通过concatenate操作融合，经过全连接层后得到最终的输出。
5.根据权利要求4所述的方法，其特征在于，步骤(4)包括：
(41)将重构的空间特征作为cnn部分的输入，以提取空间相关性；
(42)将重构的时间特征作为lstm部分的输入，以提取时间相关性；
(43)将提取的空间相关性和时间相关性相融合，通过神经网络层来预测目标站点太阳总辐射。
6.一种基于cnn-lstm的短期太阳辐射预测装置，其特征在于，包括：
数据采集模块，用于将目标站点和邻近站点作为样本站点，采集所述样本站点的气象参数和辐射数据；
空间特征重构模块，用于基于所有所述样本站点的气象参数组成二维矩阵，重构空间特征；
时间特征重构模块，用于基于所有所述样本站点的历史太阳辐射序列重构时间特征；
预测模块，用于将所述空间特征和所述时间特征分别作为cnn-lstm混合模型的cnn部分和lstm部分的输入，预测目标站点太阳总辐射。
7.根据权利要求6所述的装置，其特征在于，所述空间特征重构模块包括：
二维矩阵构建单元，用于将所有所述样本站点待预测时刻的气象参数构成二维矩阵，每列表示所有站点的一种气象参数，每行表示一个站点的所有气象参数；所述二维矩阵的元素由下式表示：
fi,j(i＝1,2,…,s；j＝1,2,…,f)
其中，s表示站点数量，f表示气象参数的数量，fi,j表示第i个站点的第j个气象参数；
相关性定义单元，用于定义第一相关性和第二相关性：特定的气象参数在不同的站点间存在第一相关性，具体表现为同一列中相邻两个元素之间的相关性，第一相关性corr1定义如下：
corr1＝g1(fi,j,fi 1,j)
其中g1(·)表示站点相关性函数；
同时，同一站点的不同气象参数之间存在第二相关性，具体表现为同一行相邻两个元素之间的相关性，第二相关性corr2定义如下：
corr2＝g2(fi,j,fi,j 1)
其中g2(·)表示气象参数相关性函数；第一相关性和第二相关性共同构成空间相关性。
8.根据权利要求7所述的装置，其特征在于，所述时间特征重构模块包括：
输入设置单元，用于将每个时间步的输入设置为所有样本站点在该时间步的辐射值；
时间特征提取单元，用于将所有样本站点的辐射序列按照时间步依次输入到lstm部分中，以提取历史太阳辐射序列中的时间特征。
9.根据权利要求8所述的装置，其特征在于，所述cnn-lstm混合模型包括cnn部分和lstm部分；将第一输入输入到cnn部分的卷积层中，经过池化层、平坦层和全连接层得到第一特征；将第二输入输入到lstm部分的lstm层中，经过全连接层后得到第二特征；将所述第一特征和所述第二特征通过concatenate操作融合，经过全连接层后得到最终的输出。
10.根据权利要求9所述的装置，其特征在于，所述预测模块包括：
空间相关性提取单元，用于将重构的空间特征作为cnn部分的输入，以提取空间相关性；
时间相关性提取单元，用于将重构的时间特征作为lstm部分的输入，以提取时间相关性；
融合预测单元，用于将提取的空间相关性和时间相关性相融合，通过神经网络层来预测目标站点太阳总辐射。
技术总结