CN114263564A - 考虑不变桨风速范围的风电机组有功功率控制方法及系统 - Google Patents

考虑不变桨风速范围的风电机组有功功率控制方法及系统 Download PDF

Info

Publication number
CN114263564A
CN114263564A CN202111543163.0A CN202111543163A CN114263564A CN 114263564 A CN114263564 A CN 114263564A CN 202111543163 A CN202111543163 A CN 202111543163A CN 114263564 A CN114263564 A CN 114263564A
Authority
CN
China
Prior art keywords
wind
speed
pitch angle
wind speed
range
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
Application number
CN202111543163.0A
Other languages
English (en)
Other versions
CN114263564B (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
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Nanjing University of Science and Technology
State Grid Jiangsu Electric Power Co Ltd
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, State Grid Jiangsu Electric Power Co Ltd, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical Nanjing University of Science and Technology
Priority to CN202111543163.0A priority Critical patent/CN114263564B/zh
Publication of CN114263564A publication Critical patent/CN114263564A/zh
Application granted granted Critical
Publication of CN114263564B publication Critical patent/CN114263564B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Abstract

本发明公开了一种考虑不变桨风速范围的风电机组有功功率控制方法及系统,方法包括:获取风机相关参数;计算风机的不变桨风速范围;分析桨距角对不变桨风速范围的影响规律;建立被动变速风轮桨距角与风速变化范围之间的适配关系;实现考虑不变桨风速范围的风机有功功率控制。本发明提出的考虑不变桨风速范围的风机有功功率控制方法,能够使被动变速风轮不变桨风速范围动态匹配风速实际波动范围,降低风轮达到转速边界的频次,缓解风轮超速和电磁功率跌落。

Description

考虑不变桨风速范围的风电机组有功功率控制方法及系统
技术领域
本发明属于风机控制领域,特别是一种考虑不变桨风速范围的风电机组有功功率控制方法及系统。
背景技术
随着风电大规模、高渗透率接入电网,逐步形成了新能源为主体的新型电力系统。在此背景下,迫切需要变速变桨风电机组从通常采用的最大功率点跟踪(MPPT)控制转变为跟踪电网功率指令的有功功率控制(APC),以参与电网的自动发电控制、支撑系统的有功功率平衡。
区别于MPPT控制,在被动变速APC控制下,风机气动功率与电磁功率完全解耦,其传动链将承受幅值更大、变化更快的不平衡功率。这不仅导致风机机电动态稳定更难维持,而且大幅增加变桨机构的调节负担。
为了更充分地利用风轮被动变速、减轻变桨机构的动作疲劳,现有研究提出了集成变速-变桨APC控制策略(IAPC)和分段桨距控制,利用任意桨距角风轮被动变速应对风速波动,有效减少了变桨动作。但是现有的APC方法中的变桨环节均是在风轮转速达到边界时启动,实现转速边界处的限转速控制,被动变速风轮的桨距角具有很强的随机性和不确定性。而且变桨执行机构动作缓慢,电磁功率不参与风机机电动态调节,风机更容易出现转速超调,造成风轮超速或电磁功率跌落问题,严重影响了风机有功功率控制性能。
发明内容
本发明的目的在于提供一种考虑不变桨风速范围的风电机组有功功率控制的方法,通过增加适配桨距角的设定环节,使得被动变速风轮的不变桨风速范围动态匹配风速波动范围,降低风轮达到转速边界的频次,进一步缓解风轮超速和电磁功率跌落。
实现本发明目的的技术解决方案为:一种考虑不变桨风速范围的风电机组有功功率控制方法,所述方法包括以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure BDA0003414937000000021
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000022
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
进一步地,步骤2所述计算风机的不变桨风速范围
Figure BDA0003414937000000023
具体过程包括:
对于功率指令为Pcmd和变速区间设定为[ωLU]的有功功率控制,当桨距角为β时,存在一个最大的风速波动范围
Figure BDA0003414937000000024
Figure BDA0003414937000000025
使得对于该风速范围内的任一风速,风机都存在一个稳定平衡点(ω00),且满足ω0∈[ωLU]与β0=β;式中,ωL与ωU分别为变速区间的下边界与上边界,vL(β)与vU(β)分别为桨距角为β时风机不变桨风速范围的下边界与上边界;vL(β)与vU(β)分别满足:
vL(β)={v|0.5ρπR2v3CPLR/v,β)=Pcmd}
vU(β)={v|0.5ρπR2v3CPUR/v,β)=Pcmd}
只要风速波动不超出
Figure BDA0003414937000000026
则可由桨距角为β的风轮在变速区间[ωLU]内的被动变速完全应对,无需启动桨距角调节;否则,风轮转速将会达到ωL或ωU而触发转速边界处的恒转速变桨;
Figure BDA0003414937000000027
为风机在桨距角β下的不变桨风速范围。
进一步地,所述ωU设定为额定转速ωN,ωU=ωN,ωL为恒定功率曲线与最优功率曲线交点对应的转速
Figure BDA0003414937000000028
Figure BDA0003414937000000029
Figure BDA00034149370000000210
式中,
Figure BDA00034149370000000211
为风机最佳叶尖速比,
Figure BDA00034149370000000212
为最大风能利用系数。
进一步地,步骤3所述分析确定桨距角β对不变桨风速范围的影响关系,具体包括:
桨距角β对
Figure BDA0003414937000000031
的影响主要体现在风速水平和波动范围两个方面;
1)
Figure BDA0003414937000000032
内的风速水平随桨距角增大呈现升高趋势,桨距角越大,被动变速风轮可应对幅值越高的风速波动;
2)
Figure BDA0003414937000000033
随桨距角增大呈现扩大趋势,当桨距角增大时,风轮捕获气动功率对转速变化的敏感度升高,相同风速下的风轮转速变化量减小;因此,桨距角越大,相同变速区间的风轮被动变速可缓冲更大范围的风速波动。
进一步地,步骤4中所述被动变速风轮桨距角与风速变化范围Uv之间的适配关系具体包括:
实际风速变化范围Uv为[vmin,vmax],vmin与vmax分别为最低风速和最高风速,其与不变桨风速范围
Figure BDA0003414937000000034
之间的关系有6种,具体描述如下:
关系1:
Figure BDA0003414937000000035
且vU1)<vmin,即Uv整体高于
Figure BDA0003414937000000036
此时,若桨距角保持在β1,则风机稳定平衡点的转速大于ωU,风轮转速势必达到ωU而启动上调桨距角;
关系2:vmin<vL2)<vU2)<vmax,即
Figure BDA0003414937000000037
此时,若桨距角保持在β2,则风机运行至ωU或ωL而启动变桨调节;
关系3:vL3)<vmin<vU3)<vmax,即
Figure BDA0003414937000000038
此时,若桨距角保持在β3,则高于vU3)风速对应的稳定平衡点的转速大于ωU,风轮转速加速到ωU而启动上调桨距角;
关系4:vmin<vL4)<vmax<vU4),即
Figure BDA0003414937000000039
此时,若桨距角保持在β4,则低于vL4)风速对应的稳定平衡点的转速小于ωL,风轮转速减速到ωL而启动下调桨距角;
关系5:vL5)<vmin<vmax<vU5),即
Figure BDA00034149370000000310
此时,风机可以仅依靠该桨距角风轮被动变速独立应对风速波动,表明风轮转速不会达到ωU和ωL而启动变桨调节;
关系6:
Figure BDA0003414937000000041
且vmax<vL6),即Uv整体低于
Figure BDA0003414937000000042
此时,若桨距角保持在β6,则风机稳定平衡点的转速小于ωL,风轮转速势必达到ωL而启动下调桨距角;
根据
Figure BDA0003414937000000043
与Uv的关系,给出适配桨距角
Figure BDA0003414937000000044
确定策略:
Figure BDA0003414937000000045
式中,[βminmax]为桨距角运行范围,通常为[0°,35°];
Figure BDA0003414937000000046
Figure BDA0003414937000000047
与Uv交集
Figure BDA0003414937000000048
对应的连续有界区间的长度,其定义如下:
Figure BDA0003414937000000049
Figure BDA00034149370000000410
特别地,若
Figure BDA00034149370000000411
Figure BDA00034149370000000412
进一步地,步骤5所述基于适配关系实现考虑不变桨风速范围的风机有功功率控制,具体包括:
步骤5-1,应用基于牛顿-拉夫逊算法的风速观测方法,根据持续观测得到的风速序列获得风速波动范围Uv
步骤5-2,预设计算Uv的周期,将当前周期结束时刻得到的Uv作为紧邻下一周期的Uv估计值;
步骤5-3,根据Uv估计值,基于所述适配桨距角
Figure BDA00034149370000000413
确定策略,搜索紧邻下一周期的适配桨距角
Figure BDA00034149370000000414
进一步地,步骤5-3中所述搜索紧邻下一周期的适配桨距角
Figure BDA00034149370000000415
具体为:在桨距角运行范围[βminmax]内采用遍历算法搜索紧邻下一周期的适配桨距角
Figure BDA00034149370000000416
遍历搜索步长为0.1°。
一种考虑不变桨风速范围的风电机组有功功率控制系统,所述系统包括:
参数采集模块,用于获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
计算模块,用于计算风机的不变桨风速范围
Figure BDA0003414937000000051
第一关系建立模块,用于分析确定桨距角β对不变桨风速范围的影响关系;
第二关系建立模块,用于基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000053
确定策略;
控制模块,用于基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure BDA0003414937000000056
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000054
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure BDA0003414937000000052
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000055
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
本发明与现有技术相比,其显著优点为:1)相比于IAPC方法,新增的桨距角设定环节克服了被动变速风轮桨距角设定的随机性和不确定性;2)桨距角β对应的不变桨风速范围与湍流风速波动范围的适配性使得大惯量风轮能更好地缓冲高频率、小幅值的风速波动分量,减小因触及转速边界而引发的风轮超速或电磁功率跌落;3)降低了风轮达到转速边界的频次,进一步提升风机APC控制在风轮转速调节、变桨机构疲劳载荷和功率指令响应三个方面的性能。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为本发明考虑不变桨风速范围的风电机组有功功率控制的流程图。
图2为一个实施例中风机不变桨风速范围示意图。
图3为一个实施例中不同桨距角下的风机不变桨风速范围图。
图4为一个实施例中不变桨风速范围与风速波动范围之间关系示意图。
图5为一个实施例中本发明APC方法的控制原理框图。
图6为一个实施例中风机在不同APC方法下的轨迹图,其中图(a)至图(d)分别为不同APC方法下风速、转速、电磁功率和桨距角的轨迹图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
需要说明,若本发明实施例中有涉及“第一”、“第二”等的描述,则该“第一”、“第二”等的描述仅用于描述目的,而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外,各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
在一个实施例中,结合图1和图2,提供了一种考虑不变桨风速范围的风电机组有功功率控制方法,所述方法包括以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure BDA0003414937000000071
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000077
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
进一步地,在其中一个实施例中,结合图3,步骤2所述计算风机的不变桨风速范围
Figure BDA0003414937000000078
具体过程包括:
对于功率指令为Pcmd和变速区间设定为[ωLU]的有功功率控制,当桨距角为β时,存在一个最大的风速波动范围
Figure BDA0003414937000000072
Figure BDA0003414937000000073
使得对于该风速范围内的任一风速,风机都存在一个稳定平衡点(ω00),且满足ω0∈[ωLU]与β0=β;式中,ωL与ωU分别为变速区间的下边界与上边界,vL(β)与vU(β)分别为桨距角为β时风机不变桨风速范围的下边界与上边界;vL(β)与vU(β)分别满足:
vL(β)={v|0.5ρπR2v3CPLR/v,β)=Pcmd}
vU(β)={v|0.5ρπR2v3CPUR/v,β)=Pcmd}
只要风速波动不超出
Figure BDA0003414937000000079
则可由桨距角为β的风轮在变速区间[ωLU]内的被动变速完全应对,无需启动桨距角调节;否则,风轮转速将会达到ωL或ωU而触发转速边界处的恒转速变桨;
Figure BDA0003414937000000074
为风机在桨距角β下的不变桨风速范围。
进一步优选地,在其中一个实施例中,所述ωU设定为额定转速ωN,ωU=ωN,ωL为恒定功率曲线与最优功率曲线交点对应的转速
Figure BDA0003414937000000075
Figure BDA0003414937000000076
Figure BDA0003414937000000081
式中,
Figure BDA0003414937000000082
为风机最佳叶尖速比,
Figure BDA0003414937000000083
为最大风能利用系数。
进一步地,在其中一个实施例中,步骤3所述分析确定桨距角β对不变桨风速范围的影响关系,具体包括:
桨距角β对
Figure BDA0003414937000000084
的影响主要体现在风速水平和波动范围两个方面;
1)
Figure BDA0003414937000000085
内的风速水平随桨距角增大呈现升高趋势,桨距角越大,被动变速风轮可应对幅值越高的风速波动;
2)
Figure BDA0003414937000000086
随桨距角增大呈现扩大趋势,当桨距角增大时,风轮捕获气动功率对转速变化的敏感度升高,相同风速下的风轮转速变化量减小;因此,桨距角越大,相同变速区间的风轮被动变速可缓冲更大范围的风速波动。
进一步地,在其中一个实施例中,结合图4和图5,步骤4中所述被动变速风轮桨距角与风速变化范围Uv之间的适配关系具体包括:
实际风速变化范围Uv为[vmin,vmax],vmin与vmax分别为最低风速和最高风速,其与不变桨风速范围
Figure BDA0003414937000000087
之间的关系有6种,具体描述如下:
关系1:
Figure BDA0003414937000000088
且vU1)<vmin,即Uv整体高于
Figure BDA0003414937000000089
此时,若桨距角保持在β1,则风机稳定平衡点的转速大于ωU,风轮转速势必达到ωU而启动上调桨距角,容易造成风轮超速;
关系2:vmin<vL2)<vU2)<vmax,即
Figure BDA00034149370000000810
此时,若桨距角保持在β2,则风机有可能运行至ωU或ωL而启动变桨调节;
关系3:vL3)<vmin<vU3)<vmax,即
Figure BDA00034149370000000811
此时,若桨距角保持在β3,则高于vU3)风速对应的稳定平衡点的转速大于ωU,风轮转速有可能加速到ωU而启动上调桨距角;
关系4:vmin<vL4)<vmax<vU4),即
Figure BDA00034149370000000812
此时,若桨距角保持在β4,则低于vL4)风速对应的稳定平衡点的转速小于ωL,风轮转速有可能减速到ωL而启动下调桨距角;
关系5:vL5)<vmin<vmax<vU5),即
Figure BDA0003414937000000091
此时,风机可以仅依靠该桨距角风轮被动变速独立应对风速波动,表明风轮转速不会达到ωU和ωL而启动变桨调节;
关系6:
Figure BDA0003414937000000092
且vmax<vL6),即Uv整体低于
Figure BDA0003414937000000093
此时,若桨距角保持在β6,则风机稳定平衡点的转速小于ωL,风轮转速势必达到ωL而启动下调桨距角,容易导致电磁功率跌落;
根据
Figure BDA0003414937000000094
与Uv的关系,给出适配桨距角
Figure BDA0003414937000000095
确定策略:
Figure BDA0003414937000000096
式中,[βminmax]为桨距角运行范围,通常为[0°,35°];
Figure BDA0003414937000000097
Figure BDA0003414937000000098
与Uv交集
Figure BDA0003414937000000099
对应的连续有界区间的长度,其定义如下:
Figure BDA00034149370000000915
Figure BDA00034149370000000910
特别地,若
Figure BDA00034149370000000911
Figure BDA00034149370000000916
进一步地,在其中一个实施例中,步骤5所述基于适配关系实现考虑不变桨风速范围的风机有功功率控制,具体包括:
步骤5-1,应用基于牛顿-拉夫逊算法的风速观测方法,根据持续观测得到的风速序列获得风速波动范围Uv
步骤5-2,预设计算Uv的周期,将当前周期结束时刻得到的Uv作为紧邻下一周期的Uv估计值;
步骤5-3,根据Uv估计值,基于所述适配桨距角
Figure BDA00034149370000000912
确定策略,搜索紧邻下一周期的适配桨距角
Figure BDA00034149370000000913
这里,优选地,所述搜索紧邻下一周期的适配桨距角
Figure BDA00034149370000000914
具体为:在桨距角运行范围[βminmax]内采用遍历算法搜索紧邻下一周期的适配桨距角
Figure BDA0003414937000000101
遍历搜索步长为0.1°。
一种考虑不变桨风速范围的风电机组有功功率控制系统,所述系统包括:
参数采集模块,用于获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
计算模块,用于计算风机的不变桨风速范围
Figure BDA0003414937000000104
第一关系建立模块,用于分析确定桨距角β对不变桨风速范围的影响关系;
第二关系建立模块,用于基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000102
确定策略;
控制模块,用于基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
关于考虑不变桨风速范围的风电机组有功功率控制系统的具体限定可以参见上文中对于考虑不变桨风速范围的风电机组有功功率控制方法的限定,在此不再赘述。上述考虑不变桨风速范围的风电机组有功功率控制系统中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure BDA0003414937000000105
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000103
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
关于每一步的具体限定可以参见上文中对于考虑不变桨风速范围的风电机组有功功率控制方法的限定,在此不再赘述。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure BDA0003414937000000113
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure BDA0003414937000000112
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
关于每一步的具体限定可以参见上文中对于考虑不变桨风速范围的风电机组有功功率控制方法的限定,在此不再赘述。
作为一种具体示例,在其中一个实施例中,对本发明进行进一步验证说明。
本实施例利用美国国家可再生能源实验室(National Renewable EnergyLaboratory,NREL)提供的开源的专业风力机仿真软件FAST(Fatigue,Aerodynamics,Structures,and Turbulence)来仿真验证效果。风力机模型采用的是容量为600kW的CART3叶片试验机型,具体参数如下表1所示。
表1容量为600kW的CART 3叶片试验机型主要参数
Figure BDA0003414937000000111
本发明考虑不变桨风速范围的风电机组有功功率控制方法包括以下内容:
1、获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
2、根据风机气动特性和变速区间设定,计算风机的不变桨风速范围
Figure BDA0003414937000000121
具体包括:
对于功率指令为Pcmd和变速区间设定为[ωLU]的有功功率控制方法,当桨距角为β时,存在一个最大的风速波动范围
Figure BDA0003414937000000122
Figure BDA0003414937000000123
使得对于该风速范围内的任一风速,风机都存在一个稳定平衡点(ω00),且满足ω0∈[ωLU]与β0=β。式中,ωL与ωU分别为变速区间的下边界与上边界,vL(β)与vU(β)分别为桨距角为β时风机不变桨风速范围的下边界与上边界。vL(β)与vU(β)分别满足:
vL(β)={v|0.5ρπR2v3CPLR/v,β)=Pcmd}
vU(β)={v|0.5ρπR2v3CPUR/v,β)=Pcmd}
3、分析不变桨风速范围的影响规律,具体包括:
桨距角对
Figure BDA0003414937000000124
的影响主要体现在风速水平和波动范围两个方面。
1)
Figure BDA0003414937000000125
内的风速水平随桨距角增大呈现升高趋势。桨距角越大,被动变速风轮可以应对幅值越高的风速波动。
2)
Figure BDA0003414937000000126
随桨距角增大呈现扩大趋势。当桨距角增大时,风轮捕获气动功率对转速变化的敏感度升高,相同风速下的风轮转速变化量减小。因此,桨距角越大,相同变速区间的风轮被动变速可以缓冲更大范围的风速波动。
4、建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,具体包括:
根据
Figure BDA0003414937000000127
与Uv的关系,给出适配桨距角
Figure BDA0003414937000000128
确定策略:
Figure BDA0003414937000000129
式中,[βminmax]为桨距角运行范围,通常为[0°,35°]。
Figure BDA00034149370000001210
Figure BDA00034149370000001211
与Uv交集
Figure BDA0003414937000000131
对应的连续有界区间的长度,其定义如下:
Figure BDA0003414937000000132
Figure BDA0003414937000000133
特别地,若
Figure BDA0003414937000000134
Figure BDA0003414937000000135
5、基于适配关系实现考虑不变桨风速范围的风机有功功率控制,具体包括:
应用基于牛顿-拉夫逊算法的风速观测方法,根据持续观测得到的风速序列获得风速波动范围Uv。本发明计算Uv的周期设定为1分钟,将当前周期结束时刻得到的Uv作为紧邻下一周期的Uv估计值。
根据Uv估计值,基于所述的适配桨距角确定策略,搜索紧邻下一周期的适配桨距角
Figure BDA0003414937000000136
由于是有限范围内的单变量寻优,在桨距角运行范围[βminmax]内采用遍历算法,遍历搜索步长为0.1°。
6、在相同的仿真风速(600s湍流风速序列,平均风速为10m/s、湍流强度为A级)和功率指令(功率指令设定为150kW)场景设定下,仿真RSC方法、IAPC方法和本发明提出的改进方法,得到风机在三种方法下的控制性能指标如表2所示。风机在不同方法下的转速、桨距角和电磁功率轨迹如图6所示。
表2不同APC方法的评价指标
Figure BDA0003414937000000137
比较三种方法的APC控制性能指标。从表2和图6可以看出,在IAPC方法的基础上增加了考虑不变桨风速范围的桨距角设定环节,根据风速信息动态更新的适配桨距角提升了风轮被动变速对湍流风速的缓冲效果,有效降低了风轮达到转速边界的频次,缓解了风机超速和避免了电磁功率跌落。被动变速风轮能更好地缓冲高频率、小幅值的风速波动分量,改善了功率指令响应性能。
以上仿真实验结果表明,采用本发明所提出的考虑不变桨风速范围的风机有功功率控制方法能有效降低风轮达到转速边界的频次,进一步缓解了风轮超速和电磁功率跌落,验证了本发明所提出改进方法的有效性和实用性。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (10)

1.一种考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,所述方法包括以下步骤:
步骤1,获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
步骤2,计算风机的不变桨风速范围
Figure FDA0003414936990000011
步骤3,分析确定桨距角β对不变桨风速范围的影响关系;
步骤4,基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure FDA0003414936990000012
确定策略;
步骤5,基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
2.根据权利要求1所述的考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,步骤2所述计算风机的不变桨风速范围
Figure FDA0003414936990000013
具体过程包括:
对于功率指令为Pcmd和变速区间设定为[ωLU]的有功功率控制,当桨距角为β时,存在一个最大的风速波动范围
Figure FDA0003414936990000014
Figure FDA0003414936990000015
使得对于该风速范围内的任一风速,风机都存在一个稳定平衡点(ω00),且满足ω0∈[ωLU]与β0=β;式中,ωL与ωU分别为变速区间的下边界与上边界,vL(β)与vU(β)分别为桨距角为β时风机不变桨风速范围的下边界与上边界;vL(β)与vU(β)分别满足:
vL(β)={v|0.5ρπR2v3CPLR/v,β)=Pcmd}
vU(β)={v|0.5ρπR2v3CPUR/v,β)=Pcmd}
只要风速波动不超出
Figure FDA0003414936990000016
则可由桨距角为β的风轮在变速区间[ωLU]内的被动变速完全应对,无需启动桨距角调节;否则,风轮转速将会达到ωL或ωU而触发转速边界处的恒转速变桨;
Figure FDA0003414936990000017
为风机在桨距角β下的不变桨风速范围。
3.根据权利要求2所述的考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,所述ωU设定为额定转速ωN,ωU=ωN,ωL为恒定功率曲线与最优功率曲线交点对应的转速
Figure FDA0003414936990000021
Figure FDA0003414936990000022
Figure FDA0003414936990000023
式中,
Figure FDA0003414936990000024
为风机最佳叶尖速比,
Figure FDA0003414936990000025
为最大风能利用系数。
4.根据权利要求1或2所述的考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,步骤3所述分析确定桨距角β对不变桨风速范围的影响关系,具体包括:
桨距角β对
Figure FDA0003414936990000026
的影响主要体现在风速水平和波动范围两个方面;
1)
Figure FDA0003414936990000027
内的风速水平随桨距角增大呈现升高趋势,桨距角越大,被动变速风轮可应对幅值越高的风速波动;
2)
Figure FDA0003414936990000028
随桨距角增大呈现扩大趋势,当桨距角增大时,风轮捕获气动功率对转速变化的敏感度升高,相同风速下的风轮转速变化量减小;因此,桨距角越大,相同变速区间的风轮被动变速可缓冲更大范围的风速波动。
5.根据权利要求4所述的考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,步骤4中所述被动变速风轮桨距角与风速变化范围Uv之间的适配关系具体包括:
实际风速变化范围Uv为[vmin,vmax],vmin与vmax分别为最低风速和最高风速,其与不变桨风速范围
Figure FDA0003414936990000029
之间的关系有6种,具体描述如下:
关系1:
Figure FDA00034149369900000210
且vU1)<vmin,即Uv整体高于
Figure FDA00034149369900000211
此时,若桨距角保持在β1,则风机稳定平衡点的转速大于ωU,风轮转速势必达到ωU而启动上调桨距角;
关系2:vmin<vL2)<vU2)<vmax,即
Figure FDA00034149369900000212
此时,若桨距角保持在β2,则风机运行至ωU或ωL而启动变桨调节;
关系3:vL3)<vmin<vU3)<vmax,即
Figure FDA0003414936990000031
此时,若桨距角保持在β3,则高于vU3)风速对应的稳定平衡点的转速大于ωU,风轮转速加速到ωU而启动上调桨距角;
关系4:vmin<vL4)<vmax<vU4),即
Figure FDA0003414936990000032
此时,若桨距角保持在β4,则低于vL4)风速对应的稳定平衡点的转速小于ωL,风轮转速减速到ωL而启动下调桨距角;
关系5:vL5)<vmin<vmax<vU5),即
Figure FDA0003414936990000033
此时,风机可以仅依靠该桨距角风轮被动变速独立应对风速波动,表明风轮转速不会达到ωU和ωL而启动变桨调节;
关系6:
Figure FDA0003414936990000034
且vmax<vL6),即Uv整体低于
Figure FDA0003414936990000035
此时,若桨距角保持在β6,则风机稳定平衡点的转速小于ωL,风轮转速势必达到ωL而启动下调桨距角;
根据
Figure FDA0003414936990000036
与Uv的关系,给出适配桨距角
Figure FDA0003414936990000037
确定策略:
Figure FDA0003414936990000038
式中,[βminmax]为桨距角运行范围,通常为[0°,35°];
Figure FDA0003414936990000039
Figure FDA00034149369900000310
与Uv交集
Figure FDA00034149369900000311
对应的连续有界区间的长度,其定义如下:
Figure FDA00034149369900000312
Figure FDA00034149369900000313
特别地,若
Figure FDA00034149369900000314
Figure FDA00034149369900000315
6.根据权利要求5所述的考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,步骤5所述基于适配关系实现考虑不变桨风速范围的风机有功功率控制,具体包括:
步骤5-1,应用基于牛顿-拉夫逊算法的风速观测方法,根据持续观测得到的风速序列获得风速波动范围Uv
步骤5-2,预设计算Uv的周期,将当前周期结束时刻得到的Uv作为紧邻下一周期的Uv估计值;
步骤5-3,根据Uv估计值,基于所述适配桨距角
Figure FDA0003414936990000041
确定策略,搜索紧邻下一周期的适配桨距角
Figure FDA0003414936990000042
7.根据权利要求6所述的考虑不变桨风速范围的风电机组有功功率控制方法,其特征在于,步骤5-3中所述搜索紧邻下一周期的适配桨距角
Figure FDA0003414936990000043
具体为:在桨距角运行范围[βminmax]内采用遍历算法搜索紧邻下一周期的适配桨距角
Figure FDA0003414936990000044
遍历搜索步长为0.1°。
8.一种考虑不变桨风速范围的风电机组有功功率控制系统,其特征在于,所述系统包括:
参数采集模块,用于获取风机相关参数,包括空气密度ρ,风轮半径R,额定转速ωN,风机气动特性CP(λ,β),功率指令Pcmd、风轮转速ω;
计算模块,用于计算风机的不变桨风速范围
Figure FDA0003414936990000045
第一关系建立模块,用于分析确定桨距角β对不变桨风速范围的影响关系;
第二关系建立模块,用于基于所述影响关系建立被动变速风轮桨距角与风速变化范围Uv之间的适配关系,并基于该关系获得适配桨距角
Figure FDA0003414936990000046
确定策略;
控制模块,用于基于所述适配关系实现考虑不变桨风速范围的风机有功功率控制。
9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述方法的步骤。
CN202111543163.0A 2021-12-16 2021-12-16 考虑不变桨风速范围的风电机组有功功率控制方法及系统 Active CN114263564B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111543163.0A CN114263564B (zh) 2021-12-16 2021-12-16 考虑不变桨风速范围的风电机组有功功率控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111543163.0A CN114263564B (zh) 2021-12-16 2021-12-16 考虑不变桨风速范围的风电机组有功功率控制方法及系统

Publications (2)

Publication Number Publication Date
CN114263564A true CN114263564A (zh) 2022-04-01
CN114263564B CN114263564B (zh) 2024-02-13

Family

ID=80827521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111543163.0A Active CN114263564B (zh) 2021-12-16 2021-12-16 考虑不变桨风速范围的风电机组有功功率控制方法及系统

Country Status (1)

Country Link
CN (1) CN114263564B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116111653A (zh) * 2022-12-21 2023-05-12 天津大学 Prsv模式下考虑平衡点优化的风电apc系统及方法
CN116591895A (zh) * 2023-05-29 2023-08-15 国网江苏省电力有限公司电力科学研究院 一种风电机组有功功率控制方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103397983A (zh) * 2013-07-18 2013-11-20 上海交通大学 变速风机有功功率控制和转速控制方法
CN108843489A (zh) * 2018-05-04 2018-11-20 南京理工大学 基于限转速平滑功率控制的风机变桨优化方法
CN108843494A (zh) * 2018-05-04 2018-11-20 南京理工大学 基于斜线平滑功率控制的风机变桨优化方法
CN112664393A (zh) * 2020-12-22 2021-04-16 南京理工大学 基于最大不变桨风速范围指标的风机有功功率控制方法
CN113315174A (zh) * 2021-05-27 2021-08-27 哈电风能有限公司 风电场机组调控方法、装置、风电场控制终端及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103397983A (zh) * 2013-07-18 2013-11-20 上海交通大学 变速风机有功功率控制和转速控制方法
CN108843489A (zh) * 2018-05-04 2018-11-20 南京理工大学 基于限转速平滑功率控制的风机变桨优化方法
CN108843494A (zh) * 2018-05-04 2018-11-20 南京理工大学 基于斜线平滑功率控制的风机变桨优化方法
CN112664393A (zh) * 2020-12-22 2021-04-16 南京理工大学 基于最大不变桨风速范围指标的风机有功功率控制方法
CN113315174A (zh) * 2021-05-27 2021-08-27 哈电风能有限公司 风电场机组调控方法、装置、风电场控制终端及存储介质

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116111653A (zh) * 2022-12-21 2023-05-12 天津大学 Prsv模式下考虑平衡点优化的风电apc系统及方法
CN116111653B (zh) * 2022-12-21 2024-04-12 天津大学 Prsv模式下考虑平衡点优化的风电apc系统及方法
CN116591895A (zh) * 2023-05-29 2023-08-15 国网江苏省电力有限公司电力科学研究院 一种风电机组有功功率控制方法及系统

Also Published As

Publication number Publication date
CN114263564B (zh) 2024-02-13

Similar Documents

Publication Publication Date Title
CN108488035B (zh) 永磁直驱风力发电机组失速和变桨混合控制方法
CN114263564A (zh) 考虑不变桨风速范围的风电机组有功功率控制方法及系统
CN110552850B (zh) 一种基于风速提前测量的风电机组有功调节方法及装置
CN105649875B (zh) 风力发电机组的变桨控制方法及装置
BRPI1015099B1 (pt) Controlador para uma turbina eólica flutuante, estrutura de turbina eólica flutuante e método de controle de uma turbina eólica flutuante
BR102015015592B1 (pt) Sistema de turbina eólica e método de operação de uma turbina eólica
CN109185054A (zh) 风力发电机桨距角控制方法、系统、装置及可读存储介质
CN110345006A (zh) 一种风力发电机组低风速区最大功率跟踪优化控制方法
CN111379665B (zh) 风力发电机组的变桨控制方法及系统
CN106065848A (zh) 控制风力涡轮机的操作的方法
CN111209638B (zh) 基于运行攻角偏差修正的低风速风力机叶片气动设计方法
CN107194122B (zh) 一种变速风机叶片的多工况气动优化的改进方法
CN106777525B (zh) 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法
CN112664393B (zh) 基于最大不变桨风速范围指标的风机有功功率控制方法
CN108537372B (zh) 一种风向预测方法及风力发电机组的偏航控制方法
CN112682258B (zh) 基于Backstepping的大型风力机最大功率点跟踪控制方法
CN110030148A (zh) 基于风速提前测量的非线性预测变桨控制方法
CN114876732A (zh) 一种风电机组变桨的控制方法及装置
CN114294155A (zh) 风电机组有功功率控制方法及装置
CN113464378A (zh) 一种基于深度强化学习的提升风能捕获的转速跟踪目标优化方法
CN112211782A (zh) 基于起始转速自适应搜索的收缩跟踪区间控制方法
CN116591895A (zh) 一种风电机组有功功率控制方法及系统
Perez-Blanco et al. Aerodynamic performance of preferred wind turbine airfoils
CN115828569B (zh) 基于鲸鱼优化算法的双馈风电机组传动系统参数辨识方法
CN112211781B (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