一种风力发电机组叶尖速比获取方法技术

本发明专利技术公开了一种风力发电机组叶尖速比获取方法，包括以下步骤：构造传动链观测器，通过使用卡尔曼观测器，从最终观测结果中获得原传动链系统的真实状态；通过传感器获得叶轮轮毂转速以及风力发电机组的发电机转矩，得到叶轮实际吸收的气动转矩；建立风力发电机的气动转矩与风力发电机的叶尖速比的映射关系；通过二次曲面的散点拟合多项式函数，得到叶尖速比。本发明专利技术通过对风力发电机物理系统本质特性的分析，直接减少了风速获取和风速估计部分，可在风力发电机控制系统实时运行环境下获得风电场风力发电机组的有效叶尖速比，有助于提高风力发电机组控制精度提升风机整机性能，增强风机控制系统的实时性、有效性和可靠性。

A Method for Obtaining Tip Speed Ratio of Wind Turbine

【技术实现步骤摘要】
一种风力发电机组叶尖速比获取方法
本专利技术涉及风电领域，特别涉及一种风力发电机组叶尖速比获取方法。
技术介绍
风能是绿色能源，风力发电对保护生态环境和改善能源结构具有重要意义。风力发电机组是一个复杂非线性多变量系统，控制系统方案的实时性与可靠性是有效保证风力发电机组安全高效运行的关键。实际运行中，不同风电场在地形地貌、机型、布置和风频分布等千差万别，发电量和安全性是决定风电场风力发电机组运行特性的最重要因素。如何在实现最大效率吸收风能的同时保证风力发电机安全平稳是风力发电机的基本设计问题。为提高风力发电机最大功率跟踪性能，学者们提出基于叶尖速比的最大功率跟踪控制器方法，该方法的基本思想是，放弃以前使用静态增益进行转矩跟踪而是将最大功率点跟踪问题转化为最优叶尖速比跟踪问题。同时，叶尖速比也可作为保证风力发电机安全性的一个基础指标。目前最佳叶尖速比法作为实现最大功率跟踪主要方法之一，传统的设计思路为要获得叶尖速比就应当先获得可靠的风速，在风速测量的基础上计算叶尖速比，再根据最佳叶尖速比得到对应的最优转距。目前风速信息的获取通常有两类方法：一类使用风速计或者使用激光雷达测风装置。由于风力机组实际运行时处于风流场环境中，风速大小方向在整个风机叶轮平面上分布不均，而且受塔影、风湍流、气温等的影响。使用风速计只能测量叶轮后某一点的风速，其无法代表有效风速，同时受到叶轮运动扰动后的测量误差较大与整个风力机旋转平面所受到的有效风速有较大差别。使用激光雷达测风装置则价格昂贵，加之使用激光雷达存在风流场演化计算等不确定因素，不利于大规模应用。另一类利用风力发电机系统的物理特性与风能利用系数Cp，设计观测器进行风速观测。这类方式的主要缺陷在于使用了风能利用系数Cp，由于风能利用系数Cp在部分区存在多值解导致求解耗时部分区域准确性不佳等问题，限制了实际可作用风速区间。因此，基于以上方式进行风力发电机组控制均不是理想的可靠方法。
技术实现思路
为了解决上述技术问题，本专利技术提供一种算法简单、计算精度高的风力发电机组叶尖速比获取方法。本专利技术解决上述问题的技术方案是：一种风力发电机组叶尖速比获取方法，包括以下步骤：步骤一：构造传动链观测器，通过使用卡尔曼观测器，从最终观测结果中获得原传动链系统的真实状态；步骤二：通过传感器获得叶轮轮毂转速ωr以及风力发电机组的发电机转矩Tg，得到叶轮实际吸收的气动转矩Ta；步骤三：建立风力发电机的气动转矩Ta与风力发电机的叶尖速比λ的映射关系；步骤四：通过二次曲面的散点拟合多项式函数，得到叶尖速比λ。上述风力发电机组叶尖速比获取方法，风机等效传动链的动力特性由如下方程组表示：其中，Je表示等效转动惯量，表示风机轮毂加速度，Ta是叶轮实际吸收的气动转矩，Tg表示风机发电机转矩，Tloss表示风机转矩损失，Pg表示风机发电机功率,ωr表示风机轮毂转速。上述风力发电机组叶尖速比获取方法，所述步骤一具体步骤为：假定风力发电机传动链线性时变系统的观测表达式为：dx(t)＝(A(t)x(t)+B(t)u(t))dt+dn(t)y(tk)＝H(tk)x(tk)+D(tk)u(tk)+v(tk)式中，dx(t)为系统状态在t时刻的一阶微分，A(t)为风机系统在t时刻的状态矩阵，B(t)为风机系统在t时刻的输入矩阵，x(t)为系统在t时刻的状态量，u(t)为在t时刻的系统输入，dt为微分时间，dn(t)为在t时刻的噪声干扰项；y(tk)为在tk时刻的系统输出，H(tk)为风机系统在tk时刻的卡尔曼观测输出矩阵，D(tk)为风机系统在tk时刻的卡尔曼观测前馈矩阵，x(tk)为系统在tk时刻的状态量，u(tk)为在tk时刻的系统输入，v(tk)为在tk时刻的观测噪声干扰项；给定其初始测量值以及实际测量数据y(t0),y(t1),y(t2),…,u(t),t≥t0上式中，y(t0),y(t1),y(t2)为系统在t0、t1、t2时刻的输出，表示系统初始状态的估计，表示在时刻的系统状态估计，表示在时刻的卡尔曼矩阵，P0表示初始状态的卡尔曼矩阵；对于实际风机系统，使用tk时刻实际传感器值更新状态更新观测结果，得到观测结果中获得原传动链系统的真实状态，其中包含叶轮轮毂转速观测值以及叶轮加速度观测值上述风力发电机组叶尖速比获取方法，所述步骤二中，使用步骤一中卡尔曼观测器并考虑工程偏差进一步得到如下表达式：其中k表示第k时刻，Je为轮毂总等效转矩，为观测得到的叶轮加速度，为叶轮实际吸收的气动转矩观测值，Tloss为系统效率损失转矩。上述风力发电机组叶尖速比获取方法，所述步骤三具体步骤为：引入风力发电机组转矩系数Cq定义如下：其中，Cq(λ,θ)是风力发电机的转矩系数，ρ是空气密度，ve是等效叶轮风速，S是叶轮等效扫风面积，R是叶轮等效半径，叶尖速比λ用来表示叶轮在不同风速下的状态，对于风力发电机的叶尖速比λ有如下定义：ve＝ωrRλ综合转矩系数Cq方程并去除其中风速相关参数有:风力发电机的转矩系数Cq(λ,θ)<1是表示在外界风流场作用下风力发电机的获得旋转扭矩的系数，其是关于叶尖速比λ和桨距角θ的一个非线性曲面函数；从上述表达式分离气动转矩Ta得到：式中，表示在条件下求解λ；明确Ta＝Ta(λ,θ,ωr)是关于叶尖速比λ、轮毂转速ωr和桨距角θ的函数；在风力发电机实际运行过程中，由于轮毂转速ωr、桨距角θ可在每个采样周期测量得到，因此当已知Cq(λ,θ)曲面的情况下，使用轮毂转速ωr、桨距角θ以及步骤三所述方式获得气动转矩Ta计算得到叶尖速比λ。上述风力发电机组叶尖速比获取方法，所述步骤四具体步骤为：对于实际运行的风力发电机，Cq(λ,θ)曲面是已知的，为使用数学方法表达Cq(λ,θ)曲面，工程使用以下表达式获得其分段拟合：式中，m为将桨距角θ离散化后的某一离散点，n为多项式曲线拟合阶数，α为各项系数，θm为第m个θ离散点，αm,n表示在第m个θ离散点的第n项λ系数；当桨距角θ为不同的确定取值时，风能利用系数Cq关于叶尖速比λ构成一簇曲线，根据不同θ值下的Cq(λ)分段曲线，拟合以叶尖速比λ为自变量的多项式，得到气动转矩Ta关于有效风速λ的表达式，最后通过迭代计算找到λ的最佳估计；考虑步骤三估计得到的风机气动转矩观测值和步骤四得到的气动转矩Ta关于叶尖速比λ的映射关系Ta＝Ta(λ,θ,ωr)，有：式中，k表示第k时刻，通过步骤二得到的气动转矩观测值和通过步骤一得到的叶轮轮毂转速观测值都是已知值，对于特定风机系统的叶轮半径R与空气密度ρ也是已知的，因此当前时刻叶尖速比λ观测值是唯一变量。本专利技术的有益效果在于：本专利技术首先构造传动链观测器，通过使用卡尔曼观测器，从最终观测结果中获得原传动链系统的真实状态；然后通过传感器获得叶轮轮毂转速ωr以及风力发电机组的发电机转矩Tg，得到叶轮实际吸收的气动转矩Ta；再建立风力发电机的气动转矩Ta与风力发电机的叶尖速比λ的映射关系；最后通过二次曲面的散点拟合多项式函数，得到叶尖速比λ，本专利技术通过对风力发电机物理系统本质特性的分析本专利技术直接减少了风速获取和风速估计部分，可在风力发电机控制系统在线实时运行环境下下，获得风电场风力发电机组的有效叶尖速比，有助于提高风力发电机组控制精度提升风机本文档来自技高网...

【技术保护点】
1.一种风力发电机组叶尖速比获取方法，包括以下步骤：步骤一：构造传动链观测器，通过使用卡尔曼观测器，从最终观测结果中获得原传动链系统的真实状态；步骤二：通过传感器获得叶轮轮毂转速ωr以及风力发电机组的发电机转矩Tg，得到叶轮实际吸收的气动转矩Ta；步骤三：建立风力发电机的气动转矩Ta与风力发电机的叶尖速比λ的映射关系；步骤四：通过二次曲面的散点拟合多项式函数，得到叶尖速比λ。

【技术特征摘要】
1.一种风力发电机组叶尖速比获取方法，包括以下步骤：步骤一：构造传动链观测器，通过使用卡尔曼观测器，从最终观测结果中获得原传动链系统的真实状态；步骤二：通过传感器获得叶轮轮毂转速ωr以及风力发电机组的发电机转矩Tg，得到叶轮实际吸收的气动转矩Ta；步骤三：建立风力发电机的气动转矩Ta与风力发电机的叶尖速比λ的映射关系；步骤四：通过二次曲面的散点拟合多项式函数，得到叶尖速比λ。2.根据权利要求1所述的风力发电机组叶尖速比获取方法，其特征在于，所述步骤一中，风机等效传动链的动力特性由如下方程组表示：其中，Je表示等效转动惯量，表示风机轮毂加速度，Ta是叶轮实际吸收的气动转矩，Tg表示风机发电机转矩，Tloss表示风机转矩损失，Pg表示风机发电机功率,ωr表示风机轮毂转速。3.根据权利要求2所述的风力发电机组叶尖速比获取方法，其特征在于，所述步骤一具体步骤为：假定风力发电机传动链线性时变系统的观测表达式为：dx(t)＝(A(t)x(t)+B(t)u(t))dt+dn(t)y(tk)＝H(tk)x(tk)+D(tk)u(tk)+v(tk)式中，dx(t)为系统状态在t时刻的一阶微分，A(t)为风机系统在t时刻的状态矩阵，B(t)为风机系统在t时刻的输入矩阵，x(t)为系统在t时刻的状态量，u(t)为在t时刻的系统输入，dt为微分时间，dn(t)为在t时刻的噪声干扰项；y(tk)为在tk时刻的系统输出，H(tk)为风机系统在tk时刻的卡尔曼观测输出矩阵，D(tk)为风机系统在tk时刻的卡尔曼观测前馈矩阵，x(tk)为系统在tk时刻的状态量，u(tk)为在tk时刻的系统输入，v(tk)为在tk时刻的观测噪声干扰项；、给定其初始测量值以及实际测量数据y(t0),y(t1),y(t2),…,u(t),t≥t0上式中，y(t0),y(t1),y(t2)为系统在t0、t1、t2时刻的输出，表示系统初始状态的估计，表示在时刻的系统状态估计，表示在时刻的卡尔曼矩阵，P0表示初始状态的卡尔曼矩阵；对于实际风机系统，使用tk时刻实际传感器值更新状态更新观测结果，得到观测结果中获得原传动链系统的真实状态，其中包含叶轮轮毂转速观测值以及叶轮加速度观测值4.根据权利要求3所述的风力发电机组叶尖速比获取方法，其特征在于，所述步骤二中，使用步骤...

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