CN116796622A - 一种基于条件重要性加权自编码器的风电出力场景生成方法 - Google Patents
一种基于条件重要性加权自编码器的风电出力场景生成方法 Download PDFInfo
- Publication number
- CN116796622A CN116796622A CN202310301809.7A CN202310301809A CN116796622A CN 116796622 A CN116796622 A CN 116796622A CN 202310301809 A CN202310301809 A CN 202310301809A CN 116796622 A CN116796622 A CN 116796622A
- Authority
- CN
- China
- Prior art keywords
- data
- wind power
- encoder
- power output
- scene
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005070 sampling Methods 0.000 claims abstract description 22
- 238000012549 training Methods 0.000 claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 238000010606 normalization Methods 0.000 claims description 18
- 230000004913 activation Effects 0.000 claims description 9
- 238000005315 distribution function Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 230000017105 transposition Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000013528 artificial neural network Methods 0.000 claims description 3
- 238000005311 autocorrelation function Methods 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013527 convolutional neural network Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000010485 coping Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
Landscapes
- Wind Motors (AREA)
Abstract
本发明公开了一种基于条件重要性加权自编码器的风电出力场景生成方法,包括以下步骤:获取所在地区的风电出力的历史实测数据并进行数据集划分;对数据进行预处理;搭建条件重要性加权自编码器场景生成模型;利用训练集数据训练条件重要性加权自编码器,直到损失函数收敛;生成新的风电出力场景,评估场景生成模型性能,验证模型的有效性。本方法在风电出力场景生成过程中,通过对每个样本进行多次隐变量采样,获得更准确的对数似然估计值,从而更准确地建模数据分布,最终提高生成场景的质量。
Description
技术领域
本发明涉及风电出力场景生成领域,具体涉及一种基于条件重要性加权自编码器的风电出力场景生成方法。
背景技术
以风能等为代表的可再生能源资源总量丰富,不产生污染。大规模的风电并网将给电网的稳定运行带来巨大的挑战,风能往往具有较强的随机性和波动性,从而导致风力发电出力不确定,因此如何描述风电出力的不确定性问题十分关键。
场景生成是一种通过构建确定性场景来分析电力***不确定性问题的方法,是应对风力发电不确定性的主要方法之一,现有方法主要是基于概率模型的方法,采用概率假设建立风电出力的数学模型,依靠抽样产生场景。但概率模型过度依赖人为假设的概率密度函数,且难以处理复杂的非线性数据,具有较大的局限性。与此相比人工智能技术无需进行统计假设,并且具有更强大的处理复杂非线性数据和容错能力。随着计算机与人工智能技术的发展,以数据驱动的场景生成模型开始广泛应用于电力***处理不确定性问题,其中典型的就是变分自编码器与生成对抗网络。这类方法能根据历史数据生成新的场景,而无需明确指定模型或拟合概率分布,其次,还能给数据制定标签生成特定情况下的场景。生成对抗网络还存在训练过程不稳定,损失函数剧烈抖动以及难以收敛的问题,这些问题目前还没有有效的解决方法。变分自编码器通过最大化变分下界函数来逼近样本的对数似然函数,但无法准确计算出实际的似然函数,因此生成效果还有待提高。重要性加权自编码器是变分自编码器的变体结构,优化目标也是最大化模型的变分下界函数,使用多个采样路径来计算潜在变量的后验分布,并应用重要性加权得出更紧凑的变分下界,更加接近真实的对数似然,最终提高生成效果,但是目前还没有应用于风电场景生成领域的技术方案。
发明内容
本发明的目的在于针对传统变分自编码器在风电历史数据中挖掘其统计规律、获取样本后验分布和泛化能力方面性能还有待提高的问题,提出一种更有效,可靠的场景生成方法,一种基于条件重要性加权自编码器的风电出力场景生成方法。
一种基于条件重要性加权自编码器的风电出力场景生成方法,其特征在于,包括以下步骤:
I.获取所在地区的风电出力的历史实测数据并进行数据集划分;
II.对数据进行预处理;
III.搭建条件重要性加权自编码器场景生成模型;
IV.利用训练集数据训练条件重要性加权自编码器场景生成模型,直到损失函数收敛;
V.生成新的风电出力场景,评估场景生成模型性能,验证模型的有效性;
优选地,所述步骤I中具体包括:
1)采用的历史实测数据结构为10min采样间隔,即每条功率曲线的时间分辨率为10min,每日样本包括144个采样点;
2)数据随机按照80%,20%的比例划分为训练集和测试集;
优选地,所述步骤II中具体包括:
1)对历史风电出力数据进行数据预处理,对数据进行清洗,剔除偏差较大的异常数据;
2)对历史数据进行归一化处理,采用归一化方法公式如下:
其中,x为原始数据,x*为归一化后的数据,xmin为历史数据最小值,xmax为历史数据最大值;
3)将每日的一维历史数据格式转换成二维矩阵格式;
优选地,所述步骤III中具体包括:
1)构建一个编码器将输入数据x映射到低维空间的隐变量z,编码器由四个卷积层和两个全连接层组成,每个卷积层后都有一个批量归一化层和一个ReLU激活函数;
其中,ReLU激活函数为:
f(x)=max(0,x) (2)
2)构建一个解码器将隐变量还原为重构样本解码器由两个全连接层和四个转置卷积层组成,每个转置卷积层中都有一个批量归一化层和一个ReLU激活函数,内部结构与编码器相互对称;
优选地,所述步骤IV中具体包括:
1)确定条件信息c,将c作为标签值赋予训练集;
2)在训练集中随机采样选取一组历史出力数据,作为编码器的输入数据;
3)构建并计算损失函数;
条件重要性加权自编码器优化目标为最大化模型的变分下界函数,用极大似然法优化目标函数可得到对数似然函数logP(x|c)为:
logP(x|c)=DKL(q(z|x,c)||p(z|x,c))+L(x,c) (3)
其中c为额外条件信息,DKL[.]表示Kullback-Leibler散度,q(z|x,c)为近似后验分布,p(z|x,c)为真实后验分布,L(x,c)为似然函数的变分下界函数;
在近似后验分布中进行K次采样得到K个平均权重的隐变量z1,z2...zk,从而适应更多形式的分布,提高生成网络的生成能力,变分下界L(x,c)表示为:
其中Eq(z|x)[.]表示数学期望,p(z|c)和p(x|z,c)表示概率分布函数;
得到损失函数LCIWAE:
LCIWAE=-E[logP(x|c)]+DKL(q(z|x,c)||p(z|x,c)) (5)
其中E[.]表示数学期望;
4)利用梯度计算计算损失函数来获取编码器和解码器参数;
梯度计算公式为:
其中wi神经网络的重要性权重;为归一化后的重要性权重;K为采样次数;
5)重复步骤2),3),4),直到损失函数收敛;
6)提取训练好的条件重要性加权自编码器用于风电出力场景生成任务;
优选地,所述步骤V中具体包括:
1)从测试集中随机挑选一组历史出力数据作为输入,生成新的风电出力场景;
2)对新生成的风电场景进行反归一化,反归一化公式为:
x=x*(xmax-xmin)+xmin (7)
其中,x为反归一化后的数据,x*为归一化后的数据,xmin为历史数据最小值,xmax为历史数据最大值;
3)将新生成的风电出力场景与测试集中实际场景进行对比,验证生成场景的有效性,采用自相关函数R(τ)来评估实际场景和生成场景的时间相关性,其数学表达式为:
其中E[.]表示期望;τ为时间延迟;σx为风电出力xt的标准差;u为风电出力xt的均值;
采用连续分级概率评分CRPS来评估生成场景的综合性能,其数学表达式为:
其中,F(x)为生成场景的累计分布函数,P(x)实际场景的累计分布函数。
附图说明
图1为本发明中条件重要性加权自编码器的风电出力场景生成方法流程示意图
图2本发明中历史风电出力曲线格式转换示意图;
图3本发明中条件重要性加权自编码器结构示意图;
图4本发明中条件重要性加权自编码器优化目标示意图;
图5本发明中条件重要性加权自编码器隐变量采样示意图。
具体实施方式
下面结合附图对本发明进行详细说明,本部分的描述仅是示范性和解释性,不应对本发明的保护范围有任何的限制作用。此外,本领域技术人员根据本文件的描述,可以对本文件中实施例中以及不同实施例中的特征进行相应组合。
一种基于条件重要性加权自编码器的风电出力场景生成方法,其流程示意图如图1所示,包括以下步骤:
I.获取所在地区的风电出力的历史实测数据并进行数据集划分;
II.对数据进行预处理;
III.搭建条件重要性加权自编码器场景生成模型;
IV.利用训练集数据训练条件重要性加权自编码器,直到损失函数收敛;
V.生成新的风电出力场景,评估场景生成模型性能,验证模型的有效性;
所述步骤I中具体包括:
1)采用的历史实测数据结构为10min采样间隔,即每条功率曲线的时间分辨率为10min,每日样本包括144个采样点;
2)数据随机按照80%,20%的比例划分为训练集和测试集;
所述步骤II中具体包括:
1)对历史风电出力数据进行数据预处理,对数据进行清洗,剔除偏差较大的异常数据;
2)对历史数据进行归一化处理,采用归一化方法公式如下:
其中,x为原始数据,x*为归一化后的数据,xmin为历史数据最小值,xmax为历史数据最大值;
3)将每日的一维历史数据格式转换成二维矩阵作为模型的输入,考虑到编码器网络使用的卷积神经网络层,而卷积神经网络更适合处理具有相同行数和列数的矩阵输入,而风电出力曲线为1×n的时间序列,可以通过格式转换为二维矩阵再作为输入数据,格式转换示意图如图2示;
所述步骤III中具体包括:
1)构建编码器由四个卷积层和两个全连接层组成,每个卷积层后都有一个批量归一化层和一个ReLU激活函数;
其中,ReLU激活函数为:
f(x)=max(0,x) (2)
2)构建解码器由两个全连接层和四个转置卷积层组成,每个转置卷积层中都有一个批量归一化层和一个ReLU激活函数,内部结构与编码器相互对称,条件重要性加权自编码器整体结构图如图3所示;
所述步骤IV中具体包括:
1)确定条件信息c,将c作为标签值赋予训练集,常用的条件信息有年份,月份,天气等;
2)在训练集中随机采样选取一组历史出力数据,作为编码器的输入数据;
3)构建并计算损失函数;
条件重要性自编码器优化目标为最大化模型的变分下界函数,优化原理示意图如图4示,到对数似然函数为:
log P(x|c)=DKL(q(z|x,c)||p(z|x,c))+L(x,c) (3)
其中c为额外条件信息,DKL[.]表示Kullback-Leibler散度,q(z|x,c)为近似后验分布,p(z|x,c)为真实后验分布,L(x,c)为似然函数的变分下界函数;
图5所示为条件重要性自编码器隐变量采样过程图,近似后验分布中进行K次采样得到K个平均权重的隐变量z1,z2...zk,从而适应更多形式的分布,提高生成网络的生成能力,L(x,c)表示为:
其中Eq(z|x)[.]表示数学期望,p(z)和p(x|z,c)表示概率分布函数;
得到损失函数:
LCIWAE=-E[logP(x|c)]+DKL(q(z|x,c)||p(z|x,c)) (5)
其中E[.]表示数学期望;
4)利用梯度计算计算损失函数来获取编码器和解码器参数;
梯度计算公式为:
其中wi神经网络的重要性权重;为归一化后的重要性权重;K为采样次数;
5)重复步骤2),3),4),直到损失函数收敛;
6)提取训练好的条件重要性加权自编码器用于风电出力场景生成任务;
所述步骤VI中具体包括:
1)从测试集中随机挑选一组历史出力数据作为输入,生成新的风电出力场景;
2)对新生成的风电场景进行反归一化,反归一化公式为:
x=x*(xmax-xmin)+xmin (7)
其中,x为反归一化后的数据,x*为归一化后的数据,xmin为历史数据最小值,xmax为历史数据最大值;
3)将新生成的风电出力场景与测试集中实际场景进行对比,验证生成场景的有效性,采用自相关函数R(τ)来评估实际场景和生成场景的时间相关性,其数学表达式为:
其中E[.]表示期望;τ为时间延迟;σx为风电出力xt的标准差;u为风电出力xt的均值;
采用连续等级概率分数CRPS来评估生成场景的综合性能,其数学表达式为:
其中,F(x)为生成场景的累计分布函数,P(x)实际场景的累计分布函数。
与现有的变分自编码器风电出力场景生成方法相比,本发明有以下优点:
(1)本发明在场景生成过程中,通过对每个样本进行多次隐变量采样,可以获得更准确的对数似然估计值,从而生成更高质量的样本,提高生成效果;
(2)更准确地建模数据分布的多模式性和异构性,增加了模型灵活性,可以捕捉更复杂的数据分布;
(3)引入重要性加权思路,对采样的隐变量进行加权平均提高了鲁棒性,更适合处理不完整或噪声数据。
这里所述的实施方式是为了帮助读者理解本发明,应被理解为本发明的保护范围并不局限于这样的特别陈述;本领域的技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其他具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (6)
1.一种基于条件重要性加权自编码器的风电出力场景生成方法,其特征在于,包括以下步骤:
I.获取所在地区的风电出力的历史实测数据并进行数据集划分;
II.对数据进行预处理;
III.搭建条件重要性加权自编码器场景生成模型;
IV.利用训练集数据训练条件重要性加权自编码器场景生成模型,直到损失函数收敛;
V.生成新的风电出力场景,评估场景生成模型性能,验证模型的有效性。
2.根据权利要求1所述,其特征在于,步骤I中具体包括:
1)采用的历史实测数据结构为10min采样间隔,即每条功率曲线的时间分辨率为10min,每日样本包括144个采样点;
2)数据随机按照80%,20%的比例划分为训练集和测试集。
3.根据权利要求1所述,其特征在于,步骤II中具体包括:
1)对历史风电出力数据进行数据预处理,对数据进行清洗,剔除偏差较大的异常数据;
2)对历史数据进行归一化处理,采用归一化方法公式如下:
其中,x为原始数据,x*为归一化后的数据,xmin为历史数据最小值,xmax为历史数据最大值;
3)将每日的一维历史数据格式转换成二维矩阵格式。
4.根据权利要求1所述,其特征在于,步骤III中具体包括:
1)构建一个编码器将输入数据x映射到低维空间的隐变量z,编码器由四个卷积层和两个全连接层组成,每个卷积层后都有一个批量归一化层和一个ReLU激活函数;
其中,ReLU激活函数为:
f(x)=max(0,x) (2)
2)构建一个解码器将隐变量还原为重构样本解码器由两个全连接层和四个转置卷积层组成,每个转置卷积层中都有一个批量归一化层和一个ReLU激活函数,内部结构与编码器相互对称。
5.根据权利要求1所述,其特征在于,步骤IV中具体包括:
1)确定条件信息c,将c作为标签值赋予训练集;
2)在训练集中随机采样选取一组历史出力数据,作为编码器的输入数据;
3)构建并计算损失函数;
条件重要性加权自编码器优化目标为最大化模型的变分下界函数,用极大似然法优化目标函数可得到对数似然函数logP(x|c)为:
logP(x|c)=DKL(q(z|x,c)||p(z|x,c))+L(x,c) (3)
其中c为额外条件信息,DKL[.]表示Kullback-Leibler散度,q(z|x,c)为近似后验分布,p(z|x,c)为真实后验分布,L(x,c)为似然函数的变分下界函数;
在近似后验分布中进行K次采样得到K个平均权重的隐变量z1,z2...zk,从而适应更多形式的分布,提高生成网络的生成能力,变分下界L(x,c)表示为:
其中Eq(z|x)[.]表示数学期望,p(z|c)和p(x|z,c)表示概率分布函数;
得到损失函数LCIWAE:
LCIWAE=-E[logP(x|c)]+DKL(q(z|x,c)||p(z||p(z|x|x,c)) (5)
其中E[.]表示数学期望;
4)利用梯度计算计算损失函数来获取编码器和解码器参数;
梯度计算公式为:
其中wi神经网络的重要性权重;为归一化后的重要性权重;K为采样次数;
5)重复步骤2),3),4),直到损失函数收敛;
6)提取训练好的条件重要性加权自编码器用于风电出力场景生成任务。
6.根据权利要求1所述,其特征在于,步骤V中具体包括:
1)从测试集中随机挑选一组历史出力数据作为输入,生成新的风电出力场景;
2)对新生成的风电场景进行反归一化,反归一化公式为:
x=x*(xmax-xmin)+xmin (7)
其中,x为反归一化后的数据,x*为归一化后的数据,xmin为历史数据最小值,xmax为历史数据最大值;
3)将新生成的风电出力场景与测试集中实际场景进行对比,验证生成场景的有效性,采用自相关函数R(τ)来评估实际场景和生成场景的时间相关性,其数学表达式为:
其中E[.]表示期望;τ为时间延迟;σx为风电出力xt的标准差;u为风电出力xt的均值;
采用连续分级概率评分CRPS来评估生成场景的综合性能,其数学表达式为:
其中,F(x)为生成场景的累计分布函数,P(x)实际场景的累计分布函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310301809.7A CN116796622A (zh) | 2023-03-27 | 2023-03-27 | 一种基于条件重要性加权自编码器的风电出力场景生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310301809.7A CN116796622A (zh) | 2023-03-27 | 2023-03-27 | 一种基于条件重要性加权自编码器的风电出力场景生成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116796622A true CN116796622A (zh) | 2023-09-22 |
Family
ID=88042836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310301809.7A Pending CN116796622A (zh) | 2023-03-27 | 2023-03-27 | 一种基于条件重要性加权自编码器的风电出力场景生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116796622A (zh) |
-
2023
- 2023-03-27 CN CN202310301809.7A patent/CN116796622A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111860982B (zh) | 一种基于vmd-fcm-gru的风电场短期风电功率预测方法 | |
CN110070226B (zh) | 基于卷积神经网络与元学习的光伏功率预测方法及*** | |
CN111784041B (zh) | 一种基于图卷积神经网络的风电功率预测方法及*** | |
CN112910690A (zh) | 基于神经网络模型的网络流量预测方法、装置及设备 | |
CN111507046B (zh) | 一种电动闸阀剩余使用寿命预测方法及*** | |
CN111242351A (zh) | 基于自编码器和gru神经网络的热带气旋轨迹预测方法 | |
CN112949904B (zh) | 一种基于特征选择与多层级深度迁移学习的风电场短期功率预测方法 | |
CN116307212B (zh) | 一种空气质量预测方法及*** | |
CN112801332B (zh) | 一种基于灰度共生矩阵的短期风速预测方法 | |
CN112232488A (zh) | 基于深度学习与数据驱动的新能源出力场景生成方法 | |
CN116739172B (zh) | 一种基于爬坡识别的海上风电功率超短期预测方法及装置 | |
CN109583588A (zh) | 一种短期风速预测方法及*** | |
CN115170874A (zh) | 一种基于解耦蒸馏损失的自蒸馏实现方法 | |
CN112418504B (zh) | 一种基于混合变量选择优化深度信念网络风速预测方法 | |
CN114186412A (zh) | 基于自注意力机制的水电站水轮机顶盖长序列水位预测***及方法 | |
CN110555566B (zh) | 一种基于b样条分位数回归的光电概率密度预测方法 | |
CN115604131B (zh) | 一种链路流量预测方法、***、电子设备及介质 | |
CN117354172A (zh) | 网络流量预测方法及*** | |
CN116796622A (zh) | 一种基于条件重要性加权自编码器的风电出力场景生成方法 | |
CN116578858A (zh) | 基于图神经网络的空压机故障预测与健康度评价方法及*** | |
CN115861779A (zh) | 一种基于有效特征表示的无偏场景图生成方法 | |
CN111967577B (zh) | 一种基于变分自编码器的能源互联网场景生成方法 | |
CN114004405A (zh) | 基于Elman神经网络和卫星云图的光伏功率预测方法及*** | |
CN113537573A (zh) | 基于双重时空特征提取的风电功率运行趋势预测方法 | |
Almohri et al. | Data Analysis and Prediction of Power Generated by Photovoltaic Systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |