Power Maximization and Control of Variable-Speed Wind Turbine System Using Extremum Seeking

Maximizing the power capture is an important issue to the turbines that are installed in low wind speed area. In this paper, we focused on the modeling and control of variable speed wind turbine that is composed of two-mass drive train, a Squirrel Cage Induction Generator (SCIG), and voltage source converter control by Space Vector Pulse Width Modulation (SPVWM). To achieve Maximum Power Point Tracking (MPPT), the reference speed to the generator is searched via Extremum Seeking Control (ESC). ESC was designed for wind turbine region II operation based on dither-modulation scheme. ESC is a model-free method that has the ability to increase the captured power in real time under turbulent wind without any requirement for wind measurements. The controller is designed in two loops. In the outer loop, ESC is used to set a desired reference speed to PI controller to regulate the speed of the generator and extract the maximum electrical power. The inner control loop is based on Indirect Field Orientation Control (IFOC) to decouple the currents. Finally, Particle Swarm Optimization (PSO) is used to obtain the optimal PI parameters. Simulation and control of the system have been accomplished using MATLAB/Simulink 2014.


Introduction
Wind energy is a fast-growing interdisciplinary field that encompasses multiple branches of engineering and science.Wind energy has gained tremendous atten-tion over the past two decades as one of the most promising renewable energies for the future, environmentally friendly solutions, pollution-free and inexhaustible sources [1].
According to speed control, WTs are classified into two types: fixed and variable speed WT.Fixed speed turbines are using SCIG directly connected to grid without using power electronic converters.The speed of the generator is constant regardless of the wind speed.Variable-speed turbines are currently the most used turbines that use different types of generators connected to grid through the power converters.Variable-speed turbines allow the rotational speed to be continuously adapted and controlled in such a manner that the turbine operates constantly at its highest level of aerodynamic efficiency then the maximum energy can be extracted relying upon the wind speed variations.
Maximizing power capture for variable speed Wind Turbine (WT) by means of control theory has been a popular topic for the researchers.According to WT operation regions that will be explained in the next section, the researchers applied different control methodologies by using generator speed control to harvest maximum energy in region II, using blade pitch control to reduce mechanical loads and maintenance cost in region III, or a combination of both.Model-based strategies have been studies for region II control [2] [3] [4].The model-based control approaches can achieve good performance for some ranges of operating conditions but these methods have some limitations and difficulties as the mathematical model of the system is directly incorporated into the controller.Thus, the controller incorporates knowledge of how the system dynamics change with the operating conditions.The Wind Energy Conversion Systems (WECS) are very complex and it's difficult to obtain the dynamics equations otherwise, model-based methods required experimental data and/or lookup tables under specific wind condition and accurate wind speed measurement.
To overcome some difficulties with model-based control methods and reduce the dependency on the mathematical model, other researchers proposed different model-free control methods for Region II operation.Johnson et al. [5] proposed Model Reference Adaptive Control (MRAC), he maximized output power by selecting a suitable control law that related the speed directly with torque then MRAC adapted the torque gain to increase power.Barakati et al. [6] proposed classical MPPT methods, like model-free Perturbation and Observation (P & O) method.The logic of this method is to perturb control input with a certain step size and observe the change in power output.
In this paper, Extremum Seeking Control (ESC) scheme [7] [8] for WT region II operation was proposed to maximize the output power capture without consideration of the load impacts.ESC is a nonmodel-based self-optimizing control strategy that aims to search for unknown input in real-time varying systems by finding the extreme point.The plant excited with some sinusoidal probing signals (dither signal).ESC is applicable in situations where there is nonlinearity in the control problem, and the nonlinearity has local minimum or a maximum.Journal of Power and Energy Engineering Creaby et al. [9] used ESC method to maximize the rotor mechanical power.
They used sinusoidal signal as a dither perturbation signal and second order filters to seek the maximum point.Ghaffari and Krstic et al. [10] used the same method with first order filters.Xiao et al. [11] presented the experimental works to the same model that used by Creaby (CART3 facility).In this work, both schemes have been considered.Moreover, maximization of the electrical power has been worked on in order to produce more accurate results than maximization of the mechanical power.
The approach in this paper focuses on the control of variable speed WT composed of two-mass drive train, SCIG, and voltage source converter control by SPVWM.The controller is designed in two loops, in the outer loop ESC has been used to set a desired reference for PI control of the generator speed and extract the maximum electrical power within high wind speed variation.The inner control loop is based on IFOC for the currents decoupling.The following Sections are organized as a follow.In Section 2, the design of wind energy conversion system has presented with concentration on the dynamic equations of SCIG in dq reference frame and two-mass drive train.In Section 3, the controller design was discussed.In Sections 4 and 5, Simulation results and conclusions are presented respectively.

Aerodynamics Model
WT is composed of different parts to achieve kinetic-to-electric energy conversion.The wind kinetic energy is converted to mechanical energy by the blades.
The mechanical energy is transmitted through the drive train to the generator, which converts mechanical energy into electric energy.This conversion is supported by power converters that deliver the power from the generator to the grid, as shown in Figure 1 [12].
The physical definition for kinetic energy is But the power definition P is the kinetic energy per unit time: where, V is the speed and m  is the air mass that passes the disc of wind tur- bine rotor in a unit length of time, e.g., m AV ρ ∞ =  Thus, the available power in the wind stream is where avail P is the wind power, ρ is the air mass density, is the swept area of turbine rotor, and V ∞ is the wind speed.The mechanical power that can be extracted from the wind power is Journal of Power and Energy Engineering where,

( )
, p C λ β is, the power coefficient, the ratio of the mechanical power to available power, β is the pitch angle of the blade, and λ is the Tip-Speed- Ratio (TSR) can be obtained by where, wtr Ω is the low-speed shaft and t R is turbine rotor radius.
Wind turbine operation regions are typically divided into four primary regions, as seen in Figure 2. Region 1 spans operation from startup to the cut-in wind speed.In Region 2 the turbine is said to be in sub-rated region when wind speeds are between cut-in and rated speed.The primary goal is to capture as much power as possible.In Region 3 wind speeds are high enough to drive the generator at its rated power output; in this case, the goal is to regulate speed and power safely at rated levels.Region 4 occurs when the turbine shuts down to prevent damage due to high wind speeds [13]. where,

Dynamic Equations for Two-Mass Model Drive Train
The model has two masses corresponding to WT rotor and generator, respectively.The dynamic equations of motion for the drive-train are written with respective to the generator mass as given below.

Squirrel Cage Induction Generator (SCIG) Model
The dynamic equations of induction machine in reference frame can be written using either fluxes, currents or both as state-variables.The four loop equations described from the equivalent circuits where the variables ds i , qs i , dr i and qr i represent the state variables and ds V , qs V , dr V and qr V represent the input vector of the state-space model in the dq reference frame.So, the state and input vectors are respectively: T For SCIG, the d and q components of the rotor voltage set to zero, i.e.
( ) since the rotor of squirrel cage machine has slot-embedded bars that are shorted by end rings.The nonlinear fourth-order state model of SCIG can be presented as in the equations below [3] [16].

( )
( ) while, the Equation ( 15) can be rewritten with rotor flux variables instead of rotor current as in Equation ( 18) to be prepared for IFOC.

Voltage Source Converters and Space Vector PWM
Electrical generators are usually controlled by power electronics converters.
Power converters are controlled using a switching strategy.Space vector mod- Sector I, it can be synthesized by 1 V , 2 V and 0 V that gives the equation be- low.The detail explanation is found in [17] and its references, where our implementation is done based on them.

Field Orientation Control (FOC)
FOC can be classified in two design approaches 1) Direct FOC and 2) Indirect FOC, according to the implementation method that is used to estimate the rotor flux r ψ and the synchronous position s θ .The rotor flux orientation is one of the most popular schemes used in AC drives control and WECS control due to its simple and easy implementation.The necessity of using FOC is the decoupled control of the rotor flux and electromagnetic torque of the machine to achieve high dynamic performance.The synchronous frequency equals to the summation of machine rotor frequency r ω and slip frequency sl ω then by the integration of them, the position s θ can be obtained.

( )
For rotor flux orientation, the mathematical model of the SCIG in d-q frame is given in Equations ( 11

Extremum Seeking Control (ESC)
A new method for driving the operating point to the energy maximum, based on Extremum Seeking (ES) was introduced to overcome challenges with the conventional MPPT algorithm and remove the dependence of the plant modeling.ESC is a class of self-optimizing control strategy which can search for unknown time varying input for optimizing a performance index of nonlinear plant.In this Work, the electrical power was used as the performance index and our objective is to maximize this function The input and the output to plant become: ( ) where, ω is the dither frequency equals to 20.6 rad/s.Then Taylor series ex- pansion for Equation (24) becomes: HPF filters out the DC term while passes the AC term then by multiplying with the demodulation signal , that yields: LPF retains DC term to estimate the gradient of the objective function with respect to the control input u while filters out AC term from Equation ( 26).
The loop is closed with an integrator to drive the gradient to zero in steady state.
Several signals were examined as dither signals for ESC.The cosine signal was applied as a harmonic signal and wind speed signal was applied as non-harmonic signal.In addition, the HPF and LPF have been implemented in both as 1 st and 2 nd order transfer functions.
The 1 st and 2 nd order transfer function for HPF is given below with cutoff frequency 3 h ω = rad/s: By feedbacking the electrical power to ES, the optimal speed value is obtain to achieve the maximum energy captured.Different types of dither signals have been examined to get the best reference speed to the system since the selection of the dither signal is an important issue in the application of ESC [19] [20].Figure 5 shows the effect of choosing different type of dither signals to generate the reference speed.The response in blue color (using 1 st order filters and cosine signal as dither signal for ES) suffers from ripples and rises above the optimal speed 157.5 rad/s in the periods (16 -20 s), while the responses in green (using 2 nd order filters and cosine signal as a dither signal for ES) and pink (using 2 nd order filters and wind signal as a dither signal for ES) have high overshoots at 0.5 s.Moreover the green response is decreasing to 156.5 rad/s at 16 s.The response in red color (using 1 st order filters and wind signal as a dither signal for ES) was found to be the best to choose as a reference speed.

Optimized PI Control Using Particle Swarm Optimization
The PI controller is the more commonly used control strategy in electrical machines for speed control.There are many different algorithms to find suitable parameters of the controller.To achieve the closed loop speed control, we firstly select one of the generated speed references from the above sub-section to be the speed command for PI controller.The PI parameters are found using trial and error method to control the speed of generator then Particle Swarm Optimization (PSO) has been presented to determine the optimal PI parameters and ensure the maximum power captured.PSO is an off-line optimization algorithm based on population.PSO algorithm gives high quality solutions and stable Figure 5.The generator reference speed using different types of dither signals.around in a multidimensional search space then during its flight each particle adjusts its position according to its own experience (P-best), and according to the experience of a neighboring particle (G-best), made the best position encountered by itself and its neighbor [21].
The velocity of each agent can be modified by the following equation: The current position can be modified by the following equation: The performance index (ITSE) is chosen as cost function for PSO and given by.
( ) As mentioned before, the closed loop control for regulating the speed of the generator had accomplished using PI controller.The gain values for K p and K i were tuned manually then improved using PSO.The presented data in Figure 6 and Table 2 had been obtained at the following values: iteration = 20, swarm size = 30, inertia = 0.9, and 1 2 2 c c = = .

Simulation Results
The modeling of variable speed wind turbine for mechanical and electrical components has been built based on Equations ( 1)- (17), the currents of the squirrel cage induction generator is decoupled using IFOC based on Equations ( 18)-( 22) and ESC was presented as explained before.MATLAB/Simulink block diagram is shown in Figure 7, also the data of 2 MW wind turbine are found in Appendix A.
The input signal to the whole model is the wind speed which has been modeled and developed based on the Kaimal spectra in [15] [22].Figure 8 shows the obtained signal.
The Simulation has been solved using fixed-step ode4 solver with sampling time ( )    The other parameters values are initially set to zero. Figure 9 and Figure 10 show that the mechanical power had been maximized as result of the approaching of power coefficient from its maximum value ( ) . ESC in the outer loop is based on the generator speed regulation via the electrical power feedback.The ESC searches for the optimal generator's speed to maximize power without accurate knowledge of power map.In the inner loop, PI controller regulates the generator speed in closed loop to attain the optimal speed that is suggested by ESC.As in Figure 6, the speed of the generator was approximately around the optimal speed 157.5 rad/s, meanwhile the speed of turbine rotor in Figure 11 is around 2.1 rad/s since the gearbox ratio equals to 75. Figure 12 shows that the electrical power attains its maximum power capture in some periods of time (2 MW in 6 -8 s, 14 -19 s), meanwhile the electromagnetic torque in Figure 13 approaches to its rated value (14,700 Nm in same periods) according to the decreasing or increasing in wind speed measurements.Finally, Figure 14 shows the three phase stator currents of SCIG.

Conclusion
In this paper, the development of MPPT with ESC to harvest maximum energy captured in the sub-rated speed region has been presented.The simulation results show the strength of ESC to get the position of the operating point neglecting the model dynamics equations as compared with the conventional MPPT methods that highly depend on the dynamic equations and speed measurements.For its being a self-optimizing model-free method, ESC has provided good convergence rate despite the time-varying in wind speed.The output electrical power was used as a performance index to the ESC scheme which shows more interest in energy capture than using the mechanical power.We examined different kinds of dither signals for ESC to achieve slow perturbation to the system and the non-harmonic wind signal with first order transfer function was the better one.The rotor flux orientation control was implemented in inner control loop to achieve the required decoupling between currents and electromagnetic torque.The PI gains values in speed control loop is tuned firstly by trial and error then improved with PSO, where the obtained parameters gains attain the optimal energy captured.
is a nonlinear function, difficult to measure and predict, has different values between each turbine and another depending on turbine parameters and design.The theoretical upper limit for p C is Betz Limit equals to 0.59.The p C is obtained based on the exponential function in[14].
Figure 3. (a) ulation (SVM) is one of the preferred real-time modulation techniques and it is widely used for digital control of voltage source converters.SVM allows reducing commutation losses and the harmonic distortion of output voltage waveform, giving higher amplitude modulation indexes.The SVM control is designed space vector combinations to get a reference voltage ref V .The summation of the voltage multiplied by the time interval of chosen space vectors equal the product of the reference voltage ref V and sampling period z T by knowing in which sector the reference voltage vector is located, we select the two adjacent basic voltage space vectors and one zero-vector.For example, when ref V falls into

dr
)-(18).The decoupling control can be achieved by aligning the stator flux component of the current ds i on the d-axis, and the torque component of the current qs i should align the q-axis, that leads to 0 by a proper selection of dither signal frequency ω , Low Pass Filter (LPF) frequency l ω and High Pass Filter (HPF) frequency h ω , and integrator gain k to set the control û in the real time.ESC is interpreted as single-input setting with dither signal (usually sinusoidal signal ( ) S t ) added to the estimated û .The schematic diagram of whole variable speed WT including control system, mechanical and electrical components, and MPPT method using ESC is shown in Figure 4 [8] [11].

Figure 4 . 8 l
Figure 4. Schematic diagram of variable speed WT control system.

.
The initial values for the speed of generator (high-speed shaft) and the speed of turbine rotor (low-speed shaft) are 156.5, 1.874 rad/s respectively.

Figure 7 .
Figure 7. Simulink block diagram of variable speed wind turbine based on SCIG.

Figure 8 .
Figure 8.The obtained wind speed signal based on kaimal spectra.
speed (m/sec) Journal of Power and Energy Engineering

Figure 14 .
Figure 14.The stator currents of the generator in three phase.

Table 2 .
PI control parameters for generator speed control.