# Control of DFIG Wind Turbines Based on Indirect Matrix Converters in Short Circuit Mode to Improve the LVRT Capability

Authors: Ahmad Khajeh; Reza Ghazi
Journal: Advances in Power Electronics
 (2013)

---

## Abstract

Nowadays, the doubly-fed induction generators (DFIGs) based wind turbines (WTs) are the dominant type of WTs. Traditionally, the back-to-back converters are used to excite the rotor circuit of DFIG. In this paper, an indirect matrix converter (IMC) is proposed to control the generator. Compared with back-to-back converters, IMCs have numerous advantages such as higher level of robustness, reliability, and reduced size and weight due to the absence of bulky electrolytic capacitor. According to the recent grid codes wind turbines must have low voltage ride-through (LVRT) capability. In this paper a new crowbar system is proposed so that along with the control system it protects the IMC from large fault currents and supports the grid voltage dips during grid faults. This crowbar system is provided using the existing converter switches to establish a short circuit mode without any extra circuitry. Even in severe fault conditions, the duration of short circuit mode is quite small so the control system will be activated shortly after the fault to inject reactive power as required by new LVRT standards.  Therefore, the new LVRT standards are well satisfied without any extra costs. PSIM simulation results confirm the efficiency of the proposed method.


---

## Body

## 1. Introduction

Wind energy has been found as the fastest growing renewable power generation resources over the world. Government support and advancement in technology of WT have accelerated this rapid growth. Based on estimation the wind power capacity becomes double every three years [1].Due to their superior characteristics, most of the grid-connected WTs operate at a variable speed. Among the different variable speed types, the DFIG is the most promising one. The stator winding of DFIG is directly connected to the grid, while the rotor winding is connected to the grid through an ac-ac power electronic converter having bidirectional switches. Traditionally the back-to-back converters are used to excite the rotors of DFIGs. The presence of the dc-link capacitor in this arrangement is a serious drawback as it increases the costs and reduces the overall lifetime of the system and also makes the system bulky [2].In this paper, the back-to-back converter arrangement is replaced by an indirect matrix converter (IMC) to control the generator as shown in Figure1. The main advantages of a matrix converter are robustness, reliability with less size and weight due to the absence of the bulky electrolytic capacitor, controllable input power factor, nearly sinusoidal input current, and output voltage with only switching frequency harmonics, along with bidirectional power flow. The direct matrix converter (DMC) encounters the commutation problems requiring a complex control circuitry. While in IMC all switches at the line side will turn on and off at zero current, so the commutation problems are eliminated [3]. Therefore, the IMCs are the most promising devices for wind energy applications regarding their robustness, reliability, and reduced size.Figure 1
Wind Turbine based on the DFIG excited by indirect matrix converter.In the past, based on most grid codes, wind turbines were allowed to be disconnected from the grid during the grid disturbances and the abnormal voltage reduction. With the increased capacity of wind powers in the power systems over the years, the sudden loss of wind turbines during the grid faults resulting from turbines disconnection could generate control problems of frequency and voltage in the system so leading to the voltage collapse in worst cases. The increased penetration of wind energy into the power system over the last decade has therefore led to a serious concern about its influence on the dynamic behavior of the power system. To handle this issue the grid codes are revised by system operators in several countries [4]. According to the new grid codes the wind turbines must be remain connected to the network on the occurrence of grid faults and tolerate the resulting voltage drops. This feature is known as the low voltage ride-through (LVRT) capability of a power plant. The LVRT requirement defines a specific voltage profile with limiting curves depending on the grid voltage dip depth. Moreover, in some grid codes such as Germany, wind turbines should also inject reactive power to help the grid voltage recovery.Because of direct connection of stator winding to the grid and the small rating of power electronic converter in the rotor side, DFIG based wind turbine is very sensitive to grid disturbances, especially to voltage dips during grid faults. Grid faults, even those occurring far from the wind turbines, can cause voltage drops at the point of wind turbine connection. The abrupt drop of the grid voltage will make an increase in the current of stator windings of the DFIG. Because of the magnetic coupling between stator and rotor, this current will also flow in the rotor circuit and the power converter. At present, the back-to-back converter is the most frequently used power electronic converter in the wind turbines industry. Therefore, the most of research works have been done to enhance the LVRT ability of DFIG utilizing this type of converters. So the grid faults cause over current in the rotor windings and over voltage in the dc-link capacitor; hence, the proper protection should be provided to safekeep the converter. Various solutions have been proposed to solve the problem. Some of them are briefly reviewed as follows.The most popular and reliable method is based on the use of a protective circuit known as crowbar. The crowbar consists of resistors connected to the rotor windings by means of electronic switches. In active crowbar the switches are controllable such as IGBTs whereas in passive one they are antiparallel thyristors. The active crowbar allows opening the circuit when the currents reach safe region. The passive crowbar allows closing the circuit but not opening it until the crowbar currents reach zero. Therefore, according to the new grid codes, the passive crowbar is not appropriate for the modern wind turbines. When a fault occurs, the rotor windings get connected to the crowbar while the converter is tripped. Several papers have discussed the implementation and control of crowbar [5–12]. Different types of crowbars realization are shown in Figure 2. The crowbar can effectively protect the rotor converter under serious grid faults. But its main drawback is that, when the crowbar is activated, the rotor converter cannot control the active and reactive power during crowbar activation, so the DFIG operates as a cage induction generator, absorbing extra reactive power from the grid and deteriorates the grid voltage profile.Different types of crowbar realization.
(a)
(b)
(c)Another approach is to modify the rotor-side converter control system to limit the fault currents without using any extra hardware circuit [13–21]. In this method at serious grid faults, the required voltage in the rotor-side converter becomes too large and exceeds its voltage capability. Therefore, the rotor current greatly increases thus leading to converter destruction. This method is beneficial as it does not require any additional hardware but it is suitable for moderate faults.Another solution for LVRT is to use an additional energy storage system (ESS) connected to the dc link [22–26]. The ESS can balance the extra power that goes through the rotor circuit during a grid fault transient, but it also requires a rotor side converter with higher current ratings having extra cost and more system complexity.Utilizing matrix converters in DFIG wind turbines is still in research phase. Most of the works have dealt with their normal operation [27–30]. The control of a DFIG using an indirect matrix converter for dynamic performance evaluation under normal grid condition is presented in [27]. Sigma-delta modulator to control IMC switches of DFIG based WT is proposed in [28]. By using this modulation method, torque pulsations and harmonic content of currents are reduced. Therefore, the power quality of WT is improved. The capability of the input converter to generate different virtual dc-link voltage levels is also proposed in [29]. This method reduces both the commutation losses of output converter and common mode voltage. Due to small voltage required around the synchronous speed, this strategy is suitable in DFIG applications.Dynamic performance of DFIG using two-stage IMC under voltage dip condition has been discussed in [30] without providing any protection. As simulation results show for 80% terminal voltage sag, the stator current reaches 2.5 pu. Without any protection method, this high current could damage the IGBT switches of the IMC. According to new LVRT standards, the reactive power injection to the grid is mandatory to improve the voltage recovery. Results of [30] indicate that, during the fault, about 0.7 pu reactive power is absorbed from the grid which makes the voltage deeper. Therefore, the existing controller which is designed for normal operating condition is ineffective during fault conditions.To solve these drawbacks and satisfy the new LVRT standards, in present paper a new method is proposed to protect the IMC with DFIG following the grid faults. In this method a new short circuit mode is proposed to control the short circuit currents of the DFIG wind turbines based on IMC. In addition to protect the IMC against large fault currents, the proposed method also makes it possible to inject reactive power during grid faults to improve voltage recovery, so the new LVRT standards are well satisfied. The PSIM simulation results confirm the efficiency of the proposed method.This paper is organized as follows. In Section2, LVRT requirement of a German transmission system operator is discussed. Modeling and control of the DFIG based wind turbine is presented in Section 3. The vector control of the DFIG is briefly discussed in Section 4. In Section 5, the proposed method is presented and analyzed. The simulation results are provided in Section 6. Finally, Section 7 concludes the paper.


## 2. LVRT Requirement

The German transmission system operator E.ON Netz was the first power system operator, which introduced grid codes for wind turbines and is followed now by many other network operators in several countries. E.ON introduces a voltage profile and the limiting curves and regions defining the LVRT requirement as shown in Figure3. Accordingly, wind turbines must stay connected even when the voltage at the point of common coupling (PCC) drops to zero. The 150 ms delay shown in Figure 3 accounts for the normal operating time of protection relays. Three-phase short circuits or fault-related symmetrical voltage dips must not lead to instability above the limit line 1 in Figure 3 or disconnection of the wind turbines from the grid. Within the shaded area and above limit line 2 in Figure 3, LVRT is also required but in case of instability a short term interruption (STI) is allowed. Below the limit line 2 in Figure 3, no LVRT is required and STI from the grid is always permissible. Here, resynchronization times of more than 2 seconds and an active power increase following fault clearance of less than 10% of the rated power per second are also possible. According to the E.ON Netz grid code, wind turbines have to provide a mandatory voltage support during voltage dips. Wind turbines have to supply 1 pu reactive current when the voltage falls below 0.5 pu.Figure 3
E.ON Netz LVRT requirement.


## 3. Modeling and Control of DFIG Based Wind Turbine

In this section, dynamic model of DFIG based wind turbine is provided. The stator winding of DFIG is directly connected to the grid, and the rotor winding is coupled via an IMC. The IMC must handle the slip power, that is, about 30% of the rated power of wind turbine. The speed range of the generator is typically ±30% of synchronous speed, thus providing flexibility to operate in both sub- and supersynchronous modes, depending on the wind conditions. The inverter of the IMC is used to control the active and reactive power of the DFIG. The rectifier of the IMC is often operated at unity power factor. Depending on the rotor speed, the IMC will either absorb power (subsynchronous) from grid or inject power (supersynchronous) to the grid. Therefore, the IMC must have the ability of bidirectional power flow to the network.


### 3.1. Turbine Model

The mechanical power extracted by a wind turbine from the wind is expressed by(1)P=12A·ρ·Cp(λ,β)·vω3,
where A is the area covered by the rotor blades, ρ is the air density, Cp is the power coefficient, representing the amount of power the turbine can extract, and vω is the wind speed. The power available in the wind cannot be extracted completely. Theoretically the maximum captured power is 59% of the power available in the wind. The power coefficient is a function of the tip-speed ratio λ and the pitch angle of the rotor blades β. The tip-speed ratio is defined by
(2)λ=RΩtvω,
where R is the radius of the rotor blades and Ωt is the angular speed of the blades. For each pitch angle of the rotor blades, there is an optimum tip-speed ratio λopt for which Cp(λopt,β) takes a maximum value. Therefore, for wind speeds below rated value, the existing power is below the rated power, so power coefficient is optimized by adjusting the tip-speed ratio to achieve the maximum power coefficient. The pitch angle is kept constant in this region which is called the maximum power point tracking (MPPT) of wind turbines. On contrary, at wind speeds above rated value the extracted wind power has to be limited by means of blade pitching. In this paper the MPPT is accomplished through a look-up table. The MPPT table ensures the plant operation at maximum power for partial loads and specifies that the rated speed is achieved at rated power. According to [4] coordination between the control of the mechanical and the electrical system exists.


### 3.2. DFIG Model

For DFIG modeling a fifth-order dynamic model of the DFIG is used in this paper [31]. The model in a two-axis d-q synchronous reference frame is given by
(3)vsdq=Rsisdq+dψsdqdt+jωsψsdq,vrdq=Rrirdq+dψrdqdt+j(ωs-ωr)ψrdq,Tem=p2(ψrqird-ψrdirq),(4)ψsdq=Lsisdq+Lmirdq,ψrdq=Lmisdq+Lrirdq,Ls=Lls+Lm,Lr=Llr+Lm.In these equations,Rs, Rr, Ls, Lr, Lls, and Llr are the resistors and inductors of the stator and rotor windings, Lm is the magnetizing inductance,vsdq, vrdq, isdq, irdq, ψsdq, and ψrdq are the space vectors of the stator and rotor voltages, currents, and fluxes, respectively, ωs is the synchronous speed of the generator, ωr is the electrical speed of the rotor, and p is the number of poles.


### 3.3.  Indirect Matrix Converter

Due to the numerous advantages of matrix converters over back-to-back converters, in this paper the IMC is used to control the DFIG. An IMC consists of a rectification part on the input side and an inversion part on the output side, connected via fictitious dc link as shown in Figure4. For purposes of analysis, we can assume that the switching frequency is far greater than the fundamental frequency of both the input voltage and output current. Thus during each switching cycle, it is assumed that both input voltage and output current are constant.Figure 4
Indirect matrix converter topology.The rectifier has six bidirectional switches with the ability of conducting current and blocking voltage in both directions. The rectifier side objective is to achieve maximum positive voltages at the fictitious dc-link and sinusoidal input currents. Usually the grid side converter of DFIG exchanges zero reactive power at the grid point. In order to obtain maximum dc-link voltage, the input phase voltage which has the highest absolute value is connected to the positive or negative rail of the dc link at 60 degree intervals depending on its polarity. To achieve sinusoidal current and unity power factor at the input side, regardless of the load type, the other two phase voltages are modulated so that the reference current space vector be in phase with the voltage space vector. Space vectors of the input current are shown in Figure5. Assume the input voltages are
(5)va=Vmcosθa=Vmcos(ωit),vb=Vmcosθb=Vmcos(ωit-2π3),vc=Vmcosθc=Vmcos(ωit+2π3),
where Vm is the maximum of the input phase voltage and ωi is the angular frequency. In order to achieve unity power factor at the input side, the input currents must be in phase with the input voltages. Therefore, input currents are as follows:
(6)Ia=Imcos(ωit),Ib=Imcos(ωit-2π3),Ic=Imcos(ωit+2π3),
where Im is the maximum of the input current. On the rectifier side always two switches, one from top and another from the bottom, are ON and the others are OFF.Figure 5
Space vectors of the input current.To calculate the duty cycle of each switch, first the angle of the input voltage space vector is obtained. This is the reference angle of the current space vector to obtain the unity power factor. Regarding the angle of current space vector, its corresponding sector is determined from Figure5. Then the current reference vector can be built using two adjacent vectors of each section. For example, when the reference current space vector is located in the first section, two adjacent vectors are “ab” and “ac.” In this section, voltage of phase “a” has the highest absolute value. Therefore, in 60 degree duration of this section, the top switch of phase “a” is ON and the bottom switches of phases “b” and “c” get modulated. Vector “ab” means that the top switch of phase “a” and the bottom switch of phase “b” are ON, and “ac” means that the top switch of phase “a” and the bottom switch of phase “c” are ON and so on.So in one period of switching frequency in this section, input currents are as follows:(7)ia=(dab+dac)idc,ib=-dabidc,ic=-dacidc,
where dab and dac are the duty cycles of I1(ab) and I2(ac), respectively, and idc is the dc-link current. As only the active vectors are used, the following relationships hold:
(8)dab+dac=1⟹ia=idc.Therefore, the duty cycles of active vectors in section one are given by(9)dab=-ibia,dac=-icia.The duty cycles in the other sections are obtained similarly. Based on this modulation method in the rectifier stage, the average fictitious dc-link voltage in each period is(10)Vdc=3·Vm2·|cos(θin)|,|cos(θin)|=max(|cos(θa)|,|cos(θb)|,|cos(θc)|),
where θin is the angle of the input voltage space vector. On the output stage, the space vector modulation (SVM) is used to generate the required rotor voltage space vector by currents controllers. Inverter voltage space vectors along with the reference vector are shown in Figure 6.Figure 6
Space vectors of the output voltage.Duty cycles of active vectors on the inverter stage are as follows:(11)dα=k·sin(60∘-θi),dβ=k·sin(θi),
where k is the modulation index and is given by
(12)k=2·|cos(θin)|·V^r3·Vm,
where Vr is the absolute of the reference voltage. Considering modulation of both the input and output sides, the final duty cycles of active vectors are
(13)dαγ=dγ·dα,dβγ=dγ·dβ,dαδ=dδ·dα,dβδ=dδ·dβ,
where dγ is the first vector of the rectifier section and dδ is the second vector of the rectifier section. Finally, the duty cycles of the zero vectors are calculated as
(14)d0=1-(dγ+dδ)·(dα+dβ),d0γ=dγ(dγ+dδ)d0,d0δ=dδ(dγ+dδ)d0.The overall switching scheme in one period is depicted in Figure7. Fictitious dc-link voltage for simulated IMC is shown in Figure 8.Figure 7
Switching scheme in one period.Figure 8
Fictitious dc-link voltage.


## 3.1. Turbine Model

The mechanical power extracted by a wind turbine from the wind is expressed by(1)P=12A·ρ·Cp(λ,β)·vω3,
where A is the area covered by the rotor blades, ρ is the air density, Cp is the power coefficient, representing the amount of power the turbine can extract, and vω is the wind speed. The power available in the wind cannot be extracted completely. Theoretically the maximum captured power is 59% of the power available in the wind. The power coefficient is a function of the tip-speed ratio λ and the pitch angle of the rotor blades β. The tip-speed ratio is defined by
(2)λ=RΩtvω,
where R is the radius of the rotor blades and Ωt is the angular speed of the blades. For each pitch angle of the rotor blades, there is an optimum tip-speed ratio λopt for which Cp(λopt,β) takes a maximum value. Therefore, for wind speeds below rated value, the existing power is below the rated power, so power coefficient is optimized by adjusting the tip-speed ratio to achieve the maximum power coefficient. The pitch angle is kept constant in this region which is called the maximum power point tracking (MPPT) of wind turbines. On contrary, at wind speeds above rated value the extracted wind power has to be limited by means of blade pitching. In this paper the MPPT is accomplished through a look-up table. The MPPT table ensures the plant operation at maximum power for partial loads and specifies that the rated speed is achieved at rated power. According to [4] coordination between the control of the mechanical and the electrical system exists.


## 3.2. DFIG Model

For DFIG modeling a fifth-order dynamic model of the DFIG is used in this paper [31]. The model in a two-axis d-q synchronous reference frame is given by
(3)vsdq=Rsisdq+dψsdqdt+jωsψsdq,vrdq=Rrirdq+dψrdqdt+j(ωs-ωr)ψrdq,Tem=p2(ψrqird-ψrdirq),(4)ψsdq=Lsisdq+Lmirdq,ψrdq=Lmisdq+Lrirdq,Ls=Lls+Lm,Lr=Llr+Lm.In these equations,Rs, Rr, Ls, Lr, Lls, and Llr are the resistors and inductors of the stator and rotor windings, Lm is the magnetizing inductance,vsdq, vrdq, isdq, irdq, ψsdq, and ψrdq are the space vectors of the stator and rotor voltages, currents, and fluxes, respectively, ωs is the synchronous speed of the generator, ωr is the electrical speed of the rotor, and p is the number of poles.


## 3.3.  Indirect Matrix Converter

Due to the numerous advantages of matrix converters over back-to-back converters, in this paper the IMC is used to control the DFIG. An IMC consists of a rectification part on the input side and an inversion part on the output side, connected via fictitious dc link as shown in Figure4. For purposes of analysis, we can assume that the switching frequency is far greater than the fundamental frequency of both the input voltage and output current. Thus during each switching cycle, it is assumed that both input voltage and output current are constant.Figure 4
Indirect matrix converter topology.The rectifier has six bidirectional switches with the ability of conducting current and blocking voltage in both directions. The rectifier side objective is to achieve maximum positive voltages at the fictitious dc-link and sinusoidal input currents. Usually the grid side converter of DFIG exchanges zero reactive power at the grid point. In order to obtain maximum dc-link voltage, the input phase voltage which has the highest absolute value is connected to the positive or negative rail of the dc link at 60 degree intervals depending on its polarity. To achieve sinusoidal current and unity power factor at the input side, regardless of the load type, the other two phase voltages are modulated so that the reference current space vector be in phase with the voltage space vector. Space vectors of the input current are shown in Figure5. Assume the input voltages are
(5)va=Vmcosθa=Vmcos(ωit),vb=Vmcosθb=Vmcos(ωit-2π3),vc=Vmcosθc=Vmcos(ωit+2π3),
where Vm is the maximum of the input phase voltage and ωi is the angular frequency. In order to achieve unity power factor at the input side, the input currents must be in phase with the input voltages. Therefore, input currents are as follows:
(6)Ia=Imcos(ωit),Ib=Imcos(ωit-2π3),Ic=Imcos(ωit+2π3),
where Im is the maximum of the input current. On the rectifier side always two switches, one from top and another from the bottom, are ON and the others are OFF.Figure 5
Space vectors of the input current.To calculate the duty cycle of each switch, first the angle of the input voltage space vector is obtained. This is the reference angle of the current space vector to obtain the unity power factor. Regarding the angle of current space vector, its corresponding sector is determined from Figure5. Then the current reference vector can be built using two adjacent vectors of each section. For example, when the reference current space vector is located in the first section, two adjacent vectors are “ab” and “ac.” In this section, voltage of phase “a” has the highest absolute value. Therefore, in 60 degree duration of this section, the top switch of phase “a” is ON and the bottom switches of phases “b” and “c” get modulated. Vector “ab” means that the top switch of phase “a” and the bottom switch of phase “b” are ON, and “ac” means that the top switch of phase “a” and the bottom switch of phase “c” are ON and so on.So in one period of switching frequency in this section, input currents are as follows:(7)ia=(dab+dac)idc,ib=-dabidc,ic=-dacidc,
where dab and dac are the duty cycles of I1(ab) and I2(ac), respectively, and idc is the dc-link current. As only the active vectors are used, the following relationships hold:
(8)dab+dac=1⟹ia=idc.Therefore, the duty cycles of active vectors in section one are given by(9)dab=-ibia,dac=-icia.The duty cycles in the other sections are obtained similarly. Based on this modulation method in the rectifier stage, the average fictitious dc-link voltage in each period is(10)Vdc=3·Vm2·|cos(θin)|,|cos(θin)|=max(|cos(θa)|,|cos(θb)|,|cos(θc)|),
where θin is the angle of the input voltage space vector. On the output stage, the space vector modulation (SVM) is used to generate the required rotor voltage space vector by currents controllers. Inverter voltage space vectors along with the reference vector are shown in Figure 6.Figure 6
Space vectors of the output voltage.Duty cycles of active vectors on the inverter stage are as follows:(11)dα=k·sin(60∘-θi),dβ=k·sin(θi),
where k is the modulation index and is given by
(12)k=2·|cos(θin)|·V^r3·Vm,
where Vr is the absolute of the reference voltage. Considering modulation of both the input and output sides, the final duty cycles of active vectors are
(13)dαγ=dγ·dα,dβγ=dγ·dβ,dαδ=dδ·dα,dβδ=dδ·dβ,
where dγ is the first vector of the rectifier section and dδ is the second vector of the rectifier section. Finally, the duty cycles of the zero vectors are calculated as
(14)d0=1-(dγ+dδ)·(dα+dβ),d0γ=dγ(dγ+dδ)d0,d0δ=dδ(dγ+dδ)d0.The overall switching scheme in one period is depicted in Figure7. Fictitious dc-link voltage for simulated IMC is shown in Figure 8.Figure 7
Switching scheme in one period.Figure 8
Fictitious dc-link voltage.


## 4. Vector Control of DFIG

The goal of the DFIG controller is the independent control of the stator active and reactive power. The active power reference is determined by MPPT algorithm and the reactive power is set in order to achieve the desired power factor. Stator fluxd-q reference frame is the most widely used DFIG vector control in the wind turbine applications. Thus, the inverter of IMC is controlled in a stator flux d-q reference frame, with the d-axis oriented along the stator flux vector position. For this reference frame selection, the DFIG model can be written as
(15)ψs=ψsd,0=ψsq.Substituting (15) into (4), we obtain
(16)ids=1Ls(ψs-Lmidr),iqs=-LmLsiqr,ψdr=σLridr+LmLsψs,ψqr=σLriqr,σ=1-Lm2LsLr,
where σ is the leakage factor. In DFIG, the rotor voltages are control variables which control the rotor currents. Substituting (16) into (3), rotor voltages can be written as
(17)vdr=Rridr+ddt(σLridr+LmLsψs)-(ωs-ωr)(σLriqr),vqr=Rriqr+ddt(σLriqr)+(ωs-ωr)(σLridr+LmLsψs).In dynamic performance analysis of the overall system, cross-coupling terms in (17) are added to control loops as feed-forward compensation terms. The stator active and reactive power can be calculated as
(18)Ps=32(vdsids+vqsiqs),Qs=32(vqsids-vdsiqs).In steady state, the stator flux is proportional to the grid voltage. Neglecting the small voltage drop in the stator resistance yields(19)Vs=vqs,0=vds,|Vs|≃ωs|ψs|.Thus, when orienting thed-axis with the stator flux, the voltage aligns with the q-axis. Combining (19) and (16) with (18), we obtain
(20)Ps=32vdsids=-32VsLmLsiqrQs=32vqsids=3Vs2Ls(ψs-Lmidr).The above equations clearly show that, under the stator flux orientation, the active and reactive powers are decoupled and can be controlled via the rotor currents. By means ofiqr, we can control the active power while the reactive power can be controlled via the idr. Using the above equations, the reference currents can be calculated from the desired powers. The schematic diagram of stator flux based vector control of DFIG along with MPPT algorithm is shown in Figure 9.Figure 9
Schematic diagram of vector control of DFIG.


## 5. New LVRT Scheme Based on the Short Circuit Mode

### 5.1. Operation and Implementation

In wind turbines based on back-to-back converters, using the crowbar is the most popular and reliable scheme to relieve the problems of over current in rotor windings and over voltage at dc link. As shown in Figure2, these schemes employ extra circuitry to do the job. However, in patent [32] for back-to-back converter a scheme is presented using only the existing converter switches. In this scheme the crowbar function is performed by handling the fault currents through applying zero vectors by inverter part of the converter. Simultaneous zero vectors are achieved by means of turning on 3 upper switches or lower switches. Here due to the presence of dc-link capacitor turning on both switches of a leg is prohibited to prevent it from short circuit.In the present paper based on this concept a new crowbar scheme is proposed and applied on DFIG excited by an indirect matrix converter (Figure10). In the indirect matrix converter due to the absence of dc-link capacitor, the limitation of simultaneous switching of two switches of each leg is removed. Therefore, in the occurrence of large voltage dips all output stage switches of IMC are turned on and those of input stage are turned off simultaneously. This mode is known as short circuit mode and could do the crowbar task. Therefore, a new low-cost crowbar is realized that does not require any additional hardware. Activation of the short circuit mode is based on the output of a comparator. The comparator block compares the maximum absolute value of the rotor currents with a threshold value. If the output of the comparator is zero, then the gate pulses generated by SVM modulator will apply to the IMC. In case of occurring faults, the rotor currents will exceed the threshold value and comparator output value will change to one. Hence, the controller switches to the short circuit mode and apply the new crowbar. Furthermore, in IMC there are six switches in this mode (Figure 10). So, each switch in this scheme should tolerate half the current of back-to-back converter case. Therefore, it can work well in DFIG based on IMC even for severe faults.Figure 10
New crowbar structure.


### 5.2. Dynamic Performance

The dynamic performance of the proposed scheme is investigated during the paper. Theoretical analysis of the dynamic behavior of the DFIG during three-phase voltage dips has been investigated in previous works, especially in [14, 15]. Both of them proposed inserting extra resistances in stator circuit along with control modification to improve the behavior of the DFIG. During faults the proposed method could be considered as a crowbar having a small resistance. In normal operating condition, the space vector of the stator flux in DFIG rotates at synchronous speed. Referenced to rotor, it rotates at the slip speed and induces small voltages in it. Due to continuity principle of the stator flux, when the DFIG experiences a voltage dip, the stator flux encompasses two components: an ac component corresponds to the remaining voltage of the DFIG terminal and a dc component to maintain flux continuity. The dc component is fixed to the stator, but it is seen by the rotor as a flux that rotates at the rotor speed. Note that the rotor speed is significantly greater than the slip speed and so the induced voltages in rotor are much greater. Therefore, it is desired to increase the damping of the dc term to improve transient performance of the DFIG. Impact of the proposed method on the stator flux dynamic is investigated to compare with disconnecting the converter as follows.If the converter is disconnected from the rotor, the dynamic equation of the stator flux based on (3) and (4) can be calculated by
(21)isdq=ψsdqLs⟹dψsdqdt=-RsLsψsdq-jωsψsdq+vsdq.The pole of the above equation iss=-Rs/Ls-jωs. The damping coefficient of this pole is α=-Rs/Ls and is equal to −0.38 for the studied DFIG. Therefore, the pole is close to imaginary axis with poor damping. The dynamic equation of the stator flux in case of applying the proposed method will change to
(22)isdq=Lrψsdq-LmψrdqLsLr-Lm2⟹dψsdqdt=-RsσLsψsdq-RsσLrψrdq-jωsψsdq+vsdq.The pole of the new equation iss=-Rs/σLs-jωs. The damping coefficient is modified to α=-Rs/σLs, and for the studied DFIG it becomes −25.6 which is quite greater than the previous damping coefficient. Therefore, in the proposed method not only the fault current distributed in more switches but the damping of the stator flux is also greatly increased.


## 5.1. Operation and Implementation

In wind turbines based on back-to-back converters, using the crowbar is the most popular and reliable scheme to relieve the problems of over current in rotor windings and over voltage at dc link. As shown in Figure2, these schemes employ extra circuitry to do the job. However, in patent [32] for back-to-back converter a scheme is presented using only the existing converter switches. In this scheme the crowbar function is performed by handling the fault currents through applying zero vectors by inverter part of the converter. Simultaneous zero vectors are achieved by means of turning on 3 upper switches or lower switches. Here due to the presence of dc-link capacitor turning on both switches of a leg is prohibited to prevent it from short circuit.In the present paper based on this concept a new crowbar scheme is proposed and applied on DFIG excited by an indirect matrix converter (Figure10). In the indirect matrix converter due to the absence of dc-link capacitor, the limitation of simultaneous switching of two switches of each leg is removed. Therefore, in the occurrence of large voltage dips all output stage switches of IMC are turned on and those of input stage are turned off simultaneously. This mode is known as short circuit mode and could do the crowbar task. Therefore, a new low-cost crowbar is realized that does not require any additional hardware. Activation of the short circuit mode is based on the output of a comparator. The comparator block compares the maximum absolute value of the rotor currents with a threshold value. If the output of the comparator is zero, then the gate pulses generated by SVM modulator will apply to the IMC. In case of occurring faults, the rotor currents will exceed the threshold value and comparator output value will change to one. Hence, the controller switches to the short circuit mode and apply the new crowbar. Furthermore, in IMC there are six switches in this mode (Figure 10). So, each switch in this scheme should tolerate half the current of back-to-back converter case. Therefore, it can work well in DFIG based on IMC even for severe faults.Figure 10
New crowbar structure.


## 5.2. Dynamic Performance

The dynamic performance of the proposed scheme is investigated during the paper. Theoretical analysis of the dynamic behavior of the DFIG during three-phase voltage dips has been investigated in previous works, especially in [14, 15]. Both of them proposed inserting extra resistances in stator circuit along with control modification to improve the behavior of the DFIG. During faults the proposed method could be considered as a crowbar having a small resistance. In normal operating condition, the space vector of the stator flux in DFIG rotates at synchronous speed. Referenced to rotor, it rotates at the slip speed and induces small voltages in it. Due to continuity principle of the stator flux, when the DFIG experiences a voltage dip, the stator flux encompasses two components: an ac component corresponds to the remaining voltage of the DFIG terminal and a dc component to maintain flux continuity. The dc component is fixed to the stator, but it is seen by the rotor as a flux that rotates at the rotor speed. Note that the rotor speed is significantly greater than the slip speed and so the induced voltages in rotor are much greater. Therefore, it is desired to increase the damping of the dc term to improve transient performance of the DFIG. Impact of the proposed method on the stator flux dynamic is investigated to compare with disconnecting the converter as follows.If the converter is disconnected from the rotor, the dynamic equation of the stator flux based on (3) and (4) can be calculated by
(21)isdq=ψsdqLs⟹dψsdqdt=-RsLsψsdq-jωsψsdq+vsdq.The pole of the above equation iss=-Rs/Ls-jωs. The damping coefficient of this pole is α=-Rs/Ls and is equal to −0.38 for the studied DFIG. Therefore, the pole is close to imaginary axis with poor damping. The dynamic equation of the stator flux in case of applying the proposed method will change to
(22)isdq=Lrψsdq-LmψrdqLsLr-Lm2⟹dψsdqdt=-RsσLsψsdq-RsσLrψrdq-jωsψsdq+vsdq.The pole of the new equation iss=-Rs/σLs-jωs. The damping coefficient is modified to α=-Rs/σLs, and for the studied DFIG it becomes −25.6 which is quite greater than the previous damping coefficient. Therefore, in the proposed method not only the fault current distributed in more switches but the damping of the stator flux is also greatly increased.


## 6. Simulation Results

In this section the performance of a DFIG wind turbine is simulated under a severe fault condition. Simulations are carried out in PSIM environment to investigate the effects of switching. As was mentioned in Section2 according to LVRT standard of Germany, wind turbines must remain connected to grid following occurrence of grid faults. Moreover, they must also inject reactive power to help the grid for voltage recovery. To simulate the fault condition, at t=0.7 s the grid voltages are dropped to 0.2 pu for 300 ms as shown in Figure 11. Figure 12 shows the rotor currents in this case without any protection. For such a severe voltage dip, the rotor currents are increased to 2.8 pu. In the absence of any type of protecting method, these high currents can damage the IGBT switches of the IMC. Figure 13 shows the rotor currents when the proposed method is in operation. Short circuit activation signal is obtained from comparison of rotor currents with a threshold value. In this simulation the threshold value is set to 1.5 pu. When the rotor current is greater than the threshold value, the short circuit mode is activated. It means that all of the output stage switches are turned on. As was mentioned in the previous section, here the high fault currents of the rotor are distributed between output switches.Figure 11
Grid voltages.Figure 12
Rotor currents in fault condition, with the controller designed for normal operation.Rotor currents in fault condition, with the proposed method.
(a)
(b)
(c)Figure14 shows the currents flowing in IGBT switches in time interval from t=0.72 s up to t=0.77 s during activation of the short circuit mode. As shown in Figure 14, the currents through all switches are in the safe operating area (SOA). Furthermore, the duration of short circuit mode is less than 20 ms, so there is a chance for reactive power injection as required by new LVRT standards. Short circuit mode activation signal is presented in Figure 13. It implies that the short circuit mode will be instantly activated if the rotor currents exceed the threshold value. The short circuit mode is active only for a short time, so, for most of the time during fault, the control system of DFIG is in operation to inject reactive power to grid. In order to meet the new LVRT standards, during fault the references of the active and reactive powers are set to 0.0 pu and 1.0 pu, respectively. As can be seen from Figure 13, injection of 1.0 pu reactive current is realized in compliance with LVRT standard.Rotor currents in fault condition, with the proposed method.
(a)
(b)


## 7. Conclusion

To enhance the LVRT capability of DFIG wind turbines based on indirect matrix converters (IMCs), in this paper a new crowbar system is proposed. This crowbar system is based on using the existing converter switches in short circuit mode to protect the IMC against large fault currents. Therefore, this method is interesting as it does not require any additional hardware. Even in severe fault conditions, the short circuit mode is active only for short time, so the control system can be in operation at most of the fault duration. Therefore, the reactive power injection is realized to help voltage recovery in order to satisfy the new LVRT standards. Simulation results confirm the efficiency of the proposed method.



---

*Source: 157431-2013-12-12.xml*