EE251 Notes0 Power System Models
EE251 Notes0 Power System Models
EE 251
Fault Studies
Review Notes
Power System Models
Prof. Rowaldo del Mundo The following notes were prepared by the
Electrical Engineering Faculty of the Electrical &
Associate Professor Electronics Engineering Institute of the University
of the Philippines
1
Transmission & Distribution Lines
Sequence Impedance
2
Geometric Mean Distance
Typical three-phase line configurations
D12 D23
D31
D12
D12 D23
D31
Dm 3 D12D23D31
Ds Dscd
Ds 1.09 4 Dscd3
Note: Dsc=GMR of a single conductor
3
Aluminum Conductors
Aluminum is preferred over copper as a material
for transmission lines due to:
• lower cost
• lighter weight
• larger diameter for the same resistance*
Stranded Conductors
Alternate layers of wire of a Steel
stranded conductor are
spiraled in opposite directions
to prevent unwinding and
make the outer radius of one
layer coincide with the inner Aluminum
radius of the next.
The number of strands depends on the number
of layers and on whether all the strands are of
the same diameter. The total number of strands
of uniform diameter in a concentrically stranded
cable is 7, 19, 37, 61, 91, etc.
4
Example: Each conductor of the 345 kV bundled
three-phase line shown is 1272 MCM ACSR 54/19
Pheasant. Assuming complete transposition and a
total length of 80 miles, determine the positive-
sequence impedance of the line. Assume d=45 cm
and D=8 meters.
d
We get
Dm 3 8(8)(16) 10.08 m
Ds 0.0142(0.45) 0.08 m
The total resistance per phase
R1 0.0821(80)(0.5) 3.28
The total reactance per phase
10.08
X1 2.02 x 10-3(60)(80) ln 46.9
0.08
The positive (or negative) sequence impedance of
the line
Z1 Z2 3.28 j46.9 / phase
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 10
5
Transmission Line Inductance
The total self-inductance of a cylindrical conductor
consists of an internal inductance and an external
inductance.
Internal Inductance:
c
L int Henry/meter
8
Note: c is the permeability of the conductor.For
a non-ferrous conductor such as copper or
aluminum,
c 4 x 10-7 Henry/meter
External Inductance:
The external inductance may be approximated by
m 2s
L ext (ln 1) Henry/meter
2 r
where
m = permeability of the medium surrounding
the conductor
r = radius of the conductor
6
Self Inductance:
c m 2s
L (ln 1) Henry/meter
8 2 r
For a non-ferrous conductor in air, we get
2s
L 2 x 10-7 (ln 1) Henry/meter
Ds
Mutual Inductance:
The mutual inductance between two conductors is
2s
M 2 x 10-7 (ln 1) Henry/meter
Dm
where Dm is the geometric mean distance
between the conductors
Note: For any unit of length,
2s
L k (ln 1) Henry/unit length
Ds
2s
M k (ln 1) Henry/unit length
Dm
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 14
7
Carson’s Line
Carson examined a single overhead conductor
whose remote end is connected to earth.
a
zaa a’
+
Local Remote
Earth
Va Ia zad Earth
-
Fictitious
REF Vd 0 Return
d d’ Conductor
Id zdd
The current returns through a fictitious earth
conductor whose GMR is assumed to be 1 foot (or
1 meter and is located a distance Dad from the
overhead conductor.
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 15
8
Primitive Impedances:
2s
zaa ra jL a ra jk(ln 1)
Dsa
2s
zdd rd jk(ln 1)
Dsd
2s
zad jM jk(ln 1)
Dad
ra, rd = resistances of overhead conductor and
fictitious ground wire, respectively
Dsa, Dsd = GMRs of overhead conductor and
fictitious ground wire, respectively
Note: Dsd is arbitrarily set to 1 unit length
Earth Resistance:
Carson derived an empirical formula for the earth
resistance.
Note : At 60 Hz,
rd 0.09528 /mile
9
Equivalent Impedance:
Substitute the primitive impedances into
zaa zaa zdd 2zad
We get 2
D
zaa (ra rd ) jk ln ad
DsaDsd
Note: 1. Dsd=1 unit length
2
D
2. Define De ad
Dsd
We get
De
zaa (ra rd ) jk ln /unit length
Dsa
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 19
10
Three-Phase Line Impedances
Extend Carson’s
analysis to a three-phase line
Ia zaa
a a’
Ib zbb zab zca b’
b
+
+ Ic zcc zbc
Va c c’
All wires
- Vb + zad grounded
- Vc zbd here
- zcd
REF
Vd 0
d d’
Id zdd
1 unit length
11
The voltage equation for the fictitious ground
conductor
is
Vd Vd' zadIa zbdIb zcdIc zddId
which can be re-written as
Vd' zadIa zbdIb zcdIc zdd(Ia Ib Ic )
Combining terms, we get
Vd' (zad zdd)Ia (zbd zdd)Ib (zcd zdd)Ic
Following a similar procedure, the voltage
equation
for phase a can be shown to be
Va Va' (zaa zad)Ia (zab zad)Ib
(zac zad)Ic
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 23
12
We can also subtract the voltage equation of the
ground conductor from the equations of phases b
and c. The resulting matrix equation is
Va zaa zab zac I
a
Vb = zab zbb zbc Ib V/unit length
Vc zac zbc zcc Ic
Self Impedances:
zaa zaa 2zad zdd /unit length
zbb zbb 2zbd zdd /unit length
zcc zcc 2zcd zdd /unit length
Mutual Impedances:
zab zab zad zbd zdd /unit length
zbc zbc zbd zcd zdd /unit length
zac zac zad zcd zdd /unit length
Primitive Impedances:
2s
zxx rx jk(ln 1) /unit length
Dsx x=a,b,c,d
2s
zxy jk(ln 1) /unit length
Dxy xy=ab,bc,ca,ad,bd,cd
13
Assumptions:
1. Identical phase conductors
Dsc=0.01668 ft @ 60 Hz a b c
De
zaa zbb zcc (ra rd ) jk ln
Ds
(0.278 0.095) j0.121 ln 0.01668
2790
14
zab zbc 0.095 j0.121 ln 2790
10
0.095 j0.683 /mile
Zab 3.81 j27.33
Line Transposition
Line transposition is used to make the mutual
impedances identical.
Ia
Phase c
Pos.1
Ib
Phase a
Pos.2
Ic
Phase b
Pos.3
s1 s2 s3
Section 1 Section 2 Section 3
15
Voltage equation for Section 1
Va Z111 Z121 Z131 I
a
Vb = Z211 Z221 Z231 Ib volts
Vc Z311 Z321 Z331 Ic
16
Va (Z111 Z222 Z333 )Ia
(Z121 Z232 Z313 )Ib
(Z131 Z212 Z323 )Ic
By a similar procedure, we get for phases b and c
Vb (Z211 Z322 Z133 )Ia
(Z221 Z332 Z113 )Ib
(Z231 Z312 Z123 )Ic
Vc (Z311 Z122 Z233 )Ia
(Z321 Z132 Z213 )Ib
(Z331 Z112 Z223 )Ic
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 33
17
Define f1, f2 and f3 as as the ratios of s1, s2 and s3
to the total length s, respectively. We get
Va (f1Z11 f2Z22 f3Z33)Ia
(f1Z12 f2Z23 f3Z31)Ib
(f1Z13 f2Z21 f3Z32 )Ic
Vb (f1Z21 f2Z32 f3Z13)Ia
(f1Z22 f2Z33 f3Z11)Ib
(f1Z23 f2Z31 f3Z12 )Ic
Vc (f1Z31 f2Z12 f3Z23)Ia
(f1Z32 f2Z13 f3Z21)Ib
(f1Z33 f2Z11 f3Z22 )Ic
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 35
Define:
18
It can be shown that
De
Zs (ra rd )s jks ln
Ds
D D D
Zk1 rds jks f1ln e f2 ln e f3 ln e
D12 D23 D31
D D D
Zk2 rds jks f1ln e f2 ln e f3 ln e
D31 D12 D23
D D D
Zk3 rds jks f1ln e f2 ln e f3 ln e
D23 D31 D12
19
Solving for the mutual impedances, we get
3.81 j25.65
Similarly, we get
20
Example: For the same line assume the
following transposition cycle:
s1=16 miles s2=24 miles s3=0
21
Example: For the same line assume a complete
transposition cycle.
10’ 10’
The GMD is
a b c
Dm 3 10(10)(20) 12.6 feet
We get the average of the mutual impedances.
Zm 3.81 j26.21
The impedance of the transposed line is
Phase- to Sequence-Impedances
Consider a transmission line that is described by
the following voltage equation:
Va Zaa Zab Zac I
a
Vb = Zab Zbb Zbc Ib volts
Vc Zac Zbc Zcc Ic
or
Vabc ZabcIabc
From symmetrical components, we have
Vabc AV012 and Iabc A I012
22
Substitution gives
AV012 ZabcA I012
or
V012 A 1ZabcA I012
which implies that
Z012 A 1Zabc A
Performing the multiplication, we get
Z0 Zs0 2Zm0 Zs2 Zm2 Zs1 Zm1
Z1 = Zs1 Zm1 Zs0 Zm0 Zs2 2Zm2
Z2 Zs2 Zm2 Zs1 2Zm1 Zs0 Zm0
Note: Z012 is not symmetric.
Zs0 1
3
(Zaa Zbb Zcc )
Zs1 1
3
(Zaa aZbb a2Zcc )
Z s2 1
3
(Zaa a2Zbb aZcc )
Zm0 1
3
(Zab Zbc Zca)
Zm1 1
3
(a2Zab Zbc aZca )
Zm2 1
3
(aZab Zbc a2Zca )
23
If the line is completely transposed,
Zs0 Zs Zm0 Zm
Zs1 Zs2 0 Zm1 Zm2 0
The sequence impedance matrix reduces to
Z0 Zs 2Zm 0 0
Z1 = 0 Zs Zm 0
Z2 0 0 Zs Zm
24
Example: For the same line and assuming a
complete transposition cycle, find the sequence
impedances of the line.
10’ 10’
25
Incomplete Transposition
The resulting circuit unbalance due to incomplete
transposition is measured in terms of the zero-
and negative-sequence line currents that arise
when positive-sequence voltages are applied at
the sending end.
For a partially-transposed transmission line,
Va zs zk1 zk 2 I
a
Vb = zk1 zs zk 3 Ib V/unit length
Vc zk 2 zk 3 zs Ic
or
Vabc zabcIabc
where
De
zs (ra rd ) jk ln
Ds
D D D
zk1 rd jk f1ln e f2 ln e f3 ln e
D12 D23 D31
D D D
zk2 rd jk f1ln e f2 ln e f3 ln e
D31 D12 D23
D D D
zk3 rd jk f1ln e f2 ln e f3 ln e
D23 D31 D12
26
Transform the phase impedances into sequence
impedances.
z012 A 1zabc A
which gives
z0 zs 2zm0 zm2 zm1
z1 = zm1 zs zm0 2zm2
z2 zm2 2zm1 zs zm0
where
zm0 1
3
(zk1 zk2 zk3 )
zm1 1
3
(a2zk1 azk 2 zk 3 )
zm2 1
3
(azk1 a2zk 2 zk 3 )
27
Assume that balanced voltages are applied at the
sending end;
that is
Va0 0 and Va2 0
We can re-write the voltage equation
Va0 z00 z01 z02 Ia0
Va1 = z10 z11 z12 I
a1
Va2 z20 z21 z22 Ia2
into
0 z00Ia0 z01Ia1 z02Ia2
Va1 z10Ia0 z11Ia1 z12Ia2
0 z20Ia0 z21Ia1 z22Ia2
The unbalanced
factor are definedas:
I I
m0 a0 and m2 a2
Ia1 Ia1
It can be shown that
z02z21 z01z22
m0
z00z22 z20z02
z20z01 z00z21
m2
z00z22 z20z02
Note: For the untransposed 40-mile line,
m0 1.05472.68o %
m2 6.624139.14o %
28
Lines with Overhead Ground Wire
Ia zaa
a
a’
Ib zbb zab zca
b
b’
+ Ic zcc zbc
c
c’
Va + Iw
zw w zad
- Vb + w w’
-
Vc +
zbd All wires
- Vw zcd grounded
- here
REF
Vd 0
d’
d
Id zdd
1 unit length
29
The last equation may be written as
Vd' zdaIa zdbIb zdcIc zdwIw
zdd(Ia Ib Ic Iw )
or
Vd' (zda zdd)Ia (zdb zdd)Ib
(zdc zdd)Ic (zdw zdd)Iw
Following the same procedure, the equation for
phase a can be expressed as
Va Va' (zaa zad)Ia (zab zad)Ib
(zac zad)Ic (zaw zad)Iw
30
We get
Va zaa zab zac zaw I
a
Vb zba zbb zbc zbw Ib
=
Vc zca zcb zcc zcw I
c
Vw zw a zw b zw c zw w Iw
where Vw 0 and
De
zxx (rx rd ) jk ln x=a,b,c,w
Dsx
De
z xy rd jk ln xy=ab,ac,aw,
Dxy bc,bw,cw
Kron’s Reduction
Consider the matrix equation
V1 Z1 Z2 I1
=
0 Z3 Z4 I2
31
From eq. 2, we get
1
I2 Z4 Z3I1
which is substituted in eq. 1 to get
1
V1 Z1I1 Z2(Z4 Z3I1)
or
1
V1 (Z1 Z2Z4 Z3 )I1
Go back to the voltage equation of the line.
Va zaa zab zac zaw I
a
Vb zba zbb zbc zbw Ib
=
Vc zca zcb zcc zcw I
c
0 zw a zw b zw c zw w Iw
We wish to eliminate Iw , where
zaa zab zac zaw
T
Z1 = zba zbb zbc Z2 = zbw = Z3
zca zcb zcc zcw
Z4 zw w
We get
zawzw a z z zawzw c
zaa zab aw w b zac
zw w zw w zw w
z z z z z z
zabc = zba bw w a zbb bw w b
zw w zw w
zbc bw w c
zw w
z z z z z z
zca cw w a zcb cw w b zcc cw w c
zw w zw w zw w
32
Example: Find the equivalent impedance of the
69-kV line shown. The phase conductors are the
same as in the previous examples. The overhead
ground wires have the following characteristics:
rw=4.0 /mile @ 25oC w
Dsw=0.001 ft @ 60 Hz
For the ground wire, we get 15’
De
zw w (rw rd ) jk ln 10’ 10’
Dsw a b c
(4.0 0.095) j0.121 ln 2790
0.001
4.095 j1.8 /mile
Zw w 163.8 j72
De
zaw zcw rd jk ln
Daw
0.095 j0.121 ln 18.03
2790
/mile
Zaw Zcw 3.81 j24.47
zbw 0.095 j0.121 ln 2790
15
/mile
Zbw 3.81 j25.36
From a previous example, we got
33
Using the ground wire impedances, we also get
3.81 j24.47
Z2 = 3.81 j25.36 = Z3T
3.81 j24.47 Z4 163.8 j72
Current Division
How much current is flowing in the overhead
ground wire? Under zero-sequence condition,
Ia Ib Ic Ia0
The primitive voltage equation reduces to
Va Va' zaa zab zac zaw zad Ia0
Vb Vb' zba zbb zbc zbw zbd Ia0
V/unit
Vc Vc' = zca zcb zcc zcw zcd Ia0
length
0 Vw' zw a zw b zw c zw w zw d Iw
0 Vd' zda zdb zdc zdw zdd Id
34
From KCL, we get
Id Iw 3Ia0
The last equation may be written as
Vd' (zda zdb zdc )Ia0 zdwIw zddId
Similarly, for the overhead ground wire
Vw' (zw a zw b zw c)Ia0 zw wIw zw dId
Subtracting the two equations, we get
0 [(zw a zw b zw c) (zda zdb zdc )]Ia0
(zw w zdw)Iw (zw d zdd)Id
Since Id (3Ia0 Iw ) , we get
0 [(zw a zw b zw c) (zda zdb zdc )]Ia0
(zw w zdw)Iw (zw d zdd)(3Ia0 Iw )
Combine
common terms and get the ratio of Iw to
3Ia0 . We get
Iw (zw a zw b zw c) (zda zdb zdc ) 3(zdd zw d)
3Ia0 3(zw w zw d zdw zdd)
Recall that
zwx zwx zwd zdx zdd x=a,b,c,w
35
Substitution gives
Iw zw a zw b zw c
(a)
3Ia0 3zw w
From Id Iw 3Ia0 , we get
Id I
1 w
3Ia0 3Ia0
Substitute eq. (a) to get
Id z zw b zw c
(b) 1 wa
3Ia0 3zw w
36
Line Transposition
Consider a transmission line with one overhead
ground wire. Assume that the line is transposed.
Ground wire
Ia
Phase c
Pos.
Ib
Phase a
Pos.
Ic
Phase b
Pos.
s1 s2 s3
Section 1 Section 2 Section 3
The impedance matrix for each section can be
found using Kron’s reduction.
37
Similarly, for Sections 3, we get
z wzw z wzw z wzw
z z z
zw w zw w zw w
z wzw z z z wzw
z f 3 = z z w w z
zw w zw w zw w
zwzw z z zwzw
z z w w z
zw w zw w zw w
It can be shown that
Va I
a
Vb s(f1Zf1 f2Zf 2 f3Zf 3 ) Ib
Vc Ic
z s zm zm
zabc = zm zs zm /unit length
zm zm z s
zs 13 (z z z )
z z z z z z
13 w w w w w w
zw w zw w zw w
38
The off-diagonal elements are
zm 13 (z z z )
z z z z z z
13 w w w w w w
zw w zw w zw w
The sequence impedance matrix will be diagonal.
z00 0 0
z012 = 0 z11 0 /unit length
0 0 z22
where
z00 zs 2zm and z11 z22 zs zm
39
The last two terms can be shown to be equal to
(k)2 M
1
3
zww
where
De 2 D D
M (ln ) (ln e )2 (ln e )2
Dw Dw Dw
D D D D
(ln e )(ln e ) (ln e )(ln e )
D w Dw Dw Dw
De D
(ln )(ln e )
Dw D w
Likewise,
De
zw w (rw rd ) jk ln
Dw w
For most transmission line configurations,
Dm
z11 ra jk ln /unit length
Ds
Note: This is the formula for the positive sequence
impedance of a transposed line without ground
wires.
40
The zero-sequence impedance can be shown to be
2
z0(m)
z00 z0(a) /unit length
z0(g)
where
3
De
z0(a) (ra 3rd ) jk ln 2
DsDm
De
z0(m) 3rd j3k ln
3 DwDwDw
z0(g) 3zww
41
Zm 13 [2(6.48 j25.0) 6.38 j21.7]
6.45 j23.9
17.57 j56.06 6.45 j23.9 6.45 j23.9
Zabc= 6.45 j23.9 17.57 j56.06 6.45 j23.9
6.45 j23.9 6.45 j23.9 17.57 j56.06
For comparison, Zabc below is for a completely
transposed line without the overhead ground wire.
Z00 Zs 2Zm
30.47 j103.86
Z11 Z22 Zs Zm
11.12 j32.17
Note: When the effect of the ground wire was not
taken into account, we got
Z0 22.55 j110.8
Z1 11.12 j32.17
42
Steel Ground Wires
The overhead Grade E.B.B. 7 Strands
3/8 inch Diameter
ground wires are 10-3 15
Resistance (/mile)
10-7 11
some alloy of GMR
steel. This poses a 10-9 9
problem since the 10-11 7
resistance and 5
10-13
GMR of these Resistance
conductors vary 10-15 3
with current. 10-17 1
10-19
0 10 20 30 40 50 60
Amperes per Cable
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 85
43
The zero-sequence impedance of the ground wire
De
z0(g) 3zw w 3(rw rd ) j3k ln
Dw w
The resistance and GMR of the ground wire is a
function of the current. For 1A, 30A and 60A, we
get for E.B.B.
rw= 3.5 /mile at 1A
= 7.8 /mile at 30A
= 6.0 /mile at 60A
Dsw= 10-3 ft at 1A
= 5 x 10-12 ft at 30A
= 1.5 x 10-10 ft at 60A
44
The zero-sequence mutual impedance is
De
z0(m) 3rd j3k ln
Dag
2790
3(0.095) j3(0.121) ln
16.96
0.29 j1.86 /mile
Z0(m) 11.43 j74.28
Substitution gives
Z00 31.0 j102.65 at 1A
Z00 18.86 j107.06 at 30A
Z00 27.21 j105.89 at 60A
Note: For comparison, we got the following zero-
sequence impedances from previous examples:
No overhead ground wire:
45
Parallel Multicircuit Lines
Consider the double-circuit transmission line
shown. One circuit is designated abc and the
other a’b’c’. Each circuit may be untransposed or
partially transposed.
m Iabc n
Zabc
Ia'b'c'
Za'b'c'
46
Using matrix notation, we get
Vabc Zaa Zaa' I
= abc
Va'b'c' Za'a Za'a' Ia'b'c'
47
Performing the transformation, we get
I A 1Y A A 1Yaa'A V012
012 = 1 aa
I0'1'2' A Ya'aA A 1Ya'a'A V0'1'2'
Assume that the voltages at both ends of the line
contain positive sequence components only. We
get
0 0
Va1 1
V012 0 0
= = Va1
V0'1'2' 0 0
Va'1 1
0 0
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 95
Note that Va1 Va'1 since the lines are bused at
both
ends. Substitution gives
Ia0 Y00 Y01 Y02 Y00' Y01' Y02' 0
I Y10 Y11 Y12 Y10' Y11' Y12' 1
a1
I Y Y21 Y22 Y20' Y21' Y22' 0
a2 = 20 Va1
Ia'0 Y0'0 Y0'1 Y0'2 Y0'0' Y0'1' Y0'2' 0
I Y1'0 Y1'1 Y1'2 Y1'0' Y1'1' Y1'2' 1
a'1
Ia'2 Y2'0 Y2'1 Y2'2 Y2'0' Y2'1' Z2'2' 0
Solving for the sequence currents, we get
Ia0 (Y01 Y01' )(Va1 )
Ia2 (Y21 Y21' )(Va1 )
Ia1 (Y11 Y11' )(Va1 )
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 96
48
Ia'0 (Y0'1 Y0'1' )(Va1 )
Ia'2 (Y2'1 Y2'1' )(Va1 )
Ia'1 (Y1'1 Y1'1' )(Va1 )
Define the net through unbalance factors as
Ia0 Ia'0
m0t per unit
Ia1 Ia'1
I I
m2t a2 a'2 per unit
Ia1 Ia'1
Define the net circulating current unbalance factors
as Ia0 Ia'0
m0c per unit
Ia1 Ia'1
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 97
Ia2 Ia'2
m2c per unit
Ia1 Ia'1
Solving for the unbalance factors, we get
m0t (Y01 Y01' Y0'1 Y0'1' ) / Y1
m2t (Y21 Y21' Y2'1 Y2'1' ) / Y1
m0c (Y01 Y01' Y0'1 Y0'1' ) / Y1
m2c (Y21 Y21' Y2'1 Y2'1' ) / Y1
where
Y1 Y11 Y11' Y1'1 Y1'1'
49
Electric Field of a Long Conductor
Consider a long, isolated cylindrical conductor with
a uniform positive charge density q and located in
a homogeneous medium.
The flux density at a distance x from the center is
q q
D âr Coul/m2 +
2x x
âr
The electric field intensity is
D q
E âr Volt/m
2x
Sequence Capacitance
50
Voltage Due to Charged Conductor
Consider two points P1 and P2 which are located at
distances D1 and D2 from the center of the
conductor.
D1 P1
The voltage drop from
P1 to P2 is P2
D2 D2
v12 D1
E dx
q
D2
D1 2x dx
q D
ln 2 Volts
2 D1
51
Due to charge qb, we also get the voltage drop vba.
qb D
vba ln
2 rb
or q D q r
vab b ln b ln b
2 rb 2 D
Applying superposition, we get the total voltage
drop from charge qa to charge qb.
qa D q r
vab ln b ln b
2 ra 2 D
Since qa+qb=0, we get
qa D2
v ab ln Volts
2 rarb
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 103
2
Can Cbn 2Cab Farad/meter
D
ln
r
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 104
52
Three-Phase Line
Consider a three-phase trans- qb
mission line with equilateral
spacing. Assume identical
phase conductors and neglect D D
the effect of ground. Assume
qa qb qc 0 qa D qc
1 D r
v ab v ac 2q ln (q q ) ln
D
a b c
2 r
Since qa=-(qb+qc), we get
3qa D
vab vac ln Volts
2 r
53
Phasor Diagram:
Vca Vcn Vab
Vab Van( 330o )
Vac Van( 3 30o ) Van
Vab Vac 3Van
Vbn Vac
We get
qa D
van ln
2 r Vbc
The capacitance to neutral is
qa 2 Farad/meter,
Can Cbn Ccn
van D to neutral
ln
r
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 107
qa Phase c
Pos.1
qb Phase a
Pos.2
qc Phase b
Pos.3
1
3
s 1
3
s 1
3
s
Section 1 Section 2 Section 3
54
For section 1, we get b
D12 D23
a D31 c
1 D12 r D
vab1 qa ln qb ln qc ln 23
2 r D12 D31
1 D31 D r
vac1 qa ln qb ln 23 qc ln
2 r D12 D31
D12 D23
c D31 b
1 D23 r D31
vab2 q
a ln qb ln qc ln
2 r D23 D12
1 D12 D31 r
vac2 q
a ln qb ln qc ln
2 r D23 D12
55
For section 3, we get c
D12 D23
b D31 a
1 D31 r D
vab3 qa ln qb ln qc ln 12
2 r D31 D23
1 D23 D r
vac3 qa ln qb ln 12 qc ln
2 r D31 D23
1 Dm r
vab q
a ln qb ln
2 r Dm
56
Similarly, we get
1 Dm r
vac q
a ln qc ln
2 r Dm
Taking summation, we get
1 Dm r
vab vac 2qa ln (qb qc ) ln
2 r Dm
3 D
3van qa ln m Volts
2 r
or
1 D
van qa ln m Volts
2 r
The capacitance of phase a to neutral is
qa 2 Farad/meter,
Can
v an D to neutral
ln m
r
Since the line is completely transposed,
57
Capacitive Reactance
1
xc
2fC
2.862 D -meter,
x 109 ln m to neutral
f r
1.779 D -mile, to
x 106 ln m neutral
f r
Note: To get the total capacitive reactance, divide
xc by the total length of the line.
Equivalent Circuit
A lumped-parameter representation is used for
short and medium-length lines. For long lines, a
distributed-parameter model is used.
Short Line: Up to 80 kilometers (50 miles)
Z=R+jL
+
+
Vs Is VR
- -
Vs Z Is VR
58
Medium-Length Line: From 80 km (50 mi) to
240 km (150 mi)
Z=R+jL
+ +
I IR
Vs s 1
2
YC jC
1
2
1
2
YC VR
- -
2(23.8)3 30 ft
Dm 3
30
x1 0.121 ln 0.812 /mi
0.0373
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 118
59
Z1 0.1284 j0.812 /mi
Z1 12.84 j0.81.16
1.779 D -mile, to
xC x 106 ln m neutral
f r
1.779 30
x 106 ln 1
60 24 (1.108)
192,034 -mile, to neutral
XC 1,920.34, to neutral
1
YC j j520.74 x 10-6 -1
XC
Mutual Capacitance
In capacitance calculations, the earth is assumed
as a perfectly conducting plane. The electric field
that results is the same if an image conductor is
used for every conductor above ground.
+qb
Dab +qc
Dac
Dan +qn
+qa
Haa
Hab Hac Han
-qa
-qn
-qc
-qb
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 120
60
The voltage drop from conductor a to ground is
va vaa'
1
2
1 H H H
(qa ln aa qb ln ab ... qn ln an
4 ra Dab Dan
r D D
qa ln a qb ln ab ... qn ln an )
Haa Hab Han
Combining common terms, we get
1 H H H
va (qa ln aa qb ln ab ... qn ln an )
2 ra Dab Dan
…
…
61
where
1 H
Pkk ln kk
2 rk
1 H
Pkj ln kj
2 Dkj
Inversion of matrix P gives
62
For a purely capacitive system, we get
63
Transposed Three-Phase Line
qa Phase c
Pos.1
qb Phase a
Pos.2
qc Phase b
Pos.3
1
3
s 1
3
s 1
3
s
Section 1 Section 2 Section 3
64
Since the line is completely transposed, the mutual
and ground capacitances are averaged values.
Average Mutual Capacitance
65
Sequence Capacitance
Using matrix notation, we have
Iabc YabcVabc
Iabc jCabcVabc
From Vabc AV012 and Iabc A I012, we get
A I012 jCabcAV012
or
I012 jA 1CabcAV012
Thus, we have
C012 A 1Cabc A
66
Example: Determine the phase and sequence
capacitances of the transmission line shown. The
phase conductors are 477 MCM ACSR 26/7 whose
radius is 0.0357 ft. The line is 50 miles long and is
completely transposed.
14’ 14’
Calculate distances
a b c
Haa=Hbb=Hcc=80 ft
Hab=Hbc=81.2 ft 40’
Hac=84.8 ft
Find the P matrix
1 H
Paa Pbb Pcc ln aa
2 0 ra
For air,
1
0 x 10-9 Farad/meter
36
Substitution gives
80
Paa 18 x 109 ln
0.0357
138.86 x 109 Meter/Farad
86.29 x 106 Mile/Farad
Similarly, we get
1 H
Pab Pbc ln ab
2 0 Dab
19.66 x 106 Mile/Farad
67
The P matrix can be shown to be
68
The self- and mutual capacitances are
69
The voltage equation for section 1 is
va Paa Pab Pac Paw qa
vb Pba Pbb Pbc Pbw qb
=
vc Pca Pcb Pcc Pcw qc
vw Pw a Pw b Pw c Pw w qw
where vw=0 and the elements of P are
1 H
Pkk ln kk k=a,b,c,w
2 rk
1 H kj=ab,ac,aw,
Pkj ln kj
2 Dkj bc,bw,cw
70
Example: Determine the phase and sequence
capacitances of the transmission line shown. The
phase conductors are 477 MCM ACSR 26/7 whose
radius is 0.0357 ft. The overhead ground wires
are 3/8 inch E.B.B. steel. The line is 50 miles long
and is completely transposed. w u
Pw a Pw b Pw c Pw w Pw u
Pua Pub Puc Puw Puu
where rw ru 0.5( 38 )(12
1
) 0.016 ft
1 110
Pw w Puu ln
2 0 0.016
99.1 x 106 Mile/Farad
71
It can be shown that the P matrix is
P1 P2
P=
P3 P4
where
86.3 19.7 12.4
P1 = 19.7 86.3 19.7 x 106 mi/F
12.4 19.7 86.3
14.8 19.6
P2 = P3T = 19.6 19.6 x 106 mi/F
19.6 14.8
99.1 23.1
P4 = x 106 mi/F
23.1 99.1
72
Using Kron’s reduction,
1
Pabc P1 P2P4 P3
or
81.30 14.15 7.69
Pabc = 14.15 80.02 14.15 x 106 mi/F
7.69 14.15 81.30
Capacitances to ground
73
Assuming a complete transposition cycle,
Electrostatic Unbalance
In the system shown, assume Van , Vbn and Vcn
are positive-sequence voltages. The line is not
transposed.
a
Ia
a’
+
Van Ib Cab
b
Vcn - -
Vbn b’
+ + Cbc Cca
c
Ic
c’
+
+ + +
Vn Va Vb Vc C cg Cbg C ag
- - -
-
74
The charging
current
at no-load
condition
Iabc jCabcVabc jBabcVabc
Transforming to sequence components, we get
I012 jA 1CabcAV012 jB012V012
From KVL, we get
Van Va Vn Vbn Vb Vn
Vcn Vc Vn
or
Va n 1
Vabc = V
bn
+ 1 Vn
Vcn 1
75
Expanding, we get
Ia0 jB00Vn jB01Van
Ia1 jB10Vn jB11Van
Ia2 jB20Vn jB21Van
Case 1: Neutral Ungrounded
If the neutral is ungrounded, Ia0 0. Electrostatic
unbalance is defined in terms of the sequence
voltages. It can be shown that
V V C
d0 a0 n 00
Va1 Van C01
76
Case 2: Neutral Grounded
If the neutral is grounded, Vn 0. Electrostatic
unbalance is defined in terms of the sequence
currents. It can be shown that
2
Ia0 Cag a Cbg aCcg
d0
Ia1 3(Cg0 3Cm0 )
Note: Electrostatic unbalance may be reduced by
Transformers
Sequence Impedance
77
Ideal Two-Winding Transformer
NH NX
+ +
I IX
VH H VX
_ _
VH N
Voltage Equation: H
VX NX
The voltage drop from the polarity-marked
terminal to the non-polarity-marked terminal of
the H winding is in phase with the voltage drop
from the polarity-marked terminal to the non-
polarity-marked terminal of the X winding.
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 155
Current Equation: NH IH NX IX
The current that enters the H winding through the
polarity-marked terminal is in phase with the
current that leaves the X winding through the
polarity-marked terminal.
Note: Balancing ampere-turns must be satisfied at
all times.
78
Ideal Autotransformer
IH nS
VH nC
IC IX VX
79
Transformer Polarity
H1 H2 H1 H2
Subtractive Additive
X1 X2 X2 X1
Practical Transformer
1. The H and X coils have a small resistance.
2. There are leakage fluxes in the H and X coils.
3. There is resistance loss in the iron core.
4. The permeability of the iron is not infinite.
iH iX
+ +
vH eH eX vX
- NH NX -
iron
core
80
Equivalent Circuit
RH jXH Iex NH NX R X jXX
+ + + +
IH R IX
VH c jXm EH EX VX
- - - -
H winding Ideal X winding
+ + + +
IH IX
VH EH EX VX
- - - -
81
Three-Winding Transformer
NX
NH +
+ IX VX
IH _
VH +
_
NY IY VY
_
VH N VH N VX N
H H X
VX NX VY NY VY NY
NH IH NX IX NY IY
Three-Winding Autotransformer
+
IY V
NY
IH _Y
VH NH
NX IX VX
VH N VH N VX N
H H X
VX NX VY NY VY NY
NH IH NX IX NY IY
82
Equivalent Circuit
From 3 short-circuit tests, get
ZHX=impedance measured at the H side when the
X winding is short-circuited and the Y winding
is open-circuited
ZHY=impedance measured at the H side when the
Y winding is short-circuited and the X winding
is open-circuited
ZXY=impedance measured at the X side when the
Y winding is short-circuited and the H winding
is open-circuited
Note: When expressed in ohms, the impedances
must be referred to the same side.
ZX
ZH +
+
ZY +
VX
VH
- VY -
-
ZHX ZH ZX
ZXY ZX ZY
ZHY ZH ZY
or
ZH 1
2
(ZHX ZHY ZXY )
ZX 12 (ZHX ZHY ZXY )
ZY 1
2
(ZHX ZHY ZXY )
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 166
83
Example: A three-winding three-phase transformer
has the following nameplate rating:
H: 30 MVA 140 kV
X: 30 MVA 48 kV Y:
10.5 MVA 4.8 kV
Short circuit tests yield the following impedances:
ZHX = 63.37 @ the H side
ZHY = 106.21 @ the H side
ZXY = 4.41 @ the X side
Find the equivalent circuit in ohms, referred to the
H side.
ZXY (140
48
)2 (4.41) 37.52
2.66
66.03 +
+
+
VX
VH 40.18
- VY -
-
84
Angular Displacement
ANSI/IEEEC57.12.00-1993: The angular displace-
ment of a three-phase transformer is the time
angle (expressed in degrees) between the line-to-
neutral voltage of the high-voltage terminal
marked H1 and the the line-to-neutral voltage of
the low-voltage terminal marked X1.
Vector Diagrams
H2 X2 H2
X2
X1
X1 X3
H1 H3 H1 H3 X3
- Connection Y- Connection
H2 X2 H2 X2
X1
X1 X3 X3
H1 H3 H1 H3
Y-Y Connection -Y Connection
85
Wye-Wye Connection
HV side LV side
H2 X2
N n
H1 H3 X1 X3
VH2N
VX2n
VX1n VX3n
VH1N VH3N
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 171
VH1N VH3N VAN VCN
VX2n Vbn
Van and VAN
are in phase.
VX1n VX3n Van Vcn
86
Let H1=A, H2=B and H3=C. Also let X1=c, X2=a
and X3=b.
VH2N VBN
VH1N VH3N VAN VCN
VX2n Van
V
Van lags AN
by 120o.
VX1n VX3n Vcn Vbn
VH1N VH3N VAN VCN
VX2n Vcn
Van lags VAN
by 240o.
VX1n VX3n Vbn Van
87
Wye-Wye Connection
HV side LV side
H2 X2
N n
H1 H3 X1 X3
VH2N
VX3n VX1n
VH1N VH3N VX2n
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 175
VH1N VH3N VAN VCN
VX3n VX1n Vcn Van
Van lags VAN
by 180o.
VX2n Vbn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 176
88
Let H1=A, H2=B and H3=C. Also let X1=c, X2=a
and X3=b.
VH2N VBN
VH1N VH3N VAN VCN
VX3n VX1n Vbn Vcn
V
Van lags AN
by 300o.
VX2n Van
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 177
VH1N VH3N VAN VCN
VX3n VX1n Van Vbn
Van lags VAN
by 60o.
VX2n Vcn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 178
89
IEC Designation
The International Electrotechnical Commission
(IEC) uses the face of a clock in identifying the
phase shift. 0
IEC Designation for -
10 2
Dd0 Dd2 Dd4
Dd6 Dd6 Dd10
8 4
IEC Designation for Y-Y
Yy0 Yy6 6
90
Consider a set of balanced 3-phase line-to-neutral
voltages with asequence
a-c-b. From
Vab Van Vnb Vcn
we get
Vab Van Vbn
Van
Similarly, weget
Vbc Vbn Vcn Vbn
Vca Vcn Van
Vab Vcn Vca
Van leads Vab by 30o
Van
Vbn leads Vbc by 30o
Vbn
Vcn leads Vca by 30o
Vbc
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 181
91
Y- and -Y Transformers
H2 H2
X2 X2
X1 X1
H1 H3 X3 X3
H1 H3
HV Y – LV
H2 X2
H2 X2
X1
X1
X3
H1 H3
H1 H3
X3
HV - LV Y
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 183
HV Wye – LV Delta
H1 X1
H2 X2
N
H3 X3
VH2N
VX1X2 VX2n
Sequence
1-2-3 VX1n
VX2X3
VX3X1
VH1N VH3N
VX3n
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 184
92
The choice of phases A, B and C is arbitrary. Let
H1=A, H2=B, H3=C. Also let X1=a, X2=b, X3=c.
VH2N VBN
VH1N VH3N VAN VCN
VX2n Vbn
VX1X2 Vab
VX1n Van lags VAN
VX2X3 Van Vbc by 30o
VX3X1 Vca
VX3n Vcn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 185
VH1N VH3N VAN VCN
VX2n Vcn
VX1X2 Vbc
VX1n
V Van lags VAN
VX2X3 bn Vca by 270o
VX3X1 Vab
VX3n Van
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 186
93
Let H1=A, H2=B, H3=C. Also let X1=c, X2=a,
X3=b.
VH2N VBN
VH1N VH3N VAN VCN
VX2n Van
VX1X2 Vca
VX1n Van lags VAN
VX2X3 Vcn Vab by 150o
VX3X1
Vbc
VX3n Vbn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 187
HV Wye – LV Delta
H1 X1
H2 X2
N
H3 X3
VH3N VX3n
Sequence VX3X2
3-2-1 VX1X3 VX2n
VH1N VH2N VX2X1
VX1n
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 188
94
The choice of phases A, B and C is arbitrary. Let
H1=C, H2=B, H3=A. Also let X1=c, X2=b, X3=a.
VH3N VAN
VH1N VH2N VCN VBN
VX3n
Van
VX3X2
Vab
VX1X3 Van leads VAN
VX2n Vca Vbn by 30o
VX2X1
VX1n Vbc
Vcn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 189
VH1N VH2N VCN VBN
VX3n Vbn
VX3X2 Vbc
VX1X3 VX2n Vab V
Vcn an leads V AN
by 150o
VX2X1 Vca
VX1n Van
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 190
95
Let H1=C, H2=B, H3=A. Also let X1=b, X2=a,
X3=c.
VH3N VAN
VH1N VH2N VCN VBN
VX3n Vcn
VX3X2 Vca
Van lags VAN
VX1X3 VX2n Vbc Van by 90o
VX2X1 Vab
VX1n Vbn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 191
C H3 X3 b A H3 X3 b
B H2 X2 a B H2 X2 c
A H1 X1 c C H1 X1 a
Van lags VAN by 150o Van lags VAN by 210o
A H3 X3 c C H3 X3 a
B H2 X2 a B H2 X2 c
C H1 X1 b A H1 X1 b
Van lags VAN by 90o Van lags VAN by 270o
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 192
96
IEC Designation
The International Electrotechnical Commission
(IEC) uses the face of a clock in identifying the
phase shift.
11 1
IEC Designation for Y-
HV Delta– LV Wye
H3 X3
H2 X2
n
H1 X1
VH2N VX2n
Sequence
VH1H2
1-2-3 VH2H3 VX n
1
VH1N
V
VH3H1 H3N VX3n
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 194
97
HV Delta– LV Wye
H3 X3
C c
H2 X2
B b
n
H1 X1
A a
VBN Vbn
VAB VBC Van Van lags VAN
by 30o
VAN VCA VCN Vcn
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 195
Positive–Sequence Voltages
Consider a Y- transformer with a 30o phase shift.
B b
H2 X2
X1
N
A H3 X3 c
H1
C a
VBN1 Vbn1
(A-B-C) Vab1
Van1 lags VAN1
Van1 Vbc1
by 30o
Vca1
VAN1 VCN1
Vcn1
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 196
98
Negative–Sequence Voltages
Consider a Y- transformer with a 30o phase shift.
B b
H2 X2
X1
N
A H3 X3 c
H1
C a
VCN2 Vcn2
Van2 leads VAN2
(A-C-B) Vcb2
by 30o
Vac2 Vbn2
Vba2
VAN2 VBN2 Van2
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 197
Positive–Sequence Currents
Consider a Y- transformer
with a 30o phase shift.
B H2
Iba1 b
X2
IB1 X1
Ib1
Icb1
IA1 Ic1
A H3 Iac1 X3
c
IC1 H1 Ia1
C a
Icb1 Ib1
IB1 (A-B-C)
Ia1 lags IA1
Ia1 by 30o
IA1 IC1 Iba1 Iac1
Electrical & Electronics Engineering Institute
Ic1 EE 251 – Fault Studies
University of the Philippines Diliman 198
99
Negative–Sequence Currents
Consider a Y- transformer
with a 30o phase shift.
B H2
Iba2 b
X2
Ib2
IB2 X1
Icb2
IA 2 Ic2
A H3 Iac2 X3
c
IC2 H1 Ia2
C a
Ic2 Iac2
IC2
(A-C-B) Ia2 leads IA 2
Ib2 by 30o
IA 2 IB2 Iba2 Icb2
Ia2
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 199
100
The Base Currents
Base MVA3 x 1000
Base I
3 Base kVL
100,000
251 A at the 230 kV side
3 (230)
100,000
2,887 A at the 20 kV side
3 (20)
The Per-Unit Phase Fault Currents
IA 3.0 p.u.
IB IC 0
2
IA1 3 (IA a IB a IC ) 1.0 p.u.
1
IA2 13 (IA a2 IB a IC ) 1.0 p.u.
The Sequence Fault Currents at the 20-kV side
Ia0 0
Ia1 IA1 30o 1.0 - 30o p.u.
Ia2 IA230o 1.030o p.u.
101
The Line Currents at the 20-kV side
Ia Ia0 Ia1 Ia2 1.732 p.u.
Ib Ia0 a2 Ia1 a Ia2 1.732 p.u.
2
Ic Ia0 a Ia1 a Ia2 0
The Line Currents in Amperes
Ia Ib 5,000 A I 0
c N
Check ampere turns: NH IH NX IX or H IH IX
NX
230 / 3
(753) 5,000 (Check)
20
Shell Type
102
Three-Legged Transformer Core
The 3-legged core type three-phase transformer
uses the minimum amount of core material. For
balanced three-phase condition, the sum of the
fluxes is zero.
30
Note: The zero-sequence
flux leaks out of the core
and returns through the
0 0 0
transformer tank.
3-Legged Core Type
103
Positive Sequence Impedance
Whether a bank of single-phase units or a three-
phase transformer unit (core type or shell type),
the equivalent impedance is the same.
Using per-unit values, the positive-sequence
equivalent circuit is
Z1 R1 jX1
+ + Note: The negative-
I I sequence impedance is
VH H X VX equal to the positive-
- - sequence impedance.
Z0 Z1
+
+
VH VX
- -
104
Zero Sequence Impedance*
Transformer Connection Zero-Sequence Network
Z0 Z1
+
+
VH VX
- -
Z0 Z1
+
+
VH VX
- -
Z0 Z1
+
+
VH VX
- -
105
Zero Sequence Impedance*
Transformer Connection Zero-Sequence Network
XH XX
+
+
VH Xm VX
- -
XH
+
+
VH Xm VX
- -
XH XX
+
+
VH Xm VX
- -
106
Zero Sequence Impedance*
Transformer Connection Zero-Sequence Network
XH XX
+
+
VH Xm VX
- -
XH XX
+
+
VH Xm VX
- -
107
Example: A three-winding three-phase transformer
has the following nameplate rating: 150/150/45
MVA 138YG-69YG-13.8 kV (Yy0d1).
H-X @ 150 MVA = 14.8%
H-Y @ 45 MVA = 21.0%
X-Y @ 45 MVA = 36.9%
Draw the positive and zero-sequence equivalent
circuits. Use 100 MVA and the transformer voltage
ratings as bases.
At the chosen MVA base,
ZHX 0.148(100 / 150) 0.10 p.u.
ZHY 0.21(100 / 45) 0.47 p.u.
ZXY 0.369(100 / 45) 0.82 p.u.
We get
ZH 12 (0.10 0.47 0.82) 0.125 p.u.
ZX 12 (0.10 0.47 0.82) 0.225 p.u.
ZY 12 (0.10 0.47 0.82) 0.595 p.u.
Positive/Negative Zero Sequence
Sequence Network Network
ZX ZX
ZH + ZH +
+ +
ZY + VX ZY + VX
VH VH
-
VY - -
VY -
- -
108
Generators
Sequence Impedance
109
Negative Sequence Impedance
IEEE Std 115-1995: Determination of the negative-
sequence reactance (Method 1)
The machine is operated at rated speed with its
field winding short-circuited. Symmetrical
sinusoidal three-phase currents of negative phase
sequence are applied to the stator. Two or more
tests should be made with current values above
and below rated current, to permit interpolation.
The line-to-line voltages, line currents and electric
power input are measured and expressed in per-
unit.
110
Zero Sequence Impedance
IEEE Std 115-1995: Determination of the zero-
sequence reactance (Method 1)
The machine is operated at rated speed with its
field winding short-circuited. A single-phase voltage
is applied between the line terminals and the
neutral point.
Measure the applied
voltage, current and V
electric power.
E A W
Field
111
Positive Sequence Impedance
The AC RMS component of the current following a
three-phase short circuit at no-load condition with
constant exciter voltage and neglecting the
armature resistance is given by
E E E t
I(t) exp
Xds Xd ' Xds '
d
E E t
exp
X
d " X '
d "
d
where E = AC RMS voltage before the short circuit.
112
IEEE Std 115-1995: Determination of the transient
and subtransient direct axis reactances (Method 1)
The direct-axis transient reactance is determined
from the current waves of a three-phase short
circuit suddenly applied to the machine operating
open-circuited at rated speed. For each test run,
oscillograms should be taken showing the short
circuit current in each phase.
The direct-axis transient reactance is equal to the
ratio of the open-circuit voltage to the value of the
armature current obtained by the extrapolation of
the envelope of the AC component of the armature
current wave, neglecting the rapid variation during
the first few cycles.
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 225
113
Example: Calculation of transient and subtransient
reactances for a synchronous machine
Phase 1 Phase 2 Phase 3 Ave
(1) Initial voltage 1.0
(2) Steady-state Current 1.4 1.4 1.4
(3) Initial Transient Current 8.3 9.1 8.6
(4) I’ = (2)+(3) 9.7 10.5 10.0 10.07
(5) Xd’ = (1)(4) 0.0993
(6) Init. Subtransient Current 3.8 5.6 4.4
(7) I” = (4)+(6) 13.5 16.1 14.4 14.67
(8) Xd” = (1)(7) 0.0682
14
12 +
10 +++ Curve B
Current in phase 1 (per unit)
+
8 ++
+ ++
++
6 + ++
Line C
++
5 ++
++
+ +
4 ++
+ +
+ + ++
3
+
2.0
+ Line A
1.5
+
+
1.0
+
0.8
0.6 +
+ Curve A
0.4 +
0 10 20 30 40 50 60
Time in half-cycles
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 228
114
Average Generator Reactances
Reactance Turbo Water-Wheel Synchronous
Generators Generators Motors
Xd 1.10 1.15 1.20
Xq 1.08 0.75 0.90
X d‘ 0.23 0.37 0.35
X q‘ 0.23 0.75 0.90
X d” 0.12 0.24 0.30
X q” 0.15 0.34 0.40
X2 0.12 0.24 0.35
Ungrounded Generator
If the generator is connected ungrounded-wye or
delta, no zero-sequence current can flow. The
sequence networks for the generator are shown
below.
F1 F2 F0
jX1
+ jX2 jX0
Eg
-
N1 N2 N0
Positive Negative Zero
Sequence Sequence Sequence
115
Grounded-Wye Generator
The sequence networks for the grounded-wye
generator are shown below.
F1 F2 F0
jX1
+ jX2 jX0
Eg
-
N1 N2 N0
F1 F2 F0
jX1 jX0
+ jX2
Eg 3Zg
-
N1 N2 N0
116
Electric Utility (Power Grid)
Thevenin’s Theorem
Any linear active network with output terminals AB can
be replaced by a single voltage source Vth in series with
a single impedance Zth
A A
Linear + Zth
Active Vth
Network -
B B
117
Utility Fault MVA
Electric Utilities conduct
Electric short circuit analysis at the
Utility Connection Point of their
Grid customers
IF
Fault Customers obtain the Fault
Data at the Connection Point
to represent the Utility Grid
for their power system
analysis
Customer
Facilities
118
Equivalent Circuit of Utility
I TPF
Vf
S TPF V f I TPF
V
f
2
Z1 Z1
Where, Z1 and Z2 are the
equivalent positive-
Z1
kV
2
sequence and
Z 2 negative-sequence
Fault MVA3 impedances of the
utility
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 237
I SLGF
3V f
S SLGF V f I SLGF
3Vf
2
Z1 Z 2 Z0 2Z 1 Z 0
Z1 Z 2
2Z 1 Z 0
3Vf
2
119
Equivalent Circuit of Utility
Example:
Z1 Z 2
kV 2
69
2
1.3603
Fault MVA3 3500
In Per Unit,
Z actual 1.3603
Z1 Z 2 0.0286 p.u.
Z base 47.61
or
100MVABASE
Z1 Z 2 0.0286 p.u.
3500MVAFAULT
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 240
120
Equivalent Circuit of Utility
Solving for the Resistance and Reactance,
2Z 1 Z 0
3Vf
2
31.0
2
0.1
S SLGF 30
Re al2 Z 1 Z 0 0.1cos tan -1 20 0.004994 p.u.
Im ag2 Z 1 Z 0 0.1sintan -1
20 0.099875 p.u.
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 242
121
Equivalent Circuit of Utility
2 Z 1 Z 0 0.004994 j0.099875
Z 0 0.004994 j0.099875 20.0013 j0.028571
0.003694 j0.042733 p.u.
+ + +
0.0013+j0.028571
+ 0.0013+j0.028571 0.003694 j0.042733
V f 10
-
- - -
Positive Negative Zero
Sequence Sequence Sequence
122
Equivalent Circuit of Utility
The equivalent sequence networks of the Electric
Utility Grid are:
+ + +
R1 +jX1
+ R2 +jX2 R0 +jX0
Eg
-
- - -
Positive Negative Zero
Sequence Sequence Sequence
Utility Thevenin Equivalent Circuits
Two-Port Network
123
Two-Port Network
I1 I2
+ RLC +
V1 V2
-
Network -
I1 I2
Impedance or Z parameters
V1 z11 z12 I1
= or V = ZI
V2 z21 z22 I2
I1 I2
+ RLC +
V1 V2
- Network -
I1 I2
Admittance or Y parameters
I1 y11 y12 V1
=
I2 y 21 y 22 V2
or
I = YV
124
I1 I2
+ RLC +
V1 V2
- Network -
I1 I2
Hybrid or H parameters
V1 h11 h12 I1
=
I2 h21 h22 V2
or
U = HW
125
V2
z21 5
I1 I2 0
Similarly, we get
V1
z12 5
I2 I1 0
V2
z22 9
I2 I1 0
Note: z12=z21
126
Example: For the electric circuit shown, find the Y
parameters.
2 4
+ +
I1 I2
V1 5 V2
- -
Let V2=0.
2 4
+ +
I1 I2
V1 5 V2=0
- -
V1 5(4) 38
zin 2
I1 54 9
or
9 1
y11
38
Using current division, we get
5 5
I2 I1 I1
54 9
127
Since I1=y11V1 when V2=0, we get
5
I2 y11V1
9
or
I2 5 5 1
y21 y11
V1 9 38
I1 I2
y12 and y22
V2 V1 0
V2 V1 0
Let V1=0.
2 4
+ +
I1 I2
V1=0 5 V2
- -
V2 5(2) 38
zout 4
I2 52 7
or
7 1
y22
38
Using current division, we get
5 5
I1 I2 I2
52 7
128
Since I2=y22V2 when V1=0, we get
5
I1 y22V2
7
or
I1 5 5 1
y12 y22
V2 7 38
We can also get the Y parameters by taking the
inverse of the matrix of the Z parameters.
7 5 1 9 5
Z= Y=Z-1= -1
5 9 38 5 7
Note: y12=y21
+ RLC & +
V1 V2
-
Sources -
I1 I2
Assume all internal sources are independent. With
all sources set to zero, the network is described by
either (1) V = ZI
(2) I = YV
(3) U = HW
129
With all sources active, the network will be
described, respectively, by
(1) V = ZI + Vs
(2) I = YV + Is
(3) U = HW + Us
where
130
+
EG1 - +
I1 EG2 I2
V1 + V2
I1 j0.27 I2
- -
+
+
I1 I2
V1 j0.27 V2
I1 I2
- -
V1 z11 I1 z12 I2 and V2 z21 I1 z22 I2
V1
z11 j(0.5 0.27 0.22) j0.99
I1 I 0
2
V
z21 2 j0.27
I1 I 0
2
V
z22 2 j0.27
I2 I 0
1
131
When the sources are active, we get
V1 z11 z12 I Vz1
= 1 +
V2 z21 z22 I2 Vz2
If I1 I2 0, then Vz1 V1 and Vz2 V2 . By
inspection, we get
Vz1 EG2 EG1 1.0 1.130o
0.552 85.08o
Vz2 EG2 1.00o
Series Connection
Two two-port networks are connected in series if
their input currents are the same and their output
currents are likewise the same.
+ Network +
I1 V1a V2a I2
A
V1 V2
I1 Network I2
- V1b V2b -
B
132
For two-port network A, we get
Remarks:
(1) For two two-port networks in series, the
over-all Z parameters is the sum of the
individual Z parameters.
(2) This result applies to any number of two-
port networks in series.
(3) Series-connected two-port networks must
be described using the Z parameters.
133
Parallel Connection
Two two-port networks are connected in parallel if
their input voltages are the same and their output
voltages are likewise the same.
+ Network +
I1 I1a I2a I2
A
V1 V2
I1 Network I2
- I1b I2b -
B
134
Adding the two-port network equations, we get
Remarks:
(1) For two two-port networks in parallel, the
over-all Y parameters is the sum of the
individual Y parameters.
(2) This result applies to any number of two-
port networks in parallel.
(3) Parallel-connected two-port networks must
be described using the Y parameters.
Hybrid Connection
Two two-port networks are connected in hybrid if
one port is connected in series and the other port
is connected in parallel.
+ Network +
I1 V1a I2a I2
A
V1 V2
I1 Network I2
- V1b I2b -
B
135
For two-port network A, we get
Remarks:
(1) For two two-port networks in hybrid
connection, the over-all H parameters is
the sum of the individual H parameters.
(2) This result applies to any number of two-
port networks in hybrid connection.
(3) Hybrid-connected two-port networks must
be described using the H parameters.
136
Electrical & Electronics Engineering Institute EE 251 – Fault Studies
University of the Philippines Diliman 273
137