[发明专利]一种基于时空神经网络的风电集群功率日前预测方法

权利要求书

1.一种基于时空神经网络的风电集群功率日前预测方法，其特征在于：它包括以下步骤：

1)风电集群的数值天气预报(NWP)特征图的构建

数值天气预报(NWP)信息中包含温度、动量通量、各高度风向、各高度风速、湿度，考虑风电集群大量数据的时空关联程度，将风电集群的同一时刻数值天气预报(NWP)信息数据排列到一起，形成空间特征图，设初始时刻为t，训练长度为n个时刻，在各个时刻形成具有丰富空间结构的特征图，n张空间特征图组成时空特征图，因此，风电场之间丰富的时空相关信息被包含在时序特征图中；

2)利用卷积神经网络(CNN)进行特征提取

①输入层：通过步骤1)对风电集群的数值天气预报(NWP)数据进行特征图构建，每一时刻的特征图水平方向代表风电集群中不同风电场的装机容量，竖直方向表示不同类型的数值天气预报(NWP)数据，将多个时刻的特征图按照时间推移顺次排列，作为时空特征输入；

②卷积层：定义一组固定大小的卷积核函数，令各个卷积核与输入数据进行卷积计算，然后加上偏置值，最后经过激活函数，便形成了这一层的特征图，多个不同的卷积核函数经过卷积运算后得到输入数据的不同特征图，体现为风电场范围内受地形、粗糙度、尾流局地因素影响的实际风波动过程的时空耦合关系特征，卷积过程用式(1)表示：

式中，表示第l层的第k个特征图，为偏置矩阵，M_k表示输入特征图的集合，表示第l-1层第i个神经元的输出，为卷积核矩阵，f表示激活函数，符号“*”表示卷积运算；

③池化层：通过使用某一位置相邻区域的总体统计特征作为网络在该位置的输出来减少特征图的尺寸和网络的参数量，同时能够有效避免网络过拟合，最大值池化过程用式(2)表示：

式中，表示第l层的第k个特征图，表示第l-1层的第k个特征图，H为卷积核宽度；

④全连接层：经过特征提取的信号进行分类回归识别，通过激活函数和偏置对输入进行非线性变换，连接层用式(3)表示：

x^l＝f(w^lx^l-1+b^l) (3)

式中，x^l-1表示第l-1层的特征图，b^l为偏置矩阵，w^l为全连接层的权重系数；

通过式(1)和式(2)进行两次训练，式(1)两次卷积核大小定义为3*3，卷积核数量定义为50,100，将得到的结果进行式(3)计算；

3)长短时神经网络(LSTM)的训练

由步骤2)得到每一时刻的特征信息，进而对每一时刻提取到的特征进行预测，首先由遗忘门的Sigmoid层决定从细胞状态中丢弃了什么信息，所述细胞状态包括了风电集群数值天气预报(NWP)的基本属性特征，因此能够预测出未来功率，当输入新的运行数据，就会希望忘记比较旧的数值天气预报(NWP)分布信息，从细胞中丢弃旧的信息，这一步决定通过遗忘门层来完成，遗忘门计算过程为式(4)：

f_t＝sigmoid(W_f·[h_t-1,x_t]+b_f) (4)

式中：f_t是遗忘门的输出，x_t是输入序列，W_f是权值矩阵，h_t-1是上一时刻细胞单元最终的输出，[h_t-1,x_t]表示两个向量连接成一个长向量，b_f是偏移项，通过sigmoid函数后输出0到1的概率；

下一步决定在细胞状态中存储了什么信息，首先，称为“输入门层”的Sigmoid层决定了更新哪些值用式(5)表示，接下来一个tanh层创建候选向量Ct，该向量将会被加到细胞的状态中用式(6)表示，在下一步中，结合这两个向量来创建更新值用式(7)表示，

i_t＝sigmoid(W_i·[h_t-1,x_t]+b_i) (5)

更新细胞状态

C_t＝i_t⊙C_t+f_t⊙C_t-1 (7)

此输出将基于细胞状态，但将是一个过滤版本，首先，运行一个Sigmoid层来决定要输出的细胞状态的部分，然后，将单元格状态通过tanh，将值规范化到-1和1之间，并将其乘以Sigmoid门的输出，两步决定在细胞状态中存储输入用式(8)、式(9)表示，

o_t＝sigmoid(W_o·[h_t-1,x_t]+b_o) (8)

h_t＝o_t·tanh(C_t) (9)

式中:i_t是输入门的输出，作为当前层的候选值可能会添加到单元状态中去，C_t是当前的记忆单元状态，整个过程是对上一时刻的记忆单元状态进行更新，即丢掉无用信息，添加新信息的过程，o_t是输出门的输出,h_t是LSTM当前时刻的最终输出；

4)仿真计算

仿真输入量：对风电集群的实测数据进行分析，确定电场的总装机容量；输入历史数据：各季度的前两个月历史功率，输入预测月份的数值天气预报(NWP)数据；数据采样间隔为15min，根据步骤1)至步骤3)，从而得到最终的功率日前预测结果；

5)评价指标

设P_mi为i时段的实际平均功率，P_pi为i时段的预测功率，C_i为i时段的开机总容量，n为所有样本个数那么平均绝对误差定义为式(10)：

均方根误差定义为式(11)：

按步骤4)，输入仿真输入量，将模型计算的预测功率同实测功率经过步骤5)中的平均绝对误差式(10)、均方根误差式(11)进行误差计算，求得预测准确率。