基于激光测风前馈的风电机组复合智能变桨控制方法与流程
本发明涉及风电机组技术领域，更具体地说，涉及一种基于激光测风前馈的风电机组复合智能变桨控制方法。
背景技术：
传统机械式测风仪不能提前测量叶轮正前方的风速，并且受叶片转动影响存在测量不准的缺陷；目前风电机组的变桨控制算法是变增益pi控制算法，利用转速误差反馈来计算桨距角的目标指令值，而转速误差滞后于风速的变化，这样桨距角并不能根据风速的变化及时动作，不利于降低机组的载荷；并且pi控制算法针对外界的风速扰动抗扰性较弱；现有的基于激光测风雷达的风速前馈控制大部分是与变增益pi反馈控制相结合的复合控制，而变增益pi反馈控制本身存在缺陷。
技术实现要素：
针对现有方法的不足，提出一种基于激光测风前馈的风电机组复合智能变桨控制方法。
本发明解决其技术问题所采用的技术方案是：构造一种基于激光测风前馈的风电机组复合智能变桨控制方法，包括：
采集风电机组叶片前方同一平面内四点的风速信息，重构风电机组轮毂高度的等效风速，同时通过获取风电机组的实时转速信息，计算实时的转速误差；
基于风电机组轮毂高度的等效风速，计算风电机组前馈控制器前后时刻的输出增量，将风电机组前馈控制器的输出增量及所述转速误差作为rbf神经网络的输入，输出结果为前馈控制器优化参数的调整变化量；
依据前馈控制器优化参数的调整变化量对前馈控制器的优化参数进行调整；
通过曲线拟合的方法，获取风电机组桨距角与状态观测器观测带宽之间的函数关系，通过调整状态观测器观测带宽，实现对风电机组自抗扰变桨控制器参数的调整；
将自抗扰控制器的输出与前馈控制器的输出相加得到对风电机组变桨执行机构的控制指令，实现对风电机组桨距角的调整。
其中，风电机组轮毂高度的等效风速表示为：
其中，vi(t)表示四点中任意一点的风速，λi表示各点的系数，取值范围在0～1。
其中，在计算风电机组前馈控制器的输出增量的步骤中，
风电机组前馈控制器在k时刻的输出信号βff(k)为从当前时刻到之前n-1时刻输入风速的线性加权和，即：
将转速误差对前馈控制器输出的敏感度
前馈控制器的输出增量为：
其中，在rbf神经网络中，网络的输入向量为x＝[x1，x2，…，xn]t，网络隐含层节点径向基向量为h＝[h1，h2，…，hm]t，其中hm为高斯函数。网络隐含层第j个节点的中心矢量为cj＝[cj1，cj2，…，cjm]t，隐含层节点基宽向量为b＝[b1，b2，…，bm]t，网络的权向量为w＝[w1，w2，…，wm]t，设第k时刻辨识网络输出为em(k)，rbf神经网络辨识输出如公式5所示：
em(k)＝w1h1+w2h2+…+wmhm(5)
rbf网络辨识性能指标函数如公式6所示：
采用梯度下降法对rbf神经网络各参数进行迭代，为了使搜索过程收敛于全局极小值，引入了包含k时刻之前两个时刻系数变化量的动量项，如公式7～11所示：
wj(k)＝wj(k-1)+η(e(k)-em(k))hj+α1(wj(k-1)-
wj(k-2))+α2(wj(k-2)-wj(k-3))(7)
bj(k)＝bj(k-1)+ηδbj+α1(bj(k-1)-bj(k-2))+
α2(bj(k-2)-bj(k-3))(9)
cji(k)＝cji(k-1)+ηδcji+α1(cji(k-1)-cji(k-2))+α2(cji(k-
2)-cji(k-3))(11)
其中动量因子α1、α2和学习速率η的取值范围为0～1；
将转速误差及前馈控制器输出增量输入rbf神经网络，输出结果为前馈控制器优化参数的调整变化量，公式表示为：
其中，η0、η1、……、ηn-1分别表示前馈滤波器各项系数的学习速率；
将前馈控制器优化参数的调整变化量对应添加至前馈控制器的待优化参数上，实现对前馈控制器的优化参数的实时调整。
其中，在通过曲线拟合的方法，获取风电机组桨距角与状态观测器观测带宽之间的函数关系，通过调整状态观测器观测带宽的步骤中，包括步骤：
桨距角与转速之间的动态特性可以用状态空间表达式表示为：
其中d(δω，x2，f)表示系统内部和外部的总扰动；
状态观测器的数学表达式如公式17所示：
设计ladrc控制器进行控制：
u0＝kp(r-y)(18)
ladrc控制器中需要整定的参数有kp、b0、β1和β2，其中，β1和β2与扩张状态观测器带宽wo存在如公式19所示的关系：
做不同幅度的风速阶跃试验，按照调整时间短，转速超调量小的原则，固定kp进行b0和wo参数的调整，b0、wo与桨距角θ的拟合函数关系如公式20所示：
通过拟合函数的方法，实现状态观测器带宽随桨距角的变化而变化，增强了观测器带宽的适应性。
区别于现有技术，本发明提供了一种基于激光测风前馈的风电机组复合智能变桨控制方法，利用激光测风雷达重构叶片前方轮毂高度的风速，根据转速误差和前馈控制量额变化，利用径向基神经网络(rbf)动态调整有限脉冲滤波器(fir)的系数，实现前馈控制算法参数的动态自适应调整；通过曲线拟合的方法，得到了桨距角与状态观测器观测带宽之间的函数关系，设计了变带宽的ladrc变桨距控制策略弥补了传统pi变桨控制算法的缺陷，将rbf-fir与ladrc两者结合形成了能够自适应调节的复合智能变桨控制算法。通过本发明，实现了控制参数的自适应调节，增强了变桨控制算法的鲁棒性和抗扰性，实现了风电机组降低转速波动和机组叶根载荷的目的。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法的流程示意图。
图2是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法的控制结构示意图。
图3是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法与现有方法的转速调节时间和超调分析对比图。
图4是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法的与现有方法的转速波动分析对比图。
图5是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法的与现有方法的功率波动分析对比图。
图6是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法的与现有方法的叶根x方向上载荷力矩分析对比图。
图7是本发明提供的一种基于激光测风前馈的风电机组复合智能变桨控制方法的与现有方法的叶根y方向上载荷力矩转速分析对比图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
参阅图1和图2，本发明提供了一种基于激光测风前馈的风电机组复合智能变桨控制方法，包括：
采集风电机组叶片前方同一平面内四点的风速信息，重构风电机组轮毂高度的等效风速，同时通过获取风电机组的实时转速信息，计算实时的转速误差；
基于风电机组轮毂高度的等效风速，计算风电机组前馈控制器的输出增量，将风电机组前馈控制器的输出增量及所述转速误差作为rbf神经网络的输入，输出结果为前馈控制器优化参数的调整变化量；
依据前馈控制器优化参数的调整变化量对前馈控制器的优化参数进行调整；
通过曲线拟合的方法，获取风电机组桨距角与状态观测器观测带宽之间的函数关系，通过调整状态观测器观测带宽，实现对风电机组自抗扰变桨控制器参数的调整；
将自抗扰控制器的输出与前馈控制器的输出相加得到对风电机组变桨执行机构的控制指令，实现对风电机组桨距角的调整。
本发明利用激光测风雷达的测量数据，设计一种抗扰性、鲁棒性强的变桨控制算法，是风电机组向着智能化发展的必然要求。本文利用激光测风雷达测量叶片前方50m处轮毂高度的风速，根据转速误差和前馈控制量额变化，利用径向基神经网络(rbf)动态调整有限脉冲滤波器(fir)的系数，实现前馈控制算法参数的动态自适应调整；通过曲线拟合的方法，得到了桨距角与状态观测器观测带宽之间的函数关系，设计了变带宽的ladrc变桨距控制策略弥补了传统pi变桨控制算法的缺陷，将rbf-fir与ladrc两者结合形成了能够自适应调节的复合智能变桨控制算法。
激光测风雷达通过采集同一平面内v1、v2、v3、v4四个点的风速信息，常用式1来重构风电机组轮毂高度的等效风速。
有限脉冲滤波器在k时刻的输出信号βff(k)为从当前时刻到之前n-1时刻输入风速的线性加权和，即：
在rbf神经网络中，网络的输入向量为x＝[x1，x2，…，xn]t，网络隐含层节点径向基向量为h＝[h1，h2，…，hm]t，其中hm为高斯函数。网络隐含层第j个节点的中心矢量为cj＝[cj1，cj2，…，cjm]t，隐含层节点基宽向量为b＝[b1，b2，…，bm]t，网络的权向量为w＝[w1，w2，…，wm]t，设第k时刻辨识网络输出为em(k)，rbf神经网络辨识输出如公式5所示：
em(k)＝w1h1+w2h2+…+wmhm(5)
rbf网络辨识性能指标函数如公式6所示：
采用梯度下降法对rbf神经网络各参数进行迭代，为了使搜索过程收敛于全局极小值，引入了包含k时刻之前两个时刻系数变化量的动量项，如公式7～11所示：
wj(k)＝wj(k-1)+η(e(k)-em(k))hj+α1(wj(k-1)-
wj(k-2))+α2(wj(k-2)-wj(k-3))(7)
bj(k)＝bj(k-1)+ηδbj+α1(bj(k-1)-bj(k-2))+
α2(bj(k-2)-bj(k-3))(9)
cji(k)＝cji(k-1)+ηδcji+α1(cji(k-1)-cji(k-2))+α2(cji(k-
2)-cji(k-3))(11)
其中动量因子α1、α2和学习速率η的取值范围为0～1；
rbf能够根据性能指标函数公式6进行在线的参数自我整定，同时能够给出前馈控制器的参数整定算法，使前馈控制器实现了参数的自适应调整。前馈控制的目标是保证在出现风速波动扰动的情况下，系统转速误差最小，所以取前馈控制的性能指标函数为：
将转速误差对前馈控制器的控制输出的敏感度
前馈控制器的输出增量为：
将转速误差及前馈控制器输出增量输入rbf神经网络，输出结果为前馈控制器优化参数的调整变化量，公式表示为：
其中，η0、η1、……、ηn-1分别表示前馈滤波器各项系数的学习速率；利用负梯度方向决定每一次迭代新的搜索方向，每次迭代都能保证前馈控制性能指标函数向着最小的方向来优化参数
桨距角与转速之间的动态特性可以用状态空间表达式表示为：
其中d(δω，x2，f)表示系统内部和外部的总扰动。在此基础上设计了扩张状态观测器(eso)，观测器的数学表达式如式17所示：
经过补偿后的对象可以看做是积分串联的对象，设计比例控制器进行控制：
u0＝kp(r-y)(18)
ladrc控制器中需要整定的参数有kp、b0、β1和β2，但是β1和β2与扩张状态观测器带宽wo存在如式19所示的关系，所以实际控制器只需要整定kp、b0和wo3个参数。
做不同幅度的风速阶跃试验，按照调整时间短，转速超调量小的原则，固定kp进行b0和wo参数的调整。式20给出了b0和wo与桨距角θ的拟合函数关系。
通过拟合函数的方法，能够实现观测器带宽随桨距角的变化而变化，增强了观测器带宽的适应性。
采用美国可再生能源实验室(nrel)开发的开源软件fast中的5mw双馈风电机组模型对所提的控制方法进行测试，对比了pi控制、rbfnnfir+pi控制、ladrc控制和rbfnnfir+ladrc控制四种变桨控制策略的控制效果。
为了验证系统的鲁棒性，分别针对1m/s，2m/s，3m/s，4m/s，5m/s，6m/s的阶跃风速，通过对系统转速超调和转速稳定的调节时间进行统计分析。
分析结果如图3所示，四种控制策略均能保证在不同阶跃风速下的控制稳定，相比较而言，本发明提出的控制策略在转速调节时间和转速超调量方面都有明显的下降趋势，故控制的鲁棒性更好。
按照iec61400标准，取中等湍流强度b级，分别采用四种控制策略，对转速、功率波动的最大值、最小值、均值、方差，以及叶根和y方向上载荷力矩的最大值、最小值、均值、方差进行统计分析，如图4-7所示。
从图4-7不难看出各个指标在四种控制策略下的变化情况，转速、功率方差以及叶根方向力矩的均值都具有明显的降低，但是方向力矩的波动无明显变化。叶根y方向力矩的均值和方差在四种控制策略下的控制效果基本持平。经过统计分析，本发明提出的复合智能变桨控制算法能够降低转速波动约47.9％，功率波动约34.3％，降低叶根方向力矩均值约8.8％。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
- 还没有人留言评论。精彩留言会获得点赞！