CN102509027A - 一种基于交叉熵理论的风电功率组合预测方法 - Google Patents
一种基于交叉熵理论的风电功率组合预测方法 Download PDFInfo
- Publication number
- CN102509027A CN102509027A CN2011104268595A CN201110426859A CN102509027A CN 102509027 A CN102509027 A CN 102509027A CN 2011104268595 A CN2011104268595 A CN 2011104268595A CN 201110426859 A CN201110426859 A CN 201110426859A CN 102509027 A CN102509027 A CN 102509027A
- Authority
- CN
- China
- Prior art keywords
- wind power
- cross entropy
- wind
- power
- individual
- 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
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于交叉熵理论的风电功率组合预测方法,包括以下步骤:(1)采用多种方法计算给定时间段内某一时刻的风电功率;(2)对历史时刻进行虚拟预测;(3)定义风电功率概率密度函数;(4)建立支持向量,并确定其与权值的关系;(5)求解最小交叉熵问题。本发明对能够对多种预测方法进行综合,利用交叉熵理论,能够通过判断不同信息源之间的两两相交程度来确定相互支持度,并依据相互支持度确定信息源的权重,为优化确定风电功率组合预测模型的权重提供了重要的参考依据。
Description
技术领域
本发明涉及一种用于风力发电功率组合预测方法,属于新能源发电技术中的风电功率预测技术领域。
背景技术
由于风的高度随机性、波动性和间歇性,大规模风电接入电网会对电力供需平衡、电力***的安全稳定以及电能质量等带来严峻挑战。风电功率预测对于电力***的调度和安全稳定运行有重要的意义,能够优化***的发电计划和备用容量配置,减少电力***运行成本,为***调度和实时运行控制提供依据,有效减轻风电接入对电网的影响,提高风电的利用率。
目前,根据预测的物理量不同,风电功率预测方法可以分为2类:第1类是对风速的预测,然后根据风电场的功率曲线得到其功率输出;第2类是直接预测风电场的输出功率。在具体实施时,通常分为三种方法,即物理方法、统计方法和物理与统计相结合的方法。其中,通过对风速、风向、气温、气压、地形和粗糙度等信息进行物理建模,再利用物理方程进行对风电功率进行预测的方法是物理方法;根据历史数据进行统计分析,找到其内在规律并用于预测的方法是统计方法;将物理方法与统计方法综合应用于风电功率预测的方法成为物理与统计相结合方法。所有预测方法都需建立一定的数学模型,根据数学模型的不同又可分为持续时间法、ARMA法、卡尔曼滤波法、人工神经网络法、小波分析法、支持向量机(SVM)回归法等。其中,持续时间法最为简单,即把最近一点的风速或功率观测值作为下一点的预测值。上述预测方法从某一个方面来说,都采用了单一的预测模型,导致某些测量点往往出现较大误差。
国内外已有专家学者开始研究多种模型组合的风电功率预测方法。虽然这些方法对提高风电功率预测精度产生了有益效果,但在确定组合模型权重时只是通过简单的数学处理,并未考虑到每种预测模型输入信息之间的相互融合特性,可能导致数据融合结果的不可靠。
随着信息科学的崛起和信息技术的不断发展,电力***逐步实现信息化,信息熵也在电力***可靠性评估中得到了应用[9-13]。其中,交叉熵[14]是其中常用的方法之一。交叉熵能够通过判断不同信息源之间的两两相交程度来确定相互支持度,并依据相互支持度确定信息源的权重,相互支持度越高,所占权重越大。这为进一步优化确定风电功率组合预测模型的权重提供了重要的参考依据。本发明将风电功率组合预测看作一个信息融合问题,利用交叉熵理论确定各个预测模型的权值,建立具有更高精度的风电功率预测组合模型。
发明内容
本发明所要解决的技术问题是为提高风电功率预测精度,提供一种基于交叉熵理论的风电功率组合预测方法。
本发明解决其技术问题所采用的技术方案如下:一种基于交叉熵理论的风电功率组合预测方法,其特征是:所述的基于交叉熵理论的风电功率组合预测方法包括以下步骤:
(1)采用N种可用方法预测给定时间段内某一时刻的风电功率,得到每种预测方法的风电功率预测结果为P tn (n=1,2,…,N);
(2)对历史时刻进行虚拟预测,选择与预测时刻风速、风向近似且功率也近似的M个历史时刻,分别用N种算法进行风电功率虚拟预测,每种方法得到M个功率值P nm (m=1,2,…,M),加上原来的P tn ,对于第n个算法,就有(M+1)个功率值,就可计算第n种算法得到的风电功率的数字特征;
(3)定义风电功率概率密度函数,设f(x)为融合各影响因素给出的风电场输出功率概率分布密度函数,X代表风电输出功率,f n (x)为第n种预测模型下的风电输出功率概率分布密度函数,满足以下关系和定义,
(1)
(4)建立支持向量,并确定其与权值的关系;
(5)假设t时刻的风电功率满足正态分布,
式中,μ为平均值,σ为方差,
在此基础上,求解最小交叉熵问题。
其中,步骤(4)的具体步骤是:
a、建立支持向量S,
式中,
进入步骤b;
b、记
为使各种预测方法的权重ω n 能够反应不同信源之间的支持程度,令
进入步骤c;
c、建立最小交叉熵优化问题的目标函数,
其中,步骤(5)的具体步骤是:
a、对于第n个算法得到的风电功率预测值P tn ,做为第n个算法得到的风电功率概率分布的均值μ n ,进入步骤b;
b、根据历史数据的虚拟预测所得到的风电功率预测值,可计算样本方差,作为σ tn ,
进入步骤c;
c、得到这N种算法预测t时刻风电功率的概率密度函数,求解最小交叉熵优化问题,确定f(x)的具体表达式,所求的t时刻风电功率即为X的数学期望。
本发明的有益效果如下:
(1) 本发明提供的基于交叉熵理论的风电功率组合预测方法,能够对多种预测方法进行综合,通过配置合适的权重降低整体预测误差;
(2) 本发明提供的基于交叉熵理论的风电功率组合预测方法,能够通过判断不同信息源之间的两两相交程度来确定相互支持度,并依据相互支持度确定信息源的权重,为优化确定风电功率组合预测模型的权重提供了重要的参考依据。
附图说明
图1为基于交叉熵理论的风电功率组合预测方法流程。
图2为单一风电功率预测方法的预测结果。
图3为单一风电功率预测方法的预测绝对误差。
图4为风电功率组合预测方法的预测结果。
图5为风电功率组合预测方法的预测绝对误差。
具体实施方式
下面参照附图并结合实施例对本发明作进一步详细描述。但是本发明不限于所给出的例子。
实例一:
本发明用于某风电场的风电功率预测。风电场含有单机容量1.5MW的风电机组134台,总容量为200MW,风电功率预测方法选择常用的三种风电功率预测方法:ARMA时序分析法、BP神经网络法和SVM支持向量机法。具体分析时,分别对单一预测方法和采用不同权重计算方法的组合预测方法进行对比分析。
图2是采用不同单一预测方法所得到的2011年4月某日风电功率预测结果,黑线为当日实际功率,蓝线为采用ARMA时序分析方法得到的风电功率预测值,红线为采用BP神经网络方法得到的风电功率预测值,绿线为采用SVM支持向量机方法得到的风电功率预测值,三种常用预测方法都能够比较好的预测风电功率的变化趋势,但具体预测值与实际情况仍存在较大偏差。图3是三种方法预测结果的绝对误差,ARMA时间序列分析方法的最大绝对误差为15.6%,BP神经网络方法的最大绝对误差为15.4%,SVM支持向量机方法的最大绝对误差高达23.12%。
在此基础上,分析采用等权重组合预测方法和基于交叉熵的组合预测方法的预测结果。采用等权重方法时,三种预测方法的权重均取为1/3,采用交叉熵方法时,三种预测方法的权重由本发明所述步骤求取。图4为两种组合预测方法的预测结果,黑线为当日实际功率,蓝线为等权重组合预测方法的风电功率预测值,红线为采用基于交叉熵理论的组合预测方法得到的风电功率预测值,与图2相比,两种组合预测方法的预测结果更加符合实际功率值。图5给出了两种组合预测方法预测结果的绝对误差,等权重组合预测方法的最大绝对误差为14.9%,基于交叉熵理论的组合预测方法的最大绝对误差约为13%。
Claims (2)
1.一种基于交叉熵理论的风电功率组合预测方法,其特征是:所述的基于交叉熵理论的风电功率组合预测方法包括以下步骤:
(1)采用N种可用方法预测给定时间段内某一时刻的风电功率,得到每种预测方法的风电功率预测结果为P tn (n=1,2,…,N);
(2)对历史时刻进行虚拟预测,选择与预测时刻风速、风向近似且功率也近似的M个历史时刻,分别用N种算法进行风电功率虚拟预测,每种方法得到M个功率值P nm (m=1,2,…,M),加上原来的P tn ,对于第n个算法,就有(M+1)个功率值,就可计算第n种算法得到的风电功率的数字特征;
(3)定义风电功率概率密度函数,设f(x)为融合各影响因素给出的风电场输出功率概率分布密度函数,X代表风电输出功率,f n (x)为第n种预测模型下的风电输出功率概率分布密度函数,满足以下关系和定义,
(2)
(4)建立支持向量,并确定其与权值的关系;
(5)假设t时刻的风电功率满足正态分布,
式中,μ为平均值,σ为方差,
在此基础上,求解最小交叉熵问题。
2.根据权利要求1所述的基于交叉熵理论的风电功率组合预测方法,其特征是:步骤(4)的具体步骤是:
a、建立支持向量S,
式中,
进入步骤b;
b、记
为使各种预测方法的权重ω n 能够反应不同信源之间的支持程度,令
进入步骤c;
c、建立最小交叉熵优化问题的目标函数,
3、根据权利要求1所述的基于交叉熵理论的风电功率组合预测方法,其特征是:其中,步骤(5)的具体步骤是:
a、对于第n个算法得到的风电功率预测值P tn ,做为第n个算法得到的风电功率概率分布的均值μ n ,进入步骤b;
b、根据历史数据的虚拟预测所得到的风电功率预测值,可计算样本方差,作为σ tn ,
进入步骤c;
c、得到这N种算法预测t时刻风电功率的概率密度函数,求解最小交叉熵优化问题,确定f(x)的具体表达式,所求的t时刻风电功率即为X的数学期望。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110426859.5A CN102509027B (zh) | 2011-12-19 | 2011-12-19 | 一种基于交叉熵理论的风电功率组合预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110426859.5A CN102509027B (zh) | 2011-12-19 | 2011-12-19 | 一种基于交叉熵理论的风电功率组合预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102509027A true CN102509027A (zh) | 2012-06-20 |
CN102509027B CN102509027B (zh) | 2014-11-26 |
Family
ID=46221112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110426859.5A Active CN102509027B (zh) | 2011-12-19 | 2011-12-19 | 一种基于交叉熵理论的风电功率组合预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102509027B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440428A (zh) * | 2013-09-12 | 2013-12-11 | 重庆大学 | 风电功率组合预测模型的自适应动态权重确定方法 |
CN103473461A (zh) * | 2013-09-16 | 2013-12-25 | 东南大学 | 基于数据特征提取的风功率预测误差估计方法 |
CN104201710A (zh) * | 2014-07-25 | 2014-12-10 | 国家电网公司 | 应用最大熵原理的风电最大装机容量预测方法 |
CN104332996A (zh) * | 2014-11-18 | 2015-02-04 | 国家电网公司 | 一种评估电力***可靠性的方法 |
CN104915736A (zh) * | 2015-06-29 | 2015-09-16 | 东北电力大学 | 基于改进熵权法提高风电功率组合预测精度的方法 |
CN106452863A (zh) * | 2016-09-30 | 2017-02-22 | 国家电网公司 | 一种信息***的资源分配及预测方法和装置 |
CN107092779A (zh) * | 2017-03-30 | 2017-08-25 | 北京空间飞行器总体设计部 | 一种基于最小交叉位移熵的传感器作动器位置优化方法 |
CN112613655A (zh) * | 2020-12-18 | 2021-04-06 | 国网新疆电力有限公司经济技术研究院 | 一种风电功率组合预测方法 |
CN113283657A (zh) * | 2021-06-01 | 2021-08-20 | 国网山东省电力公司嘉祥县供电公司 | 一种电力应急物资需求的预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1828668A (zh) * | 2006-04-10 | 2006-09-06 | 天津大学 | 基于嵌入式隐马尔可夫模型与交叉熵的台风中心定位方法 |
CN101414751A (zh) * | 2008-11-20 | 2009-04-22 | 北京方鸿溪科技有限公司 | 风电功率预测***及其方法、电网*** |
-
2011
- 2011-12-19 CN CN201110426859.5A patent/CN102509027B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1828668A (zh) * | 2006-04-10 | 2006-09-06 | 天津大学 | 基于嵌入式隐马尔可夫模型与交叉熵的台风中心定位方法 |
CN101414751A (zh) * | 2008-11-20 | 2009-04-22 | 北京方鸿溪科技有限公司 | 风电功率预测***及其方法、电网*** |
Non-Patent Citations (2)
Title |
---|
李永红等: "基于交叉熵的可靠性评估技术", 《航空标准化与质量》 * |
韩爽: "风电场功率短期预测方法研究", 《万方学位论文数据库》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440428A (zh) * | 2013-09-12 | 2013-12-11 | 重庆大学 | 风电功率组合预测模型的自适应动态权重确定方法 |
CN103473461A (zh) * | 2013-09-16 | 2013-12-25 | 东南大学 | 基于数据特征提取的风功率预测误差估计方法 |
CN103473461B (zh) * | 2013-09-16 | 2016-06-15 | 东南大学 | 基于数据特征提取的风功率预测误差估计方法 |
CN104201710A (zh) * | 2014-07-25 | 2014-12-10 | 国家电网公司 | 应用最大熵原理的风电最大装机容量预测方法 |
CN104332996A (zh) * | 2014-11-18 | 2015-02-04 | 国家电网公司 | 一种评估电力***可靠性的方法 |
CN104915736A (zh) * | 2015-06-29 | 2015-09-16 | 东北电力大学 | 基于改进熵权法提高风电功率组合预测精度的方法 |
CN106452863A (zh) * | 2016-09-30 | 2017-02-22 | 国家电网公司 | 一种信息***的资源分配及预测方法和装置 |
CN106452863B (zh) * | 2016-09-30 | 2019-09-24 | 国家电网公司 | 一种信息***的资源分配及预测方法和装置 |
CN107092779A (zh) * | 2017-03-30 | 2017-08-25 | 北京空间飞行器总体设计部 | 一种基于最小交叉位移熵的传感器作动器位置优化方法 |
CN107092779B (zh) * | 2017-03-30 | 2019-11-29 | 北京空间飞行器总体设计部 | 一种基于最小交叉位移熵的传感器作动器位置优化方法 |
CN112613655A (zh) * | 2020-12-18 | 2021-04-06 | 国网新疆电力有限公司经济技术研究院 | 一种风电功率组合预测方法 |
CN113283657A (zh) * | 2021-06-01 | 2021-08-20 | 国网山东省电力公司嘉祥县供电公司 | 一种电力应急物资需求的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102509027B (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102509027A (zh) | 一种基于交叉熵理论的风电功率组合预测方法 | |
Shi et al. | Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features | |
Li et al. | Short-term operation model and risk management for wind power penetrated system in electricity market | |
CN103049798B (zh) | 一种应用于光伏发电***的短期发电功率预测方法 | |
CN107591844A (zh) | 考虑节点注入功率不确定性的主动配电网鲁棒重构方法 | |
CN103683274B (zh) | 区域中长期风电发电量概率预测方法 | |
CN103390116A (zh) | 采用分步方式的光伏电站发电功率预测方法 | |
CN103268366A (zh) | 一种适用于分散式风电场的组合风电功率预测方法 | |
CN102184453A (zh) | 基于模糊神经网络和支持向量机的风电功率组合预测方法 | |
CN103699944A (zh) | 一种多预测模式的风电光伏发电功率预测*** | |
CN102509026A (zh) | 基于最大信息熵原理的风电场短期输出功率预测综合模型 | |
CN103904695A (zh) | 一种基于mcs-pso的临近海岛多微网动态调度方法 | |
CN106684885A (zh) | 基于多场景分析的风电机组配电网无功优化方法 | |
CN110350518A (zh) | 一种用于调峰的电网储能容量需求评估方法及*** | |
CN103971175A (zh) | 一种多级变电站的短期负荷预测方法 | |
CN104037761A (zh) | 一种agc功率多目标随机优化分配方法 | |
CN103326394B (zh) | 一种计及风电波动性的多场景概率优化调度方法 | |
Zhou et al. | Interval optimization combined with point estimate method for stochastic security-constrained unit commitment | |
CN115114854A (zh) | 一种虚拟电厂分布式资源两级自组织优化聚合方法及*** | |
CN116526469A (zh) | 一种水风光互补***长期随机动态调度方法 | |
CN103996087A (zh) | 一种新能源发电功率预测方法及*** | |
CN107358059A (zh) | 短期光伏能量预测方法和装置 | |
Zhou et al. | Distributionally robust unit commitment considering unimodality-skewness information of wind power uncertainty | |
CN113139698A (zh) | 负荷预测方法、装置及设备 | |
Tseng et al. | OPTIMIZATION OF WIND TURBINE PLACEMENT LAYOUT ON NON-FLAT TERRAINS. |
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 |