DFIG-Based Wind Turbine Robust Control Using High-Order Sliding Modes and A High Gain Observer
DFIG-Based Wind Turbine Robust Control Using High-Order Sliding Modes and A High Gain Observer
DFIG-Based Wind Turbine Robust Control Using High-Order Sliding Modes and A High Gain Observer
mg
= Generator speed;
T
g
= Generator electromagnetic torque;
J = Turbine total inertia;
K = Turbine total external damping;
s (r) = Stator (rotor) index;
d (q) = Synchronous reference frame index;
V (I) = Voltage (current);
P (Q) = Active (reactive) power;
| = Flux;
T
em
= Electromagnetic torque;
R = Resistance;
L (M) = Inductance (mutual inductance);
o = Leakage coefficient, o = 1 M
2
/L
s
L
r
;
r
= Rotor position;
r
(
s
) = Angular speed (synchronous speed);
s = Slip;
p = Pole pair number.
I. Introduction
Variable speed wind turbines are continuously
increasing their market share, since it is possible to track
the changes in wind speed by adapting shaft speed and
thus maintaining optimal power generation. The more
variable speed wind turbines are investigated, the more it
becomes obvious that their behavior is significantly
affected by the used control strategy. Typically, they use
aerodynamic controls in combination with power
electronics to regulate torque, speed, and power. The
aerodynamic control systems, usually variable-pitch
blades or trailing-edge devices, are expensive and
complex, especially when the turbines are larger. This
situation provides a motivation to consider alternative
control approaches and to introduce more intelligence [1].
Therefore, this paper deals with power generation
control in variable speed wind turbines. A control strategy
is proposed to ensure power extraction optimization and
reduce mechanical stresses in the drive train of a DFIG-
based wind turbine. The proposed control strategy
combines an MPPT using a high gain observer and
second-order sliding mode for the DFIG control. This
strategy presents attractive features such as chattering-free
behavior, finite reaching time, robustness and unmodeled
dynamics (generator and turbine).
The overall control strategy, illustrated by Fig. 1, has
been validated on a 1.5-MW three-blade wind turbine
using the NREL wind turbine simulator FAST.
II. Wind Turbine Modeling
The global scheme for a grid-connected wind turbine is
given in Fig. 2. The wind turbine modeling is inspired
from [2]. In the following, the wind turbine components
models are briefly described.
II.1. Turbine Model
In this case, the aerodynamic power P
a
captured by the
wind turbine is given by
( )
2 3
1
2
a p
P R C v = t (1)
where
R
v
e
= (2)
h
a
l
-
0
0
6
3
0
4
3
6
,
v
e
r
s
i
o
n
2
-
4
N
o
v
2
0
1
1
Author manuscript, published in "International Review on Modelling and Simulations 4, 3 (2011) 1148-1155"
DFIG
G
e
a
r
Grid 50 Hz
Stator power
Rotor
power
Rotor side
converter
Grid side
converter
Second-Order
Sliding Mode Controller
P
o
w
e
r
Speed
v
1
v
2
v
n MPPT
u v
I
r
e
r
T
em
DFIG
G
e
a
r
Grid 50 Hz
Stator power
Rotor
power
Rotor side
converter
Grid side
converter
Second-Order
Sliding Mode Controller
Second-Order
Sliding Mode Controller
P
o
w
e
r
Speed
v
1
v
2
v
n
P
o
w
e
r
Speed
v
1
v
2
v
n MPPT
u v
I
r
e
r
T
em
Fig. 1. The proposed control structure.
WT Aerodynamics
Gearbox DFIG & Converters
Grid connection
Wind
WT Aerodynamics
Gearbox DFIG & Converters
Grid connection
Wind
Fig. 2. Wind turbine global scheme.
Normally, a variable speed wind turbine follows the
C
pmax
to capture the maximum power up to the rated
speed by varying the rotor speed to keep the system at
opt
. Then it operates at the rated power with power
regulation during high wind periods by active control of
the blade pitch angle or passive regulation based on
aerodynamic stall.
The rotor power (aerodynamic power) is also defined
by
a a
P T = e (3)
According to [2], the following simplified model is
adopted for the turbine (drive train) for control purposes.
a g
J T K T e = e (4)
II.2. Generator Model
The WT adopted generator is the DFIG. DFIG-based
WT will offer several advantages including variable
speed operation (33% around the synchronous speed),
and four-quadrant active and reactive power capabilities.
Such system also results in lower converter costs
(typically 25% of total system power) and lower power
losses compared to a system based on a fully fed
synchronous generator with full-rated converter.
Moreover, the generator is robust and requires little
maintenance [3-4].
The control system is usually defined in the
synchronous d-q frame fixed to either the stator voltage
or the stator flux. For the proposed control strategy, the
generator dynamic model written in a synchronously
rotating frame d-q is given by (5).
( )
sd
sd s sd s sq
sq
sq s sq s sd
rd
rd r rd r rq
rq
rq r rq r rd
sd s sd rd
sq s sq rq
rd r rd sd
rq r rq sq
em rd sq rq sd
d
V R I
dt
d
V R I
dt
d
V R I
dt
d
V R I
dt
L I MI
L I MI
L I MI
L I MI
T pM I I I I
|
= + e |
= + + e |
= + e |
|
= + + e |
| = +
| = +
| = +
| = +
(5)
h
a
l
-
0
0
6
3
0
4
3
6
,
v
e
r
s
i
o
n
2
-
4
N
o
v
2
0
1
1
For simplification purposes, the q-axis is aligned with
the stator voltage and the stator resistance is neglected
[4]. These will lead to (6).
1
1
rd sd
rd r rd s r rq
r s
rq
rq r rq s r rd s sd
r s
em sd rq
s
dI d M
v R I s L I
dt L L dt
dI
M
v R I s L I s
dt L L
M
T p I
L
| | |
= + e o
|
o
\ .
| |
= e o e |
|
o
\ .
= |
(6)
III. The MPPT Strategy
III.1. Problem Formulation
The reference torque of the MPPT block must address
two problems: the captured power maximization and be
driven according to the three fundamental operating
regions illustrated by Fig. 3.
In practice, there are two possible regions of turbine
operation, namely the high- and low-speed regions. High
speed operation (III) is frequently bounded by the speed
limit of the machine. Conversely, regulation in the low-
speed region (II) is usually not restricted by speed
constraints. However, the system has nonlinear
nonminimum phase dynamics in this region. This adverse
behavior is an obstacle to perform the regulation task [5].
In region II, the control objective is to optimize the
capture wind energy by tracking the optimal torque T
opt
.
2
opt
T k = e , with
max
5
3
1
2
p
opt
C
k R = t
(7)
The standard strategy neglect the dynamic of the drive
train, which means that the aerodynamic torque is
supposed to be equal to the generator torque. It is obvious
that in many cases, and especially for turbulent winds,
this assumption will not be realistic.
Wind speed (m/s)
P
o
w
e
r
(
M
W
)
Wind power
C
p
= 1
Available power
C
p
= 0.593
(Betz limit)
Nonidealturbine power
C
p
= 0.412
Region I
Region II
Region III
P
o
w
e
r
(
M
W
)
Wind power
C
p
= 1
Available power
C
p
= 0.593
(Betz limit)
Nonidealturbine power
C
p
= 0.412
Region I
Region II
Region III
Wind speed (m/s)
P
o
w
e
r
(
M
W
)
Wind power
C
p
= 1
Available power
C
p
= 0.593
(Betz limit)
Nonidealturbine power
C
p
= 0.412
Region I
Region II
Region III
P
o
w
e
r
(
M
W
)
Wind power
C
p
= 1
Available power
C
p
= 0.593
(Betz limit)
Nonidealturbine power
C
p
= 0.412
Region I
Region II
Region III
Fig. 3. Steady-state power curves.
The proposed control strategy relies then on the
aerodynamic torque estimation using a high gain observer
[6-7]. This estimate is then used to derive a high-order
sliding mode controller that ensures T
opt
tracking in finite
time.
III.2. High Gain Observer
A high gain observer will therefore be used to estimate
the aerodynamic torque. A key feature of a high gain
observer is that it will reduce the chattering induced by a
sliding mode observer [6].
From (4), we have
a em
T T K
J J J
e
e = (8)
The following notations are introduced.
1
2
a
x
T
x
J
= e
(9)
Thus we have
1 2 1
2
( )
em
T K
x x x
x J J
x f t
(10)
or in matrix form
( , ) ( ) x Ax x u t
y Cx
= + + c
where
| |
1
0 1
0 0
1 0
( , )
0
0
( )
( )
A
C
Kx u
x u
J
t
f t
(
=
(
(
=
(
c =
(
.
A candidate observer could be [6]
1 1
( , ) ( )
T
x Ax x u S C C x x
u
= + uA
(11)
where
1 0
1 1
,
1
1 2 0
S
u
(
(
(
A = =
(
(
(
u
Let S be the unique solution of the algebraic Lyapunov
equation
0
T T
S A S SA C C + + = (12)
Let us define ( ) x x x
u
= A
h
a
l
-
0
0
6
3
0
4
3
6
,
v
e
r
s
i
o
n
2
-
4
N
o
v
2
0
1
1
then ( )
1
( ) ( ) ( ) ( )
T
x A S C C x x x t
u u
= u + A A c
(13)
Consider the quadratic function
T
V x Sx = (14)
then ( )
2
2 ( ) ( )
2 ( )
T T T
T
V xSx
V V x C Cx x S x x
x S t
u
u
=
= u + A
A c
Therefore
( )
max
( ) ( )
2 ( )
( )
x x
V V x S
t
u
u
| | A
| s u +
|
+ A c
\ .
(15)
We can assume that (triangular structure and the
Lipschitz assumption on )
( )
( ) ( )
( )
x x x
f t
u
A s
s o
with
K
J
=
It comes that
2
max max
2 ( ) 2 ( ) V V x S x S
o
s u + +
u
then
1 2
V V c V c V
o
s u + +
u
with
max
1
min
max
2
min
( )
2
( )
( )
2
( )
S
c
S
S
c
S
=
Now taking
{ }
( )
0 1
max
min
1
max
min 1
max 1,
( )
( )
2
( )
2
( )
c
S
S
c
S
M
S c
u
u
u =
u
=
and u > u
0
, we obtain
( ) exp( ) (0) e t t e M
u u
s u + o (16)
With
2
a
T Jx = , it comes that
exp( ) (0)
a a a
T T T J t e M
u u
= s u + o(
(17)
A practical estimate of the aerodynamic torque is then
obtained as M
= e e
(20)
then
T g
e FT G = +
(21)
Let us consider the following observer based on the
supertwisting algorithm [2], [8].
( )
( )
1
2
1
2
sgn
sgn
g T T
T
T y B e e
y B e
= +
= +
(22)
The gains B
1
and B
2
are chosen as
( )
( )
2
1
2 1 1 2
2
2
1 1
2
4
0
m
M
m
m M
B
A
B
A
G
F
u
>
u I + u
>
I u
< u
< I s < I
(23)
Thus we will guaranty the convergence of e
T
to 0 in a
finite time t
c
. The aerodynamic torque estimation is then
deduced.
T
a
= T
opt
, t > t
c
. (24)
The above proposed high gain observer principal is
illustrated by the block diagram in Fig. 4.
IV. Control of the DFIG-Based Wind Turbine
The DFIG-based WT control objective is to optimize
the extracted power by tracking the optimal torque T
opt
(7). The control is a compromise between conversion
efficiency and torque oscillation smoothing [9-11].
The reactive power is expressed as follows.
s sq sd sd sq
Q V I V I = (25)
h
a
l
-
0
0
6
3
0
4
3
6
,
v
e
r
s
i
o
n
2
-
4
N
o
v
2
0
1
1
2
nd
Order Sliding
Mode Controller
Wind Turbine
T
opt
c
T
P
g
T
g
+
_
e
r
Wind
2
opt opt r
T k = e
2
nd
Order Sliding
Mode Observer
T
a
^
2
nd
Order Sliding
Mode Controller
Wind Turbine
T
opt
c
T
P
g
T
g
+
_
e
r
Wind
2
opt opt r
T k = e
2
nd
Order Sliding
Mode Observer
2
nd
Order Sliding
Mode Observer
T
a
^
High Gain Observer
2
nd
Order Sliding
Mode Controller
Wind Turbine
T
opt
c
T
P
g
T
g
+
_
e
r
Wind
2
opt opt r
T k = e
2
nd
Order Sliding
Mode Observer
T
a
^
2
nd
Order Sliding
Mode Controller
Wind Turbine
T
opt
c
T
P
g
T
g
+
_
e
r
Wind
2
opt opt r
T k = e
2
nd
Order Sliding
Mode Observer
2
nd
Order Sliding
Mode Observer
T
a
^
High Gain Observer
Fig. 4. High gain observer principle.
Adapting (25) to our hypotheses, it comes then
s s s
s rd
s s
V V M
Q I
L L
|
= (26)
As the stator reactive power reference is zero, then
_
s s
s rd ref
s s
V V
I
M
| = =
e e
(27)
Let us consider the following tracking errors.
_
rd
em
I rd rd ref
T em ref
e I I
e T T
=
(28)
Then we will have
_
1
rd
em
sd
I rd r rd s r rq rd ref
r s
rq r rq s r rd
T s ref
s sd s r
s
d M
e V R I g L I I
L L dt
V R I g L i
M
e p T
M
g L L
L
| | |
= + e o
|
o
\ .
e o | |
|
= |
|
e | o
|
\ .
(29)
If we define the functions G
1
and G
2
as follows.
1 _
2
1
sd
s r rq rd ref
r s
s s r rd s sd ref
s r s
d M
G g L I I
L L dt
M M
G p g L I g T
L L L
| | |
= e o
|
o
\ .
| |
= | e o e |
|
o
\ .
Thus we have
1
2
1 1
rd
em
I rd r rd
r r
s rq s r rq
s r s r
e V G R I
L L
M M
e p V G p R I
L L L L
I
= +
o o
= | + + |
o o
(30)
Now, lets us consider the following high-order sliding
mode controller based on the supertwisting algorithm [2],
[10].
1
2
1 1
1 2
1
2
2 3
2 4
1
( )
( )
( )
( )
em em
em
rd rd
rd
rq T T r rd
r
T
rd I I s r rq
s r
I
V y B e Sgn e R I
L
y B Sgn e
M
V y B e Sgn e p R I
L L
y B Sgn e
= + +
o
= +
= |
(31)
Where the constants B
1
, B
2
, B
3
, B
4
, u
1
, and u
2
satisfy
(32).
Thus, we can assert that there exist finite times t
Tem
and
t
Ird
leading to (33).
This means that the control objective is achieved.
( )
( )
( )
( )
1 1
1 1 1 2
1 1 2
2 2
1 1
2 2
2 3 2 2 2
3 4
2 2
3 2
4
,
4
,
s
s r
r
r
r
G
B
M
B p B
L L
L B
G
B
B B
L
L B
< u
u + u
> | u >
o
o u
< u
u + u u
> >
o
o u
(32)
_
,
,
rd
em
rd ref rd I
ref em T
I I t t
T T t t
= >
= >
(33)
h
a
l
-
0
0
6
3
0
4
3
6
,
v
e
r
s
i
o
n
2
-
4
N
o
v
2
0
1
1
V. Validation Results with FAST
The proposed SOSM control strategy has been tested
for validation using the NREL FAST code [13]. An
interface has been developed between FAST and Matlab-
Simulink