一种结合气象因素的超短期风电场功率预测方法技术

一种结合气象因素的超短期风电场功率预测方法，属于电力系统发电功率预测技术领域，包括如下步骤：步骤一、对风电功率历史数据、NWP气象数据进行预处理，补齐缺失数据和修改异常数据；步骤二、生成风电功率预测模型；步骤三、利用训练好的模型和未来NWP气象数据进行未来的风电功率预测。本发明专利技术在迭代过程中使用混沌策略周期性地向种群中添加新的萤火虫个体，从而提升超短期风电功率的预测精度，为电网部门的调度人员对短期的调度决策安排提供有利的依据。

【技术实现步骤摘要】
一种结合气象因素的超短期风电场功率预测方法
本专利技术属于电力系统发电功率预测
，具体涉及一种结合气象因素的超短期风电场功率预测方法。
技术介绍
近年来，我国风力发电发展迅速，截止2015年底，我国风电总装机容量已经达到1.283亿kW,连续4年居世界第一，到2020年我国风电装机有望达2.4亿kW以上。与水电、火电、核电等能源相比，风力发电主要由风速、风向、气压、温度、湿度等自然条件决定，具有间歇性、波动性、随机性的特点。大规模风电并网对电网的安全稳定运行带来了挑战，风电功率预测通过气象预报数据、风电场历史运行数据、风电场运行状态数据等参数，预测风电出力变化趋势，为电网安全、电力调度、电力运营带来了积极的影响，在考虑气象因素影响的条件下，单风电场的风电功率预测大多以该风电场SCADA提供的历史风电功率数据和当地NWP系统提供的气象数据为预测模型的训练样本数据，一般会在模型的建立上做出相应的改进和优化，以达到提升风电功率预测精度的目的，超短期风电功率预测的准确性关乎着电网调度部门的工作人员的决策部署，这对提高电网运行的稳定性和经济型都有着重要的意义。目前，我国风电功率预测的准确性仍有着很大的上升空间。门控循环单元神经网络隶属于循环递归神经网络家族，是根据长短期神经网络改进出的一种适用于短时间预测模型构建的神经网络。它对比常用的BP神经网络，GRU网络的收敛速度更快且不容易发生像BP神经网络在误差反馈过程中出现的梯度消失和梯度爆炸的情况，主要原因是该网络复杂的门控单元(细胞状态)通过偏置和相关参数来实时和灵活的调整激活函数的输入，进而控制门的开启和关闭，保存和更新有用的信息，将过程中的梯度进行了抵消，使得GRU深度门控网络具有高效的性能，适合于处理时间序列中间隔和延迟非常长的事件，如预测风速，入库径流，未来天气状况、温度湿度等。萤火虫算法(FireflyAlgorithm，FA)是剑桥学者Yang于2008年提出的一种新颖的进化算法，随着对FA研究的不断深入，目前该算法已在数值优化、工程技术、资源管理等领域得到应用并取得了良好效果。相比于粒子群算法，萤火虫算法能够将种群细分为各个子群体，避免所有个体趋于一个最优值，从而陷入局部最优解。但其算法本身仍然具备其他种群优化算法的缺点。
技术实现思路
针对上述技术问题，本专利技术的提出种一种结合气象因素的超短期风电场功率预测方法，以提高当前风电场所面临的超短期风电功率预测精度不高的问题，利于其风电场参与风电并网调度，为电网部门的工作人员提供做出调度决策的准确依据。本专利技术采用如下技术方案：一种结合气象因素的超短期风电场功率预测方法，包括如下步骤：步骤一、对风电功率历史数据、NWP气象数据进行预处理，补齐缺失数据和修改异常数据；步骤二、生成风电功率预测模型，包括以下步骤：(1)利用pearson系数分析NWP数据和历史功率数据之间的相关性，对NWP数据和历史功率数据归一化处理，提取相关性强的气象数据作为模型的输入，用于训练预测模型；利用pearson系数对NWP系统提供的气象数据和风电场SCADA数据采集系统提供的历史功率数据做相关性分析，提取出的气象因素包括风速、温度、大气密度、湿度分别构成与历史功率数据具有相同时间步长的向量X和Y,然后利用pearson系数公式：对二者进行相关性分析，其中相关系数在(0.8-1.0)和(0.6-0.8)的气象因素为相关性强的气象因素，然后对提取出的气象因素与风电功率历史数据进行归一化处理，其归一化表达式为：式中xmax和xmin分别为该变量的极大值和极小值；经过预测模型得到的风电功率预测数据，再进行反归一化处理使其具有物理意义，反归一化的表达式为：x＝0.5[x'(xmax-xmin)+(xmax-xmin)](3)分析气象因素与风电功率的相关性能够更好的利用NWP数值天气预报的气象数据；(2)门控循环单元神经网络作为风电功率预测的预测模型，隐藏层信息传递是通过一个更新门u(t)和一个重置门r(t)来实现的，在t时刻GRU接收当前状态x(t)、上一时刻的隐藏状态h(t-1)，在信息传递过程中，更新门同时进行遗忘和记忆2个步骤，运算决定选择多少信息输入网络，遗忘多少过去的信息；重置门决定忘记过去的信息量；GRU网络隐藏层的输出h(t)最终由更新门和重置门的动态控制形成。GRU各个变量间的计算公式如下：u(t)＝ψ(Wuh(t-1)+Uux(t))(4)r(t)＝ψ(Wrh(t-1)+Urx(t))(5)式中：Wu、Wr、Wc为连接隐藏层输出信号的参数矩阵；Uu、Ur、Uc为连接输入信号的参数矩阵；ψ(·)为激活函数，为sigmoid激活函数；为任何非线性函数是双曲正切函数tanh；为混合了结合了x(t)和r(t)中间记忆状态；⊙为逻辑运算符，表示矩阵中对应的元素相乘。隐含层t时刻的h(t)将会作为t时刻的信息乘上相应的权值并传入GRU的输出层y'i＝σ(ωih(t))+bi(8)对yi进行反归一化处理yi＝0.5[y′i(y'max-y'min)+(y'max+y'min)](9)利用已有的历史功率数据与预测值yi构建误差目标函数(3)基于混沌理论的萤火虫群体优化算法，对由门控循环单元神经网络得到的预测值和历史值构成的误差修正函数进行寻优，得到最优解对应的神经网络参数，生成预测模型；其算法的主要步骤为：步骤1初始化每个萤火虫的最大吸引度β0，吸收系数γ，随机步长大小α，迭代次数MaxT，目标函数自变量个数8。步骤2首先随机产生Wu、Wr、Wc、Uu、Ur、Uc、ωi、bi各个参数，其矩阵维数要根据具体的预测数据维度和门控循环单元隐含层公式来设计，采用Tent映射产生混沌序列其表达式为：对萤火虫的位置进行初始化，萤火虫初始位置即为门控循环单元神经网络的初始权值，对萤火虫个体搜索空间定义上下限Ud,Dd，并把混沌变量映射到搜索空间中步骤3根据(13)确定萤火虫之间的距离β(γ)＝β0*exp(-γRj)(14)Xi(t+1)＝Xi(t)+β0*exp(-γRj)*(Xj(t)-Xi(t))+α*(r1-0.5)*XM(15)XB(t+1)＝XB(t)+α(r2-0.5)*XM(16)式(14)更新萤火虫吸引度，式(15)更新萤火虫的位置，式(16)确定最优个体的萤火虫随机飞行。步骤4利用(17)和(18)来更新萤火虫的最优个体NG＝G*(1+Cauchy(σ)*range)(17)步骤5若迭代次数等于Np，则开始生成新的萤火虫个体。1)采用式(19)产生ps个新个体；2)采用随机方式替换掉原种群中的ps个个体；3)对于每个萤火虫个体，计算其新位置对应的目标值f(yi)；4)记录全局最优值本文档来自技高网...

【技术保护点】
1.一种结合气象因素的超短期风电场功率预测方法，其特征在于：包括如下步骤：/n步骤一、对风电功率历史数据、NWP气象数据进行预处理，补齐缺失数据和修改异常数据；/n步骤二、生成风电功率预测模型，包括以下步骤：/n(1)利用pearson系数分析NWP数据和历史功率数据之间的相关性，对NWP数据和历史功率数据归一化处理，提取相关性强的气象数据作为模型的输入，用于训练预测模型；利用pearson系数对NWP系统提供的气象数据和风电场SCADA数据采集系统提供的历史功率数据做相关性分析，提取出的气象因素包括风速、温度、大气密度、湿度分别构成与历史功率数据具有相同时间步长的向量X和Y,然后利用pearson系数公式：/n

【技术特征摘要】
1.一种结合气象因素的超短期风电场功率预测方法，其特征在于：包括如下步骤：
步骤一、对风电功率历史数据、NWP气象数据进行预处理，补齐缺失数据和修改异常数据；
步骤二、生成风电功率预测模型，包括以下步骤：
(1)利用pearson系数分析NWP数据和历史功率数据之间的相关性，对NWP数据和历史功率数据归一化处理，提取相关性强的气象数据作为模型的输入，用于训练预测模型；利用pearson系数对NWP系统提供的气象数据和风电场SCADA数据采集系统提供的历史功率数据做相关性分析，提取出的气象因素包括风速、温度、大气密度、湿度分别构成与历史功率数据具有相同时间步长的向量X和Y,然后利用pearson系数公式：



对二者进行相关性分析，其中相关系数在(0.8-1.0)和(0.6-0.8)的气象因素为相关性强的气象因素，然后对提取出的气象因素与风电功率历史数据进行归一化处理，其归一化表达式为：



式中xmax和xmin分别为该变量的极大值和极小值；
经过预测模型得到的风电功率预测数据，再进行反归一化处理使其具有物理意义，反归一化的表达式为：
x＝0.5[x'(xmax-xmin)+(xmax-xmin)](3)
分析气象因素与风电功率的相关性能够更好的利用NWP数值天气预报的气象数据；
(2)基于门控循环单元神经网络对输入的气象数据处理得到初始预测功率；采用门控循环单元神经网络作为风电功率预测的预测模型，隐藏层信息传递是通过一个更新门u(t)和一个重置门r(t)来实现的，在t时刻GRU接收当前状态x(t)、上一时刻的隐藏状态h(t-1)，在信息传递过程中，更新门同时进行遗忘和记忆2个步骤，运算决定选择多少信息输入网络，遗忘多少过去的信息；重置门决定忘记过去的信息量；GRU网络隐藏层的输出h(t)最终由更新门和重置门的动态控制形成，GRU各个变量间的计算公式如下：
u(t)＝ψ(Wuh(t-1)+Uux(t))(4)
r(t)＝ψ(Wrh(t-1)+Urx(t))(5)






式中：Wu、Wr、Wc为连接隐藏层输出信号的参数矩阵；Uu、Ur、Uc为连接输入信号的参数矩阵；ψ(·)为激活函数，为sigmoid激活函数；为任何非线性函数是双曲正切函数tanh；为混合了结合了x(t)和r(t)中间记忆状态；⊙为逻辑运算符，表示矩阵中对应的元素相乘；
隐含层t时刻的h(t)将会作为t时刻的信息乘上相应的权值并传入GRU的输出层
y'i＝σ(ωih(t))+bi(8)
对yi进行反归一化处理
yi＝0.5[y′i(y'max-y'min)+(y'max+y'min)](9)
(3)基于混沌理论的萤火虫群体优化算法，对由门控循环单元神经网络得到的预测值和历史值构成的误差修正函数进行寻优，得到最优解对应的神经网络参数，生成预测模型；利用...

【专利技术属性】
技术研发人员：单锦宁，王洪哲，王荣茂，王琛淇，陈刚，王鑫，马欣慰，赵琰，宁兆秋，马艳娟，
申请(专利权)人：国网辽宁省电力有限公司阜新供电公司，国家电网有限公司，
类型：发明
国别省市：辽宁;21