发明内容
本发明的目的是提供一种提高风电机组弱风发电量的变桨控制方法,该方法在恒最低转速区域通过提前桨距角动作增加发电量,提高了风力机气动转矩,从而在所述风速区域提高发电量。
为实现以上目的,本发明提供一种提高风电机组弱风发电量的变桨控制方法,所述方法在恒最低转速区域通过提前桨距角动作增加发电量,仅需对常规控制方法的变桨控制做小幅修改,变流器控制则无需改动,变桨控制的改进通过修改软件完成,无需增加硬件。具体而言,所提高发电量的风速区间是切入风速Vcutin到MPPT段起始风速v1之间这一弱风速区间,所适用的机组是在切入风速Vcutin到MPPT段起始风速v1之间限定最低转速的风电机组。
本发明所述变桨控制方法是一种基于风速检测的变桨控制方法,包括如下步骤:
第一步,检测风速,对检测到的风速信号进行低通滤波得到风速V,若风速V处于[Vcutin,v1]风速区间,且发电机转速ω满足|ω-ωmin|<ωL,ωL设置为0.04倍的标幺值,发电机转速的基准值为发电机同步转速,ωmin为风电机组发电运行的最低转速,执行以下步骤;
第二步,根据第一步得到的风速V,按照以下公式计算最优变桨角度βopt:
βopt=c1(v1-V)2+c2(v1-V),
其中:待定系数c1和c2通过辨识得到,v1为起始风速;
第三步,对第二步得到的最优变桨角度βopt经过指令死区模块输出桨距角指令值βopt1;
第四步,将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*,由风机主控发送至变桨执行机构。
本发明在[Vcutin,v1]风速区间根据上述βopt=c1(v1-V)2+c2(v1-V)公式提前动作桨距角,从而提高[Vcutin,v1]风速区间的气动转矩,并提高发电量。
优选地,第二步中,所述待定系数c1和c2通过离线辨识方法得到,具体为:在专业风电仿真软件,如Bladed中,针对采用实际机组参数的特定机组模型,对[Vcutin,v1]风速区间内某一风速分别设定不同桨距角大小,比较不同桨距角下的功率输出,输出功率最大所对应的桨距角即为该风速下的βopt;选取[Vcutin,v1]风速区间内的其他风速点采用同样的方法分别确定出各自对应的最优变桨角度βopt,通过参数拟合得到待定系数c1和c2。
优选地,第三步中,在计算最优变桨角度βopt时,为避免桨距角频繁动作,需要对检测到的风速信号进行低通滤波处理得到V,对于计算得到的最优变桨角度βopt需要设定指令死区即通过指令死区模块。所述指令死区模块的功能由受指令死区阈值βdz控制的保持器实现,如果保持器输出βopt1与最优变桨角度βopt差值的绝对值大于βdz,则保持器输出当前的最优变桨角度βopt;否则保持器处于保持状态,输出上一时刻的βopt1。
本发明所述变桨控制也可以采用基于功率查表的方法实现,包括如下步骤:
第一步,检测发电机功率,对检测到的发电机功率信号进行低通滤波,若发电机功率大于切入风速Vcutin对应的稳态功率Pvcutin且小于风速v1对应的稳态功率Pv1,且发电机转速ω满足|ω-ωmin|<ωL,ωL设置为0.04倍的标幺值,发电机转速的基准值为发电机同步转速,执行以下步骤;
第二步,根据发电机功率查表得到最优变桨角度βopt;
第三步,对第二步得到的最优变桨角度βopt经过指令死区模块输出桨距角指令值βopt1。
第四步,将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*,由风机主控发送至变桨执行机构。
所述第二步,根据发电机功率查表得到最优变桨角度βopt,具体为:在确定待定系数c1和c2、并已知[Vcutin,v1]风速区间内的发电机稳态功率曲线的基础上,选取[Vcutin,v1]风速区间内一定数量的离散风速点,根据与这些风速点一一对应的发电机稳态功率和最优变桨角度βopt制作一张关于发电机稳态功率和最优变桨角度βopt一一对应的表,表中对于不在Vcutin对应稳态功率和v1对应稳态功率之间的功率对应的最优变桨角度βopt为0;实现时根据实时测量得到的发电机功率通过查表和线性插值获取最优变桨角度βopt。
与现有技术相比,本发明具有如下的有益效果:
本发明在现有技术的基础上提高[Vcutin,v1]风速区间的发电量。本发明仅仅需要在现有变桨控制策略基础上做出适当修改,无需增加额外硬件装置,且可以一定程度降低叶根挥舞力矩和轴向推力。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1所示,本发明是对现有常规控制方法在机组特定运行范围内的优化和补充,本发明的基本内容是在所述恒最低转速区域通过提前桨距角动作增加发电量的控制方法,仅需对常规控制方法的变桨控制做小幅修改,变流器控制则无需改动,变桨控制的改进通过修改软件完成,无需增加硬件。具体而言,所提高发电量的风速区间是切入风速Vcutin到MPPT(最大功率跟踪)段起始风速v1之间这一弱风速区间,所适用的机组是在切入风速Vcutin到MPPT段起始风速v1之间限定最低转速的风电机组。特别说明本发明提高的是风力机这一原动力的出力,适用对象是根据风力机是否在切入风速Vcutin到MPPT段起始风速v1之间限定最低转速来决定的,与发电机类型及电气拓扑无关。
以下将就本发明方案的原理进行阐述。
风电机组最大功率跟踪就是要最大限度获取风力机气动功率,由于机组在[Vcutin,v1]风速区间恒转速运行,则[Vcutin,v1]风速区间的最大气动功率跟踪问题可转化为最大气动转矩跟踪问题。水平轴风力机的气动转矩Ta可以用切向力系数Ct来描述,
ρ为空气密度,B为桨叶个数,W为叶素上的合成风速,l为桨叶弦长,r为叶素所在半径,r0为轮毂半径,R为叶尖处半径,dr是积分微元,其中,
Ct=Cl(i)sinI-Cd(i)cosI(2)
式中Cl为升力系数,Cd为阻力系数,I为入流角,Cl、Cd随攻角i变化,某典型机组升、阻力系数随攻角i变化的曲线可近似表示如下(考虑到风力机在正常发电运行状态攻角的变化范围,这里仅给出了攻角i在0~10度之间的规律):
Cl=kcli;0<i<i1(4)
kcd是阻力系数随攻角i变化的变化率,kcl是升力系数随攻角i变化的变化率。其中Cd0的典型值为0.005,kcd的典型值为0.07878,kcl的典型值是7.21927,i0的典型值是0.0384弧度,i1的典型值是0.1745弧度。上述典型值来源于某5MW机组。
i=I-β(6)
式中β为桨距角,a为大型风电机组轴向诱导因子,b为大型风电机组周向诱导因子,ω为风力机转速,λ为叶尖速比,V为风速。
由于v1是MPPT段和恒最低转速段的分界点,需要借助风速为v1时的转矩系数Ctv1进行分析,假定风速v1时采用常规控制在气动稳态下a和b分别为av1和bv1,记该状态下的入流角
若采用常规机组控制,桨距角恒定为β0(通常为0度),记风速为v1时的Ctv1如式(8)所示
式(7)、(8)考虑到入流角I很小,从而sinI≈I,cosI≈1,并据此做出了简化,式中Clv1和Cdv1分别是风速v1下的升力系数和阻力系数。
风速变化、转速变化和桨距角变化均会导致a和b变化,但是三者的小幅变化对a和b变化的影响是微弱的。则在v1附近足够小的邻域内可假定a和b保持与风速v1时相同。现假设一个略小于v1的风速v0=v1-dv,假设dv足够小,则此时入流角I可以简化为
攻角
i=I-β≈kv0-β(10)
若采用常规机组控制在[Vcutin,v1]风速区间将桨距角设为0,考虑正常运行状态中i>i0,则风速v0下的切向力系数Ct可计算如下:
根据式(2)、(6),在v1附近足够小的风速左邻域内:
由于v1对应于MPPT段,该风速下常规桨距角(0度)为最优值,则有
由于kcd、kcl及k均大于0,根据式(13),则v0略小于v1时,
kclkv0-kcd<0(14)
可见式(11)等式右边第二项在[vcutin,v1)风速区间对切向力系数有正面作用。若在[vcutin,v1)风速区间采用变桨控制,假设变桨角度Δβ,则
dv代表风速小幅增量。
可见桨距角正调一定角度可以提高式(11)等式右边第二项的正面作用,从而提高切向力系数达到增大气动转矩增加发电量的目的。根据式(10),桨距角增大则i减小,由于i<i0后Cd不再减小(kcd=0),且变桨存在非线性特性,可以得出特定风速下存在一个最优变桨角度的结论。
假设小于v1的某一特定风速v01下的最优桨距角为β01,并假设该风速下对应于式(7)中的系数k发生一定量变后变化为k’,则
则在略小于v01的风速v02,
式(16)、(17)表明,在风速v01最优桨距角β01的基础上进一步增大桨距角可以增加风速v02下的Ct,即随风速减小βopt增大。
以上对本方案的原理性陈述结果表明:1.在[vcutin,v1]风速区间内,合理的正调桨可以提高发电量,该结论表明了变桨的方向;2.[vcutin,v1]风速区间内某一特定风速下存在一个对应的最优桨距角βopt;3.随风速减小,特定风速下的βopt单调增大。
此外经过专业风电软件验证,本发明方法在增加发电量的同时可以一定程度上减小叶根挥舞弯矩和轴向推力,并对机舱和塔架水平方向的载荷和位移有改善作用。
上述结论表明了变桨的方向和最优变桨角度随风速的单调性关系,属于定性的说明,特定风速下的最优变桨角度βopt仍难以直接定量。即便如此,根据上述定性结论可以做出如下合理假定(考虑到风速为v1时最优桨距角为0):
βopt=c1(v1-V)2+c2(v1-V)(18)
式(18)中的待定系数c1和c2可以辨识得到。一种离线辨识方法如下:在专业风电仿真软件(如Bladed)中对[Vcutin,v1]风速区间内某一风速分别设定不同桨距角大小,比较不同桨距角下的功率输出,输出功率最大所对应的桨距角即为该风速下的βopt。选取[Vcutin,v1]风速区间内的其他风速点采用同样的方法分别确定出各自对应的βopt,建议获取4组以上V和βopt对应的数据,通过参数拟合得到式(18)中的c1和c2。
本发明在[Vcutin,v1]风速区间根据式(18)提前动作桨距角,从而提高[Vcutin,v1]风速区间的气动转矩,并提高发电量。根据式(18)确定βopt时,为避免桨距角频繁动作,需要对检测到的风速信号进行低通滤波处理得到V,对于式(18)计算得到的βopt需要设定指令死区。如图2所示,本方案的实现方法流程。
指令死区模块设计如下,其功能由受指令死区阈值βdz控制的保持器实现,如果保持器输出βopt1与βopt差值的绝对值大于βdz,则保持器输出当前的βopt;否则保持器处于保持状态,输出上一时刻的βopt1。如图3所示。
另外考虑到风速测量准确性问题,在对机组特性十分清楚的情况下,提供以下替代计算βopt的方法:在确定式(18)中待定系数并已知[Vcutin,v1]风速区间内的发电机稳态功率曲线的基础上,选取[Vcutin,v1]风速区间内一定数量的离散风速点(建议10组以上),根据与这些风速点一一对应的发电机稳态功率和βopt制作一张关于发电机稳态功率和βopt一一对应的表(表中对于不在Vcutin对应稳态功率和v1对应稳态功率之间的功率对应的βopt为0)。实现时根据实时测量得到的发电机功率通过查表和线性插值获取βopt。
基于上述分析,本发明方案是对常规方法在机组特定运行范围内的优化和补充。
根据所述的研究内容,在离线确定式(18)中待定系数c1和c2的基础上,本发明基于风速检测的变桨控制方法为:
第一步.检测风速,对检测到的风速信号进行低通滤波得到风速V,若风速V处于[Vcutin,v1]风速区间,执行以下步骤。
第二步.根据第一步得到的风速V,按照式(17)计算桨距角指令值βopt;
第三步.对第二步得到的桨距角指令值βopt经过指令死区模块输出βopt1。
第四步.将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*,由风机主控发送至变桨执行机构。
本发明基于功率查表的变桨控制方法为:
第一步.检测发电机功率,对检测到的发电机功率信号进行低通滤波,若发电机功率大于切入风速对应的稳态功率Pvcutin且小于风速v1对应的稳态功率Pv1,执行以下步骤;
第二步.根据发电机功率查表得到βopt;
第三步.对第二步得到的桨距角指令值βopt经过指令死区模块输出βopt1。
第四步.将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*,由风机主控发送至变桨执行机构。
运用本实施例所述方法在一台商用1.5MW机组上进行多组仿真验证,结果证明可提高1%左右发电量,所验证机组的切入风速为3m/s,v1为5m/s,风力机最低转速1.13rad/s,桨叶长度43m。
在Bladed中经过反复试凑,得出风速3.5m/s时的βopt为1.9度、风速3.8m/s时的βopt为1.6度、4.2m/s时的βopt为1.1度、4.5m/s时的βopt为0.5度。根据这四组数据离线确定出式(18)中的待定系数为c1=-0.2052,c2=1.572。
此后仿真中[Vcutin,v1]风速区间的βopt由式(18)确定,部分仿真结果如下:
仿真选取风速为3m/s的情况下,根据式(18)计算得到βopt为2.3度,采用本发明所述方法与常规方法(桨距角0度)的对比如图4所示,可见该风速下发电量相比于常规控制提高20%左右,提高的功率约为5kW。
仿真选取风速为4m/s的情况下,根据式(18)计算得到βopt为1.5度,采用本发明所述方法与常规方法(桨距角0度)的对比如图5所示,可见该风速下发电量相比于常规控制提高2.5%左右,提高的功率约为2.5kW。
本发明方案也在动态风速下进行了验证,验证中仍然采用式(18)计算βopt,采用基于风速检测的实现方法进行验证。选取恒最低转速段风速如图6(a)所示,采用本发明所提出的技术与常规控制仿真对比结果(虚线代表常规方法,实线代表本发明方法)如图6所示。
本发明在现有技术的基础上提高[Vcutin,v1]风速区间的发电量。本发明仅仅需要在现有变桨控制策略基础上做出适当修改,无需增加额外硬件装置,且可以一定程度降低叶根挥舞力矩和轴向推力。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。