CN113323804B - 解决风力发电机组塔架二阶前后振动的控制方法与模块 - Google Patents
解决风力发电机组塔架二阶前后振动的控制方法与模块 Download PDFInfo
- Publication number
- CN113323804B CN113323804B CN202110653039.3A CN202110653039A CN113323804B CN 113323804 B CN113323804 B CN 113323804B CN 202110653039 A CN202110653039 A CN 202110653039A CN 113323804 B CN113323804 B CN 113323804B
- Authority
- CN
- China
- Prior art keywords
- blade
- pitch
- angle
- tower
- azimuth
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000001133 acceleration Effects 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 23
- 230000001965 increasing effect Effects 0.000 claims abstract description 9
- 238000001914 filtration Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000011217 control strategy Methods 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims 1
- 238000005452 bending Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/022—Adjusting aerodynamic properties of the blades
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D1/00—Wind motors with rotation axis substantially parallel to the air flow entering the rotor
- F03D1/06—Rotors
- F03D1/065—Rotors characterised by their construction elements
- F03D1/0675—Rotors characterised by their construction elements of the blades
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D13/00—Assembly, mounting or commissioning of wind motors; Arrangements specially adapted for transporting wind motor components
- F03D13/20—Arrangements for mounting or supporting wind motors; Masts or towers for wind motors
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/0204—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/022—Adjusting aerodynamic properties of the blades
- F03D7/0224—Adjusting blade pitch
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/0296—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor to prevent, counteract or reduce noise emissions
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/328—Blade pitch angle
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/329—Azimuth or yaw angle
-
- 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
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种解决风力发电机组塔架二阶前后振动的控制方法与模块,基于塔架中间段的前后加速度以及叶轮方位角,分析二阶模态阵型,采用P控制器,在常规变桨控制器输出的变桨速率指令上,给三个叶片分别叠加一个不同的动态变桨速率,经过限幅得到各个叶片最终的变桨速率值,然后根据计算得出各个叶片最终桨距角给定值,即通过微调桨距角,根据叶片不同方位角来增加或减小桨距角,减小叶片受力不平衡,最终达到减小塔架在二阶频率点上的前后振动的目的,保障风力发电机组安全稳定运行。
Description
技术领域
本发明涉及风力发电机组的技术领域,尤其是指一种解决风力发电机组塔架二阶前后振动的控制方法与模块。
背景技术
随着风力发电技术的发展以及市场降低度电成本需求,风力发电机组容量越来越大,叶片越来越长,塔架越来越高,并且由于海上资源丰富,风速稳定,东南沿海一带经济发达,电量需求量大,可自行解决消纳问题,提高风电能利用率,因此海上风电将迎来快速发展。
但是无论是陆上机组还是海上机组,塔架振动问题都是一个棘手的问题。特别是海上风力发电机组的基础立于海床上,塔架上端受到风速、尾流等的影响,下端经受洋流海浪等的影响,以及塔架固有频率的降低,其塔架振动问题相较于陆上机组更严重,尤其是塔架二阶前后振动更是严重。
对于塔架前后振动问题,现采用塔架前后加阻策略,即在常规变桨控制器的基础上,增加一个基于塔架前后加速度的变桨给定值,通过微调桨距角增加塔架前后振动模态阻尼,从而减小塔架前后振动。此方案针对一阶塔架前后振动有较好效果,但针对以塔架前后二阶模态为主的振动,振动效果并不明显。
发明内容
本发明的第一目的在于克服现有技术的缺点与不足,提供了一种解决风力发电机组塔架二阶前后振动的控制方法,可减小塔架前后二阶模态为主的振动,从而解决塔架在二阶频率点上的前后振动,保障风力发电机组安全稳定运行。
本发明的第二目的在于提供一种解决风力发电机组塔架二阶前后振动的控制模块。
本发明的第一目的通过下述技术方案实现:一种解决风力发电机组塔架二阶前后振动的控制方法,该方法是基于塔架中间段的前后加速度以及叶轮方位角的独立变桨控制策略IPC(Individual Pitch Control)来解决塔架在二阶频率点的前后振动,通过减小叶轮旋转平面内三个叶片面内受力不均而造成的气动不平衡,从根本上减小塔架前后的振动,即根据叶片不同方位角来增加或减小桨距角,从而减小三个叶片受力不平衡,进而减小塔架前后二阶模态为主的振动,具体过程是:首先,通过传感器实时测量塔架中间段的前后加速度,并对其进行滤波;然后,再测量并计算三个叶片对应的方位角;根据塔架中间段的前后加速度以及方位角,采用P控制器,得到各个叶片额外的变桨速率,并将其叠加到各个叶片常规变桨速率上,即给三个叶片分别叠加一个不同的动态变桨速率,最后,经过限幅并计算出各个叶片最终的桨距角给定值输送给机组的变桨***执行,从而解决塔架在二阶频率点的前后振动问题。
进一步,所述解决风力发电机组塔架二阶前后振动的控制方法,包括以下步骤:
1)测量塔架中间段的前后加速度anacelle,并进行滤波,经过传递函数F(s)滤波后的塔架中间段的前后加速度为anacellef,即:
anacellef=anacelle*F(s)
式中,s为拉普拉斯算子,ξ为阻尼系数,ω为频率,T1和T2为时间常数;其中,定义anacellef的正方向为向前时为正;
2)计算叶轮方位角;
测量叶片1方位角,根据三个叶片均匀分布在一个叶轮平面,计算出叶片2 和叶片3的方位角,定义叶片垂直向上方位角为0°,通过传感器测量得到叶片1方位角为则叶片2的方位角是在叶片1的方位角基础上加叶片3的方位角是在叶片1的方位角基础上加即:
式中,RotorAzimuth1为叶片1的方位角,RotorAzimuth2为叶片2的方位角,RotorAzimuth3为叶片3的方位角;
3)计算各个叶片额外的变桨速率;
根据塔架中间段的前后加速度,基于叶片在不同方位角下的受力分析,分别计算三个叶片分别所需的额外变桨速率,具体如下:
式中,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率,Kp为P控制器的比例增益系数,A为叶轮方位角超前角度;
4)计算各个叶片对应变桨速率;
将步骤3)得到的各个叶片的额外变桨速率值与常规控制器输出的变桨速率值Vo(k)进行叠加,经过限幅,分别得到三个叶片对应的变桨速率,具体如下:
叶片1的变桨速率V1(k):V1(k)=Vo(k)+ΔV1(k)
叶片2的变桨速率V2(k):V2(k)=Vo(k)+ΔV2(k)
叶片3的变桨速率V3(k):V3(k)=Vo(k)+ΔV3(k)
再将各个叶片的变桨速率进行限幅,得到叶片1、叶片2、叶片3最终的变桨速率分别为Vf1(k)、Vf2(k)、Vf3(k),其中Vo(k)为常规控制器计算得到的统一变桨速率;
5)计算各个叶片对应的桨距角给定值;
根据步骤4)计算得到最终的变桨速率,分别计算三个叶片对应的桨距角给定值,并将其输送到变桨***执行,实现三个叶片桨距角的单独变桨控制,即 IPC控制;
叶片1的桨距角:
θ1(k)=θ1(k-1)+Vf1(k)*T
叶片2的桨距角:
θ2(k)=θ2(k-1)+Vf2(k)*T
叶片3的桨距角:
θ3(k)=θ3(k-1)+Vf3(k)*T
式中,θ1(k-1)为前一时刻叶片1的桨距角,θ1(k)为叶片1当前时刻的桨距角,θ2(k-1)为前一时刻叶片2的桨距角,θ2(k)为叶片2当前时刻的桨距角,θ3(k-1)为前一时刻叶片3的桨距角,θ3(k)为叶片3当前时刻的桨距角,T为Controller Cycle time控制算法循环时间常数;
6)执行桨距角变桨;
变桨***根据三个叶片对应的桨距角给定值调节叶片桨距角,减小由于气动不平衡引起的塔架在二阶频率点上的前后振动,从而实现风力发电机组基于塔架中间段的前后加速度以及叶轮方位角的振动优化控制。
本发明的第二目的通过下述技术方案实现:一种解决风力发电机组塔架二阶前后振动的控制模块,该模块是基于塔架中间段的前后加速度以及叶轮方位角的独立变桨控制策略IPC(Individual Pitch Control)来解决塔架在二阶频率点的前后振动,通过减小叶轮旋转平面内三个叶片面内受力不均而造成的气动不平衡,从根本上减小塔架前后的振动,即根据叶片不同方位角来增加或减小桨距角,从而减小三个叶片受力不平衡,进而减小塔架前后二阶模态为主的振动,具体过程是:首先,通过传感器实时测量塔架中间段的前后加速度,并对其进行滤波;然后,再测量并计算三个叶片对应的方位角;根据塔架中间段的前后加速度以及方位角,采用P控制器,得到各个叶片额外的变桨速率,并将其叠加到各个叶片常规变桨速率上,即给三个叶片分别叠加一个不同的动态变桨速率,最后,经过限幅并计算出各个叶片最终的桨距角给定值输送给机组的变桨***执行,从而解决塔架在二阶频率点的前后振动问题。
进一步,所述解决风力发电机组塔架二阶前后振动的控制模块,包括:
测量单元,用于测量塔架中间段的前后加速度,并进行滤波;
叶轮方位角计算单元,用于测量计算各叶片方位角;
叶片额外变桨速率计算单元,用于计算各个叶片额外的变桨速率,根据塔架中间段的前后加速度,基于叶片在不同方位角下的受力分析,分别计算三个叶片分别所需的额外变桨速率;
叶片对应变桨速率计算单元,用于计算各个叶片对应的变桨速率,将叶片额外变桨速率计算单元计算得到的各个叶片额外的变桨速率值与常规控制器输出的变桨速率值进行叠加,经过限幅,分别得到三个叶片对应的变桨速率;
叶片对应桨距角给定值计算单元,用于计算三个叶片对应的桨距角给定值,根据叶片对应变桨速率计算单元计算得到最终的变桨速率,分别计算三个叶片对应的桨距角给定值,并将其输送给变桨***执行调节叶片桨距角,实现三个叶片桨距角的单独变桨控制,即IPC控制。
进一步,在所述测量单元中,测量塔架中间段的前后加速度anacelle,并进行滤波,经过传递函数F(s)滤波后的塔架中间段的前后加速度为anacellef,即:
anacellef=anacelle*F(s)
式中,s为拉普拉斯算子,ξ为阻尼系数,ω为频率,T1和T2为时间常数;其中,定义anacellef的正方向为向前时为正。
进一步,在所述叶轮方位角计算单元中,测量叶片1方位角,根据三个叶片均匀分布在一个叶轮平面,计算出叶片2和叶片3的方位角,定义叶片垂直向上方位角为0°,通过传感器测量得到叶片1方位角为则叶片2的方位角是在叶片1的方位角基础上加叶片3的方位角是在叶片1的方位角基础上加即:
式中,RotorAzimuth1为叶片1的方位角,RotorAzimuth2为叶片2的方位角,RotorAzimuth3为叶片3的方位角。
进一步,在所述叶片额外变桨速率计算单元中,采用以下公式:
式中,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率,Kp为P控制器的比例增益系数,anacellef为滤波后的塔架中间段的前后加速度,A为叶轮方位角超前角度。
进一步,在所述叶片对应变桨速率计算单元中,采用以下公式:
叶片1的变桨速率V1(k):V1(k)=Vo(k)+ΔV1(k)
叶片2的变桨速率V2(k):V2(k)=Vo(k)+ΔV2(k)
叶片3的变桨速率V3(k):V3(k)=Vo(k)+ΔV3(k)
再将各个叶片的变桨速率进行限幅,得到叶片1、叶片2、叶片3最终的变桨速率分别为Vf1(k)、Vf2(k)、Vf3(k),其中Vo(k)为常规控制器计算得到的统一变桨速率,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率。
进一步,在所述叶片对应桨距角给定值计算单元中,采用以下公式:
叶片1的桨距角:
θ1(k)=θ1(k-1)+Vf1(k)*T
叶片2的桨距角:
θ2(k)=θ2(k-1)+Vf2(k)*T
叶片3的桨距角:
θ3(k)=θ3(k-1)+Vf3(k)*T
式中,θ1(k-1)为前一时刻叶片1的桨距角,θ1(k)为叶片1当前时刻的桨距角,θ2(k-1)为前一时刻叶片2的桨距角,θ2(k)为叶片2当前时刻的桨距角,θ3(k-1)为前一时刻叶片3的桨距角,θ3(k)为叶片3当前时刻的桨距角,T为 Controller Cycle time控制算法循环时间常数,Vf1(k)、Vf2(k)、Vf3(k)别为叶片1、叶片2、叶片3最终的变桨速率。
本发明与现有技术相比,具有如下优点与有益效果:
1、本发明采取基于塔架中间段的前后加速度,在常规变桨控制器的基础上叠加变桨速率的控制方式减小塔架前后振动,使得机组稳定运行。相对于现有的桨距角控制方式,能有效减小机组振动(特别是塔架在二阶频率点上的前后振动)而引起的载荷偏大问题。
2、本发明无需增加机组设备,只需在控制方法中增加相应的功能模块,就能减小塔架前后振动,减小机组载荷,从而节省成本,提高机组竞争力。
3、本发明有很强的理论依据,易于被相关技术人员接受,对后续控制优化改进和设备维护奠定基础。
综上所述,本发明可以有效可减小塔架前后二阶模态为主的振动,从而解决塔架在二阶频率点上的前后振动,降低机组载荷,保障风力发电机组安全稳定运行,本方案同时适用于一切由于气动不平衡造成载荷偏大的工况,具有十分广泛的应用前景,值得推广。
附图说明
图1为二阶模态阵型图。
图2为本发明方法的流程图。
图3为解决风力发电机组塔架二阶前后振动的控制框图。
图4为本发明模块的架构图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
实施例1
首先我们分析了风力发电机组的受力情况:
叶轮的叶片由于受到重力、风剪切、湍流、偏航误差以及上升气流的影响,三个叶片受力不均衡导致塔顶受到波动的合力与合弯矩,引起塔架的振动,特别是不平衡的前后方向倾覆弯矩会激化塔架二阶模态振动,二阶模态对塔架中间段的振动影响大,二阶模态阵型如图1所示,其表现为塔架中间段加速度大,两端小,而基于塔架中间段前后振动加速度大小,结合基于方位角的IPC控制策略,产生作用于塔架前后方向的倾覆弯矩,抵消由于湍流、风切变等产生的不平衡倾覆弯矩,进而达到降低塔架二阶模态振动的目的。
具体控制如下:当叶片在0°方位角(叶片位于垂直正上方)时,测量到塔架中间段有向前的加速度时,意味着塔顶受到正向的倾覆弯矩,此时应调大叶片的桨矩角,降低该处叶片的气动推力,进而降低作用于塔顶的倾覆弯矩;同理,当叶片在180°方位角(叶片位于垂直正下方)时,此时应调小叶片的桨矩角,增大该处叶片的气动推力,产生作用于塔顶负向的倾覆弯矩,同样起到削弱正向倾覆弯矩,降低振动的目的。如若塔架中间段有向后的加速度时,分析结论则刚好相反。
因此,基于上述分析,本实施例提供了一种解决风力发电机组塔架二阶前后振动的控制方法,其是基于塔架中间段的前后加速度以及叶轮方位角的独立变桨控制策略IPC(Individual Pitch Control)来解决塔架在二阶频率点的前后振动,通过减小叶轮旋转平面内三个叶片面内受力不均而造成的气动不平衡,从根本上减小塔架前后的振动,即根据叶片不同方位角来增加或减小桨距角,从而减小三个叶片受力不平衡,进而减小塔架前后二阶模态为主的振动,如图2 所示,具体过程是:首先,通过传感器实时测量塔架中间段的前后加速度,并对其进行滤波;然后,再测量并计算三个叶片对应的方位角;根据塔架中间段的前后加速度以及方位角,采用P控制器,得到各个叶片额外的变桨速率,并将其叠加到各个叶片常规变桨速率上,即给三个叶片分别叠加一个不同的动态变桨速率,最后,经过限幅并计算出各个叶片最终的桨距角给定值输送给机组的变桨***执行,从而解决塔架在二阶频率点的前后振动问题;其具体包括以下步骤:
1)测量塔架中间段的前后加速度anacelle,并进行滤波,经过传递函数F(s)滤波后的塔架中间段的前后加速度为anacellef,即:
anacellef=anacelle*F(s)
式中,s为拉普拉斯算子,ξ为阻尼系数,ω为频率,T1和T2为时间常数;其中,定义anacellef的正方向为向前时为正;
2)计算叶轮方位角;
测量叶片1方位角,根据三个叶片均匀分布在一个叶轮平面,计算出叶片2 和叶片3的方位角,定义叶片垂直向上方位角为0°,通过传感器测量得到叶片 1方位角为则叶片2的方位角是在叶片1的方位角基础上加叶片3的方位角是在叶片1的方位角基础上加即:
式中,RotorAzimuth1为叶片1的方位角,RotorAzimuth2为叶片2的方位角,RotorAzimuth3为叶片3的方位角;
3)计算各个叶片额外的变桨速率;
根据塔架中间段的前后加速度,基于叶片在不同方位角下的受力分析,分别计算三个叶片分别所需的额外变桨速率,具体如下:
式中,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率,Kp为P控制器的比例增益系数,A为叶轮方位角超前角度;
4)计算各个叶片对应变桨速率;
将步骤3)得到的各个叶片的额外变桨速率值与常规控制器输出的变桨速率值Vo(k)进行叠加,经过限幅,分别得到三个叶片对应的变桨速率,具体如下:
叶片1的变桨速率V1(k):V1(k)=Vo(k)+ΔV1(k)
叶片2的变桨速率V2(k):V2(k)=Vo(k)+ΔV2(k)
叶片3的变桨速率V3(k):V3(k)=Vo(k)+ΔV3(k)
再将各个叶片的变桨速率进行限幅,得到叶片1、叶片2、叶片3最终的变桨速率分别为Vf1(k)、Vf2(k)、Vf3(k),其中Vo(k)为常规控制器计算得到的统一变桨速率;
5)计算各个叶片对应的桨距角给定值;
根据步骤4)计算得到最终的变桨速率,分别计算三个叶片对应的桨距角给定值,并将其输送到变桨***执行,实现三个叶片桨距角的单独变桨控制,即IPC控制;
叶片1的桨距角:
θ1(k)=θ1(k-1)+Vf1(k)*T
叶片2的桨距角:
θ2(k)=θ2(k-1)+Vf2(k)*T
叶片3的桨距角:
θ3(k)=θ3(k-1)+Vf3(k)*T
式中,θ1(k-1)为前一时刻叶片1的桨距角,θ1(k)为叶片1当前时刻的桨距角,θ2(k-1)为前一时刻叶片2的桨距角,θ2(k)为叶片2当前时刻的桨距角,θ3(k-1)为前一时刻叶片3的桨距角,θ3(k)为叶片3当前时刻的桨距角,T为Controller Cycle time控制算法循环时间常数;
6)执行桨距角变桨;
变桨***根据三个叶片对应的桨距角给定值调节叶片桨距角,减小由于气动不平衡引起的塔架在二阶频率点上的前后振动,从而实现风力发电机组基于塔架中间段的前后加速度以及叶轮方位角的振动优化控制。
实施例2
见图3所示,风力发电机组常规控制器由转矩控制器和变桨控制器两部分组成:转矩控制器在额定风速以下时起作用,用于风能最大捕获;变桨控制器用于额定风速以上时,通过调节桨距角保持发电机转速在额定转速附近,保证机组正常运行的同时确保机组输出功率为额定功率。考虑到目前长叶片,高塔架,特别是海上风力发电机组的塔架上端受风、尾流等影响,下端经受洋流海浪等的影响,以及塔架固有频率的降低,其塔架振动问题相较于陆上机组更严重,尤其是塔架二阶前后振动更是严重。因此,在常规变桨控制器的基础上,建立了一个新的控制模块,用于解决风力发电机组塔架二阶前后振动问题,该控制模块是基于塔架中间段的前后加速度以及叶轮方位角的独立变桨控制策略 IPC(IndividualPitch Control)来解决塔架在二阶频率点的前后振动,通过减小叶轮旋转平面内三个叶片面内受力不均而造成的气动不平衡,从根本上减小塔架前后的振动,即根据叶片不同方位角来增加或减小桨距角,从而减小三个叶片受力不平衡,进而减小塔架前后二阶模态为主的振动,具体过程是:首先,通过传感器实时测量塔架中间段的前后加速度,并对其进行滤波;然后,再测量并计算三个叶片对应的方位角;根据塔架中间段的前后加速度以及方位角,采用P控制器,得到各个叶片额外的变桨速率,并将其叠加到各个叶片常规变桨速率上,即给三个叶片分别叠加一个不同的动态变桨速率,最后,经过限幅并计算出各个叶片最终的桨距角给定值输送给机组的变桨***(包含变桨控制器、变桨执行器等)执行,减小塔架前后二阶模态为主的振动,从而解决塔架在二阶频率点上的前后振动,降低机组载荷,保障风力发电机组安全稳定运行。
如图4所示,本实施例所提供的解决风力发电机组塔架二阶前后振动的控制模块,包括以下功能单元:
测量单元,用于测量塔架中间段的前后加速度anacelle,并进行滤波,经过传递函数F(s)滤波后的塔架中间段的前后加速度为anacellef,即:
anacellef=anacelle*F(s)
式中,s为拉普拉斯算子,ξ为阻尼系数,ω为频率,T1和T2为时间常数;其中,定义anacellef的正方向为向前时为正。
叶轮方位角计算单元,用于测量计算各叶片方位角,具体情况如下:
测量叶片1方位角,根据三个叶片均匀分布在一个叶轮平面,计算出叶片2 和叶片3的方位角,定义叶片垂直向上方位角为0°,通过传感器测量得到叶片 1方位角为则叶片2的方位角是在叶片1的方位角基础上加叶片3的方位角是在叶片1的方位角基础上加即:
式中,RotorAzimuth1为叶片1的方位角,RotorAzimuth2为叶片2的方位角,RotorAzimuth3为叶片3的方位角。
叶片额外变桨速率计算单元,用于计算各个叶片额外的变桨速率,根据塔架中间段的前后加速度,基于叶片在不同方位角下的受力分析,分别计算三个叶片分别所需的额外变桨速率,具体如下:
式中,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率,Kp为P控制器的比例增益系数,A为叶轮方位角超前角度。
叶片对应变桨速率计算单元,用于计算各个叶片对应的变桨速率,将叶片额外变桨速率计算单元计算得到的各个叶片额外的变桨速率值与常规控制器输出的变桨速率值Vo(k)进行叠加,经过限幅,分别得到三个叶片对应的变桨速率,具体如下:
叶片1的变桨速率V1(k):V1(k)=Vo(k)+ΔV1(k)
叶片2的变桨速率V2(k):V2(k)=Vo(k)+ΔV2(k)
叶片3的变桨速率V3(k):V3(k)=Vo(k)+ΔV3(k)
再将各个叶片的变桨速率进行限幅,得到叶片1、叶片2、叶片3最终的变桨速率分别为Vf1(k)、Vf2(k)、Vf3(k),其中Vo(k)为常规控制器计算得到的统一变桨速率。
叶片对应桨距角给定值计算单元,用于计算三个叶片对应的桨距角给定值,根据叶片对应变桨速率计算单元计算得到最终的变桨速率,分别计算三个叶片对应的桨距角给定值,并将其输送给变桨***执行调节叶片桨距角,实现三个叶片桨距角的单独变桨控制,即IPC控制;其中,采用以下公式计算:
叶片1的桨距角:
θ1(k)=θ1(k-1)+Vf1(k)*T
叶片2的桨距角:
θ2(k)=θ2(k-1)+Vf2(k)*T
叶片3的桨距角:
θ3(k)=θ3(k-1)+Vf3(k)*T
式中,θ1(k-1)为前一时刻叶片1的桨距角,θ1(k)为叶片1当前时刻的桨距角,θ2(k-1)为前一时刻叶片2的桨距角,θ2(k)为叶片2当前时刻的桨距角,θ3(k-1)为前一时刻叶片3的桨距角,θ3(k)为叶片3当前时刻的桨距角,T为 Controller Cycle time控制算法循环时间常数。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (2)
1.一种解决风力发电机组塔架二阶前后振动的控制方法,其特征在于,该方法是基于塔架中间段的前后加速度以及叶轮方位角的独立变桨控制策略IPC来解决塔架在二阶频率点的前后振动,通过减小叶轮旋转平面内三个叶片面内受力不均而造成的气动不平衡,从根本上减小塔架前后的振动,即根据叶片不同方位角来增加或减小桨距角,从而减小三个叶片受力不平衡,进而减小塔架前后二阶模态为主的振动,具体过程是:首先,通过传感器实时测量塔架中间段的前后加速度,并对其进行滤波;然后,再测量并计算三个叶片对应的方位角;根据塔架中间段的前后加速度以及方位角,采用P控制器,得到各个叶片额外的变桨速率,并将其叠加到各个叶片常规变桨速率上,即给三个叶片分别叠加一个不同的动态变桨速率,最后,经过限幅并计算出各个叶片最终的桨距角给定值输送给机组的变桨***执行,从而解决塔架在二阶频率点的前后振动问题;其包括以下步骤:
1)测量塔架中间段的前后加速度anacelle,并进行滤波,经过传递函数F(s)滤波后的塔架中间段的前后加速度为anacellef,即:
anacellef=anacelle*F(s)
式中,s为拉普拉斯算子,ξ为阻尼系数,ω为频率,T1和T2为时间常数;其中,定义anacellef的正方向为向前时为正;
2)计算叶轮方位角;
测量叶片1方位角,根据三个叶片均匀分布在一个叶轮平面,计算出叶片2和叶片3的方位角,定义叶片垂直向上方位角为0°,通过传感器测量得到叶片1方位角为则叶片2的方位角是在叶片1的方位角基础上加叶片3的方位角是在叶片1的方位角基础上加即:
式中,RotorAzimuth1为叶片1的方位角,RotorAzimuth2为叶片2的方位角,RotorAzimuth3为叶片3的方位角;
3)计算各个叶片额外的变桨速率;
根据塔架中间段的前后加速度,基于叶片在不同方位角下的受力分析,分别计算三个叶片分别所需的额外变桨速率,具体如下:
式中,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率,Kp为P控制器的比例增益系数,A为叶轮方位角超前角度;
4)计算各个叶片对应变桨速率;
将步骤3)得到的各个叶片的额外变桨速率值与常规控制器输出的变桨速率值Vo(k)进行叠加,经过限幅,分别得到三个叶片对应的变桨速率,具体如下:
叶片1的变桨速率V1(k):V1(k)=Vo(k)+ΔV1(k)
叶片2的变桨速率V2(k):V2(k)=Vo(k)+ΔV2(k)
叶片3的变桨速率V3(k):V3(k)=Vo(k)+ΔV3(k)
再将各个叶片的变桨速率进行限幅,得到叶片1、叶片2、叶片3最终的变桨速率分别为Vf1(k)、Vf2(k)、Vf3(k),其中Vo(k)为常规控制器计算得到的统一变桨速率;
5)计算各个叶片对应的桨距角给定值;
根据步骤4)计算得到最终的变桨速率,分别计算三个叶片对应的桨距角给定值,并将其输送到变桨***执行,实现三个叶片桨距角的单独变桨控制,即IPC控制;
叶片1的桨距角:
θ1(k)=θ1(k-1)+Vf1(k)*T
叶片2的桨距角:
θ2(k)=θ2(k-1)+Vf2(k)*T
叶片3的桨距角:
θ3(k)=θ3(k-1)+Vf3(k)*T
式中,θ1(k-1)为前一时刻叶片1的桨距角,θ1(k)为叶片1当前时刻的桨距角,θ2(k-1)为前一时刻叶片2的桨距角,θ2(k)为叶片2当前时刻的桨距角,θ3(k-1)为前一时刻叶片3的桨距角,θ3(k)为叶片3当前时刻的桨距角,T为Controller Cycle time控制算法循环时间常数;
6)执行桨距角变桨;
变桨***根据三个叶片对应的桨距角给定值调节叶片桨距角,减小由于气动不平衡引起的塔架在二阶频率点上的前后振动,从而实现风力发电机组基于塔架中间段的前后加速度以及叶轮方位角的振动优化控制。
2.一种解决风力发电机组塔架二阶前后振动的控制模块,其特征在于,该模块是基于塔架中间段的前后加速度以及叶轮方位角的独立变桨控制策略IPC来解决塔架在二阶频率点的前后振动,通过减小叶轮旋转平面内三个叶片面内受力不均而造成的气动不平衡,从根本上减小塔架前后的振动,即根据叶片不同方位角来增加或减小桨距角,从而减小三个叶片受力不平衡,进而减小塔架前后二阶模态为主的振动,具体过程是:首先,通过传感器实时测量塔架中间段的前后加速度,并对其进行滤波;然后,再测量并计算三个叶片对应的方位角;根据塔架中间段的前后加速度以及方位角,采用P控制器,得到各个叶片额外的变桨速率,并将其叠加到各个叶片常规变桨速率上,即给三个叶片分别叠加一个不同的动态变桨速率,最后,经过限幅并计算出各个叶片最终的桨距角给定值输送给机组的变桨***执行,从而解决塔架在二阶频率点的前后振动问题;其包括:
测量单元,用于测量塔架中间段的前后加速度,并进行滤波;
叶轮方位角计算单元,用于测量计算各叶片方位角;
叶片额外变桨速率计算单元,用于计算各个叶片额外的变桨速率,根据塔架中间段的前后加速度,基于叶片在不同方位角下的受力分析,分别计算三个叶片分别所需的额外变桨速率;
叶片对应变桨速率计算单元,用于计算各个叶片对应的变桨速率,将叶片额外变桨速率计算单元计算得到的各个叶片额外的变桨速率值与常规控制器输出的变桨速率值进行叠加,经过限幅,分别得到三个叶片对应的变桨速率;
叶片对应桨距角给定值计算单元,用于计算三个叶片对应的桨距角给定值,根据叶片对应变桨速率计算单元计算得到最终的变桨速率,分别计算三个叶片对应的桨距角给定值,并将其输送给变桨***执行调节叶片桨距角,实现三个叶片桨距角的单独变桨控制,即IPC控制;
在所述测量单元中,测量塔架中间段的前后加速度anacelle,并进行滤波,经过传递函数F(s)滤波后的塔架中间段的前后加速度为anacellef,即:
anacellef=anacelle*F(s)
式中,s为拉普拉斯算子,ξ为阻尼系数,ω为频率,T1和T2为时间常数;其中,定义anacellef的正方向为向前时为正;
在所述叶轮方位角计算单元中,测量叶片1方位角,根据三个叶片均匀分布在一个叶轮平面,计算出叶片2和叶片3的方位角,定义叶片垂直向上方位角为0°,通过传感器测量得到叶片1方位角为则叶片2的方位角是在叶片1的方位角基础上加叶片3的方位角是在叶片1的方位角基础上加即:
式中,RotorAzimuth1为叶片1的方位角,RotorAzimuth2为叶片2的方位角,RotorAzimuth3为叶片3的方位角;
在所述叶片额外变桨速率计算单元中,采用以下公式:
式中,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率,Kp为P控制器的比例增益系数,anacellef为滤波后的塔架中间段的前后加速度,A为叶轮方位角超前角度;
在所述叶片对应变桨速率计算单元中,采用以下公式:
叶片1的变桨速率V1(k):V1(k)=Vo(k)+ΔV1(k)
叶片2的变桨速率V2(k):V2(k)=Vo(k)+ΔV2(k)
叶片3的变桨速率V3(k):V3(k)=Vo(k)+ΔV3(k)
再将各个叶片的变桨速率进行限幅,得到叶片1、叶片2、叶片3最终的变桨速率分别为Vf1(k)、Vf2(k)、Vf3(k),其中Vo(k)为常规控制器计算得到的统一变桨速率,ΔV1(k)、ΔV2(k)、ΔV3(k)分别为叶片1、叶片2、叶片3所需的额外变桨速率;
在所述叶片对应桨距角给定值计算单元中,采用以下公式:
叶片1的桨距角:
θ1(k)=θ1(k-1)+Vf1(k)*T
叶片2的桨距角:
θ2(k)=θ2(k-1)+Vf2(k)*T
叶片3的桨距角:
θ3(k)=θ3(k-1)+Vf3(k)*T
式中,θ1(k-1)为前一时刻叶片1的桨距角,θ1(k)为叶片1当前时刻的桨距角,θ2(k-1)为前一时刻叶片2的桨距角,θ2(k)为叶片2当前时刻的桨距角,θ3(k-1)为前一时刻叶片3的桨距角,θ3(k)为叶片3当前时刻的桨距角,T为Controller Cycle time控制算法循环时间常数,Vf1(k)、Vf2(k)、Vf3(k)别为叶片1、叶片2、叶片3最终的变桨速率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110653039.3A CN113323804B (zh) | 2021-06-11 | 2021-06-11 | 解决风力发电机组塔架二阶前后振动的控制方法与模块 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110653039.3A CN113323804B (zh) | 2021-06-11 | 2021-06-11 | 解决风力发电机组塔架二阶前后振动的控制方法与模块 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113323804A CN113323804A (zh) | 2021-08-31 |
CN113323804B true CN113323804B (zh) | 2022-03-29 |
Family
ID=77420552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110653039.3A Active CN113323804B (zh) | 2021-06-11 | 2021-06-11 | 解决风力发电机组塔架二阶前后振动的控制方法与模块 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113323804B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114109720B (zh) * | 2021-12-31 | 2023-04-25 | 中国华能集团清洁能源技术研究院有限公司 | 一种水平轴风力机风轮增功提效装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104214045A (zh) * | 2013-05-30 | 2014-12-17 | 成都阜特科技股份有限公司 | 双馈式变速变桨风力发电机组的独立变桨距控制方法 |
KR20170054828A (ko) * | 2015-11-10 | 2017-05-18 | 대우조선해양 주식회사 | 풍력발전기 타워의 진동저감 제어시스템 |
CN110439747A (zh) * | 2019-08-02 | 2019-11-12 | 明阳智慧能源集团股份公司 | 一种降低风电塔筒左右方向振动及载荷的ipc控制方法 |
CN112128052A (zh) * | 2020-08-18 | 2020-12-25 | 明阳智慧能源集团股份公司 | 一种基于控制策略优化的风力发电机组塔架降载方法 |
CN112412700A (zh) * | 2019-01-28 | 2021-02-26 | 柳盼 | 一种风电机组塔架振动控制*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9644606B2 (en) * | 2012-06-29 | 2017-05-09 | General Electric Company | Systems and methods to reduce tower oscillations in a wind turbine |
WO2015086023A1 (en) * | 2013-12-09 | 2015-06-18 | Vestas Wind Systems A/S | Counteracting tower oscillations of an idling wind turbine |
-
2021
- 2021-06-11 CN CN202110653039.3A patent/CN113323804B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104214045A (zh) * | 2013-05-30 | 2014-12-17 | 成都阜特科技股份有限公司 | 双馈式变速变桨风力发电机组的独立变桨距控制方法 |
KR20170054828A (ko) * | 2015-11-10 | 2017-05-18 | 대우조선해양 주식회사 | 풍력발전기 타워의 진동저감 제어시스템 |
CN112412700A (zh) * | 2019-01-28 | 2021-02-26 | 柳盼 | 一种风电机组塔架振动控制*** |
CN110439747A (zh) * | 2019-08-02 | 2019-11-12 | 明阳智慧能源集团股份公司 | 一种降低风电塔筒左右方向振动及载荷的ipc控制方法 |
CN112128052A (zh) * | 2020-08-18 | 2020-12-25 | 明阳智慧能源集团股份公司 | 一种基于控制策略优化的风力发电机组塔架降载方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113323804A (zh) | 2021-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Namik et al. | Individual blade pitch control of floating offshore wind turbines | |
EP2422080B1 (en) | Extracting wave energy in a wind turbine installation | |
CN102168650B (zh) | 基于主控的兆瓦级风力机统一和独立变桨混合控制方法 | |
CN104428531B (zh) | 操作风力涡轮机的方法以及适合于所述方法的*** | |
CN105971821B (zh) | 一种风力发电机组基于风轮推力预估的控制方法 | |
CN110446853A (zh) | 用于管理风力涡轮机塔架的扭转振荡的***和方法 | |
EP2556248B1 (en) | A wind turbine | |
US20130134711A1 (en) | Wind turbine | |
CN103527405B (zh) | 一种双馈式风力发电机组变桨距控制方法 | |
CN105179168B (zh) | 一种大型风电机组塔架虚拟阻尼控制方法 | |
CN103742357B (zh) | 一种风力发电机组风轮非对称载荷控制方法 | |
CN106532739A (zh) | 风电机组分频段参与电力***一次调频方法 | |
CN112523941B (zh) | 一种防止风力发电机组超速的控制方法与模块 | |
CN102926930A (zh) | 一种风力发电***的独立变桨控制方法 | |
CN105156270A (zh) | 风力发电机独立变桨控制***及方法 | |
CN103629047A (zh) | 一种降低风电机组载荷的非线性桨距角控制方法 | |
CN112128052A (zh) | 一种基于控制策略优化的风力发电机组塔架降载方法 | |
CN113323804B (zh) | 解决风力发电机组塔架二阶前后振动的控制方法与模块 | |
CN105351144A (zh) | 一种减小风机疲劳载荷的桨叶振动反馈方法 | |
CN113224774B (zh) | 一种电网频率扰动过程风储协调控制方法及*** | |
CN113309665B (zh) | 解决风力发电机组塔架一阶侧向振动的控制方法与模块 | |
CN112610410B (zh) | 一种风力发电机组塔架净空调控方法与模块 | |
US20240159215A1 (en) | Fault-tolerant control method and apparatus of floating wind turbine | |
US11846270B2 (en) | Control method and device of a wind park | |
CN112610411B (zh) | 一种解决风力发电机组塔架净空问题的控制方法与模块 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |