CN105046383A - 一种基于集合经验模态分解和相关向量机的风电功率实时预测方法 - Google Patents
一种基于集合经验模态分解和相关向量机的风电功率实时预测方法 Download PDFInfo
- Publication number
- CN105046383A CN105046383A CN201510599478.5A CN201510599478A CN105046383A CN 105046383 A CN105046383 A CN 105046383A CN 201510599478 A CN201510599478 A CN 201510599478A CN 105046383 A CN105046383 A CN 105046383A
- Authority
- CN
- China
- Prior art keywords
- prediction
- wind power
- vector machine
- mode decomposition
- empirical mode
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明是一种基于集合经验模态分解和相关向量机的风电功率实时预测方法,其特点是包括数据获取及处理、建立多步滚动预测模式、建立集合经验模态分解和相关向量机的风电功率预测模型和预测评价指标等步骤,具有科学合理,能够降低风电功率的非线性、非平稳性带来的影响,满足预测精度要求,利于电力***调度,保证电能质量,降低运行成本等优点。
Description
技术领域
本发明涉及风电技术领域,是一种基于集合经验模态分解和相关向量机的风电功率实时预测方法。
背景技术
全球性的能源危机的变化使新能源的需求迅速增加,风能是目前最具大规模开发利用潜力的可再生能源,风力发电是大规模利用风能的有效途径,也是我国能源和电力可持续发展战略的最现实选择。然而近年来,随着风电装机容量的快速增长,风电在电网中的比例逐年增大。由于风能本身的波动性和随机性,当风电穿透功率超过一定值后,将会对电力***的调度运行和电能质量带来严峻挑战,这严重限制了风电的发展。如果能有效地对风电功率进行预测,不仅可以减少电力***的备用容量、降低***运行成本,而且还可以有效提高电力***中风电最大装机比例,提高风电竞争力。
风电功率实时预测是指自预测时刻起对未来15分钟至4小时的风电功率进行预测,采样间隔为15分钟,一次预测数据为16个,由此可知风电功率实时预测为超短期多步预测。
发明内容
本发明的目的是提供一种科学合理,能够降低风电功率的非线性、非平稳性带来的影响,满足预测精度要求,利于电力***调度,保证电能质量,降低运行成本的基于集合经验模态分解和相关向量机的风电功率实时预测方法。
本发明的目的是由以下技术方案来实现的:一种基于集合经验模态分解和相关向量机的风电功率实时预测方法,其特征在于,它包括以下步骤:
(1)数据获取及处理
选取某风电场整场机组以15分钟为时间间隔的实际风电功率数据,将历史数据分为两部分,前一部分作为训练样本,后一部分作为预测样本;
(2)建立多步滚动预测模式
进行风电功率预测时,一般假设当前时刻记为i,采样间隔记为i*,已知历史数据y(i-ni*),n=0,1,2,…,N,预测的值为y(i+mi*),m=1,2,…,M,M为多步预测的步数,则滚动多步预测的预测值可以表示为:
yG(i+mi*)=f(y(i-(N-m+1)i*),
(1)
y(i),yG(i+i*),…,yG(i+(m-1)i*))
(3)建立集合经验模态分解和相关向量机的风电功率预测模型
①采用集合经验模态对原始序列分解
集合经验模态分解(EEMD)是将多组不同的白噪声序列加入到原始序列之上,然后分别对其进行经验模态分解(EMD),之后将相应的固有模态分量(IMF)的均值看作是真实分量,每个固有模态分量(IMF)需要同时满足信号的零点数和极值点数最多相差一个;任意时刻,由局部极小值点和极大值点定义的包络线的均值为零,即上下曲线关于坐标轴局部对称,具体内容是:
a.向原始序列x(t)中加入白噪声序列,白噪声序列服从(0,(αε)2)的正态分布,得到复合序列X(t),
b.对复合序列X(t)进行经验模态分解(EMD),得到各阶固有模态分量(IMF),此时其中ci(t)为各阶固有模态分量(IMF),rn(t)为剩余分量,代表原信号的平均趋势,
c.重复步骤(a)、(b)r次,每次加入的白噪声幅值不同,
d.利用高斯白噪声频谱的零均值原理,原始信号对应的固有模态分量(IMF)为
②建立相关向量机预测模型
相关向量机(RVM)是一种稀疏概率模型,它在先验参数的结构基础上,把样本数据进行迭代,使得大部分参数的后验概率分布趋于零,将模型稀疏化,对集合经验模态分解(EEMD)得到的各分量建立相关向量机(RVM)预测模型;
③叠加对各分量得到最终的预测结果;
(4)预测评价指标
根据中华人民共和国能源局的有关规定,采用均方根误差M、准确率zq和合格率r作为评价的指标,当均方根误差越小、准确率越大、合格率越大时,预测精度越高,各指标的计算式如下:
全天预测的均方根误差M:
式中,M为均方根误差,是第i次实测功率时第k时刻的预测值,是第i次实测功率时第k时刻功率的实际值,pcap是风电总场的装机容量,
日平均预测准确率zq表示为:
式中,zqi为第i次预测的准确率,
日平均预测合格率r为:
式中,ri为第i次预测的合格率。
本发明的基于集合经验模态分解和相关向量机的风电功率实时预测方法的优点体现在:
1.集合经验模态分解很好地解决了经验模态分解算法在分解过程中出现的端点效应和模态混叠等现象,能够有效降低信号的非平稳特性;
2.相关向量机是一种稀疏概率模型,具有很好的非线性拟合能力;
3.预测精度高;
4.适用于风电功率的实时预测,将有助于电力***调度部门根据风电功率的变化及时作出调整计划,保证电能质量,降低运行成本,对实现风电规范化并网至关重要,其科学合理,实用价值高。
附图说明
图1为集合经验模态分解流程图;
图2为在本发明的预测方法下,向阳风电场风电功率超短期预测曲线与实际功率曲线的对比示意图。
具体实施方式
下面结合附图和具体实施例,对本发明的基于集合经验模态分解和相关向量机的风电功率实时预测方法作进一步说明。
本发明的基于集合经验模态分解和相关向量机的风电功率实时预测方法,包括以下步骤:
(1)数据获取及处理
本发明以吉林省向阳风电场2012年7月份的实测数据为例进行分析。该风电场装机容量为400.5MW,风机数量为267台,单台风机的额定容量为1500kW,数据采样间隔为15分钟;
(2)建立多步滚动预测模式
进行风电功率预测时,一般假设当前时刻记为i,采样间隔记为i*,已知历史数据y(i-ni*),n=0,1,2,…,N,预测的值为y(i+mi*),m=1,2,…,M,M为多步预测的步数,则滚动多步预测的预测值可以表示为:
yG(i+mi*)=f(y(i-(N-m+1)i*),
(1)
y(i),yG(i+i*),…,yG(i+(m-1)i*))
(3)建立集合经验模态分解(EEMD)和相关向量机(RVM)的风电功率预测模型
①采用集合经验模态对原始序列分解
集合经验模态分解(EEMD)是对传统经验模态分解的模态混叠等现象进行改进的一种自适应筛选方法,可以降低信号的非平稳特性。为了得到数据的实际信号,就需将多组不同的白噪声序列加入到原始序列之上,然后分别对其进行EMD分解,之后将相应的固有模态分量(IMF)的均值看作是真实分量。每个固有模态分量(IMF)需要同时满足信号的零点数和极值点数最多相差一个;任意时刻,由局部极小值点和极大值点定义的包络线的均值为零,即上下曲线关于坐标轴局部对称。集合经验模态分解(EEMD)流程图如图1,具体内容是:
a.向原始序列x(t)中加入白噪声序列,白噪声序列服从(0,(αε)2)的正态分布,得到复合序列X(t),
b.对复合序列X(t)进行EMD分解,得到各阶固有模态分量(IMF),此时其中ci(t)为各阶固有模态分量(IMF),rn(t)为剩余分量,代表原信号的平均趋势,
c.重复步骤(a)、(b)r次,每次加入的白噪声幅值不同,
d.利用高斯白噪声频谱的零均值原理,原始信号对应的IMF分量为
②建立相关向量机(RVM)预测模型
相关向量机(RVM)是一种稀疏概率模型,它在先验参数的结构基础上,把样本数据进行迭代,使得大部分参数的后验概率分布趋于零,将模型稀疏化,具有很好的非线性拟合能力,对集合经验模态分解(EEMD)得到的各分量建立相关向量机(RVM)预测模型;
③叠加对各分量得到最终的预测结果;
(4)预测评价指标
根据可中华人民共和国能源局对风电功率实时预测的管理规定提出的评价指标对误差进行分析。
表1和表2以7月10日为例,分别对87号风机和整场267台机组风电功率在各预测方法下进行的预测精度评价指标。
表187号机组风电功率预测结果评价指标
根据中华人民共和国能源局发布的《风电场功率预测预报管理暂行办法》中对全天预测的均方根误差应该小于20%,由表1和表2可以看出,无论是单台机组还是整场机组,eemd-rvm预测方法都满足要求,不仅均方根误差小于20%,准确率和合格率也都很高,说明了eemd-rvm模型的有效性。
87号机组在7月10日的第10次预测的预测值和真实值对比图如图2所示。可以发现,预测值和真实值接近,说明eemd-rvm模型的适用性。
本发明的特定实施例已对本发明的内容作出了详尽的说明,但不局限本实施例,本领域技术人员根据本发明的启示所做的任何显而易见的改动,都属于本发明权利保护的范围。
Claims (1)
1.一种基于集合经验模态分解和相关向量机的风电功率实时预测方法,其特征在于,它包括以下步骤:
(1)数据获取及处理
选取某风电场整场机组以15分钟为时间间隔的实际风电功率数据,将历史数据分为两部分,前一部分作为训练样本,后一部分作为预测样本;
(2)建立多步滚动预测模式
进行风电功率预测时,一般假设当前时刻记为i,采样间隔记为i*,已知历史数据y(i-ni*),n=0,1,2,...,N,预测的值为y(i+mi*),m=1,2,...,M,M为多步预测的步数,则滚动多步预测的预测值可以表示为:
(3)建立集合经验模态分解和相关向量机的风电功率预测模型
①采用集合经验模态对原始序列分解
集合经验模态分解(EEMD)是将多组不同的白噪声序列加入到原始序列之上,然后分别对其进行经验模态分解(EMD),之后将相应的固有模态分量(IMF)的均值看作是真实分量,每个固有模态分量(IMF)需要同时满足信号的零点数和极值点数最多相差一个;任意时刻,由局部极小值点和极大值点定义的包络线的均值为零,即上下曲线关于坐标轴局部对称,具体内容是:
a.向原始序列x(t)中加入白噪声序列,白噪声序列服从(0,(αε)2)的正态分布,得到复合序列X(t),
b.对复合序列X(t)进行经验模态分解(EMD),得到各阶固有模态分量(IMF),此时其中ci(t)为各阶固有模态分量(IMF),rn(t)为剩余分量,代表原信号的平均趋势,
c.重复步骤(a)、(b)r次,每次加入的白噪声幅值不同,
d.利用高斯白噪声频谱的零均值原理,原始信号对应的固有模态分量(IMF)为
②建立相关向量机预测模型
相关向量机(RVM)是一种稀疏概率模型,它在先验参数的结构基础上,把样本数据进行迭代,使得大部分参数的后验概率分布趋于零,将模型稀疏化,对集合经验模态分解(EEMD)得到的各分量建立相关向量机(RVM)预测模型;
③叠加对各分量得到最终的预测结果;
(4)预测评价指标
根据中华人民共和国能源局的有关规定,采用均方根误差M、准确率zq和合格率r作为评价的指标,当均方根误差越小、准确率越大、合格率越大时,预测精度越高,各指标的计算式如下:
全天预测的均方根误差M:
式中,M为均方根误差,是第i次实测功率时第k时刻的预测值,是第i次实测功率时第k时刻功率的实际值,pcap是风电总场的装机容量,
日平均预测准确率zq表示为:
式中,zqi为第i次预测的准确率,
日平均预测合格率r为:
式中,ri为第i次预测的合格率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510599478.5A CN105046383A (zh) | 2015-09-19 | 2015-09-19 | 一种基于集合经验模态分解和相关向量机的风电功率实时预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510599478.5A CN105046383A (zh) | 2015-09-19 | 2015-09-19 | 一种基于集合经验模态分解和相关向量机的风电功率实时预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105046383A true CN105046383A (zh) | 2015-11-11 |
Family
ID=54452910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510599478.5A Pending CN105046383A (zh) | 2015-09-19 | 2015-09-19 | 一种基于集合经验模态分解和相关向量机的风电功率实时预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105046383A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108876039A (zh) * | 2018-06-21 | 2018-11-23 | 浙江工业大学 | 一种基于支持向量机的含分布式电源配电网电能质量预测方法 |
CN111191398A (zh) * | 2019-12-06 | 2020-05-22 | 云南电网有限责任公司玉溪供电局 | 基于svr的变电站直流***蓄电池退化趋势预测方法 |
CN112329339A (zh) * | 2020-10-27 | 2021-02-05 | 河北工业大学 | 一种风电场短期风速预测方法 |
CN115496627A (zh) * | 2022-11-22 | 2022-12-20 | 国网天津市电力公司电力科学研究院 | 一种可调资源响应潜力的评估方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521671A (zh) * | 2011-11-29 | 2012-06-27 | 华北电力大学 | 一种风电功率超短期预测方法 |
CN102567636A (zh) * | 2011-12-27 | 2012-07-11 | 华北电力大学 | 基于相关向量机的高精度风电场功率区间预测方法 |
CN103400052A (zh) * | 2013-08-22 | 2013-11-20 | 武汉大学 | 一种风电场短期风速组合预测方法 |
CN104239689A (zh) * | 2014-08-18 | 2014-12-24 | 华北电力大学 | 基于优化相关向量机的短期风速预测方法 |
-
2015
- 2015-09-19 CN CN201510599478.5A patent/CN105046383A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521671A (zh) * | 2011-11-29 | 2012-06-27 | 华北电力大学 | 一种风电功率超短期预测方法 |
CN102567636A (zh) * | 2011-12-27 | 2012-07-11 | 华北电力大学 | 基于相关向量机的高精度风电场功率区间预测方法 |
CN103400052A (zh) * | 2013-08-22 | 2013-11-20 | 武汉大学 | 一种风电场短期风速组合预测方法 |
CN104239689A (zh) * | 2014-08-18 | 2014-12-24 | 华北电力大学 | 基于优化相关向量机的短期风速预测方法 |
Non-Patent Citations (3)
Title |
---|
严干贵等: "《两种风电功率多步预测方式的分析及评价》", 《东北电力大学学报》 * |
敖培等: "《于最小机会损失准则的超短期风电功率组合预测》", 《计算机光盘软件与应用》 * |
肖燕彩等: "《基于EMD与SVM的风电功率短期预测》", 《北京交通大学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108876039A (zh) * | 2018-06-21 | 2018-11-23 | 浙江工业大学 | 一种基于支持向量机的含分布式电源配电网电能质量预测方法 |
CN108876039B (zh) * | 2018-06-21 | 2021-07-27 | 浙江工业大学 | 一种基于支持向量机的含分布式电源配电网电能质量预测方法 |
CN111191398A (zh) * | 2019-12-06 | 2020-05-22 | 云南电网有限责任公司玉溪供电局 | 基于svr的变电站直流***蓄电池退化趋势预测方法 |
CN112329339A (zh) * | 2020-10-27 | 2021-02-05 | 河北工业大学 | 一种风电场短期风速预测方法 |
CN112329339B (zh) * | 2020-10-27 | 2022-07-19 | 河北工业大学 | 一种风电场短期风速预测方法 |
CN115496627A (zh) * | 2022-11-22 | 2022-12-20 | 国网天津市电力公司电力科学研究院 | 一种可调资源响应潜力的评估方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features | |
CN102663513B (zh) | 利用灰色关联度分析的风电场功率组合预测建模方法 | |
CN102479339B (zh) | 基于混合神经网络的风电场短期风速预测方法和*** | |
Zhang et al. | Improved Wind Speed Prediction Using Empirical Mode Decomposition. | |
CN105184391A (zh) | 基于小波分解和支持向量机的风电场风速与功率预测方法 | |
CN101793907A (zh) | 一种风电场短期风速预测方法 | |
CN105426956A (zh) | 一种超短期光伏预测方法 | |
CN102496926B (zh) | 风电场功率预测输入数据的判断及处理方法 | |
CN103208037B (zh) | 一种应用于新能源电站的基于在线校正的功率预测方法 | |
CN103683326A (zh) | 一种区域电网风电多点接入最佳接纳能力的计算方法 | |
CN102831321A (zh) | 一种基于蒙特卡洛方法的风电场风险评估方法 | |
CN107609697A (zh) | 一种风力发电功率预测方法 | |
CN105046383A (zh) | 一种基于集合经验模态分解和相关向量机的风电功率实时预测方法 | |
CN109376951A (zh) | 一种光伏概率预测方法 | |
CN103996079A (zh) | 一种基于条件概率的风电功率加权预测方法 | |
CN102479347A (zh) | 基于数据驱动的风电场短期风速预测方法和*** | |
CN106611243A (zh) | 一种基于garch模型的风速预测残差修正方法 | |
CN103440529A (zh) | 一种风电功率预测***的横向误差平移修正方法 | |
Yang et al. | Short-term PV generation system direct power prediction model on wavelet neural network and weather type clustering | |
CN105896599A (zh) | 一种考虑风机特性的风机出力多指标评估方法 | |
CN104834980A (zh) | 一种用于优化绿色数据中心电费的负载调度方法 | |
CN103927597A (zh) | 基于自回归滑动平均模型的风电功率超短期预测方法 | |
CN103530508B (zh) | 一种建立风速功率转换概率模型的方法 | |
CN105741192B (zh) | 一种风电场风电机机舱短期风速组合预报方法 | |
CN105046349A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151111 |
|
RJ01 | Rejection of invention patent application after publication |