CN104112236A - 风电场发电功率的计算方法 - Google Patents
风电场发电功率的计算方法 Download PDFInfo
- Publication number
- CN104112236A CN104112236A CN201410231145.2A CN201410231145A CN104112236A CN 104112236 A CN104112236 A CN 104112236A CN 201410231145 A CN201410231145 A CN 201410231145A CN 104112236 A CN104112236 A CN 104112236A
- Authority
- CN
- China
- Prior art keywords
- mark post
- wind energy
- turbine set
- energy turbine
- power
- 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 abstract description 6
- 230000003321 amplification Effects 0.000 claims abstract description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000017105 transposition Effects 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 abstract description 3
- 230000005611 electricity Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/04—Automatic control; Regulation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/04—Automatic control; Regulation
- F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
- F03D7/048—Automatic control; Regulation by means of an electrical or electronic controller controlling wind farms
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/028—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/82—Forecasts
- F05B2260/821—Parameter estimation or prediction
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/335—Output power or torque
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Abstract
本发明涉及一种风电场发电功率的计算方法，包括以下步骤：对所有标杆风机的功率序列进行基于经验正交函数分解，得到空间向量；对空间向量还原得到各标杆风机的典型功率序列；按照对应馈线上的开机数对各标杆风机典型功率进行比例放大，得到标杆风机所在馈线总的功率；对所有标杆风机进行叠加得到整个风电场的功率。
Description
技术领域
本发明属于电力系统规划领域，尤其针对大规模风电场（30MW及以上）理论可发电功率的计算方法，以计算风电场在自然状态下的发电量。
背景技术
近年来，我国在一系列政策的推动下，风电产业快速发展。预计到2015年，我国并网运行的风电装机容量达到1亿千瓦。但受电网输送能力和风电消纳能力等因素限制，实际风电并网运行工程中存在明显的限电情况。目前，很多电厂的弃风量是按照风机的理论年发电量笼统计算得出或者采用常规的功率扩大法得出，这与实际的弃风电量有很大的差别。
针对我国特有的风电发展模式，已有的风电发电量估算方法精度较低。由于估算弃风量与实际弃风量误差较大，就无法为风电集群控制、调度运行提供重要的参考信息，不利于风电产业的健康发展。
发明内容
综上所述，确有必要提供一种符合我国风电发展模式，并且能够精确估算风电场发电功率的计算方法。
一种风电场发电功率的计算方法，包括以下步骤：步骤S10，对所有标杆风机的功率序列进行基于经验正交函数分解，得到空间向量 ；步骤S20，对空间向量还原得到各标杆风机的典型功率序列；步骤S30，按照对应馈线上的开机数对各标杆风机典型功率进行比例放大，得到标杆风机所在馈线总的功率；步骤S40，对所有标杆风机进行叠加得到整个风电场的功率。
相对于现有技术，本发明提供的风电场发电功率的计算方法，利用EOF分解挖掘各标杆风机的共性与特性，能够克服现有风电场发电功率计算方法误差过大的不足，精确估算风电场理论发电量，为评估计算限电或检修经济损失提供了重要参考依据。
附图说明
图1为本发明提供的基于EOF分解的风电场发电功率计算方法的流程图。
具体实施方式
下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
请参阅图1，本发明提供的基于经验正交函数分解（Empirical Orthogonal Function，简称EOF）的风电场发电功率计算方法，包括以下步骤：
步骤S10，对所有标杆风机的功率序列进行EOF分解，得到空间向量；
步骤S20，基于EOF分解的空间向量还原得到各标杆风机的典型功率序列；
步骤S30，按照对应馈线上的开机数对各标杆风机典型功率进行比例放大，得到标杆风机所在馈线总的功率；
步骤S40，对所有标杆风机进行叠加得到整个风电场的功率。
在步骤S10中，以风电场标杆风机全年可用的功率数据为基础，设标杆风机数为m，每个标杆风机含有n个点的时间序列。根据EOF分解原理，将空间场抽象为一个m维随机变量X，得到容量为n的样本X1，X2，…，Xn，其中是m维向量，记为：
； (1)
根据EOF分解形式，式（1）可以表示为如下形式：
； (2)
其中，是m维待求空间向量，是与之对应的m维误差向量，是第k个空间向量表示时的权重系数，称为时间系数。
所述空间向量的获取具体包括如下步骤：
步骤S11，获取单位相量，使式：
（3）
的剩余误差平方和E1最小，即有：
（4）
其中，或。且〈〉=0，，进一步可以得到：
（5）
其中，，也即空间场的总方差，它取决于研究的空间场，与、无关。于是有：
（6）
把用和表示出来，则有：
（7）
其中，的第行第列元素具体表示为，即为第号格点序列和第号格点序列的协方差，那么就为协方差矩阵，记：
（8）
进一步有：
（9）
步骤S12，将上述公式的求解转化为求多元函数的条件极值，即在的条件下，求的极小值及相应的值。
采用拉格朗日法，构造函数：
（10）
通过对函数求导，可以得到：
（11）
令其为零向量，可得
（12）
于是所求的就是协方差矩阵的一个特征向量，而拉格朗日乘数是与之对应的特征值。带回式（9），有：
（13）
为达到使最小的目标，应取为最大的一个特征值，记，而便是与该之对应的特征向量。此时有
（14）
步骤S13，计算方差贡献率η。
空间向量的方差贡献为该空间向量加入展开式后使场的总误差方差减少的数量，可得：
（15）
又由于场的总方差为，所以的方差贡献率η为：
（16）
分量的方差贡献率反映了该空间型对原始空间场描述能力，方差贡献率越大，该空间向量还原得到的空间场与原始场越接近，通过这一指标确定典型空间向量V。
在步骤S20中，设方差贡献率η的阀值为95%，从而得到典型空间向量。令该方差贡献率所对应的典型空间向量为，时间系数为，其中，
（17）
则各标杆风机的典型功率序列为
（18）
在步骤S30中，按照对应的馈线上的开机数对典型功率进行比例放大。具体的，设m台标杆风机分别所在馈线上的开机数序列为，则开机数序列的转置为：
（19）
则m台风机所在馈线总的功率：
（20）
在步骤S40中，将公式（20）中向量的所有分量相加就是整个风电场在自然状态下理论发电量P估：
（21）。
相对现有技术，本发明提供的风电场发电功率估算方法，利用EOF分解在挖掘各标杆风机的共性与特性方面的优势，以及分解后某些随机性序列的叠加相消特点，克服了风电场场在自然状态下的理论发电量计算误差较大的现状，提高了所述风电场理论发电量的计算精度，从而为评估计算限电或检修经济损失提供了重要的参考依据。
另外，本领域技术人员还可在本发明精神内作其它变化，当然这些依据本发明精神所作的变化，都应包含在本发明所要求保护的范围内。
Claims (9)
1.一种风电场发电功率的计算方法，包括以下步骤：
步骤S10，对所有标杆风机的功率序列进行基于经验正交函数分解，得到空间向量 ；
步骤S20，对空间向量还原得到各标杆风机的典型功率序列；
步骤S30，按照对应馈线上的开机数对各标杆风机典型功率进行比例放大，得到标杆风机所在馈线总的功率；
步骤S40，对所有标杆风机进行叠加得到整个风电场的功率。
2.如权利要求1所述的风电场发电功率的计算方法，其他特征在于，所有标杆风机的功率序列以风电场标杆风机全年可用的功率数据为基础，设标杆风机数为m，每个标杆风机含有n个点的时间序列。
3.如权利要求2所述的风电场发电功率的计算方法，其特征在于，设是m维向量，代表m个标杆风机：
；
其中，是m维待求空间向量，是与之对应的m维误差向量， 是第k个空间向量表示时的权重系数，称为时间系数。
4.如权利要求3所述的风电场发电功率的计算方法，其特征在于，所述空间向量的获取具体包括：
获取单位相量，使式：的剩余误差平方和E1最小：
，
其中为空间场的总方差；
将上述公式的求解转化为求多元函数的条件极值，即在的条件下，求的极小值及相应的值；
计算方差贡献率η：
；
其中，空间向量的方差贡献为该空间向量加入展开式后使场的总误差方差减少的数量：
。
5.如权利要求4所述的风电场发电功率的计算方法，其特征在于，设设差贡献率η的阀值为95%，得到与该方差贡献率所对应的典型空间向量为，时间系数为，其中，
。
6.如权利要求5所述风电场发电功率的计算方法，其特征在于，各标杆风机的典型功率序列P为：
。
7.如权利要求6所述的风电场发电功率的计算方法，其特征在于，设m台标杆风机分别所在馈线上的开机数序列为，则开机数序列的转置为：
。
8.如权利要求7所述的风电场发电功率的计算方法，其特征在于，m台风机所在馈线总的功率：
。
9.如权利要求8所述的风电场发电功率的计算方法，其特征在于，整个风电场在自然状态下理论发电量P估：
。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410231145.2A CN104112236B (zh) | 2014-05-29 | 2014-05-29 | 风电场发电功率的计算方法 |
| US14/515,530 US9328719B2 (en) | 2014-05-29 | 2014-10-16 | Method of calculating available output power of wind farm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410231145.2A CN104112236B (zh) | 2014-05-29 | 2014-05-29 | 风电场发电功率的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104112236A true CN104112236A (zh) | 2014-10-22 |
| CN104112236B CN104112236B (zh) | 2018-04-27 |
Family
ID=51709017
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410231145.2A Active CN104112236B (zh) | 2014-05-29 | 2014-05-29 | 风电场发电功率的计算方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US9328719B2 (zh) |
| CN (1) | CN104112236B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105978041A (zh) * | 2016-03-23 | 2016-09-28 | 三重型能源装备有限公司 | 一种配置标杆风机的风电场有功功率控制方法 |
| CN107313901A (zh) * | 2017-08-17 | 2017-11-03 | 深圳市易邻益邦科技有限公司 | 一种风电机组发电量损失的评估方法及装置 |
| CN108204341A (zh) * | 2016-12-19 | 2018-06-26 | 北京金风科创风电设备有限公司 | 风电场运行状态的识别方法和装置 |
| CN112347655A (zh) * | 2020-11-17 | 2021-02-09 | 国网青海省电力公司 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020028578A1 (en) * | 2018-07-31 | 2020-02-06 | Alliance For Sustainable Energy, Llc | Distributed reinforcement learning and consensus control of energy systems |
| CN110707704A (zh) * | 2019-10-08 | 2020-01-17 | 河海大学 | 基于gmm及多点线性半不变量法的电-热互联综合能源系统概率潮流分析方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1815247A (zh) * | 2006-03-10 | 2006-08-09 | 中国航天科技集团公司第五研究院第五一四研究所 | 一种电功率测量方法及其测量装置 |
| CN103268573A (zh) * | 2013-05-09 | 2013-08-28 | 国家电网公司 | 一种基于主成分分析的风电场标杆风机选择方法 |
| CN103683274A (zh) * | 2013-07-16 | 2014-03-26 | 国家电网公司 | 区域中长期风电发电量概率预测方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1566772A (en) * | 1977-09-15 | 1980-05-08 | Standard Telephones Cables Ltd | Microstrip antenna radiators |
| US7116010B2 (en) * | 2002-09-17 | 2006-10-03 | Wisconsin Alumni Research Foundation | Control of small distributed energy resources |
| US10324132B2 (en) * | 2010-06-07 | 2019-06-18 | Abb Inc. | Systems and methods for power line event zone identification |
| EP2698870A1 (en) * | 2012-08-14 | 2014-02-19 | Alcatel-Lucent | Antenna feed |
-
2014
- 2014-05-29 CN CN201410231145.2A patent/CN104112236B/zh active Active
- 2014-10-16 US US14/515,530 patent/US9328719B2/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1815247A (zh) * | 2006-03-10 | 2006-08-09 | 中国航天科技集团公司第五研究院第五一四研究所 | 一种电功率测量方法及其测量装置 |
| CN103268573A (zh) * | 2013-05-09 | 2013-08-28 | 国家电网公司 | 一种基于主成分分析的风电场标杆风机选择方法 |
| CN103683274A (zh) * | 2013-07-16 | 2014-03-26 | 国家电网公司 | 区域中长期风电发电量概率预测方法 |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105978041A (zh) * | 2016-03-23 | 2016-09-28 | 三重型能源装备有限公司 | 一种配置标杆风机的风电场有功功率控制方法 |
| CN105978041B (zh) * | 2016-03-23 | 2019-06-18 | 三一重型能源装备有限公司 | 一种配置标杆风机的风电场有功功率控制方法 |
| CN108204341A (zh) * | 2016-12-19 | 2018-06-26 | 北京金风科创风电设备有限公司 | 风电场运行状态的识别方法和装置 |
| CN108204341B (zh) * | 2016-12-19 | 2019-12-10 | 北京金风科创风电设备有限公司 | 风电场运行状态的识别方法和装置 |
| CN107313901A (zh) * | 2017-08-17 | 2017-11-03 | 深圳市易邻益邦科技有限公司 | 一种风电机组发电量损失的评估方法及装置 |
| CN112347655A (zh) * | 2020-11-17 | 2021-02-09 | 国网青海省电力公司 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US9328719B2 (en) | 2016-05-03 |
| CN104112236B (zh) | 2018-04-27 |
| US20150345474A1 (en) | 2015-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104112236A (zh) | 风电场发电功率的计算方法 | |
| CN103020423B (zh) | 基于copula函数获取风电场出力相关特性的方法 | |
| CN105429129B (zh) | 一种考虑网络约束的间歇性能源发电容量置信度评估方法 | |
| CN104682381B (zh) | 大型风电场柔性直流输电系统可靠性计算方法 | |
| EP2505827A2 (en) | Wind power prediction method of single wind turbine generator | |
| CN114143629B (zh) | 一种建筑实时碳排放量监控系统 | |
| CN106643765A (zh) | 一种采集异常维修时长计算方法 | |
| CN104331572A (zh) | 考虑风速和风电机组故障相关性的风电场可靠性建模方法 | |
| CN110659693A (zh) | 基于k近邻分类的配电网快速拓扑识别方法、系统及可读存储介质 | |
| CN107257130A (zh) | 基于区域量测解耦的低压配电网损耗计算方法 | |
| CN102982420A (zh) | 基于多目标决策理想点法的资产全寿命周期评估方法 | |
| CN104915788B (zh) | 一种考虑多风场相关性的电力系统动态经济调度的方法 | |
| CN105529714A (zh) | 基于正态分布组合特征的快速概率潮流计算方法 | |
| CN105279575B (zh) | 基于发电预测的多能源主体分布式博弈优化方法 | |
| CN108462180B (zh) | 一种基于vine copula函数确定概率最优潮流的方法 | |
| CN105939026A (zh) | 基于混合Laplace分布的风电功率波动量概率分布模型建立方法 | |
| CN105354761B (zh) | 一种风电接入电网的安全与效能评估方法及系统 | |
| CN104200001B (zh) | 标杆风机的选取方法 | |
| CN105095674A (zh) | 一种分布式风机出力相关性多场景分析方法 | |
| CN105825435A (zh) | 一种基于时序运行模拟的电网运营效率评估方法 | |
| CN104881818B (zh) | 输电网安全风险综合评价指标体系的建立方法 | |
| CN104680263A (zh) | 基于粒子群算法的电力运输网络拓扑结构设计方法 | |
| CN108092263B (zh) | 一种dg并网选址的方法 | |
| CN104158190A (zh) | 电网潮流安全预测装置及方法 | |
| Fan et al. | A low energy algorithm of wireless sensor networks based on fractal dimension |
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 |