CN105041584B - 一种风电机组塔体倾斜度计算方法 - Google Patents
一种风电机组塔体倾斜度计算方法 Download PDFInfo
- Publication number
- CN105041584B CN105041584B CN201510299649.2A CN201510299649A CN105041584B CN 105041584 B CN105041584 B CN 105041584B CN 201510299649 A CN201510299649 A CN 201510299649A CN 105041584 B CN105041584 B CN 105041584B
- Authority
- CN
- China
- Prior art keywords
- wind speed
- tower body
- inclination
- tower
- wind
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/728—Onshore wind turbines
Landscapes
- Wind Motors (AREA)
Abstract
本发明公开了一种风电机组塔体倾斜度计算方法，建立塔体倾斜角度同风速、桨距角之间的关系式，采用叠加法计算得到的塔体偏移量；当风速v≤12m/s时，采用自适应最小二乘算法拟合风速与塔体倾斜度之间的函数关系，进而求得当风速为零或接近于零时，风机塔体的倾斜度；大于额定风速运行情况下算式的求取；当12<v≤25m/s时，采用多元非线性回归分析法求取倾斜角度同风速、桨距角之间的关系，建立二次多项式回归模型，采用最小二乘法估计多元线性回归模型的参数，得到倾斜角度。本发明的有益效果是能得到不同风速情况下风力发电机组塔体倾斜的角度，计算结果准确。
Description
技术领域
本发明属于风电技术领域，涉及一种风电机组塔体倾斜度计算方法。
背景技术
风电机组塔体在风力发电机中主要起支撑作用，塔体本身承受自身的重力、风的推力、叶轮的扭力等复杂多变的负荷，同时受气象及地质因素的影响，使得风力发电机组在运行的过程中，塔体作为一个弹性刚体会产生一定幅度的摇摆。在长期运行过程中，塔基会因塔体摇摆等因素的作用而产生沉降等现象，进而使塔体发生偏斜。塔体过大的倾斜会影响风力发电机组的正常运行，严重的还会产生安全事故。因此，研究出一种风力发电机组塔体倾斜度计算方法是非常有必要的。
在现有的技术方案中，有一种方案通过两个倾角传感器计算得到风电机组塔体的倾斜度：此方案包括在塔体顶部设置的第一倾角传感器、塔体底座设置的第二倾角传感器以及分别与两个传感器相连接的处理器。处理器根据第一倾角数据和第二倾角数据分别计算得到塔体顶部中心相对于塔体底座中心在水平方向上的第一偏移量和塔体底座中心在水平方向上的第二偏移量，进而经过处理器处理得到第一倾角数据αx、αy和第二倾角数据βx、βy，当其大于报警阈值时，便向主控机房发出报警信号。
同时，也有研究者对风电机组运行时塔筒的受力情况进行了分析。通过分析塔筒所受推力、升力力矩、重力对其产生的影响，进一步得到塔筒倾斜的偏移量与风速之间的关系。
上述采用传感器测量计算方案只是动态测量计算塔体顶部中心同塔体底部中心的相对偏移量，得不到静止条件下塔体的实际偏移量。在风电机组塔体偏移量与受力之间关系的研究中，其只是得到一个偏移量与风速之间的关系，并没有考虑不同风速下其受力情况的变化以及倾斜度与桨距角之间的关系。
发明内容
本发明的目的在于提供一种风电机组塔体倾斜度计算方法，解决了目前在风电机组塔体偏移量与受力之间关系的研究中，其只是得到一个偏移量与风速之间的关系，并没有考虑不同风速下其受力情况的变化以及倾斜度与桨距角之间的关系的问题。
本发明所采用的技术方案是按照以下步骤进行：
步骤1：塔体倾斜角度同风速、桨距角之间的关系式如下：
式中，v代表风速，β代表桨距角，θ代表塔体倾斜角度，f表示拟合函数关系；
在风电机组运行过程中，其受风的推力、自身的的重力等产生倾斜，采用叠加法计算得到的塔体偏移量f为：
式中，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)为常系数，表示回归参数；采用最小二乘法估计多元线性回归模型的参数，得到倾斜角度。
进一步，所述步骤2的方法为：
1)对得到的数据组(vi,θi)(i＝1,2,…m)，vi、θi表示一组数据中第i个风速和塔体倾斜角度，选取n＝2次多项式进行拟合，拟合函数：
用常规的最小二乘法求出拟合函数的系数即得到第一次拟合函数：
2)对拟合曲线进行上下平移△θ，△θ的取值为：
对的值，取
对的值，取
3)数据(vi,θi)(i＝1,2,…m)经过一次拟合迭代重整对再进行以上步骤，得到第二次拟合重整数据计算这两次迭代拟合的RMS均方根值之差x，若x>ε，ε为误差，则继续迭代，直到满足x≤ε，迭代停止，此时得到的拟合函数即为n＝2时的最终拟合结果，拟合次数为M2；
经过该拟合得到其拟合函数关系，通过取极限或令v＝0，最终得到风机塔体在静止条件下倾斜的角度。
进一步，所述步骤3中估计多元线性回归模型的参数方法为：
设序列x1,x2,…,xn，进行如下变换：
式中， 则y1,y2,…yn为新序列，其均值为0，方差为1；令z1＝v，z2＝v2，z3＝β，z4＝β2，z5＝vβ，将上述非线性模型转化为多元线性模型，即为：
θ＝k1+k2z1+k3z2+k4z3+k5z4+k6z5 (7)
已知n组监测数据(zi1,zi2,…,zi5；θi)，样本求得的误差为：
采用最小二乘法估计多元线性回归模型的参数，即令下式取最小值时的解。
式中，Q为观测值与估计值之差的平方和，θi表示塔体倾斜角度的第i个观测值，表示塔体倾斜角度的第i个估计值；对式(9)分别求关于回归参数k1,k2,....k6的偏导数，并令其等于0；然后联立求解即得到回归参数(k1,k2,....k6)的估计值，估计值表示为最后，将估计值代入式(5)得倾斜角度。
本发明的有益效果是提供一种风力发电机组倾斜角与风速、桨距角之间的函数关系，并计算得到静止时风力发电机组塔体倾斜的角度，对风电机组塔体进行状态监测；然后在倾斜角度变化超过预设值时发出预警信号，防止发生风电机组倒塌事故，能够精确得到不同风速下的倾斜度。
附图说明
图1是风力发电机组功率曲线示意图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
1风电机组运行工况的划分
变速变桨距风电机组的一个重要运行特性就是运行工况随风速变化切换的特性，如图1所示。根据风速情况和风力发电机的功率特性，可以将整个运行过程划分为四个典型工况，每个工况下变桨距系统控制的目标与策略均有所不同。这四个典型工况分别为：起动并网阶段、最大风能捕获阶段、恒功率控制阶段、超风速切出阶段。
由图1可以看出，在第一工况下，风速未达到切入风速，风机处于停运状态，此时塔体倾斜只与风速有关；在第二工况下，风机处于欠功率运行状态，此时桨距角维持不变，当风机转速达到额定转速时，通过调节叶尖速比λ使风机输出的机械功率继续增大至额定功率，这时整个机组在最佳的状态运行；当风机处于第三工况时，为防止其输出功率超过额定值，需要通过变桨距系统调整叶片的桨距角，降低风能利用率CP，从而使风机能够以额定功率长期稳定输出，此时叶片受到的阻力增大；当风速大于切出风速时，风机停机。
2风电机组塔体倾斜度分段计算算式的求取
从上述分析中可知，当风机处于第一、第二工况时，塔体的倾斜角度只与风速有关；当风机处于第三工况时，塔体的倾斜角度不仅与风速有关还与桨距角有关。考虑到第四工况出现频率较小，本发明不进行深入分析。设风机切入风速v＝3m/s，额定风速v＝12m/s，切出风速v＝25m/s，从而可得，塔体倾斜角度同风速、桨距角之间的关系式如下：
式中，v代表风速，β代表桨距角，θ代表塔体倾斜角度，f表示拟合函数关系。
在风电机组运行过程中，其受风的推力、自身的的重力等产生倾斜，采用叠加法计算得到的塔体偏移量为：
式中，h为塔筒高，E为弹性模量，I为塔筒弯曲时关于中性轴的惯性矩，ρ为空气密度，S为叶片面积，K1、K2为常数，v为风速。
·表示相乘。
对于一确定的风机，空气密度变化不大时，式中为一常数取为K3，则：f＝K3v2，即其与风速的平方成正比。因此，在建立多元非线性回归模型时只分析到平方项。
2.1小于额定风速运行情况下算式的求取
当v≤12m/s时，采用自适应最小二乘算法拟合风速与塔体倾斜度之间的函数关系，进而求得当风速为零或接近于零时，风机塔体的倾斜度。具体步骤如下：
对得到的数据组(vi,θi)(i＝1,2,…m)(vi、θi表示一组数据中第i个风速和塔体倾斜角度)，选取n＝2次多项式进行拟合，拟合函数：
用常规的最小二乘法求出拟合函数的系数即得到第一次拟合函数：
对拟合曲线进行上下平移△θ，△θ的取值为：
对的值，取
对的值，取
数据(vi,θi)(i＝1,2,…m)经过一次拟合迭代重整对再进行以上步骤，得到第二次拟合重整数据计算这两次迭代拟合的RMS(均方根值)之差x。若x>ε(ε为误差，如取值1×10-6)，则继续迭代，直到满足x≤ε，迭代停止，此时得到的拟合函数即为n＝2时的最终拟合结果，拟合次数为M2；
经过该拟合得到其拟合函数关系，通过取极限或令v＝0，最终得到风机塔体在静止条件下倾斜的角度。
2.2大于额定风速运行情况下算式的求取
当12<v≤25m/s时，此时风机塔体倾斜角度不仅与风速有关，而且与桨距角有关。采用多元非线性回归分析法求取倾斜角度同风速、桨距角之间的关系，建立二次多项式回归模型为：
θ＝k1+k2v+k3v2+k4β+k5β2+k6vβ (5)
式中，v为风速，β为桨距角，ki(i＝1,…,6)为常系数。
为了消除量纲影响和变量自身变异大小和数值大小的影响，利用标准化后的数据进行数据分析。此处对倾斜角θ、风速v、桨距角β采用Z-score标准化法。具体方法为：设序列x1,x2,…,xn，进行如下变换：
式中， 则y1,y2,…yn为新序列，其均值为0，方差为1。
令z1＝v，z2＝v2，z3＝β，z4＝β2，z5＝vβ，这样可将上述非线性模型转化为多元线性模型，即为：
θ＝k1+k2z1+k3z2+k4z3+k5z4+k6z5 (7)
已知n组监测数据(zi1,zi2,…,zi5；θi)，样本求得的误差为：
采用最小二乘法估计多元线性回归模型的参数，即令下式取最小值时的解。
式中，Q为观测值与估计值之差的平方和，θi表示塔体倾斜角度的第i个观测值，表示塔体倾斜角度的第i个估计值。
对式(9)分别求关于回归参数k1,k2,....k6的偏导数，并令其等于0；然后联立求解即可得到回归参数(k1,k2,....k6)的估计值，估计值表示为最后，将估计值代入式(5)便可得倾斜角度同风速、桨距角之间的函数关系式。
下面列举具体实施例对本发明进行说明。
实施例1：按照发明的公式，求得不同风速下的倾斜角度，如表1所示。
表1倾斜角、风速及桨距角对应数据
最终得到：
当v＝0时，θ＝0.1352°，即：静止时风机塔体偏斜角度为0.1352°。
以上所述仅是对本发明的较佳实施方式而已，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施方式所做的任何简单修改，等同变化与修饰，均属于本发明技术方案的范围内。
Claims (1)
1.一种风电机组塔体倾斜度计算方法，其特征在于按照以下步骤进行：
步骤1：塔体倾斜度同风速、桨距角之间的关系式如下：
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mo>&le;</mo>
<mn>12</mn>
<mi>m</mi>
<mo>/</mo>
<mi>s</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>12</mn>
<mo><</mo>
<mi>v</mi>
<mo>&le;</mo>
<mn>25</mn>
<mi>m</mi>
<mo>/</mo>
<mi>s</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，v代表风速，β代表桨距角，θ代表塔体倾斜度，f表示拟合函数关系；
在风电机组运行过程中，其受风的推力、自身的重力产生倾斜，采用叠加法计算得到的塔体偏移量f为：
<mrow>
<mi>f</mi>
<mo>=</mo>
<mfrac>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mrow>
<mi>E</mi>
<mi>I</mi>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msqrt>
<mrow>
<mfrac>
<mrow>
<mn>9</mn>
<msubsup>
<mi>K</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msup>
<mi>&rho;</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>16</mn>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>K</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<msup>
<mi>S</mi>
<mn>2</mn>
</msup>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>&CenterDot;</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，h为塔筒高，E为弹性模量，I为塔筒弯曲时关于中性轴的惯性矩，ρ为空气密度，S为叶片面积，K1、K2为常数，v为风速；
步骤2：小于额定风速运行情况下算式的求取；
当v≤12m/s时，采用自适应最小二乘算法拟合风速与塔体倾斜度之间的函数关系，进而求得当风速为零或接近于零时，风机塔体的倾斜度；
对得到的数据组(vi,θi)(i＝1,2,…m)，vi、θi表示一组数据中第i个风速和塔体倾斜度，选取n＝2次多项式进行拟合，拟合函数：
<mrow>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<msup>
<mi>v</mi>
<mi>k</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
用常规的最小二乘法求出拟合函数的系数即得到第一次拟合函数：
对拟合曲线进行上下平移Δθ，Δθ的取值为：
<mrow>
<mi>&Delta;</mi>
<mi>&theta;</mi>
<mo>=</mo>
<mi>M</mi>
<mi>A</mi>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&theta;</mi>
<mi>n</mi>
<mn>1</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2.0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
对的值，取
对的值，取
数据(vi,θi)(i＝1,2,…m)经过一次拟合迭代重整对再进行以上步骤，得到第二次拟合重整数据计算这两次迭代拟合的RMS均方根值之差x，若x>ε，ε为误差，则继续迭代，直到满足x≤ε，迭代停止，此时得到的拟合函数即为n＝2时的最终拟合结果，拟合次数为M2，通过取极限或令v＝0，最终得到风机塔体在静止条件下倾斜度；
步骤3：大于额定风速运行情况下算式的求取；
当12＜v≤25m/s时，采用多元非线性回归分析法求取倾斜度同风速、桨距角之间的关系，建立二次多项式回归模型
θ＝k1+k2v+k3v2+k4β+k5β2+k6vβ (5)
式中，v为风速，β为桨距角，ki(i＝1,…,6)为常系数，表示回归参数；采用最小二乘法估计多元线性回归模型的参数，得到倾斜度；
估计多元线性回归模型的参数方法为：
设序列x1,x2,…,xn，进行如下变换：
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mi>s</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，则y1,y2,…yn为新序列，其均值为0，方差为1；令z1＝v，z2＝v2，z3＝β，z4＝β2，z5＝vβ，将非线性模型转化为多元线性模型，即为：
θ＝k1+k2z1+k3z2+k4z3+k5z4+k6z5 (7)
已知n组监测数据(zi1,zi2,…,zi5；θi)，样本求得的误差为：
<mrow>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mo>...</mo>
<msub>
<mi>k</mi>
<mn>6</mn>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
采用最小二乘法估计多元线性回归模型的参数，即令下式取最小值时的解；
<mrow>
<mi>Q</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mo>...</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>6</mn>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Q为观测值与估计值之差的平方和，θi表示塔体倾斜度的第i个观测值，表示塔体倾斜度的第i个估计值；对式(9)分别求关于回归参数k1,k2,….k6的偏导数，并令Q等于0；然后联立求解即得到回归参数(k1,k2,….k6)的估计值，估计值表示为最后，将估计值代入式(5)得倾斜度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510299649.2A CN105041584B (zh) | 2015-06-03 | 2015-06-03 | 一种风电机组塔体倾斜度计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510299649.2A CN105041584B (zh) | 2015-06-03 | 2015-06-03 | 一种风电机组塔体倾斜度计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105041584A CN105041584A (zh) | 2015-11-11 |
| CN105041584B true CN105041584B (zh) | 2018-01-30 |
Family
ID=54448462
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510299649.2A Expired - Fee Related CN105041584B (zh) | 2015-06-03 | 2015-06-03 | 一种风电机组塔体倾斜度计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105041584B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ES2875430T3 (es) | 2016-10-07 | 2021-11-10 | Siemens Gamesa Renewable Energy As | Determinación del ángulo de inclinación de una torre de turbina eólica |
| CN113390376B (zh) * | 2020-03-13 | 2023-04-07 | 新疆金风科技股份有限公司 | 确定风力发电机组的机舱位移的方法、装置和系统 |
| CN111475951B (zh) * | 2020-04-09 | 2023-06-02 | 国网山东省电力公司电力科学研究院 | 一种热电机组工况分析方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2063110A1 (en) * | 2007-11-26 | 2009-05-27 | Siemens Aktiengesellschaft | Method of damping power vibrations of a wind turbine and inclination control system |
| CN102929217A (zh) * | 2012-09-04 | 2013-02-13 | 上海发电设备成套设计研究院 | 一种单台发电机组可靠性在线监控装置及方法 |
| CN104655133A (zh) * | 2015-02-13 | 2015-05-27 | 西安应用光学研究所 | 一种车载升降桅杆平台姿态动态测量方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070092888A1 (en) * | 2003-09-23 | 2007-04-26 | Cornelius Diamond | Diagnostic markers of hypertension and methods of use thereof |
| US20050069936A1 (en) * | 2003-09-26 | 2005-03-31 | Cornelius Diamond | Diagnostic markers of depression treatment and methods of use thereof |
-
2015
- 2015-06-03 CN CN201510299649.2A patent/CN105041584B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2063110A1 (en) * | 2007-11-26 | 2009-05-27 | Siemens Aktiengesellschaft | Method of damping power vibrations of a wind turbine and inclination control system |
| CN102929217A (zh) * | 2012-09-04 | 2013-02-13 | 上海发电设备成套设计研究院 | 一种单台发电机组可靠性在线监控装置及方法 |
| CN104655133A (zh) * | 2015-02-13 | 2015-05-27 | 西安应用光学研究所 | 一种车载升降桅杆平台姿态动态测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105041584A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111801493B (zh) | 确定用于风力涡轮机的控制设置 | |
| EP3056726B1 (en) | System and method for operating a wind turbine based on rotor blade margin | |
| US10605228B2 (en) | Method for controlling operation of a wind turbine | |
| JP5984791B2 (ja) | 風力発電装置のモニタリングシステム及びモニタリング方法 | |
| CN107709766B (zh) | 校准风力涡轮机的负载传感器的方法 | |
| DK2306003T3 (en) | System and methods for controlling a wind turbine | |
| KR101476986B1 (ko) | 풍력 터빈용 제어 장치 | |
| US10669986B2 (en) | Relating to the determination of rotor imbalances in a wind turbine | |
| CA2840441C (en) | Method and apparatus for wind turbine noise reduction | |
| CA2755154A1 (en) | Method and system for adjusting a power parameter of a wind turbine | |
| Kragh et al. | Rotor speed dependent yaw control of wind turbines based on empirical data | |
| EP2981710B1 (en) | Method for controlling thrust load on a wind turbine | |
| EP2876300A1 (en) | Methods and systems to shut down a wind turbine | |
| CN108757312A (zh) | 一种风力发电机变桨控制方法 | |
| KR101383792B1 (ko) | 라이다 측정을 이용한 나셀 풍속계 보정 방법 | |
| CN105041584B (zh) | 一种风电机组塔体倾斜度计算方法 | |
| US6993965B2 (en) | Horizontal axis wind turbine and method for measuring upflow angle | |
| CN107100795A (zh) | 一种低风速下风力发电机组mppt自适应控制方法 | |
| JP2022107523A (ja) | 風の乱流のアクティブセンシングを用いた風力タービンのための推力制御 | |
| WO2019150805A1 (ja) | 風力発電装置及び風力発電システム | |
| WO2016128003A1 (en) | Control system capable of estimating a spatial wind field of a wind turbine system having multiple rotors | |
| Pitance et al. | Experimental validation of Pharwen code using data from Vertical-axis wind turbines | |
| US10961981B2 (en) | Control for a wind turbine | |
| Ossmann et al. | Multivariable controller design verification for a Liberty wind turbine | |
| Qu et al. | Improved Data-Driven Yaw Misalignment Calibration of Wind Turbine via LiDAR Verification |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180130 Termination date: 20190603 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |