一种台区分布式光伏电站群功率预测方法技术

本发明专利技术公开了一种台区分布式光伏电站群功率预测方法，包括序列分解，低频分量时空预测，高频分量时空预测和序列相加，本发明专利技术不需要气象数据，仅利用光伏出力数据即可达到预测的目的，不需要为分布式光伏额外安装气象量测设备，节省经济成本，同时本发明专利技术的输出为台区内多个分布式光伏电站的出力预测值，能够同时预测多个分布式光伏电站的出力，较单站预测可以更好地学习到多站光伏输出功率之间的时空相关性，提高了区域光伏预测的精度并降低了训练成本，且将原始的出力序列分解为低频和高频分量分别建模，分解后得到的子分量序列能够更好地从不同维度表征原始数据的特性，通过针对性的预测建模可得到更为准确的预测结果

【技术实现步骤摘要】
一种台区分布式光伏电站群功率预测方法


[0001]本专利技术涉及分布式光伏电站进功率预测
，具体为一种台区分布式光伏电站群功率预测方法
。

技术介绍

[0002]现有分布式光伏预测技术主要有插值法
、
参考电站法和历史出力法三类，插值法主要通过对卫星气象数据进行插值，推算分布式站点所在地理位置的气象数据，再利用物理模型法或数据驱动法预测发电功率，参考电站法主要将分布式光伏站点邻近的某个集中式光伏电站作为参考站点，建立出力相关模型，借助参考电站的功率预测值推算分布式光伏出力预测值，历史出力法主要在无气象数据的条件下，仅基于分布式光伏的历史出力数据建立时序预测模型；
[0003]但是当前现有的分布式光伏预测技术中，插值法依赖于卫星气象数据，但卫星气象数据空间分辨率较低，气象插值结果误差较大，导致功率预测结果误差较大，而参考电站法依赖于邻近集中式光伏电站提供的气象数据，但实际中无法确保所有分布式光伏电站周围都有集中式光伏电站，因此适用性有限，历史出力法往往对单个分布式光伏电站进行功率预测，缺乏对多个电站出力之间空间相关性的考虑，预测误差较大，且不同电站之间的模型适应性差
、
无法共用
。

技术实现思路

[0004]本专利技术提供一种台区分布式光伏电站群功率预测方法，可以有效解决上述
技术介绍
中提出当前现有的分布式光伏预测技术中，插值法依赖于卫星气象数据，但卫星气象数据空间分辨率较低，气象插值结果误差较大，导致功率预测结果误差较大，而参考电站法依赖于邻近集中式光伏电站提供的气象数据，但实际中无法确保所有分布式光伏电站周围都有集中式光伏电站，因此适用性有限，历史出力法往往对单个分布式光伏电站进行功率预测，缺乏对多个电站出力之间空间相关性的考虑，预测误差较大，且不同电站之间的模型适应性差
、
无法共用的问题
。
[0005]为实现上述目的，本专利技术提供如下技术方案：一种台区分布式光伏电站群功率预测方法，基于序列分解和门控融合图卷积的手段，充分考虑光伏电站群出力序列不同分量的时空相关性，在无气象信息的条件下，仅基于历史出力数据实现台区内多个光伏电站功率的同时预测，具体包括如下步骤：
[0006]S1、
序列分解，利用极大重叠离散小波变换将原始光伏出力序列分解为低频序列分量和若干高频序列分量；
[0007]S2、
低频分量时空预测，基于静态图卷积网络和
Transformer
模型来预测低频分量；
[0008]S3、
高频分量时空预测，基于门控融合图卷积和
Transformer
模型来预测高频分量；
[0009]S4、
序列相加，将低频预测值和高频预测值相加，从而得到最终的预测结果
。
[0010]根据上述技术方案，所述
S1
中，序列分解具体包括如下实现步骤：
[0011]S101、
模型输入；
[0012]S102、
最大值归一化；
[0013]S103、
小波分解
。
[0014]根据上述技术方案，所述
S101
中，模型输入为台区内
N
个光伏电站
T
个历史时刻的出力
P
i,t
组成的矩阵
P
，其中，
i
＝1，2…
，
N
，
t
＝1，2…
，
T
；
[0015][0016]所述
S102
中，最大值归一化是指将台区内
N
个光伏电站出力序列进行归一化出力，则有
[0017][0018]所述
S103
中，小波分解主要指利用
matlab
小波分析工具箱中
wavedec
函数和
wrcoef
函数将光伏出力序列分解为趋势分量
a3
和波动分量
d1
，
d2
，
d3
，选择“db5”小波，具体程序如下：
[0019][C,L]＝
wavedec(P
i
,3,'db5')
；
[0020]d1
＝
wrcoef('d'C,L,wname,1)
；
[0021]d2
＝
wrcoef('d',C,L,wname,2)
；
[0022]d3
＝
wrcoef('d',C,L
，
,wname,3)
；
[0023]a3
＝
wrcoef('a',C,L,wname,3)。
[0024]根据上述技术方案，所述
S2
中，低频分量时空预测具体包括如下实现步骤：
[0025]S201、
时空特征提取层；
[0026]S202、
预测层
。
[0027]根据上述技术方案，所述
S201
中，时空特征提取层是指构建静态图卷积网络提取低频分量的时空特征，具体算法为：
[0028]首先根据
Haversine
公式求解各个光伏电站的地理距离：
[0029]d
ij
＝
Rarccos(cos
β
i
cos
β
j
cos(
α
i
‑
α
j
)+sin
β
i
sin
β
j
)(i,j
＝
1,2,...,N)
；
[0030]依据光伏电站的地理距离构建静态图结构邻接矩阵：
[0031][0032]然后构建双层
SGCN
提取光伏电站各个时刻的空间特征：
[0033][0034]其中，为自连接的静态邻接矩阵，为邻接矩阵的度矩阵，
P
tA
为低频分量在时刻
t
的值，
W1、W2为权重参数；
[0035]拼接
T
个历史时刻出力空间特征个历史时刻出力空间特征
[0036]将
SGCN
提取的各个时刻的空间特征和原始低频分量拼接后输入
Transformer
的
encoder
结构，具体算法如下：
[0037]X
A
＝
concat(H
A
,P
A
)
[0038]将
X
A
分解为其中
[0039]对应的位置编码：
[0040][0041]其中，为第
t
...

【技术保护点】

【技术特征摘要】
Networks
，
SGCN)
提取低频分量的时空特征，具体算法为：首先根据
Haversine
公式求解各个光伏电站的地理距离：
d
ij
＝
Rarccos(cos
β
i
cos
β
j
cos(
α
i
‑
α
j
)+sin
β
i
sin
β
j
)(i,j
＝
1,2,...,N)
；依据光伏电站的地理距离构建静态图结构邻接矩阵：然后构建双层
SGCN
提取光伏电站各个时刻的空间特征：其中，为自连接的静态邻接矩阵，为邻接矩阵的度矩阵，
P
tA
为低频分量在时刻
t
的值，
W1、W2为权重参数；拼接
T
个历史时刻出力空间特征个历史时刻出力空间特征将
SGCN
提取的各个时刻的空间特征和原始低频分量拼接后输入
Transformer
的
encoder
结构，具体算法如下：
X
A
＝
concat(H
A
,P
A
)
将
X
A
分解为其中其中对应的位置编码：其中，为第
t
行第
j
列对应的位置编码；将与位置编码
E
对位相加即可：则有由作线性映射得到：作线性映射得到：作线性映射得到：通过点乘得到时空特征：
所述
S202
中，预测层是指为缓解过拟合和梯度消失问题，采用残差连接将步骤
S201
中提取的时空特征与原始输入拼接后输入全连接层得到光伏电站群低频分量的预测值：
Transformer
并行输出
N
个电站的低频分量预测值
6.
根据权利要求1所述的一种台区分布式光伏电站群功率预测方法，其特征在于：所述
S3
中，高频分量时空预测具体包括如下实现步骤：
S301、
静态时空特征提取层
S302、
动态时空特征提取层；
S303、
门控融合层；
S304、
全连接层
。7.
根据权利要求6所述的一种台区分布式光伏电站群功率预测方法，其特征在于：所述
S301
中，静态时空特征提取层是指构建双层
SGCN
提取光伏电站各个时刻的空间特征，如下式所示：其中，为自连...

【专利技术属性】
技术研发人员：赵洪山，孙承妍，温开云，郭帅，宋玮琼，李海涛，王硕，李季巍，潘全成，
申请(专利权)人：国网北京市电力公司国家电网有限公司，
类型：发明
国别省市：