FaultCal Adv
FaultCal Adv
FaultCal Adv
Lecture 2
Fault Calculations
V I
V pu= I pu =
Vb Ib
MVA
MVA pu = MVA b = 3Vb ph I b ph
MVAb
Z Vb
Z pu = where Z b =
Zb Ib
IbZ
=
Vb
3Vb I b Z MVA SC = 3Vb I f
=
3Vb2 MVAb
=
MVAb Z Z T pu
=
(Vb line )2
Symmetrical Components
V a= Zs Ia + Zm Ib + Zm Ic
V b= Zm Ia + Zs Ib+ Zm Ic
V c= Zm Ia+ Zm Ib+ Zs Ic
[V ] = [Z ][I ]
⎡ ⎤ ⎡V a ⎤ ⎡ ⎤ ⎡I a⎤ ⎡ ⎤ ⎡ Z s Z m Z m⎤
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢V ⎥ = ⎢V b ⎥ ⎢ I ⎥ = ⎢ I b⎥ ⎢Z ⎥ = ⎢ Z m Z s Z m ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ ⎦ ⎢⎣V c ⎥⎦ ⎣ ⎦ ⎢⎣ I c ⎥⎦ ⎣ ⎦ ⎢⎣ Z m Z m Z s ⎥⎦
⎡ ⎤ ⎡1 1 1⎤ ⎡1 1 1⎤
-1
⎡ ⎤
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ 1⎢ ⎥
⎢ ⎥
A = ⎢1 a
2
a ⎥ ⎢ ⎥
A = ⎢ 1 a a
2
⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ 3 ⎢ ⎥
⎣ ⎦ ⎢⎣1 a a 2 ⎥⎦ ⎣ ⎦ ⎢⎣1 a 2 a ⎥⎦
a = 1∠120° 1+ a + a2 = 0
[V ] = [Z ][A][A]-1 [I ]
[A]−1 [V ] = [A]-1 [Z ][A]× [A]-1 [I ]
-1
⎡ ⎤ ⎡ ⎤ ⎡1 1 1⎤ ⎡V a ⎤ ⎡V o ⎤ ⎡ ⎤
⎢ ⎥ ⎢ ⎥ 1⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢ A ⎥ ⎢V ⎥ = 3 ⎢1 a a ⎥ ⎢V b ⎥ = ⎢V 1⎥ = ⎢ V s ⎥
2
⎢ ⎥ ⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ ⎦ ⎣ ⎦ ⎣⎢1 a
2
a ⎥⎦ ⎢⎣V c ⎦⎥ ⎣⎢V 2 ⎦⎥ ⎣⎢ _ ⎦⎥
V0 = 1 (Va + Vb + Vc )
3
(
V1 = 1 Va + aVb + a 2Vc
3
)
(
V2 = 1 Va + a 2Vb + aVc
3
)
-1
⎡ ⎤ ⎡ ⎤ ⎡1 1 1⎤ ⎡ I a ⎤ ⎡ I o ⎤ ⎡ ⎤
⎢ ⎥ ⎢ ⎥ 1⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢ ⎥ ⎢ ⎥ 3
A I = ⎢ 1 a a
2
⎥ ⎢ I b ⎥ = ⎢ I 1⎥ = ⎢ Is ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ ⎦ ⎣ ⎦ ⎣⎢1 a
2
a ⎦⎥ ⎣⎢ I c ⎦⎥ ⎣⎢ I 2 ⎦⎥ ⎣⎢ _ ⎦⎥
I 0 = 1 (I a + I b + I c )
3
I 1 = 1 (I a + aI b + a 2 I c )
3
I 2 = 1 (I a + a 2 I b + aI c )
3
[A]−1 [Z ][A] = [Z ]
⎡1 1 1⎤ ⎡ Z s Z m Z m ⎤ ⎡1 1 1⎤
1⎢ ⎥⎢ ⎥⎢ ⎥
⎢1 a a 2 ⎥ ⎢ Z m Z s Z m ⎥ ⎢1 a 2 a ⎥
3⎢ ⎥⎢ ⎥⎢ ⎥
⎢⎣1 a 2 a ⎥⎦ ⎢⎣ Z m Z m Z s ⎥⎦ ⎢⎣1 a a 2 ⎥⎦
⎡Z s + 2 Z m 0 0⎤ ⎡Z o 0 0⎤
⎢ ⎥ ⎢ ⎥
=⎢ 0 Zs - Zm 0⎥ = ⎢ 0 Z 1 0⎥
⎢ ⎥ ⎢ ⎥
⎣⎢ 0 0 Z s - Z m ⎦⎥ ⎣⎢ 0 0 Z 2 ⎦⎥
⎡V o ⎤ ⎡ Z o 0 0 ⎤ ⎡ I o ⎤
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢V 1⎥ = ⎢ 0 Z 1 0 ⎥ ⎢ I 1⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢⎣V 2 ⎥⎦ ⎢⎣ 0 0 Z 2 ⎥⎦ ⎢⎣ I 2 ⎥⎦
Hong Kong Polytechnic University 05/09/2005 11
[ I ] = [ A ] [ Is ]
⎡ I a ⎤ ⎡1 1 1⎤ ⎡ I o ⎤
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢ I b⎥ ⎢ a
= 1 2
a ⎥ ⎢ I 1⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎣⎢ I c ⎦⎥ ⎣⎢1 a a ⎦⎥ ⎣⎢ I 2 ⎦⎥
2
⎡V a ⎤ ⎡1 1 1⎤ ⎡V o ⎤
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢V b ⎥ = ⎢1 a 2 a ⎥ ⎢V 1⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎣⎢V c ⎦⎥ ⎣⎢1 a a ⎦⎥ ⎣⎢V 2 ⎦⎥
2
Ia
Ib
Ic
Va = 0 ∴, V1 + V2 + V0 = 0
I a = I f = I1 + I 2 + I 0
Ib = Ic = 0
1
I1 = I 2 = I 0 = I f
3
Z1T
+ve G E1
Z2T
-ve
Z0T
0
1 1
I1 = I 2 = I 0 = = = − j 0.8
Z 1T + Z 2T + Z 0T j1.25
I F = 3I 0 = − j 2.4
V1 = E − I 1 Z 1S = 1 − (− j 0.8)( j 0.1) = 0.92
V2 = − I 2 Z 2 S = −(− j 0.8)( j 0.15) = −0.12
V0 = − I 0 Z 0 S = −(− j 0.8)( j 0.2 ) = −0.16
V A = V1 + V2 + V0 = 0.64
V B = a 2V1 + aV2 + V0 = 0.92∠240° − 0.12∠120° − 0.16
= −0.56 − j 0.9007 = 1.06∠ − 122°
V B = aV1 + a 2V2 + V0 = 0.92∠120° − 0.12∠240° − 0.16
= −0.56 + j 0.9007 = 1.06∠122°
Hong Kong Polytechnic University 05/09/2005 17
Ib
Ic
Rf
Va = I f Zf ∴, V1 + V2 + V0 = I f R f
I a = I f = I1 + I 2 + I 0
Ib = Ic = 0
1
I1 = I 2 = I 0 = I f
3
V1 + V2 + V0 = 3I 0 R f
Z1T
+ve G E1
-ve Z2T
Z0T
0
3Rf
Example 2 – Calculate
The value of the neutral earthing resistance in
pu to limit the magnitude of the fault current
of a L-E fault to 1.2 pu
Z1=j0.1 pu
Z2=j0.15 pu
Z0=j0.2 pu
R
G •
Z1=Z2=j0.2 pu
E=1.0 pu
Z0=j0.4 pu
3 3
I F = 3I 0 = =
Z 1T + Z 2T + Z 0T + 3R N 3R N + j1.25
3
IF = = 1.2
3R N + j1.25
32
= 1.2 2
(3R N ) 2
+ 1.25 2
R N = 0.72
Ia
Ib
Ic
Ia = 0 ∴, I b = − I c = I f
I1 + I 2 = 0
Vb = Vc = V f
V1 = V2
Z1T
+ve G E1
-ve Z2T
Z0T
0
1 1
I1 = I 2 = = = − j1.538
Z1T + Z 2T j 0.65
I F = I B = a 2 I 1 + aI 2 = a 2 − a I 1( )
= j 3I 1 = 2.66∠180°
Ib
Ic
Rf
Ia = 0 I1 + I 2 + I 0 = 0
I f = Ib + Ic I b + I c = 3I 0
Vb = Vc = I f R f
V1 =
1
3
( 1
3
)
V a + aV b + a 2 V c = (V a − V b )
1
3
( 1
3
)
V 2 = V a + a 2 V b + aV c = (V a − V b )
V 0 = (V a + V b + V c ) = (V a + 2V b )
1 1
3 3
= V1 + V b = V1 + 3 I 0 R f
V1 = V 2 = V 0 − 3 I 0 R f
Z1T
+ve G E1
-ve Z2T
Z0T
0
3Rf
Ib
Ic
Rf
Z1T
+ve G E1
-ve Z2T
Z0T
0
3Rf
V A1 = C1 I 1 Z LA1 + V1
V A 2 = C 2 I 2 Z LA 2 + V2 = C1 I 1 Z LA1 + V2
V A0 = C 0 I 0 Z LA0 + V0
V1 + V2 + V0 = 3I 0 R f
V A = 2C1 I 1 Z LA1 + C 0 I 0 Z LA0 + 3I 0 R f
Earthing Methods
Solidly earthed
Resistance earthed
Reactance earthed
Resonant or Peterson Coil earthed
Insulated earthed
Earthing transformer
Solidly earthed
The zero sequence impedance is a minimum. The
earth fault current is highest. However, overvoltage
in unfaulted phases is a minimum. Earth fault
protection is simple as the fault current is usually
high. Overvoltage during earth faults is usually less
than 0.8 times phase to phase voltage.
Resistance earthed
Earthing the system by resistor reduces both the fault
current and transient overvoltages. However, high
value of earthing resistance will produce overvoltage
in unfaulted phases during earth faults approaching
to phase to phase voltage.
Hong Kong Polytechnic University 05/09/2005 48
Reactance earthed
Reactors are usually smaller and less expensive than
resistors, it can reduce fault current but this is
associated with a increase in transient overvoltages
during earth faults. Transient overvoltages are a
maximum when the value of earthing reactance is
approximately one third of the value required for
resonant earthing. In order to reduce the transient
overvoltages to an admissible level the earthing
reactance has to be reduced so that the earth fault
current approaches that for a solidly earthed system.
Overvoltages during earth faults will be between 0.8
times and full phase to phase voltage.
Hong Kong Polytechnic University 05/09/2005 49
Insulated earthed
The earth fault current is capacitive and if
small may be self extinguished. Automatic
segregation of faulty zones is extremely
difficult. The system can be run with an
earth fault for long periods. Overvoltages
during earth faults may be greater than phase
to phase voltage. Arcing earth faults are very
likely and these can result in high transient
overvoltages.
Earthing transformer
In some cases where the neutral of the power
transformer is not available, the system is earthed
via a earthing transformer of zig-zag type. A
minimum impedance is offered to the flow of zero
sequence currents. Under normal operating
conditions the currents flowing through the
windings are the magnetising currents of the
earthing transformer. The earthing transformer is
designed to carry the maximum fault current for 30
seconds.
Effectively
earthed system
Non-effectively
earthed system
A-E Fault
3K 2 K 0
Va = Ea
(K 2 + K 0 + K 2 K 0 )
3K 0
= Ea if K 2 = 1
(1 + 2K 0 )
B-C-E Fault
3V AN
VNG = − I F Z E = − ZE
Z1 (1 + K 2 + K 0 ) + 3Z E
ZE
=− V AN if K 2 = K 0 = 1
Z1 + Z E
3
=2
ZT1 + ZT 2 + Z T 0 + 3RN
3
=2
3RN + j 0.787
[(3R N ) 2 + 0.787 2 = 2.25 ]
RN = 0.4256 = 49.4Ω
1
I1 = I 2 = I 0 =
3 × 0.4256 + j 0.787
= 0.5675 − j 0.35 = 0.667∠ − 31.6°
V B / B1 = E − I 1 ( Z S 1 + Z Tx1 )
= 1 − (0.5675 − j 0.35)( j 0.221) = 0.923 − j 0.125
V B / B 2 = − I 2 ( Z S 2 + Z Tx 2 )
= −(0.5675 − j 0.35)( j 0.264) = −0.0923 − j 0.1498
VB / B0 = − I 0 ( Z Tx 0 + 3 R N )
= −(0.5675 − j 0.35)(1.277 + j 0.125) = −0.768 + j 0.376
VB / B A = VB / B1 + VB / B 2 + VB / B 0
= 0.06 + j 0.1 = 9.1kV
VB / B B = a 2VB / B1 + aVB / B 2 + VB / B 0
= (−0.5 − j 0.866)(0.923 − j 0.125) +
(−0.5 + j 0.866)(−0.0923 − j 0.1498) + (−0.768 + j 0.376)
= −1.162 − j 0.366 = 1.22∠ − 162.5° = 93kV
VB / B C = aVB / B1 + a 2VB / B 2 + VB / B 0
= (−0.5 + j 0.866)(0.923 − j 0.125) +
(−0.5 − j 0.866)(−0.0923 − j 0.1498) + (−0.768 + j 0.376)
= −1.21 + j1.393 = 1.84∠131° = 140kV
I L 0 n C = I H 0 (n S + n C )
With the tertiary winding open circuited
nC
I H 0= I L 0
n S + nC
1
= I L0
1+ N
nS
where N =
nC
EH n + nC
= S =1+ N
EL nC
⎡ 1 ⎤
I L 0 − I H 0 = I L 0 ⎢1 −
⎣ 1 + N ⎥⎦
N
= I L0
1+ N
Hong Kong Polytechnic University 05/09/2005 71
I T 0 (Z X + Z Z ) = I L 0 (Z L + Z T + 3 Z n )
Z X + Z Z = Z L + Z T + 3Z n by using p.u. values
N
As Z X = Z L + 3Z n
1+ N
1
∴, Z Z = Z T + 3Z n
1+ N
Hong Kong Polytechnic University 05/09/2005 73
E xa = E ya + I a Z a and E x0=
1
3
[
(E ya + I a Z a ) + E yb + E yc ]
E xb = E yb
[
= (E ya + E yb + E yc ) + I a Z a
1
3
]
E xc = E yc
1
E x0= E y0+ I a Z a
3
1
and similarly E x1 = E y1 + I a Z a
3
1
E x2= E y2+ I a Z a
3
Hong Kong Polytechnic University 05/09/2005 74
IaZa = Za (I1 + I2 + I0 )
1 1
3 3
Ib I b’
Ic Ic’
G
Hong Kong Polytechnic University 05/09/2005 76
Strong source
Weak Source
Further Study
Paul M. Anderson, “Analysis of Faulted
Power Systems”, IEEE Press, 1995.
Self Evaluation
Do you understand the physical meaning of positive
and negative sequence components?
What is the relationship between mutual coupling
impedance and sequence impedance in the power
system?
How do you transform between the actual system
voltages and currents into their corresponding sequence
components?
Do you understand why the actual system network can
be transformed into sequence networks?
How the sequence networks are connected under
various unsymmetrical faults?
Hong Kong Polytechnic University 05/09/2005 85
Self Evaluation
Do you know how to form the zero sequence network
of a power transformer? Now do you understand how to
form the zero sequence network of the power system?
Do you know how to solve tutorial questions from 1 to
10?
How do you connect the sequence networks for an open
circuited fault?
How to represent the simultaneous fault by sequence
networks?
Are you able to solve the open circuit and simultaneous
fault problems from Q11 to Q15?