Next Article in Journal
Pool-Based Genetic Programming Using Evospace, Local Search and Bloat Control
Next Article in Special Issue
Adaptive Control of Resistance Spot Welding Based on a Dynamic Resistance Model
Previous Article in Journal
Comparison of Splitting Methods for Deterministic/Stochastic Gross–Pitaevskii Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Direct Power Control Optimization for Doubly Fed Induction Generator Based Wind Turbine Systems

by
Mohammed Mazen Alhato
1,* and
Soufiene Bouallègue
1,2,*
1
Research Laboratory in Automatic Control (LARA), National Engineering School of Tunis (ENIT), University of Tunis, El Manar, BP 37, Le Belvédère, 1002 Tunis, Tunisia
2
High Institute of Industrial Systems of Gabès, 6011 Gabès, Tunisia
*
Authors to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(3), 77; https://doi.org/10.3390/mca24030077
Submission received: 16 July 2019 / Revised: 18 August 2019 / Accepted: 22 August 2019 / Published: 26 August 2019

Abstract

:
This study presents an intelligent metaheuristics-based design procedure for the Proportional-Integral (PI) controllers tuning in the direct power control scheme for 1.5 MW Doubly Fed Induction Generator (DFIG) based Wind Turbine (WT) systems. The PI controllers’ gains tuning is formulated as a constrained optimization problem under nonlinear and non-smooth operational constraints. Such a formulated tuning problem is efficiently solved by means of the proposed Thermal Exchange Optimization (TEO) algorithm. To evaluate the effectiveness of the introduced TEO metaheuristic, an empirical comparison study with the homologous particle swarm optimization, genetic algorithm, harmony search algorithm, water cycle algorithm, and grasshopper optimization algorithm is achieved. The proposed TEO algorithm is ensured to perform several desired operational characteristics of DFIG for the active/reactive power and DC-link voltage simultaneously. This is performed by solving a multi-objective function optimization problem through a weighted-sum approach. The proposed control strategy is investigated in MATLAB/environment and the results proved the capabilities of the proposed control system in tracking and control under different scenarios. Moreover, a statistical analysis using non-parametric Friedman and Bonferroni–Dunn’s tests demonstrates that the TEO algorithm gives very competitive results in solving global optimization problems in comparison to the other reported metaheuristic algorithms.

1. Introduction

Recently, the increasing consumption of electrical energy, depletion of fossil fuels and the environmental problems related to using the non-renewable sources have promoted a growing interest in renewable energies [1,2]. Wind power as a renewable source represents an important and a promising solution for electrical demand and has recently gained more attention for its significant merits of cleanness and resource abundances. Various wind energy configurations are produced from the intensive studies and researches that are carried out in wind systems. One of the most common configurations is the grid-connected Doubly Fed Induction Generators (DFIGs) equipped with variable speed Wind Turbines (WTs). This configuration was widely installed in wind industry due to its significant merits such as independent control of active and reactive powers, low converters’ costs and mechanical stress reduction [2].
A DFIG is a wound rotor generator in which the stator windings are directly connected to the grid. The rotor windings are connected to the grid through a back-to-back power converter that is composed from the Rotor Side Converter (RSC) and Grid Side Converter (GSC) components. This configuration allows the converter to handle the fraction of 20% to 30% of the total power. Therefore, the losses in the converter are lower compared to a system where the converter has to handle the full power leading to an improvement in the total efficiency [2,3].
However, the switching of Insulated-Gate Bipolar Transistors (IGBTs) at both the RSC and GSC components is usually performed by Sinusoidal Pulse Width Modulation (SPWM) strategy, which generates the high frequency content in the utility current. A simple L-filter or an LCL-filter is adopted to reduce the Total Harmonic Distortion (THD) of the utility current [4]. To meet the grid code requirements, the LCL-filter is predominant in reducing the utility current harmonics. Indeed, it can lead to a better attenuation of harmonics using small values of inductances. However, the resonance phenomena of the LCL-filter must be damped properly in order to prevent the possible instability of these systems [4]. One way for solving this problem is employing a passive damping circuit. The passive damping is built by inserting a resistor branch in series or in parallel with the inductors or capacitors of the LCL-filter [4].
The conventional control scheme of the grid connected DFIG wind turbine system is built based on a vector control method [5,6,7,8]. Here, the Stator Flux Orientation (SFO) scheme is selected for DFIG power regulation at the RSC. The SFO-based vector control strategy enables a decoupled regulation of the active and reactive powers that is flowing between the DFIG and the grid. The active and reactive powers’ control strategy is performed by regulating the converter currents using Proportional-Integral (PI) controllers [5,6,7,8]. On the other hand, the Voltage Oriented Control (VOC) is adopted at the GSC, which the DC-link voltage is maintained to constant value [9].
Zhou and Blaabjerg [5] proposed a frequency-domain based approach to attain the PI controllers’ gains for the inner and outer loops at both RSC and GSC circuits. The gains of the PI controllers that verify the design condition are selected for control loop. However, Hamane et al. [6] performed the direct power control of the DFIG using both PI and sliding mode controllers. The synthesis of the PI controller is based on an algebraic pole compensation method. Moreover, Hamane et al. [7] developed a comparative study of PI, RST, sliding mode and fuzzy supervisory controllers, and the synthesis of PI and polynomial RST controllers are tuned based on the pole compensation and pole assignment methods, respectively. In addition, the indirect vector control strategy is adopted, in which the matrix converter is used at the wound rotor of the DFIG instead of the conventional back-to-back converter in [8]. The active/reactive powers and current components are regulated by using the PI controllers that are designed using the pole placement technique. Although these methods have presented good performances, the main drawback for this type of control is that the performance of the DFIG system highly depends on a proper tuning of the PI controller gains. However, since the parameters of PI controllers depend on the precise mathematical representation, the control schemes become prone to error. In addition, the execution of tuning process can be time-consuming, and the optimal gains may not be obtained. Hence, proposing a systematic approach to find the best parameters setting of PI controllers is an interesting task and the metaheuristics-based hard optimization theory may be considered a promising solution.
In recent years, an enormous variety of metaheuristics optimization algorithms has been applied to solve complex and hard problems in science and engineering fields. Bekakra and Attous [10] presented the Particle Swarm Optimization (PSO) algorithm to obtain the optimal gains of PI controllers for the indirect control of the active and reactive powers’ loops of a DFIG. Integral Absolute Error (IAE), Integral Time-weighted Absolute Error (ITAE) and Integral Square Error (ISE) performance criteria were selected as objective functions. Vieira et al. [11] applied the Genetic Algorithm (GA) to obtain the gains of the PI controllers at the RSC, where the active and reactive powers loops were directly controlled. In addition, an indirect power control strategy was achieved using GA-tuned integral sliding mode controllers, where the GA was used to tune the parameters of the direct and quadrature rotor current controllers [12]. Moreover, Assareh et al. [13] proposed a hybrid GA along with a gravitational search algorithm to attain the optimal gains of the PI controller for the torque regulation of a DFIG-based WT system. In all of these studies, the optimization algorithms are applied to obtain only the PI controllers’ gains of the active and reactive powers at the RSC without discussion about the DC-link voltage regulation loop at the GSC and their effect on the overall control system.
Therefore, this work tends to apply a unified Thermal Exchange Optimization (TEO) metaheuristic algorithm to optimize the gains of PI controllers for the outer-loops in the classical vector control scheme of a DFIG-based wind energy system. In particular, this paper deals with the PI controllers tuning of the active and reactive power control loops in the RSC. It also treats the PI controller optimization-based design for the DC-link voltage loop in the GSC component.While the inner current loops at both RSC and GSC are tuned according to classical methods. The global TEO metaheuristic that has been proposed by Kaveh and Dadras in 2017 [14] is mainly adopted to tune the controllers’ parameters of the modeled DFIG-based Wind Energy Converter System (WECS). Moreover, other well-known methods such as PSO, GA, Harmony Search Algorithm (HSA), Water Cycle Algorithm (WCA) and Grasshopper Optimization Algorithm (GOA) are used for an empirical comparison study. Moreover, a statistical analysis is performed to check the significance of each algorithm by using nonparametric tests such as Friedman’s rank and Bonferroni–Dunn’s test. The main contribution of this work is to provide a systematic and less complex procedure based on an advanced TEO algorithm in order to design and tune all outer-loops PI controllers for the well-known vector control scheme of a DFIG-based wind energy system. The classical trials-errors based methods of PI controller tuning are no longer used and the design time is further reduced. The drawbacks of the classical tuning methods are significantly reduced.
The remainder of this paper is organized as follows. The mathematical model of DFIG based WT is presented in Section 2. The description of the vector control scheme for the DFIG system at both RSC and GSC is analyzed in Section 3. In addition, Section 4 is devoted to the formulation of the outer-loops PI controllers’ tuning problem given as a constrained optimization problem to be solved by the proposed TEO algorithm. In Section 5, such a TEO algorithm is described and its pseudo-code for the software implementation is given. Section 6 presents the implementation and the validation of the proposed TEO-tuned PI controller approach. Concluding remarks are given in Section 7.

2. Modelling of the DFIG Based Wind Energy Converter

The configuration of the DFIG-based WECS is depicted in Figure 1. A WT is joined to the DFIG by means of a gearbox. The DFIG is an induction machine, which the stator windings are directly connected to the grid, while the rotor windings are connected to the grid thanks to a back-to-back converter.

2.1. Modelling of the Wind Turbine

The rotor blade of the WT is responsible for catching the wind power and converting it into kinetic energy. The captured mechanical power P m is given as follows [2]:
P m = 1 2 ρ C p ( λ , β ) π R 2 V w 3
where ρ is the air density, C p is the power coefficient and it depends on both the tip-speed ratio λ and the blade pitch angle β , R is the turbine radius and V w is the wind speed.
The tip speed ratio is defined as the ratio of the blade tip speed to the wind speed. It is given by:
λ = Ω t R V w
where Ω t is the angular speed of the WT.

2.2. Modelling of the DFIG

The dynamical model of the studied DFIG can be expressed by the following electrical equations [6,7,8]:
{ V d s = R s i d s + d φ d s d t ω s φ q s V q s = R s i q s + d φ q s d t + ω s φ d s V d r = R r i d r + d φ d r d t ( ω s ω m ) φ q r V q r = R r i q r + d φ q r d t + ( ω s ω m ) φ d r
where V s and i s are the stator voltage and current, V r and i r are the rotor voltage and current, φ s and φ r are the stator and rotor flux linkages, R s and R r are the stator and rotor resistances, ω s and ω m are the stator and rotor angular frequencies, respectively. The subscripts d and q denote the direct and quadrature axis components, respectively. The stator and rotor flux linkages are defined as follows:
{ φ d s = L s i d s + L m i d r   and   φ q s = L s i q s + L m i q r φ d r = L r i d r + L m i d s   and   φ q r = L r i q r + L m i q s
where L s , L r and L m are the stator, rotor and magnetizing inductances, respectively.

2.3. Modelling of the GSC and the DC-Link

The GSC is connected to the grid through an LCL-filter. However, for better understanding the control of GSC, it is necessary to describe the model of the grid side system and DC-link parts. The mathematical formulation of the GSC in the dq synchronous frame is defined by Equation (5). In this equation, L T represents the sum of the converter L i and grid side L g inductances. In fact, the LCL-filter can be approximated to an L-filter equal to the sum of the LCL-filter inductors [4]:
{ L T d i g d d t = R T i g d + ω g L T i g q + e g d V d f L T d i g q d t = R T i g q ω g L T i g d + e g q V q f
where V d f and V q f are the components of the converter side voltage, i g d and i g q are the components of the grid currents, e g d and e g q are the components of the grid voltage, and R T is the sum of the converter and grid side resistors.
The grid active and reactive powers are expressed as [9]:
{ P g = 3 2 ( e d g i d g + e q g i q g ) Q g = 3 2 ( e q g i d g e d g i q g )

3. Vector Control of the DFIG-Based Wind Energy Converter

The classical direct power control of the DFIG system is divided into the RSC and GSC control loops, where the control structure for the RSC and GSC components consist of two cascaded control layers. The inner PI layers are adopted to regulate the components of rotor and grid currents. On the other hand, the outer PI ones in the RSC are implemented to control the active and reactive powers while the outer PI loop in the GSC is adopted to maintain the DC-link voltage to its reference value. However, in this work, the inner PI loops in the RSC and GSC circuits are designed based on a pole assignment technique [8,15], while the outer PI ones are tuned based on the proposed metaheuristics algorithms-based approach.

3.1. Control of the RSC

The RSC circuit is responsible to independently regulate the active and reactive powers in order to extract the maximum available power. The vector control by the SFO scheme is employed for the power regulation of the DFIG system.
Considering that the electrical network is stable, the stator flux φ s is constant. In addition, for the medium and high power DFIG, the value of the stator resistance is very small and therefore it can be neglected [6,7]. Hence, based on these assumptions, the stator voltages and fluxes in Equation (3) and Equation (4) can be rewritten respectively as follows:
{ V d s = 0   V q s = V s = ω s φ d s
{ φ d s = φ d s = L s i d s + L m i d r   φ q s = 0 = L r i q s + L m i q r
The active and reactive stator powers and rotor voltages are given, respectively, by:
{ P s = 3 2 ( V d s i d s + V q s i q s ) = 3 2 L m L s V s i q r Q s = 3 2 ( V q s i d s + V d s i q s ) = 3 2 ( V s 2 ω s L s L m L s V s i d r )
{ V d r = R r i d r + σ L r d i d r d t σ ω r L r i q r Δ d r V q r = R r i q r + σ L r d i q r d t + σ ω r L r i d r + L m L s ω r φ d s Δ q r
where σ = L r L m 2 / L s is the machine leakage coefficient and ω r = ω s ω m is the slip angular frequency.
As shown in Equation (9), the stator active and reactive powers are independent from each other. In addition, the components of stator power are linearly varying with the direct and quadrature rotor currents. This powers’ regulation is thus performed using PI controllers for the d and q axis components of the rotor currents.
The outer control loops for the stator active and reactive powers are tuned based on the proposed metaheuristics-based procedure. Hence, the design of PI controllers for the inner rotor current loops is detailed according to Figure 2.
Based on a pole placement method as follows:
R P I ( s ) = K p i r ( 1 + 1 T i i r s )
where K p i r and T i i r are the proportional gain and time constant of the rotor PI controllers, respectively.
The transfer function that describes the relationship between the voltages and currents dynamics is given as follows [8,15]:
R plant ( s ) = i r d q V r d q = 1 R r + σ L r s
Since the same transfer functions for the direct and quadrature rotor currents are considered, the closed-loop transfer function R c l i ( s ) between the reference rotor current i r d q and the actual one i r d q is written as:
R c l i ( s ) = i r d q i r d q = K p i r ( T i i r s + 1 ) T i i r σ L r s 2 + T i i r ( R r + K p i r ) s + K p i r
By equating the denominator of the closed-loop transfer function of Equation (13) under the general form of a second order system, the parameters of PI rotor current controllers can be found as follows [8,15]:
{ K p i r = 2 ξ ω n σ L r R r T i i r = 2 ξ ω n R r σ ω n 2 L r
where ξ and ω n are the damping coefficient and natural frequency of the desired closed-loop reference model, respectively.
For the design purpose, the discussion is made with respect to the choice of the damping coefficient ξ and natural frequency ω n . Since the inner current loop in the cascade control scheme will have a much larger bandwidth than the one used in the outer-loop, a regulation at 1 / 20 of the switching frequency f s w , R S C is retained. For the outer loop, a regulation between 1 / 20 and 1 / 10 of the inner loop bandwidth is made [5]. Then, feed forward compensations Δ d r and Δ q r should be added back to generate the desired rotor voltages V d r and V q r .

3.2. Control of the GSC

The VOC strategy is applied to control the GSC, which usually contains one outer PI control loop that regulates the DC-link voltage regardless of the magnitude and the direction of the rotor power. Two inner PI current loops that regulate the direct and quadrature gird currents are also included [16,17]. In addition, the passive damping strategy is employed to mitigate the resonance problem. The mathematical representation of GSC is built on the approximated model of the LCL-filter, which is used to tune the PI grid current controllers. Whereas, the whole transfer function of the LCL-filter with passive damping is taken into account for the stability analysis and investigation [4]. The block diagram of the grid current loops is shown in Figure 3.
To implement the VOC scheme, the d-axis is aligned with the grid vector voltage. Therefore, this leads to a d-axis grid voltage equal to its magnitude, and the q-axis voltage which is then equal to zero. Hence, the grid power expressions can be expressed as [16,17]:
{ P g = 3 2 e d g i d g Q g = 3 2 e d g i q g
Equation (15) decides the active power and consequently the DC-link voltage is controlled via the direct current whereas the quadrature current is used to regulate the reactive power that exchanges with the grid [17]. With the VOC approach, the dynamic equations of the grid currents are rewritten as follows:
e d g = R T i d g ω g L T i q g Δ d g + L T d i d g d t + V d f 0 = R T i q g + ω g L T i d g Δ q g + L T d i q g d t + V q f
By applying the same methodology of PI rotor current controllers design in the RSC, the gains of PI controllers for the grid currents dynamics are found as [15,18]:
{ K p i g = 2 ξ ω n L T R T K i i g = ω n 2 L T
where ω n = 847 . 80   rad / sec and ξ = 0.707 .

4. PI Controllers Tuning Problem Formulation

In the PI control framework, appropriate values of K p and K i gains are generally obtained by empirical methods and trials-errors based procedures [19]. These nonsystematic and challenging tasks become more difficult and time-consuming, especially for the complex and large-scale systems like the studied DFIG-based WECS. So, the idea of formulating the K p and K i gains’ selection as an optimization problem is a promising solution. Such a control problem can be nonlinear, non-smooth or even non-convex and can be effectively solved thanks to advanced metaheuristics [10,11]. In this work, three PI controllers for the outer-loops at both RSC and GSC components are considered for the optimization process. These PI controllers for the DC-link voltage, active and reactive powers’ dynamics are systematically tuned thanks to the proposed TEO metaheuristic which is described in the following. Figure 4 gives the proposed optimization-based tuning scheme of the PI controllers for the DFIG-based WECS.
The decision variables of the formulated optimization problem are obviously the proportional and integral gains of PI controllers for the active/reactive powers and DC-link loops. They are specified for the optimization process as follows:
x = [ K p P s , K i P s , K p Q s , K i Q s , K p d c , K i d c ] T   S R + 6
where S = { x R + 6 , x l o w x x u p } denotes the bounded search space for all PI parameters.
However, the defined objective functions are minimized taking into account a scope of time-domain restrictions. These are concerning to the maximum overshoot δ max , steady-state error E s s , rise time t r and/or settling time t s of the closed-loop system step-response [20]. So, the tuning issue associating the PI controllers of the DFGI-based WECS can be expressed as follows:
{ Minimize   f m ( x ) , m { I A E , I S E , I T S E , I T A E } x = [ K p P s , K i P s , K p Q s , K i Q s , K p d c , K i d c ] T S R + 6 subject   to : g 1 ( x ) = δ P s δ P s max 0 g 2 ( x ) = δ Q s δ Q s max 0 g 3 ( x ) = δ d c δ d c max 0 K p j , min K p K p j , max K i j , min K i K i j , max , j { P s , Q s , d c }
where f m : R + 6 R are the cost functions, g q : R + 6 R are the problem’s inequality constraints, δ d c , δ P s and δ Q s are the overshoots of the controlled DC-voltage, active and reactive powers dynamics, respectively, δ d c max , δ P s max and δ Q s max denote their maximum given value.
Since the optimization problem (19) is a multi-objective type, weighted sum-based aggregation functions are used as follows to handle with the multiple costs [21]:
f I A E ( x ) = w P s 0 T | e P s ( x , t ) | d t + w Q s 0 T | e Q s ( x , t ) | d t + w d c 0 T | e d c ( x , t ) | d t
f I S E ( x ) = w P s 0 T e P s 2 ( x , t ) d t + w Q s 0 T e Q s 2 ( x , t ) d t + w d c 0 T e d c 2 ( x , t ) d t
f I T S E ( x ) = w P s 0 T t e P s 2 ( x , t ) d t + w Q s 0 T t e Q s 2 ( x , t ) d t + w d c 0 T t e d c 2 ( x , t ) d t
f I T A E ( x ) = w P s 0 T t | e P s ( x , t ) | d t + w Q s 0 T t | e Q s ( x , t ) | d t + w d c 0 T t | e d c ( x , t ) | d t
where T denotes the total simulation time, w P s > 0 , w Q s > 0 and w d c > 0 are the weighting coefficients of the aggregation functions satisfying w P s + w Q s + w d c = 1 , and e j ( . ) , j { P s , Q s , d c } denote the tracking errors between the plant output and the relative set-point values, i.e., e P s ( x , t ) = P s * P s ( x , t ) , e Q s ( x , t ) = Q s * Q s ( x , t ) and e d c ( x , t ) = V d c V d c ( x , t ) .
The judgment matrix method is used to determine the weighting coefficients of the objective functions [22,23]. Such a method grades all objective functions based on the importance of each one. After calculating the eigenvalues of such a matrix, the weighting coefficients of functions (20)–(23) can be chosen as w P s = 0 . 6370 , w Q s = 0 . 2583 and w d c = 0.1047 .

5. Thermal Exchange Optimization Algorithm

The Thermal Exchange Optimization (TEO) algorithm is a novel metaheuristic inspired by Newton’s law of the cooling [3]. In this population-based metaheuristic, each agent is modeled as a cooling object. By associating another agent as a surrounding fluid, a heat transfer and thermal exchange occurs between them. The new temperature of each agent is considered as a new position of the potential solution in the search space.
In a d-dimensional search space and at the k t h iteration time, each object of the population is described by its temperature T i k = ( T i , 1 k , T i , 2 k , , T i , d k ) , ( i , k ) [ 1 , N p o p ] × [ 0 , N i t e r ] . In order to enhance the optimization performances, the TEO algorithm uses a Thermal Memory (TM) to store a number of historically best vectors as well as their related fitness. Therefore, these stored solution vectors are added to the population and the same numbers of current worst objects are deleted. After that, a growing order of solutions is retained according to their cost function values. The N p o p sorted objects are equally divided into two groups of environment and cooling objects as shown in Figure 5. The environment objects are denoted as T 1 k , T 2 k ,   ,   T N p o p / 2 k while the cooling ones are T N p o p / 2 + 1 k , T N p o p / 2 + 2 k , ,   T N p o p k .
Since the object has lower η can exchange the temperature slightly, the TEO metaheuristic proposed a similar formula to evaluate the value of η for each object as described in Equation (24):
η = Cost   ( object ) Cost ( worst   object )
Over the course of iterations, the value t = k / N i t e r increases as follows leading to improving the exploration mechanism [3]. In order to enhance the TEO exploration capacity, Equation (25) has been proposed to avoid the trapping in local optima and update the environmental temperature as follows:
T e n v , i k + 1 = ( 1 U { 0 , 1 } ( c 1 + c 2 ( 1 t ) ) ) T e n v , i k
where c 1 and c 2 are the controlling variables chosen as 0 or 1, U { 0 , 1 } is a uniformly random number, T e n v , i k is the previous environmental temperature, ( 1 t ) is proposed to reduce the randomness by closing to the last iterations. When the decreasing of the randomness is not required, the value of c 2 is set to zero.
The new temperatures for each object, i.e., either cooling objects or environmental ones, are updated as follows:
T i k + 1 = T e n v , i k + ( T i k T e n v , i k ) e η t
In the TEO formalism, another formulation is suggested to improve the exploration ability. A control parameter 0 p r o 1 is introduced and determines whether a component of each cooling object must be changed or not. For each object, such a parameter is compared to R a n d ( i ) , ( i = 1 , 2 , , N p o p ) , which is a random number uniformly distributed between 0 and 1. If R a n d ( i ) < p r o , one dimension of the i t h object is chosen randomly and its value is regenerated as follows [3,24]:
T ( i , j ) = T j , min + U ( 0 , 1 ) × ( T j , max T j , min )
where T ( i , j ) is the j t h component of the i t h object, T j , max and T j , min are the upper and lower bounds of the j t h component, respectively.
Finally, the steps of the proposed TEO algorithm are summarized as follows:
  • Step 1. Randomly initialize the temperature for all objects T i 0 , i = 1 , 2 , , N p o p .
  • Step 2. Calculate the fitness of each search object.
  • Step 3. Save some T best vectors and their related cost values in the TM.
  • Step 4. Add the saved solutions and remove the same numbers of the worst objects.
  • Step 5. Arrange the objects according to their related fitness in an ascending order.
  • Step 6. Divide the objects into two equal groups: environment and cooling objects.
  • Step 7. Calculate the parameters η and t .
  • Step 8. Change the environment temperatures by Equation (25).
  • Step 9. Update the temperatures according to Equations (26) and (27).
  • Step 10. Check the termination criterion and repeat the iterations.
To describe the proposed metaheuristics-based tuning strategy of the PI controllers in the studied DFIG-based energy conversion system, a detailed flowchart is shown in Figure 6.

6. Simulation Results and Discussions

6.1. Execution of the Metaheuristic Algorithms

The proposed metaheuristics-based direct power control of the studied DFIG-based WECS is built using the MATLAB/Simulink environment. The respective nominal parameters of the DFIG, power converters at both sides, L-filter and LCL-filter are presented in Table 1. Specifically, this work deals with the PI controllers tuning for the active and reactive powers loops in the RSC circuit, as well as it treats with the PI controller optimization-based selection of the DC-link voltage loop in the GSC component. It is worth indicating that the inner current loops in the RSC and GSC are tuned according to the pole assignment method. Since the grid current closed-loop of the LCL-filter is unstable, the passive damping method is applied as a solution to ensure the stability of such current dynamics. To investigate the effectiveness and superiority of the proposed TEO method, the well-known global algorithms GA, PSO, HSA, WCA and GOA are considered for comparison purposes.
All reported algorithms are independently run 10 times. The termination criterion is set as a maximum number of iteration reached N i t e r = 100 for a population size of the thermal agents equal to N p o p = 50 . All the values of the used common parameters are kept equal. All algorithms have been executed on a PC computer with Core TM i5-7200U CPU and 2.5 GHz/8.00 GB RAM. The specific control parameters of the reported algorithms are listed in Table 2.
The optimization problem (19) is minimized for various performance criteria, i.e., IAE, ISE, ITAE and Integral Time-weighted Square Error (ITSE), and under time-domain operational constraints. These performance indexes are calculated within two conditions, which are considered as follows:
-
14.3% step change in the reference of DC-link voltage at time t = 0.5 sec ;
-
step change in the reference stator reactive power at time t = 0.8 sec .
To find the accurate values of the decision variables x = [ K p P s , K i P s , K p Q s , K i Q s , K p d c , K i d c ] , the bounded search domain is initially set to 0 x R + 6 10 5 . Several runs of the algorithm are performed within this initial limitation. The results indicated that the search space could be reduced to 0 x R + 6 400 . Therefore, the TEO algorithm explores within a smaller search domain for the next runs and more precise solution output vector is achievable.
Table 3 gives the statistical results attained by the introduced algorithms under minimizing the cost functions described in Equations (20)–(23). The worst, average, best values and the standard deviation (STD) experimental results are summarized in the Table 3, where the optimal mean value of each index is highlighted in bold and underlined. In addition, the Elapsed Time (ET) is also listed, which is defined as the time that any algorithm requires finding the best solution. It can be clearly observed that the proposed TEO produces very competitive solutions with the reported algorithms. Appendix A summarizes the obtained gains of the PI controllers for each proposed optimization method. Indeed, these tuned PI controllers’ gains lead to the best transient and steady state responses of the entire reported algorithms.
In addition, Figure 7 shows the convergence histories of the mean objective function values of IAE, ISE, ITAE and ITSE criteria, respectively. Indeed, it is shown that the proposed TEO metaheuristic for both IAE and ISE indices outperforms the other reported methods in terms of the fastness and non-premature convergence as well as the solutions’ quality. The TEO-based method for ITAE and ITSE criteria gives the best solution as a second order after the GA one. From these results, the superiority of the TEO metaheuristic is still shown in terms of exploitation and exploration capabilities for local and global searches. This further justifies the use of such a global metaheuristic to systematic and easy design of the direct power control of the DFIG-based WECS. The Box-and-Whisker plots of the mean objective function values are shown in Figure 8. From these results, one can observe that the boxplots of the TEO algorithm for all given performance criteria are generally lower and narrower than other algorithms. This confirms the high reproductivity of the TEO algorithm toward finding the optimal values of the solution.
In the remainder of the results, the obtained gains of the optimized PI controllers are used to assess the frequency-and time-domains performances of the outer-loops at the RSC and GSC power components. Moreover, comparisons with the classical pole placement [8,15], frequency response [5] symmetrical optimum [27], Ziegler–Nichols [28] and Tyreus–Luyben [29] tuning methods are made for the DC-link voltage dynamics as shown in Table 4 and Figure 9. The closed-loop response of the DC-link voltage dynamics is investigated under variable voltage profile. The DC-link voltage reference is varied up and down at various step levels. Figure 9 describes such a response around the final set-point value for different tuned PI controllers. Although, in real systems, the DC-link voltage is required to be constant. This scenario is proposed to check the capability of the proposed TEO algorithm under different circumstances. The aim is to show the difference between the classical tuning methods and the proposed optimization-based one. Referring to this result, the ITAE-based TEO algorithm can regulate the DC-link voltage dynamics with higher performance compared to the other algorithms.
The voltage tracking is achieved with best performance in terms of response precision and fastness, i.e., the steady-state error and rise/settling time metrics are minimal as depicted in Table 4. The performance in terms of response damping is acceptable for the proposed TEO-tuned PI control as shown in Figure 9. Referring to the numerical results in Table 4, one can clearly observe that the proposed TEO algorithm showed the best rise/settling times and steady-state error except the overshoot index, which the Tyreus–Luyben method gained the minimum overshoot. Since several runs were executed to obtain the PI controllers’ gains with the classical Ziegler–Nichols and Tyreus–Luyben methods, the tuning process becomes a tedious and time-consuming task. These minor degradations of the overshoot performance do not influence the effectiveness of the proposed TEO-based tuning method with respect to the systematization of the synthesis procedure, the simplicity of implementation and the superiority in other performance indices. The proposed TEO algorithm found the optimal gains of the PI controllers within a reduced computation time and under operational constraints compared to the reported trials-errors based tuning procedures which generally give local solutions for the formulated control problem.
Regarding to the grid current dynamics, the same tuning methods with adding the technical optimum instead of the symmetrical ones are used to select and adjust the PI controllers’ gains. All obtained results are reported in Table 5.
The control performance of the introduced tuning methods is investigated by supposing that there is a step change in the q-axis grid current reference at time t = 0.2 sec . The reference q-axis grid current is changed from zero to −0.2 pu. However, the q-axis reference grid current is usually set at zero to achieve a unity power factor. Figure 10 shows the quadrature grid current response under step changes. It can be observed that the Tyreus–Luyben method gives superior performance in comparison with the other reported methods. Here, it is important to mention that the gains of the PI controllers for the inner current loops at both RSC and GSC are selected thanks to the classical methods and not with the proposed optimization algorithms. Since the plant models of the inner loops are available, the use of the pole-placement, Tyreus–Luyben and Ziegler–Nichols methods is well adapted instead of metaheuristics-based methods which do not require models of systems to be controlled. The Bode plot of the grid current loop is presented in Figure 11.
The quality of the three phase grid currents is investigated through measuring the THD of the AC grid currents. As depicted in Figure 12. It is obviously shown that the use of the LCL-filter achieves a better attenuation with a THD value about 4.05% for which the IEEE519-1992 standard limits are respected. Further analysis can be made to compare the THD of the AC grid currents in the case of using an L-filter based structure. According to the results of Figure 12, it is clear that the THD for adopting the LCL-filer is smaller than that using the L-filer type at the same value of inductance.
Referring to the RSC control results, the stator active and reactive powers can be regulated by controlling the q-axis rotor current. The reference value of stator active power is generated from a Maximum Power Point Tracking (MPPT) strategy to extract the maximum power from the wind. In addition, the value of reference reactive power is set at different step levels to check the performance of the proposed PI controllers. As shown in Figure 13, both active and reactive powers track effectively their reference values with good performance in terms of speed and damping dynamics. In addition, the decoupling between the active and reactive powers is perfectly assured and the DFIG extracts the maximum available power, which is approximately about 1.5 MW. Figure 14 demonstrates the control performance of the PI controller for the reactive power loop around its final set-point value based on different tuned PI methods. The aim is to show the difference between the optimization tuning-based methods. On the other hand, both left and right sides of Figure 15 show the tracking performance of the direct and quadrature rotor currents, respectively. Such tracking dynamics are perfectly performed and improved thanks to the proposed TEO-based tuning and control method.

6.2. Statistical Analysis and Comparison

In this part, the mean execution values related to the different optimization criteria are sorted to assess the best operating one according to its average objective function performance. Moreover, a statistical comparison based on the nonparametric Friedman and Bonferroni–Dunn is carried out by using these mean performances [30,31]. The average ranks for all the proposed methods based on the four performances indices are provided in Table 6. One can note that the proposed TEO metaheuristic has worthily attained the lowest average ranks compared to the remaining methods. A statistical analysis has been performed to highlight the importance of the TEO-based tuning method over other algorithms [32].
The Friedman test for six competitor algorithms ( m = 6 ) and four indices ( l = 4 ) provides the computed value χ F 2 = 12.14 of the χ -distribution. The critical value of such a distribution with degrees of freedom m 1 = 5 and at confidence level α = 0.05 is equal to χ 5 , 0.05 2 = 11.07 . Since the above computed score is greater than this statistical value, the null hypothesis is declined. Moreover, for the Iman–Davenport test [30,31], the statistic is distributed with m 1 = 5 and ( m 1 ) × ( l 1 ) = 15 degrees of freedom. For this test, the null hypothesis is also rejected as the calculated value F F = 4.64 is greater than the critical value F 5 , 15 , 0.05 = 2.9 at the same significant level of confidence. All these statistical results indicate that there are significant differences among the performances of the reported algorithms for the optimization problem (19). Thus, the proposed TEO is found to be the most effective one having the best average Friedman ranking as given in Table 6. However, a Bonferroni–Dunn post-hoc test is called to investigate whether or not the proposed TEO algorithm is significantly better than another algorithm at the above considered level of confidence [32]. The corresponding Critical Differences (CD) of the reported algorithms at the confidence levels α = 0.05 (95% significance level) and α = 0.1 (90% significance level) are computed as CD 0.05 = 3.17 and CD 0.1 = 2.88 , respectively. Figure 16 illustrates the graphical representation of the Bonferroni–Dunn test considering the TEO as the control algorithm [32,33].
The bar of the TEO algorithm is the lowest high among the reported algorithms and the heights of the bars corresponding to the PSO algorithm. The GOA method violates the horizontal lines of significant levels. This reveals that the TEO algorithm performs at least significantly better than these two algorithms over the solutions equality.

6.3. Computational Time Efficiency

In this subsection, the average elapsed times of Table 3 are used to assess the computational efficiency of the reported algorithms. Such an algorithmic property of resource usage is quantified by the Computational Time Efficiency (CTE) metric defined as follows [34]:
C T E = E T a l g o r i t h m E T t o t a l ( % )
where E T a l g o r i t h m and E T t o t a l denote the average and total elapsed times for the same index.
According to the ET measures of problem (19), the CTE for each algorithm over the reported optimization criterion is summarized in Table 7.
In terms of computational efficiency, it can be observed from these results that the proposed TEO algorithm attained the second rank for IAE and ITAE criteria and the third and fourth ranks for ISE and ITSE criteria, respectively. Roughly, the HSA algorithm presents the best computational efficiency in terms of time resource usage. The GOA metaheuristic achieved the best average elapsed time as shown in Figure 17.

6.4. Sensitivity Analysis

In the TEO formalism, the parameters c 1 and c 2 are introduced to reduce the randomness of the algorithm over the course of iterations and to improve the exploration capabilities. The size of thermal memory TM is considered as a predefined parameter, which can enhance the performances of the algorithm without increasing the computational cost [3,24]. However, p r o is a user parameter that improves the global search capacity. For illustration purposes, the resolution of problem (19) is tested by changing these above main parameters. As shown in Table 8, the case with the parameters set ( c 1 = 1 , c 2 = 1 ) shows the superiority of the optimization results compared to other scenarios. In addition, it can be seen from Table 9 that the user parameter T M = 10 leads to obtain the best solutions. Finally, Table 10 investigates the effect of the parameter p r o on the algorithm performance. The convergence of the TEO algorithm is always guaranteed under all these control parameters variations which shows the insensitivity of the proposed physics-inspired metaheuristic.

7. Conclusions

This paper proposes an intelligent metaheuristics-based design procedure to tune the outer-loop PI controllers for the direct power control scheme of DFIG-based energy conversion systems. In this study, only outer PI loops in the RSC and GSC circuits are optimized for the active and reactive powers and DC-link voltage dynamics. Since the reported classical techniques of PI controller tuning are tedious, time-consuming and not systematic, a TEO-based approach has been proposed and successfully implemented. The PI controllers tuning problem is firstly formulated as a constrained optimization program under nonlinear and non-smooth operational constraints. The introduced TEO algorithm is then employed separately to minimize several time-domain performance criteria such as IAE, ISE, ITAE and ITSE indices as objective functions. The proposed TEO-tuned PI controllers methodology improves the performance and robustness of the controlled DFIG-based energy converter in terms of rising time, settling time, and steady state indices. The classical trials-errors based methods of PI controllers tuning are no longer used and the design time is further reduced. To evaluate the performance superiority of the proposed TEO-based approach in finding the global minimum value of the objective function for various performance indices, a comparison study with the PSO, GA, HSA, GSO and WCA is performed. The demonstrative results exhibit that the proposed TEO gives very complete results in terms of global search capabilities, robustness and non-premature convergence. Finally, the statistical analysis is achieved by using the Friedman’s rank and Bonferroni-Dunn’s test. The corresponding results show that the proposed TEO-based method is a promising alternative approach for controlling the DFIG system by systematically tuning the unknown PI controllers’ parameters efficiently.

Author Contributions

Conceptualization, M.M.A. and S.B.; Methodology, S.B.; Software, M.M.A.; Validation, M.M.A.; Formal analysis, S.B.; Investigation, S.B.; Resources, M.M.A. and S.B.; Writing—original draft preparation, M.M.A.; Writing—review and editing, S.B.; Visualization, S.B.; Supervision, S.B.

Funding

This research received no external funding.

Acknowledgments

The authors would like to express their sincere appreciation to the anonymous reviewers for their valuable and constructive comments, which greatly improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CDCritical Difference
CTEComputational Time Efficiency
ETElapsed Time
DFIGsDoubly Fed Induction Generators
HSAHarmony Search Algorithm
IAEIntegral Absolute Error
IGBTsInsulated-Gate Bipolar Transistors
ISEIntegral Square Error
ITAEIntegral Time-weighted Absolute Error
ITSEIntegral Time-weighted Square Error
GAGenetic Algorithm
GOAGrasshopper Optimization Algorithm
GSCGrid Side Converter
MPPTMaximum Power Point Tracking
PIProportional-Integral
PSOParticle Swarm Optimization
RSCRotor Side Converter
SFOStator Flux Orientation
SPWMSinusoidal Pulse Width Modulation
STDStandard Deviation
TEOThermal Exchange Optimization
THDTotal Harmonic Distortion
TMThermal Memory
VOCVoltage Oriented Control
WECSWind Energy Converter System
WCAWater Cycle Algorithm
WTWind Turbine
WTsWind Turbines

Notations

C d c DC-link capacitance
C f Filter capacitance
C p Power conversion coefficient
e g ( d , q ) d-q axis grid voltages
i g ( d , q ) d-q axis grid currents
i r ( d , q ) d-q axis rotor currents
i s ( d , q ) d-q axis stator currents
L g Filter grid side inductance
L i Filter converter side inductance
L m Magnetizing inductance
L r , L s Rotor and stator inductances
L T Filter total inductance
P g , P s Grid and stator active powers
P m Mechanical turbine power
Q g , Q s Grid and stator reactive powers
R Turbine radius
R d Filter damping resistance
R r , R s Rotor and stator resistances
V d c DC-link voltage
V f ( d , q ) d-q axis grid converter voltage sides
V r   ( d , q ) d-q axis rotor converter voltage sides
V s   ( d , q ) d-q axis stator voltages
V w Wind speed
λ Tip speed ratio
ω r , ω s Rotor and stator angular frequencies
Ω t Mechanical rotational speed
ρ Air density
φ d r , φ q r d-q axis rotor fluxes
φ d s , φ q s d-q axis stator fluxes

Appendix A. Decision Variables of Problem (19) Relative to the Optimization Mean Case

The following gains of the optimized PI controllers (Table A1), as decision variables of problem (19), are used to perform all numerical simulations of the paper.
Table A1. Comparative results of the PI controllers’ coefficients tuning.
Table A1. Comparative results of the PI controllers’ coefficients tuning.
IndicesAlgorithmsPI Controllers’ Gains
K p P s K i P s K p Q s K i Q s K p d c K i d c
IAEPSO14.16139.37376.43381.82931.32
GA16.92394.6218.22147.4624.3662.17
HSA7.8983.96151.27227.7228.0386.11
WCA2.28329.93122.5540026.1141.53
GOA16.91394.6218.22147.4726.3659.85
TEO4.2358.4310.3715.7389.10192.71
ISEPSO11.23207.70169.3135369.65272.40
GA2.25313.46187.62398.8279.95345.49
HSA1 × 10−51.87118.84360.1633.2939.65
WCA12.88143.34360.15400191.11293
GOA12.42388.971.93371.52148.46220.78
TEO35.94165.2515.4471.591.483
ITSEPSO0.11 × 10−5238.93351.810.9925.03
GA1 × 10−51 × 10−57.27214.4224.3462.20
HSA303.71320.521.036.2527.6645.22
WCA4.13248.9810.7924.4527.8979.04
GOA9.82108.24333.2527.36131.74205.91
TEO8.565.0467.7269.4410.7131.96
ITAEPSO250.16289.53346.50142.3473.74238.53
GA1 × 10−51 × 10−545.0285.22130.5726.43
HSA8.164.83323.56384.1933.3347.41
WCA9.5463.73329.42395.36143.52400
GOA13.71399.92302.9265.497.91196
TEO4.8 × 10−611.281.6 × 10−68.3 × 10−610.2468.66

References

  1. Sahu, K.B. Wind Energy Developments and Policies in China: A Short Review. Renew. Sustain. Energy Rev. 2018, 81, 1393–1405. [Google Scholar] [CrossRef]
  2. Belmokhtar, K.; Doumbia, M.L.; Agbossou, K. Modelling and Fuzzy Logic Control of DFIG Based Wind Energy Conversion Systems. In Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE), Hangzhou, China, 28–31 May 2012; pp. 1888–1893. [Google Scholar]
  3. Hu, J.; Nian, H.; Xu, H.; He, Y. Dynamic Modeling and Improved Control of DFIG under Distorted Grid Voltage Conditions. IEEE Trans. Energy Convers. 2011, 26, 163–175. [Google Scholar] [CrossRef]
  4. Liserre, M.; Blaabjerg, F.; Hansen, S. Design and Control of an LCL-Filter-based Three-Phase Active Rectifier. IEEE Trans. Ind. Appl. 2005, 50, 1281–1291. [Google Scholar] [CrossRef]
  5. Zhou, D.; Blaabjerg, F. Bandwidth Oriented Proportional-Integral Controller Design for Back-To-Back Power Converters in DFIG Wind Turbine System. IET Renew. Power Gen. 2017, 11, 941–951. [Google Scholar] [CrossRef]
  6. Hamane, B.; Doumbia, M.L.; Bouhamida, M.; Benghanem, M. Direct Active and Reactive Power Control of DFIG based WECS using PI and Sliding Mode Controllers. In Proceedings of the 40th Annual Conference of the IEEE Industrial Electronics Society, Dallas, TX, USA, 29 October–1 November 2014; pp. 2050–2055. [Google Scholar]
  7. Hamane, B.; Doumbia, M.L.; Bouhamida, M.; Draou, A.; Chaoui, H.; Benghanem, M. Comparative Study of PI, RST, Sliding Mode and Fuzzy Supervisory Controllers for DFIG based Wind Energy Conversion System. Int. J. Renew. Energy Res. 2015, 5, 1174–1185. [Google Scholar]
  8. Hamane, B.; Doumbia, M.L.; Bouhamida, M.; Draou, A.; Chaoui, H.; Benghanem, M. Modeling and Control of a Wind Energy Conversion System Based on DFIG Driven by a Matrix Converter. In Proceedings of the Eleventh International Conference on Ecological Vehicles and Renewable Energies, Monte Carlo, Monaco, 6–8 April 2016; pp. 1–8. [Google Scholar]
  9. Tanvir, A.A.; Merabet, A.; Beguenane, R. Real-Time Control of Active and Reactive Power for Doubly Fed Induction Generator (DFIG)-Based Wind Energy Conversion System. Energies 2015, 8, 10389–10408. [Google Scholar] [CrossRef] [Green Version]
  10. Bekakra, Y.; Attous, D.B. Optimal Tuning of PI Controller Using PSO Optimization for Indirect Power Control for DFIG Based Wind Turbine with MPPT. Int. J. Syst. Assur. Eng. Manag. 2014, 5, 219–229. [Google Scholar] [CrossRef]
  11. Vieira, J.P.A.; Nunes, M.N.A.; Bezerra, U.H. Design of Optimal PI Controllers for Doubly Fed Induction Generators in Wind Turbines Using Genetic Algorithm. In Proceedings of the IEEE Power and Energy Society General Meeting—Conversion and Delivery of Electrical Energy in the 21st Century, Pittsburgh, PA, USA, 20–24 July 2008; pp. 1–7. [Google Scholar]
  12. Asgharnia, A.; Shahnazi, R.; Jamali, A. Performance and Robustness of Optimal Fractional Fuzzy PID Controllers for Pitch Control of a Wind Turbine Using Chaotic Optimization Algorithms. ISA Trans. 2018, 79, 27–44. [Google Scholar] [CrossRef]
  13. Assareh, E.I.; Poultangari, I.; Tandis, E.; Nedaei, M. Optimizing the Wind Power Generation in Low Wind Speed Areas Using an Advanced Hybrid RBF Neural Network Coupled with the HGA-GSA Optimization Method. J. Mech. Sci. Technol. 2016, 30, 4735–4745. [Google Scholar] [CrossRef]
  14. Kaveh, A.; Dadras, A. A Novel Meta-Heuristic Optimization Algorithm: Thermal Exchange Optimization. Adv. Eng. Softw. 2017, 110, 69–84. [Google Scholar] [CrossRef]
  15. Gagra, S.K.; Mishra, S.; Sing, M. Performance Analysis of Grid Integrated Doubly Fed Induction Generator for a Small Hydropower Plant. Int. J. Renew. Energy Res. 2018, 8, 2310–2323. [Google Scholar]
  16. Qiao, W. Dynamic Modeling and Control of Doubly Fed Induction Generators Driven by Wind Turbines. In Proceedings of the IEEE/PES Power Systems Conference and Exposition, Seattle, DC, USA, 15–18 March 2009; pp. 1–8. [Google Scholar]
  17. Xu, D.; Blaabjerg, F.; Chen, W.; Zhu, N. Advanced Control of Doubly Fed Induction Generator for Wind Power Systems, 1st ed.; John Wiley & Sons: New York, NY, USA, 2018. [Google Scholar]
  18. Hato, M.M.; Bouallègue, S.; Ayadi, M. Water cycle algorithm-tuned PI control of a doubly fed induction generator for wind energy conversion. In Proceedings of the 9th International Renewable Energy Congress (IREC), Hammamet, Tunisia, 20–22 March 2018; pp. 1–6. [Google Scholar]
  19. Åström, A.K.; Hägglund, T. Advanced PID Control; ISA-Instrumentation, Systems and Automation Society: Research Triangle Park, NC, USA, 2006. [Google Scholar]
  20. Bouallègue, S.; Haggège, J.; Ayadi, M.; Benrejeb, M. PID-Type Fuzzy Logic Controller Tuning based on Particle Swarm Optimization. Eng. Appl. Artif. Intel. 2012, 25, 484–493. [Google Scholar] [CrossRef]
  21. Qu, B.Y.; Zhu, Y.S.; Jiao, Y.S.; Wu, M.Y.; Suganthan, P.N.; Liang, J.J. A Survey on Multi-objective Evolutionary Algorithms for the Solution of the Environmental/Economic Dispatch Problems. Swarm Evol. Comput. 2018, 8, 1–11. [Google Scholar] [CrossRef]
  22. Gong, L.; Cao, W.; Liu, K.; Zhao, J.; Li, X. Spatial and Temporal Optimization Strategy for Plug-in Electric Vehicle Charging to Mitigate Impacts on Distribution Network. Energies 2018, 11, 1373. [Google Scholar] [CrossRef]
  23. Zhang, X.; Yang, J.; Wang, W.; Jing, T.; Zhang, M. Optimal Operation Analysis of The Distribution Network Comprising a Micro Energy Grid based on an Improved Grey Wolf Optimization Algorithm. Appl. Sci. 2018, 8, 923. [Google Scholar] [CrossRef]
  24. Kaveh, A.; Ghazaan, M.I. Enhanced Colliding Bodies Optimization for Design Problems with Continuous and Discrete Variables. Adv. Eng. Softw. 2014, 77, 66–75. [Google Scholar] [CrossRef]
  25. Geem, Z.W.; Kim, J.H.; Loganathan, G.V. A New Heuristic Optimization Algorithm: Harmony Search. Sage J. 2012, 76, 60–68. [Google Scholar]
  26. Saremi, S.; Mirjalili, S.; Lewis, A. Grasshopper Optimization Algorithm: Theory and Application. Adv. Eng. Softw. 2017, 105, 30–47. [Google Scholar] [CrossRef]
  27. Blasko, V.; Kaura, V. A Novel Control to Actively Damp Resonance in Input LC Filter of a Three-Phase Voltage Source Converter. IEEE Trans. Ind. Appl. 1997, 33, 542–550. [Google Scholar] [CrossRef]
  28. Bingi, K.; Ibrahim, R.; Karsiti, M.N.; Chung, T.D.; Hassan, S.M. Optimal PID control of pH neutralization plant. In Proceedings of the 2nd IEEE International Symposium on Robotics and Manufacturing Automation, Ipoh, Malaysia, 25–27 September 2016; pp. 1–6. [Google Scholar]
  29. Mishra, K.P.; Kumar, V.; Rana, K.P.S. Stiction combating intelligent controller tuning a comparative study. In Proceedings of the International Conference on Advances in Computer Engineering and Application, Ghaziabad, India, 19–20 March 2015; pp. 534–541. [Google Scholar]
  30. Raghuwanshi, B.S.; Shukla, S. Generalized Class-Specific Kernelized Extreme Learning Machine for Multiclass Imbalanced Learning. Expert. Syst. 2019, 121, 244–255. [Google Scholar] [CrossRef]
  31. Zeng, F.; Zhang, W.; Zhang, S.; Zheng, N. Re-Kissme: A Robust Resampling Scheme for Distance Metric Learning in the Presence of Label Noise. Neurocomputing 2019, 330, 138–150. [Google Scholar] [CrossRef]
  32. Sadhu, A.K.; Konar, A.; Bhattacharjee, A.T.; Das, S. Synergism of Firefly Algorithm and Q-Learning for Robot Arm Path Planning. Swarm Evol. Comput. 2018, 43, 50–68. [Google Scholar] [CrossRef]
  33. Liu, Y.; Chen, S.; Guan, B.; Xu, P. Layout Optimization of Large-Scale Oil–Gas Gathering System based on Combined Optimization Strategy. Neurocomputing 2019, 332, 159–183. [Google Scholar] [CrossRef]
  34. Sababha, M.; Zohdy, M.; Kafafy, M. The Enhanced Firefly Algorithm Based on Modified Exploitation and Exploration Mechanism. Electronics 2018, 7, 132. [Google Scholar] [CrossRef]
Figure 1. Configuration of the Doubly Fed Induction Generator (DFIG)-based Wind Energy Converter System (WECS).
Figure 1. Configuration of the Doubly Fed Induction Generator (DFIG)-based Wind Energy Converter System (WECS).
Mca 24 00077 g001
Figure 2. Proportional-Integral (PI) control scheme for the inner current loop of the RSC.
Figure 2. Proportional-Integral (PI) control scheme for the inner current loop of the RSC.
Mca 24 00077 g002
Figure 3. Design of the PI controller current loop of grid side converter.
Figure 3. Design of the PI controller current loop of grid side converter.
Mca 24 00077 g003
Figure 4. The proposed control algorithm applied to the DFIG system.
Figure 4. The proposed control algorithm applied to the DFIG system.
Mca 24 00077 g004
Figure 5. Pairs of environment and cooling objects.
Figure 5. Pairs of environment and cooling objects.
Mca 24 00077 g005
Figure 6. Flowchart of the proposed Thermal Exchange Optimization (TEO)-tuned PI controllers’ parameters.
Figure 6. Flowchart of the proposed Thermal Exchange Optimization (TEO)-tuned PI controllers’ parameters.
Mca 24 00077 g006
Figure 7. Convergence rates comparison: (a) IAE criterion; (b) ISE criterion; (c) ITAE criterion; (d) ITSE criterion.
Figure 7. Convergence rates comparison: (a) IAE criterion; (b) ISE criterion; (c) ITAE criterion; (d) ITSE criterion.
Mca 24 00077 g007aMca 24 00077 g007b
Figure 8. Box-and-Whisker plot of optimization problem (19): (a) IAE criterion; (b) ISE criterion; (c) ITAE criterion; (d) ITSE criterion.
Figure 8. Box-and-Whisker plot of optimization problem (19): (a) IAE criterion; (b) ISE criterion; (c) ITAE criterion; (d) ITSE criterion.
Mca 24 00077 g008aMca 24 00077 g008b
Figure 9. DC-link voltage responses against step changes for different-tuned PI controllers.
Figure 9. DC-link voltage responses against step changes for different-tuned PI controllers.
Mca 24 00077 g009
Figure 10. Quadrature grid current responses against a step change for different-tuned PI controllers.
Figure 10. Quadrature grid current responses against a step change for different-tuned PI controllers.
Mca 24 00077 g010
Figure 11. Bode plot of the open loop grid current loop: with and without passive damping based technical optimum tuning.
Figure 11. Bode plot of the open loop grid current loop: with and without passive damping based technical optimum tuning.
Mca 24 00077 g011
Figure 12. Harmonics spectrum of the AC grid currents.
Figure 12. Harmonics spectrum of the AC grid currents.
Mca 24 00077 g012
Figure 13. Time-domain variations of the TEO-tuned PI controlled active and reactive powers.
Figure 13. Time-domain variations of the TEO-tuned PI controlled active and reactive powers.
Mca 24 00077 g013
Figure 14. Performance comparison of the controlled reactive power under different PI controllers’ tuning methods.
Figure 14. Performance comparison of the controlled reactive power under different PI controllers’ tuning methods.
Mca 24 00077 g014
Figure 15. Time-domain variations of the TEO-based control of the direct and quadrature rotor currents.
Figure 15. Time-domain variations of the TEO-based control of the direct and quadrature rotor currents.
Mca 24 00077 g015
Figure 16. Graphical representation of Bonferroni–Dunn’s test for problem (19).
Figure 16. Graphical representation of Bonferroni–Dunn’s test for problem (19).
Mca 24 00077 g016
Figure 17. Average elapsed time over all criteria for the introduced algorithms.
Figure 17. Average elapsed time over all criteria for the introduced algorithms.
Mca 24 00077 g017
Table 1. Parameters for a 1.5 MW DFIG systems.
Table 1. Parameters for a 1.5 MW DFIG systems.
EquipmentParameterValueUnit
DFIGRated power1500kW
RMS grid line voltage575V
Slip range0.2-
Rated electrical frequency50Hz
Stator resistance0.023pu
Stator leakage inductance0.18pu
Rotor resistance0.016pu
Rotor leakage inductance0.16pu
Magnetizing inductance2.9pu
Power convertersRated power300kW
Switching frequency at f s w , R S C and f s w , G S C 2700Hz
DC-link voltage1050V
DC-link capacitance10000 μ F
L-filterL-filter grid inductance0.1pu
LCL-filterLCL-filter grid side inductance0.018pu
LCL-filter capacitance0.104pu
LCL-filter converter side inductance0.077pu
Passive damping resistor0.124pu
Table 2. Parameter setting of reported algorithms.
Table 2. Parameter setting of reported algorithms.
AlgorithmsParameters Setting
PSOCognitive and social coeffs. c 1 = c 2 = 2 , weights w max = 0.9 , w min = 0.2 [10].
GARoulette wheel, crossover P c r o s s = 1 , mutation P m u t = 0.01 [11].
HSAHarmony memory rate HMCR = 0.9, pitch adjusting rate PAR = 0.3 [25].
WCASummation number of rivers N s r = 8 and d m a x = 1 × 10 3 [18].
GOA c max = 1 , c min = 0.00001 [26].
TEOThermal memory T M = 10 , p r o = 0.50 , c 1 = 1 , c 2 = 1 [3].
Table 3. Statistical results of optimization problem (19) over 10 independent runs.
Table 3. Statistical results of optimization problem (19) over 10 independent runs.
Indices PSOGAHSAWCAGOATEO
IAEBest1.571.621.571.521.601.55
Mean1.731.671.601.592.011.58
Worst2.691.731.641.783.201.60
STD3.4 × 10−13.2 × 10−22.3 × 10−27.6 × 10−26.3 × 10−11.9 × 10−2
ET (sec)290222872025480329002040024640
ISEBest 33.1734.8632.6432.6634.0631.45
Mean41.4237.2535.5145.8845.0633.46
Worst61.8040.5338.7861.2464.5036.79
STD8.271.921.8010.779.081.64
ET (sec)281402742019083269001920023420
ITSEBest17.4717.4918.4517.1417.8617.79
Mean18.4617.9623.1818.3219.8318.24
Worst20.4518.2126.7719.7528.1119.20
STD0.9210.2073.6168.1E-12.964.1E-1
ET (sec)215602628018720275401950023380
ITAEBest0.6500.6320.6460.6400.6440.642
Mean0.6830.6500.6830.6590.6880.658
Worst0.7070.620.7180.7360.7420.675
STD0.0180.0190.0210.0290.0400.014
ET (sec)200602048016880258802084019680
Table 4. Time-domain performances for controlled DC-link voltage under step changes scenario.
Table 4. Time-domain performances for controlled DC-link voltage under step changes scenario.
PI Tuning MethodsUnit Step Change Response
t r   ( sec ) t s   ( sec ) t p   ( sec ) δ   ( % ) E s s
Frequency response0.04414.38134.13821.1680.3355
Pole placement0.01634.11114.04412.0060.2813
Symmetrical optimum0.01854.13244.04592.7680.3568
Ziegler–Nichols0.00374.02524.00801.2090.2843
Tyreus–Luyben0.00604.02494.01820.2990.2977
TEO0.00294.01934.00762.5370.2539
Table 5. Time-domain performances for controlled quadrature grid current under a step change scenario.
Table 5. Time-domain performances for controlled quadrature grid current under a step change scenario.
PI Tuning MethodsUnit Step Change Response
t r   ( sec ) t p   ( sec ) δ   ( % ) E s s
Frequency response0.01712.038728.18020.0035
Pole placement0.00272.009621.950.0045
Technical optimum0.01462.03838.29310.0055
Ziegler–Nichols3.582 × 10−42.002610.74000.0051
Tyreus–Luyben0.00202.01557.16570.0061
Table 6. Average Rank based statistical analysis of mean performances.
Table 6. Average Rank based statistical analysis of mean performances.
AlgorithmsIndicesAverage Rank
IAEISEITSEITAE
ScoreRankScoreRankScoreRankScoreRank
PSO1.73541.42418.4640.68354.25
GA1.67437.25317.9610.65012.25
HSA1.60335.51223.1860.68343.75
WCA1.59245.88618.3230.65933.75
GOA2.01644.06519.8350.68865.5
TEO1.58133.46118.2420.65821.5
Table 7. The computational time efficiency (CTE) of each algorithm over the four indices.
Table 7. The computational time efficiency (CTE) of each algorithm over the four indices.
IndexAlgorithms
PSOGAHSAWCAGOATEO
IAE18.01%17.82%15.81%20.41%12.66%15.29%
ISE19.52%19.02%13.24%18.66%13.32%16.25%
ITSE15.74%19.19%13.67%20.11%14.24%17.07%
ITAE16.20%16.54%13.63%20.90%16.83%15.89%
Table 8. Comparison under different values of c 1 and c 2 : IAE criterion.
Table 8. Comparison under different values of c 1 and c 2 : IAE criterion.
Optimization Scenario T M = 10 ,   p r o = 0.5   and   N i t e r = 100
WorstMeanBestSTD
c 1 = 0 , c 2 = 0 1.6691.6521.6381.9 × 10−2
c 1 = 1 , c 2 = 0 1.6641.6401.6221.4 × 10−2
c 1 = 0 , c 2 = 1 1.6541.6401.6317.5 × 10−3
c 1 = 1 , c 2 = 1 1.6081.5821.5531.9 × 10−2
Table 9. Comparison under different values of T M : IAE criterion.
Table 9. Comparison under different values of T M : IAE criterion.
Optimization Scenario c 1 = 1 ,   c 2 = 1 ,   p r o = 0.5   and   N i t e r = 100
WorstMeanBestSTD
T M = 4 1.6991.6941.6893.5 × 10−3
T M = 8 1.6541.6401.6318.4 × 10−2
T M = 10 1.6081.5821.5531.9 × 10−2
Table 10. Comparison under different values of p r o : IAE criterion.
Table 10. Comparison under different values of p r o : IAE criterion.
Optimization Scenario c 1 = 1 ,   c 2 = 1 ,   T M = 10   and   N i t e r = 100
WorstMeanBestSTD
p r o = 0.20 1.6931.6541.6331.7 × 10−2
p r o = 0.35 1.6781.6171.5755.2 × 10−2
p r o = 0.50 1.6081.5821.5531.9 × 10−2

Share and Cite

MDPI and ACS Style

Alhato, M.M.; Bouallègue, S. Direct Power Control Optimization for Doubly Fed Induction Generator Based Wind Turbine Systems. Math. Comput. Appl. 2019, 24, 77. https://doi.org/10.3390/mca24030077

AMA Style

Alhato MM, Bouallègue S. Direct Power Control Optimization for Doubly Fed Induction Generator Based Wind Turbine Systems. Mathematical and Computational Applications. 2019; 24(3):77. https://doi.org/10.3390/mca24030077

Chicago/Turabian Style

Alhato, Mohammed Mazen, and Soufiene Bouallègue. 2019. "Direct Power Control Optimization for Doubly Fed Induction Generator Based Wind Turbine Systems" Mathematical and Computational Applications 24, no. 3: 77. https://doi.org/10.3390/mca24030077

APA Style

Alhato, M. M., & Bouallègue, S. (2019). Direct Power Control Optimization for Doubly Fed Induction Generator Based Wind Turbine Systems. Mathematical and Computational Applications, 24(3), 77. https://doi.org/10.3390/mca24030077

Article Metrics

Back to TopTop