一种大型风电机组的鲸鱼群PID独立变桨控制方法技术

本发明专利技术涉及一种大型风电机组的鲸鱼群PID独立变桨控制方法，包括如下步骤：通过叶片载荷模型，计算出各叶片在垂直于叶轮面方向的弯矩；将各叶片垂直于叶轮面方向的弯矩进行MBC坐标变换，得到coleman坐标系下的俯仰弯矩和偏航弯矩；通过鲸鱼群PID控制器算法求解得到分别用于消除俯仰弯矩和偏航弯矩的影响角；从而得到各叶片的桨距角增量，再与统一变桨角度进行叠加，得到各叶片的独立变桨桨距角后，通过变桨驱动器完成独立变桨动作。通过对鲸鱼群算法进行优化，在PID参数整定精度、稳定性方面具有良好的动态性能，有效降低机组的不平衡载荷，减少机组的振动情况，保证输出功率更加稳定平滑，提升机组的寿命。

A PID independent pitch control method for large wind turbines

The invention relates to a whale group PID independent pitch control method for a large wind turbine, which comprises the following steps: calculating the bending moment of each blade in the direction perpendicular to the impeller surface through the blade load model; transforming the bending moment of each blade in the direction perpendicular to the impeller surface through MBC coordinate transformation to obtain the pitch bending moment and yaw bending moment in the Coleman coordinate system; calculating the bending moment through the whale group PID controller algorithm The solution is used to eliminate the influence angle of pitch and yaw moment respectively, so as to obtain the pitch angle increment of each blade, and then it is superposed with the unified pitch angle to obtain the independent pitch angle of each blade, and then the independent pitch action is completed by the pitch driver. Through the optimization of whale group algorithm, it has good dynamic performance in PID parameter setting accuracy and stability, effectively reduces the unbalanced load of the unit, reduces the vibration of the unit, ensures the output power is more stable and smooth, and improves the life of the unit.

【技术实现步骤摘要】
一种大型风电机组的鲸鱼群PID独立变桨控制方法
本专利技术涉及风力发电领域，具体涉及一种大型风电机组的独立变桨控制方法。
技术介绍
风能是一种安全、生态的可再生能源，变桨距风电机组对风能的利用效率高，占据了市场的主导地位。变桨距控制主要有统一变桨和独立变桨两种方式。统一变桨控制是假设风轮旋转平面内的风速是均匀的，但随着风轮直径的增加，风切变、塔影效应、偏航误差等因素使风轮面受力的不均衡度增加，附加载荷也越来越大，严重影响机组的安全运行和部件的使用寿命。目前，风能市场的快速增长要求整机厂家能够生产具有更大风轮直径且额定功率更高的风电机组。因此，减少负载成为了这种风电机组设计中具有挑战性的研究课题。近年来，提出了许多降低风机载荷，特别是风机叶片载荷的技术。其中，独立变桨控制可在限制功率的同时改善机组的不均衡载荷，是降低风机叶片载荷最有效的技术之一。由于PID控制算法结构简单，在工程控制中具有很高的实用性，因此，现有的风电机组独立变桨控制多采用PID控制算法。PID控制算法研究的重点在于PID控制器的参数整定方法，PID参数整定问题是一个典型的多模态优化问题，最优解很难求解。而现有的PID参数整定方法对于该多模态优化问题还没有很好的处理方案，不能确保解的质量。因此，亟需找到一种能尽量保证PID参数解的质量的方法，从而使控制系统有着良好的动态性能，降低机组的载荷，减少机组的振动情况，保证输出功率更加稳定平滑。WSA-IC(WhaleSwarmAlgorithmwithIterativeCounter)算法(详见参考文献：BingZeng,XinyuLi,LiangGao.etal.NeuralComput&Applic(2019).https://doi.org/10.1007/s00521-018-3949-4)是一种新型元启发式算法，有如下优点：1)不需要针对不同的问题设置小生境参数；2)可以在迭代过程中有效识别并跳出已经找到的极值点，从而节省不必要的函数评价，提高全局搜索能力；3)算法特定的参数不需要针对不同问题手动设置；4)种群大小不需要与最优解的个数精确匹配。鉴于WSA-IC算法的诸多优点，可以通过对WSA-IC算法进行改进，提高算法的求解精度、收敛速度和稳定性，并结合鲸鱼群算法用于PID控制器的参数整定。
技术实现思路
为了解决上述技术问题，本专利技术提供了一种具有自寻优、高质量解的大型风电机组的鲸鱼群PID独立变桨控制方法。本专利技术的一种大型风电机组的鲸鱼群PID独立变桨控制方法，包括以下步骤：步骤1、计算弯矩，通过叶片载荷模型，计算出各叶片在垂直于叶轮面方向的弯矩；步骤2、计算俯仰弯矩和偏航弯矩，根据步骤1得到的各叶片垂直于叶轮面方向的弯矩进行MBC(muti-bladecoordinate)坐标变换，得到coleman坐标系下的俯仰弯矩和偏航弯矩；步骤3、计算消除俯仰弯矩的影响角和消除偏航弯矩的影响角，通过鲸鱼群PID控制器求解得到coleman坐标系下分别用于消除步骤2得到的俯仰弯矩和偏航弯矩的影响角，所述鲸鱼群PID控制器中鲸鱼群算法计算每条鲸鱼个体的目标函数值f的函数模型如下：式中，比例系数kp和积分时间常数ki为待整定的两个PID参数，gm表示系统的幅值裕度，pm表示系统的相位裕度，st表示系统的调节时间，os表示系统的超调量；gth表示预先指定的最小幅值裕度，pth表示预先指定的最小相位裕度，sth表示期望的调节时间；[lp,up]表示kp的定义域，[li,ui]表示ki的定义域；步骤4、计算独立变桨桨距角，将步骤3得到的影响角进行MBC逆变换，得到各叶片的桨距角增量，再与统一变桨角度进行叠加，得到各叶片的独立变桨桨距角；步骤5、完成独立变桨，根据步骤4得到的独立变桨桨距角，各变桨驱动器完成独立变桨执行动作。优选的，所述步骤1中，通过叶片载荷模型，计算出各叶片在垂直于叶轮面方向的弯矩的公式如下所示：M1＝Mf,1·cos(θ1)-Me,1·sin(θ1)M2＝Mf,2·cos(θ2)-Me,2·sin(θ2)M3＝Mf,3·cos(θ3)-Me,3·sin(θ3)式中，M1、M2、M3分别表示1号、2号、3号叶片在垂直于叶轮面方向的弯矩，Mf,1、Mf,2、Mf,3分别表示1号、2号、3号叶片在面外方向flapwise的根部弯矩，Me,1、Me,2、Me,3分别表示1号、2号、3号叶片在面内方向edgewise的根部弯矩，θ1、θ2、θ3分别表示1号、2号、3号叶片的桨距角。优选的，所述步骤2中，风电机组在coleman坐标系下的俯仰弯矩Mtilt和偏航弯矩Myaw的计算公式如下所示：Mtilt＝M1·cos(ψ1)+M2·cos(ψ2)+M3·cos(ψ3)Myaw＝M1·sin(ψ1)+M2·sin(ψ2)+M3·sin(ψ3)式中，ψ1、ψ2和ψ3分别表示叶片1、叶片2和叶片3的方位角，其中，ψ2＝ψ1+120，ψ3＝ψ1+240。优选的，所述步骤3中，coleman坐标系下用于消除俯仰弯矩的影响角θtilt和用于消除偏航弯矩的影响角θyaw的计算方法如下所示：式中，比例系数kp和积分时间常数ki是待整定的两个PID参数。优选的，对PID控制器参数kp和ki整定的鲸鱼群算法包括以下步骤：Step1、设置鲸鱼群算法参数，其中包括种群大小p、迭代次数n、稳定性阈值Ts和邻域搜索半径r；Step2、初始化每条鲸鱼个体Ωi＝(kpi,kii)，得到初始鲸鱼种群Ω；Step3、计算每条鲸鱼个体的目标函数值f，将当前全局最优适应度值fgbest设置为0；Step4、设置编号i＝1；Step5、寻找鲸鱼Ωi的“较优且最近”的鲸鱼Y，若Y存在，则执行Step6，否则执行Step11；“较优且最近”的鲸鱼Y是当前鲸鱼的引导个体，是所有优于当前鲸鱼的鲸鱼中，离当前鲸鱼最近的鲸鱼；Step6、生成鲸鱼Ωi的副本X，X在Y的引导下根据如下位置更新公式进行移动，并执行Step7，式中，和分别指鲸鱼X的第i个元素在第t步与t+1步迭代时的位置；指鲸鱼Y的第i个元素在第t步迭代时的位置；dX,Y指X与Y之间的距离；表示产生0到之间的一个随机数；Step7、计算X的目标函数值f(X)，若f(X)大于f(Ωi)，执行Step8；否则执行Step9；Step8、将X赋值给Ωi，将鲸鱼Ωi的迭代计数器Ωi.c设置为0，执行Step17；Step9、如果鲸鱼Ωi的迭代计数器Ωi.c小于Ts，将Ωi.c加1，执行Step17；否则执行Step10；Step10、重新初始化鲸鱼Ωi，并计算鲸鱼Ωi的目标函数值f(Ωi)，执行Step17；Step11、生成鲸鱼Ωi的副本X′，对副本X′执行邻域搜索，执行Step12；Step12、计算邻域搜索后的目标函数值f(X′)，若f(X′)大于f(Ωi)，执行Step13；否则执行Step14；Step13、将X′赋值给Ωi，将鲸鱼Ωi的迭代计数器Ωi.c设置为0，执行Step17；Step14、如果鲸鱼Ωi的迭代计数器Ωi.c小于Ts，将Ωi.c加1，执行Step17；否则执行Step15；Step15、如果鲸鱼Ωi的目标函数值f(Ωi)大于fgbest，先将fgbes本文档来自技高网...

【技术保护点】
1.一种大型风电机组的鲸鱼群PID独立变桨控制方法，其特征在于包括以下步骤：步骤1、计算弯矩，通过叶片载荷模型，计算出各叶片在垂直于叶轮面方向的弯矩；步骤2、计算俯仰弯矩和偏航弯矩，根据步骤1得到的各叶片垂直于叶轮面方向的弯矩进行MBC(muti‑blade coordinate)坐标变换，得到coleman坐标系下的俯仰弯矩和偏航弯矩；步骤3、计算消除俯仰弯矩的影响角和消除偏航弯矩的影响角，通过鲸鱼群PID控制器求解得到coleman坐标系下分别用于消除步骤2得到的俯仰弯矩和偏航弯矩的影响角，所述鲸鱼群PID控制器中鲸鱼群算法计算每条鲸鱼个体的目标函数值f的函数模型如下：

【技术特征摘要】
1.一种大型风电机组的鲸鱼群PID独立变桨控制方法，其特征在于包括以下步骤：步骤1、计算弯矩，通过叶片载荷模型，计算出各叶片在垂直于叶轮面方向的弯矩；步骤2、计算俯仰弯矩和偏航弯矩，根据步骤1得到的各叶片垂直于叶轮面方向的弯矩进行MBC(muti-bladecoordinate)坐标变换，得到coleman坐标系下的俯仰弯矩和偏航弯矩；步骤3、计算消除俯仰弯矩的影响角和消除偏航弯矩的影响角，通过鲸鱼群PID控制器求解得到coleman坐标系下分别用于消除步骤2得到的俯仰弯矩和偏航弯矩的影响角，所述鲸鱼群PID控制器中鲸鱼群算法计算每条鲸鱼个体的目标函数值f的函数模型如下：式中，比例系数kp和积分时间常数ki为待整定的两个PID参数，gm表示系统的幅值裕度，pm表示系统的相位裕度，st表示系统的调节时间，os表示系统的超调量；gth表示预先指定的最小幅值裕度，pth表示预先指定的最小相位裕度，sth表示期望的调节时间；[lp,up]表示kp的定义域，[li,ui]表示ki的定义域；步骤4、计算独立变桨桨距角，将步骤3得到的影响角进行MBC逆变换，得到各叶片的桨距角增量，再与统一变桨角度进行叠加，得到各叶片的独立变桨桨距角；步骤5、完成独立变桨，根据步骤4得到的独立变桨桨距角，各变桨驱动器完成独立变桨执行动作。2.如权利要求1所述的一种大型风电机组的鲸鱼群PID独立变桨控制方法，其特征在于，所述步骤1中，通过叶片载荷模型，计算出各叶片在垂直于叶轮面方向的弯矩的公式如下所示：M1＝Mf,1·cos(θ1)-Me,1·sin(θ1)M2＝Mf,2·cos(θ2)-Me,2·sin(θ2)M3＝Mf,3·cos(θ3)-Me,3·sin(θ3)式中，M1、M2、M3分别表示1号、2号、3号叶片在垂直于叶轮面方向的弯矩，Mf,1、Mf,2、Mf,3分别表示1号、2号、3号叶片在面外方向flapwise的根部弯矩，Me,1、Me,2、Me,3分别表示1号、2号、3号叶片在面内方向edgewise的根部弯矩，θ1、θ2、θ3分别表示1号、2号、3号叶片的桨距角。3.如权利要求2所述的一种大型风电机组的鲸鱼群PID独立变桨控制方法，其特征在于，所述步骤2中，风电机组在coleman坐标系下的俯仰弯矩Mtilt和偏航弯矩Myaw的计算公式如下所示：Mtilt＝M1·cos(ψ1)+M2·cos(ψ2)+M3·cos(ψ3)Myaw＝M1·sin(ψ1)+M2·sin(ψ2)+M3·sin(ψ3)式中，ψ1、ψ2和ψ3分别表示叶片1、叶片2和叶片3的方位角，其中，ψ2＝ψ1+120，ψ3＝ψ1+240。4.如权利要求3所述的一种大型风电机组的鲸鱼群PID独立变桨控制方法，其特征在于，所述步骤3中，coleman坐标系下用于消除俯仰弯矩的影响角θtilt和用于消除偏航弯矩的影响角θyaw的计算方法如下所示：式中，比例系数kp和积分时间常数ki是待整定的两个PID参数。5.如权利要求4所述的一种大型风电机组的鲸鱼群PID独立变桨控制方法，其特征在于，对PID控制器参数kp和ki整定的鲸鱼群算法包括以下步骤：Step1、设置鲸鱼群算法参数，其中包...

【专利技术属性】
技术研发人员：曾冰，黄凌翔，张硕望，
申请(专利权)人：湘电风能有限公司，湘潭电机股份有限公司，
类型：发明
国别省市：湖南,43