Influence of the representation of the distribution transformer core configuration

on voltages developed during unbalanced operations

J. Vernieri B. Barbieri P. Arnera
Instituto de Investigaciones IITREE IITREE
Tecnolgicas para Redes y Equipos UNLP. UNLP.
Elctricos (IITREE)-Universidad Calle 48 y 116. (1900) La Plata, Calle 48 y 116. (1900) La Plata,
Nacional La Plata (UNLP). Argentina. Argentina.
Calle 48 y 116. (1900) La Plata. [email protected] [email protected]
Argentina
[email protected]

Abstract - The object of this paper is to emphasize the three-legged core type or the five-legged core type
importance of the representation of the core available models.
configuration when seeking transients developed in a
distribution transformer under an unbalanced 3-LEG CORE TYPE TRANSFORMER CIRCUIT
situation. TRANSFORMER MAIN BREAKER BREAKER

Digital time-domain simulations are performed with

the Alternative Transients Program (ATP) for the case UNDERGROUND
CABLE
of a distribution three-phase transformer using the
following two core configuration representations
3 PHASE
available in the ATP: three-legged core form and five- LOAD SOURCE
legged core.
Severe overvoltages and high saturated currents would Fig. 1 - Simplified circuit representation.
not be evidenced if the wrong core configuration
representation was used during simulations. The event performed during the simulations consisted
Keywords: Transient Analysis, Transformer, ATP. of opening two poles of the circuit breaker located at the
opposite end of the cable, as depicted in the simplified
I. INTRODUCTION circuit of Fig.1.
Transformer load was considered disconnected. Table 1
Transformer core configuration representation becomes shows the transformer specifications data, and Table 2
important especially in cases of unbalanced conditions in shows the positive sequence saturation curve obtained
which zero sequence parameters are excited. from the excitation test (LV side).
In a three-legged core type transformer zero sequence
fluxes cannot circulate exclusively in the core. They are Table 1. Transformer specification data.
forced to do it outside the windings through the air gap, HV: High Voltage
structural steel and the tank. As a consequence a higher LV: Low Voltage
reluctance is seen by this flux, and the resulting zero Voltages Rating (VHV/VLV): 132 kV/13.86 kV
Power Rating = 40 MVA
sequence magnetising inductance will be less than that for
Winding connection (HV/LV): YNyn0 grounded-grounded
the positive sequence. An appropriate representation of Positive sequence parameters:
such behaviour should be required. Excitation losses = 24 kW
On the contrary, a five-legged core type transformer Excitation current = 0.15 %
provides an iron-core path where the homopolar flux can Short-circuit losses = 176 kW
close. Therefore, for such configuration it is reasonable to Short-circuit impedance = 13.3 %
assume the zero sequence parameters equal to the positive Zero sequence parameters:
sequence ones. Excitation losses (at excitation voltage)= 2.1 kW
Excitation current LV side (3 I0)= 77 A
Excitation voltage LV side = 140 V
II. DESCRIPTION
Table 2. Positive sequence saturation curve.
The case under study consists of an unloaded U (p.u.) ULINE rms (kV) I rms (A)
distribution transformer fed through an underground cable, 0.96 13.30 1.80
and under unbalanced excitation conditions. 0.97 13.48 2.00
The transformer considered is a two winding, three- 1.00 13.86 2.50
legged core type, three-phase, 40 MVA, 132 kV/13.8 kV 1.02 14.20 3.20
transformer one. With both windings WYE-connected, and 1.05 14.60 4.30
with both STAR points grounded.
1.08 15.00 6.45
In order to reveal the influence of the core
1.09 15.10 7.40
configuration representation in the resulting voltage and
current waveforms, some simulations were performed with 1.10 15.25 9.55
the Alternative Transients Program (ATP) using either the 1.158 16.05 22.53
 A 13.5 km in length underground cable, with a positive magnetising inductance, the same represented in CASE A,
sequence capacitance of 0.404 F/km, and a three phase is indicated as SATURA in Fig. 2, and it is connected
equivalent source with its inductance, represent the rest of through the LV side winding in each phase.
the 132 kV system. Besides, the zero sequence excitation losses are very
different from those for the positive sequence. As a
III. TRANSFORMER REPRESENTATION consequence the selected model uses three linear mutually
coupled resistances (RMAG), externally connected in
The transformer was modelled with the SATURABLE parallel with the non-linear magnetising inductance, as
TRANSFORMER component, which consists of a three- indicated in Fig. 2, to represent positive and zero sequence
single phase units bank. Magnetic circuit saturation effect losses.
is represented with a non-linear magnetising inductance
connected phase to ground on the LV windings of each IV. SIMULATIONS
phase. Magnetic coupling between phases is not
represented. Simulations have been performed for the two analyzed
cases. Circuit breaker opens poles B and C at t = 0.1 sec.,
CASE A. Five-legged core type
while phase A remains energized.
For the five-legged core type configuration the zero The response of the cable-transformer system to de-
sequence parameters are supposed to be the same as the energising two phases can be determined by observing the
positive sequence ones. The model represents the non- voltage waveform at the open phases.
linear positive sequence magnetising inductance as
connected phase to ground on the LV windings of each CASE A: Results
phase. At the linear region its value is Lm = 1,242 H seen
Since the model does not represent magnetic coupling
from HV side. Excitation losses are represented internally
among phases, after the poles of the circuit breaker (Fig.1)
through a linear resistance parallel to the magnetising
are opened the cable becomes isolated from the excitation
reactance in each phase.
source in these two phases. As a consequence the trapped
CASE B. Three-legged core type charge in the cable starts to decay through the transformer
oscillating at a frequency that it is influenced by the
Since the homopolar reluctance is very high because of
saturation characteristic of the magnetizing inductance.
the air gap where the flux is forced to close, the zero
Fig. 3 shows phases B and C currents flowing from the
sequence magnetising curve is nearly linear, and a totally
cable to the transformer. As the core reaches its saturation
different model is needed. The ATP library offers a model
condition, magnetizing current climbs up to several times
which represents the equivalent circuit depicted in Fig. 2,
the transformer rated excitation current (110 A peak
where the zero sequence magnetising curve is represented
value). Anyway the current peak values are lower than the
through the linear magnetising inductance L0 (1.483 H
transformer rated peak current (247 A).
from HV side). The non-linear positive sequence

Fig. 2 - Equivalent circuit representing a two-winding, 3-leg core-type, three-phase transformer.

 Fig. 4 shows line-to-ground voltage waveforms in 60
phases B and C, which always remain below rated voltage,
[A]
and decrease as the trapped charge decays.
40

Current (A)
20
120

[A]

80 0

40 -20
Current (A)

0 -40

-40 -60
0.0 0.2 0.4 0.6 0.8 t[s] 1.0
(file CASE-A.pl4; x-var t) c:XX0028-
-80
Time (s)
Fig. 5 - CASE A - Earth current from the transformer
-120
HV side.
0.0 0.2 0.4 0.6 0.8 t[s] 1.0
(file CASE-A.pl4; x-var t) c:SBSTB -HVB c:SBSTC -HVC Hence after opening the two poles of the circuit
Time (s)
breaker, voltage applied to phase A results in magnetically
induced voltages in the two apparently de-energised phases
Fig. 3 - CASE A - Phases B and C currents. B and C.
The magnitudes of the line-to-ground voltages induced
in phases B and C depend on the circuit parameters, i.e. the
zero sequence transformer magnetising inductance, the
120 underground cable capacitance, the transformer losses and
[kV] the load.
80
Line-to-ground voltage waveforms in phases B and C
reach values that exceed the rated voltage for more than
four times (440 kV peak value > 4 p.u. value), as depicted
40
in Fig. 7. It can also be observed that the voltage
Voltage ( kV)

waveforms are modulated by a low frequency.

0 Fig. 8 shows phases A, B and C current waveforms in
the transformer. As the core reaches its saturation
-40 condition magnetizing currents climb up to several times
transformer rated excitation current, moreover they exceed
-80
for more than two and a half times the transformer rated
peak current (650 A peak value for phases B against 247 A
transformer rated peak current).
-120
0.0 0.2 0.4 0.6 0.8 1.0
Fig. 9 shows phases B and C magnetizing current
t[s]
waveform, and Fig. 10 shows the current flowing to earth
(file CASE-A.pl4; x-var t) v:SBSTB v:SBSTC Time (s)
from the neutral point, all of them considered on the HV
side of the transformer.
Fig. 4 - CASE A - Phases B and C line-to-ground
voltage waveforms.

X Eq
The current flowing to earth from the neutral point of EA
the transformer HV side can be seen in Fig. 5. Lm Cc UA
Current and voltage waveforms show how the core gets
in and out from the saturation region. 1:1

Cc UB
CASE B: Results Lm 3 L0

Although the model is not intended to represent the 1:1

existing magnetic coupling among phases the equivalent
Cc UC
electric circuit depicted in Fig. 6 provides a path for the Lm
zero sequence current that, in some way, reproduces the
coupling among phases.
Fig. 6 - CASE B - Simplified equivalent circuit.
 120
phase A 700 phase A
[A]
[kV]
525
80

350

40
Voltage ( kV)

Current ( A)
175

0 0

-175
-40

-350
-80
-525

-120 -700
0.0 0.2 0.4 0.6 0.8 t[s] 1.0 0.0 0.2 0.4 0.6 0.8 t[s] 1.0
(file case-b.pl4; x-var t) v:SBSTA
Time (s) (file case-b.pl4; x-var t) c:SBSTA -HVA Time (s)

500 phase B 700 phase B

[kV] [A]
375 525

250 350

Current ( A)
125 175
Voltage ( kV)

0 0

-125 -175

-250 -350

-375 -525

-500 -700
0.0 0.2 0.4 0.6 0.8 t[s] 1.0 0.0 0.2 0.4 0.6 0.8 t[s] 1.0
(file case-b.pl4; x-var t) v:SBSTB Time (s) (file case-b.pl4; x-var t) c:SBSTB -HVB Time (s)

500 phase C 700 phase C

[kV] [A]
375 525

250 350
Current ( A)

125 175
Voltage ( kV)

0
0

-175
-125

-350
-250

-525
-375

-700
-500 0.0 0.2 0.4 0.6 0.8 1.0
t[s]
0.0 0.2 0.4 0.6 0.8 t[s] 1.0
(file case-b.pl4; x-var t) c:SBSTC-HVC
(file case-b.pl4; x-var t) v:SBSTC Time (s) Time (s)

Fig. 7 - CASE B - Phases A, B and C line-to-ground Fig. 8 - CASE B - Phases A, B and C current
voltage waveforms. waveforms.
 250.0 phase B
[A]
To understand the cause of such high voltages and the
187.5 presence of the low frequency modulation the equivalent
circuit shown in Fig. 6 will be analyzed. Transformer
125.0 losses and transformer leakage inductance have not been
depicted in the figure, in order to make it simpler.
62.5
Current ( A)

While the voltage drop across the phase A equivalent

0.0 impedance XEq is negligible, the excitation voltage EA can
be considered as directly applied to the direct sequence
-62.5 transformer phase A non-linear magnetising inductance
(Lm) in parallel with the phase A cable capacitance (Cc),
-125.0
and in parallel to a series-circuit which includes three times
-187.5 the zero sequence transformer magnetising inductance
(3L0), as well as the two other phases shunt circuits.
-250.0 The two shunt circuits, phases B and C, constituted by
0.0 0.2 0.4 0.6 0.8 t[s] 1.0
(file CASE-B.pl4; x-var t) c:X0028B-
the direct sequence magnetising inductance and the cable
factors: 1 0,105 capacitance, have a resonance frequency at 1.94 Hz. This
offsets: 0,00E+00 0,00E+00 Time (s) frequency is responsible for the low frequency modulation
seen in voltages and current waveforms of Fig. 7 and 8.
250.0 phase C
Under the excitation frequency (50 Hz) and at the
[A]
187.5
linear region the dominant characteristic of the equivalent
shunt circuit is the cable capacitance. Consequently, under
125.0 such conditions the equivalent circuit of Fig. 6 can be
reduced to a series circuit constituted by the two
62.5
capacitance (phases B and C), three times the zero
Current ( A)

0.0
sequence transformer magnetising inductance and the
excitation voltage, considering the following fasorial
-62.5 relation (1).
-125.0
EA = U3L0 UB UC (1)
-187.5
where EA is the equivalent source terminal excitation
-250.0 voltage on phase A, U3L0 is the voltage drop across the
0.0 0.2 0.4 0.6 0.8 t[s] 1.0
three-time zero sequence transformer magnetising
(file CASE-B.pl4; x-var t) c:X0028C-
factors: 1 0,105 inductance, UB and UC are the line-to-ground voltages at
offsets: 0,00E+00 0,00E+00 Time (s) the HV transformer phase B and C terminals, respectively.
This series circuit has a resonance frequency near
Fig. 9 - CASE B - Phases B and C magnetising current 45 Hz, which can be determined using (2) and available
waveforms reflected on HV side. data.
1
f = (2)
Cc
2 3L0 2
1600
[A] where f is the resonance frequency of the series equivalent
1200
circuit.
Since at the excitation frequency, 50 Hz, this circuit is
800
close to its oscillation frequency it could be expected
important overvoltages in the middle points of the series
Current ( A)

400
LC circuit. Furthermore, these line-to ground overvoltages
0 UB and UC, applied to the non-linear magnetizing
inductance (Lm) in phases B and C cause core saturation,
-400 and consequently increase the overvoltages even more.
The core saturation condition can be observed through
-800 the current waveforms of phases B and C magnetizing
current, shown in Fig. 9.
-1200

-1600
0.0 0.2 0.4 0.6 0.8 t[s] 1.0 V. CONCLUSIONS
(file case-b.pl4; x-var t) c:XX0030- Time (s)
When the distribution transformer core configuration
Fig. 10 - CASE B - Earth current from the transformer HV influence is underestimated, choosing a five legged core
side. configuration representation instead of a three legged one,
could lead to overlook overvoltages and saturation additional protective action, may affect the transformer
conditions that in a real distribution transformer with a insulation.
three-legged core type may occur under unbalanced A real case confirms the overvoltages predicted
excitation. considering that serious damage were registered in a
Although neither the five-legged nor the three-legged lightly loaded three-phase three-legged core form
core type configuration selected models represent magnetic distribution transformer when one pole remained closed
coupling among phases properly, the equivalent electric after an anomalous three-phase circuit breaker opening
circuit employed by the second one reproduces, in some operation.
way, through the zero sequence magnetising inductance To sum up, the representation of the special behaviour
the magnetic coupling among phases. Therefore this model of the three-legged core configuration under zero sequence
is more appropriate to make evident the overvoltages that excitation conditions becomes essential for the system
could occur on the other two phases involving the cable being analysed. The three-legged representation available
capacitance and the non-linear behaviour of the core. from the ATP seems to show, in this case, a better
The five-legged core type available model definitely performance to reveal the abnormal conditions that stress
fails to reproduce voltages induced over the two other the transformer, than the five-legged representation also
phases as well as the non-linear phenomena developed on available in the ATP.
such phases.
In the case under study several different circumstances
appeared and were responsible for the overvoltages VI. REFERENCES
developed over the two apparently de-energised
transformer phases. First of all it is worth to be mentioned [1] R. Iravani et al., Modelling Guidelines for Low
the unbalance condition that excites the low zero sequence Frequency Transients, IEEE PES Special
magnetising inductance that this kind of core configuration Publication Modelling and Analysis of System
presents. Then it is also important to mention the series LC Transients 99TP133-0, 1998, pp.3-1/3-29.
circuit near resonant conditions which result in high [2] Phillipe Ferracci, Ferrorresonance, Cahier
voltages exciting the non-linear behaviour of the Technique Schneider n190, March 1998.
transformer core. Finally the lack of damping, cause of [3] Bonneville Power Administration, Electromagnetic
long permanence of high voltages that, without any Transients Program Rule Book, Portland, Oregon,
April 1982.