Open Phase Conditions in Transformers - Analysis and Protection Algorithm
Open Phase Conditions in Transformers - Analysis and Protection Algorithm
I. I NTRODUCTION
Open phase condition frequently occurs in power systems,
particularly in medium voltage distribution networks (such
as 33KV or below) in rural or remote areas where broken
conductor, loose connection or a blown fuse can lead to an
open phase. Some other cases of an open phase may involve
a defective contact of a circuit breaker. This is considered
a series fault with a different fault current calculation and
sequence networks connection than other known shunt faults
such as LG or LLG. In addition, an open phase on one side of
a transformer can be very difficult to detect on the other side
depending on the transformers core and winding configuration
as well as loading level.
In medium voltage distribution systems, industrial plants,
and rural areas and for economic reasons often smaller transformers, such as 5MVA or below, are protected mostly on the
load (low) side with the possibility of only fuses on the high
side usually provided by the local utility. As will be shown,
with an open phase on the high side of a transformer there
are cases that the three phase voltages and currents on the low
side could be fairly balanced, particularly in lower loading
conditions, such that neither overcurrent nor negative sequence
current elements can pick up a fault. In fact in lower loading
conditions, even as high as 50% , the transformer may continue
normal operation for hours and possibly longer with acceptable
balanced voltages and currents on the low side while there
is an open phase on the high side. Only in higher loading
conditions may the protection relays pick up a fault, usually
on negative sequence or sensitive ground overcurrent.
An open phase anomaly can also occur in distributed
generation (DG). Many small DG sites, including landfill,
Fig. 2.
Fig. 3.
Fig. 1.
Fig. 4.
Vm
2Vrms
=
sin t =
sin t
N1
2f N1
(2)
Therefore
magnetic flux, , is sinusoidal with an amplitude
2Vrms
of m = 2f
N1 . We can write equation 2 in the following
well-known format:
Vrms = 4.44N1 f m
(3)
Fig. 5.
(5)
1
1
VA = VC = VBC = VBC 6 180
2
2
(6)
IA = 0; IC = IB
(7)
Fig. 7.
(8)
1
VBC
2
(9)
IA = 0; IC = IB
However, in other types of transformer cores upon open
phase conditions either a fraction of B and C will pass
through winding A (4-leg, 5-leg, shell type) or there is no flux
linkage between the phases at all (three single phase bank).
In either case there is very little or no voltage induced in the
winding with lost phase.
As illustrated in Figure 6 and regardless of the core type,
upon loss of phase A the ground current at the neutral point
of the transformer is per below equation. We will use this
equation later in the discussions on the open phase condition
in transformers.
IG = 3I0 = IB + IC
(4)
(10)
Fig. 9.
Fig. 8.
(11)
Fig. 10.
Further simulation results show that on the low side II21 and VV21
are about 0.62% which demonstrate a quite balanced system.
Transformer LV Voltages- Vab, Vbc, Vca
1.05
0.036
pu
pu
0.034
1
0.032
0.95
0.03
0.05
0.1
0.15
0.2
1.05
pu
0.06
pu
0.04
0.02
0.95
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
1.1
0.4
pu
pu
1
0.3
0.9
0.8
0.2
1.1
0.6
pu
pu
0.9
0.8
0.4
0.2
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
0.3
pu
0.2
0.1
0.05
0.1
Time (Sec)
0.15
0.2
Fig. 13.
Coil Currents on the Secondary Side After Loss of Phase ALoad=600KW
IA = 0
I2
= 16 180
I1
(13)
I2
= 16 60
I1
(14)
I2
= 16 60
I1
(15)
6
I1Sec. = I1P ri. 60 as I1P ri. and I2P ri will transfer to the
secondary side with the same ratio but opposite phase angle
shifts.
I
We can now establish the following values for I21Sec. , as in
Sec.
Table A, that apply to negative and positive sequence currents
on the low side of a transformer with an ungrounded primary
winding of or Y upon loss of a single phase on the high
side. This is regardless of the core construction and secondary
winding configuration. It also applies to DG cases where there
is a generation source on the low side exporting power to the
grid. All we need to know is the phase angle shift of low side in
relation to the high side. In each case high side equations 13,
14 and 15 are re-written for the low side with appropriate
phase angle shift.
Table A: II12 Values on the Transformers Low side
with an Open Phase on the Ungrounded High side
Case
A Open
B Open
C Open
1 180
16 60
16 60
1 120
16 120
16 0
1 120
16 0
6
1 120
3) Yg -Yg Transformers: So far weve seen that, regardless of the core construction, in Yg - transformers the lost
phase is always re-generated while in transformers with
or Y primary the lost phase is never re-generated. In Yg -Yg
transformers the core construction has a major impact on the
transformers response to an open phase. Depending on the
core type a Yg -Yg transformer may or may not re-generate the
lost phase, making it a special case when dealing with open
phase conditions in transformers, although it is less commonly
used than Yg - transformers. Figure 16 illustrates a Yg -Yg
transformer with phase A open on the high side.
Transformer LV Voltages- Van, Vbn, Vcn
0.8
0.8
pu
Fig. 14.
pu
0.6
0.6
0.4
0.4
I2/I1- Magnitiude
1
-40
-60
degrees
Table A shows the common phase angle shifts in transformers. Examples of a transformer with primary and secondary
in phase are - (D1-D1) or Y -Y . For other special phase
angle shifts the corresponding II21 can be further determined
with similar exercise. Since phase sequence of the power
system (ABC or ACB) also impacts the vector rotation of
a transformer phase angle shift must be carefully determined
before using the above equations. In any case, with a or
ungrounded Y primary and upon loss of a single phase positive
and negative sequence currents will have equal magnitudes,
i.e. |I2 | = |I1 | both on primary and secondary sides. By
0.5
-80
-100
-120
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
Fig. 15. Simulation Results for a -Yg Transformer After Loss of Phase
B- Load=100KW
1.05
0.036
0.034
pu
0.032
0.95
0.03
1.05
0.06
pu
pu
Fig. 16.
pu
0.04
0.02
0.95
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
Fig. 17. Simulation Results for a Three-leg Yg -Yg Transformer After Loss
of Phase A- Load=60KW
10
1
0.3
0.8
0.6
pu
pu
0.2
0.4
0.1
0.2
0
pu
0.3
0.8
0.4
0.35
0.6
pu
0.2
0.4
pu
pu
0.9
0.1
0.3
0.2
0.8
0
0.25
0.7
0.6
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
0.2
Fig. 19. Simulation Results for a Five-leg Yg -Yg Transformer After Loss
of Phase A- Load=600KW
1.1
0.6
pu
pu
0.9
A. Yg - Transformers
0.4
0.8
0.2
0.7
0.6
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
Fig. 18. Simulation Results for a Three-leg Yg -Yg Transformer After Loss
of Phase A- Load=600KW
Fig. 20.
Figure 20 illustrates a Yg - transformer with winding configured as D1 where it lags the Y side by 30
if system sequence is ABC. With transformers ratio being
V nom.
n = | VAB
| the current ratio is: IIba
= IIcb
= IIac
= n3 .
A
B
C
abnom.
Using mathematical characteristics of fractions we will have:
11
n
3I0
n
Iba + Icb + Iac
=
=
(16)
IA + IB + IC
IG
3
3
Where I0 is the circulating zero sequence current inside
winding and IG is the measured ground current by the
high side neutral current transformer. Therefore I0 can be
calculated in the following equation:
1
n
(17)
I0 = I0Y = nIG
3
3 3
In Figure 20 we have the following set of equations for the
low side currents of the transformer:
Fig. 21.
Ia = Iba Iac
Ib = Icb Iba
Ic = Iac Icb
(18)
Ia = Ica Iab
Ib = Iab Ibc
Ia + Ib + Ic = 0
In the above set of equations Ia , Ib and Ic are available from
the low side phase CTs and I0 is available from equation 17
which comes from the neutral CT of the Yg side. We can now
solve the set of equations in (18) to obtain Iba , Icb and Iac .
By doing so the coil currents of the winding are:
1
(Ia Ib ) + I0
3
1
Icb = (Ib Ic ) + I0
(19)
3
1
Iac = (Ic Ia ) + I0
3
Then high side currents of the transformer are obtained by
V nom.
simply applying the transformer ratio, n = | VAB
|, to the
abnom.
currents which leads us to the following set of equations:
Iba =
3
Iba
n
3
IB =
Icb
(20)
n
3
IC =
Iac
n
Equations 19 and 20 provide primary side currents in a
Yg - transformer in any conditions and loading level as
long as there is no internal fault to the transformer. While
external faults dont violate these equations an internal fault
in the transformer can make the equations invalid. In our
application of open phase condition by using these equations
we can identify an open phase by observing that one of the
phase currents of the high side is zero. Considering the CTs
inaccuracy we may declare an open phase when one of the
currents obtained from equation 20 is much below the other
two phases while currents on the low side of the transformer
appear to be highly balanced.
One other way of connecting three coils in a winding is
illustrated in Figure 21. In this configuration side leads the
Y side by 30 if system sequence is ABC and is called D11.
IA =
Ic = Ibc Ica
(21)
(22)
3
Ica
n
3
Iab
IB =
n
3
IC =
Ibc
n
IA =
(23)
12
Amp.
Fig. 23.
75
70
Transformer Measured HV Currents- IA, IB, IC
Amp.
6
4
nom.
1
ICA = Ic
n 3
1
IBC = Ib
n 3
1
IAB = Ia
n 3
0
Transformer Calculated HV Currents- IA, IB, IC
6
4
2
0
0.05
0.1
Time (Sec)
0.15
(24)
0.2
Fig. 22. Simulation Results for Measured and Calculated Primary Currents
in a Yg - Transformer After Loss of Phase A- Load=60KW
With ICA , IBC and IAB obtained from (24) we can now
calculate the primary side line currents per below:
IA = IAB ICA
IB = IBC IAB
(25)
IC = ICA IBC
The above set of equations are true whether or not there is
zero sequence current on the low side (Yg ), for instance in case
of an unbalance load. This is because any I0 on the low side
will be considered in the high side coils by equation 24,
although I0 will be cancelled out in line currents by (25).
Figure 24 illustrates simulation results for the same system
in Section III-B2 with an unbalanced total load of 100KW.
Note that the unbalanced load creates zero sequence current
on the low side which will be circulating inside the winding
coils. The calculated primary side currents identify that phase
A is open.
13
(26)
Fig. 25.
IA = IAC IBA
IB = IBA ICB
(27)
IC = ICB IAC
160
120
Amp.
Amp.
120
80
80
40
40
4
Amp.
Amp.
0.05
0.1
Time (Sec)
0.15
0.2
0.05
0.1
Time (Sec)
0.15
0.2
Fig. 24. Simulation Results for Measured and Calculated Primary Currents
in a -Yg Transformer After Loss of Phase A- Unbalanced Load=100KW
With similar exercise primary side currents can be calculated from the low side CTs in other less commonly used
transformers such as Y -Yg or - which are not presented
here to stay focused on the commonly used transformers. This
method of calculating the high side current is also true and
applicable in DG cases with direct use of the above equations.
In summary, there are two reliable and effective methods
for identification and protection of open phase conditions in
transformers with or ungrounded Y primary winding: the
method of II21 values of the low side and the method of high
side current calculation. These methods may be combined
as a single algorithm to provide a robust and effective way
to identify an open phase condition in transformers with
ungrounded primary.
C. Yg -Yg Transformers
Although Yg -Yg transformers demonstrate different characteristics in open phase conditions depending on the core type
the principle of calculating the high side currents from the
low side CTs provides a general method for identification
and protection of open phase conditions in Yg -Yg transformers
regardless of the physical construction.
Figure 26 illustrates a Yg -Yg transformer with one phase
open. Positive and negative sequence currents of the transformers primary and secondary are still governed by the
transformer ratio in the new unbalanced system after loss of
V nom.
| we will have :
a phase. If n = | VAB
ab
nom.
I1Sec.
I2
= Sec. = n
(28)
I1P ri.
I2P ri.
However, the above is not necessarily true for the zero sequence currents in three-leg core transformers. Zero sequence
magnetizing current and core losses in three-leg cores are
significantly larger than that of a positive sequence model.
Therefore, the parallel branch in the transformers zero sequence model has more considerable impact, i.e. the current
in the parallel branch may not be totally ignored. Hence I0 of
the two sides will have a ratio different than n. Furthermore,
upon loss of a phase in a three-leg transformer, the primary
side has only two phases with very high I0 (IG = IB + IC ).
However, a large portion of this zero sequence current flows in
the parallel branch for zero sequence magnetization and core
loss and a small portion will be induced in the secondary side.
The secondary side is still a three-phase system with much
smaller I0 , although not necessarily highly balanced depending
on the load size. Therefore each side will have a different I0
that can freely circulate in each grounded side.
Fig. 26.
14
160
Amp.
Amp.
10
140
120
Amp.
Amp.
10
5
0.05
0.1
Time (Sec)
0.15
0.2
10
V. C ONCLUSION
0.05
0.1
Time (Sec)
0.15
0.2
[1] A.E. Fitzgerald, C. Kingsley Jr. and S.D. Umans, Electric Machinery,
4th ed. McGraw-Hill, 1983.
[2] Stephen J. Chapman, Electric Machinery Fundamentals, 4th ed.
McGraw-Hill, 2005.
[3] J. Lewis Blackburn, Thomas J. Domin, Protective Relaying- Principles
and Applications, 3rd ed. CRC Press, 2006.
[4] John Horak and Gerald F. Johnson, A Practical Guide for Detecting
Single-Phasing on a Three-Phase Power System, Western Protective
Relay Conference, 2002.
[5] Jialong Wang and Raluca Lascu, Zero Sequence Circuit of Three-legged
Core Type Transformers. (Basler Electric)
[6] E. L Harder, Sequence Network Connections for Unbalanced Load and
Fault Conditions, The Electric Journal, December 1937.
[7] Paul M. Anderson, Analysis of Faulted Power Systems, IEEE Press,
1995.
Fig. 27. Simulation Results for Measured and Calculated Primary Currents
in a Three-leg core Yg -Yg Transformer After Loss of Phase A- Load=120KW
Amir Norouzi received his B.Sc degree from Tehran Polytechnic in 1996
and his M.Sc from the University of New Brunswick, NB, in 2003 both in
electrical engineering. He has over 13 years of experience in power system
engineering including substations, industrial applications, power generation
and distribution. He joined GE Digital Energy in Markham, Ontario, in 2011
as a technical application engineer, focusing on providing system solutions in
protection and control.
Presented at MIPSYCON 2013
Previously Presented at 2012 WPRC and Georgia Tech 2013.