CN102938562B - 一种区域内风电总功率的预测方法 - Google Patents
一种区域内风电总功率的预测方法 Download PDFInfo
- Publication number
- CN102938562B CN102938562B CN201210250429.7A CN201210250429A CN102938562B CN 102938562 B CN102938562 B CN 102938562B CN 201210250429 A CN201210250429 A CN 201210250429A CN 102938562 B CN102938562 B CN 102938562B
- Authority
- CN
- China
- Prior art keywords
- wind
- power
- energy turbine
- turbine set
- output
- 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.)
- Active
Links
- 230000005611 electricity Effects 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000013507 mapping Methods 0.000 claims abstract description 17
- 238000013528 artificial neural network Methods 0.000 claims abstract description 4
- 210000002569 neuron Anatomy 0.000 claims description 16
- 238000010586 diagram Methods 0.000 claims description 6
- 206010008190 Cerebrovascular accident Diseases 0.000 claims description 4
- 208000006011 Stroke Diseases 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010248 power generation Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 8
- 230000001537 neural effect Effects 0.000 description 5
- 230000005684 electric field Effects 0.000 description 4
- 238000003066 decision tree Methods 0.000 description 3
- 238000003064 k means clustering Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
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
Landscapes
- Wind Motors (AREA)
Abstract
本发明属于风力发电功率预测技术领域,涉及一种面向区域内风电总功率的预测方法。该方法包括:A)在区域内选择代表风电场,并进一步获取各代表风电场的输出功率预测值;B)通过径向基函数神经网络建立映射模型,并根据各代表风电场的输出功率预测值计算出区域内风电总功率预测值;C)输出区域内风电总功率预测值。该预测方法具有使用操作方便、输入输出简单、预测准确度高、误差小,可扩展性好、运用灵活等优点,当代表风电场的预测功率已知,可直接通过映射模型计算出区域内风电总功率预测值;当代表风电场的预测功率未知,可将映射模型接入代表风电场,通过代表风电场的NWP数据来驱动模型,进一步计算出区域内风电总功率预测值。
Description
技术领域
本发明属于风力发电功率预测技术领域,涉及一种面向区域内风电总功率的预测方法。
背景技术
风电输出功率具有很强的随机性,目前的预报水平还不能满足电力系统实际运行的需要,准确的区域风电功率预测有利于电力调度部门安排该区域电力系统的发电计划,保证电能质量,减少系统的备用容量,这对于降低电力系统运行成本有重要的意义。目国内风电功率预测主要集中在对单个风电场功率的预测,对于某一区域内所属风电场总功率的预测相对较少。
区域风电总功率指特定地理范围或者行政隶属范围内所有并网风电场在相同时刻接入电网的功率之和,最简单的预测方式就是把单个风电场的预测结果求和便得出省级全网风电总功率。事实上,这种方式行不通。首先,很多风电场未安装风功率预测系统。其次,各个风电场预测结果准确度良莠不齐,直接相加有可能使误差放大,得出的结果不可信。
实质上,综合考虑到各个风电场之间的地理位置,风功率的间歇性,各个风电场风机运行的稳定性,随机性干扰等因素的影响,本发明通过分析各个风电场以及区域内总风电功率历史数据,建立了一个有效的映射模型以实现区域风电功率预测。
发明内容
为了克服现有技术的,本发明的目的在于提出一种误差小、准确度高、使用方便的区域内风电总功率的预测方法。
本发明是通过如下技术方案来实现的:
一种区域内风电总功率的预测方法,包括如下步骤:
A.在区域内选择代表风电场,并进一步获取各代表风电场的输出功率预测值;
B.通过径向基函数神经网络建立映射模型,并根据各代表风电场的输出功率预测值计算出区域内风电总功率预测值;
C.输出区域内风电总功率预测值。
本发明的预测方法具有如下优点:
(1)该方法使用操作方便,通过建立映射模型实现流程化操作,输入输出简单。
(2)该方法预测准确度高、误差小,可扩展性好。
(3)该方法运用灵活,当代表风电场的预测功率已知(即通过风电功率预测系统直接读取),可直接通过映射模型计算出区域内风电总功率预测值;当代表风电场的预测功率未知,可将映射模型接入代表风电场,通过代表风电场的NWP数据(包括风向和风速)来驱动模型,进一步计算出区域内风电总功率预测值。
附图说明
图1为本发明区域内风电总功率的预测方法的流程示意图;
图2为通过风向和风速的划分规则建立的决策树示意图;
图3为通过径向基函数神经网络(以下简称RBF神经网络)建立的映射模型示意图。
具体实施方式
下面结合附图对本发明的区域内风电总功率的预测方法作进一步详细的说明。
如图3所示,该预测方法可分为如下几步:
步骤1.选择代表风电场。
通过下式(1)求取各风电场的输出功率与区域内风电总功率的相关系数:
式(1)中,xt为t时刻单个风电场的输出功率,yt为t时刻区域内风电总功率,为所有风电场输出功率的平均值,为区域内风电总功率的平均值;
找出与区域内风电总功率有较好相关性(即:相关系数r≥0.8)的风电场作为“代表风电场”,并参与后面的步骤。
步骤2.获取各代表风电场的输出功率预测值,该值可通过通过风电功率预测系统直接读取,也可以通过建立风速功率决策单元、并根据各代表风电场的风向和风速计算得出,下面对第二种方式进行详细说明:
该决策单元的建立是基于把某一电场看成一个整体,由于受到地形,风机排布,风机尾流以及周边建筑物或者其他不确定因素的影响,在此假定沿着某一特定方向,某一特定风速段(大风,中风,小风),上述因素造成此方向风速功率的特定影响,建立一个如图1所示的决策树。在图1中,W(d(t),v(t))表示t时刻风速矢量,其中,
风向划分规则如下:
NEE(东北东方向):0≤d(t)<45;
NEN(东北北方向):45≤d(t)<90;
NWN(西北北方向):90≤d(t)<135;
NWW(西北西方向):135≤d(t)<180;
SWW(西南西方向):180≤d(t)<225;
SWS(西南南方向):225≤d(t)<270;
SES(东南南方向):270≤d(t)<315;
SEE(东南东方向):315≤d(t)<360;
风速划分规则如下:
B(大风):0<v(t)≤5;
M(中风):5<v(t)≤10;
S(小风):v(t)≥10;
风速功率决策单元的建立是按照上述风方向和风速大小划分规则对历史数据进行归类,总共有24类,然后对相应的时刻的风速功率时间序列进行曲线拟合,确定R1到R24的风速功率关系。
下面举例说明归类的方法:在t时刻某一代表风电场的风向为d(t)=55,风速大小为v(t)=8,此刻的功率为P(t)。那么根据上述风向及风速划分规则,t时刻的功率P(t)就应该归类到NEN方向M(中风)类中。
所述曲线拟合的具体步骤包括:
1)通过所读取的数据(P(ti),d(ti))(i=1,2,...,n)绘制散点图,并根据散点图确定曲线类型和曲线拟合方程式中的项数m;
所述曲线拟合方程式为:
P(ti)=amdm(ti)+am-1dm-1(ti)+…+a0(2)
式(2)中,P(ti)、d(ti)分别为t时刻第i个代表风电场的输出功率和风向;a0、……、am-1、am为待确定的系数;m为曲线拟合方程式中的项数,通常取值为3;;
2)建立目标函数:
3)通过最小二乘法确定目标函数中的系数a0……am,并代入曲线拟合方程式(2)中,求得各代表风电场的输出功率预测值,具体步骤包括:
31)令
32)即Da=P,划成正则化方程DTDa=DTP;
33)解正则化方程DTDa=DTP,得a=(DTD)-1DTP,即确定了系数a、并确定了曲线拟合方程式(2)。
注:由于数据量庞大的限制和计算过程的庞杂,举例篇幅巨大,不在此赘述。
当R1到R24确定后,风速功率决策单元建立完毕。
步骤3.建立功率映射模型,并根据各代表风电场的输出功率预测值计算出区域内风电总功率预测值。
为了在已知t时刻所有代表风电场预测功率的基础上对t时刻区域内总风电功率进行预测,在此采用RBF神经网络建立映射模型。
通过RBF神经网络建立的映射模型如附图2所示,包括输入层、隐含层和输出层:设代表风电场有m个,即输入层有m个神经元;设隐含层有n+1个神经元,其中一个神经元的输出值恒为1,即隐含层的阈值设输出层有一个区域总电场;
首先,输入层的每个神经元将通过其输入的一个代表风电场的输出功率预测值分别传至隐含层中除输出值恒为1以外的n个神经元,所述n个神经元的输出值通过下式计算:
式(4)中,Xi(t)为t时刻第i个代表风电场的输出功率预测值,C(t)为t时刻高斯函数的中心,σ为高斯函数的方差,j=1,2……n,i=1,2……6;C(t)可以用K-均值聚类算法求得,K-均值聚类算法在执行过程中会产生若干聚类中心。当C(t)确定后,σ通过公式:获得,其中dmax为所选取中心之间的最大距离。
其次,隐含层中各神经元的输出值与相应权重进行计算后,分别传至输出层进行求和计算,即得到区域内风电总功率预测值。隐含层中各神经元的输出值与相应权重进行计算的具体步骤如下:
本例中,将隐含层中编号为0的神经元的输出值设置为恒等于1,该神经元的输出值与权重w0相乘后得到值w0,并传至输出层;
将隐含层中其它n个神经元的输出值分别与相对应的权重w1(t)……wi(t)……wI(t)通过下式进行计算后传至输出层:
式(5)中,yj(Xi(t))为隐含层中第j个神经元传至输出层的输入值,w0为与隐含层中输出值恒为1的第0个神经元与输出层的连接权重,wj为与隐含层中第j个神经元与输出层的连接权重,权重是通过学习样本计算确定的,确定好以后,在做预测的时候直接用。在做预测的时候无需进行权重的计算;为隐含层中第j个神经元的输出值。
通过用历史数据对神经网络进行学习,确定合适的隐含层节点(即神经元)的个数N,以及合理的权重和阈值,功率映射模型建立完毕。
1)本发明中通过下式来确定隐含层中的神经元个数:
其中,m,n分别为输入层、输出层的神经元个数,m等于代表风电场的个数,n为一个区域总电场,round()意为取整运算。
2)权重是通过下式进行更新确定的:
其中,η为学习速率,为误差目标函数对权重的一阶偏导数,该模型首次运行时,先为所有权重赋值,首次输出的区域内风电总功率预测值P(t),该P(t)值与预先设置的目标值的差值即为上式中的误差值E(n),而每次权重均以上一次权重和误差值为基础进行实时更新。
3)本例中,设置阈值
步骤4.输出区域内风电总功率的预测值P(t)。
本发明综合采用统计,智能算法实现了一种以区域代表点NWP(数值天气预报)数据为输入,区域风电总预测功率为输出的预测模型。通过对区域内单个风电场输出功率与区域总输出率之间的相关性分析,选择相关性较好的风电场作为代表风电场。在此基础上利用历史数据,对风速进行分类,建立不同方向,不同风速水平下的风速功率关系,最终用决策树的思想建立每个代表风电场风速功率决策单元,实现在已知该代表风电场NWP数据的情况下,就可以预测出其输出功率的功能。最后在代表风电场功率已知的情况下,为了预测出区域风电总功率,采用RBF神经网络建立功率映射模型解决了此问题。
最后应该说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,结合上述实施例对本发明进行了详细说明,所属领域的普通技术人员应当理解到:本领域技术人员依然可以对本发明的具体实施方式进行修改或者等同替换,但这些修改或变更均在申请待批的权利要求保护范围之中。
Claims (4)
1.一种区域内风电总功率的预测方法,其特征在于,该方法包括如下步骤:
A.在区域内选择代表风电场,并进一步获取各代表风电场的输出功率预测值;
B.通过径向基函数神经网络建立映射模型,并根据各代表风电场的输出功率预测值计算出区域内风电总功率预测值;
C.输出区域内风电总功率预测值;
所述步骤A中,通过各代表风电场的风向和风速计算得出各代表风电场输出功率预测值,其具体步骤包括:
首先,通过风向和风速的划分规则对各代表风电场的输出功率进行归类;所述风向d(t)的划分规则为:
NEE(东北东方向):0≤d(t)<45;
NEN(东北北方向):45≤d(t)<90;
NWN(西北北方向):90≤d(t)<135;
NWW(西北西方向):135≤d(t)<180;
SWW(西南西方向):180≤d(t)<225;
SWS(西南南方向):225≤d(t)<270;
SES(东南南方向):270≤d(t)<315;
SEE(东南东方向):315≤d(t)<360;
所述风速v(t)的划分规则为:
B(大风):0<v(t)≤5;
M(中风):5<v(t)≤10;
S(小风):v(t)≥10;
其次,对归类的各代表风电场的输出功率进行曲线拟合后,即得到各代表风电场的输出功率预测方程;所述曲线拟合的具体步骤包括:
1)通过所读取的数据P(ti)和d(ti)绘制散点图,并根据散点图确定曲线类型和曲线拟合方程式中的项数m;所述曲线拟合方程式为:
P(ti)=amdm(ti)+am-1dm-1(ti)+…+a0
式中,P(ti)、d(ti)分别为t时刻第i个代表风电场的输出功率和风向,a0、……、am-1、am为待确定的系数,m为项数;
2)建立下式中的目标函数:
3)通过最小二乘法确定目标函数中的a0……am,并代入曲线拟合方程式求得各代表风电场的输出功率预测值。
2.如权利要求1所述的方法,其特征在于,所述步骤A中,选择代表风电场的具体方法包括:
A1.通过下式求取各风电场的输出功率与区域内风电总功率的相关系数:
式中,xt为t时刻某一风电场的输出功率,yt为t时刻区域内风电总功率,为所有风电场输出功率的平均值,为区域内风电总功率的平均值;
A2.找出相关系数r≥0.8的风电场作为代表风电场。
3.如权利要求1所述的方法,其特征在于,所述步骤B的具体方法包括:
所述映射模型包括输入层、隐含层和输出层;设代表风电场有m个,即输入层有m个神经元;设隐含层有n+1个神经元、且其中一个神经元的输出值恒为1;
首先,所述输入层的每个神经元将各代表风电场的输出功率预测值分别传至隐含层中除输出值恒为1以外的其它n个神经元,所述n个神经元的输出值通过下式计算:
式中,Xi(t)为t时刻第i个代表风电场的输出功率预测值,C(t)为t时刻高斯函数的中心,σ为高斯函数的方差,j=1,2……n,i=1,2……6;
其次,所述隐含层中各神经元的输出值与相应权重进行计算后,分别传至输出层进行求和计算,即得到区域内风电总功率预测值。
4.如权利要求3所述的方法,其特征在于,
所述隐含层中一个输出恒为1的神经元的输出值与相应权重w0进行乘积后传至输出层;
所述隐含层中其它n个神经元输出值分别与相应权重通过下式进行计算后传至输出层:
式中,yj(Xi(t))为隐含层中第j个神经元传至输出层的输入值,w0、wj分别为隐含层与输出层的连接权重。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210250429.7A CN102938562B (zh) | 2012-07-19 | 2012-07-19 | 一种区域内风电总功率的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210250429.7A CN102938562B (zh) | 2012-07-19 | 2012-07-19 | 一种区域内风电总功率的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102938562A CN102938562A (zh) | 2013-02-20 |
CN102938562B true CN102938562B (zh) | 2015-07-08 |
Family
ID=47697442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210250429.7A Active CN102938562B (zh) | 2012-07-19 | 2012-07-19 | 一种区域内风电总功率的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102938562B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473621A (zh) * | 2013-09-29 | 2013-12-25 | 中能电力科技开发有限公司 | 风电场短期功率预测方法 |
CN105023070A (zh) * | 2015-08-12 | 2015-11-04 | 河海大学常州校区 | 一种光伏系统输出功率预测方法 |
CN107194499A (zh) * | 2017-04-28 | 2017-09-22 | 国网冀北电力有限公司 | 区域风电短期功率的预测方法及预测装置 |
CN108089135B (zh) * | 2017-12-22 | 2021-03-30 | 广州市香港科大霍英东研究院 | 一种基于极限学习模型的电池状态预测系统及其实现方法 |
CN109063936A (zh) * | 2018-10-12 | 2018-12-21 | 南京千智电气科技有限公司 | 一种建立风电功率预测模型的方法及装置 |
CN109783881A (zh) * | 2018-12-21 | 2019-05-21 | 中国电力科学研究院有限公司 | 一种风电机组发电功率确定方法及装置 |
CN110400011B (zh) * | 2019-07-12 | 2022-02-22 | 国能日新科技股份有限公司 | 电力现货交易中风电场出力申报方案的确定方法及装置 |
CN110277802B (zh) * | 2019-07-24 | 2021-03-02 | 云南电力技术有限责任公司 | 一种用于风电场有功功率的实时预测方法及系统 |
CN116307307B (zh) * | 2023-05-25 | 2023-08-18 | 长江三峡集团实业发展(北京)有限公司 | 一种风电场的超短期功率预测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201413997Y (zh) * | 2009-03-30 | 2010-02-24 | 吉林省电力有限公司 | 风电功率预测系统 |
CN102496927A (zh) * | 2011-12-16 | 2012-06-13 | 中国电力科学研究院 | 基于误差统计修正的风电场功率预测方法 |
CN102570449A (zh) * | 2011-12-21 | 2012-07-11 | 西安交通大学 | 一种基于时空分布特性的区域风电功率预测方法 |
-
2012
- 2012-07-19 CN CN201210250429.7A patent/CN102938562B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201413997Y (zh) * | 2009-03-30 | 2010-02-24 | 吉林省电力有限公司 | 风电功率预测系统 |
CN102496927A (zh) * | 2011-12-16 | 2012-06-13 | 中国电力科学研究院 | 基于误差统计修正的风电场功率预测方法 |
CN102570449A (zh) * | 2011-12-21 | 2012-07-11 | 西安交通大学 | 一种基于时空分布特性的区域风电功率预测方法 |
Non-Patent Citations (2)
Title |
---|
基于RBF神经元网络的风电功率短期预测;武小梅等;《电力系统保护与控制》;20110801;第39卷(第15期);P80-83 * |
韩爽.风电场功率短期预测方法研究.《中国博士学位论文全文数据库(电子期刊)工程科技II辑》.2009,(第2期),C042-11. * |
Also Published As
Publication number | Publication date |
---|---|
CN102938562A (zh) | 2013-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102938562B (zh) | 一种区域内风电总功率的预测方法 | |
Liu et al. | Random forest solar power forecast based on classification optimization | |
Xu et al. | Data-driven configuration optimization of an off-grid wind/PV/hydrogen system based on modified NSGA-II and CRITIC-TOPSIS | |
Liu et al. | A novel hybrid methodology for short-term wind power forecasting based on adaptive neuro-fuzzy inference system | |
Hu et al. | A new clustering approach for scenario reduction in multi-stochastic variable programming | |
Sawle et al. | PV-wind hybrid system: A review with case study | |
Tascikaraoglu et al. | A review of combined approaches for prediction of short-term wind speed and power | |
Yildiz et al. | Household electricity load forecasting using historical smart meter data with clustering and classification techniques | |
Chen et al. | Online 24-h solar power forecasting based on weather type classification using artificial neural network | |
Liu et al. | Wind power plant prediction by using neural networks | |
CN103390116B (zh) | 采用分步方式的光伏电站发电功率预测方法 | |
Michael et al. | Short-term solar irradiance forecasting based on a novel Bayesian optimized deep Long Short-Term Memory neural network | |
Bhaskar et al. | Wind speed forecasting: Present status | |
CN103218674A (zh) | 基于bp神经网络模型的光伏发电系统输出功率预测方法 | |
CN103106544B (zh) | 一种基于t-s型模糊神经网络的光伏发电预测系统 | |
CN102184453A (zh) | 基于模糊神经网络和支持向量机的风电功率组合预测方法 | |
CN103268366A (zh) | 一种适用于分散式风电场的组合风电功率预测方法 | |
CN104951834A (zh) | 基于遗传算法和粒子群集成的lssvm风速预测方法 | |
CN107516145A (zh) | 一种基于加权欧氏距离模式分类的多通道光伏出力预测方法 | |
CN105069521A (zh) | 一种基于加权fcm聚类算法的光伏电站输出功率预测方法 | |
CN102509027B (zh) | 一种基于交叉熵理论的风电功率组合预测方法 | |
CN106295857A (zh) | 一种风电功率超短期预测方法 | |
CN103927460A (zh) | 一种基于rbf的风电场短期风速预测方法 | |
CN104809532A (zh) | 一种光伏系统发电量的预测方法 | |
CN105846425A (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 |