CN105545592B - 一种获取角度测量误差曲线的方法 - Google Patents
一种获取角度测量误差曲线的方法 Download PDFInfo
- Publication number
- CN105545592B CN105545592B CN201510945647.6A CN201510945647A CN105545592B CN 105545592 B CN105545592 B CN 105545592B CN 201510945647 A CN201510945647 A CN 201510945647A CN 105545592 B CN105545592 B CN 105545592B
- Authority
- CN
- China
- Prior art keywords
- measurement error
- anemoclinograph
- wind
- power
- angle measurement
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 239000007787 solid Substances 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005511 kinetic theory Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003672 processing 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/0204—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
-
- 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/32—Wind speeds
-
- 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/321—Wind directions
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种角度测量误差曲线的数据预处理方法,包括以下步骤:按风速值将电磁功率转化为叶轮功率;统计风机的叶轮功率与实测风向偏差角度之间的相互关系;基于最小二乘方法得到一条最小二乘拟合曲线,并选取该曲线最高点对应的实测风向偏差角度作为该风速值下的角度测量误差;采用最小二乘拟合法,对各风速值下的角度测量误差进行最小二乘拟合,形成风速风向仪的角度测量误差函数。本发明将从主控系统中获取风机的“电磁功率”转换为不可被主控系统监测的“叶轮功率”,进而通过“叶轮功率”与实测风向偏差角度进行统计和拟合,极大地改善了数据间的曲线相关性,从而提高了整条角度测量误差曲线δ=f(v)的计算精度。
Description
技术领域
本发明涉及一种风速风向仪,特别是一种风速风向仪的角度测量误差补偿方法。
背景技术
风力发电机将风动能转化为电能,目标是在保持风机承受较低机械载荷的同时尽可能多的发电。要成功实现以上目标,关键是要让风机叶轮精确的对准风向。
如图1-2所示,根据风机动力学理论,当风速恒定且发电机转速低于额定转速时,风机的发电功率与风向偏差角度θ的余弦的三次方成正比。设风向无偏差时,风机叶轮获得的功率为Power1;当风速不变而风向偏差角度为θ时,风机叶轮获得的功率为Power2,则二者满足如下公式:
Power2=Power1×cos3θ
因此,当风向偏差角度θ为15度时,会带来约10%的发电量损失。另外,叶轮的偏离会导致在叶轮乃至整个风机的机械载荷不平衡。这类载荷相比其他载荷会大得多,如果能降低,就能延长风机使用寿命,或者让现有风机带动更大的叶轮。
目前,在大多数风机上,风向偏差角度由安装在机舱上方的风速风向仪决定。但在风机的实际运行过工程中,风速风向仪所测量的风向偏差角度与叶轮处的实际风向偏差角度之间存在误差。
如图3所示,因为风速风向仪测量的是风机机舱尾部的风向偏差角度θ2,而风机主控系统需要的是风机叶轮处的实际风向偏差角度θ1,即两者间的角度测量误差表示为:
δθ=θ2-θ1
大量实验表明,这个误差不是固定的,而是与风速相关的,并将这个随风速变化的误差定义为如下的角度测量误差曲线或角度测量误差函数:
δ=f(v)≈a0+a1·v+a2·v2+…+an·vn
风速风向仪计算角度测量误差曲线的核心是绘制出各风速段下,角度测量误差与功率之间的曲线关系,并在曲线最大值处获得风速风向仪的角度测量误差。为此需要获得风机的风速、风向偏差角度和功率三个参数。
其中风速和风向偏差角度直接来自风速风向仪的实测数据,并经过低通滤波器滤除风速和风向偏差角度中的干扰噪声。
功率参数可以从主控系统中获取,但是从主控系统获取的功率为风机的“电磁功率”,而角度测量误差曲线需要的是“叶轮功率”,为此需要进行数据预处理才能获得更准确的角度测量误差曲线。
关于角度测量误差曲线的数据预处理方法,目前还未见报道。
发明内容
为解决现有技术存在的上述问题,本发明要设计一种角度测量误差曲线的数据预处理方法,以便能够提供更准确数据来提高角度测量误差曲线的补偿精度。
为了获取上述角度测量误差曲线,本发明的技术方案如下:一种风速风向仪的角度测量误差曲线的数据预处理方法,包括以下步骤:
A、风速风向仪以1Hz~10Hz的频率连续采集风场目标风机在1个月内的各项参数,所述的参数包括风速风向仪实测风向偏差角度、风速风向仪实测风速和主控系统实时监测的电磁功率PElectric;
B、风速风向仪从上述采集的风机参数中,按风速从小至大的顺序选取一个风速值vi;下标i为风速值的序号,设i=1;
C、将该风速值vi对应的实测风向偏差角度θ2和电磁功率PElectric进行统计,并用空心圆点描述在角度测量误差统计图中;
D、将电磁功率PElectric转化为叶轮功率PRotor,转化公式如下:
PRotor=PElectric+PGenerator+PGearbox+PInertia
式中:PRotor为叶轮功率,PElectric为电磁功率,PInertia为转动惯量功率,PGearbox为齿轮箱功率流失,PGenerator为发电机功率流失,其中,PInertia、PGearbox及PGenerator的计算方法从风机的传递函数和风机控制理论中计算获得;
E、统计风机的叶轮功率PRotor与风速风向仪实测风向偏差角度θ2之间的相互关系,并用实心圆点描述在角度测量误差统计图中;
F、对角度测量误差统计图中的实心点进行基于最小二乘方法的余弦曲线拟合,得到一条最小二乘拟合曲线,将风机的叶轮功率PRotor与风速风向仪实测风向偏差角度θ2之间的相互关系用一条余弦形式的拟合曲线表示,并选取该曲线最高点对应的风速风向仪实测风向偏差角度作为该风速值vi下的角度测量误差δi;
G、如果该风速值vi不是最大值,则选取下一个风速值,即令i=i+1,返回步骤C;如果该风速值vi是最大值,即获得风速风向仪在该风速值vi下的角度测量误差δi;
H、采用最小二乘拟合法,对步骤G中的各风速值vi下的风速风向仪的角度测量误差δi进行最小二乘拟合,形成风速风向仪的角度测量误差函数或角度测量误差曲线:
δ=f(v)≈a0+a1·v+a2·v2+…+an·vn。
与现有技术相比,本发明具有以下有益效果:
本发明通过数据预处理方法,将从主控系统中获取的风机“电磁功率”转换为不可被主控系统监测的“叶轮功率”,进而通过“叶轮功率”与风速风向仪的实测风向偏差角度进行统计和拟合,极大地改善了数据间的曲线相关性,能够更准确的表述某一风速值vi下,风速风向仪的角度测量误差δi,从而提高了整条角度测量误差曲线δ=f(v)的计算精度。
附图说明
本发明共有附图6张,其中:
图1是风机机舱正对风向示意图。
图2是风机机舱偏航示意图。
图3是实际风速和实测风速示意图。
图4是叶轮功率转化为电磁功率示意图。
图5是角度测量误差统计图。
图6是最小二乘拟合曲线示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步地说明。一种风速风向仪的角度测量误差曲线的数据预处理方法,核心是绘制出各风速值下,风向偏差角度与功率之间的曲线关系,并在曲线最大值处获得风速风向仪的角度测量误差。具体包括以下步骤:
A、风速风向仪以1Hz的频率连续采集河北某风场某风机在1个月内的各项参数,所述的参数包括风速风向仪实测风向偏差角度、风速风向仪实测风速和主控系统实时监测的电磁功率PElectric。
B、风速风向仪从上述采集的风机参数中,按风速从小至大的顺序选取一个风速值vi;下标i为风速值的序号,设i=1,v1=10m/s;
C、将该风速值vi对应的实测风向偏差角度θ2和电磁功率PElectric进行统计,并用空心圆点描述在角度测量误差统计图中,如图5所示;由图5可知,风机的电磁功率PElectric与风速风向仪实测风向偏差角度θ2之间没有很好的相关性;
D、如图4所示,将电磁功率PElectric转化为叶轮功率PRotor,转化公式如下:
PRotor=PElectric+PGenerator+PGearbox+PInertia
式中:PRotor为叶轮功率,PElectric为电磁功率,PInertia为转动惯量功率,PGearbox为齿轮箱功率流失,PGenerator为发电机功率流失,其中,PInertia、PGearbox及PGenerator的计算方法从风机的传递函数和风机控制理论中计算获得;
E、统计风机的叶轮功率PRotor与风速风向仪实测风向偏差角度θ2之间的相互关系,并用实心圆点描述在角度测量误差统计图中,如图5所示;由图5可知,实心点的分布有明显的走势变化,其曲线类似于余弦形式;
F、对角度测量误差统计图中的实心点进行基于最小二乘方法的余弦曲线拟合,得到一条最小二乘拟合曲线,将风机的叶轮功率PRotor与风速风向仪实测风向偏差角度θ2之间的相互关系用一条余弦形式的拟合曲线表示,如图6所示,并很容易的发现该曲线的最高点并没有出现在实测风向偏差角度θ2=0处,而是当风速风向仪实测风向偏差角度θ2=2°时,得到曲线最高点。即说明河北该风场该风机在风速值v1=10m/s时,风速风向仪实测风向偏差角度具有的角度测量误差为δ1=2°。
G、如果该风速值vi不是最大值,则选取下一个风速值,即令i=i+1,返回步骤C;如果该风速值vi是最大值,即获得风速风向仪在该风速值vi下的角度测量误差δi;
H、采用最小二乘拟合法,对步骤G中的各风速值vi下的风速风向仪的角度测量误差δi进行最小二乘拟合,形成风速风向仪的角度测量误差函数或角度测量误差曲线:
δ=f(v)≈a0+a1·v+a2·v2+…+an·vn。
Claims (1)
1.一种角度测量误差曲线的数据预处理方法,其特征在于:包括以下步骤:
A、风速风向仪以1Hz~10Hz的频率连续采集风场目标风机在1个月内的各项参数,所述的参数包括风速风向仪实测风向偏差角度、风速风向仪实测风速和主控系统实时监测的电磁功率PElectric;
B、风速风向仪从上述采集的风机参数中,按风速从小至大的顺序选取一个风速值vi;下标i为风速值的序号,设i=1;
C、将该风速值vi对应的实测风向偏差角度θ2和电磁功率PElectric进行统计,并用空心圆点描述在角度测量误差统计图中;
D、将电磁功率PElectric转化为叶轮功率PRotor,转化公式如下:
PRotor=PElectric+PGenerator+PGearbox+PInertia
式中:PRotor为叶轮功率,PElectric为电磁功率,PInertia为转动惯量功率,PGearbox为齿轮箱功率流失,PGenerator为发电机功率流失,其中,PInertia、PGearbox及PGenerator的计算方法从风机的传递函数和风机控制理论中计算获得;
E、统计风机的叶轮功率PRotor与风速风向仪实测风向偏差角度θ2之间的相互关系,并用实心圆点描述在角度测量误差统计图中;
F、对角度测量误差统计图中的实心圆点进行基于最小二乘方法的余弦曲线拟合,得到一条最小二乘拟合曲线,将风机的叶轮功率PRotor与风速风向仪实测风向偏差角度θ2之间的相互关系用一条余弦形式的拟合曲线表示,并选取该曲线最高点对应的风速风向仪实测风向偏差角度作为该风速值vi下的角度测量误差δi;
G、如果该风速值vi不是最大值,则选取下一个风速值,即令i=i+1,返回步骤C;如果该风速值vi是最大值,即获得风速风向仪在该风速值vi下的角度测量误差δi;
H、采用最小二乘拟合法,对步骤G中的各风速值vi下的风速风向仪的角度测量误差δi进行最小二乘拟合,形成风速风向仪的角度测量误差函数或角度测量误差曲线:
δ=f(v)≈a0+a1·v+a2·v2+…+an·vn。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510945647.6A CN105545592B (zh) | 2015-12-16 | 2015-12-16 | 一种获取角度测量误差曲线的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510945647.6A CN105545592B (zh) | 2015-12-16 | 2015-12-16 | 一种获取角度测量误差曲线的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105545592A CN105545592A (zh) | 2016-05-04 |
CN105545592B true CN105545592B (zh) | 2018-09-25 |
Family
ID=55825047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510945647.6A Active CN105545592B (zh) | 2015-12-16 | 2015-12-16 | 一种获取角度测量误差曲线的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105545592B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112648139B (zh) * | 2020-11-12 | 2022-03-04 | 北京金风慧能技术有限公司 | 风力发电机组对风不正校正方法、装置及控制器 |
CN115330092B (zh) * | 2022-10-13 | 2023-03-24 | 山东东盛澜渔业有限公司 | 基于人工智能的海洋牧场再生能源供能控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498282A (zh) * | 2008-02-01 | 2009-08-05 | 北京能高自动化技术有限公司 | 大型风力发电机组偏航控制方法 |
CN101568721A (zh) * | 2007-05-25 | 2009-10-28 | 三菱重工业株式会社 | 风力发电装置、风力发电系统以及风力发电装置的发电控制方法 |
JP5323133B2 (ja) * | 2011-06-03 | 2013-10-23 | 株式会社東芝 | 風力発電システムの制御方法 |
CN103514353A (zh) * | 2012-06-29 | 2014-01-15 | 国际商业机器公司 | 新增风机的功率预测方法和系统 |
CN103630705A (zh) * | 2013-01-16 | 2014-03-12 | 中国科学院电子学研究所 | 一种固态二维风速风向测量仪及其测量方法 |
CN104314757A (zh) * | 2014-10-15 | 2015-01-28 | 国电联合动力技术有限公司 | 一种风力发电机组偏航控制方法及系统 |
CN104481804A (zh) * | 2014-12-05 | 2015-04-01 | 北京金风科创风电设备有限公司 | 风力发电机组对风矫正控制方法、装置和系统 |
-
2015
- 2015-12-16 CN CN201510945647.6A patent/CN105545592B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101568721A (zh) * | 2007-05-25 | 2009-10-28 | 三菱重工业株式会社 | 风力发电装置、风力发电系统以及风力发电装置的发电控制方法 |
CN101498282A (zh) * | 2008-02-01 | 2009-08-05 | 北京能高自动化技术有限公司 | 大型风力发电机组偏航控制方法 |
JP5323133B2 (ja) * | 2011-06-03 | 2013-10-23 | 株式会社東芝 | 風力発電システムの制御方法 |
CN103514353A (zh) * | 2012-06-29 | 2014-01-15 | 国际商业机器公司 | 新增风机的功率预测方法和系统 |
CN103630705A (zh) * | 2013-01-16 | 2014-03-12 | 中国科学院电子学研究所 | 一种固态二维风速风向测量仪及其测量方法 |
CN104314757A (zh) * | 2014-10-15 | 2015-01-28 | 国电联合动力技术有限公司 | 一种风力发电机组偏航控制方法及系统 |
CN104481804A (zh) * | 2014-12-05 | 2015-04-01 | 北京金风科创风电设备有限公司 | 风力发电机组对风矫正控制方法、装置和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105545592A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DK1906192T3 (en) | Apparatus for evaluating sensors and / or for controlling the operation of an apparatus which includes a sensor | |
EP2169218B1 (en) | System and method for estimating wind condition for wind turbines | |
CN203685475U (zh) | 风力涡轮机控制系统及风力涡轮机系统 | |
KR101476986B1 (ko) | 풍력 터빈용 제어 장치 | |
CN103758700B (zh) | 一种校准风机对风偏差的方法 | |
US10018177B2 (en) | Control system and method for mitigating rotor imbalance on a wind turbine | |
CN105894391A (zh) | 基于scada运行数据提取的风电机组转矩控制性能评估方法 | |
EP2520793A2 (en) | Methods and apparatus for controlling wind turbine thrust | |
EP3111084A1 (en) | Method and system for improving energy capture efficiency from an energy capture device | |
CN105464903B (zh) | 一种准确获取角度测量误差曲线的方法 | |
CN102859184A (zh) | 用于风力涡轮机的监测方法 | |
EP2466133A2 (en) | System and method for measuring shaft deflection in a wind turbine | |
CN105569922A (zh) | 一种基于风速影响的风速风向仪的角度测量误差补偿方法 | |
CN105545592B (zh) | 一种获取角度测量误差曲线的方法 | |
CN105041571B (zh) | 预测风速风向的智能控制系统及其控制方法 | |
CN106951977B (zh) | 一种基于尾流效应的风速预测模型的构建方法 | |
CN105138845B (zh) | 获得风力发电机风速值的方法 | |
CN105569921B (zh) | 加入主控系统数据传递误差修正的角度测量误差补偿方法 | |
US11274656B2 (en) | Method of determining a power curve of a wind turbine | |
CN108691727B (zh) | 一种风力机导流罩 | |
CN105512416A (zh) | 一种风机尾流对角度测量误差影响关系的获取方法 | |
Pedersen et al. | Using wind speed from a blade-mounted flow sensor for power and load assessment on modern wind turbines | |
CN105512820A (zh) | 风电场风能资源监测与机组运行状态评估方法 | |
EP3677772B1 (en) | Method of operating wind turbine based on maximum thrust limit | |
CN207336561U (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
PP01 | Preservation of patent right |
Effective date of registration: 20190422 Granted publication date: 20180925 |
|
PP01 | Preservation of patent right | ||
PD01 | Discharge of preservation of patent |
Date of cancellation: 20220422 Granted publication date: 20180925 |
|
PD01 | Discharge of preservation of patent |