DirectDrive Modeling
DirectDrive Modeling
I. I NTRODUCTION
e n
c y
e r t e r
G
G
r i d
r e q
C
r e q
C
S
G
e a r
u
v
e n
c y
e r t e r
r i d
2
P
g
e n
r e f
- I
r e f
e n
r e f
- I
i d
r e f
i q
r e f
r i d
c
g
e n
- I
- I
i q
e n
e r a t o
r e c t i f i e r
r - s i d
e
b
s t e r
g
i n
/ A
- s i d
e r t e r
i d
r e f
r e f
r e f
r e f
- I
r i d
D
e
i d
r e f
i q
r e f
e n
c o
/ A
e r a t o
n
- I
i q
r e f
i d
r i d
- I
i d
i q
m
m
r e f
r e f
i d
D
D
C
r - s i d
- I
e r t e r
c o
i q
/ A
r i d
- s i d
e r t e r
According to the equivalent circuits of figures 5 and 6 the statorvoltage equations can be expressed as follows:
"
!$#%'( & *)+ ( ,
, "
!$# , %-(& , %.+ (
(1)
0/. 1
/#2/3%4( & /
5 1
67#268%9( & 6
5 1
:#2:%4( & :
(2)
; + & <>=.?@%A=.B
C & <1A+<
(3)
+GIHJ+<
C IH C <
Vdc
UVACac
(4)
G
For completing the synchronous machine model the flux linkage equations are required.
Stator flux linkage:
(
"!%#%$&2#%#"!#2/3%#"! #26
(, "! , %'%$&2# , %#"! , #2:
C
v
a l v
a l v
(5)
Fig. 8. Valve representation by equivalent switches
( /.("!
( G
6 ("!
( :.("! ,
(6)
Automatic voltage regulators (AVR) of electrically excited machines can be considered by connecting the AVR-model to / .
Permanent magnet generators can be modeled with the above equations by keeping the excitation current # / to a constant value.
This synchronous machine model of seventh order is able to represent rotor and stator transients correctly. In stability studies however,
transient phenomena of the electrical network are usually not considered [10]. Applying the principle of neglecting stator transients to
the synchronous machine model reduces stator equations to arithmetic
equations:
1 "
!$#*) + ( ,
, "
!$# , %.+ (
(7)
Equation (7) allows expressing the electrical torque (4) by voltages and
currents only:
H , ,
,
,
=.B
+ / # %I# )I"!0#21 %I#21 43*6+5 < 7B7)879
:; :
;
!=< ?
> A @ <
(9)
>
The AC-voltage phase angle is defined by the PWM converter.
The pulse-width modulation index @ is5C
the
of the
B control
@ B variable
PWM converter. Equation (9) is valid for
5 . For values
larger than 1 the converter starts saturating and the level of low order
harmonics starts increasing. The complete characteristic of the PWM
converter, including the saturated range is shown in Fig. 9. Because of
the mentioned increase of low-order harmonics, the saturated are of a
PWM-converter is usually not used or just up to
.
(8)
PWM-Converter Characteristic
0,9
0,8
0,7
Grid-side converters are usually realized by self commutated pulsewidth modulated circuits (see Figure 7). The circuit is built from six
valves with turn-off capability and six antiparallel diodes. Valves with
turn-off capability are typically realized with bipolar transistors or IGBTs (insulated gate bi-polar transistors) because they allow for higher
switching frequencies than classical GTOs.
In older wind generators, line commutated, thyristor-converters
were sometimes used, but due to their high level of harmonics and the
lack of reactive power control capabilities there are not used in modern
designs.
Depending on the application, the circuit according to Figure 7 is
actually modeled using switches according to Figure 8 for representing
the valves.
|VAC |/VDC
0,6
B. PWM Converter
0,5
0,4
0,3
0,2
0,1
0
0
0,5
1,5
2,5
tion:
;
=! < != < E 5
V dc
7!=<,7 <
%
8 5
(10)
D. DC-booster
The average values of voltage and current pulse streams generated
by the dc-booster (see Figure 3) are as follows:
;
;
; 1 @;
1 1
V ac
Rl
(14)
(15)
B @ B @
5
. The cir-
The grid-side converter (see Figures 3 and 4) operates in a statorvoltage oriented reference frame. Hence, d-axis represents the active
and q-axis the reactive current component.
E 9
C. Diode Rectifier
I
d
243 5 6
ACB=D
2
243 5 6
V
d
V
c
78
;:
@ 356
V
d
I
c
V
0
<
+>
;
<
+
;
%
A5
7 356
8 !$#% &
( )1>?,
249 8
2 9
A very fast inner control loop regulates the d- and q-axis current
components of the grid-side PWM-converter. Current references are
defined by a slower outer control-loop regulating active- and reactive
power.
The objective of the MPT-characteristic defining the active power
reference is to drive the generator into the optimum speed-power operation point (optimum GF +H ) at any time.
As mentioned before, there are possible options for the active power
regulation shown in figure 13. In many designs, the active-current
reference is given by a speed controller, instead of the power controller.
In this case, a speed reference is calculated from the actual electrical
power using the inverse of the MPT-characteristic.
B. Generator-Side Converter
(11)
(12)
C. Booster Controller
(13)
The DC-booster of Figure 3 regulates the DC-voltage at the DCside of the grid-side PWM-controller. The concept is analogous to the
PWM-converter control concept. An inner loop regulates DC-current
and an outer voltage-controller loop defines the current reference.
<
.-
<
!=<
78356
0
Fig. 12. Equivalent circuit of the fundamental frequency diode rectifier model
a c
)
%$
&
`#bc d
!#"
'
(
!"
`#a
]
_
A. Wind Fluctuations
V. T URBINE
Equation (16) shows the aerodynamic equation of a wind turbine
that relates mechanical power to wind speed and mechanical speed of
the turbine (e.g. [5]):
(16)
with:
? : Mechanical power of the wind turbine
- : Air density
: Rotor radius
H
: Tip speed ratio
+ : Blade pitch angle
GF
coefficient as a function H and +
/ :: Power
Wind speed
B. Torsional Oscillations
When the simulated applications are limited to the impact of wind
fluctuations, it is usually sufficient to consider just a single-mass shaft
model because shaft oscillations of variable speed wind generators are
not reflected to the electrical grid due to the fast active power control
[5].
In stability analysis however, when the system response to heavy
disturbances is analyzed, the shaft must be approximated by at least
a two mass model. One mass represents the turbine inertia, the other
mass is equivalent to the generator inertia.
The equations describing the mechanical coupling of turbine and
generator through the gear box by a two-mass model can be expressed
as follows (see e.g. [9]):
; ?ik+Ej ? .
= ?) =.D
i+l <
; <Si =.G
D %A= B
i?l <
+Ej ? ) +<
i'm
i=.
l D>o
n ? < m ? <7%qp ? < +Ej ?) +<
E+ ?2
/
(17)
CE FG
0
ef g h
)
CD
QSRUTVNWX
Y Z[\
39
0 678 7,:
;=< 7>@?8@A,B
24
LNM
C M
Fig. 15. Generic wind turbine model
(19)
(20)
(21)
=.?E +E7 ??
j
IKJ
1 2 3 45
(18)
C. Protection
In case of heavy disturbances, direct drive synchronous machines
quickly disconnect from the system by blocking the PWM converter.
Usually, the following protection criteria a supervised:
r Under- and overvoltage
r Overcurrent
r Overspeed
With regard to under-voltage tripping, the tendency is more and more
to enable wind generators to stay on the network as long as possible,
even for deep voltage dips (voltage ride-through). With the grid-side
converter, reactive current support can be provided during voltage dips
for lifting the voltage at the generator terminals in case of faults near
to the connection point.
STABILITY
DIgSILENT
2.00
For investigating transient stability phenomena, large system models have to be studied considering the interaction of all generators in
the system. For reducing calculation time, simplified models of direct
drive synchronous generators should be used.
In stability analysis, the model should represent correctly:
r Electrical output at the connection point (P and Q depending on
voltage and angle).
r Generator and turbine acceleration.
-0.030 s
0.965 p.u.
1.00
-0.020 s
0.967 p.u.
0.00
-1.00
P
P
g
r e f
e n
- I
i d
r e f
r e f
Q
g
r i d
-2.00
-0.04
-0.02
-0.00
[s]
PWM Grid-Side Converter: Current, phase A in p.u., switched model
PWM Grid-Side Converter: Current, phase A in p.u., fundamental frequency EMT model
PWM Grid-Side Converter: Current, phase A in p.u., stability model
i q
r e f
0.02
1.25
i d
- I
i q
1.00
V
d
m
d
50.000 Hz
0.966 p.u.
0.75
D
c
/ A
C
1100.000 Hz
0.028 p.u.
0.50
g
c o
r i d
- s i d
900.000 Hz
0.036 p.u.
e
e r t e r
0.25
A simplified model can be derived under the assumption that the DCvoltage is constant. This assumption can be justified by the large capacitance in the DC-circuit and by the DC-voltage controllers, which
can be found in both concepts, the concept according to figure 3 and
figure 4.
With this assumption, only the grid-side converter needs to be represented fed by a DC-voltage source (see Figure 17).
Considering converter losses by a constant value (assumption of
constant dc-voltage), the electrical torque can be calculated by (8) using:
B )
(22)
< < 87
2000.00
4000.00
6000.00
PWM Grid-Side Converter: Current, phase A in p.u., frequency domain
[Hz]
8000.00
DIgSILENT
1.20
1.00
0.80
0.60
0.40
0.20
DIgSILENT
Breaker
0.00
-0.00
0.25
0.50
0.75
Point of Interconnection: Voltage, Magnitude in p.u.
Point of Interconnection: Voltage Phasor Magnitude in p.u., EMT-model
[s]
1.00
-0.00
0.25
0.50
PWM Grid-Side Converter: Active Power in MW
PWM Grid-Side Converter: Active Power in MW, EMT-model
0.75
[s]
1.00
-0.00
0.25
0.50
PWM Grid-Side Converter: Reactive Power in MVAr
PWM Grid-Side Converter: Reactive Power in MVAr, EMT-model
0.75
[s]
1.00
6.00
5.00
4.00
Point of Interconnection
3.00
2.00
1.00
AC Filter
0.00
PWM DC Capacitor
2.00
DC/DC Converter
1.50
1.00
DC Bus 1
0.50
DC Inductance
0.00
DC Bus 2
DC Capacitor
Fig. 20.
concept
G
~
D. Simplified model of converter driven synchronous wind generators for transient stability simulations
According to section VI, a simplified model was represented. The
electrical model of the generator and the DC circuit dynamics are neglected. However, the mechanical behavior of the generator (shaft,
aerodynamics and blade angle control) is represented and allows for
consideration of wind turbulences [8].
Figure 24 presents the results of the simulation of a voltage sag of
80% in the connection point with the simplified model. The resulting curves are compared against the results obtained with the detailed
model of the concept according to figure 3. The results match very
well. The behavior of the active and reactive power at the point of
1.275
DIgSILENT
DIgSILENT
1.20
1.00
1.250
0.80
0.60
1.225
0.40
0.20
1.200
0.00
1.175
-0.00
0.25
0.50
Point of Interconnection: Voltage, Magnitude in p.u.
Point of Interconnection: Voltage, Magnitude in p.u., EMT-model
0.75
[s]
1.00
-0.00
0.25
0.50
PWM Grid-Side: Active Power in MW
PWM Grid-Side: Active Power in MW, EMT-model
0.75
[s]
1.00
-0.00
0.25
0.50
PWM Grid-Side: Reactive Power in MVAr
PWM Grid-Side: Reactive Power in MVAr, EMT-model
0.75
[s]
1.00
6.00
5.00
1.150
0.00
1.25
2.50
Pitch Controller: Generator Speed
Pitch Controller: Generator Speed, EMT-model
3.75
[s]
5.00
4.00
3.00
8.00
2.00
1.00
6.00
0.00
4.00
1.00
0.70
2.00
0.40
0.10
-0.00
-0.20
-2.00
0.00
Fig. 21.
concept
1.25
2.50
Generator: Active Power in MW
Generator: Active Power in MW, EMT-model
3.75
[s]
5.00
VIII. C ONCLUSIONS
The modeling of converter driven synchronous machines commonly
applied to wind generation applications was treated in this paper. The
required representation of the different components particularly for
power system stability assessment was carefully described. Detailed
representation of these wind generator concepts was implemented and
simulation results were presented. A simplified model is also proposed
for large-scale power system analysis. The simulation results of this
model are compared to those obtained with the detailed models for
severe system disturbances. The proposed simplified model provides
adequate accuracy for transient and dynamic stability analysis and is
efficient with regard to calculation time.
R EFERENCES
[1] C. Ender, Wind Energy Use in Germany - Status 31.12.2002, DEWI
Magazin Nr. 22, February 2003
[2] Verband der Elektriztatswirtschaft (VDEW), Stromerzeugung in
Deutschland 2002,17.02.2003, http://www.strom.de
[3] S. Heier Grid Integration of Wind Energy Conversion Systems John
Wiley and Sons, 1999
[4] J. L. Rodriguez-Amenedo, Automatic Generation Control of a Wind
Farm With Variable Speed Wind Turbines, IEEE Transactions on Energy Conversion, Vol 17, No. 2, June 2002
-0.50
Fig. 22. Voltage sag simulation of the wind generator concept with two PWM
converters
DIgSILENT
1.275
1.250
1.225
1.200
1.175
0.00
1.25
2.50
Speed measurement: speed
Speed measurement: speed, EMT-model
3.75
[s]
5.00
DIgSILENT
1.150
1.20
1.00
0.80
6.00
0.60
0.40
5.00
0.20
0.00
4.00
-0.00
0.50
1.00
Point of Interconnection: Voltage, Magnitude in p.u.
1.50
[s]
2.00
-0.00
0.50
1.00
PWM Grid Side: Active Power in MW (Simplified Model)
PWM Grid-Side Converter: Active Power in MW (Detailed Model)
1.50
[s]
2.00
-0.00
0.50
1.00
1.50
PWM Grid Side: Reactive Power in MVar (Simplified Model)
PWM Grid-Side Converter: Reactive Power in MVar (Detailed Model)
[s]
2.00
-0.00
0.50
Simplified Generator: speed in p.u.
Generator: speed in p.u.
[s]
2.00
6.00
3.00
5.00
4.00
3.00
2.00
2.00
1.00
1.00
0.00
0.00
0.00
1.25
2.50
Permanent Magnet Generator: Active Power in MW
Permanent Magnet Generator: Active Power in MW, EMT-model
3.75
[s]
5.00
3.50
3.00
2.50
2.00
Fig. 23. Voltage sag simulation of the wind generator concept with two PWM
converters
1.50
1.00
0.50
B IOGRAPHIES
Sebastian Achilles received the Electrical Engineer degree from
the University of La Plata, Argentina and the MS degree in Electric
Power Eng. from the Rensselaer Polytechnic Institute in Troy, NY in
1996 and 2000 respectively. He worked on power systems consulting
activities for ABB T&D Systems (Spain), GE Power Systems Energy
Consulting (NY, USA) and for the Instituto de Investigaciones
Tecnologicas de Redes y Equipos Electricos (Argentina). He currently
works in the Consulting Department of DIGSILENT GmbH carrying
out training and consulting activities on several topics including
integration of wind generation to power systems, power system
stability, network planning, test and modeling of generator controls
and insulation coordination
Markus Poller received Dipl.-Ing. degrees from the University
of Stuttgart and Ecole Nationale Superieure des Telecommunications
Paris in 1995. In 2000, he received the Dr.-Ing. degree from the
University of Hannover. Since 1995 he works with DIgSILENT
GmbH, Germany, where he holds the position of Manager Research
& Development. His activities involve the development of models
and algorithms for the power system analysis program DIgSILENT
PowerFactory, power system studies and the presentation of softwareand power system analysis courses. His current research interests
include wind power systems, optimal power flow dispatch and
probabilistic load flow analysis.
1.275
1.250
1.225
1.200
1.175
1.150
1.00
1.50