CN105023099A - 一种考虑湍流强度的风力发电机出力评估方法 - Google Patents
一种考虑湍流强度的风力发电机出力评估方法 Download PDFInfo
- Publication number
- CN105023099A CN105023099A CN201510416670.6A CN201510416670A CN105023099A CN 105023099 A CN105023099 A CN 105023099A CN 201510416670 A CN201510416670 A CN 201510416670A CN 105023099 A CN105023099 A CN 105023099A
- Authority
- CN
- China
- Prior art keywords
- wind
- wheel hub
- reference value
- blower fan
- series
- 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
Abstract
本发明公开了一种考虑湍流强度的风力发电机出力评估方法,首先获取所要评估的风力发电机的机型参数和气动参数,之后根据风力发电机机型参数和气动参数确定该风机关于轮毂处平均风速和湍流参考值变化的风功率曲面;接着将欲进行风机出力评估的风速序列按照分段周期T进行分段,确定每段风速序列的湍流参考值和轮毂处平均风速;最后通过风功率曲面确定该风机每个分段周期T的出力,合计全部时间段的风机出力得出风机总出力。本发明适用于按照机型分类对风力发电机进行出力评估,进而进行风电场的规划和风电场的年发电量预测等。该方法避免了只考虑平均风速影响带来的评估误差,较大程度的提高了风力发电机出力评估的准确度。
Description
技术领域
本发明属于风力机出力评估领域,尤其涉及一种考虑湍流强度的风力发电机出力评估方法。
背景技术
随着全球工业水平的飞速发展,石油等非可再生能源的储量已急剧减少,而人类社会对能源的需求却与日俱增。当人类尝试寻找风能等新能源代替传统能源时,日益突出的环境问题再一次加速了能源结构变革的速度。风能等新型能源共同为解决能源危机和环境问题提供了新的思路。风电的主要优势在于建设风电场的周期比较短,风电机组占地面积比较小,基本不影响风场及周边生态环境,因此风能的开发和利用在70年代以后便有了长足的进步。对于风能的利用,如在进行风力发电场的年发电量计算等场景中,都涉及到了风力发电出力的计算。
根据研究,风力发电机的出力不但与轮毂处的平均风速有关,还与湍流强度密不可分。目前传统的风力发电机的出力评估方法,一般只在某一固定湍流强度下根据风功率曲线进行风机出力的计算。由于这种方法只考虑平均风速而没有考虑湍流强度的变化,所以对风机出力的计算带来一定的误差。而风机出力的计算方法已广泛应用于风功率预测、风电场的年发电量预测等场景,较准确的风电场功率计算已成为风力发电稳定运行以及包含风力发电的电网安全调度的前提。
基于上述情况,目前迫切需要一种新的风力发电机功率计算方法,同时考虑湍流强度与平均风速对风机出力的影响。但是现有技术中尚无相关描述。
发明内容
本发明的目的在于提供一种考虑湍流强度的风力发电机出力评估方法。
实现本发明目的的技术解决方案为:一种考虑湍流强度的风力发电机出力评估方法,包括以下步骤:
步骤1、获取所要评估的风力发电机的机型参数和气动参数,包括:额定转速、额定转矩、风轮半径风机轮毂高度、塔高、机舱质量、轮毂质量、切入风速、空气密度、空气运动粘滞系数;
步骤2、根据步骤1的风力发电机机型参数和气动参数确定该风机关于轮毂处平均风速V和湍流参考值Iref变化的风功率曲面;包括以下步骤:
步骤2-1、设定湍流参考值变化范围和步长;湍流参考值变化范围为0.05到0.19,步长为0.01。
步骤2-2、在每一个湍流参考值的条件下绘制一条评估风机关于轮毂处平均风速V变化的风功率曲线,得到风功率曲线簇;
步骤2-3、将风功率曲线簇绘制成评估风机关于轮毂处平均风速V和湍流参考值Iref变化的风功率曲面。
步骤3、将欲进行风机出力评估的风速序列按照分段周期T进行分段,确定每段风速序列的湍流参考值Iref和轮毂处平均风速V;具体为:
步骤3-1、将欲进行风机出力评估的风速序列,按照分段周期T进行分段,分成每周期T为一段的风速序列;
步骤3-2、对于每个T时间的风速序列,确定该序列的轮毂处平均风速;所用公式为:
式中,V为该段风速序列轮毂处平均风速,N为该段风速序列中风速数据的个数,Vj为该段风速序列中轮毂处第j个的风速;
步骤3-3、确定每个T时间风速序列的湍流强度IT,根据湍流强度IT求出序列的湍流参考值Iref,所用公式为:
式中,Iref为该段风速序列当风速为15m/s时的湍流强度的预期值,即本文中的湍流参考值,b=5.6m/s,IT为该段风速序列的湍流强度,其计算公式为σ为该段风速序列的风速标准差,其计算公式为所述时间T为10分钟。
步骤4、根据步骤3确定的每段风速序列的湍流参考值Iref和轮毂处平均风速V,通过风功率曲面确定该风机每个分段周期T的出力Pi,最后合计全部时间段的风机出力得出风机总出力P。合计全部时间段的风机出力得出风机总出力所用公式为:
式中,Pi为第i个风速序列的指定风机的出力,Ns为风速序列个数,P为全部时间段内指定风机的全部出力。
本发明与现有技术相比,其显著优点为:本发明提出一种考虑湍流强度的风力发电机出力评估方法,通过生成多条不同湍流参考值的风功率曲线形成三维的风功率曲面,考虑湍流强度与轮毂处平均风速的共同影响进行风机出力评估。本发明能有效提高计算风力机输出功率的精确度,从而保证风电场功率输出模型的准确性,为风电场的规划、调度提供相关可靠依据。
附图说明
图1为本发明的风机出力评估方法流程图。
图2为本发明的风力发电机功率曲线族。
图3为本发明的风力发电机功率曲面。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
一种考虑湍流强度的风力发电机功率曲面计算方法,如图1所示,包括以下步骤:
步骤1、首先获取所要评估的风力发电机的机型参数和气动参数,包括:额定转速、额定转矩、风轮半径风机轮毂高度、塔高、机舱质量、轮毂质量、切入风速、空气密度、空气运动粘滞系数。如表1所示,为1.5WM级风力发电机的机型参数和气动参数。
表1 1.5MW的风力发电机机型参数和气动参数
参数 | 数值 |
切入风速 | 4(m/s) |
额定转速 | 1800(rpm) |
额定转矩 | 8376.58(Nm) |
风轮半径 | 35(m) |
空气密度 | 1.255(kg/m^3) |
风机轮毂高度 | 84.2876(m) |
空气运动粘滞系数 | 1.4639E-5(m^2/sec) |
塔高 | 82.39(m) |
机舱质量 | 51170(kg) |
轮毂质量 | 15148(kg) |
步骤2、针对指定型号的风力机,利用风电场的实时统计数据,或者专业的风力机仿真软件(如Bladed、FAST等)参照IEC64100-12-1标准中风功率曲线的绘制方法,绘制指定风机在某一湍流参考值Iref下,随轮毂处平均风速V变化的风功率曲线。
湍流参考值由0.05遍历到0.19,步长为0.01,绘制出每个湍流参考值条件下指定风机关于风速变化的风功率曲线族,如图2所示。最终将这簇风功率曲线绘制成指定风机关于轮毂处平均风速V和湍流参考值Iref变化的三维风功率曲面。
步骤3、将欲进行风机出力计算的风速序列,分成每10分钟一段的风速序列。
步骤4、对每个10分钟序列,计算每个10分钟序列的轮毂处平均风速V。
式中,V为风速序列轮毂处平均风速,N为该段风速序列中风速数据的个数,Vj为该段风速序列中轮毂处第j个的风速。
步骤5、计算每个10分钟序列的湍流强度IT,其计算公式如下:
式中,IT为该段风速序列的湍流强度,σ为该段风速序列的风速标准差,其计算公式为
根据湍流强度IT求出序列的湍流参考值Iref:
式中,Iref为该段风速序列当风速为15m/s时的湍流强度的预期值,即本文中的湍流参考值,b=5.6m/s。
步骤6、对每个10分钟序列,根据湍流参考值Iref和轮毂处平均风速V,在步骤1绘制的三维风功率曲面中得出该10分钟风机的出力Pi(Iref,V)。
步骤7、合计全部时间段的风机处理,得出风机总出力:
式中,Pi为第i个风速序列的指定风机的出力,Ns为风速序列个数,P为全部时间段内指定风机的全部出力。
下面结合实施例对本发明做进一步详细的描述:
实施例1
以1.5WM的风力发电机为例,具体参数如表1所示。首先使用TurbSim(美国国家能源部可再生能源实验室(NREL)提供的开源的湍流风模拟软件)按照自然条件,生成不同湍流强度,时长为20小时,频率为20Hz,轮毂处平均风速的变化范围从4m/s到13m/s的风速数据,共计7200000个数据点。
利用美国国家能源部可再生能源实验室(NREL)提供的开源的专业风力机仿真软件FAST(Fatigue,Aerodynamics,Structures,and Turbulence)模拟1.5WM风力发电机的出力,在该风况下1.5WM的风力发电机20小时的功率为4054.8kW·h,以此作为参考。
利用本专利的评估方法,绘制出1.5WM风力发电机关于轮毂处平均风速V和湍流参考值Iref变化的风功率曲面(如图3所示),经过计算,使用本方法计算1.5WM风力发电机20小时的功率为4128.2kW·h,而使用传统的功率曲线法计算结果为4325.5kW·h。因此,使用本发明的方法优于传统的功率曲线计算法。
Claims (6)
1.一种考虑湍流强度的风力发电机出力评估方法,其特征在于,包括以下步骤:
步骤1、获取所要评估的风力发电机的机型参数和气动参数,包括:额定转速、额定转矩、风轮半径风机轮毂高度、塔高、机舱质量、轮毂质量、切入风速、空气密度、空气运动粘滞系数;
步骤2、根据步骤1的风力发电机机型参数和气动参数确定该风机关于轮毂处平均风速V和湍流参考值Iref变化的风功率曲面;
步骤3、将欲进行风机出力评估的风速序列按照分段周期T进行分段,确定每段风速序列的湍流参考值Iref和轮毂处平均风速V;
步骤4、根据步骤3确定的每段风速序列的湍流参考值Iref和轮毂处平均风速V,通过风功率曲面确定该风机每个分段周期T的出力Pi,最后合计全部时间段的风机出力得出风机总出力P。
2.根据权利要求1所述的考虑湍流强度的风力发电机出力评估方法,其特征在于,步骤2中确定该风机关于轮毂处平均风速V和湍流参考值Iref变化的风功率曲面时,包括以下步骤:
步骤2-1、设定湍流参考值变化范围和步长;
步骤2-2、在每一个湍流参考值的条件下绘制一条评估风机关于轮毂处平均风速V变化的风功率曲线,得到风功率曲线簇;
步骤2-3、将风功率曲线簇绘制成评估风机关于轮毂处平均风速V和湍流参考值Iref变化的风功率曲面。
3.根据权利要求2所述的考虑湍流强度的风力发电机出力评估方法,其特征在于,步骤2-1中湍流参考值变化范围为0.05到0.19,步长为0.01。
4.根据权利要求1所述的考虑湍流强度的风力发电机出力评估方法,其特征在于,步骤3将欲进行风机出力评估的风速序列按照分段周期T进行分段,确定每段风速序列的湍流参考值Iref和轮毂处平均风速V具体为:
步骤3-1、将欲进行风机出力评估的风速序列,按照分段周期T进行分段,分成每周期T为一段的风速序列;
步骤3-2、对于每个T时间的风速序列,确定该序列的轮毂处平均风速;所用公式为:
式中,V为该段风速序列轮毂处平均风速,N为该段风速序列中风速数据的个数,Vj为该段风速序列中轮毂处第j个的风速;
步骤3-3、确定每个T时间风速序列的湍流强度IT,根据湍流强度IT求出序列的湍流参考值Iref,所用公式为:
式中,Iref为该段风速序列当风速为15m/s时的湍流强度的预期值,即本文中的湍流参考值,b=5.6m/s,IT为该段风速序列的湍流强度,其计算公式为σ为该段风速序列的风速标准差,其计算公式为
5.根据权利要求4所述的考虑湍流强度的风力发电机出力评估方法,其特征在于,时间T为10分钟。
6.根据权利要求1述的考虑湍流强度的风力发电机出力评估方法,其特征在于,步骤4合计全部时间段的风机出力得出风机总出力所用公式为:
式中,Pi为第i个风速序列的指定风机的出力,Ns为风速序列个数,P为全部时间段内指定风机的全部出力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510416670.6A CN105023099B (zh) | 2015-07-15 | 2015-07-15 | 一种考虑湍流强度的风力发电机出力评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510416670.6A CN105023099B (zh) | 2015-07-15 | 2015-07-15 | 一种考虑湍流强度的风力发电机出力评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105023099A true CN105023099A (zh) | 2015-11-04 |
CN105023099B CN105023099B (zh) | 2018-07-06 |
Family
ID=54413052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510416670.6A Active CN105023099B (zh) | 2015-07-15 | 2015-07-15 | 一种考虑湍流强度的风力发电机出力评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105023099B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105699894A (zh) * | 2016-01-27 | 2016-06-22 | 云南电网有限责任公司电力科学研究院 | 一种风力发电机在线监测方法及系统 |
CN106932607A (zh) * | 2017-03-09 | 2017-07-07 | 北京玖天气象科技有限公司 | 一种简便判断测风数据受测风塔体遮挡影响的方法和装置 |
CN107654342A (zh) * | 2017-09-21 | 2018-02-02 | 湘潭大学 | 一种考虑湍流的风电机组功率异常的检测方法 |
CN112628069A (zh) * | 2020-12-18 | 2021-04-09 | 南京理工大学 | 考虑湍流风速对转速恢复影响的风机惯性响应控制方法及系统 |
CN115564161A (zh) * | 2022-03-23 | 2023-01-03 | 师亚琦 | 一种智慧电网发电设施工程项目建设评估分析方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10988923B2 (en) * | 2012-11-30 | 2021-04-27 | Rensselaer Polytechnic Institute | Methods and systems of modifying air flow at building structures |
CN103745024B (zh) * | 2013-11-26 | 2018-12-04 | 沈阳工业大学 | 基于三维尾流模型修正风电机组尾部风速功率特性评估法 |
-
2015
- 2015-07-15 CN CN201510416670.6A patent/CN105023099B/zh active Active
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105699894A (zh) * | 2016-01-27 | 2016-06-22 | 云南电网有限责任公司电力科学研究院 | 一种风力发电机在线监测方法及系统 |
CN106932607A (zh) * | 2017-03-09 | 2017-07-07 | 北京玖天气象科技有限公司 | 一种简便判断测风数据受测风塔体遮挡影响的方法和装置 |
CN106932607B (zh) * | 2017-03-09 | 2020-04-24 | 北京玖天气象科技有限公司 | 一种简便判断测风数据受测风塔体遮挡影响的方法和装置 |
CN107654342A (zh) * | 2017-09-21 | 2018-02-02 | 湘潭大学 | 一种考虑湍流的风电机组功率异常的检测方法 |
CN112628069A (zh) * | 2020-12-18 | 2021-04-09 | 南京理工大学 | 考虑湍流风速对转速恢复影响的风机惯性响应控制方法及系统 |
CN112628069B (zh) * | 2020-12-18 | 2022-03-18 | 南京理工大学 | 考虑湍流风速对转速恢复影响的风机惯性响应控制方法及系统 |
CN115564161A (zh) * | 2022-03-23 | 2023-01-03 | 师亚琦 | 一种智慧电网发电设施工程项目建设评估分析方法 |
CN115564161B (zh) * | 2022-03-23 | 2023-11-17 | 甘肃电通电力工程设计咨询有限公司 | 一种智慧电网发电设施工程项目建设评估分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105023099B (zh) | 2018-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Simulation of turbulent flow inside and above wind farms: model validation and layout effects | |
Adaramola et al. | Experimental investigation of wake effects on wind turbine performance | |
CN103244348B (zh) | 变速变桨风力发电机组功率曲线优化方法 | |
CN105023099B (zh) | 一种考虑湍流强度的风力发电机出力评估方法 | |
CN103886179B (zh) | 一种基于尾流效应分群的风电场风机聚合方法 | |
CN101592538B (zh) | 一种基于实测数据风电场稳态输出功率的计算方法 | |
CN101949363A (zh) | 计及风电场输入风速和风向随机波动的风电机组分组方法 | |
CN103745024B (zh) | 基于三维尾流模型修正风电机组尾部风速功率特性评估法 | |
CN102663251A (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
CN101852172A (zh) | 一种考虑风电场内尾流影响的风电机组输入风速计算方法 | |
CN103617308B (zh) | 一种风电场频域等效模型的构建方法 | |
CN106815773A (zh) | 一种风电场功率特性评估方法 | |
CN107947228B (zh) | 基于Markov理论的含风电电力系统随机稳定性分析方法 | |
CN103699804B (zh) | 一种实时监测风电机组运行时的功率出力偏差的方法 | |
Salami et al. | Evaluation of Wind Potential for an Optimum Choice of Wind Turbine Generator on the Sites of Lomé, Accra, and Cotonou Located in the Gulf of Guinea. | |
CN106919730B (zh) | 一种采用风速衰减因子的风电场尾流确定方法 | |
CN103696913A (zh) | 一种实时监测风电机组运行时的风能利用偏差的方法 | |
CN106779202B (zh) | 一种考虑空气湿度的风电功率预测方法 | |
Deskos et al. | Development and validation of the higher-order finite-difference wind farm simulator, WInc3D | |
CN109038651B (zh) | 一种基于机组相关性的风电场功率波动频谱预测模型 | |
Hong et al. | The design and testing of a small-scale wind turbine fitted to the ventilation fan for a livestock building | |
CN103944159B (zh) | 含风电的电力系统规划中检修备用容量确定方法 | |
Tjahjana et al. | Wind potential assessment to estimate performance of selected wind turbines in Pandansimo Beach-Yogyakarta | |
CN113468767B (zh) | 一种海上风电机组发电量评估方法与系统 | |
Charhouni et al. | Analysis of wake impact on wind farm performance using two analytical models |
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 |