Model Development and Validation of Brushless Exciters
Model Development and Validation of Brushless Exciters
Exciters
Tony Bertes
DIgSILENT Pacific
Melbourne, Australia 3004
Email: [email protected]
AbstractThis paper describes the process of performing site
tests and model development of brushless exciters. Brushless
exciters are often a source for model inaccuracy when studying
synchronous machines and excitation system, which may lead to
incorrect stabiliser design and poor performance. The process described here involves on site testing of the excitation system both
with the generating unit operating open circuit and synchronised
to the network which will enable the extraction of parameters to
fit a model for the brushless exciter based on the IEEE Std 421.5.
Two case studies shall be presented to contrast the difference in
approach given the different limitations faced at each site.
KeywordsBrushless exciter, model validation, excitation system, AVR, synchronous machine, stability, field testing.
I NTRODUCTION
For a number of years, manufacturers of AVRs often provided IEEE models for use in power stability studies. However
there has been a gradual movement towards detailed models for
each type of AVR that is on the market. This has largely been
assisted by the improvement of control system implementation
in digital based hardware which enables greater transparency
of how the regulator is programmed and behaves under various
conditions.
Unfortunately this is not the case with exciter machines.
The inaccuracies of exciter models are often amplified during
a refurbishment project. For example, a generator (asset owner)
II.
T HE B RUSHLESS E XCITER
1
sTE
vE
fEX
VE,min
iFD
v X = v E SE (v E )
fEX = f (i N )
KE
K i'
i N = C FD
vE
vX
IN
i'FD
KD
Figure 1.
e'FD
Renormalization
I.
GEM ,11(s ) =
e'FD (s )
KVE
=
,
eFD (s ) K E + KVx + sTE
KVx =
v x
,
v E
KVE =
i N fEX (fEX v E )
v E i N fEX
Figure 2.
C. Comparison of Exciters
iFD = eFD
1
TE
1 + s. KE+SE
(1)
~
SM
Rotating exciter
B. Rectifier Model
The rectifier conduction introduces non-linear effects, depending on the field current (iFD ) and the applied voltage, VE.
These non-linear effects have been calculated and represented
in a three step linearised function, FEX . This function is a
simplification of a non-linear response but experience with the
widely used approximation has been such that it is seldom
modelled to a higher level of complexity. It is also not
uncommon to have KC set to zero, which introduces further
approximations into the model.
All ac sources that supply rectifier circuits have an internal
impedance that is predominately inductive. This impedance has
an effect in the commutation process and produces a non linear
characteristic as the rectifier load current and voltage varies.
SM
Static excitation
Figure 3. Comparison of the two main exciter types (a) the brushless exciter,
and (b) the static exciter
Brushless Excitation
Static Excitation
-40
-60
-80
-100
-120
-140
-160
-180
-200
-1
10
10
Frequency (Hz)
10
DIgSILENT
-20
III.
1.022
DIgSILENT
1.008
1.002
0.995
0.988
10.00
14.40
Measured: Measured Terminal Voltage [p.u.]
G6: Terminal Voltage in p.u.
Measured: Upper Band
Measured: Lower Band
18.80
23.20
27.60
[s]
32.00
A. Field Testing
At this particular installation, the excitation system comprised of:
1.045
1.030
1.015
1.000
0.985
0.000
1.015
3.200
sym_30451_1: Static - Terminal Voltage
sym_30451_1: Brushless - Terminal Voltage
6.400
9.600
12.80
[s]
16.00
Output
Input
Set Point
Error
+
-
Vr
Controller
(Regulator)
Efd
Terminal Voltage
Plant
(Generator)
Power Amplifier
(Exciter)
Terminal Voltage
Feedback Element
Figure 7.
Figure 9.
550
Figure 10 shows the composite model used for the Parameter Identification simulations. A window function was used
to ensure the Parameter Identification calculation commenced
only when a transient event occurs.
500
450
Gain
400
350
300
250
200
150
100
1
10
Frequency (Hz)
Simulated AVR Gain
DIgSILENT
0.1
40
35
Compare Block
30
25
Measured Data
Phase
20
15
in1mea
Sig
*
10
in2mea
0
-5
-10
0.1
in1sim
sym Slot
ElmSym*
10
Compare
ElmComp*
Frequency (Hz)
Simulated Data
in2sim
avr Slot
ElmAvr*,ElmVc..
4
qzpf
window
*
Window Function
B. Exciter saturation
The exciter saturation function SE(eFD ) is defined as
a multiplier of pu exciter output voltage to represent the
increase in exciter excitation requirements due to saturation.
To accurately derive the exciter saturation, the exciter output
voltage needs to be measured and plotted against exciter field
current. Unfortunately, the exciter output voltage could not be
measured and the exciter saturation could not be derived from
site testing. Typical data has been assumed from similar sized
Brush machines that allows for a small amount of saturation
when the exciter is loaded.
Figure 10.
3)
4)
5)
DIgSILENT
1.033
Due to the extensive testing, the process of model development did not rely on Parameter Identification, but rather on
deriving values for the exciter parameters from first principles.
The model was therefore developed and verified in the time
and frequency domain using the results from site.
1.026
1.019
1.012
1.005
12.00
IV.
14.00
sym_30451_1: Terminal Voltage in p.u.
Signal Model: Measured Terminal Votage
Signal Model: UgLower
Signal Model: UgUpper
16.00
18.00
20.00
[s]
22.00
With the availability of the mentioned signals, each component of the excitation system could be theoretically be tested
and verified independently. In this case, the AVR is planned
to be replaced by a new digital based system and could be
considered as a black box as a model for the existing AVR
is not required to be submitted as part of the connection of
the planned alteration. However, for closed loop simulations a
reasonably accurate model is required. To verify the adequacy
of the model for the regulator, a transfer function measurement
was made. Figure 12 shows the frequency response of the AVR
10
5
0
-5
Gain (dB)
Gain (dB)
-15
-20
Simulated
Measured
13
-25
12
-30
11
-35
-2
10
-1
10
10
10
10
9
20
8
7
Measured
Simulated
6
-1
10
10
Phase Angle
10
-20
-40
-60
-5
Phase Angle
-10
-80
-10
-100
-2
10
-1
10
10
10
Frequency (Hz)
-15
-20
Figure 13. Measured transfer function of the exciter compared with the model
-25
-30
-1
10
Figure 12.
the model
10
Frequency (Hz)
10
The frequency response of the model shows good agreement with the measured transfer function, which suggests that
the developed model is suitable for use in large scale stability
studies, including the design of a stabiliser.
B. Resistance measurements
The resistance measurements of the generator stator, rotor
and exciter field winding provide valuable information to be
able to assess the steady state performance of the model. As
the measurements were taken at standstill, the temperature of
the winding under test must be considered to be able to correct
to rated temperature.
For all three tests, the resistance was measured three
times, both in the forward and reverse direction. For the rotor
resistance, RFD the per unitisation occurs as follows:
ZF Dbase =
RF D =
3phV Abase
i02
F Dbase
M easuredResistance
ZF Dbase
(2)
(3)
The rotor resistance is of particular importance particularly when deriving steady state generator field current from
measurements of field voltage (from Ohms law), and also
understanding the affect it has on the generators (and exciters)
time constant. It is also of significance when determining the
saturation characteristic of the exciter and the synchronous
generator.
C. Frequency response
The frequency response of the exciter was measured during
unsynchronised operation. The unit was operated below 0.8 pu
terminal voltage to minimize the effects of saturation in the
exciter or in the main machine during the test. Figure 13 show
the measured and corresponding simulated frequency response
of the ac exciter. The simulated response was calculated using
the validated settings, and corrected to test temperature.
(4)
TransGrid recommended that a scaling factor to the normally calculated value of Fex may be required to
meet the steady state accuracy of exciter field current. An outcome of determining exciter field current
base was to apply a scaling factor of 0.94 to the Fex value.
Based on the data presented, the base exciter field current is determined as 8.84 A that is 8.84 A exciter
field current for 1.0 p.u generator terminal voltage on the air gap line.
10
The response of the rotor voltage to this sudden deexcitation and discharge is shown in Figure 16.
2
1.00
X = 13.370 s
DIgSILENT
X = 13.858 s
-2
0.75
-4
-6
0.50
-8
0.341
0.337
Y = 0.337
-10
0
50
100
150
200
Measured rotor voltage (V)
250
0.25
300
Figure 4: Percentage error in steady exciter field current values against measured rotor voltage
DIgSILENT
Figure 14. Percentage error in steady exciter field current values against
measured generator field voltage
1000.0
[-]
800.00
0.00
-0.25
13.00
13.40
Signals: Measured Rotor Voltage in p.u.
G3: Simulated Rotor Voltage in p.u.
13.80
14.20
14.60
[s]
15.00
Figure 16. De-excitation test of the field to determine the open loop time
constant of the exciter
600.00
400.00
200.00
6
0.0000
0.0000
8.0000
16.000
avr_T1G3: Exciter Field Current (A) / Exciter Voltage (V) Air Gap
avr_T1G3: Exciter Field Current (A) / Exciter Voltage (V) with Saturation
Figure 15.
24.000
32.000
[-]
40.000
G. Dynamic performance
The dynamic performance of the model in time domain
was verified by performing unsynchronised and synchronised
voltage step response simulations and comparing them to the
actual plant response. Figure 17 compares the responses of
the model and the actual plant response to a synchronised
voltage step response. Good agreement is obtained between
the simulated and measured terminal voltage.
In addition to the step response tests, the exciter model
was further validated by open loop simulations. This involved
injecting measured excitation voltage (treating the AVR as a
black box) to the field of the exciter model and comparing
the simulated rotor voltage with the measured data. Figure 18
compares the open loop response of the exciter model and the
actual plant response to a synchronised voltage step response.
DIgSILENT
1.095
1.087
1.079
1.071
ACKNOWLEDGMENT
The author would like to thank the staff at DIgSILENT
Pacific, in particularly Tim George and David Browne for their
assistance in the model validation work and support, and also
Ron Oosterwijk for his assistance in field testing and support.
1.063
1.055
12.00
14.00
Signals: Measured Terminal Voltage in p.u.
G3: Terminal Voltage in p.u.
Signals: UgLower
Signals: UgUpper
16.00
18.00
20.00
[s]
22.00
R EFERENCES
DIgSILENT
Figure 17. Synchronised voltage step response to a 2.5% step applied to the
AVR of Hydro, overlayed with actual plant response
1.40
[1]
[2]
1.20
[3]
1.00
[4]
0.80
0.60
[5]
0.40
10.00
12.80
15.60
iexc file: Measured Excitation Current in p.u.
avr_T1G3: Simulated Excitation Current (Vfe) in p.u.
18.40
21.20
[s]
24.00
12.80
Signals: Measured Rotor Voltage in p.u.
sym_20843_3: Simulated Rotor Voltage in p.u.
18.40
21.20
[s]
24.00
2.10
1.70
1.30
0.90
0.50
0.10
10.00
15.60
Figure 18. Open loop response to a 2.5% step applied to the AVR of Hydro,
overlayed with actual plant response
V.
C ONCLUSION