Power System Parameters and Equations
Power System Parameters and Equations
Power System Parameters and Equations
Network Calculations
1
4/10/2013
2
4/10/2013
3
4/10/2013
4
4/10/2013
5
4/10/2013
In matrix form,
6
4/10/2013
7
4/10/2013
1 .5 ∠ 0 0
I1 = I 3 = = 1 . 2 ∠ − 90 0 = 0 − j1 . 20 per unit
j1 . 25
1.5∠ − 36.87 0
I2 = = 1.2∠ − 126.87 0 = −0.72 − j0.96 per unit
j1.25
8
4/10/2013
Numerical Methods
Direct Methods
a) Cramer’s Rule
b) Matrix Inversion
c) Gauss Method
d) Gauss-Jordan Method
Iterative Methods
a) Gauss Iterative Method
b) Gauss-Seidel Method
c) Newton-Raphson Method
9
4/10/2013
Solution:
[I] = [Ybus] [V]
Therefore:
[V] = [Ybus]-1 [I]
10
4/10/2013
Therefore:
V1 1.4111 - j0.2668
V2 = 1.3830 - j0.3508
V3 1.4059 - j0.2824
V4 1.4009 - j0.2971
11
4/10/2013
IA K L VA
=
IX LT M VX
Where:
IX – the submatrix composed of the currents entering the
nodes to be eliminated. Every element here is zero.
VX – the submatrix composed of voltages of the nodes to be
eliminated.
K – self and mutual admittances identified with the nodes to
be retained.
M – self and mutual admittances identified with the nodes to
be eliminated.
L, LT – are composed of the mutual admittances common to
nodes to be retained and to be eliminated.
IA K VA + L VX
=
IX LT VA + M VX
12
4/10/2013
VX = - M-1 LT VA
IA = K VA - L M-1 LT VA
[ IA ] = [ K - L M-1 LT ] [VA ]
[ IA ] = [ Ybus ] [VA ]
Therefore the new bus admittance matrix:
IA K L VA
=
IX LT M VX
13
4/10/2013
14
4/10/2013
15
4/10/2013
16
4/10/2013
Ykn Ynj
Ykj(new) = Ykj(old) -
Ynn
Ykn Ynj
Ykj(new) = Ykj(old) -
Ynn
17
4/10/2013
E th 1.432 ∠ − 11.97 0
IC = = = 0.316 ∠ 78.03 0 per unit
Z th + Z C j0.4733 − j5.0
18
4/10/2013
19
4/10/2013
Ik
k Orig. network
with bus k and
the reference
Ip
bus extracted
p
Zb
0 0
Ik Orig. network
k with bus k
and the
Ip reference bus
p extracted
Zb
0 0
V1 0 I1
V2 0 I2
Zbus(orig)
: = : :
Vn 0 In
Vp 0 0 ⋅⋅⋅ 0 Zb Ip
Zbus(new)
20
4/10/2013
Ik Ik + Ip
k Orig. network
with bus k and
the reference
Ip
bus extracted
p
Zb
Ik Ik + Ip
k Orig. network
with bus k and
the reference
Ip
bus extracted
p
Zb
V1 Z1k I1
V2 Z2k I2
Zbus(orig)
: = : :
Vn Znk In
Vp Zk1 Zk2 ⋅ ⋅ ⋅ Zkn Zkk + Zb Ip
Zbus(new)
21
4/10/2013
Ik Ik + Ip
k Orig. network
with bus k and
the reference
Zb bus extracted
Ib
0 0
Ik Ik + Ip
k Orig. network
with bus k and
the reference
Zb bus extracted
Ib
0 0
V1 Z1k I1
V2 Z2k I2
Zbus(orig)
: = : :
Vn Znk In
0 Zk1 Zk2 ⋅ ⋅ ⋅ Zkn Zkk + Zb Ib
Zbus(new)
22
4/10/2013
V1 Z1k I1
V2 Z2k I2
Zbus(orig)
: = : :
Vn Znk In
0 Zk1 Zk2 ⋅ ⋅ ⋅ Zkn Zkk + Zb Ib
Zbus(new)
Ij Ij + Ib
j Orig. network
with buses j, k
Ib Zb and reference
bus extracted
k
Ik Ik - Ib
23
4/10/2013
Ij Ij + Ib
j Orig. network
with buses j, k
Ib Zb and reference
bus extracted
k
Ik Ik - Ib
Zbus(new) 0
V1 Z1j - Z1k I1
: : :
Vj Zbus(orig) Zjj - Zjk Ij
Vk Zkj - Zkk Ik
: = : :
Vn Znj - Znk In
0 (Zj1 - Zk1) ⋅ ⋅ ⋅ (Zkj - Zkk) ⋅ ⋅ ⋅ Zbb Ib
V1 Z1j - Z1k I1
: : :
Vj Zbus(orig) Zjj - Zjk Ij
Vk Zkj - Zkk Ik
: = : :
Vn Znj - Znk In
0 (Zj1 - Zk1) ⋅ ⋅ ⋅ (Zkj - Zkk) ⋅ ⋅ ⋅ Zbb Ib
Zbus(new)
24
4/10/2013
25
4/10/2013
1 2
1 j1.25 0
Zbus = 2
0 j1.25
1 2 3
1 j1.25 0 j1.25
Zbus = 2 0 j1.25 0
3 j1.25 0 j1.50
1 2 3
1 j1.25 0 j1.25 j1.25
Zbus = 2 0 j1.25 0 0
3 j1.25 0 j1.50 j1.50
j1.25 0 j1.50 j2.75
1 2 3
1 j0.68182 0 j0.56818
Zbus = 2 0 j1.25 0
3 j0.56818 0 j0.68182
26
4/10/2013
1 2 3
1 j0.68182 0 j0.56818 -j0.56818
2 0 j1.25 0 j1.25
Zbus = 3 j0.56818 0 j0.68182 -j0.68182
-j0.56818 j1.25 -j0.68182 j2.33182
1 2 3
1 j0.54338 j0.30458 j0.40205
Zbus = 2 j0.30458 j0.57992 j0.36550
3 j0.40205 j0.36550 j0.48246
1 2 3 4
1 j0.54338 j0.30458 j0.40205 j0.30458
2 j0.30458 j0.57992 j0.36550 j0.57992
Zbus =
3 j0.40205 j0.36550 j0.48246 j0.36550
4 j0.30458 j0.57992 j0.36550 j0.77992
27
4/10/2013
1 2 3 4
1 j0.54338 j0.30458 j0.40205 j0.30458 j0.23880
2 j0.30458 j0.57992 j0.36550 j0.57992 -j0.27534
Zbus = 3 j0.40205 j0.36550 j0.48246 j0.36550 j0.03655
4 j0.30458 j0.57992 j0.36550 j0.77992 -j0.47534
j0.23880 -j0.27534 j0.03655 -j0.47534 j0.91414
1 2 3 4
1 j0.48100 j0.37651 j0.39250 j0.42875
2 j0.37651 j0.49699 j0.37651 j0.43675
Zbus =
3 j0.39250 j0.37651 j0.48100 j0.38451
4 j0.42875 j0.43675 j0.38451 j0.53275
28
4/10/2013
1 2 3 4
1 j0.48100 j0.37651 j0.39250 j0.42875 -j0.03625
2 j0.37651 j0.49699 j0.37651 j0.43675 -j0.06024
Zbus =
3 j0.39250 j0.37651 j0.48100 j0.38451 j0.09649
4 j0.42875 j0.43675 j0.38451 j0.53275 -j0.14824
-j0.03625 -j0.06024 j0.09649 -j0.14824 j0.36973
1 2 3 4
1 j0.48100 j0.37651 j0.39250 j0.42875 -j0.03625
2 j0.37651 j0.49699 j0.37651 j0.43675 -j0.06024
Zbus =
3 j0.39250 j0.37651 j0.48100 j0.38451 j0.09649
4 j0.42875 j0.43675 j0.38451 j0.53275 -j0.14824
-j0.03625 -j0.06024 j0.09649 -j0.14824 j0.36973
1 2 3 4
1 j0.47745 j0.37060 j0.40196 j0.41422
2 j0.37060 j0.48718 j0.39223 j0.41260
Zbus =
3 j0.40196 j0.39223 j0.45582 j0.42320
4 j0.41422 j0.41260 j0.42320 j0.47331
29