Detection of blade root moment sensor failures is an important problem for faulttolerant individual pitch control, which plays a key role in reduction of uneven blade loads of large wind turbines. A new method for detection of blade root moment sensor failures which is based on variations induced by a vertical wind shear is described in this paper. The detection is associated with monitoring of statistical properties of the difference between amplitudes of the first harmonic of the blade load, which is calculated in two different ways. The first method is based on processing of the load sensor signal, which contains a number of harmonics. The first harmonic is recovered via least squares estimation of the blade load signal with harmonic regressor and strictly diagonally dominant (SDD) information matrix. The second method is a modelbased method of estimation of the first harmonic, which relies on the blade load model and upwind speed measurements provided by multibeam Light Detection and Ranging (LIDAR). This is a new application for future LIDARenabled wind turbine technologies. Moreover, adaptation of the load model in a uniform wind field is proposed. This adaptation improves accuracy of the load estimation and hence the performance of the blade load sensor failure detection method.
Individual pitch control (IPC) is the only suitable tool for mitigation of uneven blade loads [19]. IPC reduces maintenance costs and increases efficiency and lifetime of the turbine components through significant contribution to load reduction. IPC is based on the blade load measurements provided by strain gauge resistors, where the resistance changes in the case of loading. The lifetime of the strain gauge is normally not very long and the main sources of the errors are the following: (1) sensitivity to temperature variations which results in a zero load drift, temperature distribution over the material and thermal stresses; (2) transverse sensitivity; (3) humidity which results in variation in stiffness of the composite materials; (4) strain cycling which introduces a zero shift; (5) fatigue of the material which results in a zero shift, change in gauge factor, and possible gauge failure in fatigue; (6) cable effects that might change the resistance, capacitance, insulation, and screening.
The factors listed above together with many other factors are the reasons for the high failure rates of blade load sensors. The failure rates are especially high for offshore wind turbines due to severe environmental conditions. An average failure rate is one failure per year. For wind turbines with three blades, the failure rate is approximately three failures per year [10].
Performance improvement of the blade load sensor failure detection mechanism is a challenge in the IPC system, which plays a key role in reduction of uneven blade loads of large wind turbines. On the other hand, new laser sensor Light Detection and Ranging (LIDAR) technologies which are capable of measuring wind speed at a distance in front of the turbine will be widely used in future in a number of turbine control loops aiming to performance improvement via integration of proactive capabilities. This will result in LIDARassisted proactive turbine speed control, collective pitch angle control as well as yaw and individual pitch control, see [1114] for details. Preview information, provided by the LIDAR usually results in a feedforward part which is properly integrated into the existing feedback controller. Information provided by the LIDAR might also be used for fault detection purposes (so far it was used for control only). LIDAR sensors are traditionally configured with a single beam. New generation of LIDAR sensors are equipped with many beams at different angles providing new opportunities for accurate estimation of incoming wind field [15, 16]. In particular, LIDARenabled blade load estimation may result in a new class of blade load fault detection methods with improved detection capabilities. This is a new application for future wind turbine technologies based on multibeam LIDAR.
Notice that the load model driven by upwind speed measurements is suitable for estimation of the first harmonic only, whereas load sensor signal contains a number of harmonics. Therefore, a reduction of uneven loads is more efficient for IPC, which is driven by the load sensor signal, provided that a pitch actuator is fast enough to counteract higher harmonics [8]. IPC which is based on LIDAR measurements [9] is used in the system in the case of blade load sensor failure. Information provided by LIDAR ensures fault tolerance of IPC as an ability of the system to maintain control objectives, despite the occurrence of a failure in the blade root moment sensor.
A new method for detection of blade root moment sensor failures which is based on a vertical wind shear estimation is described in this paper. Wind shear is almost always present in the wind speed distribution (see Fig. 2 in [9]), which illustrates a vertical wind shear. Traditional detection methods are usually based on monitoring of a change in mean value of the sensor signal (see [10] for example), whereas turbine cycle variations act as a disturbance, leading to deterioration of the detection performance and misdetection. This detection method is based on monitoring of the turbine cycle variations and captures information, which was previously ignored. Assessment of turbine cycle variations provides the basis for fast detection methods with respect to the methods, which are based on monitoring of the mean value. Fast detection is especially efficient in the case of drifting sensor failures.
The method is based on monitoring of the evolution of the first harmonic of the signal in the turbine cycle loading, which appears due to the wind shear. Amplitude of the first harmonic is estimated in two ways, which are described as follows: (1)The first harmonic can be extracted from the blade root moment sensor signal, which contains a number of harmonics and a measurement noise. A number of algorithms can be used for estimation of the first harmonic, taking into account a property of persistence of excitation of the blade root moment sensor signal in the presence of the wind shear. This method is associated with processing of the blade load signal using leastsquares algorithm with harmonic regressor [1721] for estimation of amplitude of the first harmonic. (2) LIDARenabled blade load estimation opens new opportunities for the modelbased blade load sensor fault detection. The first harmonic of the turbine cycle loading can be estimated using upwind speed measurements in multibeam configuration which is illustrated in Figure 1. This method is referred as a modelbased blade load estimation method, driven by upwind speed measurements.

Figure 1. Multibeam LIDAR measurements of two wind speeds V1 and V2 (located at the heights which are equal to the hub height H with the half of the blade length L added/subtracted) at a distance D in front of the turbine. A periodic loading appears on the blades with the rotation of the turbine rotor in the presence of a vertical wind shear, which is associated with a change (increase) in wind speed with height, V1 > V2. This periodic loading has a number of harmonics and a first harmonic is associated with the rotational frequency.

Estimation of the first harmonic described above is based on two sensors with different measurement principles, and therefore can be combined for a high performance failure detection. The detection is associated with monitoring of statistical properties of the difference between amplitudes of the first harmonic calculated in two different ways.
Notice that load sensors (strain gauges), located on the fixed part of the nacelle, which measure tilt and yaw nacelle moments, could also be used for estimation of the blade loads via inverse Coleman transformation as is described in [10]. The same measurement principle associated with estimated and measured blade loads together with inaccuracies in Coleman transformation (see [22] for details) result in poor detection performance in this case.
The fault detection method proposed in this paper is associated with a change detection, where the change in parameters of distribution of the difference between amplitudes of the first harmonic identifies the fault.
Three types of faults can be detected: (1) abrupt jump faults with a steplike behavior, where the signal changes abruptly from the nominal value to a faulty value; (2) incipient faults with a driftlike behavior, where the signal gradually changes from the nominal value to a faulty value; (3) intermittent faults, where the signal changes from the nominal value to a faulty value, and returns to the nominal value after some time. All three types of faults can be identified via detection of the changes in the parameters of distribution of the difference between the amplitudes.
Notice that this paper does not cover all types of faults which may appear in a wind turbine. Additional fault scenarios together with a benchmark wind turbine model are presented in [23].
Two wind speeds, measured at different heights, see Figure 1, are the input to the turbine load model. The periodic individual blade wind speeds V_{i} , i = 1, 2, 3, defined at the center of each blade in rotating frame can be calculated using these two wind speeds (as boundary conditions), which together with the turbine rotor speed ω_{r} are associated with the individual blade tipspeed ratio λ_{i} as follows:

(1) 
where R is the rotor radius. Individual blade wind speeds have a fundamental frequency associated with a turbine rotational speed. The individual blade wind speed together with the pitch angle and rotor speed define the individual blade flapwise bending moment M_{f,i} , which can be presented as an average value and a periodic component induced by wind shear. This component contains the first harmonic of the turbine rotational frequency and high order harmonics and can be described as follows:

(2) 
where A_{0} = A_{0}(V, ω_{r} , β_{c} ) is the average value of the flapwise bending moment, and A_{h,i} = A_{h,i} (V_{1}, V_{2}, ω_{r}, β_{c} ) is the amplitude of the harmonic h, where d is a total number of harmonics. The average value of the flapwise bending moment depends on the hub wind speed V, turbine speed ω_{r} and collective pitch angle β_{c} . The amplitudes depend on two wind speeds V_{1} and V_{2}, which characterize wind shear, see Figure 1, turbine speed and collective pitch angle.
Notice that wind shear consists of vertical and horizontal components. A vertical component is accounted only in this model as the most pronounced component of the wind shear. A horizontal component, which is not accounted in this model is treated as stochastic variations of the amplitudes.
Discretetime measurements of individual blade flapwise bending moment M_{fk} (where index i is dropped for simplicity) with d harmonics can be written as follows:

(3) 
where _{k} is the harmonic regressor and θ_{*k} is the vector of unknown parameters defined as follows:

(4) 

(5) 
where ω_{rk} is a discretized turbine speed, ξ_{k} is a zero mean white Gaussian measurement noise, k = 1, 2, … is the step number. Notice that the average value of the flapwise bending moment can be easily estimated using a lowpass filter and subtracted from the individual blade flapwise bending moment signals. Therefore, signal (3) contains the periodic component only. Leastsquares estimate:

(6) 
of unknown parameter vector θ_{*k} in a window of a size N, which is moving in time, can be written as follows [19]:

(7) 
and the first harmonic M_{ffk} is recovered according to the relation:

(8) 
where A1k is estimated of the first harmonic.
The matrix is called as information matrix and it is an SDD matrix for a sufficiently large window size [20]. The SDD matrix can easily be inverted making algorithm computationally efficient and implementable. Moreover, estimated parameters can also be rapidly and accurately calculated without matrix inversion [21].
Estimation problem stated above is a challenging estimation problem due to timevarying amplitudes and rotational speed. The accuracy of estimation can be improved via reduction of the window size for fast varying parameters. The properties of the information matrix based on harmonic regressor have not been studied for a sufficiently small window size N. The case, where the window size is not large enough for information matrix to be an SDD matrix is considered in the Appendix. A positive definiteness of the information matrix based on harmonic regressor with four components is shown in this Appendix using the partitioning method.
A measurement noise together with nonstationary nature of estimated parameters and turbine rotational speed are the main obstacles to high performance estimation of the first harmonic of the flapwise bending moment in realtime.
Postprocessing estimation, where all the signals are saved in the buffer and future values of the measured signal are available shows significant improvement with respect to realtime estimation (see Fig. 2).

Figure 2. Reconstruction of the first harmonic from a noisy blade load signal in post processing. Measured signal of the flapwise bending moment, which contains four harmonics is plotted with a blue line. The first harmonic is plotted with a green line, and estimate of this harmonic is plotted with a red line. The flapwise bending moment is presented in normalized unit.

Two wind speeds which are measured at a distance in front of the turbine (see Fig. 1) are used for modelbased estimation of the amplitude of the first harmonic of blade load. The wind speeds which are expected to arrive to the turbine site after some time can be calculated using a classical frozen turbulence assumption [24]. Expected periodic individual blade wind speeds V_{ie} , i = 1, 2, 3 defined at the center of each blade in rotating frame can be calculated using two expected wind speeds. Expected periodic individual blade wind speeds together with turbine speed and pitch angle are the inputs to the lookup table, which estimate individual blade flapwise bending moments as follows:

(9) 
This model is suitable for estimation of first harmonic of turbine rotational frequency only. This harmonic appears in the load in the presence of the wind shear. Amplitude of this harmonic is estimated using least squares method. Notice that two techniques of estimation of the first harmonic described above are based on sensors with different measurement principles and therefore can be combined for high performance failure detection.
Blade load sensor failure detection mechanism is driven by two sensors: (1) LIDAR, which measures two wind speeds at a distance in front of the turbine, and (2) blade strain gauges, which measure a flapwise bending moment. A chart of the detection mechanism is presented in Figure 3. Two amplitudes of the first harmonic of rotational frequency are estimated and compared. A difference between the amplitudes , which is calculated in two ways is normally distributed with a zero mean value for healthy system. The parameters of such distribution can be identified using a sufficiently large number of measured points in the case of a healthy system. The fault detection is associated with a change detection, where the change of the parameters of this distribution identifies the fault.

Figure 3. Detection of blade load sensor faults.

A slow drift of the amplitude of the first harmonic of the strain gauge signal is shown as an example in Figure 4 with corresponding distributions plotted in Figure 5. Such a drift may appear due to the fatigue of the material in the cycle loading or humidity change, which results in a change of gauge factor (see section 'Individual pitch control and unreliable blade load measurements' for all types of faults).

Figure 4. A time chart of the amplitudes of the first harmonic of the healthy system (the first subplot) and a failed system (the second subplot). Amplitudes of the first harmonic calculated using load model and sensor signal are plotted with black and red lines, respectively. Incipient fault associated with a slow drift occurs in step number 700. Histograms of healthy and failed systems are plotted in Figure 5.


Figure 5. Histogram of the difference between the amplitudes of the first harmonic for a healthy system (plotted with red color) and for a failed system, plotted with blue color.

Each sample of the difference between amplitudes should belong to a parent distribution, associated with a healthy system. The fault is detected if the sample is identified as an outlier (outlier is an observation point which is distant from other observations).
An outlier is detected via a two sample ttest, where a hypothesis that a mean value of the distribution that describes a healthy system is equal to the suspected outlier, which is treated as a mean value of a virtual distribution [25], is taken as a null hypothesis. This hypothesis is tested against an alternative hypothesis that the observation point does not belong to the parent distribution that describes the healthy system. The failure is detected if the null hypothesis is rejected in favor of an alternative hypothesis.
A slow drift can also be detected via monitoring of the fluctuations of the mean value of the difference between the amplitudes. This monitoring is performed in a window of a certain size, which is moving in time. Notice that the difference in mean values should be statistically significant for reliable detection of the fault. This can be verified using the two sample ttest, where a significance level represents tradeoff between fastness and detection performance.
The performance of the detection mechanism, described above, depends on the performance of the load model, which is a mean value model that describes the flapwise blade root bending moment (see Fig. 3). The model of the flapwise bending moment is presented in the form of lookup tables (the surfaces in threedimensional space) with the tipspeed ratio and blade pitch angle as input variables. The mean value of the bending moment calculated via lookup tables should coincide with the outputs of the blade load sensors in a uniform wind field for the healthy system. Deviations between the load sensor measurements and output of the model necessitate adaptation of the model, provided that the wind field is uniform. The uniformity of the wind speed across the rotor swept area is detected by comparison of the blade load sensor signals of all three blades and/or by LIDAR measurements. Average values of the load calculated using blade load sensors in a number of working points are memorized for adaptation of lookup tables. The data should be acquired over a relatively large time segment to ensure statistical consistency. Notice that additional requirements may be imposed on input data to avoid erroneous adaptation. If new data are available in a certain operating region only (e.g., at low wind speeds and for small pitch angles), then the part of the surface parameters is adapted (e.g., gradient in pitch angle direction). Adaptation of the lookup table is associated with a motion of the surface in three dimensional space (see Fig. 6). The position and orientation of the surface in threedimensional space change only after adaptation, which in turn allows for a prediction of the bending moment for a wide range of operating variables. This adaptation method was developed first in automotive applications [25] and was successfully applied to turbine model validation with fusion of simulation and measurement data on Big Glenn wind turbine, located outside Gothenburg, Sweden [26].

Figure 6. The flapwise bending moment is presented as a surface with tipspeed ratio and pitch angle as input variables for a certain turbine speed [14]. Mean values of measured flapwise bending moment are plotted with plus signs of a red color. The surface is adapted against measured data, minimizing deviations between the surface and measured points. All variables are presented in normalized units.

Increasing demands on operational reliability, safety, and power output of wind turbines necessitate the development of new highperformance fault detection techniques. Highperformance fault detection is directly associated with predictive maintenance [27], where a component is replaced before the system breaks down, which implies significant savings as well as increases power output. Moreover, high performance fault detection methods are the basis for novel faulttolerant turbine control strategies, where a fault is predictively detected, and the turbine is switched to a safe operation mode to prevent damages, until a maintenance crew arrives at the turbine site. This extends the turbine operation time and increases the power output. Fault tolerance is associated with an ability of the turbine control system to maintain control objectives, despite the occurrence of a fault [2830]. New blade load sensor fault detection technique proposed in this paper is the basis for high performance faulttolerant IPC. This is also a new application for future multibeam LIDAR technology, which is utilized now in the previewbased control only. New detection algorithms can be easily integrated into existing IPC functionality aiming for improvement of uneven load reduction for large turbines.
This work was supported by ÅF. The author is grateful to WindVector AB and Albert Santiago from Avent Lidar Technology for interesting discussions about multibeam LIDAR technologies.
None declared.
Consider the following harmonic regressor:

(A1) 
where ω is a constant rotational frequency and hω with h = 2, 3, 4 are higher harmonics, k = 1, 2, … is the step number. Consider the following information matrix calculated over the window of a size N:

Elements of the matrix G can be evaluated explicitly using the following relation:

(A2) 
where h = 1, …, 4. The first two elements of the first row of this matrix are presented below:


(A3) 
All other elements of this matrix can be evaluated using similar arguments.
Symmetric matrix G can be partitioned as follows:

where


Notice that the diagonal elements of the matrices A and C have periodic part and average part N/2, whereas all other elements of matrices A, B and C have periodic parts only. Matrix G becomes an SDD matrix for a sufficiently large window size N. Suppose that the window size N is large enough for matrices A and C to be SDD matrices, but it is not sufficiently large for matrix G to be an SDD matrix. The matrix G is positive definite if and only if the matrices A and C − B^{T}A^{−1}B are positive definite matrices (see theorem 7.7.6 in [31]). Inverse of the matrix A can be approximated as , when neglecting periodic part in matrix A, where I is the identity matrix. The matrix also has an average part that is proportional to the window size N. This part is canceled by the diagonal elements of matrix D in the term B^{T}A^{−1}B. Therefore there exists a large enough window size N such that the matrix C − B^{T}A^{−}^{1}B is an SDD matrix. Hence the matrix G is a positive definite matrix, since matrix A is also an SDD and positive definite matrix. Notice that the matrix G is not an SDD matrix.
This technique of determination of positive definiteness of information matrix is illustrated via numerical example in the next Section.
Consider the following information matrix G based on the regressor equation (A1) with ω = 0.0228 and the window size N = 83:

This matrix G is not an SDD matrix since the diagonal elements in the second and the third rows are not larger than the sum of the magnitudes of all the other (nondiagonal) entries in the second and third rows. Nevertheless, the matrices and are positive definite and SDD matrices. Therefore the matrix G is a positive definite matrix. Notice that the matrices A and C−B^{T}A^{−1}B are still positive definite (but not SDD) matrices even for a sufficiently small window size N.
Published on 01/06/17
Submitted on 01/06/17
Licence: Other
Are you one of the authors of this document?