CN105868853A - 一种短期风电功率组合概率预测方法 - Google Patents
一种短期风电功率组合概率预测方法 Download PDFInfo
- Publication number
- CN105868853A CN105868853A CN201610182801.3A CN201610182801A CN105868853A CN 105868853 A CN105868853 A CN 105868853A CN 201610182801 A CN201610182801 A CN 201610182801A CN 105868853 A CN105868853 A CN 105868853A
- Authority
- CN
- China
- Prior art keywords
- model
- prediction
- wind power
- wind
- period
- 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 81
- 238000009826 distribution Methods 0.000 claims abstract description 63
- 238000004836 empirical method Methods 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 97
- 238000012549 training Methods 0.000 claims description 27
- 238000007476 Maximum Likelihood Methods 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 9
- 238000005315 distribution function Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 2
- 230000001186 cumulative effect Effects 0.000 claims 1
- 238000002474 experimental method Methods 0.000 claims 1
- 238000012417 linear regression Methods 0.000 abstract description 5
- 238000012706 support-vector machine Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 239000013598 vector Substances 0.000 description 8
- 238000012360 testing method Methods 0.000 description 6
- 238000013277 forecasting method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000007726 management method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004899 motility Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- 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)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Power Engineering (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种短期风电功率组合概率预测方法,包括:构建每个前瞻时段的各单项概率预测模型。将支持向量机和线性回归预测得到的点预测结果假设为符合weibull分布和gauss分布的概率结果,应用经验方法统计得到所假设各分布的形状参数。构建每个前瞻时段的扩展BMA模型。将风电功率预测所需要的数据输入所述各单项概率预测模型,得到各单项预测风电功率条件概率密度函数,利用扩展BMA模型将其组合预测得到风电功率组合概率密度函数。本模型可以根据不同风场数据的特点自适应地得到任意形式的分布形式,可有效降低模型分布形式的选择带来的模型误差,真正做到预测模型的普适性。
Description
技术领域
本发明涉及新能源发电过程中风电功率预测技术领域,特别是涉及一种短期风电功率组合概率预测方法。
背景技术
由于风电随机性和波动性的存在,大规模风电并网对电力系统运行造成很大影响,如系统的稳定性和电能质量,因此风力发电被认为是不可调度的。风功率不确定性预测是解决这些问题、提高电力系统风电并网能力的一种有效工具。
最优决策与预测分布的分位数直接相关结论的证明,为风力发电的不确定性预测被应用到电力系统运行提供了理论基础。在随机优化的帮助下,风电出力的不确定性信息已被用于满足相关决策问题的要求、风电机组的交易策略以及考虑到风力发电的不确定性后的能源存储大小和风电站优化调度。这些研究表明,不确定性预测的应用使风力发电的利用率大大增加。概率预测是最常用的不确定性表示。不同的决策情况需要不同的分位数或者预测区间,可以通过预测完整的概率密度函数(PDF)来满足需求的灵活性。
典型的风电功率概率密度预测方法多为单项预测方法,通常根据数据从多种假设模型中选择一种作为最优预测模型。然而,任何一种单项预测方法均有其固有的局限性,且只适用于部分风场,目前还没有一种概率预测方法可以适用于所有风场。因此,所选择的单项预测模型对现有数据来说不一定是最优模型。其他合理的模型也可以针对样本数据给出不同的预测结果,这些合理的模型就是预测结果不确定性的来源。然而,应用以上的典型方法,认为单一模型是“最好”的,忽略了这种预测不确定性的来源,从而低估了不确定性。而组合预测方法能够综合多种相同预测水平的单项方法的优势,得到不同情况下各种单项预测方法优势互补的预测方法,从而降低单个预测模型中随机因素的影响,以提高预测精度。
现有方法或者将多个单点预测值进行加权平均,再对其组合误差进行统计得到误差的PDF,此方法主要是对单点预测值进行组合;或者利用组合方法预测得到风速风向,再利用风机功率曲线得到风电功率的PDF。现实中,这一转换误差是相对较大的。现有的组合概率预测模型均假设了单项模型符合相同分布的参数预测,这与现实中的实际问题有所偏差。在现实问题中,预测目标量的分布可能随着时空的变化趋向于不同的分布,有些甚至不符合某种分布。
发明内容
本发明的目的就是为了解决上述难题,提供了一种短期风电功率组合概率预测方法,该方法能够解决现有风电功率预测技术中存在的或者选择单一模型作为最优预测模型、或者组合多种确定性预测模型、或者组合同分布单项预测模型等问题。
为实现上述目的,本发明的具体方案如下:
一种短期风电功率组合概率预测方法,包括:
构建每个前瞻时段的SBL(Sparse Bayesian Learning,稀疏贝叶斯学习)单项预测模型和KDE(Kernel Density Estimation,核密度估计)单项预测模型;
构建每个前瞻时段的单项gauss分布预测模型,构建SVM(Support VectorMachine,支持向量机)预测模型,利用经验法统计得到gauss分布函数的方差值;
构建每个前瞻时段的单项weibull分布预测模型,构建LR(Linear Regression,线性回归)预测模型,利用经验法统计得到weibull分布函数的形状参数k;
利用训练数据样本,将风电功率预测所需要的数据输入所述SBL预测模型,得到每个前瞻时段的单项风电功率预测条件概率密度函数正态分布函数p(y|M1);
利用训练数据样本,将风电功率预测所需要的数据输入所述KDE预测模型,得到每个前瞻时段的单项风电功率预测条件概率密度函数KDE函数p(y|M2);
利用训练数据样本,将风电功率预测所需要的数据输入所述SVM预测模型,得到每个前瞻时段的风电功率预测值E3,结合所述经验统计得到的gauss函数的方差值,得到每个前瞻时段的单项风电功率预测条件概率密度函数gauss函数p(y|M3);
利用训练数据样本,将风电功率预测所需要的数据输入所述LR预测模型,得到每个前瞻时段的风电功率预测值E4,结合经验统计法得到的形状参数k对所述比例参数λ进行修正,得到修正后的比例参数λ,从而得到每个前瞻时段的单项风电功率预测条件概率密度函数weibull函数p(y|M4);
构建扩展BMA模型,将上述SBL单项预测模型、KDE单项模型、SVM预测模型以及LR预测模型得到的概率密度函数进行组合,得到未来时段的风电功率值的目标量y的概率密度函数;
利用训练样本数据训练扩展BMA模型得到模型权值参数,以概率精度指标值最小为目标函数,对模型权值参数进行优化,得到最终的风电功率组合概率预测结果。
将所述4种单项预测方法概率预测模型得到的条件概率密度函数p(y|M1)、p(y|M2)、p(y|M3)和p(y|M4)通过加权平均得到最终的组合概率密度函数,求解各权值参数;
将风电功率的分布预测所需要的数据输入所述扩展BMA模型,将风电功率的分布预测所需要的数据输入各单项概率预测模型,得到每个单项模型在前瞻时段的预测条件概率密度函数,通过将所述权值参数和各单项预测条件概率密度函数代入公式预测得到最终的风电功率组合概率密度函数。
优选地,构建每个前瞻时段的SBL单项预测模型,包括:
利用风速历史数据和对应的风电场输出功率的历史数据作为样本对SBL进行训练,得到SBL模型权重系数,形成每个前瞻时段的SBL预测模型。
优选地,构建每个前瞻时段的KDE单项预测模型,包括:
利用风速历史数据和对应的风电场输出功率的历史数据作为样本对KDE进行训练,得到KDE模型带宽参数,形成每个前瞻时段的KDE预测模型。
优选地,构建每个前瞻时段的gauss分布单项预测模型,包括:
构建每个前瞻时段的SVM预测模型,利用风速历史数据和对应的风电场输出功率的历史数据作为样本对SVM进行训练,得到SVM模型权重系数,形成每个前瞻时段的SVM预测模型。
利用经验法对风电功率历史数据进行统计,得到gauss分布函数的方差值。
优选地,构建每个前瞻时段的weibull分布单项预测模型,包括:
构建每个前瞻时段的LR预测模型,利用风速历史数据和对应的风电场输出功率的历史数据作为样本对LR进行训练,得到LR模型权重系数,形成每个前瞻时段的LR预测模型。
利用经验法对风电功率历史数据进行统计,得到weibull分布的形状参数k。
优选地,构建每个前瞻时段的扩展BMA预测模型,包括:
利用所述4种单项预测模型,以测试样本的目标时段风速数据和距离所述目标时段最近三个时段的风电场输出功率数据作为输入,得到所述训练样本的每个前瞻时段的风电功率预测条件概率密度函数p(y|M1)、p(y|M2)、p(y|M3)和p(y|M4);
将所述4种单项模型概率预测结果进行加权平均得到组合概率密度预测结果,组合公式为:
优选地,采用极大似然估计方法初步模型的权值参数,包括:
利用用极大似然估计初步求解(w1,w2,…,wk)。似然函数为其中,θ为所有的参数集合。
利用求解最大似然函数的常用方法EM方法求解所述似然函数。
步骤一(E step):
步骤二(M step):
优选地,以24小时的概率指标连续排名概率得分(CRPS)平均值最小为目标函数,采用PSO方法进一步优化所述权值参数,包括:
24小时前瞻时段内的CRPS平均值定义为其中T=24,且
目标函数为min CRPS,约束条件为wk-α≤wk≤wk+α和由于目标函数解析式较复杂,难以解析的求解最优参数。采用粒子群算法对参数进行寻优,wk参数范围定为以极大似然估计结果为中心的小范围区间,α此处取值为0.05。
优选地,将风电功率的分布预测所需要的数据输入所述各单项概率预测模型,包括:
将所述预测目标时段风速风向数据、所述风电功率预测值和距离所述预测目标时段最近三个时段的误差数据输入所述各单项概率预测模型。
优选地,在第4种单项概率预测模型中,利用所述LR预测均值和所述经验统计法得到的形状参数计算得到weibull分布预测的另一重要参数比例参数,包括:
利用所述形状参数k和所述LR预测模型得到的均值E4,根据公式对比例参数λ进行计算。
优选地,将所述数值风速数据的横纵坐标形式转换为风速向量形式,即风速幅值和风向角度,进一步对其进行归一化,且采用五点三次平滑滤波方法对所有样本数据进行滤波处理。
由以上技术方案可见,本发明实施例公开了一种短期风电功率组合概率预测方法,包括,在贝叶斯模型平均(BMA)方法的基础上,将BMA所组合的单项模型采用不同形式的概率函数表示,包括:gauss、weibull函数以及核密度估计函数,并且以24小时前瞻时段内的CRPS平均值作为目标函数,对模型参数进行进一步优化。本文选用的稀疏贝叶斯(SBL)模型和核密度估计模型,均以目标时段的风速风向数据作为输入,对目标时段的风电功率概率密度函数直接进行预测。此外,为了达到组合预测的效果,将支持向量机和线性回归预测得到的点预测结果分别假设为符合为gauss分布和weibull分布的概率结果,应用经验方法统计得到所假设各分布的形状参数。点预测值结合经验方法得到的形状参数,得到最终的概率密度预测结果。综上,本模型所用单项模型概率密度函数符合不同分布,包括SBL模型得到的Gaussian函数、核密度估计的概率结果、经验法扩展得到的gauss函数和weibull函数四种概率函数。本模型可以根据不同风场数据的特点自适应地得到任意形式的分布形式,可有效降低模型分布形式的选择带来的模型误差,真正做到预测模型的普适性。
附图说明
图1为本发明实施例提供的一种短期风电功率组合概率密度预测方法的流程示意图;
图2为本发明实施例提供的构建每个前瞻时段的gauss分布预测模型的步骤图;
图3为本发明实施例提供的构建每个前瞻时段的weibull分布预测模型的步骤图;
图4为本发明实施例提供的前瞻24h风电功率概率密度曲线与真实值展示图;
图5(a)-图5(f)为本发明实施例提供的1:00,5:00,9:00,13:00,17:00,21:00的风电功率预测概率密度与真实值展示图;
图6为本发明实施例提供的前瞻24小时风电功率预测区间与真实值对比图。
具体实施方式:
下面结合附图对本发明进行详细说明:
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
SVM是在VC维理论和经验风险最小化原则的基础上提出的一种新型学习机器,其最大的特点是采用少数支持向量代表整个样本集,改变了传统的经验风险最小化原则,具有很好的泛化能力,另外,SVM在处理非线性问题,通过将非线性问题映射到高维空间转化为线性问题,从而巧妙的解决了反复计算问题,有效的克服了维数灾难以及局部极小问题。SVM回归的思想是通过一个非线性映射将数据映射到高维特征空间H,并在这个空间进行线性回归,具体为,给定k个样本数据,其值为{xk,yk},其中xk∈Rn为n维向量,yk∈R为相应的输出变量,SVM通过引入损失函数来解决回归问题,一般采用 ω∈Rn来估计输出值,式中,ω为权重向量,b为偏置项,可以通过解凸最优化问题得到ω和b。通常用核函数K替代内积计算,即式中ai和为不等于0的拉格朗日乘子。
SBL理论,由Tipping于2000年在SVM的基础上提出,与SVM采用结构风险最小化原则不同,SBL采用了贝叶斯推理方法,这一变化不仅使SBL具有更好的稀疏特性,而且使其具有概率预测的能力。稀疏贝叶斯学习机预测模型可表示为:式中:youtput为待预测随机变量;xinput为输入向量;xi为训练样本中的输入向量;K(·)为核函数,在回归预测中可采用高斯核函数形式;M为训练样本总数;ε为误差项,服从正态分布wi(w0同)为权重系数,在稀疏贝叶斯学习机中被看作随机变量,并假设其先验分布为容易看出,当式所述学习机训练完成时,对于任意给定的输入向量,均可得到被预测量的概率密度函数。学习机的训练过程,则是根据贝叶斯原则对参数w=[w0,w1,…,wM]T、超参α=[α0,α1,…,αM]T以及σ2的后验分布进行推断的过程,即可表示为p(w,α,σ2|y),其中,y=[y1,y2,…,yM]T。然而,由于对全部参数进行贝叶斯推断计算复杂,在实际计算过程中,超参α与σ2是通过对联合分布的极大似然估计得到的,表示为αMP与进而,可以利用贝叶斯推断得到w的后验分布为其中:
Φ=[φ(x1),φ(x2),…,φ(xM)]T;
φ(xi)=[1,K(xi,x1),…,K(xi,xM)]T;
Α=diag(α0,α1,…αM)。
在得到w及ε的后验分布后,带入式(1)即完成了稀疏贝叶斯学习机的训练过程。而实际上,由于youtput仍然服从正态分布,可直接写出其均值与方差为:与
KDE是在概率论中用来估计未知的密度函数,属于非参数检验方法之一,由Rosenblatt(1955)和Emanuel Parzen(1962)提出。对于数据x1,x2,...,xn,核密度估计的形式为:这是一个加权平均,是一个加权平均,而核函数K(·)是一个权函数,核函数的形状和值域控制着用来估计f(x)在点x的值时所用数据点的个数和利用的程度,直观来看,核密度估计的好坏依赖于核函数和带宽h的选取。
LR是利用数理统计中回归分析,来确定两种或两种以上变量间相互依赖的定量关系的一种统计分析方法,运用十分广泛。其表达形式为y=w'x+e,e为误差服从均值为0的正态分布。
经验统计方法是根据用于评估分布与数据的拟合度,估计百分位数并比较不同的样本分布。
BMA是一种组合来自不同来源的预测分布结果的方法。BMA预测得到的组合概率密度函数(PDF)是以多个单项预测概率密度函数为核心的加权平均函数,其以后验概率作为权重对所有可能的单项概率预测模型进行加权平均(后验概率表示单项预测模型真实描述实际过程的概率)。这些权重可反映每个单项预测模型对组合预测结果的相对贡献,也可以在有大量单项预测模型的情况下作为选择较优单项模型的基础。基于训练数据yT,根据BMA模型,由K个概率模型M1、M2、…、MK组合得到的目标量y的概率密度函数可描述为其中,p(Mk|yT)为有训练数据yT得到的模型Mk为正确的拟合模型的后验概率,可作为模型Mk在组合模型中的权值系数,用来反映模型Mk对训练数据的拟合程度。p(y|Mk)是由第k个模型Mk单独预测得到的目标值y的条件概率密度函数。原始BMA模型为静态预测模型,2005年Rafter等人将其扩展为动态预测模型,主要用于气象数据(如风速、温度、海平面气压等等)预测,每个气象站作为一个单项预测模型,综合多个气象站对目标量的点预测结果得到更准确的未来时段的风电功率值的目标量y以及目标量的概率密度函数,为:其中,fk表示第k个单项模型对目标量的预测值;gk(y|fk)表示基于预测值fk的目标量yn的条件概率密度函数。所有单项模型根据预测目标量的特点选取概率密度函数的形式,但是对于同一目标量采用相同的分布形式。如用于温度和海平面气压预测时,假设单项概率密度函数符合均值为线性函数ak+bkyk、标准差为σ的高斯分布。其中yk为第k个气象站预测得到的目标量点预测值。线性函数系数由各气象站预测历史数据以及历史真实值应用线性回归得到,标准差σ应用EM方法得到。
本实施例将基于上述SBL、KDE、SVM、LR、经验统计法以及BMA方法,对未来1-48h,的风电功率概率密度函数进行预测,其中,时间分辨率为1h。本实施例由构建单项模型预测模型、构建组合预测模型以及风电功率组合概率密度函数预测3部分组成,并通过分步预测对多时段的风电功率进行预测。在具体预测过程,首先,通过SBL及KDE单项预测模型对风电功率概率密度函数进行预测;其次,通过SVM预测模型对假设符合gauss分布的风电功率均值进行预测,结合经验统计法得到的gauss分布的方差;进而,通过LR预测模型对假设符合Weibull分布的风电功率均值进行预测,结合经验统计法得到的weibull分布的形状参数k,得到weibull分布的比例参数最后,应用求得权值参数的扩展BMA模型将各单项模型得到的概率密度函数进行组合得到最后的风电功率组合概率密度函数。并且本实施例中以全球能源预测比赛(GEFCOM)的十个不同风电场数据为例对方法进行验证,证明本实施例方法的有效性。
参见图1,为本发明实施例提供的一种短期风电功率组合概率预测方法的流程示意图,该方法包括:
步骤101:构建每个前瞻时段的各单项概率预测模型。
所用样本数据包含风速、风向信息,在构建每个前瞻时段的单项概率预测模型时,两者均为风电场输出功率的主要影响因素,因此本实施例中,采用风速、风向历史数据和对应的风电场输出功率的历史数据作为样本对各单项概率预测模型进行训练,形成每个前瞻时段的单项概率预测模型。
样本数据的处理:
所述风速历史数据包括风电场的风速相量横纵坐标,此处为了模型的应用方便将其转换为风速向量形式,即风速的幅值和风向角度,经验证表明将数据归一化有利于数据规律的统计,因此对风速幅值和风向角度进行归一化处理。
在样本数据的测量中,会不可避免的引入错误点、冗余点以及扫描环境所带来的测量噪声等,这些点对统计模型的建立会带来很大影响,为了更好的提取数据样本中的特征数据,进行数据滤波,将这些错误点等进行去除。此处采用五点三次滤波方法对风速、风向以及对应的风电功率样本数据进行滤波处理,以提高所建立统计模型的精确度。
在构建每个前瞻时段的SBL和KDE预测模型时,利用所述SBL和KDE预测模型,以测试样本的目标时段风速、风向数据和距离所述目标时段最近三个时段的风电场输出功率数据作为输入,得到所述测试样本的每个前瞻时段的风电功率预测概率密度函数。
在构建每个前瞻时段的gauss分布预测模型时,具体包括如下步骤:
步骤201:利用所述SVM预测模型,以测试样本的目标时段风速、风向数据和距离所述目标时段最近三个时段的风电场输出功率数据作为输入,得到所述测试样本的每个前瞻时段的风电功率预测目标值,即所假设gauss分布的均值。
步骤202:利用所述经验统计方法,以历史风电功率数据为样本,得到gauss分布的方差。
在构建每个前瞻时段的weibull分布预测模型时,具体包括如下步骤:
步骤301:利用所述LR预测模型,以测试样本的目标时段风速、风向数据和距离所述目标时段最近三个时段的风电场输出功率数据作为输入,得到所述测试样本的每个前瞻时段的风电功率预测目标值,即所假设weibull分布的均值E4。
步骤302:利用所述经验统计方法,以历史风电功率数据为样本,得到weibull分布的形状参数k。
步骤102:构建每个前瞻时段的扩展BMA模型。
针对风电功率数据进行预测,并无传统BMA模型所需的各气象站预测模型所得到的目标量预测数据。因此,此处采用其他不同的单项模型对风电功率概率密度函数进行预测,且所用单项模型概率密度函数符合不同分布,形成进一步扩展的BMA方法对符合不同分布的单项模型进行组合。
在现有的风电功率预测方法中,直接对目标时段的风电功率概率密度函数预测的并不多。因此,此方法中单项预测模型的选取是一个难点。此处选用近年来研究较多、预测效果较好且计算相对较简单的稀疏贝叶斯(SBL)模型和核密度估计模型。这两种方法均以目标时段的风速风向数据作为输入,对目标时段的风电功率概率密度函数直接进行预测。此外为了达到组合预测的效果,本文借鉴raftery的方法,将支持向量机和线性回归预测得到的点预测结果分别假设为符合为gauss分布和weibull分布的概率结果,应用经验方法统计得到所假设各分布的形状参数。针对weibull分布,结合点预测值与经验方法得到的形状参数,计算得到weibull分布的另一重要参数比例参数。
利用训练样本构建扩展BMA模型,以所述训练样本中的目标时段风速、风向数据和距离所述目标时段最近三个时段的风电场输出功率数据作为输入,预测得到所述4种单项概率模型在每个前瞻时段的风电功率条件概率密度函数数p(y|M1)、p(y|M2)、p(y|M3)和p(y|M4)。将所述单项概率预测结果代入组合公式此处K=4。
应用极大似然估计初步求解(w1,w2,…,wK)。似然函数为:
其中,θ为所有的参数集合。
利用EM方法对似然函数进行求解。
步骤一(E step):
步骤二(M step):
以24小时前瞻时段内的CRPS平均值最小为目标函数,进一步优化模型参数。
CRPS定义:
其中,T=24。
目标函数为:min CRPS。
约束条件:
wk-α≤wk≤wk+α (6)
目标函数解析式较复杂,难以解析的求解最优参数。采用粒子群算法对参数进行寻优,参数范围定为以极大似然估计结果为中心的小范围区间,α此处取值为0.05。
步骤103:将风电功率预测所需要的数据输入所述各单项概率预测模型,得到各单项预测风电功率条件概率密度函数,利用扩展BMA模型将其组合预测得到风电功率组合概率密度函数。
将风电功率预测所需要的数据输入所述SBL和KDE预测模型,得到这两个单项预测模型每个前瞻时段的风电功率概率密度函数。
将风电功率预测所需要的数据输入所述SVM预测模型,得到gauss分布单项预测模型每个前瞻时段的风电功率均值,结合所述经验分布发得到的方差值,得到符合gauss分布的风电功率概率密度函数。
将风电功率预测所需要的数据输入所述SVM预测模型,得到weibull分布单项预测模型每个前瞻时段的风电功率均值E4,结合所述经验分布发得到的形状参数,利用公式得到weibull分布的另一重要参数比例参数λ,结合所述weibull分布的比例参数和形状参数,得到符合weibull分布的风电功率概率密度函数。
将4种单项模型预测得到的风电功率条件概率密度函数代入扩展BMA模型公式得到目标时段的风电功率概率密度函数。
利用上述提供的组合预测方法,本实施例利用全球能源预测比赛的十个风场数据,包括风速、风向和风电功率历史数据,模拟进行前瞻T(T=24)小时的风电功率预测,进行有效性验证。
首先将样本分为训练集和验证集两部分,分别含有1000和500个样本,然后,利用训练集数据对每个前瞻时段训练构建4种单项概率预测模型,利用每个单项概率预测模型对训练集样本进行模拟预测,得到4种单项风电功率概率密度预测结果;其后,将训练集产生的每个前瞻时段的4种单项风电功率概率密度预测结果应用扩展BMA模型进行组合,先用EM方法对模型参数初步求解,再以24小时前瞻时段内的CRPS平均值最小为目标函数,对模型参数进一步优化,构建每个前瞻时段的扩展BMA模型;最后,对所述验证集中的数据进行预测,最终得到关于所述验证集样本的每个前瞻时段的风电功率组合概率密度函数。
下面是关于确定性预测精度和概率预测结果的预测效果分析:
首先,NMAE是最基本的风电功率预测模型的确定性结果评价指标,如下式所示:
N表示数据总数;表示i时刻的实际测量数据;表示i时刻的预测数据;YN表示风场的额定装机容量。
10个风场24小时的NMAE指标值(%)如表一所示。由表可知,确定性预测效果随着预测时段的增加而变化,且10个风场的确定性预测效果相似,表明了本模型的普适性。多方法求解的10个风场的24小时平均NMAE指标值(%)对比如表二所示。由表可知,相较传统BMA模型和所组合的每种单项预测模型,本实施方法的确定性预测效果更好。
表一10个风场24小时的NMAE指标值(%)
表二:多方法求解的10个风场的24小时平均NMAE指标值(%)对比
其次,直观的展示了本发明的方法的预测概率密度曲线与预测区间,说明本发明方法在概率预测效果方面的合理性。
为直观说明预测效果,图4直观展示了前瞻24小时(2012年2月17日)的预测风电功率概率密度曲线以及实测值和预测均值。从图中可看出,预测得到的风电功率概率密度曲线是非参数且时变的。图5(a)-图5(f)分别展示了其中1:00、5:00、9:00、13:00、17:00、21:00的预测概率密度曲线和置信水平为80%的预测区间,其中,圆圈表示风电功率实际值,星号为80%置信区间的上下限值。由图可看出,实测值均落在预测结果的80%置信区间内,并且可以直观看出组合概率密度曲线综合了所有单项预测结果的锐度和点预测准确度优势,概率密度曲线随着时间变化做出相应的调整,以达到最优。图6展示了前瞻24小时的10%-90%预测置信区间,可以看到,预测区间随着时间的变化而变化,且真实值大多数落到了80%置信区间内,表明本方法有较高的区间预测准确度。
通过上述采用确定性精度以及对预测概率密度函数和预测区间的展示,表明本实施例方法可通过组合符合不同分布的单项模型概率密度函数得到有效的风电功率概率密度函数,且可以根据不同风场数据的特点自适应地得到任意形式的分布形式,可有效降低模型分布形式的选择带来的模型误差,真正做到预测模型的普适性。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种短期风电功率组合概率预测方法,其特征在于,包括:
(1)利用训练样本数据,分别构建每个前瞻时段的SBL单项预测模型和KDE单项模型;构建每个前瞻时段的gauss分布模型和SVM预测模型;构建每个前瞻时段的weibull分布模型以及每个前瞻时段的LR预测模型;
(2)将风电功率预测所需要的数据分别输入所述SBL单项预测模型和KDE单项模型,得到两种单项预测模型每个前瞻时段的预测风电功率概率密度函数,分别为gauss函数和KDE函数;
将风电功率预测所需要的数据输入所述SVM预测模型,得到每个前瞻时段的风电功率预测均值以及预测条件概率密度函数gauss函数;
将风电功率预测所需要的数据输入所述LR预测模型,得到每个前瞻时段的风电功率预测均值以及预测条件概率密度函数weibull函数;
(3)构建扩展BMA模型,将上述SBL单项预测模型、KDE单项模型、SVM预测模型以及LR预测模型得到的概率密度函数进行组合,得到未来时段的风电功率值的目标量y的概率密度函数;
利用训练样本数据训练扩展BMA模型得到模型权值参数,以概率精度指标值最小为目标函数,对模型权值参数进行优化,得到最终的风电功率组合概率预测结果。
2.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,构建每个前瞻时段的SBL单项预测模型,包括:
利用风速历史数据和对应的风电场输出功率的历史数据作为样本对SBL进行训练,得到SBL模型权重系数,形成每个前瞻时段的SBL单项预测模型。
3.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,构建每个前瞻时段的KDE单项预测模型,包括:
利用风速历史数据和对应的风电场输出功率的历史数据作为样本对KDE模型进行训练,得到KDE模型带宽参数,形成每个前瞻时段的KDE单项预测模型。
4.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,构建每个前瞻时段的单项gauss分布预测模型,包括:
构建每个前瞻时段的单项SVM预测模型,利用风速历史数据和对应的风电场输出功率的历史数据作为样本对SVM进行训练,得到SVM模型权重系数,形成每个前瞻时段的SVM预测模型。
5.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,构建每个前瞻时段的单项weibull分布预测模型,包括:
构建每个前瞻时段的单项LR预测模型,利用风速历史数据和对应的风电场输出功率的历史数据作为样本对LR模型进行训练,得到LR模型参数,形成每个前瞻时段的LR预测模型。
6.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,所述步骤(2)中,将风电功率预测所需要的数据输入所述SVM预测模型,得到每个前瞻时段的风电功率均值预测结果E3,利用经验方法统计训练数据得到所构建gauss分布函数的方差参数得到SVM预测模型预测时段的概率密度预测结果gauss概率密度函数。
7.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,所述步骤(2)中,将风电功率预测所需要的数据输入所述LR预测模型,得到每个前瞻时段的风电功率均值预测结果E4,利用经验方法统计训练数据得到所构建weibull分布函数的形状参数k,根据预测结果E4和形状参数k计算weibull分布的比例参数λ;根据所述形状参数k和比例参数λ,得到LR预测模型在预测时段的概率密度预测结果weibull函数。
8.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,所述步骤(3)中,
由SBL单项预测模型M1、KDE单项模型M2、SVM预测模型M3和LR预测模型M4组合得到的未来时段的风电功率值的目标量y的概率密度函数具体为:
其中,p(Mk|yT)为由训练样本数据yT得到的模型Mk为正确的拟合模型的后验概率,作为模型Mk在组合模型中的权值系数wk,用来反映模型Mk对训练数据的拟合程度p(y|Mk)是由第k个模型Mk单独预测得到的目标值y的条件概率密度函数。
9.根据权利要求1所述的一种短期风电功率组合概率预测方法,其特征在于,所述步骤(3)中,利用极大似然估计方法对扩展BMA模型权值参数进行求解,利用粒子群法对扩展BMA模型权值参数进行优化。
10.根据权利要求9所述的一种短期风电功率组合概率预测方法,其特征在于,以概率评价指标CRPS值最小作为目标函数,利用粒子群法对模型权值参数进行优化:
目标函数为
其中,
Ft,n(y)为前瞻t时段的第n次预测给出的风电功率累积分布函数;yt,n为前瞻t时段的第n次实验的风电功率测量值;p代表以装机容量为基准值标幺化的风电功率随机变量;N表示验证集样本每个时段内的预测样本数。
约束条件为
其中,权值参数wk参数范围定为以极大似然估计结果为中心的小范围区间,α为常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610182801.3A CN105868853B (zh) | 2016-03-28 | 2016-03-28 | 一种短期风电功率组合概率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610182801.3A CN105868853B (zh) | 2016-03-28 | 2016-03-28 | 一种短期风电功率组合概率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105868853A true CN105868853A (zh) | 2016-08-17 |
CN105868853B CN105868853B (zh) | 2017-05-17 |
Family
ID=56624944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610182801.3A Active CN105868853B (zh) | 2016-03-28 | 2016-03-28 | 一种短期风电功率组合概率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105868853B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106530123A (zh) * | 2016-10-26 | 2017-03-22 | 广东电网有限责任公司电力科学研究院 | 一种发电负荷稳态点提取方法及其装置 |
WO2018082132A1 (zh) * | 2016-11-04 | 2018-05-11 | 江南大学 | 一种用于能源系统输出功率的短期预测方法 |
CN108022019A (zh) * | 2017-12-14 | 2018-05-11 | 西南石油大学 | 基于风电机组分类的风电场调度方法和系统 |
CN108599147A (zh) * | 2018-04-13 | 2018-09-28 | 华北电力大学 | 基于正态指数平滑法与核密度估计的组合区间预测方法 |
CN109800919A (zh) * | 2019-01-30 | 2019-05-24 | 云南电网有限责任公司电力科学研究院 | 一种区域风电场整体有功功率预测方法及装置 |
CN109840308A (zh) * | 2017-11-29 | 2019-06-04 | 中国电力科学研究院有限公司 | 一种区域风电功率概率预报方法及系统 |
CN109891235A (zh) * | 2016-10-28 | 2019-06-14 | 恩德斯+豪斯流量技术股份有限公司 | 基于统计方法,自动在线检测流体的实际状态与流体的参考状态的偏差的方法,特别是用于监测饮用水供应 |
CN109978222A (zh) * | 2019-01-04 | 2019-07-05 | 国网山东省电力公司电力科学研究院 | 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 |
CN110858715A (zh) * | 2018-08-25 | 2020-03-03 | 南京理工大学 | 基于风电随机性和相关性的风电穿透功率极限分析方法 |
CN111815039A (zh) * | 2020-06-28 | 2020-10-23 | 山东大学 | 基于天气分类的周尺度风电功率概率预测方法及系统 |
CN112269967A (zh) * | 2020-10-14 | 2021-01-26 | 清华大学 | 一种基于联合机会约束的迭代拆分方法及系统 |
CN113487105A (zh) * | 2021-07-27 | 2021-10-08 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440541A (zh) * | 2013-09-18 | 2013-12-11 | 山东大学 | 多风电场短期输出功率的联合概率密度预测方法 |
CN103903067A (zh) * | 2014-04-09 | 2014-07-02 | 上海电机学院 | 风电功率短期组合预测方法 |
CN105225006A (zh) * | 2015-09-30 | 2016-01-06 | 杨明 | 一种短期风电功率非参数概率预测方法 |
-
2016
- 2016-03-28 CN CN201610182801.3A patent/CN105868853B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440541A (zh) * | 2013-09-18 | 2013-12-11 | 山东大学 | 多风电场短期输出功率的联合概率密度预测方法 |
CN103903067A (zh) * | 2014-04-09 | 2014-07-02 | 上海电机学院 | 风电功率短期组合预测方法 |
CN105225006A (zh) * | 2015-09-30 | 2016-01-06 | 杨明 | 一种短期风电功率非参数概率预测方法 |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106530123B (zh) * | 2016-10-26 | 2020-03-27 | 广东电网有限责任公司电力科学研究院 | 一种发电负荷稳态点提取方法及其装置 |
CN106530123A (zh) * | 2016-10-26 | 2017-03-22 | 广东电网有限责任公司电力科学研究院 | 一种发电负荷稳态点提取方法及其装置 |
CN109891235A (zh) * | 2016-10-28 | 2019-06-14 | 恩德斯+豪斯流量技术股份有限公司 | 基于统计方法,自动在线检测流体的实际状态与流体的参考状态的偏差的方法,特别是用于监测饮用水供应 |
CN109891235B (zh) * | 2016-10-28 | 2021-12-31 | 恩德斯+豪斯流量技术股份有限公司 | 基于统计方法,自动在线检测流体的实际状态与流体的参考状态的偏差的方法,特别是用于监测饮用水供应 |
US11193920B2 (en) | 2016-10-28 | 2021-12-07 | Endress+Hauser Flowtec Ag | Method for the automated in-line detection of deviations of an actual state of a fluid from a reference state of the fluid on the basis of statistical methods, in particular for monitoring a drinking water supply |
WO2018082132A1 (zh) * | 2016-11-04 | 2018-05-11 | 江南大学 | 一种用于能源系统输出功率的短期预测方法 |
CN109840308A (zh) * | 2017-11-29 | 2019-06-04 | 中国电力科学研究院有限公司 | 一种区域风电功率概率预报方法及系统 |
CN109840308B (zh) * | 2017-11-29 | 2023-12-19 | 中国电力科学研究院有限公司 | 一种区域风电功率概率预报方法及系统 |
CN108022019B (zh) * | 2017-12-14 | 2022-02-11 | 西南石油大学 | 基于风电机组分类的风电场调度方法和系统 |
CN108022019A (zh) * | 2017-12-14 | 2018-05-11 | 西南石油大学 | 基于风电机组分类的风电场调度方法和系统 |
CN108599147B (zh) * | 2018-04-13 | 2021-01-05 | 华北电力大学 | 基于正态指数平滑法与核密度估计的组合区间预测方法 |
CN108599147A (zh) * | 2018-04-13 | 2018-09-28 | 华北电力大学 | 基于正态指数平滑法与核密度估计的组合区间预测方法 |
CN110858715A (zh) * | 2018-08-25 | 2020-03-03 | 南京理工大学 | 基于风电随机性和相关性的风电穿透功率极限分析方法 |
CN109978222A (zh) * | 2019-01-04 | 2019-07-05 | 国网山东省电力公司电力科学研究院 | 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 |
CN109800919A (zh) * | 2019-01-30 | 2019-05-24 | 云南电网有限责任公司电力科学研究院 | 一种区域风电场整体有功功率预测方法及装置 |
CN109800919B (zh) * | 2019-01-30 | 2023-07-11 | 云南电网有限责任公司电力科学研究院 | 一种区域风电场整体有功功率预测方法及装置 |
CN111815039A (zh) * | 2020-06-28 | 2020-10-23 | 山东大学 | 基于天气分类的周尺度风电功率概率预测方法及系统 |
CN112269967A (zh) * | 2020-10-14 | 2021-01-26 | 清华大学 | 一种基于联合机会约束的迭代拆分方法及系统 |
CN113487105B (zh) * | 2021-07-27 | 2023-04-07 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
CN113487105A (zh) * | 2021-07-27 | 2021-10-08 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN105868853B (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105868853B (zh) | 一种短期风电功率组合概率预测方法 | |
Cao et al. | Hybrid ensemble deep learning for deterministic and probabilistic low-voltage load forecasting | |
Li et al. | Photovoltaic power forecasting with a hybrid deep learning approach | |
Alessandrini et al. | A novel application of an analog ensemble for short-term wind power forecasting | |
CN103218675B (zh) | 一种基于聚类和滑动窗口的短期负荷预测方法 | |
CN108985965A (zh) | 一种结合神经网络和参数估计的光伏功率区间预测方法 | |
CN108388962B (zh) | 一种风电功率预测系统及方法 | |
Zhang et al. | GEFCom2014 probabilistic solar power forecasting based on k-nearest neighbor and kernel density estimator | |
CN107730044A (zh) | 一种可再生能源发电和负荷的混合预测方法 | |
CN105373857A (zh) | 一种光伏电站辐照度预测方法 | |
CN105225006B (zh) | 一种短期风电功率非参数概率预测方法 | |
CN105512745A (zh) | 基于粒子群—bp神经网络的风电功率区间预测方法 | |
CN111612244B (zh) | 基于qra-lstm的日前光伏功率非参数概率预测方法 | |
CN114330935B (zh) | 基于多种结合策略集成学习的新能源功率预测方法和系统 | |
Mandal et al. | A new recursive neural network algorithm to forecast electricity price for PJM day‐ahead market | |
Qiu et al. | Fuzzy time series forecasting model based on automatic clustering techniques and generalized fuzzy logical relationship | |
Nagaraja et al. | A survey on wind energy, load and price forecasting:(Forecasting methods) | |
Aliberti et al. | Forecasting Short-term Solar Radiation for Photovoltaic Energy Predictions. | |
Dong et al. | Short-term building cooling load prediction model based on DwdAdam-ILSTM algorithm: A case study of a commercial building | |
CN109087215A (zh) | 多风电场输出功率联合概率密度预测方法 | |
CN110222910A (zh) | 一种有源配电网态势预测方法及预测系统 | |
Tan et al. | Wind power scenario generation with non-separable spatio-temporal covariance function and fluctuation-based clustering | |
CN112307672A (zh) | 基于布谷鸟算法优化的bp神经网络短期风功率预测方法 | |
Zhang | Short‐term power load forecasting based on SAPSO‐CNN‐LSTM model considering autocorrelated errors | |
CN114372615A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |