CN103366225B - 风功率预测误差识别方法 - Google Patents
风功率预测误差识别方法 Download PDFInfo
- Publication number
- CN103366225B CN103366225B CN201310295998.8A CN201310295998A CN103366225B CN 103366225 B CN103366225 B CN 103366225B CN 201310295998 A CN201310295998 A CN 201310295998A CN 103366225 B CN103366225 B CN 103366225B
- Authority
- CN
- China
- Prior art keywords
- wind power
- wind speed
- forecasting
- data
- error
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000013507 mapping Methods 0.000 claims abstract description 29
- 238000013528 artificial neural network Methods 0.000 claims description 19
- 238000000342 Monte Carlo simulation Methods 0.000 description 6
- 230000005611 electricity Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种风功率预测误差识别方法,包括步骤:S1.获取风电场的风速数据、风速预测误差和风功率实测数据;S2.将步骤S1获取到的所有数据进行预处理以剔除错误的数据;S3.根据步骤S2得到的预处理后的风速数据和风功率实测数据得到风功率与风速之间的映射关系;S4.根据步骤S2得到的预处理后的风速预测误差数据得到风速预测误差的分布特性;S5.根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系以及风速预测值与风功率预测误差最小似然值的关系。本发明能够减小日前发电计划误差。
Description
技术领域
本发明涉及新能源发电与控制领域,尤其涉及一种风功率预测误差识别方法。
背景技术
进入新世纪以来,化石能源短缺和环境污染的状况越来越严重,促使电力行业寻找开发可再生的清洁能源替代现有的化学能源,优化能源结构,其中,风电作为一种大规模存在的可再生的清洁能源开始逐渐受到人们的重视。另一方面,风电的一次能源风能具有很大的波动性和间歇性,会给电力系统造成较大的干扰,因此需要对风功率进行一定的预测,将风电纳入常规发电计划,才能更好地管理和利用风电。根据国家能源局的要求,风电场应该上报日前预测数据(未来24小时,96个点)和实时预测(未来4小时,15分钟一个点)。但是实际中,风电场上报的日前预测数据误差很大,用于日前发电计划计算时可能产生较大误差。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:提供一种简单实用的基于蒙特卡洛仿真的风功率预测误差识别方法,用于解析风功率预测误差的统计特性,以便对日前发电计划进行调整,从而减小日前发电计划误差。
(二)技术方案
为解决上述问题,本发明提供了一种风功率预测误差识别方法,该方法包括步骤:
S1.获取风电场的风速数据、风速预测误差和风功率实测数据;
S2.将步骤S1获取到的所有数据进行预处理以剔除错误的数据;
S3.根据步骤S2得到的预处理后的风速数据和风功率实测数据得到风功率与风速之间的映射关系;
S4.根据步骤S2得到的预处理后的风速预测误差数据得到风速预测误差的分布特性;
S5.根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系以及风速预测值与风功率预测误差最小似然值的关系。
优选地,步骤S2包括:
S2.1对于任意时刻,若该时刻的风速数据或风功率实测数据没有数据,或显示为NaN,则剔除该时刻风速数据和风功率实测数据;
S2.2若该时刻的风速数据不在预先设定的第一范围内,或者该时刻的风功率实测数据不在预先设定的第二范围内,则剔除该时刻风速数据和风功率实测数据。
优选地,根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系包括:
S5.1记s=0,在该第一范围内随机生成一个风速预测值;
S5.2根据步骤S4得到的风速预测误差分布N(μ,σ2),随机生成一个服从该风速预测误差分布的实数
S5.3根据和步骤S3得到的映射关系计算得到步骤5.1生成的风速预测值对应的风功率误差,即令s=s+δp;
S5.4将S5.2~S5.3循环执行M次;
S5.5计算该风功率预测误差平均值E(δp)=s/M;
S5.6将步骤S5.1~S5.5执行L次,得到L组对应的数据,以为BP人工神经网络的输入,E(δp)为BP人工神经网络的输出,训练得到风速预测值与风功率预测误差平均值E(δp)之间的映射关系,记为
优选地,根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差最小似然值的关系包括:
S5.a记s=0,在该第一范围内随机生成一个风速预测值;
S5.b根据步骤S4得到的风速预测误差分布N(μ,σ2),随机生成M个服从该分布的实数i=1,2,…M;
S5.c根据和步骤S3得到的映射关系计算得到步骤5.a生成的风速预测值对应的风功率误差,即
S5.d取M'=[βM],[*]符号表示取整,该β为预先设定的概率值;
S5.e将{δp1,δp2,…,δpM}从小到大排序,第M'个记为风功率预测误差最小似然值
S5.f将步骤S5.a~S5.e执行L次,得到L组对应的数据,以为BP人工神经网络的输入,为BP人工神经网络的输出,训练得到风速预测值与风功率预测误差最小似然值之间的映射关系,记为
优选地,该M大于或等于3000,该L大于或等于2000。
(三)有益效果
本发明通过蒙特拉洛仿真方法和BP人工神经网络方法,对风功率预测误差数据进行识别,得到风功率预测误差的统计分布规律,从而能够对日前发电计划进行调整,进而减小日前发电计划误差。
附图说明
图1为依照本发明一种实施方式的风功率预测误差识别方法的流程图。
具体实施方式
本发明提出的风功率预测误差识别方法,结合附图及实施例详细说明如下。
本发明的核心思想为:通过蒙特拉洛仿真方法和BP人工神经网络方法,对风功率预测误差数据进行识别,得到风速预测值与风功率预测误差平均值的关系以及风速预测值与风功率预测误差最小似然值的关系,进而得到风速预测值与风功率预测误差的统计分布规律,从而能够对日前发电计划进行调整,进而减小日前发电计划误差。
如图1所示,依照本发明一种实施方式的基于蒙特卡洛仿真的风功率预测误差识别方法包括步骤:
S1.获取风电场的风速数据、风速预测误差和风功率实测数据。
例如,可以获取该风电场最近一个自然年内风速数据、风速预测误差和风功率实测数据,分别记为和Porigin(T),其中风速数据是矢量,不光包含风速大小信息,还包括风速方向信息,T=1,……,N表示一年内风速数据与风功率实测数据的时标数列(其中,N表示时标的总个数,N=96×D,D表示一年内的总天数,根据实际情况可取值365或366)。按照中国能源局的规定,风电场上报的风功率实测数据为15min一个数据点,因此一个自然天有96个数据点,一个自然年有96×D个数据点。
S2.将步骤S1获取到的所有数据进行预处理以剔除错误的数据。
可以将步骤S1得到的数据进行预处理,剔除数据采集、传输、保存过程中产生错误的数据,输出预处理之后的风速、风速预测误差和风功率实测数据,分别记为和P(t),其中t=1,……,n表示预处理之后合格的风速数据与风功率实测数据的时标数列。
S3.根据步骤S2得到的预处理后的风速数据和风功率实测数据得到风功率与风速之间的映射关系。
根据步骤S2得到的预处理后的风速数据和风功率实测数据,形成BP人工神经网络运算的输入、输出数据组。其中,每一组输入数据包括两个数据单元,即某一时刻的风速大小数据V(t)和风速方向数据每一组输出数据包括一个数据单元,是同一时刻的风功率实测数据P(t)。随机选取数据组的2/3作为训练集,剩下的1/3作为检验集,利用BP人工神经网络训练得到任意时刻风电场风功率与风速之间的映射关系,记为
S4.根据步骤S2得到的预处理后的风速预测误差数据得到风速预测误差的分布特性。
可以根据步骤S2得到的预处理后的风速预测误差数据,计算得到风速预测误差数据的期望μ和方差σ2,拟合得到风速预测误差的分布特性,记为N(μ,σ2)。
S5.根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系以及风速预测值与风功率预测误差最小似然值的关系。
例如,根据步骤S3得到的映射关系和步骤S4得到的风速预测误差的分布特性,利用蒙特卡洛仿真方法和BP人工神经网络方法得到风速预测值与风功率预测误差平均值的关系,记为根据步骤S3得到的映射关系和步骤S4.得到的风速预测误差分布特性,利用蒙特卡洛仿真方法和BP人工神经网络方法得到风速预测值与风功率预测误差最小似然值的关系,记为其中β为事先设定好的概率值。
本发明通过蒙特拉洛仿真方法和BP人工神经网络方法,对风功率预测误差数据进行识别,得到风速预测值与风功率预测误差平均值的关系以及风速预测值与风功率预测误差最小似然值的关系,进而得到风速预测值与风功率预测误差的统计分布规律,从而能够对日前发电计划进行调整,进而减小日前发电计划误差。
优选地,步骤S2可以具体包括:
S2.1对于任意时刻,若该时刻的风速数据或风功率实测数据没有数据,或显示为NaN,则剔除该时刻风速数据和风功率实测数据。
例如,对于任意时刻,若该时刻的风速数据或风功率实测数据没有数据,或显示为NaN(NotaNumber,没有数据),则剔除该时刻风速数据和风功率实测数据。
S2.2若该时刻的风速数据不在预先设定的第一范围内,或者该时刻的风功率实测数据不在预先设定的第二范围内,则剔除该时刻风速数据和风功率实测数据。
优选地,根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系具体包括:
S5.1记s=0,在所述第一范围内随机生成一个风速预测值。
S5.2根据步骤S4得到的风速预测误差分布N(μ,σ2),随机生成一个服从所述风速预测误差分布的实数
S5.3根据和步骤S3得到的映射关系计算得到步骤5.1生成的风速预测值对应的风功率误差,即令s=s+δp。
S5.4将S5.2~S5.3循环执行M次。
优选地,M大于或等于3000,例如,M可以为3000或者3500。
S5.5计算所述风功率预测误差平均值E(δp)=s/M。
S5.6将步骤S5.1~S5.5执行L次,得到L组对应的数据,以为BP人工神经网络的输入,E(δp)为BP人工神经网络的输出,训练得到风速预测值与风功率预测误差平均值E(δp)之间的映射关系,记为
优选地,L大于或等于2000,例如L可以为2000或者2500。
优选地,根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差最小似然值的关系具体包括:
S5.a记s=0,在所述第一范围内随机生成一个风速预测值。
S5.b根据步骤S4得到的风速预测误差分布N(μ,σ2),随机生成M个服从该分布的实数i=1,2,…M。
S5.c根据和步骤S3得到的映射关系计算得到步骤5.a生成的风速预测值对应的风功率误差,即
S5.d取M'=[βM],[*]符号表示取整,所述β为预先设定的概率值。
优选地,M大于或等于3000,例如,M可以为3000或者3500。
S5.e将{δp1,δp2,…,δpM}从小到大排序,第M'个记为风功率预测误差最小似然值
S5.f将步骤S5.a~S5.e执行L次,得到L组对应的数据,以为BP人工神经网络的输入,为BP人工神经网络的输出,训练得到风速预测值与风功率预测误差最小似然值之间的映射关系,记为
优选地,L大于或等于2000,例如L可以为2000或者2500。
本发明的基于蒙特卡洛仿真的风功率预测误差识别方法提高电网的控制技术水平,克服现有风功率预测手段不足、预测精度相对较差等问题,提高了电网的调度和运行质量,前瞻电网未来的发展,可很大限度地改善含高风电渗透率的绿色电力系统运行的经济性和电能质量;此外,本发明提供的方法计算速度快,能够满足在线计算的要求。
本发明的基于蒙特卡洛仿真的风功率预测误差识别方法可以用于我国各省级电力系统和地区风电场的经济调度自动化系统之中,能够提高风电场出力预测的安全性和准确性,具有重大的经济和社会效益。
以上实施方式仅用于说明本发明,而并非对本发明的限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行各种组合、修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种风功率预测误差识别方法,其特征在于,该方法包括步骤:
S1.获取风电场的风速数据、风速预测误差和风功率实测数据;
S2.将步骤S1获取到的所有数据进行预处理以剔除错误的数据;
S3.根据步骤S2得到的预处理后的风速数据和风功率实测数据得到风功率与风速之间的映射关系;
S4.根据步骤S2得到的预处理后的风速预测误差数据得到风速预测误差的分布特性;
S5.根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系以及风速预测值与风功率预测误差最小似然值的关系;
其中,步骤S2包括:
S2.1对于任意时刻,若该时刻的风速数据或风功率实测数据没有数据,或显示为NaN,则剔除该时刻风速数据和风功率实测数据;
S2.2若该时刻的风速数据不在预先设定的第一范围内,或者该时刻的风功率实测数据不在预先设定的第二范围内,则剔除该时刻风速数据和风功率实测数据;
其中,根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差平均值的关系包括:
S5.1记s=0,在所述第一范围内随机生成一个风速预测值;
S5.2根据步骤S4得到的风速预测误差分布N(μ,σ2),随机生成一个服从所述风速预测误差分布的实数
S5.3根据和步骤S3得到的映射关系fNN计算得到步骤5.1生成的风速预测值对应的风功率误差,即令s=s+δp;
S5.4将S5.2~S5.3循环执行M次;
S5.5计算所述风功率预测误差平均值E(δp)=s/M;
S5.6将步骤S5.1~S5.5执行L次,得到L组对应的数据,以为BP人工神经网络的输入,E(δp)为BP人工神经网络的输出,训练得到风速预测值与风功率预测误差平均值E(δp)之间的映射关系,记为
其中,根据步骤S3得到的风功率与风速之间的映射关系以及步骤S4得到的风速预测误差的分布特性得到风速预测值与风功率预测误差最小似然值的关系包括:
S5.a记s=0,在所述第一范围内随机生成一个风速预测值;
S5.b根据步骤S4得到的风速预测误差分布N(μ,σ2),随机生成M个服从该分布的实数i=1,2,…M;
S5.c根据和步骤S3得到的映射关系fNN计算得到步骤5.a生成的风速预测值对应的风功率误差,即
S5.d取M'=[βM],[*]符号表示取整,所述β为预先设定的概率值;
S5.e将{δp1,δp2,…,δpM}从小到大排序,第M'个记为风功率预测误差最小似然值δp为{δp1,δp2,…,δpM}中的任意一个元素;
S5.f将步骤S5.a~S5.e执行L次,得到L组对应的数据,以为BP人工神经网络的输入,为BP人工神经网络的输出,训练得到风速预测值与风功率预测误差最小似然值之间的映射关系,记为
2.如权利要求1所述的风功率预测误差识别方法,其特征在于,所述M大于或等于3000,所述L大于或等于2000。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310295998.8A CN103366225B (zh) | 2013-07-15 | 2013-07-15 | 风功率预测误差识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310295998.8A CN103366225B (zh) | 2013-07-15 | 2013-07-15 | 风功率预测误差识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103366225A CN103366225A (zh) | 2013-10-23 |
CN103366225B true CN103366225B (zh) | 2016-03-30 |
Family
ID=49367518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310295998.8A Active CN103366225B (zh) | 2013-07-15 | 2013-07-15 | 风功率预测误差识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103366225B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103925155B (zh) * | 2014-04-09 | 2016-10-05 | 中国水利水电科学研究院 | 一种风电机组输出功率异常的自适应检测方法 |
CN104484704A (zh) * | 2014-12-22 | 2015-04-01 | 国家电网公司 | 一种基于bp神经网络的风电场理论功率确定方法 |
EP3262732B1 (en) | 2015-02-23 | 2020-06-24 | Okinawa Institute of Science and Technology School Corporation | System and method of determining forecast error for renewable energy fluctuations |
CN105719002A (zh) * | 2016-01-18 | 2016-06-29 | 重庆大学 | 一种基于组合预测的风电机组状态参数异常辨识方法 |
CN106056254B (zh) * | 2016-06-15 | 2020-10-02 | 中国能源建设集团甘肃省电力设计院有限公司 | 一种计及出力水平影响的风电预测误差模拟方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102479339A (zh) * | 2010-11-24 | 2012-05-30 | 香港理工大学 | 基于混合神经网络的风电场短期风速预测方法和系统 |
-
2013
- 2013-07-15 CN CN201310295998.8A patent/CN103366225B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102479339A (zh) * | 2010-11-24 | 2012-05-30 | 香港理工大学 | 基于混合神经网络的风电场短期风速预测方法和系统 |
Non-Patent Citations (3)
Title |
---|
一种改进的风电功率预测误差分布模型;刘斌 等;《华东电力》;20120224;第40卷(第2期);286-291 * |
含风电场的电力系统优化调度研究;张成相;《中国优秀硕士学位论文全文数据库工程科技II辑》;20130615(第06期);C042-596 * |
基于预测误差分布特性统计分析的概率性短期负荷预测;杨文佳 等;《电力系统自动化》;20061010;第30卷(第19期);47-52 * |
Also Published As
Publication number | Publication date |
---|---|
CN103366225A (zh) | 2013-10-23 |
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 | |
CN103683326B (zh) | 一种区域电网风电多点接入最佳接纳能力的计算方法 | |
CN103762589B (zh) | 一种电网中新能源容量配比分层优化方法 | |
CN103366225B (zh) | 风功率预测误差识别方法 | |
Zhu et al. | Renewable energy consumption technology under energy internet environment | |
Ma et al. | Economic analysis of real-time large scale PEVs network power flow control algorithm with the consideration of V2G services | |
CN109103926A (zh) | 基于多辐照特性年气象场景的光伏发电接纳能力计算方法 | |
CN106816886B (zh) | 一种大规模风电并网系统调峰需求确定方法 | |
CN102496968A (zh) | 间歇式能源与常规能源协调调度模式下发电计划优化方法 | |
CN104779611A (zh) | 基于集中式和分布式双层优化策略的微电网经济调度方法 | |
CN103049798A (zh) | 一种应用于光伏发电系统的短期发电功率预测方法 | |
CN108039737A (zh) | 一种源网荷协调运行模拟系统 | |
CN102129511A (zh) | 一种基于matlab的风电场短期风速预测系统 | |
Lu et al. | Initial exploration of wind farm cluster hierarchical coordinated dispatch based on virtual power generator concept | |
CN103997039A (zh) | 基于概率区间预测的计及风电接纳的旋转备用区间预测方法 | |
CN103345585A (zh) | 基于支持向量机的风功率预测校正方法及系统 | |
CN102738792B (zh) | 一种风电功率预测方法 | |
CN104143839A (zh) | 基于功率预测的风电场集群限出力有功功率分配方法 | |
Zhang et al. | Assessing the integration potential of new energy in river basin clean energy corridors considering energy-power coupled complementary operation modes | |
CN108767855B (zh) | 一种时序持续混合的电力系统随机生产模拟方法 | |
CN105279575B (zh) | 基于发电预测的多能源主体分布式博弈优化方法 | |
CN104182808A (zh) | 一种基于等比例限电的新能源场站发电计划制定方法 | |
Hong-Tao et al. | Local Consumption model of wind and photovoltaic power based on demand side response | |
CN103530508B (zh) | 一种建立风速功率转换概率模型的方法 | |
CN105633950A (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 |