1.一种基于多尺度时空图变换网络的海上风电功率预测方法,其特征在于,包括以下步骤:(1)采集海上多个风电场的历史数据,并对历史数据进行小样本增强,获取高分辨率历史数据;
(2)构建多尺度时空图变换网络模型MSSTGTN,所述MSSTGTN模型包括多时间尺度自相关特征提取模块,时空特征提取模块以及预测输出模块;
(3)利用多时间尺度自相关特征提取模块对经过数据增强的高分辨率历史数据构造不同预测时间长度的卷积算子,自适应地提取高分辨率历史数据中不同时间尺度的自相关特征;
(4)通过时空特征提取模块中的混合膨胀图卷积网络和变体Transformer网络来提取复杂的时空特征;并引入门控融合机制模块来增强MSSTGTN模型对数据中时空交融性的学习能力;
(5)在MSSTGTN模型中引入无导数优化算法对模型超参数进行动态调整,引入时间与形状失真损失函数作为MSSTGTN模型的附加损失函数,采用多目标贝叶斯算法建立最小化时间与形状失真损失函数的优化目标函数,优化MSSTGTN模型的超参数以及时间与形状失真损失函数的权重;
(6)将所提取到的自相关特征和时空特征通过预测输出模块输出预测风电功率,得到预测结果;
步骤(1)所述历史数据包括历史风电功率数据、数值天气预报系统数据、地理潜在依赖数据序列、邻近数据序列、周期性重复数据序列以及历史气象要素数据序列;
所述步骤(3)实现过程如下:
通过多时间尺度自相关特征提取模块的一维卷积算子自适应提取不同时间尺度的风速数据的自相关特征,作为后续时空特征提取模块的输入特征,其中,一维卷积算子数学公式描述如下:其中,F表示激活函数, 表示输入的高分辨历史数据向量, 表示特征映射的偏差, 是第k个卷积核的权重,L是卷积算子的长度, 是由L和步长S确定的高分辨率历史数据的元素;
所述步骤(4)包括以下步骤:
(41)时空特征提取模块包括四个时空卷积块,分别对地理潜在依赖数据序列、邻近数据序列、周期性重复数据序列以及历史气象要素数据序列进行特征提取,每个时空卷积块都包含一个空间特征提取子模块、一个时间特征提取子模块和一个门控融合机制模块;
(42)利用空间特征提取子模块中的混合膨胀图卷积网络自适应捕获区域内风电机组之间的空间依赖性;
(43)利用时间特征提取子模块中的变体Transformer网络有效捕捉时间序列中不同时间尺度下的时间依赖性;
(44)在时空特征提取模块的门控融合机制模块中,第一层的输入为上一层的输出,然(l)后依次叠加;输入矩阵H 分别进入空间特征提取子模块、时间特征提取子模块以及两个Linear层,然后将从空间和时间两个特征提取子模块输出的结果分别和权重矩阵进行点(l+1)乘,将两个点乘后的结果求和得到输出H :
(l+1) ξ τ T
H =z⊙X+(1‑z)⊙(X) (2)(l) (l) T
z=σ(Linear(H )+Linear((H ))) (3)其中,z表示门控融合机制模块中两个Linear层的比例;
(45)对于每个子模块,在每个网络层的输入和输出之间添加一个恒等映射,即不对原始输入做任何改变,只对输出求和,残差块输出的表达式如下:H(x)=F(x)+x (4)
其中,x表示残差块的输入,H(x)表示残差块的输出,F(x)表示添加的恒等映射;
步骤(42)中所述空间特征提取子模块的具体步骤如下:(421)将区域内风电机组之间的位置信息抽象为图G=(V,E,A),其中V和E分别代表风电机组的节点和节点之间的连接,A是风电机组之间具有欧氏距离的邻接矩阵,矩阵中元素的具体计算式如下:其中,vi和vj代表区域内任意两个节点i和j,(0,1)表示节点i和节点j之间的关联程度值Ai,j的大小在0~1范围内;
N×C
(422)关于t时刻的图G中的输入数据,整个图的信号表示为x∈R ,并且卷积核定义为gθ(L),gθ(L)在傅里叶域中卷积的公式为:T
gθ(L)*x=Ugθ(Ux) (6)T
其中,θ为模型的参数,U为归一化拉普拉斯矩阵的特征向量,Ux表示对x的傅里叶变换;
通过卷积核gθ提取图中每个节点的0到(K-1)阶邻域的信息;最后的卷积公式表示为:其中,为递归公式Tk(x)=2Tk‑1(x)‑Tk‑2(x)的最大特征值,其中T0=1,T1=x;
(423)将得到的图卷积核作用在图G的每一个节点上,捕获到图中的一阶相邻信息;混合膨胀图卷积网络由一层混合膨胀卷积层和一层Add‑Norm共同构成,在混合膨胀图卷积网络之后添加一个前馈神经网络层,空间特征提取子模块的前馈神经网络层表示为:FFN(xm)=max(0,xmW1+b1)W2+b2 (8)N×T
其中,xm∈R 为输入矩阵;
所述步骤(43)包括以下步骤:
(431)在时间维度上进行建模,变体Transformer网络由位置和时间嵌入模块与多头注意力模块两个部分组成;
p
(432)对输入序列使用相对位置嵌入法,生成一个位置矩阵D ,并将其整合到节点特征中,通过与相对位置编码类似的方式将时间嵌入为节点的特征,对采样的绝对时间进行编d h t τ码;将日、小时和分钟分别编码为D、D和D,然后求和作为最终的绝对时间嵌入D:τ d h t
D=D+D+D (9)
p τ
(433)将时间序列数据的编码X、相对位置编码D和绝对时间嵌入D这3种编码求和得到τ一个含有位置信息的矩阵X:
τ p τ
X=X+D+D (10)
(434)将位置和时间嵌入模块的输出作为多头注意力模块的输入;在多头注意力模块中,将一个时间片的空间的所有信息组成一个词向量,学习不同时间片之间的相关性;将输入特征矩阵按照T×N×C的形式输入,通过前馈神经网络映射到具有可学习性的潜在子空间中;多头注意力机制将输入分成同源的3个部分,分别为查询Query、键Key、值Value;并对它们进行线性变换:τ Query τ Key τ Value
Query=XW ,Key=XW ,Value=XW (11)其中, 都为权重矩阵,且
在自注意力机制中,
(435)将3个同源矩阵输入到放缩点积注意力模块中,得到自注意力矩阵:(436)重复(434)到(435),每次只进行一个自注意力操作,每次变换时参数W不同,则整个多头注意力模块表示为:O=MultiH(Query,Key,Value)=Concat(h1,…,hn)W (14)其中,O为多头注意力输出矩阵,n为注意力的头数,hi表示每次进行注意力运算的结果。
2.根据权利要求1所述的基于多尺度时空图变换网络的海上风电功率预测方法,其特征在于,所述步骤(5)实现过程如下:T p
(51)训练样本Train={(Xi,Yi)|i=1,2,..,N},其中Xi=[Xi1,Xi2,...,Xip] ∈R表示T qMSSTGTN模型的输入数据集,Yi=[Yi1,Yi2,...,Yiq] ∈R为预测风电功率数据集的真实值;
时间与形状失真损失函数DILATE表示为:
其中, 为预测风电功率数据集的预测值,γ 为平滑参数,为规整路径,α∈[0,1]为控制平衡的参数, 为欧氏距离方阵,Ω为损失方阵;考虑可靠性和准确性要求的多目标MSSTGTN模型表示为:式中,LMSE和LDILATE表示多目标MSSTGTN模型的目标函数;λ是由多目标MSSTGTN模型的超参数组成的决策变量;S为多目标优化问题的可行域;
(52)初始化算法参数:设置样本点个数n,最大迭代次数K,代理模型GPR(x),采集函数UCB(x|D);其中,代理模型采用高效的高斯过程回归模型GPs,观测到的函数值yl通过如下的形式产生:得到高斯似然函数形式如下:
其中,f(x)服从如下分布:
f(x):GP(m(x),kθ(x,x')) (20)式中,m(x)为均值函数,kθ(x,x')为核函数,θ为核函数的超参数;
采集函数是根据后验分布Pθ(f(·)|D)来构造,选取GPs作为代理模型条件下,其后验分布仍然是高斯分布:Pθ(f(x1:q)|D)=N(μθ(x1:q),∑θ(x1:q)) (21)则采集函数UCB(x|D)表示如下:
其中,γθ(xj)=f(xj)‑μθ(xj),μθ(xj)是后续点的后验均值;
(53)初始化样本点:在可行域S中使用拉丁超立方采样生成n个样本点的集合Xinit={x1,x2,...,xn};
(54)计算适应度值:对Xinit={x1,x2,...,xn}进行适应度值计算得到f(Xinit),初始化评估点集D0={Xinit,f(Xinit)},令t=n,Dt‑1=D0;
(55)根据当前获得的评估点集,构建代理模型GPR(x);
(56)更新评估点:基于代理模型GPR(x),最大化采集函数UCB(x|Dt‑1),获得下一个评估点:xt=argminUCB(x|Dt‑1) (23)(57)更新适应度值:计算评估点xt的适应度值f(xt),并将其添加到当前评估点集合中:Dt=Dt‑1∪{xt,f(xt)} (24)(58)输出最优解:通过非支配排序算法得到当前Pareto最优解集,若t<N,则跳转步骤(55);否则停止迭代输出优化问题的Pareto最优解集;利用最优解集对MSSTGTN模型的超参数进行动态调整。
3.根据权利要求1所述的基于多尺度时空图变换网络的海上风电功率预测方法,其特征在于,步骤(6)所述的预测输出模块为一个两层的全连接前馈神经网络。
4.一种电子设备,其特征在于,包括存储器和处理器,其中:存储器,用于存储能够在处理器上运行的计算机程序;
处理器,用于在运行所述计算机程序时,执行如权利要求1‑3任一项所述的基于多尺度时空图变换网络的海上风电功率预测方法步骤。
5.一种存储介质,其特征在于,所述存储介质上存储有计算机程序,所述计算机程序被至少一个处理器执行时实现如权利要求1‑3任一项所述的基于多尺度时空图变换网络的海上风电功率预测方法步骤。