State of The Art of Wind Turbines Modelling With Induction Generator (DFIG)
State of The Art of Wind Turbines Modelling With Induction Generator (DFIG)
State of The Art of Wind Turbines Modelling With Induction Generator (DFIG)
Juan Alvaro
Fuentes and Angel
Molina are with the Department of Electrical Engineering, Universidad Politecnica de Cartagena, Cartagena 30202
Spain; (emails: [email protected], [email protected]). (phones:
+34 968 325 604, +34 968 325 462, fax: +34 968 325 356)
Fig. 1.
SCIG Squirrel Cage IG; WRIG Wound Rotor IG; DFIG Double
Fed IG; PMSG Permanent Magnet SG; WRSG Wound Rotor SG.
TABLE I
Type
Type
Type
Type
A
B
C
D
Stall
Type A0
Type B0
Type C0
Type D0
Power control
Pitch
Active stall
Type A1
Type A2
Type B1
Type B2
Type C1
Type C2
Type D1
Type D2
(1)
where = ej2/3 , 1 is a phase shift and is the synchronous angle corresponding to the synchronous frequency,
Fig. 2.
(2)
(3)
Fig. 3.
v~s = Rs i~s jd ~s
v~r = Rr i~r jdm ~r + d~r /dt
(6)
= (Ls + Lm )/Lm , LM = Lm
LL = Ls 2 Lr , RR = 2 Rr
vR = vr , iR = ir /, R = r
(7)
Fig. 4.
harmonic frequency components. In the steady state, the analysis is simplified by resolving each variable into a component
at fundamental frequency s , plus the main harmonics [4].
Note that, when the rotor voltage (Vr or VR ) is short
circuited, the equivalent circuit for the DFIG becomes the
ordinary equivalent circuit for a SFIG.
B. Mechanical modelling
The drive-train of a wind turbine generator system (WTGS)
in general consists of a blade pitching mechanism with a
spinner, a hub with blades, a rotor shaft relatively long, and a
heavy gearbox with breaker and generator. The inertia moment
of the wind wheel (hub with blades) is about 90% of the drivetrain total moment, while the generator rotor inertia moment
is about 6-8%. The remaining parts of the drive train comprise
the rest (2-4%) of the total inertia moment. At the same
time, the generator represents the biggest torsional stiffness.
The rotor shaft stiffness is about 100 times less and the
stiffness of the hub with blades is about 50 times less than the
generator stiffness [6]. Because the torsional vibration of the
drive-train elements are inevitable in the low-speed side, and
moreover their character in the high-speed side (e.g. frequency,
amplitude) can highly influence the WTGS performance.
Regarding to the mechanical part of the WT, the general
way for modelling the drive-train treats the rotor as a number
g
g
0
00 dw
Jwt
dt = Twt D1,2 (w NGB ) K2M (w NGB )
D
d
1,2
g
g
g
Jg0 dt = Te NGB (w NGB ) K2M (w NGB )
(8)
00
being Jwt
, Jg0 the equivalent inertia constants, K2M the
equivalent shaft spring constant; D1,2 the equivalent damping
0
, Te the mechanical and electrical torques; w ,
coefficient; Twt
g the mechanical angular speeds (elect. rad/s); w , g the
mechanical angular positions (elect. rad); NGB = gn /wn
the gearbox transformation ratio, and wn , gn the windwheel and generator rated speeds, respectively.
The WTGS drive-train is not advisable to be modelled as
a single lumped mass in transient analysis because it is too
simple for representing the dynamics of a WT, and stability
analysis based on this model may give significant errors,
especially in FSWT [7]. However, in the steady-state operation
analysis the minimal realization of the drive-train model is
based on the one-lumped mass [8], whose model can be
described by the following simple equation:
Fig. 5.
(10)
CP max
3opt
(11)
000
00
Jwt
dr /dt = Twt
Te D1M r
(9)
Popt =
000
Jwt
where
is the wind turbine equivalent inertia constant;
00
D1M is a friction loss coefficient (damping), Twt
is the equivalent mechanical torque; and r is the mechanical angular
speed (elect. rad/s).
C. Wind Turbine Modelling
The mechanical power generated for the WT is given by:
p
Popt
0
3
= Kopt
Popt
Kopt
s
p
Topt
00
=
= Kopt
Topt
Kopt
(12)
(13)
(15)
II) 1 < s < Sr , sPgs > 0, Pgr > 0, sPgs = Pgs + Pcur
III. OVERVIEW OF POWER CONTROLLABILITY
Power flow in DFIG: a) sub-synchronous, b) supersynchronous mode (Pcur copper losses, Sr rated slip)
Fig. 6.
(14)
Fig. 8.
around the rated mean value of the power. The third option
is the active stall control, i.e. the stall of the blade is actively
controlled by pitching the blades. At low wind speeds the
blades are pitched similar to a pitch-controlled WT, in order
to achieve maximum efficiency. At high wind speeds the
blades go into a deeper stall by being pitched slightly into the
direction opposite to that of a pitch-controlled WT. The active
stall WT achieves a smoother limited power, without high
fluctuations. This control type has the advantage of being able
to compensate variations in air density. The combination with
the pitch mechanism makes it easier to carry out emergency
stops and to start up the WT.
The main tasks of the wind wheel (turbine) control system
are to keep the value of the rotor shaft speed and the active
power at a set levels. Therefore, the control system usually
consists of speed, power and blade pitch angle regulators, as
can be seen in the example in figure 8.
B. Speed control subsystem
When the wind speed is under its rated value, the DFIG
is working at partial load, being this state the most regular
operating mode of the WT, because the wind speed is usally
lower than the rated value. In high power WT systems is
very important to implement a specific control system to
improve the energy efficiency of the wind power. There is
three methods, at least, to carry out the task of the control
system. Regarding the optimal power curve, the rotor speed is
related to the maximum power by a optimal constant, as can
be seen in subsection II-C.
QG
/t
P/g
1/s
Logic
PG
Calcul.
IqG
VabcG
Fault
TRF - L
3 --> 2
VdG
R wire
pcc
RL filter
PR, QR
ref
Defuzzification
RSC
RL filter
DFIG
PI
Tg
IdR
S_C
Torque
3 --> 2
IqR
Calqul.
Generator
Reference
IabcR
Crowbar
DC link
GSC
PWM
PWM
VabcR
Equivalent
control
input
ref
eq
+
+
Alternate
Control
input
P*
-
PS, QS
g*
b)
1/s
Control
rules
evaluation
Fuzzification
/t
Tg
PI
PG, QG
IabcG
3 --> 2
VqG
a)
Grid
TRF - H
IdG
PQ
RL wire
Gradient
detection
1/s
T_C
Q_C
VdR
Id_C
+
Iq_C
VqR
2 --> 3
Vdc
GCC
IRn
CVC
Qref
Tg
Fig. 11.
Fig. 12.
c)
Fig. 9.
Fig. 10.
Fig. 13.
when the wind speed changes suddenly either when the grid
connection of the WT is affected by some type of electrical
fault too. Regarding the order of the mechanical model, it is
proved that as wind speed is intermittent and stochastic in
nature, the torques acting on blades of wind turbines are not
always equal. This effect can only be analysed through the
6-mass drive train model. But this mechanical model slows
the simulations due to complex and lengthy mathematical
computation with small time-steps. Moreover, the unequal
blade torque distributions do not have effect on the transient
stability of WTGS, which encourages the consideration of the
reduced order 3-mass model [7], although 2-mass model is
the most used in the references. Finally, control systems of
the WT and DFIG are included, but they usually vary, since
their structures and parameters depend on the WTGS specific
model, the rated power and the rotor speed range variation [6].
R EFERENCES
[1] T. Ackermann, Wind Power in Power Systems. John Wiley & Sons,
Ltd,, 2005.
[2] A. Hansen, Generators and power electronics for wind turbines, John
Wiley & Sons, Ltd, p. 24, 2004.
[3] A. Feijoo, J. Cidras, and C. Carrillo, A third order model for the doublyfed induction machine, Electric Power Systems Research, vol. 56, pp.
121127, 2000.
[4] G. Slemon, Modelling of induction machines for electric drives, IEEEJ-IA, vol. 25, no. 6, pp. 11261131, 1989.
[5] A. Petersson, L. Harnefors, and T. Thiringer, Evaluation of current
control methods for wind turbines using doubly-fed induction machines,
IEEE Transactions on Power Electronics, vol. 20, pp. 227 235, 2005.
[6] Z. Lubosny, Wind Turbine Operation in Electric Power Systems. Advanced modeling. Berlin, Germany, Springer-Verlag, 2003.
[7] S. Muyeen, M. Hasan Ali, and R. Takahashi, Comparative study
on transient stability analysis of wind turbine generator system using
different drive train models, Renewable Power Generation, IET, vol. 1,
pp. 131141, 2007.
[8] J. Ekanayake, L. Holdsworth, and W. XueGuang, Dynamic modeling
of doubly fed induction generator wind turbines, IEEE Transactions on
Power systems, vol. 18, pp. 803809, 2003.
[9] R. Datta and V. Ranganathan, Variable-speed wind power generation
using doubly fed wound rotor induction machine-a comparison with
alternative schemes, IEEE Transaction on Energy Conversion, vol. 17,
no. 3, pp. 414421, Sept. 2002.
[10] A. Pujante, E. Gomez, and J. Fuentes, A dfig dynamic model for its
use on wind power generation studies, International Conference on
Sustainable Energy Technologies (ICSET), pp. 1028 1033, 2008.
[11] I. Cadirci and M. Errniq, Double-output induction generator operating at subsynchronous and supersynchronous speeds: steady-state
performance optimisation and wind-energy recovery, IEE ProceedingsElectric Power Applications, vol. 139, pp. 429 442, 1992.
[12] I. Munteanu, A. I. Bratcu, N.-A. Cutululis, and E. Ceanga, Optimal
Control of Wind Energy Systems: Towards a Global Approach, 1st ed.
Springer, March 2008, iSBN-13: 978-1848000797.
[13] E. Hau, Wind turbines. fundamentals, technologies, application, economics, Springer, 2006.
[14] D. Novotny and T. Lipo, Vector control and dynamics of ac drives,
Clarendon Press, Oxford,, 1996.
[15] N. Mohan, T. Undeland, and W. Robbins, Power electronics: Converters,
Applications, and Design. 2nd Edition, John Wiley & Sons, Inc., 1995.
[16] M. Rodriguez, G. Abad, I. Sarasola, and A. Gilabert, Crowbar control
algorithms for doubly fed induction generator during voltage dips, EPE
2005 Dresden, pp. 110, 2005.
[17] J. Morren, W. Sjoerd, and H. de Haan, Short-circuit current of wind
turbines with doubly fed induction generator, IEEE transactions on
Energy conversion, 2007.