1. Introduction
With the development of wind turbines (WT) in the direction of large-scale, the towers are getting higher and higher, and the blade radius is getting longer and longer, which makes the influence of wind shear and tower shadow effects on the aerodynamic load of the WT more obvious. Wind shear describes the variation in wind speed with vertical elevation, whereas tower shadow reflects the reassignment of wind speed due to the presence of a tower. Due to the factors such as wind shear and tower shadow, the actual wind speed is different everywhere on the rotor disk area, and the difference will change greatly with the WT geometry such as the blade radius and the tower outer diameter, as well as the vegetation and other geography [
1]. Therefore, the wind speed in nature has been changing with time and space, which can be called the spatiotemporal distribution characteristic of wind speed.
As shown in
Figure 1, an anemometer is usually installed near the hub center point O to detect the real-time wind speed, and then the average value over a period of time is taken as the average wind speed (AWS) at the center of the hub. However, the wind speed at this point does not completely reflect the wind speed on the whole rotor disk area. In particular, the rotor dimension of the modern WT is very large, and the wind speed differs across the rotor disk area and changes randomly. Thus, the wind turbine drive torque, the pitch bending moment, and the yaw moment, and so on, will change accordingly. The wind speed spatiotemporal distribution characteristics directly affect the aerodynamic load resulting in load fluctuation of the WT [
2,
3,
4]. Besides, it will also affect the power output of the wind turbine generator system (WTGS) [
1,
5,
6,
7]. The influence of wind shear and tower shadow effects on power in terms of power fluctuation [
1], power loss [
5], and average output power [
6] have been investigated in detail and a frequency domain approach for evaluating the impact of tower shadow and wind shear on the tie-line power oscillations was described in [
7]. They concluded that both wind shear and tower shadow were sources of periodic power fluctuations and average power loss.
From the above researches, it can be inferred that the WTGS fault characteristic may be influenced by wind speed spatiotemporal distribution condition. However, this problem is seldom discussed in previous studies. In [
8], the impact of blade mass imbalance fault on the power characteristics of a doubly-fed induction generator (DFIG) was analyzed considering the effects of wind shear and tower shadow. In [
9], the influence factors causing voltage flicker of WTGS were analyzed, including wind shear, tower shadow, gearbox tooth, and blade break down. And a classifier algorithm which can detect different causes of flicker was proposed. In [
10], the torque and vibration characteristics of the wind turbines drive train during voltage dips were investigated considering the wind shear and tower shadow effects. The other similar researches about this topic are truly seldom reported by now.
Currently, most of the mechanical and electrical fault diagnoses for WT with DFIG are usually based on the average wind speed at the center of the hub. Since the spatiotemporal distribution of wind speed is not considered, it may cause the analytical method or the diagnosis result to deviate from the actual situation, thus affecting the accuracy of the diagnosis. In order to better understand the fault characteristics on the condition of the actual wind speed considering the spatiotemporal distribution, two faults are selected from the common mechanical and electrical faults of the WTGS for comparative analysis. These two faults are the rotor aerodynamic imbalance of WT and the stator winding asymmetry fault of DFIG.
Aerodynamic imbalance means the aerodynamic torques of the three blades are unevenly distributed, which is different from the blade mass imbalance in [
8]. The reason may be a pitch angle difference among the three blades, caused by manufacturing or control errors, or blade airfoil changes due to icing on the blade and other factors. Aerodynamic imbalance can cause the main shaft vibration and further aggravate the fatigue of the blade, bearing, gear, and other parts. At present, the aerodynamic imbalance fault study is mainly focused on two aspects. The first is to analyze the vibration characteristics of the WT to diagnose the aerodynamic imbalance fault [
11,
12], and the second is to study the electrical characteristics under the aerodynamic imbalance fault [
13,
14]. These researches on the aerodynamic imbalance fault are mostly based on the AWS at the hub center, and the spatiotemporal distribution characteristics of the actual wind speed are less considered.
It is reported that about 38% of the WTGS faults are related to the stator according to the failure statistics [
15], and the stator faults is very critical because its important role in the WTGS. To avoid severe damage to stator and the WTGS, the early fault detecting such as the winding asymmetry and resistance variations may have important significance. Because the early fault characteristic is little evident, the diagnostic accuracy becomes very important. Stefani et al. [
16,
17] researched on the stator winding asymmetrical fault based on the frequency analysis of the rotor modulating signals. Dai [
18] used the modified Hilbert–Huang Transform method to analyze the stator current in order to make the fault characteristics more obvious. Williamson et al. [
19] derived simple expressions for the frequencies of the harmonic components in the steady state stator line current of a DFIG operating under various conditions of winding asymmetry. These studies have helped us understand the characteristics of stator winding asymmetry faults. However, they do not consider the spatiotemporal distribution of actual wind speed. Although Gritli et al. [
20] researched the stator winding fault by using wavelet analysis under time-varying conditions, the wind shear and tower shadow were not considered.
A wind speed model named equivalent wind speed (EWS) considering the effects of wind shear and tower shadow is established in our researches [
21,
22] based on the relevant researches [
23,
24,
25]. However, these researches only study the establishment of the wind speed model. The rear drive trains system, the generator and the control system are not taken into account in these researches. To this end, the EWS model is added to the DFIG and its control system, and a complete WTGS model with DFIG is built in this paper. Then the influence of the EWS on the mechanical and electrical fault of the WTGS is studied. Since the motor current signature analysis (MCSA) method is widely used and easier to measure, the DFIG current characteristics of the faults are analyzed mainly in the paper. The main objective of the research is to provide theoretical support for optimizing the fault monitoring system of the WTGS in the future.
2. Equivalent Wind Speed
Considering the effects of wind shear and tower shadow, the EWS of the 3-blade WT can be expressed as the sum of the average wind speed at the hub, the wind shear component and the tower shadow component:
where
Veq is the equivalent wind speed (m/s);
VH is the wind speed at the height of hub center (m/s);
Vws is the wind shear fluctuation component (m/s); and
Vts is the tower shadow fluctuation component (m/s).
According to the former research basis, the expression corresponding to the
Vws and
Vts can be respectively expanded as [
21,
22,
23]:
where
R is the blade radius (m);
r is the distance from the blade element to the hub center (m), as shown in
Figure 1, range from 0 to
R;
H is the hub height (m);
α is the wind shear exponent;
A is the tower radius (m);
x is the distance from the rotor disk plane to the tower center line (m);
β is the azimuth angle (rad);
βb (b = 1,2,3) is the azimuth angle corresponding to the three blades (rad), and
β1 =
β,
β2 =
β1 + 2π/3,
β3 =
β1 + 4π/3.
Substitute the Equations (2) and (3) into (1) and extract
VH from the Equations, and we can get the EWS expressions as follows:
Weq can be called the equivalent wind speed transform coefficient. The EWS model established by Equations (5) and (6) can calculate the wind speed value at any point on the entire rotor disk area. It is an accurate and practical wind speed calculation model. Applying it to the study of the fault diagnosis, it can reflect the wind load conditions of the WT and the fault characteristics.
According to the parameters of a 3-blade 1.5 MW WT, the distribution map of the equivalent wind speed transform coefficient is plotted, as shown in
Figure 2. It can be seen from the figure that due to the effects of wind shear and tower shadow, the wind speed on the rotor disk area is not uniform, and there is obvious periodic fluctuation on the time axis, and the main frequency is three times of the rotor rotating frequency.
In order to facilitate the subsequent analysis, the EWS calculation equation is simplified, and the Fourier fitting is performed on the EWS curve corresponding to the WT whose blade radius is 35 m, and
VH is 11 m/s. The fitting curve and the original EWS curve are shown in
Figure 3.
The fitting curve is highly coincident with the original EWS curve, and the equation corresponding to the fitting curve is:
The equation of the fitting curve contains only the constant term and the trigonometric term, where the constant term
a0 is close to the hub average wind speed
VH = 11 m/s; the
ω in trigonometric term has a value of 9.42, which is close to three times of the rotor rotating angular frequency (3
ωw), so the trigonometric term can be written as cos(3
kωwt), and the more the number of trigonometric terms, the higher the fitting accuracy of the fitted curve and the original curve. In the above figure, when the number of trigonometric terms reaches 8, the fitting accuracy has reached 0.99. The Fourier fitting is also performed on the other EWS curve corresponding to the WT of different parameters, which is basically consistent with the above conclusion, but the number of trigonometric terms should be adjusted according to the actual situation. Therefore,
Veq can be approximated as:
where
Vk and
φk are respectively the amplitudes (m/s) and the phase angles (rad) corresponding to the
k-th trigonometric term; and
ωw is the rotor rotating angular frequency (rad/s).
The Equation (8) is only used to facilitate the following analysis, and the final simulation adopts the Equations (1)–(6) to build the equivalent wind speed model.
According to the wind turbine aerodynamics theory, there is a proportional relationship between the mechanical torque output by the wind turbine and the quadratic wind speed. So, combining the Equation (8), the output mechanical torque
Tm can be obtained as follows:
where
Tm0 is the fundamental component of aerodynamic torque (N·m);
Tm0 = 0.5
ρπ
R3Cp/λ;
ρ is the air density (kg/m
3);
Cp is the optimum power coefficient; and λ is the optimum tip speed ratio.
The third term in Equation (9) is also a polynomial containing cos(3
kωwt) after expansion. Then the expression of the mechanical torque obtained by combining the second and third terms is:
where
Tk and
ϕk are respectively the amplitudes (N·m) and the phase angles (rad) of the torque oscillation components caused by the effects of wind shear and tower shadow in the EWS.
4. Simulation Analysis
In order to verify the correctness of the theoretical analysis described above, a simulation platform of 1.5 MW WTGS with DFIG is built in the MATLAB/Simulink environment (MATLAB R2017a, MathWorks Company, Natick, MA, USA). The sketch of the simulation platform is shown in
Figure 4. The simulation parameters are shown in
Table 1.
The simulation platform mainly includes the EWS model, WT aerodynamics model, gearbox model, DFIG model, and vector control model. Among them, the WT aerodynamics model includes the blade element moment theory, the tower vibration model, and the coordinate transformation model, etc. This aerodynamics model can calculate the aerodynamic torque and other parameters of the WT output under the aerodynamic imbalance caused by the inconsistent pitch angle among three blades, and the detail can be found in our article Reference [
13].
As far as the simulation of stator winding faults, two methods are usually used. One method is to simulate the stator inter-turn short circuit fault by changing the number of shorted turns, and the other is to add an additional resistor or inductor in series with the stator winding to simulate the winding asymmetry fault. Since the second method is simple and easy to implement, it is applied in many situations. For example, References [
16,
17,
18,
20] all adopt the second method. So, in this paper, the stator winding asymmetry fault is simulated by connecting an additional resistor in series with the stator phase A. This method is only used for the simulation of stator winding asymmetry originating from the resistance variations because of unreasonable structure design and electromagnetism thermal field, et al. The stator inter-turn short circuit fault is not considered here.
4.1. Normal Condition Simulation
Firstly, the simulation of normal running condition is performed under two different kinds of wind speed—AWS and EWS. The average constant wind speed at the hub centre is 12 m/s, and the rotor rotating speed in the normal condition is 30 r/min (the corresponding rotating frequency P is 0.5 Hz). The parameters using in the EWS model are: wind shear exponent is 0.4; the distance from the tower middle line to the blade is 4.5 m; tower radius is 1.7 m; blade radius is 35 m; the hub height is 70 m; the wind speed at the hub centre is 12 m/s. The comparison results are shown in
Figure 5,
Figure 6,
Figure 7,
Figure 8.
Figure 5 is the fast Fourier transform (FFT) spectrum of the mechanical torque output from the WT. It can be seen from the figure that there are obvious harmonic fluctuation components in the torque due to the influence of wind shear and tower shadow effects in the EWS, and the harmonic frequencies are 1.5 Hz, 3 Hz…3
kP (
k takes positive number). Besides, the harmonic with frequency of 3P (1.5 Hz) is the main component.
Figure 6 shows the FFT spectrum of the rotor rotating speed, which has the similar result as the torque, and also has distinct harmonic components of 3
kP.
Figure 7 shows the power spectral density (PSD) analysis results of the stator current. On the normal condition with EWS, there exists distinct harmonics in the stator current besides the fundamental wave. The harmonic frequencies in red line around the fundamental frequency are 50 ± 3
kP (1.5
k). That is to say, the modulation frequency is 3
kP which is three times of the rotor rotating frequency. However, these harmonics are not present on the normal condition with AWS as shown in the blue line in
Figure 7.
Figure 8 shows the PSD analysis results of the rotor current. Similar to the analysis of the stator current, the harmonic with modulation frequency of 3
kP appear on both sides of the fundamental frequency (10 Hz) of the rotor current. However, unlike the stator current results, the rotor current harmonic amplitudes are relatively high compared with the fundamental amplitude. In the stator current, the harmonic amplitude at frequency 50 + 3P is about 0.004 p.u., which is about 0.58% of the fundamental amplitude. However, in the rotor current, the harmonic amplitude at frequency 10 + 3P is about 0.0335 p.u., which is about 4.55% of the fundamental amplitude. The reason may be that the rotor side converter mainly controls the stator current, so the harmonic performance in the rotor current is relatively obvious than that of the stator current. The similar effect can be observed in the PSD spectrum. In
Figure 7, the amplitude at the 51.5 Hz is about −24 dB. However, the amplitude at 11.5 Hz is about −14.75 dB in
Figure 8, which is higher than that of the stator current.
4.2. Rotor Aerodynamic Imbalance Simulation
In this part, the rotor aerodynamic imbalance simulation is performed under AWS and EWS respectively. Under each case, three aerodynamic imbalance scenarios are simulated with the pitch angle of one blade adjusted by +1°, +2°, and +3°, respectively, while the other two blades are kept constant. The other running parameters are the same as the former mentioned in normal condition.
Figure 9 and
Figure 10 show the simulation results of the stator current and the rotor current when the pitch angle is adjusted by +3°.
Figure 9 shows the analysis results of the stator current PSD. In
Figure 9, the comparison of four simulation results is given, including the aerodynamic imbalance under the EWS, the aerodynamic imbalance under the AWS, the normal operation under the EWS and the normal operation under the AWS.
It can be seen from
Figure 9 that, on the condition of the aerodynamic imbalance under the AWS (in blue line), the fault harmonic frequencies in the stator current are mainly 49.5 and 50.5 Hz. That is to say, the modulation frequency is mainly the rotor rotating frequency P. The amplitudes at the modulation frequencies of 2P, 3P, 4P, 5P, and 6P are relatively small and the other frequencies are not observed in the figure. However, on the condition of aerodynamic imbalance under the EWS (in red line), the main modulation frequencies in the stator current are 3P (1.5 Hz), 6P (3 Hz) … 3
kP in addition to P. The modulation frequencies such as 2P, 4P, 5P, 7P, 8P, etc. also can be observed, however their amplitudes are relatively small compared with those of the main modulation frequencies. On the normal condition under EWS (in black line), there are only the modulation frequencies of 3
kP, and on the normal condition under AWS (in green line), there is not harmonic frequency except the fundamental frequency of the stator current.
Figure 10 shows a comparison of DFIG rotor current under the four simulation conditions as mentioned above. It can be seen that, on the condition of the aerodynamic imbalance under the AWS, the harmonic frequencies of the rotor current are mainly 9.5 and 10.5 Hz, i.e., the modulation frequency is P. On the condition of aerodynamic imbalance under the EWS, the main modulation frequencies of the rotor current include 3P (1.5 Hz), 6P (3 Hz) ... 3
kP in addition to P. The modulation frequencies of 2P, 4P, 5P, 7P, 8P, etc. can also be observed.
However, unlike the stator current, the main harmonic amplitudes in the rotor current are relatively high. On the condition of aerodynamic imbalance under EWS, the maximum harmonic amplitude (at 51.5 Hz) in the stator current is about 0.0044 p.u., which is about 0.66% of the fundamental amplitude. However, the maximum harmonic amplitude (at 8.5 Hz) in rotor current is about 0.036 p.u., which is about 5.7% of the fundamental amplitude. In the PSD spectrum of
Figure 10, the harmonics amplitudes with the modulation frequencies of P and 3
kP are higher than those of the stator current in
Figure 9. The maximum harmonic at
f + 3P is −23.56 dB in
Figure 9, and the maximum amplitude at
sf − 3P is −16.51 dB in
Figure 10. The comparison data is shown in
Table 2. Consequently, the harmonic performance in rotor current is more obvious than that of the stator current on the aerodynamic fault condition.
Figure 11 shows a comparison of the stator current PSD for three degrees of aerodynamic imbalance under EWS. It can be seen from the figure that the curves under the three imbalance degrees are basically consistent. The main difference is that the amplitude at the frequency
f ± P increases as the imbalance degree increases. However, the amplitudes at the frequency
f ± 3
kP do not change much, mainly because they are caused by the effects of wind shear and tower shadow, and less affected by the aerodynamic imbalance. The rotor current comparison has the similar characteristic as the stator current under the three imbalance degrees.
4.3. Stator Winding Asymmetry Simulation
Through the rotor aerodynamic imbalance fault simulation, the influence of the EWS on the mechanical fault characteristics of the WT is analyzed. Next, in order to analyze the influence on the electrical fault characteristics, the DFIG stator winding asymmetry fault is simulated.
The stator winding asymmetry has been simulated by means of an additional resistor connected in series with one stator phase winding and equal to the rated phase resistance
Rs, and the resistances of the other two phases remain unchanged. Then the simulations under the two kinds of wind speed are respectively performed, which are the stator winding asymmetrical faults under the EWS and the conventional AWS. The simulation results are shown in
Figure 12 and
Figure 13 below.
In
Figure 12, the stator current PSD comparison of four cases is shown, including stator winding asymmetry fault under EWS, stator winding asymmetry fault under AWS, normal operation under EWS and normal operation under AWS.
Figure 12a shows the result of the stator current PSD with frequency in the range of 0–300 Hz. But since the harmonic frequencies are not clear in this figure, the partial enlargement views of stator currents PSD are shown in
Figure 12b,c with frequency ranges of 0–100 Hz and 120–180 Hz respectively.
There are two characteristics in the stator current PSD. Firstly, there is a clear difference between the fault and the normal condition, i.e., the odd-numbered harmonic frequencies of 150 Hz, 250 Hz, etc. exist in the stator current under the winding asymmetry fault, which is consistent with the theoretical analysis as mentioned above. Although only the fault frequencies of 150 Hz and 250 Hz are given in
Figure 12a, there are other odd-numbered harmonic frequencies, which are not given because their amplitudes are relatively small. However, there are not these fault harmonics on the normal conditions. Secondly, there exists distinct difference between EWS and AWS on the condition of stator winding asymmetry fault. In the case of winding asymmetry fault with EWS, there are harmonics with modulation frequency of 3
kP on both sides of the odd-numbered frequencies, while no such harmonics appear around the odd-numbered frequencies under stator winding asymmetry fault with AWS.
Figure 13 shows the rotor current PSD comparison of the four cases as mentioned above.
Figure 13a is the result of the rotor current PSD with the frequency in the range of 0 to 220 Hz.
Figure 13b,c shows the partial enlargement views of the rotor current PSD with the frequency from 0 to 30 Hz and 80 to 120 Hz respectively.
The fault characteristic in the rotor current is basically the same as the stator current. From
Figure 13a, it can be seen that the fault frequencies 90 Hz, 110 Hz, 190 Hz, 210 Hz … (2
k ±
s)
f (
f is 50 Hz, and
s is the slip ratio of −0.2) appear in the rotor current in the case of stator winding asymmetry fault except for the fundamental frequency of 10 Hz. However, there are not these fault frequencies under the normal conditions. In addition, it also can be seen from
Figure 13 that there are harmonics with modulation frequency of 3
kP on both sides of the fundamental frequency and the fault frequencies (2
k ±
s)
f in the case of stator winding asymmetry fault with EWS. However, in the case of stator winding asymmetry with AWS, these modulation harmonics do not appear, and there are only fundamental frequency and the fault frequencies of (2
k ±
s)
f.
5. Discussion
According to the analysis above mentioned, both the mechanical torque and the rotating speed of WT include the periodic fluctuation with the frequencies of 3kP (P is the rotor rotating frequency, k stands for positive integer) due to the effects of wind shear and tower shadow. Two typical mechanical and electrical faults under the EWS and the AWS are mainly studied, and the comparison analysis results are as follows:
Comparison results of the rotor aerodynamic imbalance simulation: in the case of AWS, the fault modulation frequency in stator and rotor currents is mainly P; and in the case of EWS, the fault modulation frequencies in stator and rotor currents are mainly P and 3kP. In addition, there exists frequency modulation between P and 3kP. The harmonic performance in rotor current is more obvious than that of the stator current on the aerodynamic fault condition.
Comparison results of the stator winding asymmetrical fault simulation: in the case of AWS, the fault frequency in the stator current is mainly k1f (k1 stands for odd number, f is the grid frequency), and the fault frequency in the rotor current is (2k2 ± s)f (k2 stands for positive integer); while, in the case of EWS, the fault frequencies in the stator current include k1f ± 3kP and k1f, and the fault frequencies in the rotor current include (2k2 ± s)f ± 3kP and (2k2 ± s)f. Similarly, the harmonic performance in rotor current is more obvious than that of the stator current.