Book Chapter SenroySir
Book Chapter SenroySir
Book Chapter SenroySir
8.1 Introduction
Wind turbine generators are increasingly being employed world-wide to
harness wind energy for the purpose of generating electricity. Worldwide
initiatives in smart grids have led to increasing interest in wind energy
conversion systems, not only because the primary resource is wind
which is free, but also because the environmental impact of conventional
generation systems is unsustainable. But the introduction of wind power
into the existing grids at bulk power levels has opened up new challenges
in system operation and control. Therefore, it is important to understand
the dynamics of wind turbine generators so as to enable high penetration
levels of wind power to be achieved in existing grids.
Most popular designs of wind turbine generators employ induction
generators or permanent magnet synchronous generators (PMSG).
The PMSG requires a power electronic converter interface, through
which the generated power flows into the network. The rating of the
converter, must match that of the machine, making these designs more
expensive. The induction generators are relatively cheaper owing to the
absence of power electronic converters, and more popular. For greater
controllability, at higher ratings the induction generators used are of
doubly fed type i.e. both stator and rotor are connected to a voltage
source. While the stator is directly connected to a three-phase network
(typically at 500-700 V), the rotor is also connected to the network
(usually through a transformer) via a back-to-back power electronic
converter with a capacitor connected in between as shown in Figure 8.1.
Such generators are called doubly fed induction generators (DFIG). This
chapter is concerned with the dynamics of DFIG systems.
274 Wind Turbine Generators in Power Systems
rAv 3 C pMAX
where kopt = BASE
indicates the maximum power (p.u.) that
2S BASE
can be extracted by the turbine at the base wind speed. For optimum
l v
power extraction, the turbine is rotated at a speed, wturbine = n om wind .
R
The normalised tip-speed ratio is defined as the ratio of the per-unit
turbine speed and the per-unit wind speed.
l wturb
lp.u. = = (8.5)
ln om vw
Equation (6) assumes that the base rotor speed is the speed that
enables optimum power extraction at base wind speed. A generator
278 Wind Turbine Generators in Power Systems
8.3 Modelling of DFIG
The equivalent circuit of a DFIG is given in Figure 8.5. The rest of the
discussion concerns with two pole machines where the electrical speed
equals the mechanical speed. In general, the electrical speed and the
mechanical speed are related through the number of pole pairs, np, as
w
wmech = elec . Figure 5 shows the rotor circuit referred to the stator
np
which is stationary while in reality the rotor circuit rotates at slip speed,
wr = s ws, where ws is the stator synchronous frequency.
280 Wind Turbine Generators in Power Systems
8.4 Equivalent Circuit
The electrical equivalent circuit of a DFIG consists of the stator
resistance, rs, in series with the stator and rotor leakage reactance,
xr and xr. The shunt branch reactance, xm, is significant in magnitude
and represents the magnetising inductance caused due to the presence of
the low-reluctance air gap within the machine. The rotor reactance, rr,
is also included in the rotor circuit in series with a controlled voltage
V
source, r , indicating the rotor injected voltage.
s
Pm 1 (1 – s ) (1 – s )
Tm = = Vr I r cos q + I r2rr (8.8)
wr wr s s
The torque-slip characteristics of the DFIG is given at different rotor
injected voltage in Figure 8.6. With respect to the stator voltage, the
injected rotor voltage, Vr, may be resolved into the in-phase component,
Vqr, and the out-of-phase component, Vdr . Independent control of these
two components can significantly affect the speed of the DFIG, as
evident from Figure 8.6. In general, the operation of the DFIG may be
summarized as follows: At sub-synchronous speeds, (s > 0), the DFIG
absorbs power through the rotor circuit, but generates power through
the stator. At super-synchronous speeds, the rotor circuit generates
power along with the stator.
d wt 1
= (T – Tshaft ) (8.9)
dt 2H t turb
dqtw
= w B (wturb – wr ) (8.10)
dt
dwr 1
= (T – Te ) (8.11)
dt 2H g shaft
¢
deqs Ê ¢
e qs ˆ
= ws w B Á K 2rr ids – + s e ds
¢ – Kvdr ˜ (8.12)
dt Ë ws t r ¯
¢
deds Ê e¢ ˆ
= ws w B Á – K 2 rriqs – s eqs¢ – ds + Kvqr ˜ (8.13)
dt Ë ws t r ¯
dids wB 2
e ′qs w
= − ws LSC iqs − (rs + K rr ) ids + + r eds
′ + Kvdr − vds
dt LSC ws t r ws
(8.15)
1 Lr + Lm w
where K = , tr = , LSC = Ls + (Lr ||Lm ), and s = 1 – r .
Lr rr ws
1+
Lm
Open Loop Operation 283
Q
Vs2 – VsV∞ cos g e
stator
Qgrid = = vds iqs – vqs ids (8.17)
Xe
e ′ds
iqr = – – Kiqs (8.19)
Xm
The complete model of the DFIG connected to an infinite bus consists
of differential equations (9)-(15) and algebraic equations (16)-(17) and
is termed as system of differential algebraic equations (DAE). The
parameters used for simulation are given in Table 9.1. The equations
may be written in the form
X = F(X, Z, U)
0 = G(X, Z, U) (8.20)
where X = [wturb, qtw, wr, eqs¢ , e¢ds, iqs, ids]T is the vector of system states,
Z = [Vs, ge]T is the vector of algebraic variables, and U = [vw, vqr, vdr]T
is the vector of inputs (wind speed and rotor injected voltage). If the
DFIG is to be connected to the infinite bus directly i.e. Xe = 0, its
terminal voltage magnitude and phase no longer depend on the power
injected into the infinite bus. In that case, (16)-(17) become irrelevant
and the system of DAE of (20) reduces to a set of seven differential
equations only.
284 Wind Turbine Generators in Power Systems
wB ksh Ht Hg Lm Ls Lr Rs Rr Xe
314.16 0.3 4 0.4 5.419 0.167 0.132 0.0084 0.0083 0.06
Figure 8.7: Open loop operation of DFIG. Wind speed is 12 m/s from
0-10 s, 13 m/s from 10-20 s and 11 m/s from 20-30 s.
Linearization of DFIG – Open Loop 285
The variation of C ppu with changes in wind speed for the same
rotor voltage injections as in Figure 8.7 Open loop operation of DFIG.
is shown in Figure 8.8. It may be observed that the SCIG develops
optimum power (Cppu = 1) at 12 m/s, but the reduced value at other
wind speeds indicate that the MPPT operation is not achieved. On the
other hand, vdr = – 0.01, vqr = 0.15, helps the DFIG achieve MPPT
at 11 m/s, while vdr = – 0.01, vqr = – 0.03 helps attain MPPT at
13 m/s.
vdr = 0 49.96 Hz, 0.028 5.70 Hz, 0.114 0.52 Hz, 0.148
vqr = 0
vdr = – 0.01 49.96 Hz, 0.028 6.71 Hz, 0.115 0.55 Hz, 0.213
vqr = – 0.03
vdr = – 0.01 49.96 Hz, 0.028 9.30 Hz, 0.061 0.59 Hz, 0.365
vqr = 0.15
du4 – 1 e′ (8.28)
= dREF – tan ds
dt eqs′
du4
d r = ki 4u 4 + k p 4 (8.29)
dt
Closed Loop Operation 289
′ iqs + eds
eqs ′ ids
The electric torque, Te, is calculated as Te = . The rotor
ws
voltage thus obtained, may be resolved into the d-q components as
vqr = Vr cos dr; vdr = Vr sin dr (8.30)
The dc-link voltage is governed by the energy balance between the
rotor side converter and the grid side converter.
dVdc Protor Pc
CVdc = – (8.31)
dt vdr idr + vqr iqr vd g idg + vqg iqg
The rotor side and grid side converter currents are calculated as
eqs′ eds
′
idr = – Kids ; iqr = – – Kiqs (8.32)
ws Lm ws Lm
(vqg – Vs ) vdg
idg = – ; iqg = (8.33)
xTFR xTFR
1
where K = . The GSC controller regulates the dc-link voltage
Lr
1+
Lm
and the reactive power exchanged through the GSC.
du5
= VdcREF – Vdc (8.34)
dt
vdg = – xTFR k u + k du5 (8.35)
i5 5 p5
dt
Qc
du6
= QcREF – vdg iqg – vqg idg (8.36)
dt
du
vqg = – xTFR k16u6 + k p 6 6 – (VsREF – VS ) (8.37)
dt
Z = [Vs, ge, vdg, vqg]T, corresponding to the equations (35) and (37)
and the following two equations modified from equations (16) and (17)
to include the effect of the converter and the dc link dynamics.
P
VV
stator
s ∞
Pgrid = sin g e = vds idg + vqs iqg + vds ids + vqs iqs (8.38)
Ve
Q
Vs2 – VsV∞ cos g e
stator
Qgrid = = vds iqs – vqs ids + vds iqg – vqs idg (8.39)
Xe
In order to solve for the behaviour of the DFIG, all the differential
and algebraic equations have to be solved simultaneously using an
appropriate solver. The plot of Cp in Figure 8.13 shows that the action
of the controllers restores the DFIG operation to its MPPT level after
every change in wind speed occurs.
The tuning of the controllers is key to the performance of the
FMAC strategy. There are six proportional-integral (PI) controllers,
and they affect the DFIG performance in different ways. The
sub-optimal values selected for the simulations in this chapter are
presented in Table 8.4.
vw = 11 m/s vw = 13 m/s
Pgrid 0.49 0.85
Qgrid 0.01 0.02
wr 0.92 1.08
E¢ 1.55–5.8° 1.57–8.4°
Vr 0.14–– 1.7° 0.13–– 163.6°
Is 1.91–– 72.2° 2.02–– 65.9°
Ir 2.15–106.2° 2.25–112.1°
ys 1.01–89.1° 1.00–89.1°
yr 1.59–95.8° 1.61–98.4°
292 Wind Turbine Generators in Power Systems
The key modes of the system are shown in Table 8.6 at 11 m/s.
The electromechanical mode for this system has a frequency of 1.81 Hz,
and this does not change appreciably with the wind speed.
dqtw
F2 = = w B (wt – wr ) (8.42)
dt
Te
K f fD
dwr 1 2
F3 = = k q – k w – (8.43)
dt 2H g sh tw opt r wr
k C v3
∂F1 ∂ 1 opt ppu w
= – ksh qtw
∂wturb ∂wturb 2H t wturb
k C v3 kopt vw3 ∂C p
1 opt pREF w 0
= – + 0 pu
2H t w 2 w ∂ w turb
turb0 turb0
= – + 0 pu
2H t w 2 wturb ∂lpu ∂wturb
turb0 0
The inputs to the linearized model are the deviations in the wind
speed and the network frequency from their initial operating conditions.
Accordingly, the sensitivity of F1, F2 and F3 to these inputs is
296 Wind Turbine Generators in Power Systems
∂F1 ∂C p
1
= 3koptC pREF vw2 + kopt vw3 pu
∂vw 2H t wturb 0 0 ∂vw
0
∂C p ∂l
1 pu
= 3koptC pREF vw2 + kopt vw3 pu
2H t wturb 0 0 ∂lpu ∂vw
0
1
= (3koptC pREF vw2 – kopt vw K(C , l ) wturb ) (8.48)
2H t wturb 0 0 p 0
0
3
kopt vw3 ∂C p kopt vw3 ∂C p ∂b
∂F1 1 ∂ koptC ppu vw
= = 0 pu
= 0 pu
kopt vw3
=
0
K(C Kb
, b) (8.49)
2H t wturb p
0
The output of the system, being the deviation in the electrical power
output of the wind turbine generator, may be considered equal to the
change in PeREF if the generator dynamics are ignored.
PeD = 3kopt wr 20 wrD – Kf fD (8.52)
The initial values of the states and algebraic variables in (51) may
be obtained by setting the left hand side of (41)-(43) to zero and
solving simultaneously. At 10 m/s, the values obtained for 10% deloading
are v0 = 0.833, wturb
0 = w r 0 = 0.8054, kopt = 0.73, b = 1°, Kf = 15,
Primary Frequency Response from Dfig-Wind Turbines 297
P1 = kdel w3 (8.53)
3
lnom
where kdel = kopt . The value of ldel should be such that the
3
ldel
corresponding Cp(l, b) obtained from equation (1) is deloaded as per
the requirements. For instance, setting ldel = 9.5909 instead of its
nominal value of 8.1 yields Cp(l, b) = 0.432, kdel = 0.44. Consequently,
the power is deloaded by 10% at every wind speed, when deloaded to
the right. Figure 8.17 shows the deloaded power reference curve along
with the MPPT curve, superimposed upon the turbine characteristics
at a wind speed of 10m/s. The operating point shifts to the right due
to deloading i.e. instead of operating at w0 to generate P0, it operates
at a higher w1 to generate P1, which is 90% of P0. During primary
frequency response, the operating point is shifted along the straight line
to {w, Pdel} as indicated; the new power reference is calculated as
P0 – P1
Pdel = P1 + (w1 – w ) (8.54)
w1 – w0
At high wind speeds, and non-zero pitch (owing to pitch based
deloading), the speed increase associated with deloading can become
a limiting factor in deciding the extent of deloading because the rotor
side converter rating depends on the maximum slip.
8.13 Remarks
Increasing the penetration of wind generation into bulk power grids,
can deteriorate the system frequency response in many ways. Firstly,
replacing conventional synchronous generators (employed in thermal,
hydro and nuclear power plants) with induction generators (associated
with wind generation systems) can lead to significant reduction in system
inertia. This is because the speed of induction generators interfaced to
the grid via power electronic converters, is de-linked from the system
frequency. Hence, any small disturbance can lead to increased frequency
excursions. Further, speed governors and automatic voltage regulators
employed with conventional generators provide primary and secondary
frequency response as well as power oscillation damping, and aid the
recovery of the system during disturbances. These controller technologies
have not matured enough in the context of induction generators, so their
impact is reduced. Moreover, an inherent uncertainty is associated with
wind as a resource, and the spinning reserve available at any point of
time due to deloading of wind generators is unreliable. Thus, the overall
system stiffness decreases with increase in wind generation.
The annual capacity factor (CF) of any generator is defined as the
total annual energy produced divided by the product of the rating of
the generator and the number of hours in a year (= 8760). For a wind
farm, the CF is typically in the range of 20% – 30%. Considering that
all the wind turbines in a wind farm are deloaded by 10%, the CF
may reduce to 27% from 30%. Due to various frequency events, and
opportunities for the wind turbine to contribute its primary frequency
reserve, this CF may ultimately improve from 27% to somewhere
close to 30%. Therefore, the loss of revenue due to deloading, must
be compensated by appropriate ancillary service pricing mechanism to
actively encourage wind farm operators to opt for deloading strategies.
The selection of the frequency droop factor, Kf, also has an important
impact on the CF of the wind turbine. Generally, at locations with
high prevailing wind speeds, the spinning reserve is likely to higher and
therefore Kf should be chosen to enable maximum frequency response.
The CF accordingly would improve substantially.
8.14 Multi-Machine Simulation
The need for simulation of multiple wind turbine generators may
arise when wind farms connected to a grid have to be studied. If the
dynamic behaviour occurring at the high-voltage grid is of interest, it
is frequently convenient to represent the entire farm as a dynamically
equivalent machine with the aggregated injection. This is because the
300 Wind Turbine Generators in Power Systems
8.15 Conclusion
In this chapter detailed modelling of DFIG based wind turbine generators
is discussed and its linearization presented. Small signal analysis of the
DFIG was carried out at different wind speeds. The ability to respond
to network frequency disturbances was also analysed using mathematical
models.
8.16 References
1. O. Anaya-Lara, N. Jenkins, J. Ekanayake, P. Cartwright, and M.
Hughes, Wind Energy Generation – Modelling and Control, John
Wiley and Sons Ltd., United Kingdom, 2009
2. B. Wu, Y. Lang, N. Zargari, and Z. Kuo, Power Conversion and
Control of Wind Energy Systems, Wiley-IEEE Press, New Jersey,
2011
3. Y. Mishra, S. Mishra, M. Tripathy, N. Senroy, Z. Y. Dong, “Improving
stability of a DFIG-based wind power system with tuned damping
controller,” IEEE Trans. Energy Conv., vol. 24, No. 3, pp. 650-660,
2009
4. S. Ghosh and N. Senroy, “Electromechanical dynamics of controlled
variable speed wind turbines,” IEEE Sys. J., vol. 9, No., pp. 639-
646, 2015
5. F. M. Hughes, O. Anaya-Lara, N. Jenkins, and G. Strbac, “A power
system stabilizer for DFIG-based wind generation,” IEEE Trans.
Power Sys., vol. 21, No. 2, pp. 763-772, May 2006
6. R. Datta and V. T. Ranganathan, “A method of tracking the peak
power points for a variable speed wind energy conversion system,”
IEEE Trans. Energy Conv., vol. 18, No. 1, pp. 163-168, 2003
References 301