CN109753759B - 一种基于等效功率的风轮等效风速计算方法 - Google Patents
一种基于等效功率的风轮等效风速计算方法 Download PDFInfo
- Publication number
- CN109753759B CN109753759B CN201910136172.4A CN201910136172A CN109753759B CN 109753759 B CN109753759 B CN 109753759B CN 201910136172 A CN201910136172 A CN 201910136172A CN 109753759 B CN109753759 B CN 109753759B
- Authority
- CN
- China
- Prior art keywords
- wind
- wind wheel
- equivalent
- wind speed
- speed
- 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
Images
Landscapes
- Wind Motors (AREA)
Abstract
本发明属于风电机组控制技术领域,尤其涉及一种基于等效功率的风轮等效风速计算方法,包括:考虑风切变效应对风轮扫略平面内风速的影响来确定叶片某方位角处的风速计算式;采用三阶泰勒级数展开简化风速计算式,并利用微元面积法,求取风轮空间平均风速其对应的叶尖速比;定义风速产生的气动转矩与风轮等效风速、叶片气动转矩影响系数的关系;求取风轮机械转矩在风轮空间平均风速所对应的叶尖速比处的偏微分线性化,综合得到风轮转矩公式;通过风轮瞬时机械功率、通过风轮平均机械功率和基于等效风速的风轮机械功率求得基于等效功率的风轮瞬时等效风速和风轮平均等效风速。本发明为风电机组设计、控制及输出功率特性分析等领域优化了理论基础。
Description
技术领域
本发明属于风电机组控制技术领域,尤其涉及一种基于等效功率的风轮等效风速计算方法。
背景技术
低风速地区和海上风电是未来风电持续规模化开发的重要增量市场。低风速地区和海上风电机组具有高塔筒、长叶片的特点,风速在整个风轮扫略平面内波动明显增加,仅用轮毂高度风速进行风电机组设计、控制及输出功率特性分析将会带来更大的误差,应该采用能够代表整个风轮扫略平面内风速效应的等效风速。但是,现有的基于等效转矩的风轮等效风速计算方法在计算过程中假设风轮是静止的,而风电机组实际运行过程中,风轮是不断转动的,且随着风速的变化,风轮转速、风能利用系数、有效转矩及有效功率也随之变化,现有等效风速计算方法并不符合风电机组实际运行状况。
发明内容
针对上述技术问题,本发明提出了一种基于等效功率的风轮等效风速计算方法,包括:
步骤1:考虑风切变效应对风轮扫略平面内风速的影响来确定叶片某方位角处的风速计算式;
步骤2:采用三阶泰勒级数展开简化风速计算式,并利用微元面积法,求取风轮空间平均风速其对应的叶尖速比;
步骤3:定义风速产生的气动转矩与风轮等效风速、叶片气动转矩影响系数的关系;
步骤4:求取风轮机械转矩在风轮空间平均风速所对应的叶尖速比处的偏微分线性化,综合得到风轮转矩公式;
步骤5:通过风轮瞬时机械功率和基于等效风速的风轮机械功率求得基于等效功率的风轮瞬时等效风速;
步骤6:通过风轮平均机械功率和基于等效风速的风轮机械功率求得基于等效功率的风轮平均等效风速。
所述基于等效功率的风轮瞬时等效风速为Veq(θ);
式中,R为风轮半径,H为轮毂高度,V0为风轮空间平均风速,Ω为风轮转速,θ为叶片方位角,CP(λ,β)为轮毂高度风速对应风能利用系数,CP(λ0,β0)为风轮空间平均风速对应风能利用系数,λ0为风轮空间平均风速V0对应叶尖速比,α为风切变系数,VH为来流风速。
所述基于等效功率的风轮平均等效风速为Veq;
式中,R为风轮半径,H为轮毂高度,Ω为风轮转速,θ为叶片方位角,CP(λ,β)为轮毂高度风速对应风能利用系数,CP(λ0,β0)为风轮空间平均风速对应风能利用系数,λ0为风轮空间平均风速V0对应叶尖速比,α为风切变系数,VH为来流风速。
本发明的有益效果:
本发明基于功率等效原则,综合考虑风电机组不同控制阶段的风轮转速和风能利用系数对输出功率的影响,构建了基于等效功率的风轮等效风速数学模型,并利用风电机组设计和运行数据,对比分析基于等效功率模型和基于等效转矩模型所得等效风速随风切变、叶片方位角及来流风速等不同参数的变化规律,本发明所提模型能够有效地反映风轮等效风速在风电机组不同控制阶段的变化规律,更符合风电机组实际运行状况,为低风速地区和海上风电机组设计、控制及输出功率特性分析等领域的研究优化了理论基础。
附图说明
图1是风轮等效风速随方位角变化规律图;
图2是风轮平均等效风速随风切变变化规律图;
图3是风轮平均等效风速随来流风速变化规律图。
具体实施方式
下面结合附图,对实施例作详细说明。
本发明涉及一种基于等效功率的风轮等效风速计算方法,数学模型构建具体步骤如下:
考虑风切变效应对风轮扫略平面内风速的影响,第i个叶片,半径为r的叶素,方位角为θi处的风速可表示为:
式(1)中,H为轮毂高度,VH为轮毂高度风速,α为风切变系数。
对于三叶片风电机组,风切变效应引起气动转矩3p振荡,三阶泰勒级数展开式能够很好的反映风切变效应,之后各项影响较小,可以忽略,风速采用3阶泰勒级数展开可简化为:
利用微元面积法,求取风轮空间平均风速V0及其对应的叶尖速比λ0,如式(3)和式(4)所示。
对于三叶片风电机组,风速V(r,θ)产生的气动转矩可定义为:
式(5)中,Veq为方位角为θ时的风轮等效风速,ψ(r)为叶片气动转矩影响系数,定义ψ(r)=kr,则:
定义ΔV=Veq(θ)-V0,并将机械转矩计算公式(7)在(V0,λ0)处求偏微分线性化(λ0为风轮空间平均风速V0对应的叶尖速比)如式(8)所示。
结合式(6)和式(8)可得:
将式(9)和(10)代入式(5)可得转矩公式为:
相对于叶片长度,轮毂半径可以忽略,即r0=0,并将式(2)、式(3)以及式(12)代入(11)可得:
风轮吸收的瞬时机械功率为:
风轮瞬时功率随着方位角的变化而周期性变化,在0到2π内求取平均机械功率,如式(15)所示。
基于等效风速的风轮机械功率可表示为:
由式(14)和式(16)求得风轮瞬时等效风速为:
由式(15)和式(16)求得风轮平均等效风速为:
下面以风电机组设计和运行数据,对比分析基于等效功率模型和基于等效转矩模型所得等效风速随风切变、叶片方位角及来流风速等不同参数的变化规律,验证本发明具有的有益效果。
基于等效转矩的风轮等效风速计算模型所得瞬时等效风速:
基于等效转矩的风轮等效风速计算模型所得平均等效风速:
风电机组额定功率为2MW,额定风速为11m/s,额定转速为18转/分,塔架高度为61.5m,轮毂直径为80m,利用Bladed软件仿真计算风电机组模型的静态功率曲线数据以及不同桨距角下对应的性能系数数据。
第一步:假设风轮轮毂高度处来流风速为11m/s,风切变系数为0.3,利用公式(3)和公式(4)计算风轮空间平均风速及其对应的叶尖速比,并采用查表法确定风轮空间平均风速对应的风能利用系数,轮毂高度风速对应的理论风能利用系数以及风轮转速。通过公式(17)和公式(19)计算分析风轮等效风速随方位角的变化规律,如附图1所示。
第二步:假设风轮轮毂高度处来流风速为11m/s,采用同上的方法获取等效风速计算所需参数,并利用公式(18)和公式(20)计算分析不同风切变下风轮平均等效风速变化规律,如附图2所示。
第三步:假设风切变系数为0.3,风轮轮毂高度处来流风速范围为4~25m/s,采用式(18)和式(20)计算不同来流风速下的风轮等效风速,并计算分析等效风速与来流风速偏差ΔV的变化规律,如附图3所示。
ΔV=Veq-VH (21)
通过两种等效风速计算模型的对比分析,可以看出本发明所提模型能够有效地反映风轮等效风速在风电机组不同控制阶段的变化规律,更符合风电机组实际运行状况。
此实施例仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (1)
1.一种基于等效功率的风轮等效风速计算方法,其特征在于,包括:
步骤1:考虑风切变效应对风轮扫略平面内风速的影响来确定叶片某方位角处的风速计算式;
步骤2:采用三阶泰勒级数展开简化风速计算式,并利用微元面积法,求取风轮空间平均风速其对应的叶尖速比;
步骤3:定义风速产生的气动转矩与风轮等效风速、叶片气动转矩影响系数的关系;
步骤4:求取风轮机械转矩在风轮空间平均风速所对应的叶尖速比处的偏微分线性化,综合得到风轮转矩公式;
步骤5:通过风轮瞬时机械功率和基于等效风速的风轮机械功率求得基于等效功率的风轮瞬时等效风速;
步骤6:通过风轮平均机械功率和基于等效风速的风轮机械功率求得基于等效功率的风轮平均等效风速;
所述基于等效功率的风轮瞬时等效风速为Veq(θ);
式中,R为风轮半径,H为轮毂高度,V0为风轮空间平均风速,Ω为风轮转速,θ为叶片方位角,CP(λ,β)为轮毂高度风速对应风能利用系数,CP(λ0,β0)为风轮空间平均风速对应风能利用系数,λ0为风轮空间平均风速V0对应叶尖速比,α为风切变系数,VH为来流风速;
所述基于等效功率的风轮平均等效风速为Veq;
式中,R为风轮半径,H为轮毂高度,Ω为风轮转速,θ为叶片方位角,CP(λ,β)为轮毂高度风速对应风能利用系数,CP(λ0,β0)为风轮空间平均风速对应风能利用系数,λ0为风轮空间平均风速V0对应叶尖速比,α为风切变系数,VH为来流风速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910136172.4A CN109753759B (zh) | 2019-02-25 | 2019-02-25 | 一种基于等效功率的风轮等效风速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910136172.4A CN109753759B (zh) | 2019-02-25 | 2019-02-25 | 一种基于等效功率的风轮等效风速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109753759A CN109753759A (zh) | 2019-05-14 |
CN109753759B true CN109753759B (zh) | 2020-10-27 |
Family
ID=66407495
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910136172.4A Active CN109753759B (zh) | 2019-02-25 | 2019-02-25 | 一种基于等效功率的风轮等效风速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109753759B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111666716B (zh) * | 2020-06-05 | 2022-05-20 | 中南大学 | 大型风电机组叶轮面等效风速预测方法 |
CN112949075B (zh) * | 2021-03-15 | 2024-03-15 | 华北电力大学 | 一种计及风切变和塔影效应的风轮等效风速计算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108700038A (zh) * | 2016-02-24 | 2018-10-23 | 乌本产权有限公司 | 用于确定等效风速的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592026A (zh) * | 2012-01-12 | 2012-07-18 | 甘肃省电力公司风电技术中心 | 一种双馈风机等效模拟的仿真建模方法 |
DE102012221823A1 (de) * | 2012-11-29 | 2014-06-05 | Zf Friedrichshafen Ag | Getriebeeinheit mit Plusgetriebesatz |
CN103617308B (zh) * | 2013-10-30 | 2016-06-08 | 河海大学 | 一种风电场频域等效模型的构建方法 |
CN105587474A (zh) * | 2015-12-10 | 2016-05-18 | 中国电力科学研究院 | 一种风电机组风轮等效风速的测试方法和装置 |
-
2019
- 2019-02-25 CN CN201910136172.4A patent/CN109753759B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108700038A (zh) * | 2016-02-24 | 2018-10-23 | 乌本产权有限公司 | 用于确定等效风速的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109753759A (zh) | 2019-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103850876B (zh) | 一种适用于无载荷测量的风电机组独立变桨控制方法 | |
Duque et al. | Navier-Stokes simulations of the NREL combined experiment phase II rotor | |
US8257040B2 (en) | System and method for wind condition estimation | |
CN108488035B (zh) | 永磁直驱风力发电机组失速和变桨混合控制方法 | |
Castaignet et al. | Frequency-weighted model predictive control of trailing edge flaps on a wind turbine blade | |
KR101734757B1 (ko) | 윈드쉬어에 의한 풍력터빈 블레이드의 추력 모델링 방법 및 시스템 | |
Ye et al. | Unsteady aerodynamic characteristics of a horizontal wind turbine under yaw and dynamic yawing | |
CN102182633B (zh) | 一种大型风电机组独立变桨控制方法及装置 | |
CN103061980A (zh) | 基于激光测风雷达的风力发电机组的前馈控制***及其控制方法 | |
CN202117846U (zh) | 一种大型风电机组独立变桨控制装置 | |
CN109753759B (zh) | 一种基于等效功率的风轮等效风速计算方法 | |
Guerri et al. | Simulations of the fluid flow around a rotating vertical axis wind turbine | |
CN112283026A (zh) | 一种基于空气密度跟踪最优模态增益的动态转矩控制方法 | |
EP3514373B1 (en) | Twist correction factor for aerodynamic performance map used in wind turbine control | |
CN105138845B (zh) | 获得风力发电机风速值的方法 | |
CN112949075B (zh) | 一种计及风切变和塔影效应的风轮等效风速计算方法 | |
CN115807734A (zh) | 一种基于尾流跟踪的海上风电场场级协同控制策略 | |
Hamlaoui et al. | Improved bem method for hawt performance predictions | |
US20210279389A1 (en) | Method of determining the wind speed in the rotor plane of a wind turbine | |
Liu et al. | Load reduction for two-bladed horizontal-axis tidal current turbines based on individual pitch control | |
CN109611268B (zh) | 一种双叶轮水平轴风力机设计优化方法 | |
Li et al. | Effects of offset blade on aerodynamic characteristics of small-scale vertical axis wind turbine | |
CN111120204A (zh) | 风力发电机组独立变桨四象限运行控制方法 | |
CN101252334A (zh) | 变速恒频风力发电机动态最优能量捕获方法 | |
Al-Abadi et al. | Development of an experimental setup for double rotor HAWT investigation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |