CN109488525B - 基于提高转速下限的转速跟踪目标优化方法 - Google Patents

基于提高转速下限的转速跟踪目标优化方法 Download PDF

Info

Publication number
CN109488525B
CN109488525B CN201811335726.5A CN201811335726A CN109488525B CN 109488525 B CN109488525 B CN 109488525B CN 201811335726 A CN201811335726 A CN 201811335726A CN 109488525 B CN109488525 B CN 109488525B
Authority
CN
China
Prior art keywords
rotating speed
speed
omega
tracking target
wind
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
Application number
CN201811335726.5A
Other languages
English (en)
Other versions
CN109488525A (zh
Inventor
陈载宇
杨晶
瞿佳刘
殷明慧
汪成根
刘建坤
李群
周前
卜京
邹云
张宁宇
卫鹏
吴盛军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Nanjing University of Science and Technology
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical Nanjing University of Science and Technology
Priority to CN201811335726.5A priority Critical patent/CN109488525B/zh
Publication of CN109488525A publication Critical patent/CN109488525A/zh
Application granted granted Critical
Publication of CN109488525B publication Critical patent/CN109488525B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • F05B2270/304Spool rotational speed
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • F05B2270/32Wind speeds
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind 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)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种基于提高转速下限的转速跟踪目标优化方法,针对大转动惯量风机的慢动态特性使得风机在转速过低时无法及时加速而导致的风能损失问题,该方法在传统叶尖速比法的基础上,引入一个低通滤波器,对理论最优转速进行滤波,然后对风机的运行转速下限进行限制,同时将参考转速整体上移,从而得到新的转速跟踪目标。通过设置转速下限,使得当风速突增时,风机能够更快地追踪上当前风速。本发明虽然放弃了低风速下的转速跟踪,但是可以有效提高风机在高风速区域的动态性能,使得风机能够捕获更多的风能。

Description

基于提高转速下限的转速跟踪目标优化方法
技术领域
本发明属于风机控制领域,具体涉及一种基于提高转速下限的转速跟踪目标优化方法。
背景技术
变速恒频风力发电机组一般采用最大功率点跟踪,通过控制风机的运行转速,使得风电机组始终处于最大功率点,从而实现效率最大化,其实现方法主要为叶尖速比、最优转矩法、爬山法。
叶尖速比法实现起来较为简单,其通过估计算法计算得到风速信息,在此基础上计算最优转速,即当风机以此转速运行时,可实现最大的风能捕获。叶尖速比法的控制器通过判断实际转速与最优转速,调节风机的电磁转矩,实现对最优转速的跟踪。
然而,最优转速根据风速计算而来,具备较强波动性的特点,对于大转动惯量风机而言,因其固有的慢动态特性,使得风机转速很难跟踪上最优转速。另一方面,对于低风速风机,用以驱动风机加速的气动转矩较小,故发电机电磁转矩的下限对风机最大功率点跟踪控制的实现会产生更大的影响。
发明内容
本发明的目的在于提供一种基于提高转速下限的转速跟踪目标优化方法,通过对最优转速进行滤波,限制转速跟踪目标的下限,同时将转速跟踪目标整体上移,从而实现对转速跟踪目标的优化。
实现本发明目的的技术解决方案为:一种基于提高转速下限的转速跟踪目标优化方法,包括以下步骤:
步骤1、获取风机的相关参数,包括风速、风轮转速、风轮半径;
步骤2、根据所得风速,计算理论最优转速;
步骤3、选取滤波参数,对理论最优转速进行滤波,得到滤波后转速;
步骤4、计算用于最大功率点跟踪控制的起始转速;
步骤5、当滤波后转速ωf大于起始转速ωbgn时,参考转速取为ωf;否则,转步骤6;
步骤6、当滤波后转速ωf小于或等于起始转速ωbgn时,参考转速取为ωbgn
步骤7、计算参考转速的上移量,将参考转速整体上移,得到转速跟踪目标。
本发明与现有技术相比,其显著优点为:(1)本发明提出了一种基于提高转速下限的转速跟踪目标优化方法,解决了风机在转速过低时无法及时加速而导致的风能损失问题;(2)本发明在保证风机跟踪效果的同时,有效提高了风能捕获效率。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为本发明所提基于提高转速下限的转速跟踪目标优化方法流程图。
图2本发明有效性的仿真结果,其中图2(a)是传统叶尖速比法与优化方法的转速跟踪目标仿真结果图,图2(b)是最优转速、传统方法与优化方法的实际风轮转速仿真结果图。
具体实施方式
对于低风速大转动惯量风机,传统的以最优转速为跟踪目标的叶尖速比难以跟踪最大功率点。本发明中从两方面对叶尖速比法的转速跟踪目标进行优化:一方面,平滑理论最优转速,使风机更易跟踪;另一方面,借鉴收缩跟踪区间的思想,通过限制转速跟踪目标的下限,同时整体上移转速跟踪目标,以提高高风速区域的风能捕获。
如图1所示,一种基于提高转速下限的转速跟踪目标优化方法,包括以下步骤:
步骤1、获取风机的相关参数,包括风速v、风轮转速ωr、风轮半径R;
步骤2、根据所得风速,计算理论最优转速ωopt
步骤3、选取滤波参数αf,对理论最优转速进行滤波,得到滤波后转速ωf
步骤4、计算用于最大功率点跟踪控制的起始转速ωbgn
步骤5、当滤波后转速ωf大于起始转速ωbgn时,参考转速取为ωf;否则,转步骤6;
步骤6、当滤波后转速ωf小于或等于起始转速ωbgn时,参考转速取为ωbgn
步骤7、计算参考转速的上移量ωu,将参考转速整体上移,得到转速跟踪目标ωref
进一步的,步骤2中的理论最优转速ωopt具体计算方式如下:
Figure BDA0001861192370000021
其中:λopt是最优叶尖速比,对固定的风机来说是定值;v为风机轮毂高度处的实时风速;R为风轮半径。
进一步的,步骤3中对最优转速进行滤波的具体形式:
理论最优转速是根据风速计算而来,所以具有风速不确定性且波动大的特点,实际上,风轮很难跟踪上这种剧烈波动的跟踪目标,因此采用一阶低通滤波器对理论最优转速进行平滑滤波,得到滤波后转速ωf
ωf(k)=(1-αfopt(k)+αfωf(k-1) (2)
其中:αf为滤波参数,取值在0.98以上,本发明中取为0.99;k为离散时间周期,ωopt(k)表示本周期的理论最优转速值,ωf(k-1)为上一周期滤波后转速值。
进一步的,步骤4中计算用于最大功率点跟踪控制的起始转速ωbgn的具体形式:
Figure BDA0001861192370000031
其中:vbgn为开始进行最大功率点跟踪的起始风速,可由下式获得:
Figure BDA0001861192370000032
其中:
Figure BDA0001861192370000033
为一段时间内的平均风速,σ为该短时间内风速的标准差。
进一步的,步骤5中获取参考转速的具体形式:
当滤波后转速ωf大于起始转速ωbgn,即ωf>ωbgn时,参考转速ωlim取为:
ωlim=ωf (5)
进一步的,步骤6中获取参考转速的具体形式:
当滤波后转速ωf小于或等于起始转速ωbgn,即ωf≤ωbgn时,参考转速ωlim取为:
ωlim=ωbgn (6)
进一步的,步骤7中转速跟踪目标的具体形式:
理论上,参考转速最优上移量同时与风速变化、风力机动态、CP-λ曲线形态等多方面因素相关,难以直接计算最优上移量。在本发明,取一段时间最优转速标准差的30%作为参考转速的上移量ωu
Figure BDA0001861192370000034
其中:λopt是最优叶尖速比;R为风轮半径;σ为一段时间内最优转速的标准差;本发明取10分钟内最优转速标准差的30%作为参考转速的上移量以及最优转速的标准差。
转速跟踪目标取为步骤5、6中的参考转速加上上移量:
ωref=ωlimu (8)
下面结合实施例对本发明做进一步详细的描述:
实施例
利用美国国家能源部可再生能源实验室(National Renewable EnergyLaboratory,NREL)提供的开源的专业风力机仿真软件FAST(Fatigue,Aerodynamics,Structures,and Turbulence)来模拟控制效果。风力机模型采用NREL开发的600kW CART3试验机型,具体参数如表1所示。
表1 NREL 600kW CART3风力机主要参数
Figure BDA0001861192370000041
本发明中将基于提高转速下限的转速跟踪目标优化方法与传统叶尖速比法进行对比,传统叶尖速比法所用转速跟踪目标即为理论最优转速,本发明所提优化转速跟踪目标可由图1所示流程获取。
选用平均风速600m/s、A级湍流强度、积风尺度为150m的湍流风速,分别用两种方法进行仿真,仿真结果如图2所示:图2(a)为传统叶尖速比法与本发明所提优化方法的转速跟踪目标对比图,可以看出经过滤波之后,转速跟踪目标变得更为平滑,且转速跟踪目标被限制在1.5rad/s以上;图2(b)为两种方法的实时风轮转速对比,因转速跟踪目标不同,所以优化方法的实时转速更为平滑,波动较小;得益于提高转速跟踪目标的下限,在低转速风速突增的情况下,优化的转速跟踪目标使得风机更快地追踪上当前风速,大大提高了风机在渐强阵风下的动态性能。
该仿真算例下,传统叶尖速比法的风能捕获效率为42.85%,本发明中所提优化方法的风能捕获效率为43.63%。可见,本发明所提的一种基于提高转速下限的转速跟踪目标优化方法有效提高了高风速区域的风能捕获性能。

Claims (7)

1.一种基于提高转速下限的转速跟踪目标优化方法,其特征在于,包括以下步骤:
步骤1、获取风机的相关参数,包括风速、风轮转速、风轮半径;
步骤2、根据所得风速,计算理论最优转速;
步骤3、选取滤波参数,对理论最优转速进行滤波,得到滤波后转速;
步骤4、计算用于最大功率点跟踪控制的起始转速;
步骤5、当滤波后转速ωf大于起始转速ωbgn时,参考转速取为ωf;否则,转步骤6;
步骤6、当滤波后转速ωf小于或等于起始转速ωbgn时,参考转速取为ωbgn
步骤7、计算参考转速的上移量,将参考转速整体上移,得到转速跟踪目标。
2.根据权利要求1所述的基于提高转速下限的转速跟踪目标优化方法,其特征在于,步骤2中的理论最优转速ωopt计算方法如下:
Figure FDA0002427321250000011
其中,λopt为最优叶尖速比,v为风机轮毂高度处的实时风速,R为风轮半径。
3.根据权利要求1所述的基于提高转速下限的转速跟踪目标优化方法,其特征在于,步骤3中对最优转速进行滤波的具体方法为:
采用一阶低通滤波器对理论最优转速进行平滑滤波,得到滤波后转速ωf
ωf(k)=(1-αfopt(k)+αfωf(k-1) (2)
其中,αf为滤波参数,k为离散时间周期,ωopt(k)表示本周期的理论最优转速值,ωf(k-1)为上一周期滤波后转速值。
4.根据权利要求3所述的基于提高转速下限的转速跟踪目标优化方法,其特征在于,αf取值在0.98以上。
5.根据权利要求1所述的基于提高转速下限的转速跟踪目标优化方法,其特征在于,步骤4中计算用于最大功率点跟踪控制的起始转速ωbgn的具体方法为:
Figure FDA0002427321250000012
其中,vbgn为开始进行最大功率点跟踪的起始风速,由下式获得:
Figure FDA0002427321250000021
其中,
Figure FDA0002427321250000022
为一段时间内的平均风速,σ为该段时间内风速的标准差。
6.根据权利要求1所述的基于提高转速下限的转速跟踪目标优化方法,其特征在于,步骤7中转速跟踪目标的计算方法为:
取一段时间内最优转速标准差的30%作为参考转速的上移量ωu
Figure FDA0002427321250000023
其中,λopt是最优叶尖速比,R为风轮半径,σ为一段时间内最优转速的标准差;
转速跟踪目标取为步骤5、6中的参考转速加上上移量,即:
ωref=ωlimu (8)
其中,ωlim为参考转速。
7.根据权利要求6所述的基于提高转速下限的转速跟踪目标优化方法,其特征在于,步骤7中取10分钟最优转速标准差的30%作为参考转速的上移量ωu
CN201811335726.5A 2018-11-11 2018-11-11 基于提高转速下限的转速跟踪目标优化方法 Active CN109488525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811335726.5A CN109488525B (zh) 2018-11-11 2018-11-11 基于提高转速下限的转速跟踪目标优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811335726.5A CN109488525B (zh) 2018-11-11 2018-11-11 基于提高转速下限的转速跟踪目标优化方法

Publications (2)

Publication Number Publication Date
CN109488525A CN109488525A (zh) 2019-03-19
CN109488525B true CN109488525B (zh) 2020-07-03

Family

ID=65694240

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811335726.5A Active CN109488525B (zh) 2018-11-11 2018-11-11 基于提高转速下限的转速跟踪目标优化方法

Country Status (1)

Country Link
CN (1) CN109488525B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112211781B (zh) * 2020-08-31 2023-03-31 南京理工大学 综合收缩跟踪区间和减小转矩增益的风机转矩曲线法
CN112377366B (zh) * 2020-11-23 2022-05-10 中国船舶重工集团海装风电股份有限公司 一种切入风速的优化方法、装置、设备及可读存储介质
CN113464378A (zh) * 2021-07-13 2021-10-01 南京理工大学 一种基于深度强化学习的提升风能捕获的转速跟踪目标优化方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102434391A (zh) * 2011-12-27 2012-05-02 南京理工大学 一种基于起始转速调整的改进最大功率点跟踪控制方法
CN102635499A (zh) * 2012-04-18 2012-08-15 中船重工(重庆)海装风电设备有限公司 一种风力发电机组转速转矩控制装置及方法
CN103758697A (zh) * 2014-01-14 2014-04-30 南京理工大学 一种基于风机有效跟踪区间的改进最大功率跟踪控制方法
CN103835878A (zh) * 2013-04-07 2014-06-04 南京理工大学 基于神经网络优化起始转速的最大功率点跟踪控制方法
CN104074679A (zh) * 2014-07-02 2014-10-01 国电联合动力技术有限公司 一种变速变桨距风电机组全风速限功率优化控制方法
CN105134487A (zh) * 2015-08-24 2015-12-09 南京理工大学 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法
CN105673322A (zh) * 2016-01-28 2016-06-15 南京理工大学 实现风力机mppt控制的变参数非线性反馈控制方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102434391A (zh) * 2011-12-27 2012-05-02 南京理工大学 一种基于起始转速调整的改进最大功率点跟踪控制方法
CN102635499A (zh) * 2012-04-18 2012-08-15 中船重工(重庆)海装风电设备有限公司 一种风力发电机组转速转矩控制装置及方法
CN103835878A (zh) * 2013-04-07 2014-06-04 南京理工大学 基于神经网络优化起始转速的最大功率点跟踪控制方法
CN103758697A (zh) * 2014-01-14 2014-04-30 南京理工大学 一种基于风机有效跟踪区间的改进最大功率跟踪控制方法
CN104074679A (zh) * 2014-07-02 2014-10-01 国电联合动力技术有限公司 一种变速变桨距风电机组全风速限功率优化控制方法
CN105134487A (zh) * 2015-08-24 2015-12-09 南京理工大学 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法
CN105673322A (zh) * 2016-01-28 2016-06-15 南京理工大学 实现风力机mppt控制的变参数非线性反馈控制方法

Also Published As

Publication number Publication date
CN109488525A (zh) 2019-03-19

Similar Documents

Publication Publication Date Title
CN109488525B (zh) 基于提高转速下限的转速跟踪目标优化方法
CN105986961B (zh) 一种变速变桨风力机功率优化控制方法
CN108533451B (zh) 一种风力发电机组的变桨控制方法
CN108843489B (zh) 基于限转速平滑功率控制的风机变桨优化方法
CN105134487B (zh) 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法
CN108825434B (zh) 基于风轮动能平滑功率控制的风机变桨优化方法
CN111209638B (zh) 基于运行攻角偏差修正的低风速风力机叶片气动设计方法
CN105508134A (zh) 一种提高风力发电机组风能利用效率的方法及装置
CN107218175B (zh) 一种实现风力机最大化风能捕获效率的转速跟踪目标优化方法
CN108843494B (zh) 基于斜线平滑功率控制的风机变桨优化方法
CN111786424B (zh) 一种风电场惯量响应及一次调频潜力的量化评估方法
CN106777525B (zh) 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法
CN114294155A (zh) 风电机组有功功率控制方法及装置
CN109139363B (zh) 一种提升多机型风力机性能的最大功率点跟踪控制方法
CN111502914B (zh) 基于线性变参数系统的风机变桨控制器设计方法
CN107194122B (zh) 一种变速风机叶片的多工况气动优化的改进方法
CN112664393B (zh) 基于最大不变桨风速范围指标的风机有功功率控制方法
CN112211782B (zh) 基于起始转速自适应搜索的收缩跟踪区间控制方法
CN115182844B (zh) 一种变速风力发电机组有界ude转矩控制方法
CN112682258B (zh) 基于Backstepping的大型风力机最大功率点跟踪控制方法
US11846270B2 (en) Control method and device of a wind park
CN109707565B (zh) 一种利用大型风力机惯性的最大风能捕获控制方法
CN105986962B (zh) 一种风力发电机组的最大风能捕获方法
CN113464378A (zh) 一种基于深度强化学习的提升风能捕获的转速跟踪目标优化方法
CN109944822B (zh) 基于限制参考转速变化率的转速跟踪目标优化方法

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