Advanced
Multi-objective Fuzzy-optimization of Crowbar Resistances for the Low-Voltage Ride-through of Doubly Fed Induction Wind Turbine Generation Systems
Multi-objective Fuzzy-optimization of Crowbar Resistances for the Low-Voltage Ride-through of Doubly Fed Induction Wind Turbine Generation Systems
Journal of Power Electronics. 2015. Jul, 15(4): 1119-1130
Copyright © 2015, The Korean Institute Of Power Electronics
  • Received : November 06, 2014
  • Accepted : March 04, 2015
  • Published : July 20, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Wenjuan Zhang
Department of Electrical Engineering, Baoji University of Arts and Sciences, Shaanxi, China
zhangwj19811130@163.com
Haomiao Ma
College of Automation and Information Engineering, Xi’an University of Technology, Shaanxi, China
Junli Zhang
Department of Electrical Engineering, Baoji University of Arts and Sciences, Shaanxi, China
Lingling Chen
94816 troop, China
Yang Qu
Department of Electrical Engineering, Baoji University of Arts and Sciences, Shaanxi, China

Abstract
This study investigates the multi-objective fuzzy optimization of crowbar resistance for the doubly fed induction generator (DFIG) low-voltage ride-through (LVRT). By integrating the crowbar resistance of the crowbar circuit as a decision variable, a multi-objective model for crowbar resistance value optimization has been established to minimize rotor overcurrent and to simultaneously reduce the DFIG reactive power absorbed from the grid during the process of LVRT. A multi-objective genetic algorithm (MOGA) is applied to solve this optimization problem. In the proposed GA, the value of the crowbar resistance is represented by floating-point numbers in the GA population. The MOGA emphasizes the non-dominated solutions and simultaneously maintains diversity in the non-dominated solutions. A fuzzy-set-theory-based is employed to obtain the best solution. The proposed approach has been evaluated on a 3 MW DFIG LVRT. Simulation results show the effectiveness of the proposed approach for solving the crowbar resistance multi-objective optimization problem in the DFIG LVRT.
Keywords
I. INTRODUCTION
In the past few decades, the utilization of doubly fed induction generators (DFIG) has increased significantly due to their advantages such as variable-speed-constant-frequency (VSCF) operation, decoupled control of the active and reactive powers, and availability of partial-scale converters [1] . However, since the stator of a DFIG is directly connected to the grid, this promising wind generator technology is vulnerable to grid disturbances, especially to voltage sag faults [2] . A grid voltage sag imposed at the connection point of a DFIG to the grid could easily induce a DFIG rotor overcurrent, DC-link overvoltage, and torque oscillations that could lead to the destruction of the rotor-side converter, the DC-link capacitor, and other mechanical parts, respectively [3] . Thus, it is desirable that wind turbines remain connected to the grid and actively contribute to the system stability during and after grid voltage sag faults and disturbances. This ability of wind turbines is termed low-voltage ride-through (LVRT) [4] , [5] .
One effective way to make a DFIG capable of LVRT is to equip it with an appropriate crowbar protection circuit in the rotor side of the DFIG. However, since the transient process of the DFIG is owe damped; this will result in the grid voltage becoming lower, because when the crowbar is in operation, the insulated-gate bipolar transistors (IGBTs) of the converters will quit, and the DFIG behaves like a squirrel-cage induction machine, absorbing a large amount of reactive power from the grid [6] . In addition, to select a reasonable value for the crowbar resistance is important. The larger the resistance, the faster the attenuation of the transient components. However, a larger resistance may cause rotor-side overvoltage and may also damage the rotor-side converter.
A good amount of literature is available on the transient characteristics and crowbar resistance design aspects of DFIGs. In [7] , Niiranen analyzes the DFIG equivalent circuit and points out that an appropriate crowbar resistance is helpful for LVRT. In addition, [8] demonstrates that a detailed dynamic representation of the converter circuitry is necessary to investigate the impact of various controls on DFIG performance. Morren reports the short-circuit maximum current and short capacity of the DFIG, and the value range of the crowbar resistance is approximately calculated in [6] . In [9] , the short circuit current analytical expression is given and the decay time constant is analyzed after the crowbar protection is put into operation. In [10] , based on a linear model of the maximum short-circuit current, a crowbar resistance design method is put forward that considers the DC bus clamping effect.
The aforementioned studies provide references for the design of the crowbar circuit. However, many of these studies have two limitations. First, the crowbar resistance is strictly confined by the safety range of the rotor current and voltage. In the case of optimization design by the common method, the constraints are deemed as “rigid.” In addition, considering that the constraints are obtained from the point of view of safe service of the double fed converter, when a crowbar resistance value is selected on the constraint's boundary, the system works on the boundary between the safe area and the dangerous area. As a result of such a design, changes to the DFIG parameters caused by heating are likely to trigger the protection operation of the DFIG, or even lead to the destruction of the converter. Second, the existing design of the crowbar resistance merely considers how to keep rotor current and voltage within reasonable ranges, while overlooking many other properties of LVRT. For example, when the crowbar is in operation, the DFIG absorbs reactive power from the grid, leading to a further drop in the grid voltage. This is harmful to the recovery of the grid voltage. It is worth noting that the ability of the DFIG to absorb reactive power depends on the value of the crowbar resistance. Hence, the design of the crowbar resistance should take into account LVRT, which is a multi-objective fuzzy optimization problem.
Based on this fact, this paper proposes a method of using fuzzy sets to describe the crowbar resistance value range, and adopts a multi-objective genetic algorithm in order to optimize the value of the crowbar resistance. The optimized crowbar resistance can provide the best compromise solution during the process of LVRT. This also avoids that the value of the crowbar resistance be taken on the constraint boundary. These considerations will improve the ability of the DFIG in terms of LVRT and enhance the safety of the system. This paper is organized as follows. In Section 2, the behavior of the DFIG before and after grid faults is analyzed. In Section 3, the proposed multi-objective fuzzy-optimization method is discussed, and in Section 4, system simulation results are shown. Finally, some conclusion are presented in Section 5.
II. BEHAVIOR OF THE DFIG BEFORE AND AFTER GRID FAULTS
Fig. 1 gives the DFIG wind turbine topological structure of the crowbar protection circuit. In this structure, a three-phase non-controlled rectifier bridge is paralleled on the rotor side, and the DC side of the rectifier bridge is connected with the crowbar resistance. When the units work normally, the rotor voltage is rectified on the crowbar, and it maintains a fixed DC voltage level at both ends of the crowbar resistance. When the grid voltage sags and the rotor current increases to the protection threshold, the crowbar circuit operates. At the same time, the rotor side converter is locked to provide an unloading path for the rotor-side fault current for the purpose of converter protection. The following sections analyze the performance of the DFIG before and after grid faults.
PPT Slide
Lager Image
DFIG crowbar protection in the rotor circuit.
- A. Mathematical Model of the DFIG before Grid Faults
The stator and rotor voltages and fluxes in a synchronously rotating reference frame are given in [11] , where the rotor variables and parameters are referred to the stator.
PPT Slide
Lager Image
PPT Slide
Lager Image
In these equations, all of the parameters are reduced to the stator side, where us = [ usd usq ] T and ur = [ urd urq ] T are the stator and rotor voltage vectors, is = [ isd isq ] T and ir = [ ird irq ] T are the stator and rotor current vectors, ψs = [ ψsd ψsq ] T and ψr = [ ψrd ψrq ] T are the stator and rotor flux vectors, while L s and L r are the stator and rotor self-inductances, respectively. L m is the mutual inductance, ω s is the stator frequency, and ω r is the rotor slip frequency. In addtion, R s and R r are the stator and rotor resistances.
Based on these equations, the equivalent circuit of the DFIG is shown in Fig. 2 . It can be used for transient analysis of an induction machine.
PPT Slide
Lager Image
Equivalent circuit of induction machine for transient analysis.
Before the occurrence of a voltage dip, at t = t 0 , the stator voltage space vector is a vector with a constant amplitude u 0 that rotates at the synchronous electrical angular frequency ω s . Thus, the stator voltage space vector can be described as:
PPT Slide
Lager Image
By substituting Equation (2) into (1), the stator voltage orientation is used, per unit ω s = 1, the slip s = ( ω s ω r )/ ω s = 1 − ω r , the steady rotor current before a voltage dip is obtained as:
PPT Slide
Lager Image
where ω b is the base angular frequency.
The stator resistance can be neglected in the steady state. Based on (2) and (3), the expressions for the stator flux and rotor flux in the steady state are obtained as:
PPT Slide
Lager Image
- B. Behavior of the DFIG after Grid Faults
1) Rotor Current Analysis: The rotor current after a voltage dip consists of two components: the first is a transient attenuation component that is caused by the grid voltage dip, and the second is a steady state current attenuation component.
For the transient attenuation component at time t = t 0 , a three-phase dip fault is assumed to have occurred at the stator of the DFIG. According to Equation (3), the stator voltage expressions are:
PPT Slide
Lager Image
where p is the voltage dip ratio. This implies that the flux in both of the windings does not change. The stator and rotor flux are given by:
PPT Slide
Lager Image
This transient flux has an exponential decay. The damping time constant of the dc components in the stator and rotor can be defined as:
PPT Slide
Lager Image
PPT Slide
Lager Image
In these equations, according to [12] , the transient stator inductance can be derived as:
PPT Slide
Lager Image
Similarly, according to [12] , the transient rotor inductance can be introduced as:
PPT Slide
Lager Image
In addition, as per [13] , the change of the stator and rotor flux after a voltage dip can be expressed mathematically in the rotor reference frame as follows:
PPT Slide
Lager Image
Since a detailed analysis of the transient rotor current,
PPT Slide
Lager Image
induced by the transient flux
PPT Slide
Lager Image
can be found in [14] it will not be described in this paper. Only the resulting dynamic current expression is given as follows:
PPT Slide
Lager Image
By substituting Equation (12) into (13), the general transient rotor current is given by:
PPT Slide
Lager Image
For the steady state current attenuation component i r(0++) , based on (4), the term can be written as:
PPT Slide
Lager Image
Thus, by adding Equations (14) and (15), the final expression of the rotor currents during a voltage dip can be written as follows:
PPT Slide
Lager Image
2) Analysis of the DFIG Absorbing Reactive Power: When a voltage dip has occurred and the rotor overcurrent is large enough to trigger the crowbar protection, the DFIG may lose its controllability and behaves as a squirrel-cag induction machine. The impedance equivalent circuit is shown in Fig. 3 where R cb is the crowbar reactance.
PPT Slide
Lager Image
Impedance equivalent circuit of the DFIG with crowbar.
Based on Fig. 3 , the equivalent impedance of the DFIG with the crowbar can be described as:
PPT Slide
Lager Image
where:
PPT Slide
Lager Image
PPT Slide
Lager Image
If the influence of the transient current for the steady state power of the DFIG is ignored, based on the relationship between the active power and reactive power, it is possible to write:
PPT Slide
Lager Image
The reactive power that the DFIG absorbs from the grid during a fault can be obtained as:
PPT Slide
Lager Image
From Eq. (21), it can be seen that the reactive power is determined by the stator voltage, the value of the crowbar resistance, and the rotor speed of the DFIG. This absorbed reactive power will lead to a further voltage drop in the grid. Therefore, the absorbed reactive power should be as little as possible.
III.MULTI-OBJECTIVE FUZZY OPTIMIZATION STRATEGY OF CROWBAR RESISTANCES
The crowbar protection of the LVRT in a doubly fed wind power generator should have the following optimization goals. First, it should achieve the safest generator operation performance. Second, it should possess the fastest fault recovery. Based on the squirrel-cage induction machine operating characteristic of the DFIG after the crowbar is applied, the optimization objective function is established to achieve the lowest rotor overcurrent after voltage dips, and the lowest amount of reactive power absorbed from the grid. In order to do that, the optimization objective function of the crowbar resistance can be constructed as follows.
- A. Objective Functions
1) Rotor Current Objective Function (FC): On the one hand, according to Equation (16), although the current vector does not reach the maximum value exactly at t = T/2, the current after half a period gives a good approximation of the maximum current [15] . On the other hand, after the crowbar is in operation, the rotor current damping time constant
PPT Slide
Lager Image
PPT Slide
Lager Image
Thus, the maximum rotor current I rmax after the grid voltage dip can be obtained by substituting
PPT Slide
Lager Image
and p = 1 in (16) for:
PPT Slide
Lager Image
2) Absorbed Reactive Power Objective Function (FR):
According to Equation (21), when the crowbar is in operation, the objective function of the reactive power that the DFIG absorbs from the grid can be expressed as follows:
PPT Slide
Lager Image
- B. Problem Constraints
The constraints are:
1) Rotor Overcurrent: As shown by Eqn. (23), the value of the crowbar resistance must be large enough to limit the DFIG rotor overcurrent and to make the time constant attenuate quickly. According to Equation (22), if
PPT Slide
Lager Image
is defined as the maximum damping time constant, the value of the crowbar resistance can be described as:
PPT Slide
Lager Image
2) Rotor Overvoltage: If U rmax is defined as the maximum rotor overvoltage during voltage sag faults, it is a function related to R cb .
PPT Slide
Lager Image
Thus, the maximum rotor overvoltage monotonically increases with an increase of the crowbar reactance after the crowbar is in operation. For example, for a 2 MW DFIG, Fig. 4 shows the relationship between U rmax and R cb .
PPT Slide
Lager Image
Relationship curve between the maximum rotor voltage and crowbar resistances.
In practical applications, in order to prevent the rotor overvoltage from reaching values beyond the safety limit of the rotor voltage, and to avoid damage to both the rotor side converter and the dc-capacitor, a safety margin coefficient λ is adopted, to which the upper limit of the crowbar reactance value is connected as follows:
PPT Slide
Lager Image
where U rlim is the rotor voltage limit.
3) DFIG Absorbing Reactive Power: According to Equation (24), for example, for a 2 MW DFIG Fig. 5 shows the relationship between F R and R cb . As can be seen, with the increase of R cb , the reactive power Q R that DFIG absorbs from the grid gradually decreases.
PPT Slide
Lager Image
Relationship curve between the absorbed reactive power and crowbar resistances.
During the LVRT process of an actual wind power plant, the reactive power absorbed by the DFIG should be as little as possible to benefit the quick recovery of the dropped grid voltage. Assuming that Q s is the reactive power which the DFIG should provide to the grid when the crowbar is in operation, according to the LVRT test requirements of national wind turbine generators, if the reactive power that the DFIG absorbs from the grid is less than β times Q s , then the objective function of the reactive power will meet the following constraint conditions:
PPT Slide
Lager Image
- C. Fuzzy Processing of Multi-Objective Optimization Model
During actual engineering design, many parameters are not the determined values. There is a transition interval between the absolute allowable values and the absolute non-allowable values of parameters. This produces the so-called fuzziness. However, the problem regarding the transition interval is not taken into account in common optimization design methods. In turn, the optimal optimization effect cannot be achieved in most cases. In this paper, considering the fuzzy factors of the crowbar resistance parameters involved in the LVRT, the fuzzy set theory is introduced into the design so that the optimization results are much more suitable for actual engineering applications. The fuzzy set theory generalizes the classical set theory to allow for partial membership with a smooth boundary. The degree of membership in a set is expressed by a number between 0 and 1, where 0 means entirely not in the set, 1 means completely in the set, and a number in between 0 and 1 means partially in the set [16] . The details of the fuzzy optimization algorithm are described below.
1) Fuzzy Processing of Constraint Conditions: The fuzzy processing of constraint conditions is performed to establish the membership functions for each allowable fuzzy interval. The distribution of the linear membership function is shown in Fig. 6 .
PPT Slide
Lager Image
Membership function distribution of constraint conditions.
From Fig. 6 , the constraint condition membership function can be written as follows [17] :
PPT Slide
Lager Image
PPT Slide
Lager Image
The upper and lower limits of the constraint conditions are determined by the method of the amplification coefficient, that is a = εb , d = δc . In general the coefficient of the upper limit is ε = 1.05~1.30 , and the coefficient of the lower limit is δ = 0.75~0.90. The membership function of the fuzzy constraint with the loosest constraint conditions can be obtained by processing the constraint conditions as membership functions.
In this section, the constraint conditions of the crowbar resistance values in Eqs. (25), (27), and (28) show fuzziness and need for fuzzing. A membership function with a negative slope is shown in Fig. 7 and is adopted as the fuzzy constraint set of the crowbar resistance value. When λ and β are less than 1, that is a converter without overcurrent and overvoltage and where the reactive power the DFIG absorbed from the grid is less than the one, which the DFIG should provide to the grid when the crowbar is in operation. In these conditions, the satisfaction of the crowbar resistance value is equal to 1. Otherwise, when λ and β are larger than 1, and are gradually increased within the specified safety margin, the satisfaction of the crowbar resistance value is dropped down from 1 to 0.
PPT Slide
Lager Image
Membership function with negative slope.
2)Fuzzy Processing of the Objective Functions:
Step 1 : Construct a fuzzy objective function for every sub-objective function.
During the fuzzy processing of the objective function, the extremum method is often employed to construct the fuzzy object set. Equations (23) and (24) are the objective functions of the crowbar resistance value when a grid fault takes place. These are then transferred into the fuzzy objective functions within the closed interval [0,1]. The upper limit
PPT Slide
Lager Image
and the lower limit
PPT Slide
Lager Image
of each constraint condition are estimated according to Equations (25), (27), and (28), as well as the coefficients of the safety margin λ and β . Finally, in accordance with [18] , the fuzzy objective function for every sub-objective function Fi (ith objective function Fi ) is:
PPT Slide
Lager Image
Step 2 : Construct a membership function for every sub-objective function fuzzy objective set.
The i th objective function Fi is constructed by a membership function μi defined by [19] .
PPT Slide
Lager Image
Step 3 : Fuzzy conditional extremum estimation
According to the level set method of the fuzzy optimization theory, a variable γ is used to separate the fuzzy sub-set into several ordinary sets Fγ . Then, the extremum of the objective function in Fγ is obtained. Finally, the extremum of the fuzzy condition in Fγ can be derived as:
PPT Slide
Lager Image
where x is an optimized design variable.
3) Multi-Objective Fuzzy Optimization Mathematical Model: Based on the above processing, the multi-objective optimized mathematical model described in this paper is applied to solve a single-objective fuzzy optimization problem. The mathematical description of this problem is as follows:
  • SolveX= {x1,x2,...xn,γ]T
  • maxγ
PPT Slide
Lager Image
where μi ( Fi ) and μi ( gi ) are the objective function and the membership function of the constraint conditions after the fuzzy processing.
4) Genetic Algorithm Implementation: The optimum solution of the fuzzy optimization problem as described in Equation (34) can be calculated by an ordinary optimization algorithm. However, this optimization method usually falls into a local optimization, preventing it from achieving a globally optimal solution. Since genetic algorithms have the advantages of a strong global search and of a high convergence rate, they are employed to solve the multi-objective fuzzy optimization model in this paper.
The use of genetic algorithms (GA), which are a class of evolutionary algorithms (EA), has become prominent in recent years [20] . First introduced in Holland in 1975, the GA is a model in which an initial population develops into a new generation, through procreation, mutation, and crossover processes among all the members. A GA offers several advantages over traditional optimization algorithms [21] - [23] .
  • 1. A GA operates only with the coding chromosomes, which are optimized in accordance with objective functions and proper constraints. Therefore, the need for initial system information to begin the optimization process is eliminated.
  • 2. Another advantage of a GA, and notably one of its essential distinctions from traditional optimization algorithms, is its ability to explore numerous regions of the search space, by iteratively analyzing multiple points of a population.
  • 3. Finally, in a GA only the best candidates among the members are chosen, leading to global minimum convergence.
When solving the crowbar resistances fuzzy optimization problem using a GA, the following steps should be followed:
Step 1 : Encoding
The application of the GA begins with parameter encoding. Each individual in the population is a candidate solution, the components of which comprise all the decision variables in the system. The decision variables in the case of the crowbar resistance fuzzy optimization include the stator voltage u 0 , the voltage dip ratio p, and the crowbar resistance value R cb . The solution variables are expressed as floating point numbers and integers. This representation is more efficient than binary coding, because it does not require converting the resulting variables into the binary type. Moreover, less computer memory is needed to store the population.
Step 2 : Evaluation
Each individual in the population is given a goodness value according to the Fitness Function (FF). This value is used to drive the evolution process. At the same time, each chromosome is assigned a fitness value, which determines both the individual’s capacity to adapt to the environment, and its estimated survival probability in the subsequent generation.
Step 3 : Selection
Those individuals with the highest fitness values are selected as potential parents, since they are more likely to survive. Consequently, the larger the chromosome’s fitness, the greater the probability of survival in the next generation. Therefore, a roulette wheel selection is adopted to select the most suitable individuals to be reproduced.
Step 4 : Crossover
A one-point crossover operation is used, in which one crossover point is required to create new offspring chromosomes. This is arbitrarily chosen after the random selection of two chromosome parents.
Step 5 : Mutation
The mutation operator’s final objective is to add diversity to the population. Therefore, the operator arbitrarily selects individuals and changes some of their characteristics. The mutation is achieved with a small probability after crossover, and it creates new individuals with information not included in previous generations.
According to the description given in section 3, the complete multi-objective fuzzy optimization process of the crowbar resistance value is shown in Fig. 8 .
PPT Slide
Lager Image
Multi-objective fuzzy genetic algorithm process.
IV. SIMULATION RESULTS
- A. Crowbar Resistance Fuzzy Optimization Design
The proposed multi-objective fuzzy optimization genetic algorithm has been applied to solve the crowbar resistances optimization problem for the DFIG LVRT. The optimization algorithm has been conducted on a 3 MW DFIG. The parameters of the generators are given in the Appendix. The DFIG experiences a voltage sag of 80% at the grid and it is sustained for a period of 625 ms.
According to Section 2.2, the transient stator and rotor inductances are:
PPT Slide
Lager Image
The maximum rotor current, calculated from Equation (23), is 5.2 p.u. The value of the crowbar resistance is chosen so that the rotor voltage limit U r-limit is 1.5 p.u, the stator and rotor current limit I r-limit is 1.5 p.u, and the reactive power absorbed by the DFIG from the grid is 0.5 p.u.
The best result of the GA was obtained with the following algorithm parameters:
Number of generations: 100; population size: 55; crossover probability: 0.85; mutation probability: 0.01.
After 100 generations it was found that all of the individuals had reached almost the same fitness value. This shows that GA has reached the optimal solution. Table I lists the optimal solutions of the crowbar resistance and the individual objective function values corresponding to each solution. The algorithm took 126 s to reach the optimal solution. The crowbar resistance design effect of both the common method and the multi-objective fuzzy genetic algorithm are compared. This comparison is given in Table II . It is obvious that the multi-objective fuzzy genetic algorithm obtained the best compromise effect.
SINGLE OBJECTIVE OPTIMAL SOLUTION NOMINAL PARAMETERS
PPT Slide
Lager Image
SINGLE OBJECTIVE OPTIMAL SOLUTION NOMINAL PARAMETERS
COMMON DESIGN AND MULTI-OBJECTIVE GENETIC FUZZY ALGORITHM COMPARED
PPT Slide
Lager Image
COMMON DESIGN AND MULTI-OBJECTIVE GENETIC FUZZY ALGORITHM COMPARED
- B. The Optimization Results Validation
The optimized value of the crowbar resistance in the above is applied to a crowbar circuit under the conditions of DFIG grid faults. The transient performances are compared with the crowbar resistance which is designed by the common method. The common method selects a value of the crowbar resistance according to reference [15] , the crowbar resistance is selected based on the permissible value of the rotor overvoltage. A 3 MW DFIG connected to a power system network was simulated using MATLAB/SIMULINK software. The simulation model is shown in Fig. 9 . The DFIG parameters are given in the Appendix. The behavior of the DFIG during a voltage sag of 80% at the generator bus for 625 ms is simulated at t = 4s. At the moment of the voltage dip the generator speed is 0.8 p.u, the wind speed is 12m/s, and the generator is operated at a unity power factor. When the grid voltage sags because the rotor overcurrent crowbar protection is in operation, the crowbar action lasts for around 70 ms. After the crowbar is cut off, the rotor converter is immediately restarted and provides reactive support to the grid. To prevent converter overcurrent, after the rotor side converter restarts, no active power is generated.
PPT Slide
Lager Image
Simulation model of the system.
The stator flux during the voltage sag is shown in Fig. 10 . The stator flux has the same change trend for the crowbar resistance regardless of whether it was designed by the common method or by the optimization method. However, for the optimization method, at the beginning and end of the voltage sag the transient peak is slightly increased. This is due to the fact that when the rotor current is not equal to zero, the dynamic response of the stator flux is affected by the flux coupling between the stator and the rotor. The smaller the value of the crowbar resistance, the better of the stator flux attenuation. The rotor flux during the voltage sag is shown in Fig. 11 . When compared with the common design method, the optimization method that designed the crowbar resistance can make the rotor flux at the beginning of the voltage sag have a smaller peak and it can make the end of the voltage sag have a faster recovery speed. The rotor flux oscillations are attenuated. The reason for this is that from Equation (23), after the crowbar is in operation, the rotor current damping time constant is increased, which accelerates the attenuation of the rotor flux.
PPT Slide
Lager Image
Stator flux.
PPT Slide
Lager Image
Rotor flux.
The stator and rotor current responses are shown in Figs. 12 and 13 . As can be noticed, if the common design method is employed, the stator and rotor overcurrents reach maximum values of 1.3 and 1.15 p.u. during the fault, respectively. It can also be seen that the stator overcurrents are near the constraint boundary. At the same time, by using the proposed optimization algorithm, the maximum values of stator and rotor overcurrents decrease to 0.95 and 0.75 p.u., respectively. These are far beyond the safety limit of the converter. In addition, the transient time constant is smaller. The rotor voltage response is shown in Fig. 14 . The multi-objective fuzzy optimization algorithm obtains the compromise solution of the crowbar resistance. Furthermore, the transient peak is slightly increased. However, it is less than the safety limit of the rotor voltage.
PPT Slide
Lager Image
Stator current.
PPT Slide
Lager Image
Rotor current.
PPT Slide
Lager Image
Rotor voltage.
The active and reactive power responses are shown in Figs. 15 and 16 . The stator active power due to the occurrence of the voltage sag decreases to nearly zero. Then, after clearing the fault it starts to oscillate until it reaches its steady state value. It is evident that when compared with the common design method, the proposed optimization algorithm has a smaller transient peak during the fault. In addition, the optimized value of the crowbar resistance makes the DFIG absorb comparatively less reactive power from the grid. This is a benefit caused by the quick recovery of the dropped terminal voltage. However, for the common design method, the reactive power magnitude is increased to almost 0.6 p.u, which is beyond the design target.
PPT Slide
Lager Image
Active power.
PPT Slide
Lager Image
Reactive power.
The dynamic behavior of the electromagnetic torque is shown in Fig. 17 . It is clear that if the common design method is employed, the electromagnetic torque magnitude is increased to 1.3 p.u. This tends to damage the mechanical parts in practical applications.
PPT Slide
Lager Image
Electromagnetic torque.
The rotor speed is shown in Fig. 18 . During the voltage sag the active power decreases and the rotor speed increases. The speed peak is decreased in the optimization design. The recovery process of the DFIG terminal voltage is shown in Fig. 19 . If the common design is employed, the terminal voltage experiences a dip below 80%, the recovery voltage transient peak is increased to 1.17 p.u., and the recovery process is unduly delayed (at t = 4.85s). The proposed optimization algorithm leads to a rapid recovery of the terminal voltage of the DFIG to its nominal voltage, at t = 4.72s.
PPT Slide
Lager Image
Rotor speed.
PPT Slide
Lager Image
Grid terminal voltage.
V. CONCLUSION
This study investigated the multi-objective fuzzy optimization of crowbar resistances for the LVRT of DFIG under voltage dips with a genetic algorithm. To effectively improve the security of the LVRT problem and to provide designers with diverse alternative options, a multi-objective fuzzy optimization model with technical and operational constraints has been constructed to minimize rotor overcurrent and to simultaneously reduce the reactive power that DFIGs absorb from the grid. A genetic algorithm has been successfully applied to the multi-objective fuzzy optimization problem and it derived a satisfactory decision by finding optimal solutions. The proposed algorithm has been applied to a 3 MW DFIG LVRT. A comparison between the results obtained with the crowbar resistance designed by the proposed optimization algorithm and by the common design method proves that the optimized crowbar resistance makes the LVRT have a higher safety performance and a faster recovery speed from grid faults. The multi-objective fuzzy optimization emphasizes non-dominated solutions and simultaneously maintains diversity in the non-dominated solutions. In the future, the proposed approach can be widely applied to solve safety-constrained optimum problems in power systems.
Acknowledgements
This work was supported by National Nature Science Foundation of China(51207002), the Industrial Science and Technology Research of Shaanxi Province(2013K07-12 and 2012K06-46), the Scientific Research Program Funded by Shaanxi Provincial Education Department(2013JK1017), the Technology transfer to promote engineering project of Xi 'an bureau of science and technology (CXY1347(4)) and the key research project Funded by Baoji University of Arts and Sciences(ZK11068).
BIO
Wenjuan Zhang was born in Shaanxi, China. She received her B.S. degree in Automation Engineering, and her M.S. and Ph.D. degrees in Power Electronics and Electrical Drives from the Xi’an University of Technology, Shaanxi, China, in 2004, 2007 and 2011, respectively. From 2007 to 2011, she was with Wind Power Generation Technology, Xi’an, China. In 2012, she joined the Baoji University of Arts and Sciences, Baoji, China, as an Associate Professor. Her current research interests include wind power generation, AC-AC converters, and PWM converter/inverter systems.
Haomiao Ma was born in Shaanxi, China. She received her B.S. degree in Electrical Engineering from the Hefei Technological University, Hefei, China, in 2000, her M.S. degree in Electrical Engineering from Xi’an Jiaotong University, Xi’an, China, in 2007, and her Ph.D. degree in Electrical Engineering from the Xi’an University of Technology, Shaanxi, China, in 2013. Since 2000, she has been with the School of Computer Science, Shaanxi Normal University, Xi’an, China, where she is presently working as an Associate Professor. Her current research interests include electronics, energy systems and signals.
Junli Zhang was born in Shaanxi, China. He received his B.S. degree in Electrical Engineering from the Baoji University of Arts and Sciences, Baoji, China, in 2001, and his M.S. degree in Electrical Engineering from the Xi’an University of Technology, Shaanxi, China, in 2007. Since 2001, he has been with the Department of Electronic and Electrical Engineering, Baoji University of Arts and Sciences. He was a Teaching Assistant from 2001 to 2006, an Associate Professor from 2006 to 2011, and he is presently working as an Assistant Professor. His current research interests include new power electronic equipment and its application.
Lingling Chen was born in Sichuan, China. She received her M.S. degree in Electrical Engineering from the Xi'an University of Technology, Xi'an, China, in 2005, and her M.S. degree in Power Electronics and Electrical Drives from the Xi’an University of Technology, Shaanxi, China, in 2007. In 2008, she joined 94816 troop with the People's Liberation Army, China. Her current research interests include PWM converter/inverter systems.
Yang Qu born in Shaanxi, China. He is presently working towards his B.S. degree in Automation Engineering at the Baoji University of Arts and Sciences, Baoji, China.
References
Li H. , Yang C. , Hu Y. G. , Zhao B. , Zhao M. , Chen Z. 2014 “Fault-tolerant control for current sensors of doubly fed induction generators based on an improved fault detection method,” Measurement 47 (2) 929 - 937    DOI : 10.1016/j.measurement.2013.10.021
Rahimi M. , Parniani M. 2010 “Grid-fault ride-through analysis and control of wind turbines with doubly fed induction generators,” Electric Power Systems Research 80 (2) 184 - 195    DOI : 10.1016/j.epsr.2009.08.019
Alberdi M. , Amundarain M. , Garrido A. , Garrido I. 2012 “Neural control for voltage dips ride-through of oscillating water column-based wave energy converter equipped with doubly-fed induction generator,” Renewable Energy 48 16 - 26    DOI : 10.1016/j.renene.2012.04.014
Vinothkumar K. , Selvan M. P. 2011 “Novel scheme for enhancement of fault ride-through capability of doubly fed induction generator based wind farms,” Energy Conversion and Management 52 (7) 2651 - 2658    DOI : 10.1016/j.enconman.2011.01.003
Ouyang J. , Xiong X. 2014 “Research on short-circuit current of doubly fed induction generator under non-deep voltage drop,” Electric Power Systems Research 107 158 - 166    DOI : 10.1016/j.epsr.2013.10.008
Mohammadi J. , Afsharnia S. , Vaez-Zadeh S. 2014 “Efficient fault-ride-through control strategy of DFIG-based wind turbines during the grid faults,” Energy Conversion and Management 78 88 - 95    DOI : 10.1016/j.enconman.2013.10.029
Niiranen J. “Voltage dip ride through of doubly-fed generator equipped with active crowbar,” in the Nordic Wind Power Conference 2004 1501 - 1507
Morren J. , Haan S. W. H. 2012 “Ride-through of wind turbines with doubly-fed induction generator during a voltage dip,” IEEE Trans. Energy Convers. 47 (2) 435 - 441
Sulla F. , Svensson J. , Samuelsson O. 2011 “Symmetrical and unsymmetrical short-circuit current of squirrel-cage and doubly-fed induction generators,” Electric Power Systems Research 81 (7) 1610 - 1618    DOI : 10.1016/j.epsr.2011.03.016
Honglin Z. , Geng Y. “Short circuit current characteristics of doubly fed induction generator with crowbar protection under different voltage dips,” in Proc. the CSEE Vol. 47, No. 2, pp. 184-191, Jan. 2009. (in Chinese)
Mohseni M. , Masoum M. A. S. , Islam S. M. 2011 “Low and high voltage ride-through of DFIG wind turbines using hybrid current controlled converters,” Electric Power Systems Research 81 (7) 1456 - 1465    DOI : 10.1016/j.epsr.2011.02.010
Vinothkumar K. , Selvan M. P. 2011 “Novel scheme for enhancement of fault ride-through capability of doubly fed induction generator based wind farms,” Energy Conversion and Management 52 (7) 2651 - 2658    DOI : 10.1016/j.enconman.2011.01.003
Lopez J. , Gubia E. , Sanchis P. 2008 “Wind turbines based on doubly fed induction generator under asymmetrical voltage dips,” IEEE Trans. Energy Convers. 23 (1) 321 - 329    DOI : 10.1109/TEC.2007.914317
Lopez J. , Sanchis P. , Roboam X. , Marroyo L. 2007 “Dynamic behavior of the doubly fed induction generator during three-phase voltage dips,” IEEE Trans. Energy Convers. 22 (3) 9 - 17    DOI : 10.1109/TEC.2006.878241
Morren J. , de Haan S. W. H. 2007 “Short-circuit current of wind turbines with doubly fed induction generator,” IEEE Trans. Energy Convers. 22 (1) 174 - 180    DOI : 10.1109/TEC.2006.889615
Banu R. N. , Devaraj D. 2012 “Multi-objective GA with fuzzy decision making for security enhancement in power system,” Applied Soft Computing 12 (9) 2756 - 2764    DOI : 10.1016/j.asoc.2012.03.057
Zhang W. , Liu Y. 2008 “Multi-objective reactive power and voltage control based on fuzzy optimization strategy and fuzzy adaptive particle swarm,” Electr Power Energy System 30 (9) 525 - 532    DOI : 10.1016/j.ijepes.2008.04.005
Perera A. T. D. , Attalage R. A. , Perera K. K. C. .K. 2013 “A hybrid tool to combine multi-objective optimization and multi-criterion decision making in designing standalone hybrid energy systems,” Applied Energy 107 412 - 425    DOI : 10.1016/j.apenergy.2013.02.049
Lee K. Y. , Yome A. S. , Park J. H. 1998 “Adaptive Hopfield neural networks for economic load dispatch,” IEEE Trans. Power Syst 13 (2) 519 - 526    DOI : 10.1109/59.667377
Cheng S. , Chen M.-Y. 2014 “Multi-objective reactive power optimization strategy for distribution system with penetration of distributed generation,” International Journal of Electrical Power & Energy Systems 62 221 - 228    DOI : 10.1016/j.ijepes.2014.04.040
Khorram E. , Zarei H. 2009 “Multi-objective optimization problems with Fuzzy relation equation constraints regarding max-average composition,” Mathematical and Computer Modelling 49 (5-6) 856 - 867    DOI : 10.1016/j.mcm.2008.10.018
Thapara A. , Pandey D. , Gaur S. K. 2012 “Satisficing solutions of multi-objective fuzzy optimization problems using genetic algorithm,” Applied Soft Computing 12 (8) 2178 - 2187    DOI : 10.1016/j.asoc.2012.03.002
Wang L. , Singh C. 2008 “Balancing risk and cost in fuzzy economic dispatch including wind power penetration based on particle swarm optimization,” Electric Power Systems Research 78 (8) 1361 - 1368    DOI : 10.1016/j.epsr.2007.12.005