CN102562469A - 基于校正算法的短期风力发电机输出功率预测方法 - Google Patents
基于校正算法的短期风力发电机输出功率预测方法 Download PDFInfo
- Publication number
- CN102562469A CN102562469A CN2011104427036A CN201110442703A CN102562469A CN 102562469 A CN102562469 A CN 102562469A CN 2011104427036 A CN2011104427036 A CN 2011104427036A CN 201110442703 A CN201110442703 A CN 201110442703A CN 102562469 A CN102562469 A CN 102562469A
- Authority
- CN
- China
- Prior art keywords
- wind speed
- weather forecast
- numerical weather
- historical data
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000012937 correction Methods 0.000 title claims abstract description 43
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 27
- 238000013528 artificial neural network Methods 0.000 claims abstract description 58
- 238000012549 training Methods 0.000 claims abstract description 22
- 238000000354 decomposition reaction Methods 0.000 claims description 44
- 230000006870 function Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 13
- 238000001514 detection method Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013398 bayesian method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Wind Motors (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了发电机输出功率管理与控制技术领域中的一种基于校正算法的短期风力发电机输出功率预测方法。包括:设定第一时间长度T1和第二时间长度T2;对第一时间长度T1的数值天气预报风速历史数据进行校正;以校正后的数值天气预报历史数据为输入值,以所述数值天气预报历史数据对应的风力发电机实测输出功率为输出值,训练BP神经网络;对第二时间长度T2的数值天气预报未来数据进行校正;将校正后的数值天气预报未来数据作为输入值输入到步骤3训练的BP神经网络中进行计算,得到的输出值即为风力发电机输出功率预测结果。本发明提高了短期风力发电机输出功率预测的准确性。
Description
技术领域
本发明属于发电机输出功率管理与控制技术领域,尤其涉及一种基于校正算法的短期风力发电机输出功率预测方法。
背景技术
风电固有的随机性和间歇性使得大规模并网风电场成为电力***安全运行的变数之一。因此,提高风电场功率的预测精度成为了电力公司和风电企业的迫切要求。运用数据校正算法,可以利用预测中的时间冗余信息,综合小波分析、模式识别和统计技术对预测数据进行处理,消除预测数据中包含的随机误差和显著误差,从而提高风电场功率预测结果的精度。
数据校正的工作主要包括显著误差检测和数据协调,显著误差检测技术的主要目的是对数据中的显著误差进行检测和识别或估算其真实值;而数据协调技术的主要目的是消除数据中随机误差,给出其协调值,并设法用已测数据对未测数据和未知参数予以估算。
显著误差检测分为基于测量残差和基于约束残差的检测方法。前者主要包括测量检验法(MT)、迭代测量检验法(IMT)、修正的迭代测量检验法(MIMT);后者主要有整体检验法(GT)、节点检验法(NT)、广义似然法(GLR)、主元分析法(PCA)等。另外神经网络、贝叶斯法、鲁棒估计法也被引入显著误差检测技术中,其中神经网络是一种非参数模型估计的方法,只用历史数据进行模型训练,不需要建立精确的数学模型,因此避免了过程建模不准确带来的估计偏差。此外,神经网络对测量数据中的随机误差的分布不作特殊要求,处理非线性问题的能力较强,计算量小,适于在线应用。
数据协调由Kuehn和Davidson提出(Kuehn D R,Davidson H.Computercontrol II.Mathematics of control[J].Chemical Engineering Progress,1961,57(6):44-47),其发展经历了从稳态数据协调、动态数据协调到非线性数据协调的过程。Darouach将基于Kalman滤波的算法用于广义线性动态***(Darouach M,Zasadzinski M.Data reconciliation in generalized linear dynamicsystems[J].AICHE J.,1991,37(2):193-210)。Liebman等提出用非线性规划法(NLP)解决非线性约束和不等式约束的动态数据协调问题(Liebman M J,EdgarT F,Lasdon L S.Efficient data reconciliation estimation for dynamic processesusing nonlinear programming technique[J].Computers Chem Engng J,1992,16(10-11):963-986)。Karjala等用Elman神经元网络来求解非线性的数据协调问题(Karjala T W,Himmelblau D M,Miikkulainen R.Data rectification usingrecurrent(Elman)neural networks[C].Intemational Joint Conference on NeuralNetwork,1992,2:901-906)。
在短期风电功率预测中,实测数据和数值天气预报(Numerical WeatherPrediction,NWP)是最重要的输入参数之一,也是最主要的误差源之一。对于风电场实测数据,文献(Zhiling Yang,Yongqian Liu,Chengrong Li.Interpolationof missing wind data based on ANFIS[C].Renewable Energy,2011,36(30):993-998)提出了采用ANFIS补充缺失数据和无效数据,减小了实测数据的误差。
理论上风力发电机的输出功率可以下面的公式表达:
其中,P为输出功率,单位为千瓦(KW);CP为风机的功率系数;ρ为空气密度,单位为千克/立方米(kg/m3);A为风轮的扫掠面积,单位为平方米(m2);v为风速,单位为米/秒(m/s)。同时,实际的输出功率还要受到湍流以及风机偏航装置对风向改变的滞后性的影响。
在风电功率预测中,最常用的方法是建立输出功率和风速之间的映射模型,而两者之间的分散性越小,模型越准确。如图1所示,风机的实测输出功率和实测风速之间的分散性较小,根据实测风速很容易建立输出功率和实测风速之间的映射模型。
在短期风电功率预测中,使用未来24h~72h(h为小时)的数值天气预报风速作为预测模型的输入,来预测未来24h~72h的风电输出功率。其思路是:由气象部门提供的风电场范围内24h~72h的数值天气预报服务(风速、风向等),通过插值得到每台风机轮毂高度的24h~72h风速、风向预报值,然后根据每台风机的风速-功率曲线模型,或者人工神经网络、支持向量机等人工智能算法,映射为每台风机输出功率的预测值,最后得到整个风电场的输出功率的预测值。
其中,数值天气预报是根据大气实际情况,在一定的初值和边值条件下,通过大型计算机作数值计算,求解描述天气演变过程的流体力学和热力学方程组,预测未来一定时段的大气运动状态和天气现象的方法。图2为数值天气预报风速与实测的风机输出功率的散点图,从该图中可以看出,由于数值天气预报风速和测风塔实测风速之间误差较大,很难从数据中准确地提炼出数值天气预报风速和风电场总输出功率之间的关系曲线。因此,需要对数值天气预报风速进行进一步的校正,减小预测误差。
数值天气预报产品来源于数值模式,是大型计算机通过数值计算,求解描写天气演变过程的流体力学和热力学方程组,作出未来一段时间内的定量和客观的预报。数值模式能给出未来大气的气压、温度、风、湿度等各种状态。但是,大气运动是一种混沌现象,在数值计算的过程中,许多难以预料的细微因素往往会影响最终的预报结果。图3显示了某风电场的某风机轮毂高度的1h平均风速的数值天气预报数据和实测数据时间序列。从图中可以看出两者之间存在着较大的误差,这也是风电功率预报的主要误差来源。
发明内容
本发明的目的在于,针对现有的短期风力发电机输出功率预测方法预测结果不准确的问题,提出一种基于校正算法的短期风力发电机输出功率预测方法。
为实现上述目的,本发明提供的技术方案是,一种基于校正算法的短期风力发电机输出功率预测方法,其特征是所述方法包括:
步骤1:设定第一时间长度T1和第二时间长度T2;
步骤2:对第一时间长度T1的数值天气预报风速历史数据进行校正;
步骤3:以校正后的数值天气预报历史数据为输入值,以所述数值天气预报历史数据对应的风力发电机实测输出功率为输出值,训练BP神经网络;
步骤4:对第二时间长度T2的数值天气预报未来数据进行校正;
步骤5:将校正后的数值天气预报未来数据作为输入值输入到步骤3训练的BP神经网络中进行计算,得到的输出值即为风力发电机输出功率预测结果。
所述对第一时间长度T1的数值天气预报风速历史数据进行校正具体包括:
步骤102:设定第三时间长度T3并令i=1;
步骤103:以第i组数值天气预报风速历史数据之前的第三时间长度T3的数值天气预报风速历史数据作为既有数值天气预报风速历史数据,以第i组数值天气预报风速历史数据之前的第三时间长度T3所对应的实测风速历史数据作为既有实测风速历史数据;
步骤104:对既有数值天气预报风速历史数据进行小波包分解,得到m个既有数值天气预报风速历史数据子序列;对既有实测风速历史数据进行小波包分解,得到m个既有实测风速历史数据子序列;
步骤105:以既有数值天气预报风速历史数据为输入层数据,以既有实测风速历史数据为输出层数据,以第j个既有数值天气预报风速历史数据子序列和第j个既有实测风速历史数据子序列为训练样本集训练BP神经网络,得到第j个BP神经网络及其对应的隐层节点权值和阈值;其中,j=1,2,…,m;
步骤106:将第i组数值天气预报风速历史数据记作当前数值天气预报风速历史数据,对当前数值天气预报风速历史数据进行小波包分解,得到m个当前数值天气预报风速历史数据子序列;
步骤107:将第j个当前数值天气预报风速历史数据子序列输入到第j个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值,将该输出值作为修正后的第j个当前数值天气预报风速历史数据子序列;其中,j=1,2,…,m;
步骤108:对修正后的m个当前数值天气预报风速历史数据子序列进行小波包重构,得到校正后的当前数值天气预报风速历史数据,即校正后的第i组数值天气预报风速历史数据;
步骤109:判断i≤n是否成立,如果成立,则令i=i+1,返回步骤103;否则,结束。
所述校正第二时间长度T2的数值天气预报未来数据具体包括:
步骤201:以24小时为长度单位,将第二时间长度T2的数值天气预报风速未来数据划分为n组;其中,
步骤202:设定第三时间长度T3并令i=1;
步骤203:以第i组数值天气预报风速未来数据之前的第三时间长度T3的数值天气预报风速历史数据作为既有数值天气预报风速历史数据,以第i组数值天气预报风速未来数据之前的第三时间长度T3所对应的实测风速历史数据作为既有实测风速历史数据;
步骤204:对既有数值天气预报风速历史数据进行小波包分解,得到m个既有数值天气预报风速历史数据子序列;对既有实测风速历史数据进行小波包分解,得到m个既有实测风速历史数据子序列;
步骤205:以既有数值天气预报风速历史数据为输入层数据,以既有实测风速历史数据为输出层数据,以第j个既有数值天气预报风速历史数据子序列和第j个既有实测风速历史数据子序列为训练样本集训练BP神经网络,得到第j个BP神经网络及其对应的隐层节点权值和阈值;其中,j=1,2,…,m;
步骤206:将第i组数值天气预报风速未来数据记作当前数值天气预报风速数据,对当前数值天气预报风速数据进行小波包分解,得到m个当前数值天气预报风速数据子序列;
步骤207:将第j个当前数值天气预报风速数据子序列输入到第j个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值,将该输出值作为修正后的第j个当前数值天气预报风速数据子序列;其中,j=1,2,…,m;
步骤208:对m个当前数值天气预报风速数据子序列进行小波包重构,得到校正后的当前数值天气预报风速数据,即校正后的第i组数值天气预报风速未来数据;
步骤209:判断i≤n是否成立,如果成立,则令i=i+1,返回步骤203;否则,结束。
所述对既有数值天气预报风速历史数据/既有实测风速历史数据进行小波包分解,具体采用Daubechies小波对既有数值天气预报风速历史数据/既有实测风速历史数据进行三层小波包分解。
所述小波包的递归表达式为
其中,hk为正交尺度函数对应的正交低通实系数滤波器,gk为正交小波函数对应的高通滤波器,gn=(-1)nhn-1,n∈Z,且满足尺度方程和小波方程
u0(t)是正交尺度函数,u1(t)是正交小波函数,t是自变量,n∈Z,k∈Z。
本发明提高了短期风力发电机输出功率预测的准确性。
附图说明
图1是2006年4月实测风速与某风机实测输出功率的散点图;
图2是2006年4月数值天气预报风速与某风机实测输出功率的散点图;
图3是2006年4月某风机轮毂高度实测风速与数值天气预报风速对比图;
图4是基于校正算法的短期风力发电机输出功率预测方法流程图;
图5是2006年4月整月的数值天气预报风速历史数据归一化数据图;
图6是小波包分解示意图;
图7是三层BP神经网络结构图;
图8是对既有数值天气预报风速历史数据和既有实测风速历史数据进行小波包分解后示意图;其中,(a)是小波包分解后的D(3,0)示意图;,(b)是小波包分解后的D(3,1)示意图;(c)是小波包分解后的D(3,2)示意图;(d)是小波包分解后的D(3,3)示意图;(e)是小波包分解后的D(3,4)示意图;(f)是小波包分解后的D(3,5)示意图;(g)是小波包分解后的D(3,6)示意图;(h)是小波包分解后的D(3,7)示意图;
图9是当前数值天气预报风速数据子序列校正前后的对比图;其中,(a)是校正前后的D(3,0)对比图;,(b)是校正前后的D(3,1)对比图;(c)是校正前后的D(3,2)对比图;(d)是校正前后的D(3,3)对比图;(e)是校正前后的D(3,4)对比图;(f)是校正前后的D(3,5)对比图;(g)是校正前后的D(3,6)对比图;(h)是校正前后的D(3,7)对比图;
图10是校正前后的数值天气预报风速未来数据与同期实测风速数据的对比图;
图11是采用校正算法和未采用校正算法的短期风力发电机输出功率预测值与风力发电机输出功率实测值的对比图。
图12是校正前后的预测输出功率和实测输出功率对比图;
图13是使用常用的预测方法和本发明的预测方法得到的输出功率与实测输出功率的对比图;
图14是根据2006年3月的数值气象预测风速、实测风速和风机实测功率数值图。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
图4是基于校正算法的短期风力发电机输出功率预测方法流程图。图4中,本发明提供的方法包括:
步骤1:设定第一时间长度T1=720小时和第二时间长度T2=24小时。
步骤2:对第一时间长度T1=720小时的数值天气预报风速历史数据进行校正。
本实施例以2006年4月整月的数值天气预报风速历史数据作为第一时间长度的数值天气预报风速历史数据并对其进行校正,如图5所示,具体过程是:
步骤102:设定第三时间长度T3=720小时并令i=1。
步骤103:以第i组数值天气预报风速历史数据之前的第三时间长度T3=720小时的数值天气预报风速历史数据作为既有数值天气预报风速历史数据。
对于2006年4月1日的数值天气预报风速历史数据(即第1组数值天气预报风速历史数据),其之前的720小时的数值天气预报风速历史数据为2006年3月2日-2006年3月31日的数值天气预报风速历史数据。将2006年3月2日-2006年3月31日的数值天气预报风速历史数据作为既有数值天气预报风速历史数据。同时,以2006年3月2日-2006年3月31日的实测风速历史数据作为既有实测风速历史数据。
对于2006年4月2日的数值天气预报风速历史数据(即第2组数值天气预报风速历史数据),其之前的720小时的数值天气预报风速历史数据为2006年3月3日-2006年4月1日的数值天气预报风速历史数据。将2006年3月3日-2006年4月1日的数值天气预报风速历史数据作为既有数值天气预报风速历史数据。同时,以2006年3月3日-2006年4月1日的实测风速历史数据作为既有实测风速历史数据。
其他组别(2006年4月其他日期)数值天气预报风速历史数据之前的第三时间长度T3=720小时的数值天气预报风速历史数据和实测风速历史数据以此类推。
步骤104:对既有数值天气预报风速历史数据进行小波包分解,得到m个既有数值天气预报风速历史数据子序列;对既有实测风速历史数据进行小波包分解,得到m个既有实测风速历史数据子序列。
以第i=30组为例,即2006年4月30日的数值天气预报风速历史数据为例进行说明。其之前的720小时的数值天气预报风速历史数据为2006年3月31日-2006年4月29日的数值天气预报风速历史数据;其之前的720小时所对应的实测风速历史数据为2006年3月31日-2006年4月29日的实测风速历史数据,将其作为既有实测风速历史数据。分别对既有数值天气预报风速历史数据和既有实测风速历史数据进行小波包分解。
小波包分解具体采用Daubechies小波进行三层小波包分解。Daubechies小波系是由法国学者Daubechies提出的一系列二进制小波的总称,在Matlab中记为dbN,N为小波的序号,N值取2,3,…,10。本步骤实际是通过小波变换将信号分解成多个具有不同时间分辨率和频率分辨率的信号,从而揭示信号在不同尺度上的时域行为特征。但是由于正交小波变换只对信号的近似部分(低频部分)进行进一步分解,而对细节部分(高频部分)不再分解,因此它不能有效地处理存在大量细节特征的风速序列数据。而小波包变换可以对信号的细节部分进行更精细的分解,更适合被用于分析风速序列数据。
小波包的递归表达式为
其中,hk为正交尺度函数对应的正交低通实系数滤波器,gk为正交小波函数对应的高通滤波器,gn=(-1)nhn-1,n∈Z,且满足尺度方程和小波方程
u0(t)是正交尺度函数,u1(t)是正交小波函数,t是自变量,n∈Z,k∈Z。
小波包分解过程可以用小波分解树来表示,如图6所示。图中的D(0,0)是一个原始的一维离散信号,如一个既有数值天气预报风速历史数据;D(t,k)是分解得到的子序列,t是第t层分解,k是第k+1个子序列的编号,如D(0,1)是经过第一层小波包分解后得到的第1个子序列。
既有数值天气预报风速历史数据和既有实测风速历史数据经过上述三层小波包分解后,分别得到8个既有数值天气预报风速历史数据子序列和8个既有实测风速历史数据子序列,如图8所示。
步骤105:以既有数值天气预报风速历史数据为输入层数据,以既有实测风速历史数据为输出层数据,以第j个既有数值天气预报风速历史数据子序列和第j个既有实测风速历史数据子序列为训练样本集训练BP神经网络,得到第j个BP神经网络及其对应的隐层节点权值和阈值;其中,j=1,2,…,8。
BP神经网络的结构如图7所示。8个既有数值天气预报风速历史数据子序列和8个既有实测风速历史数据子序列一共可以训练得到8个BP神经网络及其与之对应的隐层节点权值和阈值。
在训练BP神经网络时,可以选用自适应的学习率递归调整方法来确定权值的变化率。其公式为:
其中,η(t)为学***方和,在本发明中为BP神经网络输出结果与既有既有实测风速历史数据的误差平方和。选用该方法可以避免学习率过大导致算法的不收敛或者学习率过小将导致算法收敛很慢以致于训练时间过长的问题。
步骤106:将第i组数值天气预报风速历史数据记作当前数值天气预报风速历史数据,对当前数值天气预报风速历史数据进行小波包分解,得到m个当前数值天气预报风速历史数据子序列。
由于步骤104以第i=30组为例进行说明,因此本步骤第i=30组数值天气预报风速历史数据即2006年4月30日的数值天气预报风速历史数据,被记作当前数值天气预报风速历史数据。对当前数值天气预报风速历史数据进行小波包分解的过程同步骤105,分解后得到8个当前数值天气预报风速历史数据子序列。
步骤107:将第j个当前数值天气预报风速历史数据子序列输入到第j个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值,将该输出值作为修正后的第j个当前数值天气预报风速历史数据子序列;其中,j=1,2,…,8。在训练BP神经网络时,可以如步骤105那样选用自适应的学习率递归调整方法来确定权值的变化率。
8个当前数值天气预报风速历史数据子序列分别输入到8个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值。计算结果有8个输出值,这8个输出值就作为修正后的当前数值天气预报风速历史数据子序列。
步骤108:对8个当前数值天气预报风速历史数据子序列进行小波包重构,得到校正后的当前数值天气预报风速历史数据,即校正后的第i组数值天气预报风速历史数据。
小波包重构实际上是小波包分解的逆过程,目的在于将分解的子序列重组。重组后的数据即为修正后的数据。
步骤109:判断i≤n是否成立,如果成立,则令i=i+1,返回步骤103;否则,结束。
到此,已经对2006年4月每天的数值天气预报风速历史数据进行了修正,从而得到了修正后的2006年4月的数值天气预报风速历史数据。
步骤3:以校正后的数值天气预报历史数据为输入值,以所述数值天气预报历史数据对应的风力发电机实测输出功率为输出值,训练BP神经网络。
本步骤以修正后的2006年4月的数值天气预报风速历史数据为输入值,以2006年4月的风力发电机实测输出功率为输出值,训练BP神经网络。
步骤4:对第二时间长度T2=24小时的数值天气预报未来数据进行校正。
以2006年5月1日的数值天气预报数据作为数值天气预报未来数据。本步骤的校正过程与步骤2类似,具体地:
步骤201:以24小时为长度单位,将第二时间长度T2的数值天气预报风速未来数据划分为n组;其中,
步骤202:设定第三时间长度T3=720小时并令i=1。
步骤203:由于本实施例只划分出1组数值天气预报风速未来数据,下面就以这组数值天气预报风速未来数据为例进行说明。
以2006年5月1日的数值天气预报风速未来数据之前的第三时间长度720小时的数值天气预报风速历史数据作为既有数值天气预报风速历史数据,即以2006年4月1日-4月30日的数值天气预报风速历史数据作为既有数值天气预报风速历史数据。
以第1组数值天气预报风速未来数据之前的第三时间长度720小时所对应的实测风速历史数据作为既有实测风速历史数据,即以2006年4月1日-4月30日的的实测风速历史数据作为既有实测风速历史数据。
步骤204:对既有数值天气预报风速历史数据进行三层Daubechies小波包分解,得到8个既有数值天气预报风速历史数据子序列;对既有实测风速历史数据进行三层Daubechies小波包分解,得到8个既有实测风速历史数据子序列。
步骤205:以既有数值天气预报风速历史数据为输入层数据,以既有实测风速历史数据为输出层数据,以第j个既有数值天气预报风速历史数据子序列和第j个既有实测风速历史数据子序列为训练样本集训练BP神经网络,得到第j个BP神经网络及其对应的隐层节点权值和阈值;其中,j=1,2,…,8。在训练BP神经网络时,可以如步骤105那样选用自适应的学习率递归调整方法来确定权值的变化率。
步骤206:将第1组数值天气预报风速未来数据记作当前数值天气预报风速数据,对当前数值天气预报风速数据进行三层Daubechies小波包分解,得到8个当前数值天气预报风速数据子序列。
步骤207:将第j个当前数值天气预报风速数据子序列输入到第j个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值,将该输出值作为校正后的第j个当前数值天气预报风速数据子序列;其中,j=1,2,…,8。在训练BP神经网络时,可以如步骤105那样选用自适应的学习率递归调整方法来确定权值的变化率。得到的校正后的子序列与校正前的子序列的对比图如图9所示。
步骤208:对8个当前数值天气预报风速数据子序列进行小波包重构,得到校正后的当前数值天气预报风速数据,即校正后的第1组数值天气预报风速未来数据。
步骤209:判断i≤n是否成立,如果成立,则令i=i+1,返回步骤203;否则,结束。
到此,得到校正后的2006年5月1日数值天气预报风速未来数据。图10是校正前后的数值天气预报风速未来数据与同期实测风速数据的对比图,将校正前后的数据与同期的实测风速进行比较,其平均绝对误差下降了32.24%,均方根误差下降了38.75%。
步骤5:将校正后的数值天气预报未来数据作为输入值输入到步骤3训练的BP神经网络中进行计算,得到的输出值即为风力发电机输出功率预测结果。
由于选定的数值天气预报未来数据为2006年5月1日数值天气预报风速数据,因此本步骤得到的风力发电机输出功率预测结果为2006年5月1日的风力发电机输出功率预测结果。
图11是采用校正算法和未采用校正算法的短期风力发电机输出功率预测值与风力发电机输出功率实测值的对比图。通过比较,其平均绝对误差下降了37.85%,均方根误差下降了24.52%。
当然,本发明还可以进一步改进,与数值天气预报风速的校正相似,也可以使用实测输出功率对预测输出功率进行进一步校正。图12显示了校正前后的4月29日9时至4月30日8时的预测输出功率以及实测输出功率。通过比较,其平均绝对误差下降了40.56%,均方根误差下降了31.65%。
图13显示了直接采用数值天气预报风速数据作为输入量的BP神经网络预测模型经过训练后,得到的4月29日9时至4月30日8时的预测输出功率,以及采用了数据校正算法后的预测输出功率以及同期的实测功率。采用数据校正算法后,预测输出功率的平均绝对误差下降了63.06%,均方根误差下降了48.04%。结果显示,数据校正算法取得了良好的效果。
图14显示了根据2006年3月的数值气象预测风速、实测风速和风机实测功率,采用数据校正算法前后的3月30日9时至3月31日8时的预测输出功率和实测功率。采用数据校正算法后,预测输出功率的平均绝对误差增加了19.37%,均方根误差增加了38.35%。结果显示,采用数据校正算法反而降低了预测准确性。这是因为2006年3月的数值天气预报风速与实测风速的数据相关性为0.831353,比较准确的数值天气预报风速得到了比较准确的预测输出功率。而2006年4月的数值天气预报风速与实测风速的数据相关性为0.578375,较大的数值天气预报误差使得预测输出功率的误差也较大。这表明在数值天气预报误差较大的情况下,数值校正算法能够提高预测准确性,但是当数值天气预报误差较小的情况下,数值校正算法反而会降低预测准确性。
本发明提出了一种数据校正算法,首次将数据校正算法应用于数值天气预报风速和风电功率短期预测结果的校正。通过对中国北方某风场的2006年3月和4月的数据进行计算,结果表明:随数值天气预报风速的准确性的不同,数据校正算法的实际效果也不同。在数值天气预报误差较大的情况下,数据校正算法能够取得比较理想的效果,较好地提高了预测的准确性。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (5)
1.一种基于校正算法的短期风力发电机输出功率预测方法,其特征是所述方法包括:
步骤1:设定第一时间长度T1和第二时间长度T2;
步骤2:对第一时间长度T1的数值天气预报风速历史数据进行校正;
步骤3:以校正后的数值天气预报历史数据为输入值,以所述数值天气预报历史数据对应的风力发电机实测输出功率为输出值,训练BP神经网络;
步骤4:对第二时间长度T2的数值天气预报未来数据进行校正;
步骤5:将校正后的数值天气预报未来数据作为输入值输入到步骤3训练的BP神经网络中进行计算,得到的输出值即为风力发电机输出功率预测结果。
2.根据权利要求1所述的方法,其特征是所述对第一时间长度T1的数值天气预报风速历史数据进行校正具体包括:
步骤102:设定第三时间长度T3并令i=1;
步骤103:以第i组数值天气预报风速历史数据之前的第三时间长度T3的数值天气预报风速历史数据作为既有数值天气预报风速历史数据,以第i组数值天气预报风速历史数据之前的第三时间长度T3所对应的实测风速历史数据作为既有实测风速历史数据;
步骤104:对既有数值天气预报风速历史数据进行小波包分解,得到m个既有数值天气预报风速历史数据子序列;对既有实测风速历史数据进行小波包分解,得到m个既有实测风速历史数据子序列;
步骤105:以既有数值天气预报风速历史数据为输入层数据,以既有实测风速历史数据为输出层数据,以第j个既有数值天气预报风速历史数据子序列和第j个既有实测风速历史数据子序列为训练样本集训练BP神经网络,得到第j个BP神经网络及其对应的隐层节点权值和阈值;其中,j=1,2,…,m;
步骤106:将第i组数值天气预报风速历史数据记作当前数值天气预报风速历史数据,对当前数值天气预报风速历史数据进行小波包分解,得到m个当前数值天气预报风速历史数据子序列;
步骤107:将第j个当前数值天气预报风速历史数据子序列输入到第j个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值,将该输出值作为修正后的第j个当前数值天气预报风速历史数据子序列;其中,j=1,2,…,m;
步骤108:对修正后的m个当前数值天气预报风速历史数据子序列进行小波包重构,得到校正后的当前数值天气预报风速历史数据,即校正后的第i组数值天气预报风速历史数据;
步骤109:判断i≤n是否成立,如果成立,则令i=i+1,返回步骤103;否则,结束。
3.根据权利要求1所述的方法,其特征是所述校正第二时间长度T2的数值天气预报未来数据具体包括:
步骤202:设定第三时间长度T3并令i=1;
步骤203:以第i组数值天气预报风速未来数据之前的第三时间长度T3的数值天气预报风速历史数据作为既有数值天气预报风速历史数据,以第i组数值天气预报风速未来数据之前的第三时间长度T3所对应的实测风速历史数据作为既有实测风速历史数据;
步骤204:对既有数值天气预报风速历史数据进行小波包分解,得到m个既有数值天气预报风速历史数据子序列;对既有实测风速历史数据进行小波包分解,得到m个既有实测风速历史数据子序列;
步骤205:以既有数值天气预报风速历史数据为输入层数据,以既有实测风速历史数据为输出层数据,以第j个既有数值天气预报风速历史数据子序列和第j个既有实测风速历史数据子序列为训练样本集训练BP神经网络,得到第j个BP神经网络及其对应的隐层节点权值和阈值;其中,j=1,2,…,m;
步骤206:将第i组数值天气预报风速未来数据记作当前数值天气预报风速数据,对当前数值天气预报风速数据进行小波包分解,得到m个当前数值天气预报风速数据子序列;
步骤207:将第j个当前数值天气预报风速数据子序列输入到第j个BP神经网络并利用该BP神经网络对应的隐层节点权值和阈值,计算得到该BP神经网络的输出值,将该输出值作为修正后的第j个当前数值天气预报风速数据子序列;其中,j=1,2,…,m;
步骤208:对m个当前数值天气预报风速数据子序列进行小波包重构,得到校正后的当前数值天气预报风速数据,即校正后的第i组数值天气预报风速未来数据;
步骤209:判断i≤n是否成立,如果成立,则令i=i+1,返回步骤203;否则,结束。
4.根据权利要求2或3所述的方法,其特征是所述对既有数值天气预报风速历史数据/既有实测风速历史数据进行小波包分解,具体采用Daubechies小波对既有数值天气预报风速历史数据/既有实测风速历史数据进行三层小波包分解。
5.根据权利要求4所述的方法,其特征是所述小波包的递归表达式为
其中,hk为正交尺度函数对应的正交低通实系数滤波器,gk为正交小波函数对应的高通滤波器,gn=(-1)nhn-1,n∈Z,且满足尺度方程和小波方程
u0(t)是正交尺度函数,u1(t)是正交小波函数,t是自变量,n∈Z,k∈Z。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110442703.6A CN102562469B (zh) | 2011-12-27 | 2011-12-27 | 基于校正算法的短期风力发电机输出功率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110442703.6A CN102562469B (zh) | 2011-12-27 | 2011-12-27 | 基于校正算法的短期风力发电机输出功率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102562469A true CN102562469A (zh) | 2012-07-11 |
CN102562469B CN102562469B (zh) | 2014-01-22 |
Family
ID=46408804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110442703.6A Active CN102562469B (zh) | 2011-12-27 | 2011-12-27 | 基于校正算法的短期风力发电机输出功率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102562469B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930356A (zh) * | 2012-11-16 | 2013-02-13 | 广东电网公司电力调度控制中心 | 基于气象因素敏感度的短期负荷预测方法 |
CN103001249A (zh) * | 2012-11-15 | 2013-03-27 | 中国电力科学研究院 | 基于bp神经网络的风电场短期功率预测方法 |
CN103473618A (zh) * | 2013-09-18 | 2013-12-25 | 上海电机学院 | 一种风电场短期风功率预测***及方法 |
CN103485977A (zh) * | 2013-09-06 | 2014-01-01 | 河海大学 | 风力发电***实时功率预测的修正方法 |
CN105591407A (zh) * | 2016-03-03 | 2016-05-18 | 国家电网公司 | 可再生能源电站有功出力预测误差间相关性的研究方法 |
CN107609697A (zh) * | 2017-09-06 | 2018-01-19 | 南京邮电大学 | 一种风力发电功率预测方法 |
CN107942410A (zh) * | 2017-09-19 | 2018-04-20 | 武汉船用机械有限责任公司 | 一种极地服务温度的预测方法和装置 |
CN109154281A (zh) * | 2016-05-23 | 2019-01-04 | 通用电气公司 | 用于预测风场的功率输出的***和方法 |
CN110188967A (zh) * | 2019-06-21 | 2019-08-30 | 合肥工业大学 | 基于混沌人群算法和贝叶斯网络的电力负荷概率性预测方法 |
CN110264006A (zh) * | 2019-06-21 | 2019-09-20 | 合肥工业大学 | 基于混沌萤火虫算法和贝叶斯网络的风电功率概率性预测方法 |
CN110365041A (zh) * | 2019-06-05 | 2019-10-22 | 华南理工大学 | 基于gan场景模拟的风电多场景鲁棒调度方法 |
CN110809672A (zh) * | 2017-07-07 | 2020-02-18 | 西门子歌美飒可再生能源公司 | 确定风速值 |
CN111336062A (zh) * | 2020-03-05 | 2020-06-26 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 一种基于测量风速的风力发电机组最大风能捕获方法 |
CN111476679A (zh) * | 2020-04-14 | 2020-07-31 | 四川北控清洁能源工程有限公司 | 根据大气稳定度修正风电机组功率曲线的方法 |
CN113153657A (zh) * | 2021-05-28 | 2021-07-23 | 上海扩博智能技术有限公司 | 风机发电率损失预测方法、***、设备及介质 |
CN113969873A (zh) * | 2021-10-18 | 2022-01-25 | 国网湖南省电力有限公司 | 一种风电机组出力超短期预测方法、***及存储介质 |
CN114169585A (zh) * | 2021-11-19 | 2022-03-11 | 西安理工大学 | 基于马尔可夫链和组合模型的风力发电功率预测方法 |
CN116757101A (zh) * | 2023-08-21 | 2023-09-15 | 湖南科技大学 | 一种基于机理模型和神经网络的机舱风速修正方法和*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004019583A (ja) * | 2002-06-18 | 2004-01-22 | Tohoku Electric Power Co Inc | 風力発電における発電出力予測方法、発電出力予測装置及び発電出力予測システム |
CN102063641A (zh) * | 2010-10-14 | 2011-05-18 | 北京大学 | 高速铁路沿线风速预测方法 |
WO2011110404A1 (de) * | 2010-03-12 | 2011-09-15 | Siemens Aktiengesellschaft | Verfahren zur rechnergestützten steuerung und/oder regelung eines technischen systems |
-
2011
- 2011-12-27 CN CN201110442703.6A patent/CN102562469B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004019583A (ja) * | 2002-06-18 | 2004-01-22 | Tohoku Electric Power Co Inc | 風力発電における発電出力予測方法、発電出力予測装置及び発電出力予測システム |
WO2011110404A1 (de) * | 2010-03-12 | 2011-09-15 | Siemens Aktiengesellschaft | Verfahren zur rechnergestützten steuerung und/oder regelung eines technischen systems |
CN102063641A (zh) * | 2010-10-14 | 2011-05-18 | 北京大学 | 高速铁路沿线风速预测方法 |
Non-Patent Citations (3)
Title |
---|
刘永前,朴金姬,韩爽: "风电场输出功率预测中两种神经网络算法的研究", 《现代电力》, vol. 28, no. 2, 10 April 2011 (2011-04-10) * |
刘永前,韩爽,胡永生: "风电场处理短期预报研究综述", 《现代电力》, vol. 24, no. 5, 10 October 2007 (2007-10-10) * |
赵攀,戴义平,夏俊荣,盛迎新: "卡尔曼滤波修正的风电场短期功率预测模型", 《西安交通大学学报》, vol. 45, no. 5, 17 January 2011 (2011-01-17) * |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103001249A (zh) * | 2012-11-15 | 2013-03-27 | 中国电力科学研究院 | 基于bp神经网络的风电场短期功率预测方法 |
CN103001249B (zh) * | 2012-11-15 | 2016-01-20 | 中国电力科学研究院 | 基于bp神经网络的风电场短期功率预测方法 |
CN102930356A (zh) * | 2012-11-16 | 2013-02-13 | 广东电网公司电力调度控制中心 | 基于气象因素敏感度的短期负荷预测方法 |
CN102930356B (zh) * | 2012-11-16 | 2016-06-22 | 广东电网公司电力调度控制中心 | 基于气象因素敏感度的短期负荷预测方法 |
CN103485977A (zh) * | 2013-09-06 | 2014-01-01 | 河海大学 | 风力发电***实时功率预测的修正方法 |
CN103485977B (zh) * | 2013-09-06 | 2015-08-05 | 河海大学 | 风力发电***实时功率预测的修正方法 |
CN103473618A (zh) * | 2013-09-18 | 2013-12-25 | 上海电机学院 | 一种风电场短期风功率预测***及方法 |
CN105591407A (zh) * | 2016-03-03 | 2016-05-18 | 国家电网公司 | 可再生能源电站有功出力预测误差间相关性的研究方法 |
US11242842B2 (en) | 2016-05-23 | 2022-02-08 | General Electric Company | System and method for forecasting power output of a wind farm |
CN109154281B (zh) * | 2016-05-23 | 2021-09-24 | 通用电气公司 | 用于预测风场的功率输出的***和方法 |
CN109154281A (zh) * | 2016-05-23 | 2019-01-04 | 通用电气公司 | 用于预测风场的功率输出的***和方法 |
CN110809672A (zh) * | 2017-07-07 | 2020-02-18 | 西门子歌美飒可再生能源公司 | 确定风速值 |
CN110809672B (zh) * | 2017-07-07 | 2021-07-23 | 西门子歌美飒可再生能源公司 | 确定风速值 |
CN107609697A (zh) * | 2017-09-06 | 2018-01-19 | 南京邮电大学 | 一种风力发电功率预测方法 |
CN107942410B (zh) * | 2017-09-19 | 2020-04-21 | 武汉船用机械有限责任公司 | 一种极地服务温度的预测方法和装置 |
CN107942410A (zh) * | 2017-09-19 | 2018-04-20 | 武汉船用机械有限责任公司 | 一种极地服务温度的预测方法和装置 |
CN110365041A (zh) * | 2019-06-05 | 2019-10-22 | 华南理工大学 | 基于gan场景模拟的风电多场景鲁棒调度方法 |
CN110365041B (zh) * | 2019-06-05 | 2021-05-14 | 华南理工大学 | 基于gan场景模拟的风电多场景鲁棒调度方法 |
CN110188967A (zh) * | 2019-06-21 | 2019-08-30 | 合肥工业大学 | 基于混沌人群算法和贝叶斯网络的电力负荷概率性预测方法 |
CN110264006A (zh) * | 2019-06-21 | 2019-09-20 | 合肥工业大学 | 基于混沌萤火虫算法和贝叶斯网络的风电功率概率性预测方法 |
CN110264006B (zh) * | 2019-06-21 | 2021-01-15 | 合肥工业大学 | 基于混沌萤火虫算法和贝叶斯网络的风电功率概率性预测方法 |
CN110188967B (zh) * | 2019-06-21 | 2021-01-15 | 合肥工业大学 | 基于混沌人群算法和贝叶斯网络的电力负荷概率性预测方法 |
CN111336062A (zh) * | 2020-03-05 | 2020-06-26 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 一种基于测量风速的风力发电机组最大风能捕获方法 |
CN111336062B (zh) * | 2020-03-05 | 2021-11-09 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 一种基于测量风速的风力发电机组最大风能捕获方法 |
CN111476679A (zh) * | 2020-04-14 | 2020-07-31 | 四川北控清洁能源工程有限公司 | 根据大气稳定度修正风电机组功率曲线的方法 |
CN113153657A (zh) * | 2021-05-28 | 2021-07-23 | 上海扩博智能技术有限公司 | 风机发电率损失预测方法、***、设备及介质 |
CN113969873A (zh) * | 2021-10-18 | 2022-01-25 | 国网湖南省电力有限公司 | 一种风电机组出力超短期预测方法、***及存储介质 |
CN113969873B (zh) * | 2021-10-18 | 2023-04-07 | 国网湖南省电力有限公司 | 一种风电机组出力超短期预测方法、***及存储介质 |
CN114169585A (zh) * | 2021-11-19 | 2022-03-11 | 西安理工大学 | 基于马尔可夫链和组合模型的风力发电功率预测方法 |
CN116757101A (zh) * | 2023-08-21 | 2023-09-15 | 湖南科技大学 | 一种基于机理模型和神经网络的机舱风速修正方法和*** |
CN116757101B (zh) * | 2023-08-21 | 2023-11-07 | 湖南科技大学 | 一种基于机理模型和神经网络的机舱风速修正方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN102562469B (zh) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102562469B (zh) | 基于校正算法的短期风力发电机输出功率预测方法 | |
Zhao et al. | Research and application based on the swarm intelligence algorithm and artificial intelligence for wind farm decision system | |
Tascikaraoglu et al. | A review of combined approaches for prediction of short-term wind speed and power | |
CN103996084B (zh) | 基于纵向时刻Markov链模型的风功率概率预测方法 | |
CN105354620A (zh) | 一种风机发电功率的预测方法 | |
CN104573876A (zh) | 基于时序长记忆模型的风电场短期风速预测方法 | |
CN112669168B (zh) | 一种风电功率短期预测方法 | |
CN113505938A (zh) | 一种超短期风电功率组合预测方法及*** | |
CN106526710A (zh) | 一种雾霾预测方法及装置 | |
CN103984986A (zh) | 实时校正的自学习arma模型风电功率超短期预测方法 | |
CN103984987B (zh) | 一种测风网络实时校正的arma模型风电功率超短期预测方法 | |
Zhang et al. | Assessing long-term wind conditions by combining different measure-correlate-predict algorithms | |
Wang et al. | Probabilistic power curve estimation based on meteorological factors and density LSTM | |
CN113095547B (zh) | 一种基于gra-lstm-ice模型的短期风功率预测方法 | |
CN115310648A (zh) | 一种基于多气象变量模型识别的中长期风电功率组合预测方法 | |
Karami et al. | Probabilistic neural network to quantify uncertainty of wind power estimation | |
Su et al. | A LSTM based wind power forecasting method considering wind frequency components and the wind turbine states | |
CN112070320A (zh) | 一种基于动态谐波回归的超短期风电功率预测方法及*** | |
Yang et al. | Day‐ahead wind power combination forecasting based on corrected numerical weather prediction and entropy method | |
CN115577854A (zh) | 一种基于eemd-rbf组合的分位数回归风速区间预测方法 | |
Liang et al. | PM2. 5 concentration forecasting based on data preprocessing strategy and LSTM neural network | |
CN114565186A (zh) | 一种基于广义回归神经网络的预测风速校正方法和*** | |
CN111027816B (zh) | 一种基于数据包络分析的光伏发电效率计算方法 | |
Gao et al. | Very-short-term prediction of wind speed based on chaos phase space reconstruction and NWP | |
CN113761023A (zh) | 一种基于改进广义神经网络的光伏发电短期功率预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |