Patent Application
20140003939
Wind farm owners and operators would benefit from knowing the dynamic performance of each wind turbine rotor given the local wind state. If the wind loading of each wind turbine could be ascertained, the maintenance, operation, and control of that turbine could be tailored to maximize uptime (mean time between inspections) and potentially better tracking of the optimal tip speed ratio for maximum energy capture (capacity). For example, wind farm operators could compare stored historical loads estimates of an individual turbine, particularly severe loads, to the assumed design loads to schedule condition-based maintenance for that turbine. In addition, operators could temporarily suppress loads that are primarily responsible for causing the growth of fatigue damage to allow for parts to be ordered in a timely manner (performance-based logistics).
Inflow wind characteristics are currently gathered from anemometers and wind vanes located on the nacelle behind the rotor disk. This position does not lead to an awareness of the inflow due to interference from the rotor blades and nacelle. Wind direction data gathered at this location is used to control yaw position, but this data does not capture real-time wind events. A more accurate measurement of the inflow would allow for more advanced control algorithms to increase energy capture. In particular, increased power production could be achieved using precision yaw control systems. Furthermore, the force distribution on the rotor is of interest from an energy production point of view and the current measurements may provide a perspective of the wind state. Because wind turbines often operate between the cut-in and rated wind speeds, leading to less generated power than their nameplate rating, room for improvement exists to increase power output towards their maximum capability.
Turbines often suffer from yaw error, or a lack of perpendicularity to the oncoming wind flow. Yaw error leads to a decrease of energy capture and subjects the turbine to large fatigue loads. Improper yaw alignment has been shown to create a local rotor blade angle of attack near stall. The addition of tower shadow effects and the cyclic separation and reattachment of flow over the rotor blades create large aerodynamic loads beyond static stall values. Yaw position errors can be contributed to several factors such as hysteresis in drive components, improper yaw sensor mesh, development of backlash in yaw position sensors and long time constants between sensor responses and drive action. Yaw error is almost unavoidable if winds have high directional variance, but an increase in yaw response can improve energy capture by decreasing the time constant of traditional methods. The development of an active yaw control algorithm has been explored using the maximum generated power as a means to determine maximum wind speed direction based on a Hill Climbing Control (HCC) program. This method relied on continuous changes in yaw angle resulting in large time constants for correct perpendicular yaw alignment and an inability to align in non steady state wind flow. One new method uses a Laser Wind Sensor for 3D mapping of the inflow (
Modern utility-scale wind turbines are equipped with several types of sensors for monitoring the wind resource as well as mechanical and electrical variables of interest in the turbine. The performance and reliability of wind turbines are largely governed by the nature of the wind loads that act on the turbine rotor blades. The ability to characterize these loads in real time is advantageous for implementing control algorithms that increase energy capture, which is the primary short-term measure of wind turbine performance. An equally important measure of performance is long-term wind turbine reliability. Reliability is related to the structural integrity of the blade root, low speed shaft, yaw joint, and other load bearing components. A more complete awareness of the temporal and spatial variations in the wind loads that act upon the turbine rotor would allow operators to implement advanced control and maintenance strategies. Such information could also provide an understanding of why similar turbines in the same wind farm experience different failure patterns. Various embodiments of the present invention pertain to how the rotor forced response changes due to yaw and pitch set-point errors and how these changes affect the sensitivity of blade measurements to damage mechanisms. This understanding can be applied to improve energy capture while simultaneously facilitating turbine health management. Utility-scale wind turbines are equipped with anemometers and wind vanes that are located on the nacelle behind the rotor disk to characterize inflow wind conditions. However, this measurement position is not ideal for sensing wind speed due to interference that is created by the blades and nacelle. Furthermore, the wind speed often varies across the rotor disk due to vertical and horizontal wind shear, as a result of the atmospheric boundary layer and other phenomena such as wake flow, but cup anemometers and wind vanes are incapable of measuring these variations. Wind direction data that is gathered by these sensors is used in part to control the yaw position, but this data yields an incomplete perspective of the force distribution on the rotor, which is of most interest for increasing energy capture and monitoring the structural integrity of the turbine.
Because of these sensor errors, turbines often suffer from a lack of perpendicularity to the oncoming wind flow; this is known as yaw error. This condition leads to a decrease in energy capture and large fatigue loads on the turbine. Improper yaw alignment has been shown to create a local rotor blade angle of attack near stall. The addition of tower shadow effects and the cyclic separation and reattachment of flow over the rotor blades create large aerodynamic loads beyond static stall values. The source of yaw control errors can be contributed to several different factors such as hysteresis in drive components, improper yaw sensor mesh, development of backlash in yaw position sensors and long time constants between sensor responses and drive action. Wind turbines can experience yaw error routinely due to the dynamic nature of the wind as it continuously changes direction. If winds have high directional variance, yaw error is almost unavoidable especially on short time scales; however, a small increase in yaw response time could significantly improve energy capture by decreasing the turbine response time constant to variations in wind direction.
The development of an active yaw control algorithm has been explored in the literature using the maximum generated power as a means of determining the maximum wind speed direction based on a Hill Climbing Control (HCC) program. This method relies on continuous changes in yaw angle resulting in large time constants for correct perpendicular yaw alignment and an inability to align in non steady state wind flow. Other methods utilize costly Doppler LIDAR (Light Detection And Ranging) systems, which apply lasers for three-dimensional mapping of the wind inflow to a turbine. These systems focus on the wind state and do not identify the loads experienced by the rotor blades.
In addition to yaw control errors, wind turbines can also suffer from pitch control errors. The use of active pitch control is common for power regulation. A pitch controlled wind turbine uses the power output of the generator to determine the useful for pitch action. When power output increases beyond the rated capacity of the turbine's generator, a command is sent from the Supervisory Command And Data Acquisition (SCADA) system to the blade pitch mechanism to immediately turn (pitch) the rotor blades out of the wind. Conversely, as power output drops (i.e. wind speed decreases) the blades are turned back into the wind. This turning action along the longitudinal axis (pitch axis) is controlled by an electronic or hydraulic pitch motor and gear system. Generally, blades are pitched a few degrees every time the wind changes in order to maintain the optimum angle of attack to maximize the power output for all wind speeds. There are also active stall controlled and passive stall controlled wind turbines with different mechanisms for manipulating the blade pitch set point.
Pitch control is the predominant means of control in current utility-scale wind turbines and is sensitive to power output. A change in pitch angle (on the order of 2 degrees of less) has an impact on the performance and health of a wind turbine. The sensitivity to turbine performance is well illustrated by Burton 2001. The sources of pitch error are similar to those for yaw error and include hysteresis in pitch drive components, development of backlash in pitch position sensors and improper installation of the rotor blade in the field.
Power deficits, however, are not solely the result of incorrectly pitched blades. Yaw error, pitch error, rotor or drivetrain damage, and wind and weather conditions can all lead to reduced power levels.
Another reliability concern for wind turbine rotor blades is blade ice accretion. As the number of wind turbines installed in cold weather climates or regions with cold and wet winters increases, blade ice accretion becomes a major reliability concern. Winter often brings the most favorable wind conditions, but downtime due to blade icing or related damage due to blade icing must be avoided. Without effective ice detection and removal, wind turbines can suffer power reductions as high as 30% per year. In addition to the economic loss, blade icing increases the wind loading experienced by the rotor blade due to a decrease in aerodynamic performance and rotor imbalance that can severely affect the drivetrain. An equally important concern of ice accretion is ice throw. As ice accumulates on the rotor blade in operation, shear and centrifugal forces act on the amassed ice, eventually leading to ice throw. It is most critical to detect ice accretion when the turbine is operating in an idle-speed condition. Here, the aerodynamic and centrifugal forces are small, as the rotational speed of the turbine is near zero. This permits the leading edge of the rotor blade to accumulate more ice at the stagnation point of the airfoil. When wind conditions become favorable for power production the rotational speed increases, the angle of attack changes and large masses of ice are shed. Ice throw fragments up to 16 inches in length and weighing up to two pounds have been recorded over 300 feet from the nearest wind turbine. Ice throw of this magnitude can pose a risk to civil structures and human life.
There are currently no standard solutions on the market for reliable ice detection that can be used as a control input for the turbine's supervisory system. Some measures to prevent icing have been successfully used and deicing methods are in development. However, not all methods operate continuously and there is a need for reliable ice detection to facilitate the activation of the deicing system. Various sensors have been tested but have not performed satisfactorily. There are four methods used for ice detection, but with limitations: 1) Infrared spectroscopy is limited to monitoring one section of the blade and requires the installation of fiber optic cables in the blade; 2) a flexible resonating diaphragm was shown to be effective but requires installation at multiple points inside the blade; 3) ultrasound has been proven effective at detecting ice on aircraft but has not been implemented in wind turbine blades and is not well suited for retro-fitting; 4) a change in capacitance was also measured using wires mounted in the surface of the blade.
What is needed are new methods and apparatus that improve the control of wind turbines. Various embodiments of the present invention do this in novel and unobvious ways.
One aspect of the present invention pertains to a method for control of a wind turbine. Some embodiments include providing a wind turbine including a plurality of blades. Other embodiments include a plurality of sensors, each blade having at least one sensor. Yet other embodiments include observing the signals during operation of the wind turbine. Still others include determining a modal response of at least one blade; and modifying operation of the wind turbine at least in part to change the modal response.
Another aspect of the present invention pertains to a method for control of a wind turbine with blades. Some embodiments include providing a control system for the wind turbine and a sensor attached to at least one of the blades, the sensor providing a signal corresponding to the response of the blade. Other embodiments include reading the signal during operation of the wind turbine. Yet other embodiments include removing the mean value of the signal. Still other embodiments include identifying a blade vibratory mode from the signal; and preparing a variable for the control system and using the variable in control of the wind turbine, the value of the variable being at least partly dependent upon a characteristic of the vibratory mode.
Yet another aspect of the present invention pertains to a method for control of a wind turbine with a non-rotating structure and blades. Some embodiments include providing a control system for the wind turbine, a first sensor on a blade and providing a first signal corresponding to the response of the blade, and a second sensor attached to non-rotating structure of the wind turbine and providing a second signal corresponding to the response of the non-rotating structure. Other embodiments include measuring the first signal and the second signal during operation of the wind turbine. Yet other embodiments include cross-correlating the first signal and the second signal; and preparing a variable for use in the control system, the value of the variable being at least partly dependent upon cross-correlating.
Still another aspect of the present invention pertains to a method for control of a wind turbine having a plurality of blades. Some embodiments include providing a control system for the wind turbine, a sensor attached to each blade and providing a signal corresponding to the response of the blade. Other embodiments include converting into the frequency domain each of the plurality of signals. Yet other embodiments include comparing the frequency content of each blade to the frequency content of each other blade, and automatically controlling the wind turbine based on comparing.
It will be appreciated that the various apparatus and methods described in this summary section, as well as elsewhere in this application, can be expressed as a large number of different combinations and subcombinations. All such useful, novel, and inventive combinations and subcombinations are contemplated herein, it being recognized that the explicit expression of each of these combinations is unnecessary
a) is a photograph of a Micon 65/13 wind turbine with a sensored rotor blade mounted.
b) is a photograph of the testing of a wind turbine rotor blade using modal impact measurements.
a) is a graphical representation of wind speed in turbine axis direction over 180 minute period.
b) is a graphical representation of an autopower spectra of wind speed using 50% overlap processing and Hanning window for 180 minute period.
a) is a graphical representation of autopower spectra of 8m blade flap acceleration using 50% overlap processing and Hanning window for 180 min period.
b) is a graphical representation of corresponding frequency response function magnitudes using along wind speed for shorter period of time.
a) is block diagram representation pertaining to the identification and use of rotor structural response for control and maintenance decision-making according to one embodiment of the present invention.
b) shows a block diagram of a control system according to one embodiment of the present invention.
a) is a photograph of a portion of a wind turbine according to one embodiment of the present invention.
b) is a schematic representation of a wind turbine system according to one embodiment of the present invention.
a): OMA data processing flowchart according to one embodiment of the present invention.
b): a block diagram of a operational modal identification method according to another embodiment of the present invention.
c): a block diagram of a yaw control method according to another embodiment of the present invention.
d): a block diagram of a pitch control method according to another embodiment of the present invention.
e): a block diagram of a damage identification method according to another embodiment of the present invention.
f): a block diagram of an ice accumulation detection method according to another embodiment of the present invention.
For the purposes of promoting an understanding of the principles of the invention, reference will now be made to the embodiments illustrated in the drawings and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the invention is thereby intended, such alterations and further modifications in the illustrated device, and such further applications of the principles of the invention as illustrated therein being contemplated as would normally occur to one skilled in the art to which the invention relates. At least one embodiment of the present invention will be described and shown, and this application may show and/or describe other embodiments of the present invention. It is understood that any reference to “the invention” is a reference to an embodiment of a family of inventions, with no single embodiment including an apparatus, process, or composition that should be included in all embodiments, unless otherwise stated. Further, although there may be discussion with regards to “advantages” provided by some embodiments of the present invention, it is understood that yet other embodiments may not include those same advantages, or may include yet different advantages. Any advantages described herein are not to be construed as limiting to any of the claims.
The use of an N-series prefix for an element number (N______.______) refers to an element that is the same as the non-prefixed element (______.______), except as shown and described thereafter The usage of words indicating preference, such as “preferably,” refers to features and aspects that are present in at least one embodiment, but which are optional for some embodiments. As an example, an element 1110 would be the same as element 110, except for those different features of element 1110 shown and described. Further, common elements and common features of related elements are drawn in the same manner in different figures, and/or use the same symbology in different figures. As such, it is not necessary to describe the features of 1110 and 110 that are the same, since these common features are apparent to a person of ordinary skill in the related field of technology. This description convention also applies to the use of prime (′), double prime (″), and triple prime (′″) suffixed element numbers. Therefore, it is not necessary to describe the features of 20.1, 20.1′, 20.1″, and 20.1′″ that are the same, since these common features are apparent to persons of ordinary skill in the related field of technology. Although various specific quantities (spatial dimensions, temperatures, pressures, times, force, resistance, current, voltage, concentrations, wavelengths, frequencies, heat transfer coefficients, dimensionless parameters, etc.) may be stated herein, such specific quantities are presented as examples only, and further, unless otherwise noted, are approximate values, and should be considered as if the word “about” prefaced each quantity. Further, with discussion pertaining to a specific composition of matter, that description is by example only, and does not limit the applicability of other species of that composition, nor does it limit the applicability of other compositions unrelated to the cited composition.
Various embodiments of the present invention pertain to methods and apparatus for improved control of wind turbines. It has been found that the modal response of wind turbine blades can be used to improve the overall operation of the wind turbine, such as with respect to power generation, detection of the health of the wind turbine, reduction in stresses during operation, and detection of accumulations of ice. Although what will be shown and described are methods and apparatus for acquiring data, processing the data, and controlling a wind turbine, it is appreciated that such methods and apparatus are applicable to other systems other than wind turbines.
In one embodiment, the modal response of the blades is detected by acquiring a signal corresponding to motion of the blades. The signal includes a steady, DC component related to the prevailing average wind velocity. However, the signal also includes modal vibrational data, since the prevailing winds are non-steady. These non-steady flows (or steady flows having a nonuniform distribution across the face plane of the wind turbine) can be considered as impact-type loads occurring once per revolution on the rotor system. In some embodiments the blade motion data is analyzed on a per revolution basis, referred to herein as the order domain.
In another embodiment of the present invention, the frequency response function of a blade is determined in real-time based on motion data for that blade. Preferably, this frequency response function is prepared in the order domain, and then the average frequency response function for a mode is integrated in the order domain from about one half per rotation to 1½ per rotation. It has been found that this integrated value is sensitive to errors in the yaw angle of the wind turbine. Therefore, a control system using this integrated value in a feedback loop can achieve higher power levels by reducing the yaw error.
In another embodiment, the modal response to the blade is compared to the modal response of another blade. In some embodiments, these responses are averaged over a number of revolutions, or over a period of time. If the modal response of one blade is sufficiently different from that of another blade, than one of the blades may be flagged as having damage to it. In yet other embodiments the modal response blade is compared to the historical modal response of the same blade. Damage of the blade can be detected as changes in the modal frequency, changes in the magnitude of the modal response, changes to the half-power band of frequencies for that mode, or changes in the measurement of the phase angle related to that mode.
In still other rent embodiments, the modal response of one or more blades can be used to detect potential accumulation of ice on the blade. For example, the presence of ice on a blade will shift the modal frequencies, especially in consideration of the additional mass of the ice. This additional mass can have a tendency to reduce a modal frequency, and that reduction in frequency can be used as an indicator of ice.
In yet another embodiment, the response of a blade is compared to the response of the nacelle or other structure of the wind turbine. A power spectrum for the blade can be cross compared to the power spectrum of the nacelle. The blade having the highest magnitude, especially within a predetermined frequency band, may be a blade whose pitch angle is too great relative to the other blades.
It has been shown that a decrease in the accuracy of a wind turbine's yaw control due to hysteresis and other factors results in loss in energy capture by the rotor. As a cost-effective means of characterizing the turbine's dynamic response when such losses are experienced, it has been demonstrated that sensors (as one example, inertial sensors) can be used either in the rotor or the nacelle to measure static and dynamic variables that are correlated with swept wind loads on the blades. It was shown using inertial sensors in the rotor that a 5 degree error in the yaw set point of the nacelle resulted in a 1% decrease in the power output and a 70% increase in the dynamic loads to the drive train. Based on these measurements, one embodiment of the present invention pertains to a method in which the yaw set point can be tuned to shape the loads on the rotor to maximize energy capture and maximize the reliability of the turbine components. Various embodiments of the present invention pertain to a data measurement and analysis methodology using integrated blade and nacelle inertial sensors for measurement of a wind turbine's yaw angle error.
One embodiment of the present invention pertains to a methodology for developing online dynamic models for wind turbines and wind farms for use in operational decision-making and automatic control. By using models that relate the wind states at upstream wind turbines to the dynamic performance of rotors at downstream wind turbines, the bandwidth of control actions could be increased.
The potential uses extend beyond condition monitoring of wind turbines and performance monitoring of wind farms. For example, assuming a wind speed of 15 m/s, a modern utility scale wind turbine using rate feedback based on low speed shaft speed would have at least about 3 rotor diameters (or 10 sec) lead time to command its pitch and yaw actuators to suppress gusts. In some embodiments, the lag in response of the turbine control laws due to the flexing of a 25 m long blade with its first natural frequency near 0.3 Hz could be reduced by several seconds assuming (in some embodiments) a single degree of freedom flap of the blade. These speed-ups in anticipatory control would result in a reduction in dynamic loads to the drive train resulting in enhanced reliability in the gearbox, bearings, and generator over the twenty-year life of the turbine.
One methodology of online system identification proposed here involves the creation of potentially multiple-input multiple-output frequency response functions relating wind speed (and possibly direction) variables that are measured at upstream wind turbines or anemometer towers to blade responses or drive train responses of downstream turbines. One embodiment of this method has been illustrated in the context of a Micon 65/13 fixed speed, pitch, and yaw turbine with CX-100 rotor blades shown in
a) shows the wind speed measurement taken 1.5 rotor diameters upstream from the turbine at hub height in the along-wind (rotor axis) direction. The data indicates that over a 3 hr period the mean speed decreases steadily. The autopower spectra shown in
a) shows a set of autopower spectra for the flap direction accelerometer at the 8m span position along the blade. These spectra were calculated at different points in time during the 3 hr data acquisition period using 22 min block sizes of data, which were processed in the same manner as previously described. Note that there are variations in the autopower amplitudes for each of the sequential 22 min datasets.
To capture the effects of these variations on the turbine response for use in implementing control algorithms, frequency response function amplitudes were estimated by normalizing the response spectra with wind speed spectra (
b) is a plot of this frequency response function magnitude, and it suggests that the frequency response functions vary less than the autopowers over some of the frequency range (such as at the rotor speed harmonics). For example, the variation in the frequency response amplitude at 0.92 Hz in
One embodiment of the present invention pertains to apparatus and methods for controlling the yaw angle of a wind turbine that is located downstream of a fluctuating disturbance, such as a second wind turbine. It is understood that the term “upstream” refers to a location receiving energy from the wind before those same atmospheric conditions reach a second (“downstream”) wind turbine. As one example, the two wind turbines could be substantially side by side, yet both be within atmospheric conditions in which the direction of the wind (either steady state, or gusting) is changing. One example would include wind turbines pointed in substantially the same direction for wind energy capture, but in which the direction of the wind changes such that one wind turbine receives a side gust before the adjacent wind turbine receives that same side gust.
Yet another example of a fluctuating disturbance could be a stationary object (such as a building) located in fixed relationship to the downstream wind turbine. In such cases, the fluctuation occurs when the direction of wind changes, such that the fixed object becomes upstream (or upwind) of the wind turbine because of the change in direction of the wind. In some of these embodiments, the downstream wind turbine includes a sensor providing a signal as to the angular orientation of the downstream wind turbine relative to the earth, and further preprogrammed software that recognizes the location of the fixed object relative to the wind turbine and locations on the Earth. Further, yet other embodiments of the present invention contemplate a downstream wind turbine capable of recognizing its angular orientation relative to a field of adjacent wind turbines.
In one embodiment the method includes one or more sensors on at least one blade of the upstream turbine. In yet other embodiments, there are one or more sensors providing information such as the angular orientation of the upstream wind turbine, the rotational speed of the upstream wind turbine, or other information pertaining to the current state of the upstream wind turbine. These sensors are in electrical communication with an electronic controller. The sensors measure the dynamic response of an upstream blade during operation, or other information pertaining to operation of the upstream wind turbine. Examples of the sensors include accelerometers, strain gages, position sensors, velocity sensors, or other sensors capable of providing signals corresponding to acceleration, strain, velocity, position, or the like. Yet other embodiments of the present invention include the use of low frequency accelerometers providing signals that can be used to sense a local initial reference frame from the rotating mass, such as those described in U.S. patent application Ser. No. 12/992,804, incorporated herein by reference.
The electronic controller receives the signal from the sensors and controls a variable of the downstream wind turbine to be controlled. In some embodiments, the controlled variable can be the yaw angle of the downstream turbine, the vibration of any one of the blades of the downstream turbine, the pitch angle of the blades of the downstream turbine, or other actuatable aspects of the downstream turbine.
From any of the control variables discussed above, an error signal can be determined by software within the controller. This error signal is then appropriately filtered, and the filtered results can be used to drive an actuator that repositions the wind turbine in terms of its yaw angle, as one example. As another example, a downstream turbine can begin to change its yaw angle based on movement of an upstream wind turbine, or a change in the relationship between the downstream wind turbine and a fixed object (such that the fixed object begins to aerodynamically “shadow” the downstream wind turbine). As yet another example, a change in speed in an upstream wind turbine, indicating an increase in wind speed, can be used to increase the pitch angle of the blades of a downstream wind turbine in anticipation of higher wind speeds reaching the downstream wind turbine.
Various other embodiments of the present invention include the achievement of greater power output using integrated blade and nacelle inertial sensors to characterize the response associated with swept wind loading on the rotor blades. This information can be used in real-time to observe an increase in dynamic excitation, which can be used to tune the yaw angle and reshape the rotor loads. By applying this methodology, it is shown that the wind turbine power output increases by 12% and the asymmetric rotor fatigue loads decrease. With better loading observations, enhanced maintenance scheduling can also be implemented for individual turbines, thereby decreasing down time and maximizing the life of wind turbine components.
An experimental apparatus according to one embodiment was designed to control the wind state and wind turbine degrees of freedom (rotor pitch, yaw, and longitudinal position) and to measure the rotor response. As shown in
Each of the turbine blades was instrumented with a triaxial, DC-coupled accelerometer and a uniaxial AC-coupled accelerometer. The triaxial and uniaxial accelerometers were attached on the low-pressure side of the blade at a location approximately 14 in and 15 in, respectively, from each of the blade tips (see
The free dynamic response was first investigated by means of a multi-reference modal impact test.
[H(jω)]27×12H[H((jω)]27×12=[V(jω)]12×12[CMIF(jω)]12×12[V(jω)]12×12H (1)
The CMIF was plotted in
Upon establishing an understanding of the basic free response dynamics of the rotor, the axial fans were then utilized to investigate the forced dynamic response and its correlation with the yaw angle. The fans were set at a steady-state speed and the accelerometer response measurements were collected wirelessly while the rotor freely rotated. The yaw position of the turbine was then changed from 0° to 5°, 10°, 20°, and 30° angles by locking the tail of the turbine in a prescribed location with straps.
Response data was first synchronously averaged, as a means of reducing background noise while preserving the rotordynamics of the turbine in the signal. This technique adequately reduced leakage, eliminating the need for the more intricate signal processing windowing techniques. However, it is understood that yet other embodiments of the present invention include the use of windowing techniques or other methods of reducing background noise. Averaging was accomplished by using the signal from the optical tachometer to establish the block size of each revolution. Note that modes in the low frequency range examined in Table 1 include flap-wise motion due to the increased compliance as compared to the edgewise direction and the input force direction normal to the blade chord. Thus, the flap-wise acceleration channels were the DOFs of interest when analyzing the operational data.
To summarize and condense the flap response of the whole rotor at each yaw angle, the magnitudes of the 3 DC-coupled response spectra were summed to produce one datum pertaining to each yaw position. As can be seen from the result shown in
Finally, the turbine's performance was investigated in relationship to the yaw angle. The unrectified output voltage was acquired with the previously mentioned data sets and was plotted in
There are implications of these embodiments on the turbine's performance both from an energy capture and reliability perspective. The data in
A 1.23% loss was measured for a 5° yaw misalignment. The cost influence of such a loss on a 3 MW turbine's revenue production can be estimated. Assuming electricity sales of $100/MWh (10¢/kWh) and continuous 3MW production, a 1.23% loss will decrease the revenue of the energy company by more than $32,000 per turbine per year. In an array of 30 turbines, this loss becomes nearly 1 million dollars. Assuming a more realistic 1.5MW continuous output, it still amounts to a half million dollar loss. In addition, there will be reductions in reliability of the drive train due to the more severe asymmetric loads introduced by this small yaw point error.
One embodiment of the present invention pertains to apparatus and methods for controlling the yaw angle of a wind turbine. In one embodiment the method includes one or more sensors on at least one blade of the turbine. These sensors are in electrical communication with an electronic controller. The sensors measure the dynamic response of the blade during operation. Examples of the sensors include accelerometers, strain gauges, or other sensors capable of providing signals corresponding to acceleration or strain.
The electronic controller receives the signal from the sensors and determines a variable of the wind turbine. In some embodiments, the controlled variable is the vibrational acceleration of the blade, either in broadband terms, or in predetermined frequency ranges. In some embodiments, the frequency ranges are adapted and configured to each include one of more known modal frequencies of the blade. In yet other embodiments, the controlled variable may be the strain measured on the blade, as analyzed in broadband terms or predetermined frequency ranges similar to that as discussed above.
The electronic controller uses the information provided by the sensor to change the yaw angle of the wind turbine in response to the measurement of the controlled variable. In some embodiments, the measurements from the sensors are compared to a preexisting model of the dynamic response of the blade. In yet other embodiments, the currently measured response is compared to a predetermined model that is periodically updated to account for the age or other characteristics of the wind turbine. In yet other embodiments, the controlled variable is compared to other measurements of the blade made at the same moment in time. As one example, the magnitude and/or phase responses in a particular range of frequencies can be compared to the phase and/or magnitude response in other predetermined ranges or frequencies. One example of this would be comparison of response of the blade at the Nth mode to the responses as measured at the N−1 and N+1 modes.
From any of the comparisons discussed above, an error signal can be determined by software within the controller. This error signal is then appropriately filtered, and the filtered results can be used to drive an actuator that repositions the wind turbine in terms of its yaw angle, as one example.
To characterize how the rotor forced response of the turbine changed due to yaw and pitch set-point errors and how these changes affected the sensitivity of blade measurements to damage mechanisms, wind speed measurements were taken on a wind turbine 30 as shown in
Hub 34 and blades 32 are coupled to a nacelle 36 that includes the gear reduction mechanism and electrical generating machinery of turbine 30. Hub 34 and nacelle 36 are coupled to a support beam 31 by a yaw actuator 37. Actuator 37 receives commands from an electronic controller, as does pitch actuator 33, to control the operation of wind turbine 30.
In various embodiments, wind turbine 30 includes various types of sensors. In some embodiments, at least one blade 32 includes a sensor 50 for detection of blade motion. As shown in
In still other embodiments, there are other sensors that provide their signal to an electronic controller 80. In some embodiments, wind turbine 30 includes an anemometer 55 providing a signal corresponding to wind velocity and a tachometer 54 providing a signal corresponding to rotor speed. In still further embodiments, the nonrotating structure of the wind turbine (including the nacelle, portions of the hub, and the support beam 31) have mounted on them a motion sensor 50 which provides its signal to an electronic controller.
A model of wind turbine 30 was placed at the inlet of the test-bed enclosure to quantify the wind profile for horizontal and vertical wind shear as well as a uniform wind condition. A 10 by 10 matrix of wind speeds totaling 100 discrete wind speed data points were sampled using the cup anemometer placed 6 inches from the honeycomb. The data was linearly interpolated between the sampled points and a color bar was added to better highlight the variation in wind speed. The resulting profile for a non-shear wind condition is plotted in
To create a vertical and horizontal wind shear condition, screening material was used to shape the wind inflow at the inlet of the test chamber. For the horizontal (side) shear condition, the screening material was used to restrict airflow on the right side enclosure inlet and create a horizontal gradient in wind speed. For the vertical wind shear condition, the screening material was used to restrict airflow on the lower half of the enclosure inlet and create a vertical gradient in wind speed. The resulting wind speed contours for the negative side shear condition and vertical shear condition are plotted in
Wind turbine rotor blade modes of vibration are excited by wind loads to a greater or lesser degree depending on both the frequency of wind loading and the spatial distribution of that loading. Various embodiments of the present invention identify changes in the operating state of the wind turbine rotor by including a characterization of the free dynamic response of the three-bladed rotor system. Knowledge of the free response of the rotor is useful when attempting to monitor the health of the blades because the operational response is a convolution of the blade free response characteristics and the blade forcing function. The mode shapes and resonant frequencies of the rotor can be the modes of vibration that will be most sensitive to changes in the turbine's operating state (i.e. changes in wind state and/or changes in yaw and pitch set point) and blade damage condition.
To identify the free response characteristics of the wind turbine rotor system, a modal impact test was performed. A multi-reference modal impact test was conducted using three measurement degrees of freedom (DOF) that were local to each rotor blade's coordinate system. The measurement degrees of freedom are denoted as follows: X: span, Y: lead-lag and Z: flap (flap is measured perpendicular to the blade surface).
The multiple-input, multiple-output frequency response function data were analyzed using the Complex Mode Indicator Function, or CMIF. The CMIF is mathematically suited to identify closely spaced modal frequencies by using the singular value decomposition of the normal matrix that was produced using the matrix of frequency response function data. The results of this computation are illustrated in FIG. 3-2. The top three spectral lines of this plot were used to identify degenerate modes of vibration with repeated roots and pseudo-repeated roots. At frequencies shared by these three spectral lines, multiple roots are said to exist. One of the more dominant modes was identified at 8.59 Hz for which there were three roots corresponding to this single mode of vibration; consequently, there were three modal deflection shapes associated with this mode of vibration.
Table 3-1 provides a summary of the damped natural frequencies represented by each peak and multiple peaks found in the CMIF plot. The mode at 8.59 Hz is a first bending mode with a large magnitude of response as seen in
To better categorize the free response behavior, the mode shapes were animated and the behavior of each mode shape was observed. The descriptions provided in Table 3-1 were based on the animation.
Accelerometer measurement data will contain some level of random noise. To improve the accuracy with which the time and frequency data can be analyzed, various embodiments include time-synchronous averaging. Blocks of time-sampled data are synchronized using the optical tachometer pulses so that they each begin at the same angular position of the rotor to eliminate the randomness associated with the differences in phase of the blocks of data that are averaged. Components of the signal that are synchronized with the trigger, which is the optical tachometer pulse in this case, are retained, while random noise is averaged out. A block size of three rotations of the rotor is used throughout this analysis so that the averaged time history is long enough to achieve sufficient frequency resolution, which is equal to the inverse of the time history length. However, yet other embodiments of the present invention are not so constrained, and include a block size comprising a single rotation, and yet other embodiments do not include any time-synchronous averaging. As one example, some embodiments utilize low pass, high pass, and band pass filters to eliminate unwanted signal content.
Further, it is recognized that the synchronous capture of material (by which information is placed in the order domain) can be based on various representations of a complete revolution of the hub and blades. As one example, for purposes of synchronization a once per revolution signal can be acquired from an inertial sensor 50 mounted to a blade. As another example, a revolution can be determined by the repetitive voltage signal produced by the generator of the wind turbine. As yet another example, a revolution of the hub and blades can be established by a magnetic pick-up or Hall effect sensor that is responsive to rotation of the electrical machinery (including gear box), or by a photocell that measures the reflection of a light source (such as a laser) reflecting from a rotating surface.
In order to generate frequency response functions (FRFs) from the accelerometer data that is acquired from the rotor blades, operational modal analysis (OMA) is applied in some embodiments. Two assumptions in using OMA are: (1) the power spectrum of the input force is broadband and smooth, i.e. has no poles or zeroes in the frequency range of interest, and (2) the forcing function is spatially distributed in a uniform manner. Assumption (1) may not particularly applicable to wind-excited structures because the power spectrum of wind is generally dominated by low-frequency components. In general, rotating machinery is self-excited at harmonics of the operating speed, and these harmonics are seen in the OMA frequency response functions that are measured from the turbine that was tested. Assumption (2) is reasonable for a rotating wind turbine, and the measured wind planes in
An OMA method 100 according to one embodiment of the present invention includes first computing the autocorrelation of the synchronously averaged time response for each channel of data. The result is two-sided, containing a positive exponential portion corresponding to the negative poles of the power spectrum, and a decaying exponential portion corresponding to the positive poles. Since both parts of the autocorrelation contain the same information, the positive exponential part is set to zero, essentially zero-padding the time signal. The resulting function is treated as an impulse response function, the discrete Fourier transform of which is the OMA frequency response function. In dynamic systems, the impulse response and frequency response functions are related through the Fourier transform in this manner.
a) schematically describes the OMA process.
Method 100 includes acquiring 110 time-based data. In some embodiments this data is synchronously averaged. Preferably, this information corresponds to motion of the blades, such as acceleration data. Method 100 can be accomplished in a variety of ways, including in the time domain, the frequency domain, or the order domain.
In some embodiments method 100 also includes removing 120 the value of the second is in (demeaning). The mean value being removed in analog fashion by use of a high pass filter on the acceleration data. However, you can also be removed mathematically from the acceleration data during signal processing within controller 90.
Method 100 further includes using a statistical method to determine the repetitive character restricts of the same. In one embodiment, autocorrelation is used for this processing. However, the present invention is not so constrained, and contemplates any method of identifying the vibratory response data corresponding to modal vibration of the blades. In some embodiments the October correlated data is processed in any manner similar to the FFT shift command of MATLAB. This command provides a symmetrical signal as seen in the upper right corner of
Method 100 further includes forcing the second half of the symmetric signal 20. This time response, shown in the lower left for
In various embodiments the present invention include the use of modal vibratory data from a component of the system in the control of the system. Method 100 includes one manner of producing modal information. However, the present invention is not so constrained, and contemplates any manner of generating modal response data, especially in those environments where the vibrating component, such as the blades in a wind turbine, also has applied to it a reasonably steady load, as well as a fluctuating load, such as the load from a non-uniform wind pattern.
A discrete Fourier transform according to one embodiment is carried out in the order domain, which measures frequency relative to the rotational position rather than absolute time (rotations−1 rather than seconds−1). This order domain approach is convenient for analyzing the dynamics of rotating machinery for two reasons: 1) as stated previously, these systems are typically self-excited at harmonics of the rotational frequency, and 2) variation of the independent parameters throughout testing, such as yaw and pitch error, affects the rotational speed of the rotor; therefore, analysis in the order domain allows for comparisons to be easily made between these unequal-speed data sets. However, it is understood that the present invention contemplates application of the filtering and signal processing described herein in any of the time, frequency, or order domains.
In some embodiments, a feature 98 is extracted from the processed data to indicate a change in the parameter being studied. In practice, such features are calculated from raw data and then used by the operator to implement control or maintenance decisions. It is understood that these features can be used to close a control loop in an automatic control system. As examples, a feature 98 can represent a feedback variable, or represent a signal to be further processed into a feedback variable, or can be one of several variables used as feedback in a control system, or a variable used in feedback of automatic control system having multiple closed loops (such as inner and outer loops). Further, in some embodiments, a feature 98 is used to set logical flags in an automatic control system or an automatic monitoring system. Such flags can be used to change the mode of control of the wind turbine, or used to generate information (such as warnings) sent to operators of the wind turbine. Some examples of these features 98 that are described below.
Due to fluctuations in aerodynamic loading, yaw error 98.1 occurs in accelerometer data at a once-per-revolution frequency: 1 rot−1 in the order domain. The result of integrating the OMA FRF from 0.5 to 1.5 rot−1 has been shown to be an indicator of the feature yaw error 98.1, even in the presence of horizontal and vertical wind shear. For analyzing the yaw variation data, a method according to one embodiment of the present invention includes some or all of the following acts:
c) depicts a block diagram for a method 200 according to another embodiment of the present invention. This feature 98.1 can be used for closed loop control and/or monitoring of a wind turbine system, such as system 90. Changes in this feature as a result of controlling the yaw angle of the wind turbine result in changes to the modal response of the blade being measured (or of the average of the blades, if multiple measurements are being taken). Feature 98.1 can be used in a control system according to some embodiments in order to correct the yaw angle of the turbine, and thereby capture a greater amount of energy from the wind stream. However, it is understood that the measured modal response data can be used in algorithms other than method 200 to improve the yaw angle of the wind turbine, based upon modal responses of one or more blades.
The blade flap responses couple with the fore-aft acceleration response of the nacelle, and that coupling varies depending on the pitch angle of each blade. In order to calculate a feature value 98.2 for pitch error, the following procedure is used:
where p3=peak in pitched blade cross power spectrum magnitude at 1 rot−1
p1=peak in non-pitched blade one's cross power spectrum magnitude at 1 rot−1
p2=peak in non-pitched blade one's cross power spectrum magnitude at 1 rot−1
d) identifies a method 300 according to another embodiment of the present invention. Method 300 in one embodiment of the present invention prepares the value of a feature 98.2 useful in correction of pitch error. It is understood that the specific method described in
For damage detection in some embodiments, the time-synchronously averaged operational response data is analyzed in the frequency domain. The averaged linear spectra are calculated using the discrete Fourier transform.
Mode shape analysis is used to identify frequencies in the operational data that correspond to asymmetric bending modes such as modes at 3.1 Hz and 8.6 Hz, but any frequency with a known mode of interest can be evaluated. Trends in the dynamic response at these frequencies can be exploited to reveal a change in the rotor blade condition and how the ability to observe that condition varies with yaw angle. Various embodiments include control systems that implement a change in yaw angle or pitch angle as a result of calculations performed on the dynamic response of the frequencies of interest. Various other embodiments include control systems that set a logical flag or indicator in response to the detection of damage or ice accumulation.
To evaluate these trends in one embodiment, the mode (frequency) of interest selected in the processed operating frequency spectra and the maximum value of the magnitude of the DC acceleration is recorded for a single degree of freedom for each blade within the frequency band. This procedure is carried out for each minute of the 10 min data set and the 10 values are averaged to provide one value for each blade at each yaw angle. The
The magnitude of the acceleration for each blade is used to compare the change in response due to a damage condition and a yaw angle in various wind regimes. In the case of blade root damage, the responses of the blades are compared by finding the difference in response between blade pairs, such as when one blade experiences a reduction in stiffness at the root boundary condition. The reduction in stiffness should cause a change in operational response of this blade and is measured against the healthy blades. For this damage case, the difference in response between healthy and unhealthy blades should, therefore, be greater than the difference between healthy blades.
In some embodiments of the present invention, there is a feature 98.3 (also useful as a variable used in a control system) which is an indicia of the health of the blade (such as whether or not there is any damage present). In some embodiments, the magnitudes of modal responses at a particular mode of interest are compared among the blades during operation. In comparing the magnitude of a modal response among blades, some embodiments include identifying the damaged blade as the blade having the greatest magnitude of response, relative to the other two blades. In yet other embodiments, the damage may also be identified by a blade that has the greatest relative difference in phase angle relative to the other blades. In yet other embodiments, the magnitude of response (or relative phase angles) are compared to baseline data, especially baseline data that includes historical data, including historical trends. In still further embodiments, the damaged blade may be identified by the width of the magnitude of response, such as by identifying the width of the half-power points in a power spectrum, which can indicate a damaged blade by its greater width, and higher damping during operation.
For the case of ice accretion, the magnitude of acceleration is used to track changes in operating response of each blade at a particular frequency of interest and to prepare a feature (or control system variable) 98.4. In some embodiments, a record is made (such as by the electronic controller) of the response of the blades, especially when new, for subsequent use as a baseline for comparison. The magnitude of the response when ice is present on the rotor blades is compared to the magnitude of the response for un-iced blades by calculating the percent error between a historical baseline response and the iced blade response.
As best seen in the bottom of
f) shows a method 500 according to one embodiment of the present invention for detection of ice on a blade using modal information about the blade. Method of 500 is similar to method 400, except the accumulation of ice can be detected as a downward shift in the frequency of a blade mode of vibration.
A block diagram of a control system according to one embodiment of the present invention is shown in
Controller 80 includes an input section 82 that includes the hardware and software 82 that performs signal processing on sensors 50. In some embodiments, signal processing software 82 includes method 100 for identifying modal information about the blades 32 during operation of wind turbine 30. However, it is understood that input software 82 further includes various types of analog and digital filtering, and that the output of software 82 includes multiple parameters provided to signal processor 84, including frequency response functions, real time data, peak magnitudes of selected modes, phase information about selected modes, and the like. The various information gleaned from sensors of the wind turbine are provided to signal processing software 84. Processing software 84 further receives external commands X from the operator, other wind turbines, or the like.
Signal processing software 84 receives inputs from both modal identification software 82 and from the operator. These inputs are further processed. Preferably one or more control variables 98 are provided as outputs to control software 86. It is understood that although reference is being made to different sections of software, such categorization is for explanation only, and is not indicative of any requirements on the software. Any of the various signal processing methods described in this document can be performed anywhere within the software of controller 80, or further, in certain cases, in analog circuitry of controller 80.
As examples of signal processing software 84, in some embodiments there is software performing some or all of the functions of method 500. In such embodiments a flag Y can be set and sent to the operator for his attention. In yet other embodiments, the output of this software can be used later in one or more control algorithms of actuator controller 86.
In one embodiment, signal processing software includes some or all of methods 200, 300, and 400, as described herein. It is appreciated that the descriptions of these methods are not restrictive in nature, and these methods can be expressed in many different software codings.
In various embodiments, control variables 98.1, 98.2, and 98.3, are prepared and provided to controller 86. Generally, variables 98.1, 98.2, and 98.3 are representative of numbers related to blade pitch, power optimization, damage detection, and the like. It is understood that the features 98 previously discussed can be directly used in feedback control (such as being provided to a summing junction), or used in the preparation of one or more variables provided to a summing junction. It is further noted that the reference to a summing junction is indicative of a number used in feedback control, and does not have to represent digital or analog summation. Instead, these features 98 include information in them that is used in control of the wind turbine 30.
Further processing of control variables 98.1, 98.2, and 98.3 is performed in control software 86 by, as examples, software modules 299, 399, or 499, respectively. It is in these software modules that the control variable or feature 98 is used to prepare the command signals that appear as outputs to pitch actuator 33 or yaw actuator 37. It is understood that various types of closed loop controls are envisioned. As one example, controller 80 implements a proportional-integral control methodology and state-space terms. As another example, controller 80 implements command signals to actuators 33 and 37 in terms of a proportional-integral-derivative controller implemented in classical control theory. Various embodiments of the present invention are not restricted to any particular manner of closing the loop to provide stable operation of wind turbine 30.
Statistical analysis is applied in some embodiments to determine how variations in the wind load affect the sensitivity to yaw or pitch error and the sensitivity to rotor damage/condition. The feature 98, x, is normalized in some embodiments to obtain its standard score according to
where ρ and σ are the estimated mean and standard deviation of x, respectively. The standard score is how many standard deviations x is above or below its mean. Because Equation 3 is normalized by σ, changes in x are more readily detected when the feature has little variance for a given operating condition. If x is assumed to be equal to μ±3σ, then the level of change in x from the mean value corresponds to a 99% confidence interval, i.e. the analyst is 99% certain that the feature has undergone a biased change even in the presence of natural variations in x. When this value of x is substituted into Equation 1 and then normalized, the following result is obtained,
expressing how much of a shift in x is required to detect a statistically helpful change in the feature mean with 99% confidence. For example, if x in the flap direction for yaw error exhibits a alp of 0.02, then a 6% change in x (3*0.02) is required useful to achieve 99% confidence that a change due to yaw error has occurred. Therefore, the result of Equation 2 represents a measurement of the sensitivity of x: if it is small, then the sensitivity is high, and that is desirable in terms of condition monitoring because it means that even amidst variations in wind loading and other factors, a change in x indicates a change in the feature.
The yaw error feature is calculated in each of the blade degrees of freedom for the three different wind conditions. These feature values are first plotted versus the yaw error angle to produce one curve for each wind condition and rotor blade measurement degree of freedom in
Several aspects in the presence of yaw error can be extracted from the plots of the flap measurement degree of freedom response and sensitivity seen in
The power extracted from the wind in the test cases involving wind shear is lower than the corresponding case involving no wind shear (refer to
It is observed that wind shear does not hinder but rather enhances the ability to detect yaw error. This increase in sensitivity in the presence of wind shear is seen in
Under horizontal wind shear the sensitivity decreases (the % change in x used to detect a change in the feature value with 99% confidence increases) as the turbine is yawed into the side of the inlet plane with higher velocity wind (positive yaw error). Horizontal wind shear can be particularly prevalent in wind farms due to the presence of wake flows from upstream turbine rotors, which produce velocity deficits on downstream turbine rotors
The span DOF exhibits similar feature curves as the lead-lag DOF in the vertical and horizontal shear conditions, but the uniform flow condition produces a nearly-flat feature curve between ±20°, which is not as desirable because it means that a wide range of yaw errors can have nearly the same feature value. Furthermore, while the feature curves are similar for the vertical and horizontal shear cases,
A method described herein is effective at detecting pitch error 98.2 in the cases of uniform and vertical wind shear conditions, as shown in
The other blade DOFs, lead-lag and span, did not exhibit as much trend in their cross power spectra with any of the nacelle DOFs. The sensitivity of the pitch feature is low for 0° and 5°, especially in the vertical and horizontal wind shear cases (see
By using an approach such as method 400 the change in response due to blade ice accretion was readily identified. The frequency of most interest for damage detection was found to be 3.1 Hz; therefore, damage results correspond to the change in the magnitude of acceleration in the 2 to 4 Hz frequency band of the linear spectra. In the order domain this frequency corresponds to 1 rot−1 response of the turbine rotor-dynamics. Additionally, the leadlag DOF revealed the greatest sensitivity to changes in the rotor-dynamics due to ice accretion.
A similar trend is observed when operating in vertical shear flow.
For the case of pitch error, similar results and trends observed in yaw error were revealed.
For the case of a turbine operating in vertical and horizontal shear regimes, the trends observed in
In the case of blade root damage, a method described herein was applied and revealed that the flap degree of freedom near 7 Hz emphasized the change in response due to the damage condition. In the order domain this frequency corresponds to 2 rot−1 dynamics of the turbine. Recalling
Trends are observed in the blade-to-blade response ratios to identify damage in the presence of changing pitch angle. The change in response amplitude for the blade pairs with the damaged blade (Blade 3) increased significantly beyond a 15° pitch angle. The blade-to-blade response ratios were plotted for the case of no damage and only the pitch angle of Blade 3 was altered.
As one example, the system can be a ship, such as a submarine. As the ship moves in the water and is subjected to wave motion or tidal motion, some component or subsystem of the ship may have a vibratory response to the unsteady or non-uniform excitation presented by the water surrounding the ship. In such cases it may be possible to change the vibratory response of the component or subsystem by use of the writer, or in the case of a submarine, the diving planes. This is because the particular component or subsystem was sensitive to the forcing function, and the movement of the writer or diving planes may be able to alter the environment of the component or subsystem, such that the modal response is beneficially change.
In yet another embodiment the system can be a rocket having engines that produce thrust and which are actually edible to different positions, or different levels of thrust. As the engine thrust changes, the entire vehicle, or a component or subsystem, will respond in one of its vibratory modes. In such cases, it is possible to alter the thrust vector of the engine, or change the angle of a control Finn, so as to provide a compensatory input to the vehicle, component, or subsystem, that provides a beneficial change in the modal response.
In a still further embodiment the system can be a spacecraft having one or more components that are deployable, or having engines used in control of the attitude of the spacecraft. As these components are deployed, or as the engines generate brief periods of thrust, one or more systems or components of the spacecraft will respond with a vibration in one of its modes. In such cases, it is possible to identify the modal response, and implement a beneficial change in that modal response, by use of an actuator on the spacecraft. This compensating actuator may be the same engines or another engine, or an electromechanical actuator such as a piezoelectric actuator mounted to the spacecraft.
Various aspects of different embodiments of the present invention are expressed in paragraphs X1, X2, X3, X4, and X5, as follows:
X1. One aspect of the present invention pertains to a method for control of a wind turbine. The method preferably includes providing a wind turbine including a plurality of blades coupled to a rotatable hub, a plurality of sensors, each blade having at least one sensor, and a controller receiving a signal from each of the sensors. The method preferably includes measuring the signals by the controller during operation of the wind turbine, and determining a modal response of at least one blade. The method preferably includes modifying operation of the wind turbine at least in part to change the modal response.
X2. Another aspect of the present invention pertains to a method for control of a wind turbine with blades. The method preferably includes providing a control system for the wind turbine and a sensor attached to at least one of the blades, the sensor providing a signal corresponding to the vibratory response of the blade. The method preferably includes removing the mean value of the signal and identifying a blade vibratory mode from the demeaned signal. The method preferably includes preparing a variable for the control system and using the variable in control of the wind turbine, the value of the variable being at least partly dependent upon a characteristic of the vibratory mode.
X3. Another aspect of the present invention pertains to a method for control of a wind turbine with a non-rotating structure and blades. The method preferably includes providing a control system for the wind turbine, a first sensor attached to a blade and providing a first signal corresponding to the vibratory response of the blade, and a second sensor attached to non-rotating structure of the wind turbine and providing a second signal corresponding to the vibratory response of the non-rotating structure. The method preferably includes cross-correlating the first signal and the second signal, and preparing a variable for use in the control system, the value of the variable being at least partly dependent upon the cross-correlating.
X4. Another aspect of the present invention pertains to a method for control of a wind turbine having a plurality of blades. The method preferably includes providing a control system for the wind turbine, a sensor attached to each blade and providing a signal corresponding to the vibratory response of the blade. The method preferably includes converting into the frequency domain each of the plurality of signals. The method preferably includes comparing the frequency content of each blade to the frequency content of each other blade. The method preferably includes automatically controlling the wind turbine based on comparing.
X5. A method for controlling a mechanical system, including providing a mechanical system including a plurality of components, an actuator for changing the orientation of a first component relative to a second component, the system being mechanically excited by a non-uniform forcing function, a sensor for providing a signal responsive to the motion of one of the first component or the second component, and a controller receiving the signal and providing commands to the actuator. In some embodiments the method includes identifying in real-time a modal response of one of the first component or the second component. The method can also include determining by the controller a change in the orientation that modifies the modal response. Preferably, the method also includes commanding the actuator by the controller to implement the change.
Yet other embodiments pertain to any of the previous statements X1, X2, X3, X4 or X5, which are combined with one or more of the following other aspects:
Wherein the hub can be yawed relative to the earth modifying includes changing the yaw angle of the hub.
Wherein the blades are coupled to the hub by a pitch control actuator, and modifying includes changing the pitch angle of at least one blade.
Wherein modifying includes identifying a condition of the wind turbine to an operator.
Wherein the condition is ice on a blade, damage to a blade, or load on a blade.
Which further comprises statistically comparing a signal before determining.
Wherein statically comparing includes an autocorrelation of the signal.
Wherein statistically comparing includes cross-correlating the signal with another signal.
Wherein the sensors each provide a signal responsive to at least one of strain, stress, displacement, velocity, or acceleration of the blade.
Wherein the modal response is one of the flap, lead-lag, or span modes.
Wherein determining is in the order domain, frequency domain, or time domain.
Wherein modifying includes a control algorithm having a control loop closed with a characteristic of the modal response.
Wherein the characteristic is a magnitude of the response, phase angle of the response, frequency of the response, or includes a comparison of the modal response with another modal response.
Wherein the characteristic is an integration of the response.
Wherein the vibratory mode is the current vibratory mode, which further comprises providing a historical baseline of the vibratory mode, and preparing includes comparing the current vibratory mode to the baseline vibratory mode.
Which further comprises integrating the blade vibratory mode, and the characteristic is the integrated value.
Wherein the sensor has at least two axes of providing two signals, recording is recording of each signal, removing the mean value is for each signal, and identifying includes averaging the two signals for the mode.
Wherein recording is of the time-domain response of the blade, or the frequency-domain response of the blade.
Wherein recording is for a single complete revolution the wind turbine, or a synchronous average over multiple revolutions of the wind turbine.
Wherein wind turbine generates electricity at a frequency corresponding to revolution of the wind turbine, and which further comprises using the frequency to synchronize recording.
Which further comprises preparing an autocorrelation of the signal before identifying.
Which further comprises ignoring a portion of the autocorrelated signal.
Wherein identifying is by preparing a Fourier transform of the demeaned signal.
Wherein cross-correlating includes preparing the cross-power spectrum of the first signal relative to the second signal, and the variable depends in part upon the cross-power spectrum.
Wherein the non-rotating structure is the nacelle, or the non-rotating structure is the beam supporting the nacelle.
Wherein cross-correlating includes determining a peak response.
Wherein cross-correlating is in the order domain, and the peak is proximate to once per revolution of the wind turbine
Which further comprises removing the mean value of the first recorded signal before cross-correlating.
Wherein cross-correlating includes preparing a Fourier transform of the cross-correlation.
Wherein preparing a Fourier transform is in the order domain.
Wherein the first signal corresponds to the flap response of the blade.
Wherein the second signal corresponds to the fore-aft response of the non-rotating structure.
Wherein automatically controlling includes setting a flag for the attention of the operator, or shutting down operation of the wind turbine.
Wherein comparing is in a predetermined range of frequencies.
Wherein comparing includes comparing the peak magnitude of each blade to the peak magnitude of each other blade at a predetermined frequency.
Wherein automatically controlling includes identifying the blade that has the greatest peak magnitude.
Wherein comparing includes comparing the current frequency of the peak magnitude of each blade to a predetermined frequency for each blade.
Wherein automatically controlling includes identifying an icing condition if there is a reduction in the current frequency for more than one blade.
wherein the forcing function is the wind, or a rocket engine, or waves in a body of water, or a change in the orientation of third component relative to one of the first or second components.
wherein the system is an aircraft, or a building, a ship, or a rocket, or a spacecraft.
wherein one of the components is one of a solar panel, instrument boom, or antenna
wherein the actuator is a piezoelectric actuator.
While the embodiments have been illustrated and described in detail in the drawings and foregoing description, the same is to be considered as illustrative and not restrictive in character, it being understood that only certain embodiments have been shown and described and that all changes and modifications that come within the spirit of some embodiments of the invention are desired to be protected.
This application claims the benefit of priority to U.S. Provisional Patent Application Ser. No. 61/452,891, filed Mar. 15, 2011, entitled LOAD SHAPE CONTROL OF WIND TURBINES, incorporated herein by reference.
This invention was made with government support under DE-EE0003265 awarded by the U.S. Department of Energy. The government has certain rights in the invention.
| Filing Document | Filing Date | Country | Kind | 371c Date |
|---|---|---|---|---|
| PCT/US2012/029254 | 3/15/2012 | WO | 00 | 9/13/2013 |
| Number | Date | Country | |
|---|---|---|---|
| 61452891 | Mar 2011 | US |
