WO2010016764A1 - System and method for compensating rotor imbalance in a wind turbine - Google Patents
System and method for compensating rotor imbalance in a wind turbine Download PDFInfo
- Publication number
- WO2010016764A1 WO2010016764A1 PCT/NL2009/050482 NL2009050482W WO2010016764A1 WO 2010016764 A1 WO2010016764 A1 WO 2010016764A1 NL 2009050482 W NL2009050482 W NL 2009050482W WO 2010016764 A1 WO2010016764 A1 WO 2010016764A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- rotor
- wind turbine
- blade
- imbalance
- rotor blades
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000006243 chemical reaction Methods 0.000 claims abstract description 38
- 238000012544 monitoring process Methods 0.000 claims abstract description 30
- 238000012937 correction Methods 0.000 claims abstract description 25
- 238000005452 bending Methods 0.000 claims description 102
- 230000001133 acceleration Effects 0.000 claims description 39
- 238000005259 measurement Methods 0.000 claims description 39
- 239000013598 vector Substances 0.000 claims description 35
- 238000012545 processing Methods 0.000 claims description 29
- 230000000750 progressive effect Effects 0.000 claims description 22
- 238000001914 filtration Methods 0.000 claims description 18
- 230000015654 memory Effects 0.000 claims description 13
- 230000010354 integration Effects 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 239000011295 pitch Substances 0.000 description 86
- 230000000694 effects Effects 0.000 description 11
- 230000001364 causal effect Effects 0.000 description 7
- 238000013507 mapping Methods 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
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/04—Automatic control; Regulation
- F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
- F03D7/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
- F03D7/045—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with model-based controls
-
- 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/30—Commissioning, e.g. inspection, testing or final adjustment before releasing for production
- F03D13/35—Balancing static or dynamic imbalances
-
- 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/022—Adjusting aerodynamic properties of the blades
- F03D7/024—Adjusting aerodynamic properties of the blades of individual 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
- 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/028—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power
- F03D7/0292—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power to reduce fatigue
-
- 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/10—Purpose of the control system
- F05B2270/109—Purpose of the control system to prolong engine life
- F05B2270/1095—Purpose of the control system to prolong engine life by limiting mechanical stresses
-
- 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/331—Mechanical loads
-
- 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
Definitions
- the present invention relates to a method for compensating rotor imbalance in a wind turbine.
- the invention also relates to a system for compensating rotor imbalance in a wind turbine.
- the invention relates to computer software for implementing said method.
- a wind turbine comprises a rotor provided with a number of rotor blades.
- the rotor is coupled by means of a shaft to an electric generator, possibly via some form of transmission.
- the electric generator is provided with an output for the output of electrical energy, for example through a power converter to an electricity network.
- the assembly of rotor, transmission, generator and power conversion unit are arranged in a nacelle on a mast.
- Wind turbines always show some rotor imbalance due to deviations in the profile and mass properties of the rotor blades, as well as differences in the blade pitch when the rotor is assembled.
- the rotor imbalance results in alternating loads being exerted on the static components. These varying loads are exerted periodically in the azimuth angle of the rotor and are designated as Ip variations that have a period of once every rotation, i.e. they have the same frequency (Ip) as the rotation of the rotor. This therefore leads to variations in the tilting and yaw moments and in the vertical and horizontal forces exerted on the nacelle. In order for the wind turbine to resist these loads, the wind turbine needs to be designed "heavier'.
- variations in tilting Ip and yaw moments often interfere with control circuits, the aim of which is to reduce the blade moments generated by excitations from the flow, for example caused by tower sweeping and by windshear and turbulence, to frequencies in the range of Ip.
- These control circuits reduce blade loads at about Ip by making adjustments to blade pitch variations almost periodically to the usual course of the blade pitch; in a 3-bladed wind turbine these variations are each offset by 120 degrees.
- these blade pitch variations are achieved by Ip modulation of the arithmetically processed tilting and yaw moment. This operation, called a feedback or control law, comprises low-pass filtering, time integration and scaling.
- This object is achieved by applying a method for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor and a mast, wherein the rotor comprises a rotor shaft that is provided with an n number of blades, wherein the rotor shaft is connected to a top section of the mast and wherein a pitch of each rotor blade is individually adjustable by means of a respective actuator; wherein the method comprises the execution in real-time of the following steps repeatedly when in use:
- the determination of at least one correction value for influencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein the feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling;
- the forces exerted on each of the rotor blades can be equalized as much as possible, so that a reduction of the rotor imbalance can be achieved under dynamic operating conditions.
- the invention also enables the exclusion of static balancing or even for it to be performed in a simplified manner, as a result of which the installation of the wind turbine can be carried out with less effort.
- the method according to the invention comprises measuring the bending moments as averages without an absolute offset.
- the method according to the invention comprises the feedback of the quantity indicative of rotor imbalance, comprising the execution of mathematical operations in accordance with a control law according to a method applied in modem control theory, selected from a group comprising Robust Control, Linear Quadratic Regulation (LQR), Linear Quadratic Gaussian control (LQG), H-infinite control (H ⁇ ) and ⁇ -synthesis.
- a control law selected from a group comprising Robust Control, Linear Quadratic Regulation (LQR), Linear Quadratic Gaussian control (LQG), H-infinite control (H ⁇ ) and ⁇ -synthesis.
- the method according to the invention comprises at least one correction value for influencing the aerodynamic conversion that relates to a control signal for adjusting a blade pitch value, and effectuating the at least one specific correction value for the transmission of the control signal to the actuator of the respective rotor blade.
- the method according to the invention comprises determining and monitoring loads, when said wind turbine is in use, comprising measuring bending moment values in a cross-section of the rotor shaft along two non-parallel vectors.
- the method according to the invention includes the wind turbine, further comprising a yaw bearing for pivotally connecting at least the rotor and the mast, wherein the determination of imbalance forces, when said wind turbine is in use, comprises determining a tilt moment and a yaw moment at the yaw bearing.
- the method according to the invention comprises the wind turbine, further comprising a yaw bearing for pivotally connecting at least the rotor and the mast, wherein the determination of imbalance forces, when said wind turbine is in use, comprises determining a tilt moment at the yaw bearing. In one embodiment the method according to the invention comprises measuring a vertical and laterally exerted force in the top of the mast.
- the method according to the invention comprises the measurement of at least a rotational acceleration along a direction vector in a plane defined by a vertical and lateral direction vector. In one embodiment the method according to the invention comprises the measurement of at least a translational acceleration along a direction vector in a plane defined by a vertical and lateral direction vector at some distance from the centre of the top of the mast.
- the method according to the invention includes the wind turbine comprising a nacelle, and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the determination of the tilt moment and the yaw moment at die yaw bearing comprises measuring accelerations at the nacelle in the vertical and horizontal plane at some distance from the centre of the top of the mast.
- the method according to the invention includes the wind turbine comprising a nacelle, and the nacelle comprises the rotor shaft which is coupled to the mast via the yaw bearing and wherein the determination of the tilt moment at the yaw bearing comprises measuring accelerations at the nacelle in the vertical and horizontal plane at some distance from the centre of the top of the mast.
- the method according to the invention includes the wind turbine comprising a nacelle, and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the bending moments are determined as averages without absolute offset, based upon a measured nacelle acceleration perpendicular to the rotor plane and upon measured blade deflection moments that are not necessarily offset-free.
- the method according to the invention comprises that a sensor circuit is selected from a group of at least one sensor and a combination of sensors.
- the invention relates to a computer system for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor and a mast, wherein the rotor comprises a rotor shaft that is provided with an n number of rotor blades, wherein the rotor shaft is coupled to a top section of the mast and wherein a blade pitch of each rotor blade is individually adjustable; wherein the wind turbine is provided with at least one sensor for monitoring forces and blade pitch sensors for monitoring a set blade pitch for each of the rotor blades; wherein the computer is provided with a central processing unit and memory, wherein the memory is linked to the processing unit and wherein the processing unit is linked to the sensors, wherein the computer is arranged, when the wind turbine is in use, to repeat the following operations:
- the determination of at least one correction value for influencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling; - effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
- the invention further relates to computer software stored on a computer-readable medium for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor and a mast, wherein the rotor comprises a rotor shaft that is provided with an number of rotor blades, wherein the rotor shaft is coupled to a top section of the mast and wherein a blade pitch of each rotor blade is individually adjustable; wherein the wind turbine is provided with at least one sensor for monitoring forces and blade pitch sensors for monitoring a set blade pitch for each of the rotor blades; wherein the computer is provided with a central processing unit and memory, wherein the memory is linked to the processing unit and wherein the central processing unit is linked to the sensors, wherein the computer software comprises executable code which, when loaded on the computer, enables the computer to repeatedly execute in real-time the following operations when the wind turbine is in use:
- the determination of at least one correction value for the aerodynamic effect of conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling; - effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1
- the computer software according to the invention comprises determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades comprising the use of either a method of approximation or method of reconstruction.
- the invention further relates to a computer-readable medium that comprises computer executable code which, when loaded on the computer system, enables the computer, as described above, to execute the method as described above.
- figure 1 shows a schematic view of a wind turbine provided with a system in a first embodiment
- figure 2 shows a schematic view of a wind turbine provided with a system in a second embodiment
- figure 3 shows a schematic view of a wind turbine provided with a system in a third embodiment
- figure 4 shows schematic orientations of loads possibly occurring in the wind turbine
- figure 5 shows a flow diagram for a method according to one embodiment
- figure 6 shows a schematic view of a computer for executing computer software in accordance with the method according to one embodiment.
- FIG. 1 shows a schematic view of a wind turbine Ia.
- Wind turbine Ia comprises a rotor R with a rotor shaft provided with a number of rotor blades B.
- the rotor R is coupled to a transmission T for attaching the rotor shaft to an electric generator G by means of the rotor shaft.
- the electric generator G is provided with an output for the output of electrical energy, for example through a power converter (not shown) to an electricity network.
- the assembly of rotor R, transmission T and generator G is located in a nacelle L on a mast M.
- the assembly of rotor R, transmission T and generator G is pivotally connected to the mast M by means of a yaw bearing KS that constitutes part of a yaw system KS.
- the transmission T may be excluded when the rotor shaft is coupled directly to the generator, i.e. in the case of a low-speed generator.
- the wind turbine 1 is adjustable by means of a number of parameters in order to achieve an optimal efficiency under certain arbitrary conditions, wherein the parameters are related to the parts R, B, K, G shown and mentioned here in figure 1. Other parameters could possibly also be of importance, such as would be known to a person skilled in the field.
- the present invention provides the insight that Ip variations in the static components of the wind turbine as a result of rotor imbalance can be eliminated by influencing the aerodynamic conversion (the conversion at the rotor blades of energy from the area of sweep to rotational energy of the rotor) of each individual rotor blade.
- a condition is that the aerodynamic conversion of each rotor blade can be influenced seperately.
- the aerodynamic conversion can be influenced by rotating the rotor blade along a longitudinal axis and by controlling the blade pitch.
- the blade pitch of each rotor blade B is adjustable so that a power absorbed from the flow of air by the wind turbine can be controlled.
- the adjustment of blade pitch by means of an adjustment control system for example by means of actuators (shown schematically for the upper rotor blade by arrow Al) can be implemented in the blade root of each rotor blade, but also halfway along the longitudinal axis.
- any method for influencing the aerodynamic conversion is suitable, as long as the deflection moment in the blade root can be significantly adjusted.
- the deflection moment is defined as the moment of force in a rotor blade, the vector of which is directed perpendicular to the plane in which the rotor shaft and the longitudinal axis of the rotor blade lie.
- a sensor SO on each rotor blade is linked to the blade pitch adjustment system for monitoring a set blade pitch value.
- the yaw system K is arranged pivotally around a yaw bearing KS in order to adjust the position of the rotor R relative to the horizontal direction of the flow of air.
- the power converter is provided with an electrical control system that is arranged to control the electric power and rotational speed characteristic of the generator G.
- each rotor blade is provided with a sensor Sl attached to the blade root in close proximity to the rotor shaft for monitoring a value of the deflection moment in the respective rotor blade. Additionally, in a further embodiment the radial and tangential force in the root of the respective rotor blade can also be measured.
- FIG. 1 shows a schematic view of a wind turbine Ib provided with a system in a second embodiment.
- the rotor shaft of the wind turbine Ib is provided with a sensor S2 for detecting bending moments in the rotor shaft.
- Figure 3 shows a schematic view of a wind turbine Ic provided with a system in a third embodiment.
- the yaw bearing KS of the wind turbine Ic is equipped with a sensor S3 for monitoring the tilting and yaw moment, wherein the tilt moment vector lies on the horizontal plane and is perpendicular to the shaft and the yaw moment vector in the vertical plane parallel to the mast.
- the system can be provided with one or more acceleration sensors S4 for conducting acceleration measurements from which tilting and yaw moments can be deduced.
- An acceleration measurement for example, can be performed by measuring the vertical and lateral nacelle acceleration at a certain distance, in the direction of the rotor R, from the centre of the top of the mast. Both here and in the following description, 'lateral' is designated as a direction in the horizontal plane, perpendicular to the rotor shaft.
- the third embodiment is provided with a position sensor S5 for monitoring the rotation angle of the rotor (rotor pitch), which is the angle between the longitudinal axis of a selected rotor blade and the horizontal as a result of the rotation of the rotor.
- the vertical and lateral force on or close to the yaw bearing can be measured (for example, by sensor S3), or the vertical and lateral acceleration in a second position on the nacelle (for example, by sensor S4).
- a condition in the first embodiment is that the signals recorded from the rotor blades are measured without any individual offset.
- An offset in the measurement results in the minimization of the sum of the actual imbalance and offset. This means that this offset, when compensation is activated, ultimately becomes effective as rotor imbalance.
- the bending moments are measured as averages without applying any absolute offset.
- the amplitudes of the Ip variations are both measured offset-free.
- Figure 4 shows, schematically, orientations in a wind turbine of any loads that may occur.
- the mast M has a vertical orientation Ol .
- the KRM vector of the yaw moment is given in the top of the mast, which is directed along the vertical orientation Ol .
- the vector for the tilt moment TM is directed in the lateral direction (designated by 02, which is perpendicular to the vertical orientation Ol and the rotor shaft RA).
- the rotation vector RV is directed along the rotor shaft RA.
- one rotor blade B is shown. In rotor blade B, the radial force FR (exerted on the rotor centre RC) is directed along the longitudinal axis LA of the rotor blade B.
- FIG. 5 shows a flow diagram for a method 100 according to one embodiment.
- a cyclic algorithm is activated with an indicative cycle time or period time of 0.1s. This then initiates a first step 105 for initializing, among other things, the adjustments to the control signals for the blade pitches.
- step 110 a measurement is conducted of the deflection moments and possibly the radial and tangential forces in the blade roots (first embodiment Ia; use of sensor Sl), or the bending moments in the rotor shaft (second embodiment Ib; use of sensor (S2), or the tilting and yaw moment and the rotor pitch and possibly the vertical and lateral force exerted on the yaw bearing (third embodiment Ic; use of a combination of sensors (S3, S5) or (S4, S5)).
- the bending moment values are determined along two vectors perpendicular to each other in the cross-section of the rotor shaft; one of these vectors is parallel to the longitudinal axis of the rotor blade selected for the rotor pitch.
- These bending moments in embodiment Ia are derived from the measured loads and the rotor geometry, as well as the distance from the rotor centre to the center of the cross-section where the radial and tangential forces also be measured.
- the desired bending moments then result from the measurements and the orientation of the vectors along which measurements are conducted in relation to the vectors of the desired bending moments.
- the bending moments are deduced from the tilting and yaw moment and the value of the rotor pitch, as well as the distance from the yaw bearing to the cross-section if the vertical and lateral force exerted on the yaw bearing is also measured.
- the bending moment values thus obtained are converted to desired adjustments to the blade pitches from the geometry of the two direction vectors for the bending moments in relation to those of the rotor. This is done with the aim of adjusting the aerodynamic load of the individual rotor blades in such a manner that the (progressive) mean of the bending moments is minimized whilst retaining the quantity of energy generated by the wind.
- the invention is based upon the principle that a relatively small change to a blade pitch of a rotor blade results in a substantially proportional change of the load exerted on the rotor blade by the wind.
- the conversion of bending moments to adjustments to the blade pitches comprises the inverse of the trigonometric relationships for deducing shaft bending moments from blade deflection moments, as incorporated in the later presented example of an algorithm.
- the feedback comprises low-pass filtering, time integration and scaling for minimizing the (progressive) mean value of the bending moments via low-frequency adjustment of the blade pitches.
- Minimum mean bending moments in the shaft automatically imply minimum rotor imbalance.
- the blade pitches are adjusted to a time-scale that is greater than approx. three times the cycle time of the rotor.
- step 140 the desired adjustments to the blade pitches derived from step 130 are converted to control signals for the appropriate adjustment of the blade pitches. If a desired adjustment is not possible in connection with the range of a blade pitch, all adjustments are amended with the same, smallest possible offset that does not lie beyond this range. As a result of this offset, the imbalance compensation will have some slight effect on the energy extracted from the wind. It is also possible not to adjust the pitch of one of the blades. Here, the imbalance compensation will generally have a greater effect on the energy extracted from the wind.
- step 150 one first waits till the cycle time has lapsed, after which step 110 is repeated. This is not done after a re-initialization signal 152, after which step 105 is executed, and after a stop signal 154, after which the method is completed by step 160.
- the objective of the imbalance compensation is to minimize the Ip variations in the tilting and yaw moment at the position of the yaw bearing.
- the rotor centre is located at some distance from the yaw bearing KS.
- the rotor centre is understood to be the point of the rotor R where the rotor blades B meet and are coupled to each other.
- the method provides for feedback control for compensation by adjusting the blade pitches based on measurements of load signals in the root of each rotor blade. Based on the measurements of the deflection moments and radial and tangential forces, the bending moment values are determined along two vectors in the (fictive) cross-section of the rotor shaft at the position of the yaw bearing.
- the rotor imbalance results in mean values of the bending moments that are unequal to zero due to the deflection moments via the trigonometric relationships arising from the rotor geometry and the momentum that arises from the radial and tangential forces as a result of not coinciding with rotor centre and yaw bearing.
- Moment is the exerted force multiplied by the arm.
- the arm is defined as the distance between rotor centre and yaw bearing).
- these bending moments are averaged to zero, as described in step 130.
- the Ip variations in the tilt and yaw moment are eliminated at the yaw bearing.
- a third embodiment is also conceivable, wherein the bending moments in the rotor centre are measured directly on the shaft.
- the bending moments at the rotor centre are composed from measurements of the tilt and yaw moment and the lateral and vertical force at the yaw bearing via the rotor pitch and the distance between the yaw bearing and rotor centre.
- a wind turbine is considered with three rotor blades all of which are arranged in a single rotor plane.
- the rotor blades are each offset at 120 degrees within the rotor plane.
- a choice was made to influence the aerodynamic conversion by rotating each of the three rotor blades along the respective longitudinal axis of the rotor blade and to regulate the blade pitch.
- the blade pitches are designated as ⁇ l, ⁇ 2 and ⁇ 3 for the first, second and third rotor blades respectively.
- the bending moment values on the rotor shaft are calculated on the basis of measurements of the deflection moment in the root of each of the rotor blades with the aid of the measuring sensor Sl for measuring the deflection moment of the respective rotor blade.
- a wind load exerted on a rotor blade that is directed perpendicular to the rotor plane causes a deflection moment in the rotor blade.
- the deflection moment is designated as MkI, Mk2 and Mk3 for the first, second and third rotor blade respectively.
- MbI MkI x cos(0°) + Mk2 x cos(120°) + Mk3 x cos(240°) [Ia]
- Mb2 MkI x sin(0°) + Mk2 x sin(120°) + Mk3 x sin(240°) [Ib],
- MbI is the bending moment along a direction vector that coincides with the rotation of the rotor that points in the direction of the longitudinal axis of the first rotor blade
- Mb2 is the bending moment along a direction vector perpendicular to the above-mentioned direction vector of bending moment MbL
- Variations in the blade pitches result in variations in the deflection moments, which are approximately proportional hereto.
- ⁇ 2 ⁇ bl x cos(120°) + ⁇ b2 x sin(120°) [3b]
- ⁇ Mb2 3/2 x C x ⁇ b2 [4b]
- control law for the mapping (feedback) of MbI, Mb2 to ⁇ bl, ⁇ b2 basically comprises an algorithm based on a time integer:
- ⁇ bl(t) ⁇ bl(t - ⁇ t) + K x Mbl(t) x ⁇ t [5a]
- ⁇ b2(t) ⁇ b2(t - ⁇ t) + K x Mb2(t) x ⁇ t [5b]
- K is a 'controller gain factor' and ⁇ t is a cycle time for the algorithm.
- the control applied has a low frequency: a frequency for adjusting the set blade pitch of each rotor blade according to the specific blade pitch adjustment is lower than the frequency of the component varying in time.
- the choice of the controller gain K in relation to the process gain C, in conjunction with low-pass filtering ensures that the control according to the invention is active on a time scale much larger than the rotation time of the rotor.
- the algorithm comprises at least: (i) a control law for mapping the bending moment values to virtual blade pitch variations and (ii) a trigonometric translation of these virtual blade pitch variations into true variations on the three blade pitches (i.e. the respective blade pitch of each of the rotor blades).
- the feedback of bending moments to virtual blade pitch variations can basically be subjected to any arithmetic operation (in the form of a control law).
- control laws that follow the modem control theory are applied, such as regulator design techniques based upon minimizing a quadratic criterion such as LQR (Linear Quadratic Regulator) and LQG (Linear Quadratic Gaussian), and further to this, the Robust Control technique and techniques based on minimizing the so- called Hoc -standard and refinements thereof, such as ⁇ -synthesis.
- LQR Linear Quadratic Regulator
- LQG Linear Quadratic Gaussian
- control law of the algorithm can be implemented in a processing computer or similar digital technology, but it is not restricted thereto; it is sometimes possible to implement the control law using analogue electronics via capacitors, coils, resistors etc..
- the bending moment values are reconstructed from measurements of deflection moments, possibly in conjunction with measurements of radial and tangential forces on the blade root of each rotor blade using a respective measuring sensor Sl.
- the bending moment values can be measured directly on the rotor shaft by using one or more measuring sensors S2.
- the bending moment values can be deduced by measurements of the tilting and yaw moment in the top of the mast by applying one or more sensors S3; no measurement values of the vertical and lateral force on the yaw bearing are used.
- the tilt moment is designated as Mm t and the yaw moment as whereas the vertical and lateral force is designated by F ve rt and Fi at .
- the horizontal distance from the top of the mast to the cross section of the shaft in which the bending moments are determined, is designated as d.
- MbI cos( ⁇ ) x (Mtii, + d x F vert ) - sin( ⁇ ) x (M yaw - d x FLaO [6a]
- Mb2 sin( ⁇ ) x (M ⁇ u + d x F ve ft) + cos( ⁇ ) x (M yaw - d x FuO [6b].
- acceleration measurements can be applied by using acceleration sensors S4 to compensate these loads for consumed portions of the imbalance forces by the inertia of the rotor and the nacelle, or by directly deducing the loads from these; this is explained in more detail below.
- Measured values of the tilt and yaw moment can also be derived from the measurement of at least one moment along an axis which lies in a plane defined by the direction vectors of the tilting and yaw moment, but which do not coincide with any of these direction vectors.
- the desired measurement values are then derived from the projections from the shaft along which measurements are conducted on shafts parallel to the direction vectors of the tilting and yaw moment.
- Values of the vertical and laterally directed force can also be deduced from the measurement of at least a force along an axis which lies in the plane defined by the vertical and lateral direction vector but which does not coincide with any of these direction vectors.
- MdHjBC Ma, + Jtiu x a vert / r [8]
- the vertical causal force F vm _ rac is then determined by measuring the vertical force F tra ⁇ si and the mass MW of the nacelle and rotor in conjunction with the translational acceleration of the centre of gravity of the assembly of nacelle and rotor, wherein the nacelle also comprises the rotor shaft and generator. Assuming that this lies at a horizontal distance Xc from the yaw bearing and that a ve ⁇ t is measured again at a distance r, then the following methods of determination apply to F ven e ⁇ C :
- M yaw exc Myaw + Jyaw X datat " ajat.top) / t [12];
- Fla t _exc Flat + MW X (ai at _ top + XC X ⁇ yaw) [14]
- M yaw _exc Jyaw X (3
- Myaw_exc Syaw X W Olyaw Xdt2xdtl + J yaw X Otyaw [ 17] or
- Myaw.exc Syaw X W - ai at _u,p) / 1 Xdt2xdtl + Jyaw X (afot - aiatjop) / r [18].
- tl is a first integration variable and t2 is a second integration variable.
- v3 The rigidity of the tower for lateral displacement is known; set this equal to si at .
- the lateral force in the top of the tower is then equal to the product of sj at and a cumulative lateral distortion across the length of the mast.
- the latter can be approximated by the double time interval of the lateral translational acceleration of the top of the tower.
- Fi at exc si a , x U ai at xdt2xdtl + MW x a tat _top + Xc x (a ⁇ at - a ⁇ aUoP )/r [19], or
- an estimator can also be based upon a fully dynamical turbine model that is included in a Kalman filter.
- rotor imbalance is compensated on the basis of a determination method that is insensitive to measurement offset for two bending moment coordinates in the rotor shaft and a decoding method for determining the adjustments to the blade pitches derived from the geometry of the direction vectors of the two bending moments in relation to those of the rotor.
- the intermediate feedback includes filtering, time integration and a scaling for minimizing the (progressive) mean value of the bending moments through low-frequency adjustment of the blade pitches.
- the decoding operation comprises the inverse of the trigonometric relationships for determining shaft bending moments derived from blade deflection moments, as is later presented in an example of an algorithm.
- the determination of 2 shaft bending moment coordinates consecutively comprises (i) an approximation of the axial force from the scaled sum of the deflection moments in the blade roots, (ii) an approximation of the tilt moment deduced from the collective effect on the measured forward-reverse movement of the tower of the previously determined axial force and the tilt moment itself, (iii) low-pass filtering of die tilt moment and (iv) determining two shaft bending moment coordinates along two direction vectors perpendicular to each other on the rotor shaft from the filtered tilt moment via the measured rotor position.
- the bending moment coordinates thus determined are then subjected to at least low-pass filtering, time integration and scaling, after which adjustments of the blade pitches are determined according to the above-described decoding method.
- an axial force measurement can be applied directly.
- bending moments are easier to measure and the rotor blades are generally already equipped with measuring sensors to do this, in view of the condition monitoring and reduction of Ip load variations on the rotor blades.
- the measurement of the blade deflection moments are only used to abstract turbulence effects via axial force from the measured forward-reverse tower movement from which the tilt moment is derived.
- the axial force affects this tower movement, especially at frequencies well below Ip, around the natural frequency of the first bending form, and at frequencies around multiples of the number of rotor blades.
- the arithmetic origin of the Ip component in the approximated tilt moment now only lies in the measured forward-reverse movement of the tower. Acceleration sensors that can be used to conduct these measurements show no offset in the Ip component.
- the lp-component contributes in both the tilting and the yaw moment to the average values of the shaft bending moments, which are a measure of rotor imbalance.
- the approximation of shaft bending moments from the Ip component of only the tilt moment results in average values which amount to half of the values obtained from both rotor moments.
- disregarding the yaw moment results in a 2p component in the shaft bending moments. Both effects on the approximated shaft bending moments do not undermine the application of the method described for imbalance compensation if this is taken into account in the aforementioned scaling and filtering in the feedback.
- a wind turbine is considered with three rotor blades in one rotor plane.
- the rotor blades are individually offset at 120 degrees in the rotor plane.
- a choice was made to influence the aerodynamic conversion by rotating each of the three rotor blades along the respective longitudinal axis of the rotor blade and to regulate the blade pitch.
- the blade pitches are designated as ThI, Th2 and Th3 for the first, second and third rotor blades respectively.
- the bending moment values on the rotor shaft are calculated on the basis of the approximated tilting and bending moment on the rotor with the aid of measuring sensor S4 for the forward-reverse acceleration of the nacelle, measuring sensor S 5 for monitoring the blade pitch and measuring sensor Sl for measuring the deflection moment of the respective rotor blade.
- the tilt moment is approximated from the collective effect on the forward-reverse tower displacement of the approximated axial force and of the tilt moment itself. This is done on the basis of the displacement equation corresponding to the forward-reverse nacelle movement in the first bending form of the tower.
- the deflection moment is designated as MkI, Mk2 and Mk3 for the first, second and third rotor blade respectively; the forward-reverse acceleration of the nacelle is designated as afa; the rotor position as Psi.
- k is a scaling factor.
- the scaling factor k can be determined by a person skilled in the field of aerodynamic conversion using software customary in the wind turbine industry and depends on the progressive average values of wind speed, blade pitches of the rotor blades and the rotational speed of the rotor.
- H is the tower height and m, d, and s are the so-called tower top equivalent mass, damper constant and elasticity constant respectively.
- these tower top equivalent parameters m, d, s can be determined by using software customary in the wind turbine industry. Only the damper constant d depends on the progressive average values of the wind speed, blade pitches of the rotor blades and rotational speed of the rotor. The remaining parameters may be assumed as being constant.
- ⁇ xfa is the progressive average forward-reverse nacelle position with time window T. Therefore, no progressive average value is obtained for the tilt moment. However, this is not a drawback because only the lp-component of the rotor tilt moment is important for the compensation of imbalance.
- Mb I ⁇ Mtilt x cos(Psi) [24a]
- Mb2 - ⁇ Mtilt x sin(Psi) [24b],
- MbI is the bending moment in the horizontal plane and Mb2 is the bending moment in the vertical plane for the rotor pitch of 0 or 180 degrees. It is evident that Ip component in ⁇ Mtilt results in the same mean value in MbI and Mb2, which will be half of the actual mean bending moment values. In addition, a 2 p component in MbI and Mb2 will manifest as a 90 degree shift.
- Mb2 sin(0°)*Mkl + sin(120 o )*Mk2 + sin(240 o )*Mk3 [24d]
- C is a proportional constant, also designated as a process gain factor.
- the decoding method in the compounded blade pitch variations ⁇ Th_l, ⁇ Th_2 and ⁇ Th_3 consists of 'virtual' variations ⁇ Thbl and ⁇ Thb2 according to the inverse of the trigonometric relationships for determining MbI, Mb2 from MkI, Mk2, Mk3:
- ⁇ Th_l ⁇ Thbl x cos(0°) + ⁇ Thb2 x sin(0°) [26a]
- ⁇ Th_2 ⁇ Thbl x cos(120°) + ⁇ Thb2 x sin(120°) [26b]
- ⁇ Th_3 ⁇ Thb 1 x cos(240°) + ⁇ Thb2 x sin(240°) [26c]
- the algorithm for the method comprises:
- ⁇ Thb2(t) ⁇ Thb2(t - ⁇ t) + K x Mb2_f2p(t) x ⁇ t [28b], where K is a 'controller gain factor' and ⁇ t is a cycle time for the algorithm and Mbl_f2p, Mb2_f2p are the proposed approximated shaft bending moment values after filtering the 2p component. Over time, this average value is essentially equal to zero due to the feedback of the measurement signal.
- the 'controller gain factor 1 K should be adjusted to the process gain factor C in equation [25] and the under-approximation of the average shaft bending moments, as indicated in equations [24a] and [24b].
- the controller has a low frequency: The choice of the controller gain factor K in relation to the process gain factor C, in conjunction with low-pass filtering and amplified filtering at about a 2p-frequency ensures that the system according to the invention operates on a time scale much larger than the rotation time of the rotor.
- the algorithm comprises at least: (i) a method for the robust determination of the Ip component in the tilt moment on the rotor, (ii) a control law for mapping bending moment values to virtual blade pitch variations and (iii) a trigonometric translation of the virtual blade pitch variations into real-time variations on the three blade pitches (i.e. the respective blade pitch of each of the rotor blades).
- the bending moments are determined by average without absolute offset, based on a measured nacelle acceleration perpendicular to the rotor plane and measured blade deflection moments that are not necessarily offset-free.
- Figure 6 shows schematic view of a computer for executing computer software in accordance with the method according to the invention.
- the computer 8 comprises a central processing unit 21 with peripheral equipment.
- the central processing unit 21 is connected to memory means 18, 19, 22, 23, 24 which store instructions and data.
- the computer may have one or more reading units 30 (for example, to read floppy disks, CDROMs, DVDs, portable non- volatile memories, etc.), a keyboard 26, and a mouse 27 as input devices and output devices, a display 28 and a printer 29.
- Other input units such as a trackball, a scanner and a touch screen, as well as other output devices, can be provided.
- the central processing unit 21 can be provided with a network adapter 32 for data communications with a network 33.
- the network adaptor 32 is connected to the network 33.
- the network can be any network suitable for data communications.
- the network can be a Local Area Network (LAN) or a Wide Area Network (WAN).
- Other computer systems can be linked to the network 33, which can communicate with the computer 8 via that network connection 32.
- the memory means shown in Figure 6 comprise one or more means selected from RAM 22 (E) EPROM 23, ROM 24, tape unit 19, and hard disk 18. However, there may be more and/or other storage units provided, as will be clear to an expert in the field. Moreover, if necessary, one or more of the memory resources may be placed at a distance from the central processing unit 21.
- the central processing unit 21 is shown as a single unit, but may also comprise various secondary processing units that operate in parallel or are controlled by a single central processing unit. These secondary processing units can be arranged at some distance from each other, as will be known to those skilled in this field.
- the computer 8 comprises an interface 34 for receiving signals from one or more measuring sensors Sl, S2, S3, S4, S5, which are arranged for measuring signals of the deflection moment and possibly the radial and tangential force in the root of each rotor blade (Sl), two bending moments in a cross-section of rotor shaft (S2), the tilt and yaw moment and possibly any vertical and lateral force in the yaw bearing (S3) or the vertical and lateral acceleration in one or more locations at some distance from the yaw bearing (S4) and/or rotational acceleration in the tilting and yaw direction (S4), the rotation angle of the rotor (S5).
- the interface 34 is connected to the central processing unit 21.
- the computer 8 also comprises an interface 35 for receiving a measuring signal of the selected blade pitch of each rotor blade from measuring sensor SO and for transmitting control signals to each of the rotor blades for defining the blade pitches, as visualized in Figure 1 via actuator Al, or for controlling other devices for influencing the aerodynamic conversion, such as micro-tabs, flaps and synthetic jets. A choice can be made not to influence for just one of the rotor blades the aerodynamic conversion for compensating imbalance, as is explained above.
- the interface 35 is connected to the central processing unit 21.
- the computer 8 comprises functionality in hardware and/or software in order to execute the above method.
- the computer is equipped, or is operable in the way of computer software, to perform calculations in accordance with one or more of the aforementioned methods.
- Such computer software stored in or on a computer-readable medium enables the computer, after being loaded from the computer-readable storage into the memory of the computer, to determine the first variable(s) and the second variable(s) according to the present invention.
- the computer is a SCADA system (SCADA: supervisory command and data acquisition) that is suitable for data processing and analysis.
- the computer is arranged to execute the following method for compensating the rotor imbalance in a wind turbine when in use.
- the processing unit of the computer is designed to receive signals with the aid of a measuring sensor (s) in order to translate rotor imbalance into controllable quantities, namely the (progressive) average bending moment values on the rotor shaft, as explained above.
- the measuring sensor(s) record loads in the rotor blades and/or the yaw bearing.
- the processing unit is arranged to record measuring signals of a set blade pitch for each of the rotor blades by means of a respective blade pitch sensor SO of the wind turbine.
- the processing unit is arranged to convert the bending moment values in the rotor shaft to desired adjustments to the blade pitches so that the aerodynamic conversion can be influences so that, essentially, rotor imbalance no longer occurs. This is achieved by minimizing the (progressive) average of the bending moment values.
- die processing unit is designed to accommodate the actuators of each rotor blade to transmit a control signal for adjusting the pitch of each rotor blade according to the specified blade pitch adjustment. It will be evident that the method and system according to the present invention can also be applied to wind turbines which have two or more than three rotor blades. To this end, only the trigonometric relationships need to be adjusted.
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
Method for compensating rotor imbalance in a wind turbine. The wind turbine comprises a rotor and a mast, wherein said rotor comprises a rotor shaft that is provided with an n number of rotor blades. A blade pitch angle of each rotor blade can be individually adjusted by a respective actuator. The method comprises closed-loop repetition of the following steps: - determining and monitoring loads of the wind turbine with the use of a sensor circuit; monitoring a set blade pitch angle for each of the rotor blades by means of a blade pitch sensor of the wind turbine; - determining a quantity indicative of rotor imbalance, based upon the recorded loads; - determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades so that the contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is minimized; effectuating at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
Description
System and method for compensating rotor imbalance in a wind turbine
The present invention relates to a method for compensating rotor imbalance in a wind turbine. The invention also relates to a system for compensating rotor imbalance in a wind turbine. In addition, the invention relates to computer software for implementing said method.
A wind turbine comprises a rotor provided with a number of rotor blades. The rotor is coupled by means of a shaft to an electric generator, possibly via some form of transmission. The electric generator is provided with an output for the output of electrical energy, for example through a power converter to an electricity network. The assembly of rotor, transmission, generator and power conversion unit are arranged in a nacelle on a mast.
Wind turbines always show some rotor imbalance due to deviations in the profile and mass properties of the rotor blades, as well as differences in the blade pitch when the rotor is assembled. The differences in profile characteristics and differences in blade pitch when assembled, result in an aerodynamic imbalance. The differences in mass properties lead to a mass imbalance. The forces and moments resulting from imbalance turn along with the rotation of the rotor.
Among other things, the rotor imbalance results in alternating loads being exerted on the static components. These varying loads are exerted periodically in the azimuth angle of the rotor and are designated as Ip variations that have a period of once every rotation, i.e. they have the same frequency (Ip) as the rotation of the rotor. This therefore leads to variations in the tilting and yaw moments and in the vertical and horizontal forces exerted on the nacelle. In order for the wind turbine to resist these loads, the wind turbine needs to be designed "heavier'.
In addition, variations in tilting Ip and yaw moments often interfere with control circuits, the aim of which is to reduce the blade moments generated by excitations from the flow, for example caused by tower sweeping and by windshear and turbulence, to frequencies in the range of Ip. These control circuits reduce blade loads at about Ip by making adjustments to blade pitch variations almost periodically to the usual course of the blade pitch; in a 3-bladed wind turbine these variations are each offset by 120 degrees. In fact, these blade pitch variations are achieved by Ip modulation of the arithmetically processed tilting and yaw moment. This operation, called a feedback or control law, comprises low-pass filtering, time integration and scaling. It now appears that Ip variations in the tilting and yaw moment can lead to amplified 2p variations (with the double frequency of Ip variations) in the ultimate blade moments. This is then associated with increased 3p variations (with the triple frequency of the Ip- variation) in the tilting and yaw moment.
It is known that the mass imbalance of the rotor can be decreased by static balancing. This is difficult and never fully effective in the dynamic situation during operation, for example due to the occurrence of blade distortions under the effects of the forces of air flow. Also, there is no method available for effectively compensating varying loads in real-time conditions that are caused by aerodynamic imbalance due to the combined effects of variable rotations and fluctuating wind speeds.
It is an objective of the invention to provide a system that mitigates or eliminates one or more of the drawbacks related to the compensation of the rotor imbalance, as known form the prior art. This object is achieved by applying a method for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor and a mast, wherein the rotor comprises a rotor shaft that is provided with an n number of blades, wherein the rotor shaft is connected to a top section of the mast and wherein a pitch of each rotor blade is individually adjustable by means of a respective actuator; wherein the method comprises the execution in real-time of the following steps repeatedly when in use:
- determining and monitoring loads of the wind turbine with the use of a sensor circuit for sensing bending moments in the rotor shaft;
- monitoring a set blade pitch for each of the rotor blades by means of a blade pitch sensor of the wind turbine. - determining the (progressive) average bending moments on the rotor shaft as a quantity indicative of rotor imbalance based upon the monitored loads;
- determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades by minimizing the (progressive) average of the bending moment values so that the contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is kept to a minimum, wherein the determination of at least one correction value for influencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein the feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling;
- effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n- 1 rotor blades.
In an advantageous manner, the forces exerted on each of the rotor blades can be equalized as much as possible, so that a reduction of the rotor imbalance can be achieved under dynamic operating conditions. The invention also enables the exclusion of static balancing or even for it to be performed in a simplified manner, as a result of which the installation of the wind turbine can be carried out with less effort.
In one embodiment the method according to the invention comprises measuring the bending moments as averages without an absolute offset.
In one embodiment the method according to the invention comprises the feedback of the quantity indicative of rotor imbalance, comprising the execution of mathematical operations in accordance with a control law according to a method applied in modem control theory, selected from a group comprising Robust Control, Linear Quadratic Regulation (LQR), Linear Quadratic Gaussian control (LQG), H-infinite control (H∞) and μ-synthesis.
In one embodiment, the method according to the invention comprises at least one correction value for influencing the aerodynamic conversion that relates to a control signal for adjusting a blade pitch value, and effectuating the at least one specific correction value for the transmission of the control signal to the actuator of the respective rotor blade.
In one embodiment the method according to the invention comprises determining and monitoring loads, when said wind turbine is in use, comprising measuring bending moment values in a cross-section of the rotor shaft along two non-parallel vectors.
In one embodiment the method according to the invention includes the wind turbine, further comprising a yaw bearing for pivotally connecting at least the rotor and the mast, wherein the determination of imbalance forces, when said wind turbine is in use, comprises determining a tilt moment and a yaw moment at the yaw bearing.
In one embodiment the method according to the invention comprises the wind turbine, further comprising a yaw bearing for pivotally connecting at least the rotor and the mast, wherein the determination of imbalance forces, when said wind turbine is in use, comprises determining a tilt moment at the yaw bearing. In one embodiment the method according to the invention comprises measuring a vertical and laterally exerted force in the top of the mast.
In one embodiment the method according to the invention comprises the measurement of at least a rotational acceleration along a direction vector in a plane defined by a vertical and lateral direction vector. In one embodiment the method according to the invention comprises the measurement of at least a translational acceleration along a direction vector in a plane defined by a vertical and lateral direction vector at some distance from the centre of the top of the mast.
In one embodiment, the method according to the invention includes the wind turbine comprising a nacelle, and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the determination of the tilt moment and the
yaw moment at die yaw bearing comprises measuring accelerations at the nacelle in the vertical and horizontal plane at some distance from the centre of the top of the mast.
In one embodiment, the method according to the invention includes the wind turbine comprising a nacelle, and the nacelle comprises the rotor shaft which is coupled to the mast via the yaw bearing and wherein the determination of the tilt moment at the yaw bearing comprises measuring accelerations at the nacelle in the vertical and horizontal plane at some distance from the centre of the top of the mast.
In one embodiment the method according to the invention includes the wind turbine comprising a nacelle, and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the bending moments are determined as averages without absolute offset, based upon a measured nacelle acceleration perpendicular to the rotor plane and upon measured blade deflection moments that are not necessarily offset-free.
In one embodiment the method according to the invention comprises that a sensor circuit is selected from a group of at least one sensor and a combination of sensors.
Furthermore, the invention relates to a computer system for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor and a mast, wherein the rotor comprises a rotor shaft that is provided with an n number of rotor blades, wherein the rotor shaft is coupled to a top section of the mast and wherein a blade pitch of each rotor blade is individually adjustable; wherein the wind turbine is provided with at least one sensor for monitoring forces and blade pitch sensors for monitoring a set blade pitch for each of the rotor blades; wherein the computer is provided with a central processing unit and memory, wherein the memory is linked to the processing unit and wherein the processing unit is linked to the sensors, wherein the computer is arranged, when the wind turbine is in use, to repeat the following operations:
- determining and monitoring loads of the wind turbine with the use of a sensor circuit for sensing bending moments in the rotor shaft;
- monitoring a set blade pitch for each of the rotor blades by means of a blade pitch sensor of the wind turbine.
- determining the (progressive) average bending moments on the rotor shaft as a quantity indicative of rotor imbalance based upon the monitored loads;
- determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades by minimizing the (progressive) average of the bending moment values so that die contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is kept to a minimum, wherein the determination of at least one correction value for influencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment
of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling; - effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
The invention further relates to computer software stored on a computer-readable medium for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor and a mast, wherein the rotor comprises a rotor shaft that is provided with an number of rotor blades, wherein the rotor shaft is coupled to a top section of the mast and wherein a blade pitch of each rotor blade is individually adjustable; wherein the wind turbine is provided with at least one sensor for monitoring forces and blade pitch sensors for monitoring a set blade pitch for each of the rotor blades; wherein the computer is provided with a central processing unit and memory, wherein the memory is linked to the processing unit and wherein the central processing unit is linked to the sensors, wherein the computer software comprises executable code which, when loaded on the computer, enables the computer to repeatedly execute in real-time the following operations when the wind turbine is in use:
- determining and monitoring loads of the wind turbine with the use of a sensor circuit for sensing bending moments in the rotor shaft;
- monitoring a set blade pitch for each of the rotor blades by means of a blade pitch sensor of the wind turbine;
- determining the (progressive) average bending moments on the rotor shaft as a quantity indicative of rotor imbalance, based upon the monitored loads; - determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades by minimizing the (progressive) average of the bending moment values so that the contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is kept to a minimum, wherein the determination of at least one correction value for the aerodynamic effect of conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling; - effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
In one embodiment, the computer software according to the invention comprises determining at least one correction value for influencing the aerodynamic conversion of at
least n-1 rotor blades comprising the use of either a method of approximation or method of reconstruction.
The invention further relates to a computer-readable medium that comprises computer executable code which, when loaded on the computer system, enables the computer, as described above, to execute the method as described above.
Further embodiments according to the present invention are described in the subsequent claims.
The invention will now be described in further detail with reference to several drawings illustrating exemplary embodiments thereof. The drawings are intended exclusively to illustrate the objectives of the invention and not to restrict the scope of the invention, which is defined by the appended claims.
In the drawings: figure 1 shows a schematic view of a wind turbine provided with a system in a first embodiment; figure 2 shows a schematic view of a wind turbine provided with a system in a second embodiment; figure 3 shows a schematic view of a wind turbine provided with a system in a third embodiment; figure 4 shows schematic orientations of loads possibly occurring in the wind turbine; figure 5 shows a flow diagram for a method according to one embodiment, and figure 6 shows a schematic view of a computer for executing computer software in accordance with the method according to one embodiment.
In the following figures, the same reference numerals refer to corresponding components in each of the figures. Figure 1 shows a schematic view of a wind turbine Ia. Wind turbine Ia comprises a rotor R with a rotor shaft provided with a number of rotor blades B. The rotor R is coupled to a transmission T for attaching the rotor shaft to an electric generator G by means of the rotor shaft. The electric generator G is provided with an output for the output of electrical energy, for example through a power converter (not shown) to an electricity network. The assembly of rotor R, transmission T and generator G is located in a nacelle L on a mast M. The assembly of rotor R, transmission T and generator G is pivotally connected to the mast M by means of a yaw bearing KS that constitutes part of a yaw system KS. It will be understood that the transmission T may be excluded when the rotor shaft is coupled directly to the generator, i.e. in the case of a low-speed generator. The wind turbine 1 is adjustable by means of a number of parameters in order to achieve an optimal efficiency under certain arbitrary conditions, wherein the parameters are related to the parts R, B, K, G shown and mentioned here in figure 1. Other
parameters could possibly also be of importance, such as would be known to a person skilled in the field.
The present invention provides the insight that Ip variations in the static components of the wind turbine as a result of rotor imbalance can be eliminated by influencing the aerodynamic conversion (the conversion at the rotor blades of energy from the area of sweep to rotational energy of the rotor) of each individual rotor blade. A condition is that the aerodynamic conversion of each rotor blade can be influenced seperately. The aerodynamic conversion can be influenced by rotating the rotor blade along a longitudinal axis and by controlling the blade pitch. The blade pitch of each rotor blade B is adjustable so that a power absorbed from the flow of air by the wind turbine can be controlled. The adjustment of blade pitch by means of an adjustment control system, for example by means of actuators (shown schematically for the upper rotor blade by arrow Al) can be implemented in the blade root of each rotor blade, but also halfway along the longitudinal axis. In essence, any method for influencing the aerodynamic conversion is suitable, as long as the deflection moment in the blade root can be significantly adjusted. The deflection moment is defined as the moment of force in a rotor blade, the vector of which is directed perpendicular to the plane in which the rotor shaft and the longitudinal axis of the rotor blade lie. A sensor SO on each rotor blade is linked to the blade pitch adjustment system for monitoring a set blade pitch value. The yaw system K is arranged pivotally around a yaw bearing KS in order to adjust the position of the rotor R relative to the horizontal direction of the flow of air.
The power converter is provided with an electrical control system that is arranged to control the electric power and rotational speed characteristic of the generator G.
In the first embodiment of the system, each rotor blade is provided with a sensor Sl attached to the blade root in close proximity to the rotor shaft for monitoring a value of the deflection moment in the respective rotor blade. Additionally, in a further embodiment the radial and tangential force in the root of the respective rotor blade can also be measured.
An anemometer W is positioned next to the wind turbine Ia. Figure 2 shows a schematic view of a wind turbine Ib provided with a system in a second embodiment.
In the second embodiment of the system, the rotor shaft of the wind turbine Ib is provided with a sensor S2 for detecting bending moments in the rotor shaft.
Figure 3 shows a schematic view of a wind turbine Ic provided with a system in a third embodiment.
In the third embodiment of the system, the yaw bearing KS of the wind turbine Ic is equipped with a sensor S3 for monitoring the tilting and yaw moment, wherein the tilt
moment vector lies on the horizontal plane and is perpendicular to the shaft and the yaw moment vector in the vertical plane parallel to the mast.
In addition, or alternatively, the system can be provided with one or more acceleration sensors S4 for conducting acceleration measurements from which tilting and yaw moments can be deduced. An acceleration measurement, for example, can be performed by measuring the vertical and lateral nacelle acceleration at a certain distance, in the direction of the rotor R, from the centre of the top of the mast. Both here and in the following description, 'lateral' is designated as a direction in the horizontal plane, perpendicular to the rotor shaft. Also, the third embodiment is provided with a position sensor S5 for monitoring the rotation angle of the rotor (rotor pitch), which is the angle between the longitudinal axis of a selected rotor blade and the horizontal as a result of the rotation of the rotor. Additionally, in a further embodiment, the vertical and lateral force on or close to the yaw bearing can be measured (for example, by sensor S3), or the vertical and lateral acceleration in a second position on the nacelle (for example, by sensor S4).
A condition in the first embodiment is that the signals recorded from the rotor blades are measured without any individual offset. An offset in the measurement results in the minimization of the sum of the actual imbalance and offset. This means that this offset, when compensation is activated, ultimately becomes effective as rotor imbalance. For the measured values recorded on the rotor shaft in the second embodiment, the bending moments are measured as averages without applying any absolute offset. For the measured values recorded at the yaw bearing in the third embodiment, the amplitudes of the Ip variations are both measured offset-free.
If the measuring signals in the second and third embodiment do not meet the above-mentioned respective conditions, then an effect will occur in both embodiments, during active compensation, similar to that described in the first embodiment.
Figure 4 shows, schematically, orientations in a wind turbine of any loads that may occur.
The mast M has a vertical orientation Ol . The KRM vector of the yaw moment is given in the top of the mast, which is directed along the vertical orientation Ol . The vector for the tilt moment TM is directed in the lateral direction (designated by 02, which is perpendicular to the vertical orientation Ol and the rotor shaft RA). The rotation vector RV is directed along the rotor shaft RA. In addition, one rotor blade B is shown. In rotor blade B, the radial force FR (exerted on the rotor centre RC) is directed along the longitudinal axis LA of the rotor blade B. The tangential force FT (exerted on the rotor centre RC) is directed perpendicular to the radial force FR in the plane of the rotor. The vector KL of the deflection moment in the blade root of the rotor blade B is directed perpendicular to the longitudinal axis LA of the rotor blade B.
Figure 5 shows a flow diagram for a method 100 according to one embodiment. When the turbine is started, or after receiving an initialization signal, a cyclic algorithm is activated with an indicative cycle time or period time of 0.1s. This then initiates a first step 105 for initializing, among other things, the adjustments to the control signals for the blade pitches.
Subsequently, in step 110, a measurement is conducted of the deflection moments and possibly the radial and tangential forces in the blade roots (first embodiment Ia; use of sensor Sl), or the bending moments in the rotor shaft (second embodiment Ib; use of sensor (S2), or the tilting and yaw moment and the rotor pitch and possibly the vertical and lateral force exerted on the yaw bearing (third embodiment Ic; use of a combination of sensors (S3, S5) or (S4, S5)).
Then, in step 120, the bending moment values are determined along two vectors perpendicular to each other in the cross-section of the rotor shaft; one of these vectors is parallel to the longitudinal axis of the rotor blade selected for the rotor pitch. These bending moments in embodiment Ia are derived from the measured loads and the rotor geometry, as well as the distance from the rotor centre to the center of the cross-section where the radial and tangential forces also be measured. In embodiment Ib the desired bending moments then result from the measurements and the orientation of the vectors along which measurements are conducted in relation to the vectors of the desired bending moments.
In the third embodiment Ic, the bending moments are deduced from the tilting and yaw moment and the value of the rotor pitch, as well as the distance from the yaw bearing to the cross-section if the vertical and lateral force exerted on the yaw bearing is also measured. In step 130 the bending moment values thus obtained are converted to desired adjustments to the blade pitches from the geometry of the two direction vectors for the bending moments in relation to those of the rotor. This is done with the aim of adjusting the aerodynamic load of the individual rotor blades in such a manner that the (progressive) mean of the bending moments is minimized whilst retaining the quantity of energy generated by the wind. In this respect, the invention is based upon the principle that a relatively small change to a blade pitch of a rotor blade results in a substantially proportional change of the load exerted on the rotor blade by the wind. The conversion of bending moments to adjustments to the blade pitches (feedback) comprises the inverse of the trigonometric relationships for deducing shaft bending moments from blade deflection moments, as incorporated in the later presented example of an algorithm. In addition, the feedback comprises low-pass filtering, time integration and scaling for minimizing the (progressive) mean value of the bending moments via low-frequency adjustment of the blade pitches. Minimum mean bending moments in the shaft automatically imply
minimum rotor imbalance. The blade pitches are adjusted to a time-scale that is greater than approx. three times the cycle time of the rotor.
In step 140, the desired adjustments to the blade pitches derived from step 130 are converted to control signals for the appropriate adjustment of the blade pitches. If a desired adjustment is not possible in connection with the range of a blade pitch, all adjustments are amended with the same, smallest possible offset that does not lie beyond this range. As a result of this offset, the imbalance compensation will have some slight effect on the energy extracted from the wind. It is also possible not to adjust the pitch of one of the blades. Here, the imbalance compensation will generally have a greater effect on the energy extracted from the wind.
In step 150, one first waits till the cycle time has lapsed, after which step 110 is repeated. This is not done after a re-initialization signal 152, after which step 105 is executed, and after a stop signal 154, after which the method is completed by step 160.
Below is an explanation of the method applied to compensate rotor imbalance of the wind turbine, before an example of an algorithm is given.
The objective of the imbalance compensation, for example, is to minimize the Ip variations in the tilting and yaw moment at the position of the yaw bearing. The rotor centre is located at some distance from the yaw bearing KS. The rotor centre is understood to be the point of the rotor R where the rotor blades B meet and are coupled to each other. In one embodiment, the method provides for feedback control for compensation by adjusting the blade pitches based on measurements of load signals in the root of each rotor blade. Based on the measurements of the deflection moments and radial and tangential forces, the bending moment values are determined along two vectors in the (fictive) cross-section of the rotor shaft at the position of the yaw bearing. The rotor imbalance results in mean values of the bending moments that are unequal to zero due to the deflection moments via the trigonometric relationships arising from the rotor geometry and the momentum that arises from the radial and tangential forces as a result of not coinciding with rotor centre and yaw bearing. (Moment is the exerted force multiplied by the arm. The arm is defined as the distance between rotor centre and yaw bearing). In the method described, these bending moments are averaged to zero, as described in step 130. As a result, the Ip variations in the tilt and yaw moment are eliminated at the yaw bearing. The deflection moments and the radial and tangential forces on the three blades will not generally be precisely equal; this will not lead to any Ip variations in the tilt and yaw moment. It will be clear that an embodiment in which measurements of the tilt and yaw moment are translated into bending moment values through the rotor pitch will result in virtually the same compensation of rotor imbalance.
When measurements of deflection moments are applied for reducing lp-blade variations in the blade load, it is then preferable that the compensation control makes the deflection moments on average equal to each other. By now considering the bending moment values along two vectors in the cross section of the rotor shaft at the position of the rotor centre, the compensation is no longer 'distributed' over the deflection moments and radial and tangential forces, but instead the deflection moments are made purely equal to each other whilst the radial and tangential forces are an (uncontrolled) effect of this 'imbalance of deflection moments'; this is due to the fact that blade forces have no effect on the bending moments in the rotor shaft in the rotor centre. This method of compensating imbalance is achieved in an embodiment wherein the bending moments are only deduced from the deflection moments via the trigonometric relationships that result from the rotor geometry.
An embodiment is also conceivable, wherein the bending moments in the rotor centre are measured directly on the shaft. A third embodiment is also conceivable, wherein the bending moments at the rotor centre are composed from measurements of the tilt and yaw moment and the lateral and vertical force at the yaw bearing via the rotor pitch and the distance between the yaw bearing and rotor centre.
An example will now be discussed of an algorithm for compensating rotor imbalance that corresponds to the method according to the present invention. In the example of an algorithm that now follows, a wind turbine is considered with three rotor blades all of which are arranged in a single rotor plane. The rotor blades are each offset at 120 degrees within the rotor plane. Furthermore, a choice was made to influence the aerodynamic conversion by rotating each of the three rotor blades along the respective longitudinal axis of the rotor blade and to regulate the blade pitch. The blade pitches are designated as θl, Θ2 and Θ3 for the first, second and third rotor blades respectively.
In this particular embodiment the bending moment values on the rotor shaft are calculated on the basis of measurements of the deflection moment in the root of each of the rotor blades with the aid of the measuring sensor Sl for measuring the deflection moment of the respective rotor blade.
A wind load exerted on a rotor blade that is directed perpendicular to the rotor plane causes a deflection moment in the rotor blade. The deflection moment is designated as MkI, Mk2 and Mk3 for the first, second and third rotor blade respectively.
The following trigonometric relationships apply to the bending moment values MbI and Mb2 with the deflection moments MkI, Mk2 Mk3 (° stands for degrees):
MbI = MkI x cos(0°) + Mk2 x cos(120°) + Mk3 x cos(240°) [Ia]
Mb2 = MkI x sin(0°) + Mk2 x sin(120°) + Mk3 x sin(240°) [Ib],
wherein MbI is the bending moment along a direction vector that coincides with the rotation of the rotor that points in the direction of the longitudinal axis of the first rotor blade, and wherein Mb2 is the bending moment along a direction vector perpendicular to the above-mentioned direction vector of bending moment MbL
Variations in the blade pitches result in variations in the deflection moments, which are approximately proportional hereto. The following then applies to variations the deflection moment ΔMki and variations in the blade pitch Δθi (i=l ,2,3): ΔMki = C x Δθi [2], wherein C is a proportional constant, also designated as a process gain factor.
The blade pitch variations ΔΘI, ΔΘ2 and ΔΘ3 are composed of 'virtual' variations Δθbl en Δθb2: ΔΘI = Δθbl x cos(0°) + Δθb2 x sin(0°) [3a]
ΔΘ2 = Δθbl x cos(120°) + Δθb2 x sin(120°) [3b]
ΔΘ3 = Δθbl x cos(240°) + Δθb2 x sin(240°) [3c]
Via equations [2], [3a], [3b, [3c] the variations ΔMbl, ΔMb2 in the bending moments MbI, Mb2 then follow: ΔMbl = 3/2 x C x Δθbl [4a]
ΔMb2 = 3/2 x C x Δθb2 [4b]
In this embodiment the algorithm for the method comprises:
- determining the bending moments MbI, Mb2 derived from deflection moments MkI, Mk2, Mk3 with the aid of the equation [ Ia], [Ib] ;
- mapping of MbI, Mb2 to virtual blade pitch variations Δθbl, Δθb2 via a control law;
- applying blade pitch variations ΔΘI, ΔΘ2, ΔΘ3 based upon Δθbl, Δθb2 on the blade pitch θl, Θ2, Θ3 of the respective rotor blade ([3a], [3b], [3c])
In one embodiment the control law for the mapping (feedback) of MbI, Mb2 to Δθbl, Δθb2 basically comprises an algorithm based on a time integer:
Δθbl(t) = Δθbl(t - Δt) + K x Mbl(t) x Δt [5a]
Δθb2(t) = Δθb2(t - Δt) + K x Mb2(t) x Δt [5b], where K is a 'controller gain factor' and Δ t is a cycle time for the algorithm.
Over time, this average value is essentially equal to zero due to the feedback of the time integer from the measurement signal. This is a result familiar to experts in the field. It will be clear that the 'controller gain factor K' will need to be adjusted to the process gain factor C in equation 4.
The control applied has a low frequency: a frequency for adjusting the set blade pitch of each rotor blade according to the specific blade pitch adjustment is lower than the frequency of the component varying in time. The choice of the controller gain K in relation to the process gain C, in conjunction with low-pass filtering ensures that the control according to the invention is active on a time scale much larger than the rotation time of the rotor.
It should be noted that the algorithm comprises at least: (i) a control law for mapping the bending moment values to virtual blade pitch variations and (ii) a trigonometric translation of these virtual blade pitch variations into true variations on the three blade pitches (i.e. the respective blade pitch of each of the rotor blades).
Alternatively, instead of adjusting the three blade pitches, similar options can be applied to influence the aerodynamic conversion of the three rotor blades, such as. synthetic jets, micro-tabs, flaps.
In the algorithm, the feedback of bending moments to virtual blade pitch variations, such as those according to [5a], [5b], can basically be subjected to any arithmetic operation (in the form of a control law).
In these cases, preference is given to certain classes of arithmetic operations. A proper function is achieved, for example, by allowing the time integer to be preceded by low-pass filtering. In addition, control laws that follow the modem control theory are applied, such as regulator design techniques based upon minimizing a quadratic criterion such as LQR (Linear Quadratic Regulator) and LQG (Linear Quadratic Gaussian), and further to this, the Robust Control technique and techniques based on minimizing the so- called Hoc -standard and refinements thereof, such as μ-synthesis.
The control law of the algorithm can be implemented in a processing computer or similar digital technology, but it is not restricted thereto; it is sometimes possible to implement the control law using analogue electronics via capacitors, coils, resistors etc..
An embodiment of a processing computer will be described below with reference to figure 6.
The bending moment values are reconstructed from measurements of deflection moments, possibly in conjunction with measurements of radial and tangential forces on
the blade root of each rotor blade using a respective measuring sensor Sl. Alternatively, the bending moment values can be measured directly on the rotor shaft by using one or more measuring sensors S2.
In one embodiment the bending moment values can be deduced by measurements of the tilting and yaw moment in the top of the mast by applying one or more sensors S3; no measurement values of the vertical and lateral force on the yaw bearing are used. The tilt moment is designated as Mmt and the yaw moment as
whereas the vertical and lateral force is designated by Fvert and Fiat . The horizontal distance from the top of the mast to the cross section of the shaft in which the bending moments are determined, is designated as d.
As a result of rotor imbalance M^, M^w. Fvεπ and Fiat show Ip variations. Determination of the bending moment values of the rotor shaft from Man and Myaw, Fvert and Fiat requires a trigonometric operation to be applied to the rotor pitch ψ in relation to the horizontal plane, i.e. the time integer of the speed of rotation. The following then applies:
MbI = cos(ψ) x (Mtii, + d x Fvert) - sin(ψ) x (Myaw - d x FLaO [6a]
Mb2 = sin(ψ) x (M^u + d x Fveft) + cos(ψ) x (Myaw - d x FuO [6b].
Because the Ip variations are of significance to M,at and Myaw, FV6rt and Fiat as regards compensation of the rotor imbalance, acceleration measurements can be applied by using acceleration sensors S4 to compensate these loads for consumed portions of the imbalance forces by the inertia of the rotor and the nacelle, or by directly deducing the loads from these; this is explained in more detail below.
Measured values of the tilt and yaw moment can also be derived from the measurement of at least one moment along an axis which lies in a plane defined by the direction vectors of the tilting and yaw moment, but which do not coincide with any of these direction vectors. The desired measurement values are then derived from the projections from the shaft along which measurements are conducted on shafts parallel to the direction vectors of the tilting and yaw moment.
Values of the vertical and laterally directed force can also be deduced from the measurement of at least a force along an axis which lies in the plane defined by the vertical and lateral direction vector but which does not coincide with any of these direction vectors.
It is clear that measurements of the vertical and laterally directed (translational) acceleration and the rotational acceleration in the tilting and yaw direction can be determined in a similar manner from acceleration measurements along one axis in such a vertical plane which does not coincide with the vertical or lateral direction vector.
By combining the moment and force measurements in the top of the mast (M) or the yaw bearing KS with acceleration measurements, a more relatively precise approximation can be made of the imbalance forces in the rotor. For example, the inclined tilting movement of the rotor and the nacelle, the 'tilting', consumes part of the imbalance forces through the moment of inertia Jt;i, and the rotational acceleration αtii,. The causal tilt moment Muu.exc is then deduced as follows from the measured tilt moment Mtiit:
Mtiit_exc = Mult + Jtilt X Cltili [7]
This causal tilt moment Mtitt_eχC is then directly related to the imbalance forces.
It will be clear that the causal tilt moment can be determined in a similar manner when measuring a vertical translational acceleration a^si at a horizontal distance r from the yaw bearing:
MdHjBC = Ma, + Jtiu x avert / r [8]
The vertical causal force Fvm_rac is then determined by measuring the vertical force Ftraπsi and the mass MW of the nacelle and rotor in conjunction with the translational acceleration of the centre of gravity of the assembly of nacelle and rotor, wherein the nacelle also comprises the rotor shaft and generator. Assuming that this lies at a horizontal distance Xc from the yaw bearing and that aveιt is measured again at a distance r, then the following methods of determination apply to Fven eλC:
Fv«t_exc = Fvert + MW x Xc x avert / r [9] ; or
Fvert_eXc = Fvert + MW x Xc x αtii, [10]
This causal vertical force Fvert_exc is again directly related to the imbalance forces.
Similar methods of determination apply to the values of the causal yaw moment My3w exc, the causal lateral force F|at_exc, albeit that the lateral acceleration of the top of the mast now also comes into play. Let's assume that this is also measured. Designate this lateral mast acceleration as aiaUop and a second measured value of the translational acceleration as aiat (at distance r). A measurement of the rotational acceleration in the yaw direction is designated as αyaw. The following methods of determination can then be applied:
"lyaw_eχc = Myaw + J yaw X <*yaW
[H] or
Myaw exc = Myaw + Jyaw X fatat " ajat.top) / t [12];
Flat exc = Fiat + MW X (aiaUop + XC X (aIat - ajat top) / r) [13] or
Flat_exc = Flat + MW X (aiat_top + XC X αyaw) [14]
The loads exerted in the yaw bearing can be determined by using methods of approximation based only on acceleration measurements. Three examples vl-v3 are given below: v 1. In the case of a freely yawing wind turbine, no yaw moment can be measured in the yaw bearing; this is always 0. The following then applies: Myaw exc = Jyaw X αyaw 1>5J» or
Myaw_exc = Jyaw X (3|at - ΛhΛjop) / f [ 16] ■
v2. There is no free yawing behaviour but the torsional rigidity of the mast is known; set this equal to syaw. In good approximation the yaw moment in the top of the mast is equal to the product of syaw and a cumulative torsional rotation across the length of the mast. The latter can be approximated by the double time integer of the rotor acceleration in the direction of yaw. The following then applies:
Myaw_exc = Syaw X W Olyaw Xdt2xdtl + Jyaw X Otyaw [ 17] or
wherein tl is a first integration variable and t2 is a second integration variable. v3. The rigidity of the tower for lateral displacement is known; set this equal to siat.
By approximation, the lateral force in the top of the tower is then equal to the product of sjat and a cumulative lateral distortion across the length of the mast. The latter can be approximated by the double time interval of the lateral translational acceleration of the top of the tower. The following then applies:
Fiat exc = sia, x U aiat xdt2xdtl + MW x atat_top + Xc x (aιat - aιaUoP)/r [19], or
Flat_exc = Slat X W *\X Xdt2xdtl + MW X aiat_top + Xc X αyaw [20]
In an embodiment, an estimator can also be based upon a fully dynamical turbine model that is included in a Kalman filter.
In one embodiment rotor imbalance is compensated on the basis of a determination method that is insensitive to measurement offset for two bending moment coordinates in the rotor shaft and a decoding method for determining the adjustments to the blade pitches derived from the geometry of the direction vectors of the two bending moments in relation to those of the rotor.
The intermediate feedback includes filtering, time integration and a scaling for minimizing the (progressive) mean value of the bending moments through low-frequency adjustment of the blade pitches. The decoding operation comprises the inverse of the trigonometric relationships for determining shaft bending moments derived from blade deflection moments, as is later presented in an example of an algorithm. As included in this algorithm, the determination of 2 shaft bending moment coordinates consecutively comprises (i) an approximation of the axial force from the scaled sum of the deflection moments in the blade roots, (ii) an approximation of the tilt moment deduced from the collective effect on the measured forward-reverse movement of the tower of the previously determined axial force and the tilt moment itself, (iii) low-pass filtering of die tilt moment and (iv) determining two shaft bending moment coordinates along two direction vectors perpendicular to each other on the rotor shaft from the filtered tilt moment via the measured rotor position. The bending moment coordinates thus determined are then subjected to at least low-pass filtering, time integration and scaling, after which adjustments of the blade pitches are determined according to the above-described decoding method.
Instead of performing an approximation of the axial force from deflection bending moments, an axial force measurement can be applied directly. However, bending moments are easier to measure and the rotor blades are generally already equipped with measuring sensors to do this, in view of the condition monitoring and reduction of Ip load variations on the rotor blades.
The measurement of the blade deflection moments are only used to abstract turbulence effects via axial force from the measured forward-reverse tower movement from which the tilt moment is derived. The axial force affects this tower movement, especially at frequencies well below Ip, around the natural frequency of the first bending form, and at frequencies around multiples of the number of rotor blades. The arithmetic
origin of the Ip component in the approximated tilt moment now only lies in the measured forward-reverse movement of the tower. Acceleration sensors that can be used to conduct these measurements show no offset in the Ip component.
Furthermore, an offset in a measurement of a deflection moment is therefore not included in the blade pitch adjustments. This therefore prevents under-compensation or over-compensation.
The lp-component contributes in both the tilting and the yaw moment to the average values of the shaft bending moments, which are a measure of rotor imbalance. The approximation of shaft bending moments from the Ip component of only the tilt moment results in average values which amount to half of the values obtained from both rotor moments. Moreover, disregarding the yaw moment results in a 2p component in the shaft bending moments. Both effects on the approximated shaft bending moments do not undermine the application of the method described for imbalance compensation if this is taken into account in the aforementioned scaling and filtering in the feedback. An example will now be discussed of an algorithm for reconstructing or approximating the loads that cause rotor imbalance that corresponds to the method according to the present invention.
The following example of an algorithm a wind turbine is considered with three rotor blades in one rotor plane. The rotor blades are individually offset at 120 degrees in the rotor plane. Furthermore, a choice was made to influence the aerodynamic conversion by rotating each of the three rotor blades along the respective longitudinal axis of the rotor blade and to regulate the blade pitch. The blade pitches are designated as ThI, Th2 and Th3 for the first, second and third rotor blades respectively.
The bending moment values on the rotor shaft are calculated on the basis of the approximated tilting and bending moment on the rotor with the aid of measuring sensor S4 for the forward-reverse acceleration of the nacelle, measuring sensor S 5 for monitoring the blade pitch and measuring sensor Sl for measuring the deflection moment of the respective rotor blade.
The tilt moment is approximated from the collective effect on the forward-reverse tower displacement of the approximated axial force and of the tilt moment itself. This is done on the basis of the displacement equation corresponding to the forward-reverse nacelle movement in the first bending form of the tower.
The deflection moment is designated as MkI, Mk2 and Mk3 for the first, second and third rotor blade respectively; the forward-reverse acceleration of the nacelle is designated as afa; the rotor position as Psi. A variation in the approximated axial force, designated as ΔFax, is proportional to the sum of variations of the three deflection moments according to
ΔFax = k*(ΔMM + ΔMk2 + ΔMk3) [21]
Here, ΔMki is the difference between the momentary measured value for Mki and a progressive average of a third of the sum of MkI, Mk2 and Mk3, for i = 1 ,2,3. k is a scaling factor.
The scaling factor k can be determined by a person skilled in the field of aerodynamic conversion using software customary in the wind turbine industry and depends on the progressive average values of wind speed, blade pitches of the rotor blades and the rotational speed of the rotor.
A customary displacement equation for the forward-reverse nacelle movement in the first bending form of the tower is as follows: m*afa + d*vfa + s*xfa = ΔFax - 3/2 * MtiMT [22]
Here, H is the tower height and m, d, and s are the so-called tower top equivalent mass, damper constant and elasticity constant respectively.
For a skilled person in the field of constructional dynamics, these tower top equivalent parameters m, d, s can be determined by using software customary in the wind turbine industry. Only the damper constant d depends on the progressive average values of the wind speed, blade pitches of the rotor blades and rotational speed of the rotor. The remaining parameters may be assumed as being constant.
The variation ΔMtilt in the rotor tilt moment is approximated on the basis of the measured nacelle acceleration afa and the approximated axial force variation ΔFax based upon the aforementioned displacement equation according to
Here, ~xfa is the progressive average forward-reverse nacelle position with time window T. Therefore, no progressive average value is obtained for the tilt moment. However, this is not a drawback because only the lp-component of the rotor tilt moment is important for the compensation of imbalance.
The following trigonometric relationships apply to the bending moment values MbI and Mb2 with the approximated tilt moment variation ΔMtilt via the rotor pitch Psi:
Mb I = ΔMtilt x cos(Psi) [24a]
Mb2 = -ΔMtilt x sin(Psi) [24b],
Here, MbI is the bending moment in the horizontal plane and Mb2 is the bending moment in the vertical plane for the rotor pitch of 0 or 180 degrees. It is evident that Ip component in ΔMtilt results in the same mean value in MbI and Mb2, which will be half of the actual mean bending moment values. In addition, a 2 p component in MbI and Mb2 will manifest as a 90 degree shift.
In order to understand this, the trigonometric relationships for the correlation between the shaft bending moments MbI and Mb2 and the deflection moments MkI, Mk2 and Mk3 are also given:
MbI = cos(0°)*Mkl + cos( 120°)*Mk2 + cos(240o)*Mk3 [24c]
Mb2 = sin(0°)*Mkl + sin(120o)*Mk2 + sin(240o)*Mk3 [24d]
Variations in the blade pitches lead to variations in the deflection moments which are proportional hereto. The following then applies to the variations of the deflection moment ΔMk_i and variations of the blade pitch ΔTh_i (i=l,2,3):
ΔMk_i = C x ΔTh_i [25],
where C is a proportional constant, also designated as a process gain factor.
The decoding method in the compounded blade pitch variations ΔTh_l, ΔTh_2 and ΔTh_3 consists of 'virtual' variations ΔThbl and ΔThb2 according to the inverse of the trigonometric relationships for determining MbI, Mb2 from MkI, Mk2, Mk3:
ΔTh_l = ΔThbl x cos(0°) + ΔThb2 x sin(0°) [26a]
ΔTh_2 = ΔThbl x cos(120°) + ΔThb2 x sin(120°) [26b]
ΔTh_3 = ΔThb 1 x cos(240°) + ΔThb2 x sin(240°) [26c]
Via equations [25], [26a], [26b, [26c] the variations ΔMbl, ΔMb2 in the bending moments MbI, Mb2 then follow:
ΔMbl = 3/2 x C x ΔThbl [27a]
ΔMb2 = 3/2 x C x ΔThb2 [27b]
In this embodiment the algorithm for the method comprises:
- approximation of the axial variation ΔFax derived from deflection moment measurements MkI, Mk2, Mk3 with the aid of equation [21[;
- approximation of the tilt moment variation ΔMtilt derived from acceleration measurement afa and axial force approximation ΔFax with the aid of equation [23];
- determining the bending moments MbI, Mb2 derived from the approximated tilt moment variation with the aid of the equation [24a], [24b];
- mapping of MbI, Mb2 to virtual blade pitch variations ΔThbl, ΔThb2 via a control law
- applying blade pitch variations ΔTh_l, ΔTh_2, ΔTh_3 based upon ΔThbl, ΔThb2 on blade pitch ThI, Th2, Th3 of the respective rotor blade ([26a], [26b], [26c]);
hi one embodiment the control law for the mapping (feedback) of MbI, Mb2 to Δ Thbl , Δ Thb2 essentially comprises an algorithm based on a time integer, but also at least a filtering of the 2p component in the approximated shaft bending moments: ΔThbl(t) = ΔThbl(t - Δt) + K x Mbl_f2p(t) x Δt [28a]
ΔThb2(t) = ΔThb2(t - Δt) + K x Mb2_f2p(t) x Δt [28b], where K is a 'controller gain factor' and Δ t is a cycle time for the algorithm and Mbl_f2p, Mb2_f2p are the proposed approximated shaft bending moment values after filtering the 2p component. Over time, this average value is essentially equal to zero due to the feedback of the measurement signal.
It will be clear that the 'controller gain factor1 K should be adjusted to the process gain factor C in equation [25] and the under-approximation of the average shaft bending moments, as indicated in equations [24a] and [24b]. The controller has a low frequency: The choice of the controller gain factor K in relation to the process gain factor C, in conjunction with low-pass filtering and amplified filtering at about a 2p-frequency ensures that the system according to the invention operates on a time scale much larger than the rotation time of the rotor.
It should be noted that the algorithm comprises at least: (i) a method for the robust determination of the Ip component in the tilt moment on the rotor, (ii) a control law for mapping bending moment values to virtual blade pitch variations and (iii) a trigonometric translation of the virtual blade pitch variations into real-time variations on the three blade pitches (i.e. the respective blade pitch of each of the rotor blades). hi the above described embodiment the bending moments are determined by average without absolute offset, based on a measured nacelle acceleration perpendicular to the rotor plane and measured blade deflection moments that are not necessarily offset-free.
Figure 6 shows schematic view of a computer for executing computer software in accordance with the method according to the invention.
The computer 8 comprises a central processing unit 21 with peripheral equipment. The central processing unit 21 is connected to memory means 18, 19, 22, 23, 24 which store instructions and data. Furthermore, the computer may have one or more reading units 30 (for example, to read floppy disks, CDROMs, DVDs, portable non- volatile memories, etc.), a keyboard 26, and a mouse 27 as input devices and output devices, a display 28 and a printer 29. Other input units, such as a trackball, a scanner and a touch screen, as well as other output devices, can be provided. Furthermore, the central processing unit 21 can be provided with a network adapter 32 for data communications with a network 33. The network adaptor 32 is connected to the network 33. The network can be any network suitable for data communications. For example, the network can be a Local Area Network (LAN) or a Wide Area Network (WAN). Other computer systems can be linked to the network 33, which can communicate with the computer 8 via that network connection 32.
The memory means shown in Figure 6 comprise one or more means selected from RAM 22 (E) EPROM 23, ROM 24, tape unit 19, and hard disk 18. However, there may be more and/or other storage units provided, as will be clear to an expert in the field. Moreover, if necessary, one or more of the memory resources may be placed at a distance from the central processing unit 21.
The central processing unit 21 is shown as a single unit, but may also comprise various secondary processing units that operate in parallel or are controlled by a single central processing unit. These secondary processing units can be arranged at some distance from each other, as will be known to those skilled in this field. The computer 8 comprises an interface 34 for receiving signals from one or more measuring sensors Sl, S2, S3, S4, S5, which are arranged for measuring signals of the deflection moment and possibly the radial and tangential force in the root of each rotor blade (Sl), two bending moments in a cross-section of rotor shaft (S2), the tilt and yaw moment and possibly any vertical and lateral force in the yaw bearing (S3) or the vertical and lateral acceleration in one or more locations at some distance from the yaw bearing (S4) and/or rotational acceleration in the tilting and yaw direction (S4), the rotation angle of the rotor (S5). The interface 34 is connected to the central processing unit 21.
The computer 8 also comprises an interface 35 for receiving a measuring signal of the selected blade pitch of each rotor blade from measuring sensor SO and for transmitting control signals to each of the rotor blades for defining the blade pitches, as visualized in Figure 1 via actuator Al, or for controlling other devices for influencing the aerodynamic conversion, such as micro-tabs, flaps and synthetic jets. A choice can be made not to influence for just one of the rotor blades the aerodynamic conversion for compensating
imbalance, as is explained above. The interface 35 is connected to the central processing unit 21.
The computer 8 comprises functionality in hardware and/or software in order to execute the above method. The computer is equipped, or is operable in the way of computer software, to perform calculations in accordance with one or more of the aforementioned methods. Such computer software stored in or on a computer-readable medium, enables the computer, after being loaded from the computer-readable storage into the memory of the computer, to determine the first variable(s) and the second variable(s) according to the present invention. In one embodiment the computer is a SCADA system (SCADA: supervisory command and data acquisition) that is suitable for data processing and analysis. hi one embodiment the computer is arranged to execute the following method for compensating the rotor imbalance in a wind turbine when in use.
The processing unit of the computer is designed to receive signals with the aid of a measuring sensor (s) in order to translate rotor imbalance into controllable quantities, namely the (progressive) average bending moment values on the rotor shaft, as explained above. In doing so, the measuring sensor(s) record loads in the rotor blades and/or the yaw bearing.
Also, the processing unit is arranged to record measuring signals of a set blade pitch for each of the rotor blades by means of a respective blade pitch sensor SO of the wind turbine.
Subsequently, the processing unit is arranged to convert the bending moment values in the rotor shaft to desired adjustments to the blade pitches so that the aerodynamic conversion can be influences so that, essentially, rotor imbalance no longer occurs. This is achieved by minimizing the (progressive) average of the bending moment values.
Finally, die processing unit is designed to accommodate the actuators of each rotor blade to transmit a control signal for adjusting the pitch of each rotor blade according to the specified blade pitch adjustment. It will be evident that the method and system according to the present invention can also be applied to wind turbines which have two or more than three rotor blades. To this end, only the trigonometric relationships need to be adjusted.
Alternative and equivalent embodiments of the present invention are conceivable within the scope of the invention, as will be apparent to those skilled in the art. The aim and scope of the invention is restricted only by the appended claims.
Claims
1. Method of rotor imbalance compensation in a wind turbine, wherein the wind turbine comprises a rotor (R) and a mast (M), wherein the rotor comprises rotor shaft that is provided with an n number of blades (B), wherein the rotor shaft is connected to a top section of the mast (M) and wherein a pitch of each rotor blade is individually adjustable by means of a respective actuator (Al), and wherein the method comprises the execution repeated in time of the following steps when said wind turbine is in use:
- determining and monitoring loads of the wind turbine (Ia; Ib; Ic) with the use of the sensor circuit (Sl; S2; S3, S5; S4, S5) for sensing bending moments in the rotor shaft;
- monitoring a blade pitch for each of the rotor blades by means of a blade pitch sensor (SO) of the wind turbine (Ia; !b; !c);
- determining the (progressive) average bending moments on the rotor shaft as a quantity indicative of rotor imbalance, based upon the monitored loads;
- determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades by minimizing the (progressive) average of the bending moment values so that the contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is kept to a minimum, wherein the determination of at least one correction value for infuencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling;
- effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
2. Method according to claim 1 , wherein the bending moments are measured as averages without absolute offset,
3. Method according to claim 1, wherein the feedback of the quantity indicative of rotor imbalance comprises the execution of mathematical operations in accordance with a control law according to a method applied in modern control theory, selected from a group comprising Robust Control, Linear Quadratic Regulation (LQR), Linear Quadratic Gaussian control (LQG), H-infmite control (H00) and μ- synthesis.
4. Method according to any of the claims 1-3, wherein the at least one correction value for influencing the aerodynamic conversion relates to a control signal for adjusting a blade pitch value, and effectuating the at least one specific correction value comprises the transmission of the control signal to the actuator of the respective rotor blade.
5. Method according to claim 1, wherein determining and monitoring loads, when said wind turbine is in use, comprises measuring bending moment values in a cross-section of the rotor shaft along two non-parallel vectors.
6. Method according to claim 1, wherein the wind turbine further comprises a yaw bearing for pivotally connecting at least the rotor and the mast, wherein the determination of imbalance forces, when said wind turbine is in use, comprises determining a tilt moment and a yaw moment at the yaw bearing.
7. Method according to claim 1, wherein the wind turbine further comprises a yaw bearing for pivotally connecting at least the rotor and the mast, wherein the determination of imbalance forces, when said wind turbine is in use, comprises determining a tilt moment at the yaw bearing.
8. Method according to claim 6 or 7, further comprising the measuring of a vertical and laterally exerted force in the top of the mast.
9. Method according to claim 6 or 7 or 8, further comprising the measuring of at least a rotational acceleration along a direction vector in a plane defined by a vertical and lateral direction vector.
10. Method according to any of the claims 6-9, further comprising het measurement of at least a translational acceleration along a direction vector in a plane defined by a vertical and lateral direction vector at some distance from the centre of the top of the mast.
11. Method according to claim 1, wherein the wind turbine comprises a nacelle (L), and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the determination of the tilt moment and the yaw moment at the yaw bearing comprises measuring accelerations at the nacelle in the vertical and horizontal plane at some distance from the centre of the top of the mast.
12. Method according to claim 1, wherein the wind turbine comprises a nacelle (L), and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the determination of the tilt moment at the yaw bearing comprises measuring accelerations at the nacelle in the vertical and horizontal plane at some distance from the centre of the top of the mast.
13. Method according to claim 1 , wherein the wind turbine comprises a nacelle (L), and the nacelle comprises the rotor shaft and is coupled to the mast via the yaw bearing and wherein the bending moments are determined by average without absolute offset, based on a measured nacelle acceleration perpendicular to the rotor plane and measured blade deflection moments that do not need to be offset- free.
14. Method according to claim 1, wherein said sensor circuit is selected from a group of at least one sensor (Sl; S2) and a combination of sensors (S3, S5; S4, S5).
15. Computer system for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor (R) and a mast (M), wherein the rotor comprises a rotor shaft that is provided withan n number of rotor blades (B), wherein the rotor shaft is coupled to a top section of the mast (M) and wherein a blade pitch of each rotor blade is individually adjustable (Al); wherein the wind turbine is provided with at least one sensor (Sl; S2; S3, S5; S4, S5) for monitoring forces and blade pitch sensors (SO) for monitoring a set blade pitch for each of the rotor blades; wherein the computer is provided with a central processing unit (21) and memory (18, 19, 22, 23, 24), wherein the memory is connected to the central processing unit, and wherein the central processing unit is linked to the sensors (SO, Sl ; SO, S2; SO, S3, S5; SO, S4, S5), wherein the computer is arranged, when the wind turbine is in use, to repeat in time the following actions:
- determining and monitoring loads of the wind turbine (Ia; Ib; Ic) with the use of the sensor circuit (Sl; S2; S3, S5; S4, S5) for sensing bending moments in the rotor shaft; - monitoring a blade pitch for each of the rotor blades by means of a blade pitch sensor (SO) of the wind turbine (Ia; !b; !c);
- determining the (progressive) average bending moments on the rotor shaft as a quantity indicative of rotor imbalance, based upon the monitored loads;
- determining at least one correction value for influencing the aerodynamic conversion of at least n- 1 rotor blades by minimizing the (progressive) average of the bending moment values so that the contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is kept to a minimum, wherein the determination of at least one correction value for influencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling;
- effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
16. Computer software stored on a computer-readable medium for compensating rotor imbalance in a wind turbine, wherein the wind turbine comprises a rotor (R) and a mast (M), wherein the rotor comprises a rotor shaft that is provided withan n number of rotor blades (B), wherein the rotor shaft is coupled to a top section of the mast (M) and wherein a blade pitch of each rotor blade is individually adjustable (Al); wherein the wind turbine is provided with at least one sensor (Sl;
S2; S3, S5; S4, S5) for monitoring forces and blade pitch sensors (SO) for monitoring a set blade pitch for each of the rotor blades; wherein the computer is provided with a central processing unit (21) and memory (18, 19, 22, 23, 24), wherein the memory is connected to the central processing unit, and wherein the central processing unit is linked to the sensors (SO, Sl; SO,
S2; SO, S3, S5; SO, S4, S5), wherein the computer software comprises executable code which, when loaded on the computer, enables the computer, to repeatedly in time execute the following operations when the wind turbine is in use:
- determining and monitoring loads of the wind turbine (Ia; Ib; Ic) with the use of the sensor circuit (Sl; S2; S3, S5; S4, S5) for sensing bending moments in the rotor shaft;
- monitoring a blade pitch for each of the rotor blades by means of a blade pitch sensor (SO) of the wind turbine (Ia; Ib; Ic);
- determining the (progressive) average bending moments on the rotor shaft as a quantity indicative of rotor imbalance, based upon the monitored loads;
- determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades by minimizing the (progressive) average of the bending moment values so that the contribution of loads on all of the rotor blades collectively on the quantity indicative of rotor imbalance is kept to a minimum, wherein the determination of at least one correction value for influencing the aerodynamic conversion of at least n-1 blade comprises feedback of the quantity indicative of rotor imbalance as an adjustment of said at least one correction value of at least n-1 rotor blades via an inverse of trigonometric relationships between deflection moments in the rotor blades and bending moments in the rotor shaft, wherein said feedback from the quantity indicative of rotor imbalance comprises low-pass filtering, time integration and scaling; - effectuating the at least one specific corrective value for influencing the aerodynamic conversion of the at least n-1 rotor blades.
17. Computer software according to claim 16, wherein determining at least one correction value for influencing the aerodynamic conversion of at least n-1 rotor blades comprises the use of either a method of approximation or method of reconstruction.
18. A computer-readable medium that comprises computer executable code which, when loaded on the computer system according to claim 15, enables the computer to execute the method according to any of the claims 1-14.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NL2001878A NL2001878C2 (en) | 2008-08-07 | 2008-08-07 | System and method for compensating rotor imbalance for a wind turbine. |
NL2001878 | 2008-08-07 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2010016764A1 true WO2010016764A1 (en) | 2010-02-11 |
Family
ID=40436386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/NL2009/050482 WO2010016764A1 (en) | 2008-08-07 | 2009-08-06 | System and method for compensating rotor imbalance in a wind turbine |
Country Status (2)
Country | Link |
---|---|
NL (1) | NL2001878C2 (en) |
WO (1) | WO2010016764A1 (en) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2482009A (en) * | 2010-07-14 | 2012-01-18 | Vestas Wind Sys As | Ice and fault detection method for wind turbine blades |
CN102338034A (en) * | 2010-05-28 | 2012-02-01 | 通用电气公司 | Method and system for validating wind turbine |
NL2005400C2 (en) * | 2010-09-27 | 2012-03-28 | Stichting Energie | Method and system for wind gust detection in a wind turbine. |
EP2535571A1 (en) * | 2011-06-17 | 2012-12-19 | IFP Energies Nouvelles | Method for optimising the power recovered by wind turbine by reducing the mechanical impact on the structure. |
WO2013182200A1 (en) * | 2012-06-06 | 2013-12-12 | Vestas Wind Systems A/S | A wind turbine with a load controller |
US9032807B2 (en) | 2010-07-14 | 2015-05-19 | Vestas Wind Systems A/S | Method and system for monitoring bending strains of wind turbine blades |
WO2015132187A1 (en) * | 2014-03-05 | 2015-09-11 | Robert Bosch Gmbh | Rotor blade setting method and device for a wind turbine |
EP3015705A1 (en) * | 2014-10-27 | 2016-05-04 | General Electric Company | System and method for adaptive rotor imbalance control |
WO2017000948A1 (en) * | 2015-06-30 | 2017-01-05 | Vestas Wind Systems A/S | Method of calibrating load sensors of a wind turbine |
EP3249217A1 (en) | 2016-05-23 | 2017-11-29 | Envision Energy (Jiangsu) Co., Ltd. | Method of identifying a wind distribution pattern over the rotor plane and a wind turbine thereof |
CN107429663A (en) * | 2015-03-27 | 2017-12-01 | 西门子公司 | The control of wind turbine |
US9970415B2 (en) | 2014-06-12 | 2018-05-15 | General Electric Company | Method and system for managing loads on a wind turbine |
CN110145444A (en) * | 2019-05-09 | 2019-08-20 | 东南大学 | Wind measuring device before a kind of impeller |
WO2020011323A1 (en) * | 2018-07-11 | 2020-01-16 | Vestas Wind Systems A/S | Method and system for controlling a wind turbine to reduce nacelle vibration |
WO2020074331A1 (en) * | 2018-10-10 | 2020-04-16 | Senvion Gmbh | Method and system for operating a wind turbine |
US10669986B2 (en) | 2014-09-01 | 2020-06-02 | Vestas Wind Systems A/S | Relating to the determination of rotor imbalances in a wind turbine |
WO2020125885A1 (en) * | 2018-12-20 | 2020-06-25 | Vestas Wind Systems A/S | Correcting pitch angle |
US10781792B2 (en) | 2017-05-18 | 2020-09-22 | General Electric Company | System and method for controlling a pitch angle of a wind turbine rotor blade |
WO2021129907A1 (en) * | 2019-12-23 | 2021-07-01 | Vestas Wind Systems A/S | A method for controlling tilt moment of a wind turbine with hinged wind turbine blades |
WO2021129906A1 (en) * | 2019-12-23 | 2021-07-01 | Vestas Wind Systems A/S | A method for controlling yaw of a wind turbine with hinged wind turbine blades |
US11608811B2 (en) | 2020-04-08 | 2023-03-21 | General Electric Renovables Espana, S.L. | System and method for mitigating loads acting on a rotor blade of a wind turbine |
EP4224011A1 (en) * | 2022-02-02 | 2023-08-09 | Siemens Gamesa Renewable Energy Innovation & Technology S.L. | Method and device of controlling a wind turbine by individual pitch control |
US11879431B2 (en) | 2019-07-30 | 2024-01-23 | Vestas Wind Systems A/S | Correcting blade pitch in a wind turbine |
US12098704B2 (en) | 2019-10-22 | 2024-09-24 | Ge Infrastructure Technology Llc | System and method for mitigating loads acting on a rotor blade of a wind turbine |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999005414A1 (en) * | 1997-07-25 | 1999-02-04 | Aloys Wobben | Wind energy installation |
US20040151577A1 (en) * | 2003-02-03 | 2004-08-05 | Kirk Pierce | Method and apparatus for wind turbine rotor load control |
JP2005098181A (en) * | 2003-09-24 | 2005-04-14 | Electric Power Dev Co Ltd | Windmill power generation system, windmill direction change program, and recording medium for information readable by computer in which the windmill direction change program is recorded |
US20050276696A1 (en) * | 2004-06-10 | 2005-12-15 | Lemieux David L | Methods and apparatus for rotor blade ice detection |
US20060033338A1 (en) * | 2004-05-11 | 2006-02-16 | Wilson Kitchener C | Wind flow estimation and tracking using tower dynamics |
EP1674724A2 (en) * | 2004-12-23 | 2006-06-28 | General Electric Company | Methods and apparatus for wind turbine fatigue load measurement and assessment |
EP1816347A1 (en) * | 2006-02-01 | 2007-08-08 | Fuji Jukogyo Kabushiki Kaisha | Wind turbine generator |
EP1870596A2 (en) * | 2006-06-19 | 2007-12-26 | General Electric Company | Methods and apparatus for balancing a rotor |
WO2008119351A2 (en) * | 2007-03-30 | 2008-10-09 | Vestas Wind Systems A/S | Wind turbine with pitch control arranged to reduce life shortening loads on components thereof |
EP2060785A1 (en) * | 2007-11-15 | 2009-05-20 | Gamesa Innovation & Technology, S.L. | Method and system for operating sensor |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DK200300670A (en) * | 2003-05-05 | 2004-11-06 | Lm Glasfiber As | Wind turbine with buoyancy regulating organs |
EP1978246A1 (en) * | 2007-04-04 | 2008-10-08 | Siemens Aktiengesellschaft | Method of reducing an unbalance in a wind turbine rotor and device for performing the method |
-
2008
- 2008-08-07 NL NL2001878A patent/NL2001878C2/en not_active IP Right Cessation
-
2009
- 2009-08-06 WO PCT/NL2009/050482 patent/WO2010016764A1/en active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999005414A1 (en) * | 1997-07-25 | 1999-02-04 | Aloys Wobben | Wind energy installation |
US20040151577A1 (en) * | 2003-02-03 | 2004-08-05 | Kirk Pierce | Method and apparatus for wind turbine rotor load control |
JP2005098181A (en) * | 2003-09-24 | 2005-04-14 | Electric Power Dev Co Ltd | Windmill power generation system, windmill direction change program, and recording medium for information readable by computer in which the windmill direction change program is recorded |
US20060033338A1 (en) * | 2004-05-11 | 2006-02-16 | Wilson Kitchener C | Wind flow estimation and tracking using tower dynamics |
US20050276696A1 (en) * | 2004-06-10 | 2005-12-15 | Lemieux David L | Methods and apparatus for rotor blade ice detection |
EP1674724A2 (en) * | 2004-12-23 | 2006-06-28 | General Electric Company | Methods and apparatus for wind turbine fatigue load measurement and assessment |
EP1816347A1 (en) * | 2006-02-01 | 2007-08-08 | Fuji Jukogyo Kabushiki Kaisha | Wind turbine generator |
EP1870596A2 (en) * | 2006-06-19 | 2007-12-26 | General Electric Company | Methods and apparatus for balancing a rotor |
WO2008119351A2 (en) * | 2007-03-30 | 2008-10-09 | Vestas Wind Systems A/S | Wind turbine with pitch control arranged to reduce life shortening loads on components thereof |
EP2060785A1 (en) * | 2007-11-15 | 2009-05-20 | Gamesa Innovation & Technology, S.L. | Method and system for operating sensor |
Non-Patent Citations (5)
Title |
---|
BOSSANYI E A: "INDIVIDUAL BLADE PITCH CONTROL FOR LOAD REDUCTION", WIND ENERGY, WILEY, CHICHESTER, GB, vol. 6, 8 October 2002 (2002-10-08), pages 119 - 128, XP008031928, ISSN: 1099-1824 * |
BOSSANYI E A: "The Design of Closed Loop Controllers for Wind Turbines", WIND ENERGY, WILEY, vol. 3, no. 3, 1 January 2000 (2000-01-01), pages 149 - 163, XP007908706, ISSN: 1095-4244, [retrieved on 20010719] * |
CASELITZ P ET AL: "REDUCTION OF FATIGUE LOADS ON WIND ENERGY CONVERTERS BY ADVANCED CONTROL METHODS", EUROPEAN WIND ENERGY CONFERENCE, XX, XX, 1 October 1997 (1997-10-01), pages 555 - 558, XP008031751 * |
CASELITZ P ET AL: "Rotor condition monitoring for improved operational safety of offshore wind energy converters", JOURNAL OF SOLAR ENERGY ENGINEERING, NEW YORK, NY, vol. 127, no. 2, 1 May 2005 (2005-05-01), pages 253 - 261, XP008084934, ISSN: 0199-6231 * |
TRUDNOWSKI D ET AL: "Independent pitch control using rotor position feedback for wind-shear and gravity fatigue reduction in a wind turbine", PROCEEDINGS OF THE 2002 AMERICAN CONTROL CONFERENCE. ACC. ANCHORAGE, AL, MAY 8 - 10, 2002; [AMERICAN CONTROL CONFERENCE], NEW YORK, NY : IEEE, US, vol. 6, 8 May 2002 (2002-05-08), pages 4335 - 4340, XP010597604, ISBN: 978-0-7803-7298-6 * |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102338034A (en) * | 2010-05-28 | 2012-02-01 | 通用电气公司 | Method and system for validating wind turbine |
CN102338034B (en) * | 2010-05-28 | 2015-02-25 | 通用电气公司 | Wind turbine, method and load control system for validating wind turbine |
GB2482009B (en) * | 2010-07-14 | 2014-07-23 | Vestas Wind Sys As | Ice detection and system for wind turbine blades |
US9523354B2 (en) | 2010-07-14 | 2016-12-20 | Vestas Wind Systems A/S | Ice detection method and system for wind turbine blades |
GB2482009A (en) * | 2010-07-14 | 2012-01-18 | Vestas Wind Sys As | Ice and fault detection method for wind turbine blades |
US9032807B2 (en) | 2010-07-14 | 2015-05-19 | Vestas Wind Systems A/S | Method and system for monitoring bending strains of wind turbine blades |
WO2012044161A3 (en) * | 2010-09-27 | 2012-05-18 | Stichting Energieonderzoek Centrum Nederland | Method and system for wind gust detection in a wind turbine |
NL2005400C2 (en) * | 2010-09-27 | 2012-03-28 | Stichting Energie | Method and system for wind gust detection in a wind turbine. |
FR2976630A1 (en) * | 2011-06-17 | 2012-12-21 | IFP Energies Nouvelles | METHOD FOR OPTIMIZING THE POWER RECOVERED BY A WIND TURBINE BY REDUCING THE MECHANICAL IMPACT ON THE STRUCTURE |
EP2535571A1 (en) * | 2011-06-17 | 2012-12-19 | IFP Energies Nouvelles | Method for optimising the power recovered by wind turbine by reducing the mechanical impact on the structure. |
US10041473B2 (en) | 2011-06-17 | 2018-08-07 | IFP Energies Nouvelles | Method of optimizing the power recovered by a wind turbine by reducing the mechanical impact on the structure |
WO2013182200A1 (en) * | 2012-06-06 | 2013-12-12 | Vestas Wind Systems A/S | A wind turbine with a load controller |
US9970413B2 (en) | 2012-06-06 | 2018-05-15 | Vestas Wind Systems A/S | Wind turbine with a load controller |
WO2015132187A1 (en) * | 2014-03-05 | 2015-09-11 | Robert Bosch Gmbh | Rotor blade setting method and device for a wind turbine |
US9970415B2 (en) | 2014-06-12 | 2018-05-15 | General Electric Company | Method and system for managing loads on a wind turbine |
US10669986B2 (en) | 2014-09-01 | 2020-06-02 | Vestas Wind Systems A/S | Relating to the determination of rotor imbalances in a wind turbine |
EP3015705A1 (en) * | 2014-10-27 | 2016-05-04 | General Electric Company | System and method for adaptive rotor imbalance control |
US9567978B2 (en) | 2014-10-27 | 2017-02-14 | General Electric Company | System and method for adaptive rotor imbalance control |
US10961981B2 (en) | 2015-03-27 | 2021-03-30 | Siemens Gamesa Renewable Energy A/S | Control for a wind turbine |
CN107429663A (en) * | 2015-03-27 | 2017-12-01 | 西门子公司 | The control of wind turbine |
US10677225B2 (en) | 2015-06-30 | 2020-06-09 | Vestas Wind Systems A/S | Method of calibrating load sensors of a wind turbine |
CN107709766B (en) * | 2015-06-30 | 2020-02-07 | 维斯塔斯风力系统集团公司 | Method of calibrating a load sensor of a wind turbine |
WO2017000948A1 (en) * | 2015-06-30 | 2017-01-05 | Vestas Wind Systems A/S | Method of calibrating load sensors of a wind turbine |
CN107709766A (en) * | 2015-06-30 | 2018-02-16 | 维斯塔斯风力系统集团公司 | The method for calibrating the load cell of wind turbine |
EP3249217A1 (en) | 2016-05-23 | 2017-11-29 | Envision Energy (Jiangsu) Co., Ltd. | Method of identifying a wind distribution pattern over the rotor plane and a wind turbine thereof |
US10781792B2 (en) | 2017-05-18 | 2020-09-22 | General Electric Company | System and method for controlling a pitch angle of a wind turbine rotor blade |
WO2020011323A1 (en) * | 2018-07-11 | 2020-01-16 | Vestas Wind Systems A/S | Method and system for controlling a wind turbine to reduce nacelle vibration |
US11572865B2 (en) | 2018-07-11 | 2023-02-07 | Vestas Wind Systems A/S | Method and system for controlling a wind turbine to reduce nacelle vibration |
WO2020074331A1 (en) * | 2018-10-10 | 2020-04-16 | Senvion Gmbh | Method and system for operating a wind turbine |
CN113454335A (en) * | 2018-12-20 | 2021-09-28 | 维斯塔斯风力系统集团公司 | Correcting pitch angle |
WO2020125885A1 (en) * | 2018-12-20 | 2020-06-25 | Vestas Wind Systems A/S | Correcting pitch angle |
CN110145444B (en) * | 2019-05-09 | 2020-12-29 | 东南大学 | Wind measuring device before impeller |
CN110145444A (en) * | 2019-05-09 | 2019-08-20 | 东南大学 | Wind measuring device before a kind of impeller |
US11879431B2 (en) | 2019-07-30 | 2024-01-23 | Vestas Wind Systems A/S | Correcting blade pitch in a wind turbine |
US12098704B2 (en) | 2019-10-22 | 2024-09-24 | Ge Infrastructure Technology Llc | System and method for mitigating loads acting on a rotor blade of a wind turbine |
WO2021129907A1 (en) * | 2019-12-23 | 2021-07-01 | Vestas Wind Systems A/S | A method for controlling tilt moment of a wind turbine with hinged wind turbine blades |
WO2021129906A1 (en) * | 2019-12-23 | 2021-07-01 | Vestas Wind Systems A/S | A method for controlling yaw of a wind turbine with hinged wind turbine blades |
US11608811B2 (en) | 2020-04-08 | 2023-03-21 | General Electric Renovables Espana, S.L. | System and method for mitigating loads acting on a rotor blade of a wind turbine |
EP4224011A1 (en) * | 2022-02-02 | 2023-08-09 | Siemens Gamesa Renewable Energy Innovation & Technology S.L. | Method and device of controlling a wind turbine by individual pitch control |
WO2023148032A1 (en) * | 2022-02-02 | 2023-08-10 | Siemens Gamesa Renewable Energy Innovation & Technology S.L. | Method and device of controlling a wind turbine by individual pitch control |
Also Published As
Publication number | Publication date |
---|---|
NL2001878C2 (en) | 2010-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2010016764A1 (en) | System and method for compensating rotor imbalance in a wind turbine | |
CN107850050B (en) | Wind turbine and method of operating a wind turbine with reduced edgewise vibrations | |
EP3180513B1 (en) | Improvements relating to the determination of rotor imbalances in a wind turbine | |
US10669986B2 (en) | Relating to the determination of rotor imbalances in a wind turbine | |
Bossanyi | The design of closed loop controllers for wind turbines | |
DK2110551T3 (en) | Method and device for forecast-based wind turbine management | |
Boukhezzar et al. | Nonlinear control of variable speed wind turbines without wind speed measurement | |
EP1646786B1 (en) | Method and apparatus for wind turbine rotor load control based on shaft radial displacement | |
CN110678646B (en) | Blade tip clearance, estimation and control of wind turbines | |
EP1907695B8 (en) | Wind flow estimation and tracking using tower dynamics | |
US9863402B2 (en) | System and method for operating a wind turbine based on rotor blade margin | |
US9714085B2 (en) | Monitoring of wind turbines | |
AU2004316333B2 (en) | Wind turbine generator, active damping method thereof, and windmill tower | |
EP2489872B1 (en) | Blade load reduction for wind turbine | |
EP2048562A1 (en) | Method and device for providing at least one input sensor signal for a control and/or monitoring application and control device | |
Kumar et al. | Simulating feedback linearization control of wind turbines using high‐order models | |
CA2707472C (en) | Compensation system for a rotor | |
Andersen et al. | Load alleviation on wind turbine blades using variable airfoil geometry | |
EP3899250B1 (en) | Control of side-side and fore-aft vibrational movement of a wind turbine | |
CN112696313A (en) | System and method for mitigating loads acting on rotor blades of a wind turbine | |
Bottasso et al. | Cascading kalman observers of structural flexible and wind states for wind turbine control | |
Wright et al. | Designing and testing controls to mitigate tower dynamic loads in the controls advanced research turbine | |
Bottasso et al. | Advanced control laws for variable-speed wind turbines and supporting enabling technologies | |
Couchman et al. | Active load reduction by means of trailing edge flaps on a wind turbine blade | |
Elmaati et al. | An integrator based wind speed estimator for wind turbine control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09788266 Country of ref document: EP Kind code of ref document: A1 |
|
DPE1 | Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101) | ||
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 09788266 Country of ref document: EP Kind code of ref document: A1 |