CN107230977A - 基于误差修正和提升小波组合预测模型的风电功率预测方法 - Google Patents
基于误差修正和提升小波组合预测模型的风电功率预测方法 Download PDFInfo
- Publication number
- CN107230977A CN107230977A CN201710313241.5A CN201710313241A CN107230977A CN 107230977 A CN107230977 A CN 107230977A CN 201710313241 A CN201710313241 A CN 201710313241A CN 107230977 A CN107230977 A CN 107230977A
- Authority
- CN
- China
- Prior art keywords
- mrow
- error
- msubsup
- wind power
- value
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 17
- 238000013277 forecasting method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 83
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 239000010410 layer Substances 0.000 claims description 60
- 230000008569 process Effects 0.000 claims description 19
- 238000005096 rolling process Methods 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 10
- 239000011229 interlayer Substances 0.000 claims description 7
- 241001123248 Arma Species 0.000 claims description 6
- 230000002146 bilateral effect Effects 0.000 claims description 6
- 238000009434 installation Methods 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 230000032798 delamination Effects 0.000 claims description 3
- 238000003909 pattern recognition Methods 0.000 claims description 3
- 238000002203 pretreatment Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 5
- 230000007812 deficiency Effects 0.000 abstract description 4
- 230000006978 adaptation Effects 0.000 abstract description 2
- 238000013459 approach Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001447 compensatory effect Effects 0.000 description 2
- 230000001373 regressive effect Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于误差修正和提升小波组合预测模型的风电功率预测方法,本发明利用提升小波分解技术处理风电功率历史数据,不仅可以提取功率数据序列的主要特性,得到特征更加明显的各频率分量,还能达到消除噪声的效果,使其更适用于各种预测算法。而根据各高低频分量的特性来选择适应的预测模型,可以消除单一预测方法的不足,并能大大的提高预测精度。本发明利用误差分层分析方法来处理误差的修正,相对于直接用预测模型来获得下一时刻误差预测值,在误差产生大幅波动情况下,该方法更能精确的分析下一时刻误差情况及补偿力度,降低在误差的预测过程中带来的误差,从而降低整体预测方法的误差。
Description
技术领域
本发明属于电力系统预测与控制技术领域,尤其涉及一种基于误差修正和提升小波组合预测模型的风电功率预测方法。
背景技术
随着风电技术的不断发展和风电场的规模不断增大,为了保证电力系统的稳定运行和供电可靠性,必须对风电系统进行有效的规划和调度。风电本身所特有的间歇性和不确定性,增加了电网调度的难度,增加了电力企业安排电网发电机组的起停和制订机组检修计划的难度,所以需要对风电场的输出功率进行预测。只有通过对风电场发电功率进行准确的预测,才可以有效降低风力发电系统的运行成本,为电网调度运行提供可靠依据。
风电功率预测方法根据不同的预测物理量可分为两大类:(1)先预测风速,再依据风电机组或者风电场的功率曲线从而获得风电场的输出功率;(2)直接预测风场的输出功率。这里直接选择输出功率为预测目标,该类方法可以简单的分为两大类。第Ⅰ类是基于确定性时序模型的预测方法,该类方法通过找出风电功率历史数据本身在时间上的相关性来进行风电功率预测,常用方法有:卡尔曼滤波法、时间序列法(ARMA)、指数平滑法等。第Ⅱ类是基于智能模型的预测方法,其实质是根据人工智能方法提取风电功率变化特性,进而进行风电功率预测。常用的方法有:人工神经网络法、小波分析法、最小二乘支持向量机(LSSVM)回归法和模糊逻辑法等。上述方法都有各自的优点,但也有很多局限性。在预测的过程中,不同的预测对象往往具有不同的特点,在选择预测方法时也是根据不同的特点选择最适合的预测方法,以提高预测精度。当预测对象受随机性太大,单一的预测方法已经不能满足其预测精度,这就可以用多种预测方法同时预测。
而现有的风电功率预测误差分析和补偿方法主要分为两类:一类是直接通过模型预测获得下一时刻误差预测值,进而对误差进行补偿和修正;另一类是应用统计方法对某一时期风电功率预测误差的概率密度曲线进行模型拟合,根据误差概率密度特征对未来误差进行估计。在第1类方法中,研究人员通过模型预测对误差数据时间及幅值特性进行分析。这类方法在结合误差数据特征的基础上进行未来误差模型预测而得到的结果精度比较有限,尤其是在误差产生大幅波动的情况下,其分析能力具有较大局限性。在第2类方法中,通常情况下风电功率预测误差幅值概率密度被假设为服从正态分布,这种假设是根据大多数预测方法和时间尺度而得到的。但是在很多情况下,尤其是在风电并网容量较大时,正态分布并不能很好的描述误差分布。而现有研究的概率分布模型拟合精度方面存在不足,并且单纯通过概率分布拟合估计预测误差的效果有限。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于误差修正和提升小波组合预测模型的风电功率预测方法,运用提升小波来处理历史数据,不仅可以提取数据的主要特性,还能达到消除噪声的效果,使其适用于各种预测算法,而组合预测模型算法是根据数据特征来选取预测模型,则可以消除单一预测方法的不足。本发明还对预测误差进行分析来修正预测结果,消除大的误差偏差,达到风电功率预测结果更加精确的效果。
为了实现上述目的,本发明提供了一种基于误差修正和提升小波组合预测模型的风电功率预测方法,具体包括以下步骤:
步骤一,获取风电功率历史数据,对风电功率历史数据进行提升小波分解的预处理,将数据信号分解为高频部分和低频部分;
步骤二,根据高低频部分信号的特性分别选择灰色预测模型、差分自回归滑动平均(ARIMA)模型和最小二乘支持向量机(LSSVM)回归模型中的一个模型来建立其对应的预测模型;
步骤三,将预测模型的预测数据进行数据重构,得到风电功率的预测初值;
步骤四,利用风电场风电功率历史数据的预测初值和实际值得到历史预测误差值,进而得到预测误差概率密度曲线,采用改进广义误差分布模型对预测误差概率密度曲线进行拟合,得到概率密度拟合模型;
步骤五,利用概率密度拟合模型计算出在不同置信水平下的风电功率预测值置信区间,并对误差进行分层;
步骤六,根据历史预测误差值,采用差分自回归移动平均模型对历史数据的最后一个时刻t的下一个时刻t+1的误差值进行预测;
步骤七,根据t+1时刻的误差预测值及t时刻的误差值在误差层中所处位置,选择不同的补偿力度对t+1时刻的误差值进行补偿,从而修正t+1时刻的风电功率预测初值,得到t+1时刻的预测结果;
步骤八,获取最新风电功率实际值,对误差和风电功率值进行滚动预测。
进一步地,所述步骤一中,首先获取风电场风电功率的实际值Pactual和预测初值Pforecast,而实际值Pactual的数据序列可用Pa表示,设k=1,2,...,N,N表示时刻;历史预测初值Pforecast的数据序列可用Pf表示,并设
再通过提升小波分解风电功率的实际值数据,获取数据特性,使得数据更有利于建模预测。而提升小波算法通过构造双正交小波函数,使用线性、非线性或空间变化的预测和更新算子进行提升变换,而且确保了变换的可逆性。
所述步骤一中,提升小波分解的预处理步骤具体如下:
1)分裂:将风电功率实际值Pa分割成相互关联的奇偶两部分,即(偶部分)和(奇部分);
其中k=1,2,...,[N/2],[N/2]为取N/2的整数部分;
2)预测:用预测获取预测值实际值与预测值的差值d1(即高频部分)反应了两者之间的逼近程度,称为细节系数或小波系数,对应于元信号Pa的高频部分。预测过程为:
式中,预测算子P可用预测函数来表示,函数可取为中的对应数据本身,即:
则
3)更新:经过分裂步骤产生子集的某些整体特征可能与原始数据并不一致,为了保持风电功率数据的这些整体特征,需要一个更新的过程。将更新过程用算子U来代替,其过程为:
式中,s1是Pa的低频部分,更新算子U可以用函数Uk()表示,即:
Uk(d1)={d1,1/2,d1,2/2,...,d1,k/2},k=1,2,...,[N/2] (5)
经提升小波分解,可将风电功率数据Pa分解为低频部分s1和高频部分d1,对于低频数据子集s1可以再进行相同的分裂、预测和更新,把s1进一步分解成s2和d2;…;如此下去,经过n次分解后,风电功率数据Pa的小波表示为{sn,dn,dn-1,…,d1}。其中sn代表了风电功率数据的低频部分,而{dn,dn-1,…,d1}则是功率数据从低到高的高频部分系列。
进一步地,所述步骤二中,根据步骤一所得每一个高频分量和低频分量的特性分别选择灰色预测模型、差分自回归滑动平均(ARIMA)模型和最小二乘支持向量机(LSSVM)回归模型中的一个模型来建立其对应的预测模型。低频分量变化较平缓且波动小,适合选用灰色预测法预测。高频分量代表原始信号中随机性最强的突然波动且无规律可循,而最小二乘支持向量机回归模型具有很强的泛化性,可对高频分量进行训练预测。对于高频信号中的周期序列则可选用ARIMA模型进行预测。
进一步地,所述步骤三中,数据重构,即提升小波的反变换过程,可以用替代的方式来计算:
Merge()即将奇序列和偶序列合并。
进一步地,所述步骤四中,首先根据风电场历史风电功率数据中的风电功率预测初值Pforecast和实际值Pactual求出预测误差Δp,即:
再算出风电功率预测误差的标幺值x(也是相对于总装机容量的相对误差),即x=(Pforecast-Pactual)/Pbase=Δp/Pbase,其中Pbase是样本数据系统所接入的风电装机总容量。再根据每个预测点的误差标幺值的概率密度来得到其概率密度曲线。
使用改进广义误差分布模型来拟合预测误差概率密度曲线,改进广义误差分布模型概率密度函数为:
其中v和λ为形状参数;Γ(·)为伽马函数。公式(8)中计算得到参数λ决定了曲线总体形状的平坦与陡峭程度,斜度参数α剥离了曲线斜度和峰度的关系使得曲线变化更具有灵活性,位置参数μ可以使模型具有拟合带偏度曲线的能力。
该模型的参数估计的方法确定:由于本文所研究的概率分布模型大部分为指数形式,故采用极大似然估计(MLE)法对模型参数进行估计。这样可以通过最大化对数似然值(log-likelihood,LL),将样本数据与模型参数的关系进行转变,便于参数估计。
进一步地,所述步骤五中,在确定误差概率密度的拟合模型后,根据分析对象数据特征选取一低一高两个不同的置信度水平(如85%和95%)作为分层标准,并根据预测误差计算风电功率预测初值的置信区间。根据风电功率实际值处于这两个置信区间中的位置,实现对误差进行分层。当实际值处于置信度小的置信区间内时,说明此时误差较小,将该误差层称为小误差层;当实际值处于置信度大的置信区间以外时,说明此时误差较大,将该层称为大误差层;当实际值处于两个置信区间临界值之间时,误差处于中等水平,将该层称为中误差层。由此,可以得到根据历史预测误差构建的预测误差分层体系。
在大多数情况下,预测误差概率密度曲线是可以看作对称的,所以计算单侧概率密度累积值F可以判断总体误差水平。
进一步地,所述步骤六中,根据历史预测误差值,采用差分自回归移动平均(ARIMA)模型对历史数据的最后一个时刻t的下一个时刻t+1的误差值进行预测;ARIMA是由三部分组成:自回归项(AR)、差分项(I)和滑动平均模型(MA),是在自回归滑动平均(ARMA)模型的基础上提出来的。
ARMA数学表达式如下:
式中,Δpt+1是t+1时刻风电功率预测误差值的预测值;表示自回归项AR,为过去观测值的线性组合;aj为常数,Δpt+1-j为t+1-j时刻的观测值;bl为常数,ht+1-l为白噪声序列(一个随机序列,期望为0,方差为常数);表示白噪声序列的滑动平均项MA。
一个风电功率预测误差值时间序列在某时刻可以用p个历史观测值的线性组合加上一个白噪声序列的q项滑动平均来表示,则该时间序列为ARMA(p,q)过程。
使用自回归滑动平均模型时,首先判断风电功率预测误差值序列Δp是否为平稳序列,若非平稳序列,则通常使用差分使之变为平稳序列,即ARIMA过程。之后根据数据的自相关和偏相关系数进行模式识别,再由最小信息准则(AIC)进行模型定阶。最后由ARMA模型的表达公式,即可估计出模型的参数值。
进一步地,所述步骤七中,由于在t+1刻误差值的预测过程中,已经考虑过去几个时刻的风电功率预测历史误差值。因而,在进行分层补偿时只需要根据t+1时刻的误差预测值及t时刻的误差值在误差层中所处位置,选择不同的补偿力度对t+1时刻的误差值进行补偿。
t+1时刻误差预测值相对于t时刻误差历史值可能会在同侧(正误差或负误差)的同一层之内或者在多误差层之间波动。为了防止产生补偿方法的误判,对误差预测值与误差历史值在不同层和不同侧做相应的补偿办法。补偿情况基于如下:
1、当误差预测值在单侧波动
当t时刻误差历史值Δpt和误差预测值Δpt+1均处于单侧同一误差层内时,此时的误差波动较小,只需要对误差预测值进行等幅反向补偿即可。即当两者都处于小误差层时,说明误差较小,则不进行补偿;当两者都处于中误差层或是大误差层时则对误差预测值进行等幅反向补偿。
2、误差预测值在层间波动
误差预测值可能会在正负误差两侧各3个误差层之间任意波动,因此对于误差预测值的层间波动不能仅考虑误差预测值(t+1时刻)的幅值大小,还应该考虑未来误差(t+2时刻)的发展趋势。此时,本文引入Δpt和Δpt+1连线斜率作为误差发展趋势的衡量标准,的定义如下:
式中|ψ1-ψ2|为由拟合模型和置信度水平确定的单侧置信区间临界值之差的绝对值。
当Δpt和Δpt+1均在双侧的小误差层中时,对风电功率预测误差不进行补偿。当Δpt+1处于双侧任一小误差层而Δpt不在小误差层中时,虽然误差仍具有较大的变化趋势,但是无法排除未来误差持续处于小误差层的情况,因此也不再进行补偿。除了以上两种特殊情况外,误差预测值均按照表1中的补偿方式和补偿幅度进行补偿。
表1误差预测值处于层间波动时的补偿方法
当值较大时,说明误差变化幅度较大。此时t+2时刻的超短期风电功率预测误差进一步增大的可能性较大,因此采用了误差预测值的反向过补偿对功率预测初值进行修正,以防出现补偿力度不足的情况。过补偿倍数应根据风电功率预测误差的波动情况来确定,并且可以通过调节过补偿倍数来对波动较大误差的补偿力度进行调控。由于以上原因,并结合本研究所用数据特点,选取1.5倍过补偿来进行方法的说明。通过以上方法对误差进行滚动分析和补偿,实现对由预测模型得到的风电功率预测初值的进一步修正,从而可以提高超短期风电功率预测值的精度。
进一步地,所述步骤八中,获取下一时刻的风电功率实际值,将其和预测值一起当作新的历史数据,并丢弃历史数据中的第一个历史数据。则第i+1次滚动预测所利用的数据序列为其中为历史数据,为前i次滚动预测所得到的预测值。
本发明的有益效果为:
(1)该发明利用提升小波分解技术处理风电功率历史数据,不仅可以提取功率数据序列的主要特性,得到特征更加明显的各频率分量,还能达到消除噪声的效果,使其更适用于各种预测算法。而根据各高低频分量的特性来选择适应的预测模型,可以消除单一预测方法的不足,并能大大的提高预测精度。
(2)该发明利用误差分层分析方法来处理误差的修正,相对于直接用预测模型来获得下一时刻误差预测值,在误差产生大幅波动情况下,该方法更能精确的分析下一时刻误差情况及补偿力度,降低在误差的预测过程中带来的误差,从而降低整体预测方法的误差。
附图说明
图1为短期风电功率组合预测方法流程图;
图2为风电功率序列提升小波分解的流程图;
具体实施方式
下面结合附图与实施例对本发明的具体实施方式作进一步详细说明。
本发明提供的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,具体包括以下步骤:
步骤一,获取风电功率历史数据,对风电功率历史数据进行提升小波分解的预处理,将数据信号分解为高频部分和低频部分;
步骤二,根据高低频部分信号的特性分别选择灰色预测模型、差分自回归滑动平均(ARIMA)模型和最小二乘支持向量机(LSSVM)回归模型中的一个模型来建立其对应的预测模型;
步骤三,将预测模型的预测数据进行数据重构,得到风电功率的预测初值;
步骤四,利用风电场风电功率历史数据的预测初值和实际值得到历史预测误差值,进而得到预测误差概率密度曲线,采用改进广义误差分布模型对预测误差概率密度曲线进行拟合,得到概率密度拟合模型;
步骤五,利用概率密度拟合模型计算出在不同置信水平下的风电功率预测值置信区间,并对误差进行分层;
步骤六,根据历史预测误差值,采用差分自回归移动平均模型对历史数据的最后一个时刻t的下一个时刻t+1的误差值进行预测;
步骤七,根据t+1时刻的误差预测值及t时刻的误差值在误差层中所处位置,选择不同的补偿力度对t+1时刻的误差值进行补偿,从而修正t+1时刻的风电功率预测初值,得到t+1时刻的预测结果;
步骤八,获取最新风电功率实际值,对误差和风电功率值进行滚动预测。
进一步地,所述步骤一中,首先获取风电场风电功率的实际值Pactual和预测初值Pforecast,而实际值Pactual的数据序列可用Pa表示,设k=1,2,...,N,N表示时刻;历史预测初值Pforecast的数据序列可用Pf表示,并设
再通过提升小波分解风电功率的实际值数据,获取数据特性,使得数据更有利于建模预测。而提升小波算法通过构造双正交小波函数,使用线性、非线性或空间变化的预测和更新算子进行提升变换,而且确保了变换的可逆性。
所述步骤一中,提升小波分解的预处理步骤具体如下:
1)分裂:将风电功率实际值Pa分割成相互关联的奇偶两部分,即(偶部分)和(奇部分);
其中k=1,2,...,[N/2],[N/2]为取N/2的整数部分;
2)预测:用预测获取预测值实际值与预测值的差值d1(即高频部分)反应了两者之间的逼近程度,称为细节系数或小波系数,对应于元信号Pa的高频部分。预测过程为:
式中,预测算子P可用预测函数来表示,函数可取为中的对应数据本身,即:
则
3)更新:经过分裂步骤产生子集的某些整体特征可能与原始数据并不一致,为了保持风电功率数据的这些整体特征,需要一个更新的过程。将更新过程用算子U来代替,其过程为:
式中,s1是Pa的低频部分,更新算子U可以用函数Uk()表示,即:
Uk(d1)={d1,1/2,d1,2/2,...,d1,k/2},k=1,2,...,[N/2] (5)
经提升小波分解,可将风电功率数据Pa分解为低频部分s1和高频部分d1,对于低频数据子集s1可以再进行相同的分裂、预测和更新,把s1进一步分解成s2和d2;…;如此下去,经过n次分解后,风电功率数据Pa的小波表示为{sn,dn,dn-1,…,d1}。其中sn代表了风电功率数据的低频部分,而{dn,dn-1,…,d1}则是功率数据从低到高的高频部分系列。
进一步地,所述步骤二中,根据步骤一所得每一个高频分量和低频分量的特性分别选择灰色预测模型、差分自回归滑动平均(ARIMA)模型和最小二乘支持向量机(LSSVM)回归模型中的一个模型来建立其对应的预测模型。低频分量变化较平缓且波动小,适合选用灰色预测法预测。高频分量代表原始信号中随机性最强的突然波动且无规律可循,而最小二乘支持向量机回归模型具有很强的泛化性,可对高频分量进行训练预测。对于高频信号中的周期序列则可选用ARIMA模型进行预测。
进一步地,所述步骤三中,数据重构,即提升小波的反变换过程,可以用替代的方式来计算:
Merge()即将奇序列和偶序列合并。
进一步地,所述步骤四中,首先根据风电场历史风电功率数据中的风电功率预测初值Pforecast和实际值Pactual求出预测误差Δp,即:
再算出风电功率预测误差的标幺值x(也是相对于总装机容量的相对误差),即x=(Pforecast-Pactual)/Pbase=Δp/Pbase,其中Pbase是样本数据系统所接入的风电装机总容量。再根据每个预测点的误差标幺值的概率密度来得到其概率密度曲线。
使用改进广义误差分布模型来拟合预测误差概率密度曲线,改进广义误差分布模型概率密度函数为:
其中v和λ为形状参数;Γ(·)为伽马函数。公式(8)中计算得到参数λ决定了曲线总体形状的平坦与陡峭程度,斜度参数α剥离了曲线斜度和峰度的关系使得曲线变化更具有灵活性,位置参数μ可以使模型具有拟合带偏度曲线的能力。
该模型的参数估计的方法确定:由于本文所研究的概率分布模型大部分为指数形式,故采用极大似然估计(MLE)法对模型参数进行估计。这样可以通过最大化对数似然值(log-likelihood,LL),将样本数据与模型参数的关系进行转变,便于参数估计。
进一步地,所述步骤五中,在确定误差概率密度的拟合模型后,根据分析对象数据特征选取一低一高两个不同的置信度水平(如85%和95%)作为分层标准,并根据预测误差计算风电功率预测初值的置信区间。根据风电功率实际值处于这两个置信区间中的位置,实现对误差进行分层。当实际值处于置信度小的置信区间内时,说明此时误差较小,将该误差层称为小误差层;当实际值处于置信度大的置信区间以外时,说明此时误差较大,将该层称为大误差层;当实际值处于两个置信区间临界值之间时,误差处于中等水平,将该层称为中误差层。由此,可以得到根据历史预测误差构建的预测误差分层体系。
在大多数情况下,预测误差概率密度曲线是可以看作对称的,所以计算单侧概率密度累积值F可以判断总体误差水平。
进一步地,所述步骤六中,根据历史预测误差值,采用差分自回归移动平均(ARIMA)模型对历史数据的最后一个时刻t的下一个时刻t+1的误差值进行预测;ARIMA是由三部分组成:自回归项(AR)、差分项(I)和滑动平均模型(MA),是在自回归滑动平均(ARMA)模型的基础上提出来的。
ARMA数学表达式如下:
式中,Δpt+1是t+1时刻风电功率预测误差值的预测值;表示自回归项AR,为过去观测值的线性组合;aj为常数,Δpt+1-j为t+1-j时刻的观测值;bl为常数,ht+1-l为白噪声序列(一个随机序列,期望为0,方差为常数);表示白噪声序列的滑动平均项MA。
一个风电功率预测误差值时间序列在某时刻可以用p个历史观测值的线性组合加上一个白噪声序列的q项滑动平均来表示,则该时间序列为ARMA(p,q)过程。
使用自回归滑动平均模型时,首先判断风电功率预测误差值序列Δp是否为平稳序列,若非平稳序列,则通常使用差分使之变为平稳序列,即ARIMA过程。之后根据数据的自相关和偏相关系数进行模式识别,再由最小信息准则(AIC)进行模型定阶。最后由ARMA模型的表达公式,即可估计出模型的参数值。
进一步地,所述步骤七中,由于在t+1刻误差值的预测过程中,已经考虑过去几个时刻的风电功率预测历史误差值。因而,在进行分层补偿时只需要根据t+1时刻的误差预测值及t时刻的误差值在误差层中所处位置,选择不同的补偿力度对t+1时刻的误差值进行补偿。
t+1时刻误差预测值相对于t时刻误差历史值可能会在同侧(正误差或负误差)的同一层之内或者在多误差层之间波动。为了防止产生补偿方法的误判,对误差预测值与误差历史值在不同层和不同侧做相应的补偿办法。补偿情况基于如下:
1、当误差预测值在单侧波动
当t时刻误差历史值Δpt和误差预测值Δpt+1均处于单侧同一误差层内时,此时的误差波动较小,只需要对误差预测值进行等幅反向补偿即可。即当两者都处于小误差层时,说明误差较小,则不进行补偿;当两者都处于中误差层或是大误差层时则对误差预测值进行等幅反向补偿。
2、误差预测值在层间波动
误差预测值可能会在正负误差两侧各3个误差层之间任意波动,因此对于误差预测值的层间波动不能仅考虑误差预测值(t+1时刻)的幅值大小,还应该考虑未来误差(t+2时刻)的发展趋势。此时,本文引入Δpt和Δpt+1连线斜率作为误差发展趋势的衡量标准,的定义如下:
式中|ψ1-ψ2|为由拟合模型和置信度水平确定的单侧置信区间临界值之差的绝对值。
当Δpt和Δpt+1均在双侧的小误差层中时,对风电功率预测误差不进行补偿。当Δpt+1处于双侧任一小误差层而Δpt不在小误差层中时,虽然误差仍具有较大的变化趋势,但是无法排除未来误差持续处于小误差层的情况,因此也不再进行补偿。除了以上两种特殊情况外,误差预测值均按照表1中的补偿方式和补偿幅度进行补偿。
表1误差预测值处于层间波动时的补偿方法
当值较大时,说明误差变化幅度较大。此时t+2时刻的超短期风电功率预测误差进一步增大的可能性较大,因此采用了误差预测值的反向过补偿对功率预测初值进行修正,以防出现补偿力度不足的情况。过补偿倍数应根据风电功率预测误差的波动情况来确定,并且可以通过调节过补偿倍数来对波动较大误差的补偿力度进行调控。由于以上原因,并结合本研究所用数据特点,选取1.5倍过补偿来进行方法的说明。通过以上方法对误差进行滚动分析和补偿,实现对由预测模型得到的风电功率预测初值的进一步修正,从而可以提高超短期风电功率预测值的精度。
进一步地,所述步骤八中,获取下一时刻的风电功率实际值,将其和预测值一起当作新的历史数据,并丢弃历史数据中的第一个历史数据。则第i+1次滚动预测所利用的数据序列为其中为历史数据,为前i次滚动预测所得到的预测值。
实施例:选用某风电场5天的风电功率数据作为样本,利用第6天的风电功率数据进行虚拟预测对比分析,风功率数据的时间分辨率是每15分钟一个采样数据点(一天96点)。对风电场功率序列进行基于提升小波的组合模型预测及误差修正,如图1所示,包括以下步骤:
步骤(1):获取风电功率历史数据序列:将选取以上的采样数据点,获得历史风电功率序列其中为训练数据,为测试数据。
步骤(2):利用提升小波分解历史风电功率序列得到历史风电功率序列的低频分量和高频分量。
根据仿真实验,对历史风电功率序列进行三层提升小波的分解比较适合预测方法的训练和预测。所得到的频率分量为{sn,dn,dn-1,…,d1},其中sn代表了风电功率数据的低频部分,而{dn,dn-1,…,d1}则是功率数据从低到高的高频部分系列。历史风电功率时间序列的提升小波分解如图2所示。
步骤(3):对每一个频率分量进行分析,并分别选取相适应的预测模型分别对风电功率序列每一个高低频分量进行预测:低频分量变化较平缓且波动小,适合选用灰色预测法预测。高频分量代表原始信号中随机性最强的突然波动且无规律可循,而最小二乘支持向量机具有很强的泛化性,可对高频分量进行训练预测。对于高频信号中的周期序列则可选用ARIMA模型进行预测。
步骤(4):进行提升小波逆运算即数据重构得到下一时刻的风电功率预测初值
步骤(5):通过对历史风电功率实际值的前480个数据的预测,按照步骤(3)和步骤(4)获得480个历史风电功率的预测初值,并根据这480个历史风电功率的实际值和预测初值求出预测误差序列Δp,Δp={Δp1,Δp2,...,Δpt},t=480;
且根据预测误差算出预测误差的标幺值x,且x=Δp/pbase,其中Pbase是样本数据系统所接入的风电装机总容量。并根据每个预测点的误差标幺值的概率密度来得到其概率密度曲线。
步骤(6):使用改进广义误差分布模型来拟合预测误差概率密度曲线。改进广义误差分布模型概率密度函数为:
其中v和λ为形状参数;Γ(·)为伽马函数。
步骤(7):利用概率密度拟合模型计算出在不同置信水平下的风电功率预测初值置信区间,并利用风电功率预测初值置信区间与实际值之间的关系来对误差进行分层;置信度水平参考值如表2所示。
表2置信度水平参考值
根据样本数据预测误差水平,本例选取95%和80%置信度水平作为分层依据,利用改进误差分布模型进行计算,可得95%置信度水平所对应的风电功率预测误差标幺值的绝对值及其对应误差区间,和80%置信度水平所对应的风电功率预测误差标幺值的绝对值及其对应误差区间。根据风电功率误差标幺值所对应的功率预测初值及误差区间,可得出风电功率预测初值在一定置信水平的置信区间。再根据风电功率误差标幺值所对应的功率实际值,依照步骤五中描述的误差分层方法实现对误差的分层,依次分为大误差层、中误差层和小误差层。
步骤(8):选取步骤(5)中所获取的历史预测误差序列Δp,利用ARIMA模型对其预测获得t+1时刻误差预测值Δpt+1,即第6天的第一个时间点的风电功率预测误差的初步预测值。
步骤(9):根据风电功率预测误差历史值的最后一个数据Δpt和步骤(8)所得的t+1时刻风电功率预测的误差预测值Δpt+1在误差层级之间所处位置的情况对误差预测值Δpt+1进行补偿。具体补偿按照步骤七中所述的情况进行对应补偿。
步骤(10):将步骤(9)中所得的补偿后的误差对步骤(4)中所得的下一时刻风电功率预测值进行修正,则可得到修正后的风电功率预测值
步骤(11):依照步骤八所述,获得下一时刻风电功率的实际值,按上述步骤及步骤八的方法进行滚动预测,获得第6天的96个预测点的风电功率预测值,并与历史实际值中第6天的数据进行对比,可得到预测方法的预测精度。
最后应当说的是,以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (9)
1.一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,该方法包括以下步骤:
步骤一,获取风电功率历史数据,对风电功率历史数据进行提升小波分解的预处理,将数据信号分解为高频部分和低频部分;
步骤二,根据高低频部分信号的特性分别选择灰色预测模型、差分自回归滑动平均(ARIMA)模型和最小二乘支持向量机(LSSVM)回归模型中的一个模型来建立其对应的预测模型;
步骤三,将预测模型的预测数据进行数据重构,得到风电功率的预测初值;
步骤四,利用风电场风电功率历史数据的预测初值和实际值得到历史预测误差值,进而得到预测误差概率密度曲线,采用改进广义误差分布模型对预测误差概率密度曲线进行拟合,得到概率密度拟合模型;
步骤五,利用概率密度拟合模型计算出在不同置信水平下的风电功率预测值置信区间,并对误差进行分层;
步骤六,根据历史预测误差值,采用差分自回归移动平均模型对历史数据的最后一个时刻t的下一个时刻t+1的误差值进行预测;
步骤七,根据t+1时刻的误差预测值及t时刻的误差值在误差层中所处位置,选择不同的补偿力度对t+1时刻的误差值进行补偿,从而修正t+1时刻的风电功率预测初值,得到t+1时刻的预测结果;
步骤八,获取最新风电功率实际值,对误差和风电功率值进行滚动预测。
2.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤一中,提升小波分解的预处理步骤具体如下:
1)分裂:将风电功率实际值Pa分割成相互关联的奇偶两部分,即偶部分和奇部分
<mrow>
<mi>S</mi>
<mi>p</mi>
<mi>l</mi>
<mi>i</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mi>a</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>o</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<mo>{</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mn>4</mn>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
<mi>a</mi>
</msubsup>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>o</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<mo>{</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mn>3</mn>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中k=1,2,...,[N/2],[N/2]为取N/2的整数部分;
2)预测:用预测获取预测值实际值与预测值的差值d1反应了两者之间的逼近程度,对应于元信号Pa的高频部分。预测过程为:
<mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mi>o</mi>
<mi>a</mi>
</msubsup>
<mo>-</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,预测算子P可用预测函数来表示,函数可取为中的对应数据本身,即:
<mrow>
<msub>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>{</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>a</mi>
</msubsup>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
则
3)更新:将更新过程用算子U来代替,其过程为:
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>+</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
1
式中,s1是Pa的低频部分,更新算子U可以用函数Uk()表示,即:
Uk(d1)={d1,1/2,d1,2/2,...,d1,k/2},k=1,2,...,[N/2] (5)
经过n次分解后,风电功率数据Pa的小波表示为{sn,dn,dn-1,…,d1}。其中,sn代表了风电功率数据的低频部分,而{dn,dn-1,…,d1}则是功率数据从低到高的高频部分系列。
3.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤二中,低频分量适合选用灰色预测法预测。高频分量适合最小二乘支持向量机回归模型。对于高频信号中的周期序列则可选用ARIMA模型进行预测。
4.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤三中,数据重构,即提升小波的反变换过程,可以用替代的方式来计算:
<mrow>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>o</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>P</mi>
<mi>a</mi>
</msup>
<mo>=</mo>
<mi>M</mi>
<mi>e</mi>
<mi>r</mi>
<mi>g</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>e</mi>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>o</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Merge()即将奇序列和偶序列合并。
5.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤四中,首先根据风电场历史风电功率数据中的风电功率预测初值Pforecast和实际值Pactual求出预测误差Δp,即:
<mrow>
<mi>&Delta;</mi>
<mi>p</mi>
<mo>=</mo>
<msup>
<mi>P</mi>
<mi>f</mi>
</msup>
<mo>-</mo>
<msup>
<mi>P</mi>
<mi>a</mi>
</msup>
<mo>=</mo>
<mo>{</mo>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mi>f</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mi>f</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mi>a</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>t</mi>
<mi>f</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>t</mi>
<mi>a</mi>
</msubsup>
<mo>}</mo>
<mo>=</mo>
<mo>{</mo>
<msub>
<mi>&Delta;p</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>&Delta;p</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>&Delta;p</mi>
<mi>t</mi>
</msub>
<mo>}</mo>
<mo>;</mo>
</mrow>
再算出风电功率预测误差的标幺值x,即x=(Pforecast-Pactual)/Pbase=Δp/Pbase,其中Pbase是样本数据系统所接入的风电装机总容量;再根据每个预测点的误差标幺值的概率密度来得到其概率密度曲线。
使用改进广义误差分布模型来拟合预测误差概率密度曲线,改进广义误差分布模型概率密度函数为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>;</mo>
<mi>v</mi>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>v</mi>
<mrow>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>v</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mo>|</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<mi>&mu;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<msup>
<mo>|</mo>
<mi>&alpha;</mi>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>{</mo>
<mfrac>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mo>&lsqb;</mo>
<mn>2</mn>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<mi>v</mi>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>3</mn>
<mi>v</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>3</mn>
<mi>v</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中v和λ为形状参数;Γ(·)为伽马函数,α为斜度参数,μ为位置参数。
采用极大似然估计(MLE)法对模型参数进行估计。
6.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤五中,根据分析对象数据特征选取一低一高两个不同的置信度水平作为分层标准,并根据预测误差计算风电功率预测初值的置信区间。根据风电功率实际值处于这两个置信区间中的位置,实现对误差进行分层。当实际值处于置信度小的置信区间内时,说明此时误差较小,将该误差层称为小误差层;当实际值处于置信度大的置信区间以外时,说明此时误差较大,将该层称为大误差层;当实际值处于两个置信区间临界值之间时,误差处于中等水平,将该层称为中误差层。由此,可以得到根据历史预测误差构建的预测误差分层体系。
7.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤六中,根据历史预测误差值,采用差分自回归移动平均模型对历史数据的最后一个时刻t的下一个时刻t+1的误差值进行预测;ARIMA是由三部分组成:自回归项(AR)、差分项(I)和滑动平均模型(MA),ARMA数学表达式如下:
<mrow>
<msub>
<mi>&Delta;p</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<msub>
<mi>&Delta;p</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>-</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>q</mi>
</munderover>
<msub>
<mi>b</mi>
<mi>l</mi>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>-</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Δpt+1是t+1时刻风电功率预测误差值的预测值;表示自回归项AR,为过去观测值的线性组合;aj为常数,Δpt+1-j为t+1-j时刻的观测值;bl为常数,ht+1-l为白噪声序列;表示白噪声序列的滑动平均项MA。
使用自回归滑动平均模型时,首先判断风电功率预测误差值序列Δp是否为平稳序列,若非平稳序列,则通常使用差分使之变为平稳序列,即ARIMA过程。之后根据数据的自相关和偏相关系数进行模式识别,再由最小信息准则(AIC)进行模型定阶。最后由ARMA模型的表达公式,即可估计出模型的参数值。
8.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤七中,在进行分层补偿时只需要根据t+1时刻的误差预测值及t时刻的误差值在误差层中所处位置,选择不同的补偿力度对t+1时刻的误差值进行补偿,补偿情况基于如下:
1)当误差预测值在单侧波动:当t时刻误差历史值Δpt和误差预测值Δpt+1均处于单侧同一误差层内时,此时的误差波动较小,只需要对误差预测值进行等幅反向补偿即可。即当两者都处于小误差层时,说明误差较小,则不进行补偿;当两者都处于中误差层或是大误差层时则对误差预测值进行等幅反向补偿。
2)误差预测值在层间波动:误差预测值可能会在正负误差两侧各3个误差层之间任意波动,引入Δpt和Δpt+1连线斜率作为误差发展趋势的衡量标准,的定义如下:
式中|ψ1-ψ2|为由拟合模型和置信度水平确定的单侧置信区间临界值之差的绝对值。
当Δpt和Δpt+1均在双侧的小误差层中时,对风电功率预测误差不进行补偿。当Δpt+1处于双侧任一小误差层而Δpt不在小误差层中时,不再进行补偿。
除了以上两种特殊情况外,误差预测值按下表中的补偿方式和补偿幅度进行补偿:
。
9.根据权利要求1所述的一种基于误差修正和提升小波组合预测模型的风电功率预测方法,其特征在于,所述步骤八中,获取下一时刻的风电功率实际值,将其和预测值一起当作新的历史数据,并丢弃历史数据中的第一个历史数据。则第i+1次滚动预测所利用的数据序列为其中为历史数据,为前i次滚动预测所得到的预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710313241.5A CN107230977A (zh) | 2017-05-05 | 2017-05-05 | 基于误差修正和提升小波组合预测模型的风电功率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710313241.5A CN107230977A (zh) | 2017-05-05 | 2017-05-05 | 基于误差修正和提升小波组合预测模型的风电功率预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107230977A true CN107230977A (zh) | 2017-10-03 |
Family
ID=59933194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710313241.5A Pending CN107230977A (zh) | 2017-05-05 | 2017-05-05 | 基于误差修正和提升小波组合预测模型的风电功率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107230977A (zh) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108173686A (zh) * | 2017-12-26 | 2018-06-15 | 北京工业大学 | 一种基于arima和小波变换的云数据中心请求流预测方法 |
CN108802535A (zh) * | 2018-06-27 | 2018-11-13 | 全球能源互联网研究院有限公司 | 筛选方法、主干扰源识别方法及装置、服务器及存储介质 |
CN109640351A (zh) * | 2019-01-25 | 2019-04-16 | 南京邮电大学 | 一种基站流量的联合预测方法 |
CN109740111A (zh) * | 2018-12-24 | 2019-05-10 | 华北科技学院 | 对地电场值预测算法 |
CN109816165A (zh) * | 2019-01-16 | 2019-05-28 | 国能日新科技股份有限公司 | 风电超短期功率预测方法及系统 |
CN109981327A (zh) * | 2017-12-28 | 2019-07-05 | 中移信息技术有限公司 | 一种业务量的预测方法及系统 |
CN110571850A (zh) * | 2019-08-28 | 2019-12-13 | 中国农业大学 | 一种风电场功率波动轨迹预测和校正控制方法 |
CN110707744A (zh) * | 2019-10-09 | 2020-01-17 | 广东电网有限责任公司电网规划研究中心 | 风电场集群发电状态的监测方法、装置及存储介质 |
CN110717610A (zh) * | 2018-07-11 | 2020-01-21 | 华北电力大学(保定) | 一种基于数据挖掘的风电功率预测方法 |
CN110782059A (zh) * | 2018-07-31 | 2020-02-11 | 北京金风科创风电设备有限公司 | 风力发电机组关键性能指标预测方法、装置、设备及介质 |
CN111105005A (zh) * | 2019-12-03 | 2020-05-05 | 广东电网有限责任公司 | 一种风电功率预测方法 |
CN112489418A (zh) * | 2020-10-22 | 2021-03-12 | 浙江交通职业技术学院 | 一种基于道路路段行程时间预测模型的路段行程时间动态误差修正方法 |
CN113033904A (zh) * | 2021-04-02 | 2021-06-25 | 合肥工业大学 | 基于s变换的风电功率预测误差分析与分类方法 |
CN113258565A (zh) * | 2021-05-11 | 2021-08-13 | 广东电网有限责任公司韶关供电局 | 一种频率调节方法、装置、设备和存储介质 |
CN113505158A (zh) * | 2021-07-16 | 2021-10-15 | 瑞幸咖啡信息技术(厦门)有限公司 | 一种时间序列的异常检测方法、装置、设备及存储介质 |
CN113792032A (zh) * | 2021-08-09 | 2021-12-14 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于正态分布误差修正的测风数据塔影效应分析方法 |
CN113809772A (zh) * | 2021-09-17 | 2021-12-17 | 国网河南省电力公司电力科学研究院 | 提升应对次时间尺度风电不确定性的安全性的方法及装置 |
CN113991638A (zh) * | 2021-08-31 | 2022-01-28 | 华能(福建)能源开发有限公司 | 一种针对不同地点的新能源场站发电功率的预测方法 |
CN114897245A (zh) * | 2022-05-12 | 2022-08-12 | 国网湖北省电力有限公司电力科学研究院 | 一种基于长短期记忆神经网络的规模风电超短期功率预测误差修正方法 |
CN116448263A (zh) * | 2023-06-16 | 2023-07-18 | 山东德圣源新材料有限公司 | 一种勃姆石生产设备运行状态检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102102626A (zh) * | 2011-01-30 | 2011-06-22 | 华北电力大学 | 风电场短期功率预测方法 |
CN102855385A (zh) * | 2012-07-31 | 2013-01-02 | 上海交通大学 | 一种风力发电短期负荷预测方法 |
CN104376368A (zh) * | 2014-08-19 | 2015-02-25 | 上海交通大学 | 一种基于频域分解的风力发电短期负荷预测方法及装置 |
CN105787606A (zh) * | 2016-03-24 | 2016-07-20 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于超短期负荷预测的电力调度在线趋势预警系统 |
US10135253B2 (en) * | 2000-12-29 | 2018-11-20 | Abb Schweiz Ag | System, method and computer program product for enhancing commercial value of electrical power produced from a renewable energy power production facility |
-
2017
- 2017-05-05 CN CN201710313241.5A patent/CN107230977A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10135253B2 (en) * | 2000-12-29 | 2018-11-20 | Abb Schweiz Ag | System, method and computer program product for enhancing commercial value of electrical power produced from a renewable energy power production facility |
CN102102626A (zh) * | 2011-01-30 | 2011-06-22 | 华北电力大学 | 风电场短期功率预测方法 |
CN102855385A (zh) * | 2012-07-31 | 2013-01-02 | 上海交通大学 | 一种风力发电短期负荷预测方法 |
CN104376368A (zh) * | 2014-08-19 | 2015-02-25 | 上海交通大学 | 一种基于频域分解的风力发电短期负荷预测方法及装置 |
CN105787606A (zh) * | 2016-03-24 | 2016-07-20 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于超短期负荷预测的电力调度在线趋势预警系统 |
Non-Patent Citations (2)
Title |
---|
叶林 等: "超短期风电功率预测误差数值特性分层分析方法", 《中国电机工程学报》 * |
李霄: "基于提升小波和最小二乘支持向量机的风电功率预测", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108173686A (zh) * | 2017-12-26 | 2018-06-15 | 北京工业大学 | 一种基于arima和小波变换的云数据中心请求流预测方法 |
CN109981327A (zh) * | 2017-12-28 | 2019-07-05 | 中移信息技术有限公司 | 一种业务量的预测方法及系统 |
CN108802535A (zh) * | 2018-06-27 | 2018-11-13 | 全球能源互联网研究院有限公司 | 筛选方法、主干扰源识别方法及装置、服务器及存储介质 |
CN108802535B (zh) * | 2018-06-27 | 2020-06-26 | 全球能源互联网研究院有限公司 | 筛选方法、主干扰源识别方法及装置、服务器及存储介质 |
CN110717610A (zh) * | 2018-07-11 | 2020-01-21 | 华北电力大学(保定) | 一种基于数据挖掘的风电功率预测方法 |
CN110717610B (zh) * | 2018-07-11 | 2023-10-31 | 华北电力大学(保定) | 一种基于数据挖掘的风电功率预测方法 |
CN110782059A (zh) * | 2018-07-31 | 2020-02-11 | 北京金风科创风电设备有限公司 | 风力发电机组关键性能指标预测方法、装置、设备及介质 |
CN109740111A (zh) * | 2018-12-24 | 2019-05-10 | 华北科技学院 | 对地电场值预测算法 |
CN109740111B (zh) * | 2018-12-24 | 2023-09-22 | 华北科技学院 | 对地电场值预测方法 |
CN109816165A (zh) * | 2019-01-16 | 2019-05-28 | 国能日新科技股份有限公司 | 风电超短期功率预测方法及系统 |
CN109816165B (zh) * | 2019-01-16 | 2020-09-25 | 国能日新科技股份有限公司 | 风电超短期功率预测方法及系统 |
CN109640351A (zh) * | 2019-01-25 | 2019-04-16 | 南京邮电大学 | 一种基站流量的联合预测方法 |
CN110571850A (zh) * | 2019-08-28 | 2019-12-13 | 中国农业大学 | 一种风电场功率波动轨迹预测和校正控制方法 |
CN110571850B (zh) * | 2019-08-28 | 2020-11-24 | 中国农业大学 | 一种风电场功率波动轨迹预测和校正控制方法 |
CN110707744A (zh) * | 2019-10-09 | 2020-01-17 | 广东电网有限责任公司电网规划研究中心 | 风电场集群发电状态的监测方法、装置及存储介质 |
CN111105005B (zh) * | 2019-12-03 | 2023-04-07 | 广东电网有限责任公司 | 一种风电功率预测方法 |
CN111105005A (zh) * | 2019-12-03 | 2020-05-05 | 广东电网有限责任公司 | 一种风电功率预测方法 |
CN112489418B (zh) * | 2020-10-22 | 2022-10-14 | 浙江交通职业技术学院 | 一种基于道路路段行程时间预测模型的路段行程时间动态误差修正方法 |
CN112489418A (zh) * | 2020-10-22 | 2021-03-12 | 浙江交通职业技术学院 | 一种基于道路路段行程时间预测模型的路段行程时间动态误差修正方法 |
CN113033904A (zh) * | 2021-04-02 | 2021-06-25 | 合肥工业大学 | 基于s变换的风电功率预测误差分析与分类方法 |
CN113033904B (zh) * | 2021-04-02 | 2022-09-13 | 合肥工业大学 | 基于s变换的风电功率预测误差分析与分类方法 |
CN113258565A (zh) * | 2021-05-11 | 2021-08-13 | 广东电网有限责任公司韶关供电局 | 一种频率调节方法、装置、设备和存储介质 |
CN113505158B (zh) * | 2021-07-16 | 2024-02-06 | 瑞幸咖啡信息技术(厦门)有限公司 | 一种时间序列的异常检测方法、装置、设备及存储介质 |
CN113505158A (zh) * | 2021-07-16 | 2021-10-15 | 瑞幸咖啡信息技术(厦门)有限公司 | 一种时间序列的异常检测方法、装置、设备及存储介质 |
CN113792032A (zh) * | 2021-08-09 | 2021-12-14 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于正态分布误差修正的测风数据塔影效应分析方法 |
CN113792032B (zh) * | 2021-08-09 | 2024-01-23 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于正态分布误差修正的测风数据塔影效应分析方法 |
CN113991638A (zh) * | 2021-08-31 | 2022-01-28 | 华能(福建)能源开发有限公司 | 一种针对不同地点的新能源场站发电功率的预测方法 |
CN113809772A (zh) * | 2021-09-17 | 2021-12-17 | 国网河南省电力公司电力科学研究院 | 提升应对次时间尺度风电不确定性的安全性的方法及装置 |
CN113809772B (zh) * | 2021-09-17 | 2023-09-08 | 国网河南省电力公司电力科学研究院 | 提升应对次时间尺度风电不确定性的安全性的方法及装置 |
CN114897245A (zh) * | 2022-05-12 | 2022-08-12 | 国网湖北省电力有限公司电力科学研究院 | 一种基于长短期记忆神经网络的规模风电超短期功率预测误差修正方法 |
CN114897245B (zh) * | 2022-05-12 | 2024-05-28 | 国网湖北省电力有限公司电力科学研究院 | 一种基于长短期记忆神经网络的规模风电超短期功率预测误差修正方法 |
CN116448263B (zh) * | 2023-06-16 | 2023-09-05 | 山东德圣源新材料有限公司 | 一种勃姆石生产设备运行状态检测方法 |
CN116448263A (zh) * | 2023-06-16 | 2023-07-18 | 山东德圣源新材料有限公司 | 一种勃姆石生产设备运行状态检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107230977A (zh) | 基于误差修正和提升小波组合预测模型的风电功率预测方法 | |
Catalao et al. | An artificial neural network approach for short-term wind power forecasting in Portugal | |
Singh et al. | An efficient time series forecasting model based on fuzzy time series | |
CN102102626B (zh) | 风电场短期功率预测方法 | |
CN108256697A (zh) | 一种用于电力系统短期负荷的预测方法 | |
Celik et al. | Generalized feed-forward based method for wind energy prediction | |
Panapongpakorn et al. | Short-term load forecast for energy management systems using time series analysis and neural network method with average true range | |
CN109242143A (zh) | 一种神经网络风电功率预测方法及系统 | |
CN111861039A (zh) | 基于lstm和广义预测控制算法的电力负荷预测方法、系统、设备及存储介质 | |
CN105447594A (zh) | 一种基于指数平滑的电力系统灰色负荷预测方法 | |
CN109242212A (zh) | 一种基于变模式分解与长短记忆网络的风电预测方法 | |
CN105117593A (zh) | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 | |
CN103559563A (zh) | 一种风电场短期风速预测方法 | |
CN104036328A (zh) | 自适应风电功率预测系统及预测方法 | |
CN104933489A (zh) | 基于自适应神经模糊推理系统的风电功率实时高精度预测方法 | |
Chen et al. | Air-conditioning load forecasting for prosumer based on meta ensemble learning | |
Rizwan et al. | Artificial intelligence based approach for short term load forecasting for selected feeders at madina saudi arabia | |
CN104463356A (zh) | 一种基于多维信息人工神经网络算法的光伏发电功率预测方法 | |
CN113869795B (zh) | 一种工业副产煤气系统长期调度方法 | |
CN110135634A (zh) | 中长期电力负荷预测装置 | |
Nagaraja et al. | A survey on wind energy, load and price forecasting:(Forecasting methods) | |
CN111489038B (zh) | 一种提升风电场超短期功率预测准确率的方法 | |
CN110222910A (zh) | 一种有源配电网态势预测方法及预测系统 | |
Lin et al. | Hybrid RNN-LSTM deep learning model applied to a fuzzy based wind turbine data uncertainty quantization method | |
Gorbatiuk et al. | Application of fuzzy time series forecasting approach for predicting an enterprise net income level |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171003 |