Power Systems K-Notes
Power Systems K-Notes
Transmission Lines
Skin Effect
It is tendency of AC current to be concentrated on the surface of conductor.
Cause: Non-uniform distribution of magnitude flux linkages Due to skin effect, the effective
area of cross section of conductor decrease and hence resistances increases.
In case of DC, There is no skin effect so
RDC R AC
0 r
8
External inductance from distance d1 to d2
Internal inductance
Lex
0r d2
ln
2 d1
0 r 0 r d2
ln
8
2 r
d d
0 r ln 1 0 r ln
2 re 4 2 r
r 0.7788r = Geometric mean radius (GMR)
Total inductance
0 r
d
ln
2 r
Total inductance = L1 L2
Lsys
0 d
ln
r
If radius of both wire is not same, assume radius of 1st wire ra & that of second wire is rb
Lsys
0 d
ln
r r
a b
ra 0.7788ra &
rb 0.7788rb
If instead of a single conductor per phase we use multiple conductor, then GMR is
replaced by self GND (Geometric Mean Distance) and d by mutual GMD.
Self GMD
self GMDfwd D11 D12 ........D1n D21 D22 .........D2n ....... Dn1 Dn2 ........Dnn
Where
Dii ri 0.7788ri
Self GMDbwd Di 1
' '
Where
m2
Dii ri 0.7788ri
Mutual GMD
Mutual GMD D11' D12' ...........D1m' ................. Dn1' Dn2' .................Dnm'
Now, with these terms all the inductance expressions change to
Single wire :
1 , 2 wire:
0
d
ln
2 self GMD
0
mutual GMD
ln
2 Self GMD
Lph
0r D
ln
2 r
r 0.7788r
mn
n2
Asymmetrical configurations
If conductors are placed horizontally or vertically.
Lph
0r Deq
ln
2 r
Dbc Dbceq
Dca Dcaeq
Dab eq Da1b1
3 12 12 3
6m
4.24m
7.348m
Dbc eq 3 6 6 3 4
Dca eq 6 9 9 6 4
= 5.71m
5
Self GMD r D
a
a1a2
Da2a1 r
0.341m
Self GMD r D
Db2b1 r
b1b2
0.2467m
Self GMD r D
c
c1c2
Dc2c1 r
0.1528m
=
L
0.2398m
0 GMD
ln
2 GMD Self
5.71
2 10 7 ln
0.634mH / km
0.2398
Capacitance
Single Phase 2 Wire System
Cab
0 r
D
ln
rr
12
Can
20 r
D
ln
r1
Cbn
20 r
D
ln
r2
20 r
GMD
ln
20 0
GMD
ln
Self GMD
or
Where f = frequency
2) Medium Line
80 km < l < 200 km
4000 < l*f < 10000
3) Long Line
l > 200 km
l*f > 10000
Modeling of transmission lines
Transmission lines are modeled as 2 port network
Vs AVR BIR
Is CVR DIR
Under no load
IR 0 , Vs AVR , VR
Vs
A
No Load Voltage
Vs
A
and IR 0 , Is CVR
Vs
VR
A
Voltage Regulation
100%
VR
This current is called as line charging current and is responsible for as effect is called as Ferranti
Effect.
Ferranti Effect
Under no-load or light load conditions receiving end voltage becomes more than sending end
voltage due to presence of line charging current.
Short transmission line
Vs VR IR R jwL
VR IR Z
Vs 1 z VR
I 0 1 I
s
R
A = D (symmetrical)
AD BC = 1 (reciprocal)
Approximate Voltage Regulation
For lagging pf
VR=
IR
VR
R cos
X sin R
R cos R
X sin R
load pf=cosR
For leading pf
VR
IR
VR
YZ
YZ
1
Z 1
Vs
2
4 VR
YZ IR
Is
Y
1
2
Here all problems are in actual values & not per unit length.
Nominal--Model
YZ
1
Z
Vs
V
2
R
Is Y 1 YZ 1 YZ IR
4
2
Vx VR cosh X IR Z c sinh x
IR IR cosh x
VR
sinh x
Zc
Where Vx & I x are voltage and current at distance x from receiving end.
Zc
R jL
= surge impedance
G jC
VR
sinh l
Zc
cosh l Z c sinh l
Vs
VR
1
cosh l IR
Is Z
A = D (symmetric)
AD - BC = 1 (reciprocal)
Power Transfer Equation
Pr
Vs Vr
B
Vs Vr
Qr
2
A
Vr cos
B
cos
2
A
Vr sin
B
sin
For Short TL
B Z
B Z ;
A 10 A
Pr
Vs Vr
Qr
A 1, =0
cos
Vs Vr
Z
sin
Vr
cos
Z
Vr
sin
Vs Vr
X
sin ; QR =
Vs Vr
X
cos
Vr
Remember, the last expression can be applied between any two bases in a power system as long
as transmission line connecting them is loss less.
Wave Propagation
Due to continuous energy transfer between L & C elements of a transmission line we consider
energy propagation from sending to receiving end & hence wave propagation.
Z c = Surge impedance or characteristic impedance.
10
= Propagation constant.
R jL G jC
j LC
j
= attenuation constant
= phase constant
LC
Velocity of wave
Wavelength
1
LC
2
2
LC
PL
2
Vr(L
L)
ZC
If ZL ZC
VR IR Z C
11
Surge Traversal
Theremin equivalent circuit
When surge voltage V is induced on the line & line can be represented as Theremin
equivalent circuit shown.
ZC = Characteristic impedance of line
Incident Voltage V1 V
Reflected voltage V
V V V2
[Voltage continuity]
Z ZC
V V L
ZL Z C
Reflected current
V
V
V
, Refracted current 2 , Incident current
ZC
ZC
ZL
Reflection coefficient:
V ZL Z C
V ZL Z C
Refraction coefficient:
V2
2ZL
V ZL Z C
Voltage Control
Usually in case of lagging loads, the voltage at receiving end falls below sending end voltage
and to boost the receiving end voltage we connect a shunt capacitor at receiving end.
Similarly, in case of leading loads, receiving end voltage is higher than sending end voltage
so we connect a shunt reactor to avoid over-voltage.
Usually in GATE, we need to calculate rating of capacitor for voltage control & it is illustrated
through a question shown below:
12
Example: A three phase overhead lines has a resistance & reactance of 5 & 20 respectively.
The load at receiving end is 30MW, 0.85 pf lagging at 33kv & we connect a compensating
equipment at receiving end to maintain voltage at each end equal to 33 Kv. Find rating of
compensating equipment?
Solution:
Assuming base (MVA) = 30 MVA
Base voltage = 33 kv
pu power =
30Mw
1pu
30MVA
Base impedance =
V2
36.3
S
pu impedance =
5 j20
0.56875.960
36.3
PR
VS VR
Z
cos
VR
Z
cos
11
12
cos 75.96
cos75.96
0.568
0.568
Z 0.568 ; =75.96 0
40.110
QR
VS VR
Z
sin
VR
Z
sin
1 1
12
sin 75.96 40.11
sin75.96
0.568
0.568
QR = - 0.645 pu
QL PL tan 1 tan cos1 0.85
QL = 0.6197 pu
QR QL QC QC 1.2647pu
QC 1.2647 30 = 37.9 MVAR
QC P1 tan 1 P2 tan 2
In case of capacitor bank, PC 0 P1 P2
If we wish to calculate capacitance per phase
(in both voltage control & pf correction )
2
QC 3CVph
QC
2
3Vph
Underground Cables
Insulation resistance
R
R
ln
2l r
L = length of cable
R = Outer radius (sheath radius)
R = conductor radius
Capacitance Model
CC : Core capacitance
14
Calculating CS & CC
1)
C1 CS 2CC
CS
C2
3
C2
3
C
C
CC 1 2
6
2
2CC C1
Cph CS 3CC
3C1 C2
2
6
C3
3CC CS
2
2
Cph 2C3
15
1
cphR
R = Insulation resistance
Overhead Insulators
For suspension type string insulator, the model for 3-discs looks like as shown.
C
Let m m
CS
V2 V1 1 m
V3 V1 1 m2 3m
string voltage
No. of discs voltage across bottom disc
V2 V3
3 V3
100%
Distribution Systems
VA VB I A I1 r1 I A I1 I2 r2 I A I1 I2 I3 r3
16
pu value=
Out of these, 2 value must be known, to convert entire system into pu system.
Sbase Vbase Ibase
Ibase
Sbase
Vbase
Z base
Vbase
Ibase
2
Vbase
Sbase
3 - System
Ibase ph
Ibase
Sbase
base
Ibase ph Ibase / 3
Sbase
In per unit system, equivalent impedance of transformer referred to primary or
secondary in same.
17
Change of base
If base of system is changed from Vbase old , Sbase old to Vbase new ,Sbase new
2
YBUS matrix
YBUS
y10 y12
y12
0
y12
y 20 y12 y 23
y 23
y 23
y 23 y 30
OBSERVATIONS
1)
2)
3)
4)
The diagonal elements are sum of all admittance connected to that particular bus.
The off-diagonal elements are negative of admittance connected between two buses.
If two buses are not connected to each other than that elements is zero.
YBUS Matrix is a symmetrical matrix.
ZBUS matrix
ZBUS YBUS
18
V1
I1
.
.
.
.
Vk
Ik
. Z .
BUS
.
.
Vn
In
V
BUS
0
0
.
.
ZBUS
I
f
.
.
0
Vk If Zkk
If
Vf
Z 0
Zkk f
Vj Z jk I"f
Vj
Z jk
Zkk
Vf
Vjf Vj
Z jk
V
Zkk f
Eg Vjf
jX"d
Classification of buses
At each bus, there are 4 parameter: V , ,P,Q .
At any bus, out of these 4 quantities any 2 are specified.
1) Slack Bus / Swing Bus/Reference Bus
P, V is specified
Q, are unspecified
19
2n m 2 2n m 2
Incremental cost
If is cost required to generate an additional unit of energy.
IC i = Incremental cost of ith generator
IC i
Ci
PGi
Transmission Loss
m
PL PP
i jBij
i 1 j 1
th
th
Pi , Pj : Real power injection at i & j buses
B ij = loss coefficient
P
1 L
PGi
IC i Li cons tant
20
From this expression, for m generator we get ( m 1) equation and mth equation is
m
P
i 1
Gi
PD PLOSS
Fault Analysis
Symmetrical Components
For an unbalanced 3 phase system, the analysis is done better by means of symmetrical
components.
Va0 Zero sequence components
Va1 Positive sequence components
Va2 Negative sequence components
Va0
1 1
1
Va1 3 1
V
1 2
a2
1
Vs A Vp
1 Va
2 Vb
Vc
Where e j120
Vp = phase voltage
Vs = Symmetrical component
1 1
A 1 2
1
Z1 jX"d
21
Transformers
Depending upon scheme of connection, we close series or shunt connection & method
of grounding.
22
Shunt connection are closed for delta connection & series connection are closed for star
connection with grounded neutral. If primary & secondary are inter changed then circuit
becomes mirror image.
Case 1
If neutral is grounded through an impedance Z n
Z 0 Z T0 3Zn
Case 2
Case 3
Case 4
23
Case 5
Transmission Lines
Z0 Zs 2Zm
Z s = Self impedance
Zm = mutual impedance
Remembers, all sequence networks are always drawn in per unit & never in actual values.
Fault Analysis
The following short circuit faults are considered
1.
2.
3.
4.
24
3-phase short circuit fault comes under the category of symmetrical SC fault whereas other 3
faults are called as unsymmetrical SC fault.
Order of severity
LG < LL < LLG < 3 - SC
But if faults occur at terminal of alternator then LG fault is most severe.
Occurrence of SC fault
LG > LL > LLG > 3 - SC
Transient on a Transmission line
Equivalent Circuit
i t it iss
Rt
Vm
V
sin e L m sin wt
Z
Z
L
2
; Z R L
R
tan1
imm
2Vm
sin
Z
If resistance is neglected,
imm
900
2Vm
cos
Z
25
Finally, when current in field winding is also reduced to zero, we enter steady state period &
reactance is called as steady state reactance Xd .
Xd " X d ' X d
We replace alternators by an emf source in series with sub transient reactance and emf
source under no-load is usually 100 pu or terminal voltage in pu.
V 00
Ef t
pu
Vbase
Ef
Z eq Z f
Z f = fault impedance.
SC MVA
SC MVA
1
Zeq Z f
MVA base
Zeq Z f
pu
MVA
Usually, circuit breakers are selected on the basis of most severe fault which is 3 phase
SC fault.
Three ratings of circuit Breakers are important.
2) Making current
Making current = 2.54 Isc
26
Factor
1.0
1.1
1.2
1.3
Ia1
Ea
Z1 Z 2 Z 0 3Z f
Fault current I a
Ia 3Ia1
Ia
3 Ea
Z1 Z 2 Z 0 3Z f
SC MVA 3Ea1 , I a1 *
27
Z1 , Z 2 & Z 0
In pu SC MVA
Z1 Z 2 Z 0 3Z f
3
pu
Z1 Z 2 Z 0 3Z f
3 Ea1 I2
3 MVA base
Z1 Z 2 Z 0 3Z f
MVA
Ia1
Ea
Z1 Z 2 Z f
Fault current
Ib
j 3 Ea
Z1 Z 2 Z f
3
pu
Z1 Z 2 Z f
3 MVA base
Z1 Z 2 Z f
MVA
Ia0
Ea
Z1 Z 2 Z 0 3Z f
I Z2
a1
Z2 Z0
28
3
Z1 Z 2 1 3Z f Z 0
pu
3 MVA base
Z1 Z 2
Z 0 3Z f
MVA
Md2
Pm Pe
dt2
M = inertia constant ( MJ-S / elect - rad)
Pm= mechanical input (MW)
Pe = electrical output (MW)
= rotor angle
Another Form
H d2
Pm Pe
f dt2
H = inertia constant ( MJ / MVA)
Pm & Pe both are in pu
29
GH
(MJ S / elect - deg)
180f
GH
M
(MJ S / elect - rad)
f
G = machine rating (MVA)
M
If two alternators are swinging coherently. Then they can be replaced by a single
alternator having
Meq M1 M2
But H cannot be added directly, they must first be on same base.
Accelerating Power,
Pa Pm Pe
In steady state Pm Pe
In transient, Pm Pe so rotor accelerate or decelerate.
Before Fault
We say maximum power transferrable is Pmax,1
& Pe Pmax,1 sin
During fault
We say maximum power transferrable is Pmax,2
Pe Pmax,2 sin
After Fault
We say maximum power transferrable is Pmax,3
Pe Pmax,3 sin
30
Pmax,2 0
Pmax,3 Pmax,1
Cr clearing angle
By equal area criteria
2
Pm Pmax,1 sin d 0
P
0 sin1 m
Pmax,1
For critical clearing 2 max
max 0
tCr
2H Cr 0
= Critical Clearing Time
f
Pm
Pmax,1
E V
Xg X1 X2
During Fault
Pmax,2
E V
Xeq
After Fault
Pmax,3
E V
Xg X1
31
P
m
Pmax,1
0 sin1
P
m
Pmax,3
max sin1
Pad 0
0
Pm 0 d Pm Pmax,3 sin d 0
tCr
2H Cr 0
f
Pm
Pad 0
32
cr cos 1
This is a generic formula and can be applied to other two cases as well after substituting
value of Pmax,1 , Pmax,2 & Pmax,3 .
But tcr can only be calculated from cr in previous two cases using expression written
before.
Differential Relays
Nr I1 I2
K I1 I2 K
Ipu
N0
2
The ratio operating coil current & restraining coil current is called as bias of differential
relay.
33
Differential relays are used for protection of large transformers and CT are always
connected in configuration opposite to power transformer.
Example: 1) If power transformer is Y then CT is Y
: 2) If power transformer is then CT is Y Y
Buccholz relay used to prevent any incipient fault below oil level in a transformer of small
KVA.
34