提高风电机组弱风发电量的变桨控制方法

摘要

本发明提供了一种提高风电机组弱风发电量的变桨控制方法，所述方法是在恒最低转速区域通过提前桨距角动作以增加发电量，所提高发电量的恒最低转速区域用风速区间来描述指的是切入风速到最大功率点跟踪(MPPT)段起始风速之间这一弱风速区间。本发明在现有技术的基础上提高[Vcutin，v1]风速区间的发电量。本发明仅仅需要在现有变桨控制策略基础上做出适当修改，无需增加额外硬件装置，且可以一定程度降低叶根挥舞力矩和轴向推力。

权利要求

1.一种提高风电机组弱风发电量的变桨控制方法，其特征在于，所述方法在恒最低转速区域通过提前桨距角动作增加发电量，所提高发电量的风速区间是切入风速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加上常规控制桨距角指令作为优化后的桨距角指令β*，由风机主控发送至变桨执行机构。

2.根据权利要求1所述的一种提高风电机组弱风发电量的变桨控制方法，其特征在于，第二步中，所述待定系数c1和c2通过离线辨识方法得到，具体为：在专业风电仿真软件中对[Vcutin，v1]风速区间内某一风速分别设定不同桨距角大小，比较不同桨距角下的功率输出，输出功率最大所对应的桨距角即为该风速下的βopt；选取[Vcutin，v1]风速区间内的其他风速点采用同样的方法分别确定出各自对应的最优变桨角度βopt，通过参数拟合得到待定系数c1和c2。

3.根据权利要求1所述的一种提高风电机组弱风发电量的变桨控制方法，其特征在于，第三步中，在计算最优变桨角度βopt时，为避免桨距角频繁动作，需要对检测到的风速信号进行低通滤波处理得到V，对于计算得到的最优变桨角度βopt需要设定指令死区即通过指令死区模块。

4.根据权利要求1-3任一项所述的一种提高风电机组弱风发电量的变桨控制方法，其特征在于，所述指令死区模块的功能由受指令死区阈值βdz控制的保持器实现，如果保持器输出βopt1与最优变桨角度βopt差值的绝对值大于βdz，则保持器输出当前的最优变桨角度βopt；否则保持器处于保持状态，输出上一时刻的βopt1。

5.根据权利要求1所述的一种提高风电机组弱风发电量的变桨控制方法，其特征在于，所述变桨控制方法是一种基于功率查表的变桨控制方法，包括如下步骤：第一步，检测发电机功率，对检测到的发电机功率信号进行低通滤波，若发电机功率大于切入风速对应的稳态功率Pvcutin且小于风速v1对应的稳态功率Pv1，且发电机转速ω满足|ω-ωmin|<ωL，ωL设置为0.04倍的标幺值，发电机转速的基准值为发电机同步转速，ωmin为风电机组发电运行的最低转速，执行以下步骤；

第二步，根据发电机功率查表得到最优变桨角度βopt；

第三步，对第二步得到的最优变桨角度βopt经过指令死区模块输出桨距角指令值βopt1；

第四步，将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*，由风机主控发送至变桨执行机构。

6.根据权利要求5所述的一种提高风电机组弱风发电量的变桨控制方法，其特征在于，所述最优变桨角度βopt，采用以下方法确定：检测风速，对检测到的风速信号进行低通滤波得到风速V，根据公式βopt＝c1(v1-V)2+c2(v1-V)，在确定待定系数c1和c2、并已知[Vcutin，v1]风速区间内的发电机稳态功率曲线的基础上，选取[Vcutin，v1]风速区间内一定数量的离散风速点，根据与这些风速点一一对应的发电机稳态功率和最优变桨角度βopt制作一张关于发电机稳态功率和最优变桨角度βopt一一对应的表；假设Vcutin风速下稳态功率为Pvcutin，v1风速下稳态功率为Pv1，则表中对于不在[Pvcutin Pv1]区间内的功率所对应的最优变桨角度βopt为0；实现时根据实时测量得到的发电机功率通过查表和线性插值获取最优变桨角度βopt。

7.根据权利要求6所述的一种提高风电机组弱风发电量的变桨控制方法，其特征在于，第二步中，所述待定系数c1和c2通过离线辨识方法得到，具体为：在专业风电仿真软件中对[Vcutin，v1]风速区间内某一风速分别设定不同桨距角大小，比较不同桨距角下的功率输出，输出功率最大所对应的桨距角即为该风速下的βopt；选取[Vcutin，v1]风速区间内的其他风速点采用同样的方法分别确定出各自对应的最优变桨角度βopt，通过参数拟合得到待定系数c1和c2。

说明书

提高风电机组弱风发电量的变桨控制方法

技术领域

[0001] 本发明涉及风力发电技术领域，具体地，涉及一种恒最低转速区域通过提前桨距角动作增加发电量的控制方法。

背景技术

[0002] 双馈风电机组受限于变流器电流可控性，通常将机组转差率限定在+0.3之内，因此在弱风速段存在一个最低转速ωmin的限制。现有风电运营商普遍采用这一方案，风力机在切入风速Vcutin到MPPT(最大功率点跟踪)段起始风速v1这一风速区间运行在恒定的最低转速ωmin。现有技术在[Vcutin，v1]风速区间将桨距角设为某一恒定值(通常为0度)，桨距角在到达额定风速前不予动作，这种运行控制策略导致了机组发电量在[Vcutin，v1]风速区间内存在一定提升空间。

[0003] 检索中发现现有技术有不少关于风电机组的变桨控制方法，比如中国专利申请号为201410171828.3的发明专利，该专利公开一种变速恒频风电机组变速变桨距联合控制方法，“该控制方法引入最优转速估计值以及转速参考值分别作为转速控制与变桨距控制的参考量，转速PID控制器根据最优转速估计值与风电机组额定转速的偏差进行转矩指令调节，变桨距PID控制器根据转速参考值与机组额定转速的偏差进行桨距角指令调节，实现转速控制与变桨距控制之间的解耦控制，避免转速控制与变桨距控制之间的干扰。当风速在额定值附近波动、风电机组捕获的气动功率大于额定功率时，气动功率估计值对应的转速参考值大于风电机组额定转速，变桨距PID控制器在转速PID控制器输出的转矩指令达到额定值之前调节桨距角指令，变速变桨超前联合控制，抑制机组转速和输出功率的超调。”[0004] 但是上述技术仍旧没有考虑到机组发电量在[Vcutin，v1]风速区间内的提升问题。

发明内容

[0005] 本发明的目的是提供一种提高风电机组弱风发电量的变桨控制方法，该方法在恒最低转速区域通过提前桨距角动作增加发电量，提高了风力机气动转矩，从而在所述风速区域提高发电量。

[0006] 为实现以上目的，本发明提供一种提高风电机组弱风发电量的变桨控制方法，所述方法在恒最低转速区域通过提前桨距角动作增加发电量，仅需对常规控制方法的变桨控制做小幅修改，变流器控制则无需改动，变桨控制的改进通过修改软件完成，无需增加硬件。具体而言，所提高发电量的风速区间是切入风速Vcutin到MPPT段起始风速v1之间这一弱风速区间，所适用的机组是在切入风速Vcutin到MPPT段起始风速v1之间限定最低转速的风电机组。

[0007] 本发明所述变桨控制方法是一种基于风速检测的变桨控制方法，包括如下步骤：

[0008] 第一步，检测风速，对检测到的风速信号进行低通滤波得到风速V，若风速V处于[Vcutin，v1]风速区间，且发电机转速ω满足|ω-ωmin|<ωL，ωL设置为0.04倍的标幺值，发电机转速的基准值为发电机同步转速，ωmin为风电机组发电运行的最低转速，执行以下步骤；

[0009] 第二步，根据第一步得到的风速V，按照以下公式计算最优变桨角度βopt：

[0010] βopt＝c1(v1-V)2+c2(v1-V)，

[0011] 其中：待定系数c1和c2通过辨识得到，v1为起始风速；

[0012] 第三步，对第二步得到的最优变桨角度βopt经过指令死区模块输出桨距角指令值βopt1；

[0013] 第四步，将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*，由风机主控发送至变桨执行机构。

[0014] 本发明在[Vcutin，v1]风速区间根据上述βopt＝c1(v1-V)2+c2(v1-V)公式提前动作桨距角，从而提高[Vcutin，v1]风速区间的气动转矩，并提高发电量。

[0015] 优选地，第二步中，所述待定系数c1和c2通过离线辨识方法得到，具体为：在专业风电仿真软件，如Bladed中，针对采用实际机组参数的特定机组模型，对[Vcutin，v1]风速区间内某一风速分别设定不同桨距角大小，比较不同桨距角下的功率输出，输出功率最大所对应的桨距角即为该风速下的βopt；选取[Vcutin，v1]风速区间内的其他风速点采用同样的方法分别确定出各自对应的最优变桨角度βopt，通过参数拟合得到待定系数c1和c2。

[0016] 优选地，第三步中，在计算最优变桨角度βopt时，为避免桨距角频繁动作，需要对检测到的风速信号进行低通滤波处理得到V，对于计算得到的最优变桨角度βopt需要设定指令死区即通过指令死区模块。所述指令死区模块的功能由受指令死区阈值βdz控制的保持器实现，如果保持器输出βopt1与最优变桨角度βopt差值的绝对值大于βdz，则保持器输出当前的最优变桨角度βopt；否则保持器处于保持状态，输出上一时刻的βopt1。

[0017] 本发明所述变桨控制也可以采用基于功率查表的方法实现，包括如下步骤：

[0018] 第一步，检测发电机功率，对检测到的发电机功率信号进行低通滤波，若发电机功率大于切入风速Vcutin对应的稳态功率Pvcutin且小于风速v1对应的稳态功率Pv1，且发电机转速ω满足|ω-ωmin|<ωL，ωL设置为0.04倍的标幺值，发电机转速的基准值为发电机同步转速，执行以下步骤；

[0019] 第二步，根据发电机功率查表得到最优变桨角度βopt；

[0020] 第三步，对第二步得到的最优变桨角度βopt经过指令死区模块输出桨距角指令值βopt1。

[0021] 第四步，将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*，由风机主控发送至变桨执行机构。

[0022] 所述第二步，根据发电机功率查表得到最优变桨角度βopt，具体为：在确定待定系数c1和c2、并已知[Vcutin，v1]风速区间内的发电机稳态功率曲线的基础上，选取[Vcutin，v1]风速区间内一定数量的离散风速点，根据与这些风速点一一对应的发电机稳态功率和最优变桨角度βopt制作一张关于发电机稳态功率和最优变桨角度βopt一一对应的表，表中对于不在Vcutin对应稳态功率和v1对应稳态功率之间的功率对应的最优变桨角度βopt为0；实现时根据实时测量得到的发电机功率通过查表和线性插值获取最优变桨角度βopt。

[0023] 与现有技术相比，本发明具有如下的有益效果：

[0024] 本发明在现有技术的基础上提高[Vcutin，v1]风速区间的发电量。本发明仅仅需要在现有变桨控制策略基础上做出适当修改，无需增加额外硬件装置，且可以一定程度降低叶根挥舞力矩和轴向推力。

附图说明

[0025] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：

[0026] 图1为本发明一实施例所适用机组在弱风速段的典型运行控制策略；

[0027] 图2为本发明一实施例实现方法框图；

[0028] 图3为本发明一实施例指令死区模块示意图；

[0029] 图4为本发明一实施例在3m/s风速下的验证；

[0030] 图5为本发明一实施例在4m/s风速下的验证；

[0031] 图6为本发明一实施例在动态风速下的验证，其中：图(a)为实施案例所采用的风速，图(b)为桨距角，图(c)为轴向推力，图(d)为电功率。

具体实施方式

[0032] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。

[0033] 如图1所示，本发明是对现有常规控制方法在机组特定运行范围内的优化和补充，本发明的基本内容是在所述恒最低转速区域通过提前桨距角动作增加发电量的控制方法，仅需对常规控制方法的变桨控制做小幅修改，变流器控制则无需改动，变桨控制的改进通过修改软件完成，无需增加硬件。具体而言，所提高发电量的风速区间是切入风速Vcutin到MPPT(最大功率跟踪)段起始风速v1之间这一弱风速区间，所适用的机组是在切入风速Vcutin到MPPT段起始风速v1之间限定最低转速的风电机组。特别说明本发明提高的是风力机这一原动力的出力，适用对象是根据风力机是否在切入风速Vcutin到MPPT段起始风速v1之间限定最低转速来决定的，与发电机类型及电气拓扑无关。

[0034] 以下将就本发明方案的原理进行阐述。

[0035] 风电机组最大功率跟踪就是要最大限度获取风力机气动功率，由于机组在[Vcutin，v1]风速区间恒转速运行，则[Vcutin，v1]风速区间的最大气动功率跟踪问题可转化为最大气动转矩跟踪问题。水平轴风力机的气动转矩Ta可以用切向力系数Ct来描述，

[0036]

[0037] ρ为空气密度，B为桨叶个数，W为叶素上的合成风速，l为桨叶弦长，r为叶素所在半径，r0为轮毂半径，R为叶尖处半径，dr是积分微元，其中，

[0038] Ct＝Cl(i)sinI-Cd(i)cos I  (2)

[0039] 式中Cl为升力系数，Cd为阻力系数，I为入流角，Cl、Cd随攻角i变化，某典型机组升、阻力系数随攻角i变化的曲线可近似表示如下(考虑到风力机在正常发电运行状态攻角的变化范围，这里仅给出了攻角i在0～10度之间的规律)：

[0040]

[0041] Cl＝kcli；0＜i＜i1  (4)

[0042] kcd是阻力系数随攻角i变化的变化率，kcl是升力系数随攻角i变化的变化率。其中Cd0的典型值为0.005，kcd的典型值为0.07878，kcl的典型值是7.21927，i0的典型值是0.0384弧度，i1的典型值是0.1745弧度。上述典型值来源于某5MW机组。

[0043]

[0044] i＝I-β  (6)

[0045] 式中β为桨距角，a为大型风电机组轴向诱导因子，b为大型风电机组周向诱导因子，ω为风力机转速，λ为叶尖速比，V为风速。

[0046] 由于v1是MPPT段和恒最低转速段的分界点，需要借助风速为v1时的转矩系数Ctv1进行分析，假定风速v1时采用常规控制在气动稳态下a和b分别为av1和b v1，记该状态下的入流角

[0047]

[0048] 若采用常规机组控制，桨距角恒定为β0(通常为0度)，记风速为v1时的Ctv1如式(8)所示

[0049]

[0050] 式(7)、(8)考虑到入流角I很小，从而sinI≈I，cosI≈1，并据此做出了简化，式中Clv1和Cdv1分别是风速v1下的升力系数和阻力系数。

[0051] 风速变化、转速变化和桨距角变化均会导致a和b变化，但是三者的小幅变化对a和b变化的影响是微弱的。则在v1附近足够小的邻域内可假定a和b保持与风速v1时相同。现假设一个略小于v1的风速v0＝v1-dv，假设dv足够小，则此时入流角I可以简化为[0052]

[0053] 攻角

[0054] i＝I-β≈kv0-β  (10)

[0055] 若采用常规机组控制在[Vcutin，v1]风速区间将桨距角设为0，考虑正常运行状态中i>i0，则风速v0下的切向力系数Ct可计算如下：

[0056]

[0057] 根据式(2)、(6)，在v1附近足够小的风速左邻域内：

[0058]

[0059] 由于v1对应于MPPT段，该风速下常规桨距角(0度)为最优值，则有

[0060]

[0061] 由于kcd、kcl及k均大于0，根据式(13)，则v0略小于v1时，

[0062] kclkv0-kcd＜0  (14)

[0063] 可见式(11)等式右边第二项在[vcutin,v1)风速区间对切向力系数有正面作用。若在[vcutin,v1)风速区间采用变桨控制，假设变桨角度Δβ，则

[0064]

[0065] dv代表风速小幅增量。

[0066] 可见桨距角正调一定角度可以提高式(11)等式右边第二项的正面作用，从而提高切向力系数达到增大气动转矩增加发电量的目的。根据式(10)，桨距角增大则i减小，由于i

[0067] 假设小于v1的某一特定风速v01下的最优桨距角为β01，并假设该风速下对应于式(7)中的系数k发生一定量变后变化为k’，则

[0068]

[0069] 则在略小于v01的风速v02，

[0070]

[0071] 式(16)、(17)表明，在风速v01最优桨距角β01的基础上进一步增大桨距角可以增加风速v02下的Ct，即随风速减小βopt增大。

[0072] 以上对本方案的原理性陈述结果表明：1.在[vcutin,v1]风速区间内，合理的正调桨可以提高发电量，该结论表明了变桨的方向；2.[vcutin,v1]风速区间内某一特定风速下存在一个对应的最优桨距角βopt；3.随风速减小，特定风速下的βopt单调增大。

[0073] 此外经过专业风电软件验证，本发明方法在增加发电量的同时可以一定程度上减小叶根挥舞弯矩和轴向推力，并对机舱和塔架水平方向的载荷和位移有改善作用。

[0074] 上述结论表明了变桨的方向和最优变桨角度随风速的单调性关系，属于定性的说明，特定风速下的最优变桨角度βopt仍难以直接定量。即便如此，根据上述定性结论可以做出如下合理假定(考虑到风速为v1时最优桨距角为0)：

[0075] βopt＝c1(v1-V)2+c2(v1-V)  (18)

[0076] 式(18)中的待定系数c1和c2可以辨识得到。一种离线辨识方法如下：在专业风电仿真软件(如Bladed)中对[Vcutin，v1]风速区间内某一风速分别设定不同桨距角大小，比较不同桨距角下的功率输出，输出功率最大所对应的桨距角即为该风速下的βopt。选取[Vcutin，v1]风速区间内的其他风速点采用同样的方法分别确定出各自对应的βopt，建议获取4组以上V和βopt对应的数据，通过参数拟合得到式(18)中的c1和c2。

[0077] 本发明在[Vcutin，v1]风速区间根据式(18)提前动作桨距角，从而提高[Vcutin，v1]风速区间的气动转矩，并提高发电量。根据式(18)确定βopt时，为避免桨距角频繁动作，需要对检测到的风速信号进行低通滤波处理得到V，对于式(18)计算得到的βopt需要设定指令死区。如图2所示，本方案的实现方法流程。

[0078] 指令死区模块设计如下，其功能由受指令死区阈值βdz控制的保持器实现，如果保持器输出βopt1与βopt差值的绝对值大于βdz，则保持器输出当前的βopt；否则保持器处于保持状态，输出上一时刻的βopt1。如图3所示。

[0079] 另外考虑到风速测量准确性问题，在对机组特性十分清楚的情况下，提供以下替代计算βopt的方法：在确定式(18)中待定系数并已知[Vcutin，v1]风速区间内的发电机稳态功率曲线的基础上，选取[Vcutin，v1]风速区间内一定数量的离散风速点(建议10组以上)，根据与这些风速点一一对应的发电机稳态功率和βopt制作一张关于发电机稳态功率和βopt一一对应的表(表中对于不在Vcutin对应稳态功率和v1对应稳态功率之间的功率对应的βopt为0)。实现时根据实时测量得到的发电机功率通过查表和线性插值获取βopt。

[0080] 基于上述分析，本发明方案是对常规方法在机组特定运行范围内的优化和补充。

[0081] 根据所述的研究内容，在离线确定式(18)中待定系数c1和c2的基础上，本发明基于风速检测的变桨控制方法为：

[0082] 第一步.检测风速，对检测到的风速信号进行低通滤波得到风速V，若风速V处于[Vcutin，v1]风速区间，执行以下步骤。

[0083] 第二步.根据第一步得到的风速V，按照式(17)计算桨距角指令值βopt；

[0084] 第三步.对第二步得到的桨距角指令值βopt经过指令死区模块输出βopt1。

[0085] 第四步.将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*，由风机主控发送至变桨执行机构。

[0086] 本发明基于功率查表的变桨控制方法为：

[0087] 第一步.检测发电机功率，对检测到的发电机功率信号进行低通滤波，若发电机功率大于切入风速对应的稳态功率Pvcutin且小于风速v1对应的稳态功率Pv1，执行以下步骤；

[0088] 第二步.根据发电机功率查表得到βopt；

[0089] 第三步.对第二步得到的桨距角指令值βopt经过指令死区模块输出βopt1。

[0090] 第四步.将第三步得到的桨距角指令值βopt1加上常规控制桨距角指令作为优化后的桨距角指令β*，由风机主控发送至变桨执行机构。

[0091] 运用本实施例所述方法在一台商用1.5MW机组上进行多组仿真验证，结果证明可提高1％左右发电量，所验证机组的切入风速为3m/s，v1为5m/s，风力机最低转速1.13rad/s，桨叶长度43m。

[0092] 在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。

[0093] 此后仿真中[Vcutin，v1]风速区间的βopt由式(18)确定，部分仿真结果如下：

[0094] 仿真选取风速为3m/s的情况下，根据式(18)计算得到βopt为2.3度，采用本发明所述方法与常规方法(桨距角0度)的对比如图4所示，可见该风速下发电量相比于常规控制提高20％左右，提高的功率约为5kW。

[0095] 仿真选取风速为4m/s的情况下，根据式(18)计算得到βopt为1.5度，采用本发明所述方法与常规方法(桨距角0度)的对比如图5所示，可见该风速下发电量相比于常规控制提高2.5％左右，提高的功率约为2.5kW。

[0096] 本发明方案也在动态风速下进行了验证，验证中仍然采用式(18)计算βopt，采用基于风速检测的实现方法进行验证。选取恒最低转速段风速如图6(a)所示，采用本发明所提出的技术与常规控制仿真对比结果(虚线代表常规方法，实线代表本发明方法)如图6所示。

[0097] 本发明在现有技术的基础上提高[Vcutin，v1]风速区间的发电量。本发明仅仅需要在现有变桨控制策略基础上做出适当修改，无需增加额外硬件装置，且可以一定程度降低叶根挥舞力矩和轴向推力。

[0098] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。