一种风电机组塔体倾斜度计算方法技术

本发明专利技术公开了一种风电机组塔体倾斜度计算方法，建立塔体倾斜角度同风速、桨距角之间的关系式，采用叠加法计算得到的塔体偏移量；当风速v≤12m/s时，采用自适应最小二乘算法拟合风速与塔体倾斜度之间的函数关系，进而求得当风速为零或接近于零时，风机塔体的倾斜度；大于额定风速运行情况下算式的求取；当12<v≤25m/s时，采用多元非线性回归分析法求取倾斜角度同风速、桨距角之间的关系，建立二次多项式回归模型，采用最小二乘法估计多元线性回归模型的参数，得到倾斜角度。本发明专利技术的有益效果是能得到不同风速情况下风力发电机组塔体倾斜的角度，计算结果准确。

【技术实现步骤摘要】

本专利技术属于风电
，涉及。
技术介绍
风电机组塔体在风力发电机中主要起支撑作用，塔体本身承受自身的重力、风的 推力、叶轮的扭力等复杂多变的负荷，同时受气象及地质因素的影响，使得风力发电机组在 运行的过程中，塔体作为一个弹性刚体会产生一定幅度的摇摆。在长期运行过程中，塔基会 因塔体摇摆等因素的作用而产生沉降等现象，进而使塔体发生偏斜。塔体过大的倾斜会影 响风力发电机组的正常运行，严重的还会产生安全事故。因此，研究出一种风力发电机组塔 体倾斜度计算方法是非常有必要的。 在现有的技术方案中，有一种方案通过两个倾角传感器计算得到风电机组塔体的 倾斜度：此方案包括在塔体顶部设置的第一倾角传感器、塔体底座设置的第二倾角传感器 以及分别与两个传感器相连接的处理器。处理器根据第一倾角数据和第二倾角数据分别计 算得到塔体顶部中心相对于塔体底座中心在水平方向上的第一偏移量和塔体底座中心在 水平方向上的第二偏移量，进而经过处理器处理得到第一倾角数据ax、ay和第二倾角数 据Px、Py，当其大于报警阈值时，便向主控机房发出报警信号。 同时，也有研究者对风电机组运行时塔筒的受力情况进行了分析。通过分析塔筒 所受推力、升力力矩、重力对其产生的影响，进一步得到塔筒倾斜的偏移量与风速之间的关 系。 上述采用传感器测量计算方案只是动态测量计算塔体顶部中心同塔体底部中心 的相对偏移量，得不到静止条件下塔体的实际偏移量。在风电机组塔体偏移量与受力之间 关系的研究中，其只是得到一个偏移量与风速之间的关系，并没有考虑不同风速下其受力 情况的变化以及倾斜度与桨距角之间的关系。【
技术实现思路
】 本专利技术的目的在于提供，解决了目前在风电机 组塔体偏移量与受力之间关系的研究中，其只是得到一个偏移量与风速之间的关系，并没 有考虑不同风速下其受力情况的变化以及倾斜度与桨距角之间的关系的问题。 本专利技术所采用的技术方案是按照以下步骤进行： 步骤1 :塔体倾斜角度同风速、桨距角之间的关系式如下： 式中，v代表风速，e代表桨距角，0代表塔体倾斜角度，f表示拟合函数关系； 在风电机组运行过程中，其受风的推力、自身的的重力等产生倾斜，采用叠加法计 算得到的塔体偏移量f?为：式中，h为塔筒高，E为弹性模量，I为塔筒弯曲时关于中性轴的惯性矩，P为空气 密度，S为叶片面积，I、1(2为常数，v为风速； 步骤2 :小于额定风速运行情况下算式的求取； 当v< 12m/s时，采用自适应最小二乘算法拟合风速与塔体倾斜度之间的函数关 系，进而求得当风速为零或接近于零时，风机塔体的倾斜度； 步骤3 :大于额定风速运行情况下算式的求取； 当12〈v彡25m/s时，采用多元非线性回归分析法求取倾斜角度同风速、桨距角之 间的关系，建立二次多项式回归模型 9 =k1+k2v+k3v2+k4 0 +k5 0 2+k6v0 (3) 式中，v为风速，0为桨距角，kji= 1，…，6)为常系数，表示回归参数；采用最小 二乘法估计多元线性回归模型的参数，得到倾斜角度。 进一步，所述步骤2的方法为：1)对得到的数据组（Vl，0J(i= 1，2,…m)，Vl、0i表示一组数据中第i个风速 和塔体倾斜角度，选取n= 2次多项式进行拟合，拟合函数： 用常规的最小二乘法求出拟合函数的系数= (U…即得到第一次拟合函 数 3)数据（Vl，0i)(i=1，2,…m)经过一次拟合迭代重整= 对2)对拟合曲线进行上下平移A0，A0的取值为： 卜,6^1' = 1义，)再进行以上步骤，得到第二次拟合重整数据^42^'=1爲___爾)，计算这 两次迭代拟合的RMS均方根值之差X，若x>e，e为误差，则继续迭代，直到满足x彡e， 迭代停止，此时得到的拟合函数…即为n= 2时的最终拟合结果，拟合次数为 k式I M2； 经过该拟合得到其拟合函数关系，通过取极限把C:(v)或令v = 〇,最终得到风 机塔体在静止条件下倾斜的角度。 进一步，所述步骤3中估计多元线性回归模型的参数方法为： 设序列Xl，x2，…，xn，进行如下变换：(6) 式中，则又。y2，…yj新序列，其均值为0,方差 为1 ;令Zi=V，Z2=V2,z3= 0，z4= 0 2,z5=V0，将上述非线性模型转化为多元线性模 型，即为： 9 =k1+k2z1+k3z2+k4z3+k5Z4+k6Z5 (7) 已知n组监测数据（Zll，zl2，…，zl5; 0J，样本求得的误差为： 采用最小二乘法估计多元线性回归模型的参数，即令下式取最小值时的解。 式中，Q为观测值与估计值之差的平方和，0 1表示塔体倾斜角度的第i个观测值， |表示塔体倾斜角度的第i个估计值；对式（9)分别求关于回归参数Lk2, ....k6的偏导 数，并令其等于〇 ;然后联立求解即得到回归参数k2, . . . .k6)的估计值，估计值表示为 H??i,，；最后，将估计值代入式（5)得倾斜角度。 本专利技术的有益效果是提供一种风力发电机组倾斜角与风速、桨距角之间的函数关 系，并计算得到静止时风力发电机组塔体倾斜的角度，对风电机组塔体进行状态监测；然后 在倾斜角度变化超过预设值时发出预警信号，防止发生风电机组倒塌事故，能够精确得到 不同风速下的倾斜度。【附图说明】 图1是风力发电机组功率曲线示意图。【具体实施方式】 下面结合【具体实施方式】对本专利技术进行详细说明。 1风电机组运行工况的划分 变速变桨距风电机组的一个重要运行特性就是运行工况随风速变化切换的特性， 如图1所示。根据风速情况和风力发电机的功率特性，可以将整个运行过程划分为四个典 型工况，每个工况下变桨距系统控制的目标与策略均有所不同。这四个典型工况分别为：起 动并网阶段、最大风能捕获阶段、恒功率控制阶段、超风速切出阶段。 由图1可以看出，在第一工况下，风速未达到切入风速，风机处于停运状态，此时 塔体倾斜只与风速有关；在第二工况下，风机处于欠功率运行状态，此时桨距角维持不变， 当风机转速达到额定转速时，通过调节叶尖速比A使风机输出的机械功率继续增大至额 定功率，这时整个机组在最佳的状态运行；当风机处于第三工况时，为防止其输出功率超过 额定值，需要通过变桨距系统调整叶片的桨距角，降低风能利用率CP，从而使风机能够以额 定功率长期稳定输出，此时叶片受到的阻力增大；当风速大于切出风速时，风机停机。 2风电机组塔体倾斜度分段计算算式的求取 从上述分析中可知，当风机处于第一、第二工况时，塔体的倾斜角度只与风速有 关；当风机处于第三工况时，塔体的倾斜角度不仅与风速有关还与桨距角有关。考虑到第四 工况出现频率较小，本专利技术不进行深入分析。设风机切入风速v= 3m/s，额定风速v= 12m/ s，切出风速v= 25m/s，从而可得，塔体倾斜角度同风速、桨距角之间的关系式如下： 式中，v代表风速，0代表桨距角，0代表塔体倾斜角度，f表示拟合函数关系。 在风电机组运行过程中，其受风的推力、自身的的重力等产生倾斜，采用叠加法计 算得到的塔体偏移量为：[当前第1页1 2 本文档来自技高网...

【技术保护点】
一种风电机组塔体倾斜度计算方法，其特征在于按照以下步骤进行：步骤1：塔体倾斜角度同风速、桨距角之间的关系式如下：θ=f(v)v≤12m/sf(v,β)12<v≤25m/s---(1)]]>式中，v代表风速，β代表桨距角，θ代表塔体倾斜角度，f表示拟合函数关系；在风电机组运行过程中，其受风的推力、自身的的重力等产生倾斜，采用叠加法计算得到的塔体偏移量f为：f=h2EI·9K12ρ216+K22S2h2·v2---(2)]]>式中，h为塔筒高，E为弹性模量，I为塔筒弯曲时关于中性轴的惯性矩，ρ为空气密度，S为叶片面积，K1、K2为常数，v为风速；步骤2：小于额定风速运行情况下算式的求取；当v≤12m/s时，采用自适应最小二乘算法拟合风速与塔体倾斜度之间的函数关系，进而求得当风速为零或接近于零时，风机塔体的倾斜度；步骤3：大于额定风速运行情况下算式的求取；当12<v≤25m/s时，采用多元非线性回归分析法求取倾斜角度同风速、桨距角之间的关系，建立二次多项式回归模型θ＝k1+k2v+k3v2+k4β+k5β2+k6vβ   (3)式中，v为风速，β为桨距角，ki(i＝1,…,6)为常系数，表示回归参数；采用最小二乘法估计多元线性回归模型的参数，得到倾斜角度。...

【技术特征摘要】

【专利技术属性】
技术研发人员：赵洪山，徐樊浩，邓嵩，徐文岐，
申请(专利权)人：华北电力大学保定，
类型：发明
国别省市：河北;13