CPT3 - Load Flow Analysis-6th Batch
CPT3 - Load Flow Analysis-6th Batch
Certificate in
Course Outline
1. The Load Flow Problem 2. Power System Models for Load Flow Analysis 3. Gauss-Seidel Load Flow 4. Newton-Raphson Load Flow 5. Backward/Forward Sweep Load Flow 6. Principles of Load Flow Control 7. Uses of Load Flow Studies
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
Basic Electrical Engineering Solution Load Flow of Distribution System Load Flow of Transmission and Subtransmission System Load Flow of a Contemplated System Load Flow of a Single Line
Competency Training & Certification Program in Electric Power Distribution System Engineering
Solve for:
1) ISR = (SR/VR )* 2) VD = ISRZL 3) VS = VR + VD 4) SS = VSx(ISR)*
VS = ?
VOLTAGE DROP = VS - VR
Competency Training & Certification Program in Electric Power Distribution System Engineering
Receiving End
Solve for:
1) ISR = (SR/VR )* 2) VD = ISRZL 3) VS = VR + VD 4) SS = VSx(ISR)*
VS = ?
VD = ( 87.48 31.79 )( 1.1034 + j2.0856 ) = 178.15 + j104.23 V VS = (7621.02 + j0 ) + ( 178.15 + j104.23 ) = 7,799.87 0.77 V VS = 7,799.87 0.77 /1000* 3 = 13.51 k V
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
VS = 13.2 kVLL
Competency Training & Certification Program in Electric Power Distribution System Engineering
Utility Grid
V1 = 67 kV P1 , Q1 = ?
How do you solve for the Voltages, Currents, Power and Losses?
I24 , Loss24 = ? V4 = ? P4 , Q4 = ? Bus4 V2 = ? Lumped Load A P2 , Q2 = ? 2 MVA 85%PF Lumped Load B 1 MVA 85%PF
Competency Training & Certification Program in Electric Power Distribution System Engineering
Competency Training & Certification Program in Electric Power Distribution System Engineering
10
Flow
Analysis
Competency Training & Certification Program in Electric Power Distribution System Engineering
11
VS = 13.2 kVLL
VR = ?
12
VS = 13.2 kVLL
VR = ?
Converting Quantities in Per Unit Base Power = 1 MVA Base Voltage = 13.2 kV Base Impedance = [13.2]2/1 = 174.24 ohms
U. P. National Engineering Center National Electrification Administration
VS(pu) = 13.2 /13.2 = 1/0 SR(pu) = 2/cos-1(0.85) / 1 Zpu = (1.1034 + j2.0856)/174.24 = 0.00633 + j0.01197
Competency Training & Certification Program in Electric Power Distribution System Engineering
13
VS = 13.2 kVLL
VR = ?
VR(0) = 1/0
For k = 2
For k = 1
14
VS = 13.2 kVLL
VR = ?
15
VS = 13.2 kVLL
VR = ?
16
VS = 13.2 kVLL
VR = ?
17
Bus Admittance Matrix, Ybus Network Models Generator Models Bus Types for Load Flow Analysis
Competency Training & Certification Program in Electric Power Distribution System Engineering
18
Competency Training & Certification Program in Electric Power Distribution System Engineering
19
[YBUS] =
The number of buses (excluding the neutral bus) determines the dimension of the bus admittance, Ybus.
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
20
Line 1 1 Line 2 3
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
21
22
23
24
25
voltage magnitude, Vp voltage phase angle, p injected active power, Pp injected reactive power, Qp
Competency Training & Certification Program in Electric Power Distribution System Engineering
26
Specify: V, Unknown: P, Q
Competency Training & Certification Program in Electric Power Distribution System Engineering
27
Specify: P, V Unknown: Q,
Competency Training & Certification Program in Electric Power Distribution System Engineering
28
V
-
29
Competency Training & Certification Program in Electric Power Distribution System Engineering
30
Bus Types
Line 3
Competency Training & Certification Program in Electric Power Distribution System Engineering
31
Linear Formulation of Load Flow Equations Gauss-Seidel Load Flow Solution Numerical Example
Competency Training & Certification Program in Electric Power Distribution System Engineering
32
The real and reactive power into any bus P is: or where
(1)
33
(2)
(3)
34
(4)
35
Y22V2 =
36
(5)
Competency Training & Certification Program in Electric Power Distribution System Engineering
37
(6)
where, = k =k+1
U. P. National Engineering Center National Electrification Administration
if p < q if p > q
Competency Training & Certification Program in Electric Power Distribution System Engineering
38
Gauss-Seidel Voltage Equations of the form shown in (6) are written for all buses except for the swing bus. The solution proceeds iteratively from an estimate of all bus voltages For a Load Bus (Type 3) whose real power and reactive power are specified, the G-S voltage equation is used directly to compute the next estimate of the bus voltage. For a Generator Bus (Type 2) where the voltage magnitude is specified, an estimate of Qp must be determined first. This estimate is then compared with the reactive power limits of the generator. If it falls within the limits, the specified voltage is maintained and the computed Qp is inputted, in the Gauss-Seidel equation. Otherwise, the reactive power is set to an appropriate limit (Qmin or Qmax) and the bus is treated as a load bus in the current iteration.
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
39
Line 1
1 Line 2 3
U. P. National Engineering Center National Electrification Administration
2 Line 3
Competency Training & Certification Program in Electric Power Distribution System Engineering
40
Bus Data Bus Voltage Generation No. V (p.u.) P Q 1 1.0 0.0 * * 2 1.0 * 0.20 * 3 * * 0 0
U. P. National Engineering Center National Electrification Administration
Competency Training & Certification Program in Electric Power Distribution System Engineering
41
42
43
= 1 0
Bus 2: Generator Bus Q2 must first be determined from: P2 - jQ2(k+1) = (V2(k))* [Y21V1(k+1) + Y22V2(k) + Y23V3(k)] then substitute it to:
V2
(k + 1 )
(k + 1 ) 1 P2 jQ 2 (k + 1 ) (k ) = Y 21 V 1 Y 23 V 3 * ( ) k Y 22 V2
Competency Training & Certification Program in Electric Power Distribution System Engineering
44
V3
1 P3 jQ 3 (k + 1 ) (k + 1 ) = Y 31 V 1 Y 32 V 2 * Y 33 V (k ) 3
Iteration 1 (k = 0): V1 (1) = 1.00 P2 - jQ2(1) = (1.00) [(-1.25 + j3.75)(1.00) + (2.9167 - j8.75)(1.00) + (-1.6667 + j5)(1.00) = 0.0 + j0.0 Q2(1) = 0.0 [This value is within the limits.]
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
45
(k + 1 ) 1 P2 jQ 2 (k + 1 ) (k ) = Y 21V 1 Y 23 V 3 * (k ) Y 22 V 2
V2
(1)
P2 jQ2
( k +1 )
Y22
(V )
(k )
2
= 1.00711.1705
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
46
1 P3 jQ 3 (k + 1 ) (k + 1 ) = Y 31V 1 Y 32 V 2 * ( ) k Y 33 V 3 P jQ 3 3
V31
(V )
3
k *
( k +1 ) Y32 V2 - (5.2705108.4349)(1.00711.1705)
= 0.9816 -1.0570
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
47
48
49
(k + 1 ) 1 P2 jQ 2 (k + 1 ) (k ) = Y 21V 1 Y 23 V 3 * k ( ) Y 22 V2
V2
(2)
50
1 P3 jQ 3 (k + 1 ) (k + 1 ) = Y 31 V 1 Y 32 V 2 * Y 33 V (k ) 3
51
52
53
(k + 1 ) 1 P2 jQ 2 (k + 1 ) (k ) = Y 21V 1 Y 23 V 3 * k ( ) Y 22 V2
V23
Competency Training & Certification Program in Electric Power Distribution System Engineering
54
1 P3 jQ 3 (k + 1 ) (k + 1 ) = Y 31 V 1 Y 32 V 2 * Y 33 V (k ) 3
- (5.2705108.4349)(0.99900.4129) = 0.9788-1.2560
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
55
56
Non-Linear Formulation of Load Flow Equations Newton-Raphson Load Flow Solution Numerical Example
Competency Training & Certification Program in Electric Power Distribution System Engineering
57
Pp jQ p = E* pIp
and the current equation may be written as
Ip =
(1)
Y
q =1
pq
Eq
n
(2)
P p jQ
= E
* p
Y E
q =1 pq
(3)
q
Competency Training & Certification Program in Electric Power Distribution System Engineering
58
E q = Vq q
Ypq = Ypq pq
Substituting into equation (3),
Pp jQ p =
Pp =
n
VVY
q =1 p q
q pq
pq
( pq + q p )
(4)
VVY
q =1 p
n q =1
co s( pq + q p )
(5) (6)
Q p = V pV q Y pq sin( pq + q p )
U. P. National Engineering Center National Electrification Administration
Competency Training & Certification Program in Electric Power Distribution System Engineering
59
Competency Training & Certification Program in Electric Power Distribution System Engineering
60
(7)
61
Pi spec
M M M
= Pi calc = Q1calc
calc = Q2
+ + +
M P
i
1 +
M P
i
+ + +
+ + +
Q1spec Q2spec
Q1 1 1 Q2 1 1
+ +
Q1 2 2 Q2 2 2
Q1 i i Q2 i i Q j i
Q1 V1 V1 Q2 V1 V1 Q j V1
Q1 V2 V2 Q2 V2 V2 Q j V2
M M M
spec j
= Q
calc j
M Q
1 +
M Q
+... +
V1 +
V2 +... +
V j
Competency Training & Certification Program in Electric Power Distribution System Engineering
62
P1 2 P2 2
L L
P1 i P2 i
P1 V1 P2 V1
P1 V2
P2 V2
L
L
M
Pi 2 Q1 2 Q2 2
M L L L M
j
M
Pi V1 Q1 V1 Q 2 V1
M
Pi V2 Q1 V2 Q 2 V2
Pi 2 Q1 i Q 2 i
L L L
M
Q
M
Q i
j
M
Q
j
M
Q
j
V1
V2
P1 1 V j P2 2 V j M M Pi i V j Q1 V1 V j Q2 V 2 V j M M Q j V j V j
Competency Training & Certification Program in Electric Power Distribution System Engineering
63
P P Q = Q
P V Q V V
P P V P V V Q = Q Q V V V
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
64
n Pp = V p V q Y p q s in ( p q + q q = 1 ,q p p J1 P p = V V Y s in ( + ) p q pq pq q p q
U. P. National Engineering Center National Electrification Administration
Competency Training & Certification Program in Electric Power Distribution System Engineering
65
J2
J3
Competency Training & Certification Program in Electric Power Distribution System Engineering
66
The solution of the load flow equations proceeds iteratively from the set of initial estimates. These estimates are updated after evaluating the Jacobian matrix.
Competency Training & Certification Program in Electric Power Distribution System Engineering
67
V p( k +1 ) = V p( k ) + V p( k )
The process is terminated once convergence is achieved whrein
MAX P( k ) p
and
MAX Q ( k ) q
Competency Training & Certification Program in Electric Power Distribution System Engineering
68
G 1 Line 2
Line 1
G 2 Line 3 3
Competency Training & Certification Program in Electric Power Distribution System Engineering
69
Bus Data Bus Voltage Generation No. V (p.u.) P Q 1 1.0 0.0 * * 2 1.0 * 0.20 * 3 * * 0 0
U. P. National Engineering Center National Electrification Administration
Competency Training & Certification Program in Electric Power Distribution System Engineering
70
Y11 = 6.25 j18.75 = 19.7642 71.5651 Y12 = 1, 25 + j 3.75 = 3.9528 108.4349 Y13 = 5 + j15 = 15.8114 108.4349 Y 21 = 1.25 + j .375 = 3.9528 108.4349 Y 22 = 2.9167 j 8.75 = 9.2233 71.5649 Y 23 = 1.6667 + j 5 = 5.2705 108.4349 Y31 = 5 + j15 = 15.811 4 108.4349 Y32 = 1.6667 + j 5 = 5.2705 108.4349 Y33 = 6.6667 j 20 = 21.0819 71.5650
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
71
P2 P2 2 P3 3= P 2 Q3 Q3 2
P2 3 P 3 3 Q3 3
P2 V3 2 V3 P V3 3 3 V3 Q3 V3 V3 V3 V3
Competency Training & Certification Program in Electric Power Distribution System Engineering
72
73
P20 = V2V1Y21 cos( 21 + 1 2 ) + V2V2Y22 cos 22 + V2V3Y23 cos( 23 + 3 2 ) = ( 1.0 )( 1.0 )( 3.9528 )cos( 108.4349 + 0.0 + 0.0 ) + ( 1.0 )( 1.0 )( 9.2233 )cos( 71.5649 ) + ( 1.0 )( 1.0 )( 5.2705 )cos( 108.4349 ) = 0.0
Competency Training & Certification Program in Electric Power Distribution System Engineering
74
Competency Training & Certification Program in Electric Power Distribution System Engineering
75
+ V 3V 2Y32 sin( 32 + 2 3 ) + V 3V 3Y33 sin 33 = ( 1.0 )( 1.0 )( 15.8114 ) sin( 108.4349 + 0.0 0.0 ) + ( 1.0 )( 1.0 )( 5.2705 ) sin( 108.4349 + 0.0 + 0.0 ) + ( 1.0 )( 1.0 )( 21.0819 ) sin( 71.5650 ) = 0.0
Competency Training & Certification Program in Electric Power Distribution System Engineering
76
P20 = 0 . 2 0 . 0 = 0 . 2 P30 = 0 .6 0 .0 = 0 .6 Q 30 = 0 .2 5 0 .0 = 0 .2 5
Evaluate elements of Jacobian Matrix
P J = Q P Q V V V
Competency Training & Certification Program in Electric Power Distribution System Engineering
77
P2 = V2V1Y21 sin( 21 + 1 2 ) 2 + V2V3Y23 sin( 23 + 3 2 ) = ( 1.0 )( 1.0 )( 3.9528 )sin( 108.4349 + 0.0 0.0 ) + ( 1.0 )( 1.0 )( 5.2705 )sin( 108.4349 + 0.0 0.0 ) = 8.75
Competency Training & Certification Program in Electric Power Distribution System Engineering
78
79
P3 = V3V1Y31 sin( 31 + 1 3 ) 3 + V3V2Y23 sin( 32 + 2 3 ) = ( 1.0 )( 1.0 )( 15.8114 )sin( 108.4349 + 0.0 0.0 ) + ( 1.0 )( 1.0 )( 5.2705 )sin( 108.4349 + 0.0 0.0 ) = 20
Competency Training & Certification Program in Electric Power Distribution System Engineering
80
P2 = V2V3Y23 cos( 23 + 3 + 2 ) V3 V3
81
Q3 = V3V2Y32 cos( 32 + 2 3 ) 2 = ( 1.0 )( 1.0 )( 5.2705 )cos( 108.4349 + 0.0 0.0 ) = 1.6667
+ V3Y2Y32 cos( 32 + 2 3 ) = ( 1.0 )( 1.0 )( 15.8114 )cos( 108.4349 + 0.0 0.0 ) + ( 1.0 )( 1.0 )( 5.2705 )cos( 108.4349 + 0.0 0.0 ) = 6.6667
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
Q3 = V3V1Y31 cos( 31 + 1 3 ) 3
82
8.75 5 1.6667
83
84
21 = 20 + 20
21 = 0.0 + 0.2283 = 0.2283 31 = 30 + 30 31 = 0.0 1.4822 = 1.4822
V = V + V
1 3 0 3 0 3
Specified Variables
V11 = 1.0
11 = 0.0
V21 = 1.0
Competency Training & Certification Program in Electric Power Distribution System Engineering
85
+ V2V2Y22 cos(22 )
+ 1.4822 = (1.0 )(1.0 )(3.9528 )cos(108.4349 + 0.0 0.2283) + (1.0 )(1.0 )(9.2233)cos( 71.5649 ) + (1.0 )(0.97855)(5.2705)cos(108.4349 1.4822 0.2283) = 0.1975
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
+ V2V3Y23 cos(23 + 3 2 )
86
P31 = V3V1Y31 cos( 31 + 1 3 ) + V3V2Y32 cos( 32 + 2 3 ) + V3V3Y33 cos( 33 ) = (0.97855 )(1.0 )(15.8114 )cos(108.4349 + 0.0 + 1.4822 ) + (0.97855 )(1.0 )( 5.2705 )cos(108.4349 + 0.2283 + 1.4822 ) + (0.97855 )(0.97855 )( 21.0819 )cos( 71.5650 ) = 0.66633
Competency Training & Certification Program in Electric Power Distribution System Engineering
87
+V3V3Y33 sin( 33 )]
+ V3V2Y32 sin( 32 + 2 3 )
= - [(0.97855 )(1.0 )(15.8114 )sin(108.4349 + 0.0 + 1.4822 ) + (0.97855 )(1.0 )( 5.2705 )sin(108.4349 + 0.2283 + 1.4822 ) +(0.97855 )(0.97855 )( 21.0819 )sin( 71.5650 )] = 0.2375
Competency Training & Certification Program in Electric Power Distribution System Engineering
88
P3 = P3 ,s p P3 ,c a lc
= 0 .6 + 0 .6 6 3 3 = 0 .0 6 3 3
Q 3 = Q 3 ,s p Q 3 ,c a lc
= .0 .2 5 + 0 .2 3 7 5 = 0 .0 1 2 5
U. P. National Engineering Center National Electrification Administration
Proceed to Iteration 2
Competency Training & Certification Program in Electric Power Distribution System Engineering
89
P 2 = VVY 2 1 21 sin( 21 + 1 2 ) 2 + VV 2 3Y 23 sin(23 + 3 2 ) = (1.0)(1.0)(3.9528)sin(108.4349 + 0.0 0.2283) + (1.0)(0.97855)(5.2705)sin(108.4349 1.4822 0.2283) = 8.6942
Competency Training & Certification Program in Electric Power Distribution System Engineering
90
Competency Training & Certification Program in Electric Power Distribution System Engineering
91
P3 = V3V1Y31 sin( 31 + 1 3 ) 3 + V3V2Y23 sin( 32 + 2 3 ) = (0.97855 )(1.0 )(15.8114 )sin(108.4349 + 0.0 + 1.4822 ) + (0.97855 )(1.0 )( 5.2705 )sin(108.4349 + 0.2283 + 1.4822 ) = 19.3887
Competency Training & Certification Program in Electric Power Distribution System Engineering
92
P2 V3 = V2V3Y23 cos( 23 + 3 + 2 ) V3
93
Q3 = V3V1Y31 cos( 31 + 1 3 ) 3
+ V3Y2Y32 cos( 32 + 2 3 ) = ( 0.97855 )( 1.0 )( 15.8114 )cos( 108.4349 + 0.0 + 1.4822 ) + ( 0.97855 )( 1.0 )( 5.2705 )cos( 108.4349 + 0.2283 + 1.4822 ) = 7.0470
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
94
95
0.0025 8.6942 4.9393 1.4842 2 0.0633 = 4.8419 19.3887 5.7205 3 0.0125 1.7762 7.0470 18.9137 V3 / V3
V30
V3
= 0.0005
Competency Training & Certification Program in Electric Power Distribution System Engineering
96
22 = 21 + 21
= 0.2283 + 0.1458 = 0.37410
32 = 31 + 31
= 1.4822 + 0.2150 = 1.2672 0 V32 = V31 + V31 = 0.97855 + 0.0005 = 0.9791
Specified Variables
V11 = 1.0
11 = 0.0
V21 = 1.0
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
97
P22 = V2V1Y21 cos( 21 + 1 2 ) + V2V2Y22 cos( 22 ) + V2V3Y23 cos( 23 + 3 2 ) = ( 1.0 )(1.0 )( 3.9528 )cos(108.4349 + 0.0 0.3741) + ( 1.0 )(1.0 )( 9.2233 )cos( 71.5649 ) + ( 1.0 )(0.9791)( 5.2705 )cos( 108.4349 1.2672 0.3741) = 0.2018
Competency Training & Certification Program in Electric Power Distribution System Engineering
98
P32 = V3V1Y31 cos( 31 + 1 3 ) + V3V2Y32 cos( 22 + 2 3 ) + V3V3Y33 cos 33 = ( 0.9791 )( 1.0 )( 15.8114 )cos( 108.4349 + 0.0 1.2672 ) + ( 0.9791 )( 1.0 )( 5.2705 )cos( 108.4349 + 0.3741 + 1.2672 ) + ( 0.9791 )( 0.9791 )( 21.0819 )cos( 71.5650 ) = 0.5995
Competency Training & Certification Program in Electric Power Distribution System Engineering
99
+V3V2Y32 sin( 32 + 2 3 ) + V3V3Y33 sin33 ] = [(0.9791)(1.0 )(15.8114 )sin(108.4349 + 0.0 1.2672 ) + (0.9791)(1.0 )( 5.2705 )sin(108.4349 + 0.37411 + 1.2672 ) +(0.9791)(0.9791)( 21.0819 )sin( 71.5650 )] = 0.2487
Competency Training & Certification Program in Electric Power Distribution System Engineering
100
P3 = P3 ,s p P3 ,c a lc
= 0 .6 + 0 .5 9 9 5 = 0 .0 0 0 5
Q 3 = Q 3 ,s p Q 3 ,c a lc
= .0 .2 5 + 0 .2 4 8 7 = 0 .0 0 1 3
U. P. National Engineering Center National Electrification Administration
101
V1 = 1.00 0 V2 = 1.00.3741
0 0
V3 = 0.9791 1.2672
Competency Training & Certification Program in Electric Power Distribution System Engineering
102
The power injected into the buses are: P1 - jQ1 = V1* [Y11V1 + Y12V2 + Y13V3 ] P1 - jQ1 = (1.00) [(19.7642-71.5651)(1.00) + (3.9528108.4349)(0.99900.4129) + (15.8114108.4349) (0.9788-1.25560) = 0.4033 - j0.2272
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
103
P2 - jQ2 = (0.999-0.4129)[(3.9528108.4349)(1.00) + (9.2233-71.5649)(0.99900.4129) + (5.2705108.4349)(0.9788-1.25560) = 0.2025 - j0.04286 P3 - jQ3 = V3* [Y31V1 + Y32V2 + Y33V3 ] P3 - jQ3 = (0.97881.256) [(15.8114108.4349)(1.00) + (5.2705108.4349)(0.99900.4129) + (21.0819 -71.5650)(0.9788-1.25560 ) = -0.600 + j0.2498
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
104
I pq = I line = y pq ( V p Vq )
I12 = y12 [V1 - V2]
I qp = I line = y pq ( Vq V p )
I21 = y12 [V2 V1]
105
Iline ypo
ypq
Vq
q
Iqp
yqo
I pq = I line + I po = y pq ( V p Vq ) + y poV p
Similarly, the line current Iqp, measured at bus q is I = I + I = y (V V ) + y V
qp line qo pq q p qo q
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
106
Competency Training & Certification Program in Electric Power Distribution System Engineering
107
S pq = Ppq jQ pq = V I
* p pq
Competency Training & Certification Program in Electric Power Distribution System Engineering
108
Competency Training & Certification Program in Electric Power Distribution System Engineering
109
= (V p + V q )I
U. P. National Engineering Center National Electrification Administration
= VpI
* pq
VqI
* pq
* pq
Competency Training & Certification Program in Electric Power Distribution System Engineering
110
Competency Training & Certification Program in Electric Power Distribution System Engineering
111
Competency Training & Certification Program in Electric Power Distribution System Engineering
112
Load Flow for Radial Distribution System Procedure: Iterative Solution Initialization Solving for Injected Currents through the nodes Backward Sweep Forward Sweep Solving for Injected Power Solving for Voltage Mismatch
Competency Training & Certification Program in Electric Power Distribution System Engineering
113
Utility Grid
0.131 + j1.595
V1 = 67 kV P1 , Q1 = ?
I24 , Loss24 = ? V4 = ? P4 , Q4 = ? 0.4223 + j0.7980 Bus4 V2 = ? Lumped Load A P2 , Q2 = ? 2 MVA 85%PF Lumped Load B 1 MVA 85%PF
Bus3 V3 = ? P3 , Q3 = ?
Competency Training & Certification Program in Electric Power Distribution System Engineering
114
0.0364 +j 0.1131 pu
Utility Grid
Bus3 V3
V1
0.0242+j0.0458pu
Bus4 V4
1 + j0 pu
V2
Base Values
Sbase = 10 MVA Vbase1 = 67 kV Vbase2 = 13.2 kV
U. P. National Engineering Center National Electrification Administration
0.085 + j0.05267 pu
0.17 + j0.10536 pu
115
116
0.0364 +j 0.1131 pu
Utility Grid
Bus3 V3
V1
0.0242+j0.0458 pu
Bus4 V4
1 + j0 pu
V2
Initialize,
0.085 + j0.05267 pu
0.17 + j0.10536 pu
Competency Training & Certification Program in Electric Power Distribution System Engineering
117
0.0364 +j 0.1131 pu
Utility Grid
Bus3 V3
0.17 + j0.10536 pu
V1
0.0242+j0.0458 pu
Bus4 V4
0.085 + j0.05267 pu
1 + j0 pu
V2
I1(0) = 0 I2(0) = 0
Competency Training & Certification Program in Electric Power Distribution System Engineering
118
0.0364 +j 0.1131 pu
Utility Grid
Bus3 V3
0.17 + j0.10536 pu
V1
0.0242+j0.0458 pu
Bus4 V4
0.085 + j0.05267 pu
1 + j0 pu
V2
I24(0) = I4(0) = _______ I23(0) = I3(0) = _______ I12(0) = 0 + I23(0) + I24(0) = _______
Competency Training & Certification Program in Electric Power Distribution System Engineering
119
Utility Grid
V1
Bus3 V3
0.0075+j 0.0915 pu
0.0242+j0.0458 pu
Bus4 V4
1 + j0 pu
V2
V1(1) = 1/0
0.085 + j0.05267 pu
0.17 + j0.10536 pu
V2(1) = V1(0) [I12(0)][Z12] = ________ V3(1) = V2(1) [I23(0)][Z23] = ________ V4(1) = V2(1) [I24(0)][Z24] = ________
Competency Training & Certification Program in Electric Power Distribution System Engineering
120
0.0364 +j 0.1131 pu
Utility Grid
V1
0.0242+j0.0458 pu
Bus4 V4
1 + j0 pu
V2
S1(1) = [V1(1)][I1(0)]* = ___________ S2(1) = [V2(1)][I2(0)]* = ___________ S3(1) = [V3(1)][I3(0)]* = ___________ S4(1) = [V4(1)][I4(0)]* = ___________
0.085+ j0.05267 pu
Competency Training & Certification Program in Electric Power Distribution System Engineering
121
0.0364 +j 0.1131 pu
Utility Grid
Bus3 V3
V1
0.0242+j0.0458 pu
Bus4 V4
1 + j0 pu
V2
Solve S1(1) = S1(sp) - S1(calc) = ____________ Power Mismatch S2(1) = S2(sp) S2(calc) = ____________ S3(1) = S3(sp) S3(calc) = ____________ S4(1) = S4(sp) S4(calc) = ____________
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
0.085 + j0.05267 pu
0.17 + j0.10536 pu
122
123
124
Competency Training & Certification Program in Electric Power Distribution System Engineering
125
126
127
Competency Training & Certification Program in Electric Power Distribution System Engineering
128
Bus3
Utility Grid
0.0075+j 0.0915 pu
0.0242+j0.0458
Bus4
INJECTED POWER P1 + jQ1 = ________ + j ________ P2 + jQ2 = ________ + j ________ P3 + jQ3 = ________ + j ________ P4 + jQ4 = ________ + j ________
Competency Training & Certification Program in Electric Power Distribution System Engineering
129
Bus3
Utility Grid
0.0075+j 0.0915 pu
0.0242+j0.0458 pu
P23 + jQ23 = ________ + j ________ 1 MVA 85%PF P24 + jQ24 = ________ + j ________ POWER FLOW (Q-P)
P21 + jQ21 = ________ + j ________ P32 + jQ32 = ________ + j ________ P42 + jQ42 = ________ + j ________
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
130
Bus3
Utility Grid
0.0075+j 0.0915 pu
0.0242+j0.0458 pu
Branch Currents I12 = ________ I23 = ________ I24 = ________ POWER LOSSES
I2R12 + jI2X12 = ________ + j ________ I2R23 + jI2X24 = ________ + j ________ I2R24 + jI2X24 = ________ + j ________
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
131
19
132
I ia I ib I ic
(k )
S / V ( k 1 ) ia ia ( k 1 ) = Sib / Vib S / V ( k 1 ) ic ic
( ( (
) ) )
Where,
I ia , I ib , I ic V ia ,V ib ,V ic Y ia ,Y ib ,Y ic
Current injections at node i Voltages at node i Admittances of all shunt elements at node i
Competency Training & Certification Program in Electric Power Distribution System Engineering
133
J la J lb J lc
Where,
(k )
I ja = I jb I jc
(k )
J ma J + mb mM J mc
(k )
J la , J lb , J jc
l and M
are the current flows on line section l Is the set of line sections connected to node j
Competency Training & Certification Program in Electric Power Distribution System Engineering
134
V ja V jb V jc
(k )
Via V = ib Vic
(k )
(k )
Competency Training & Certification Program in Electric Power Distribution System Engineering
135
S S S
(k ) ia (k ) ib (k ) ic
=V =V =V
(k ) ia (k ) ib (k ) ic
(I ) (I ) (I )
(k ) ia (k ) ib (k ) ic
Y Via
* ia * ia * ic
S ia S ib S ic
Y Vib Y Vic
If the real and imaginary part (real and reactive power) of any of these power mismatches is greater than a convergence criterion, steps 1, 2 & 3 are repeated until convergence is achieved.
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
136
Prime mover and excitation control of generators Reactive Var Compensation (e.g., Capacitors) Control of tap-changing and voltage regulating transformers
Competency Training & Certification Program in Electric Power Distribution System Engineering
137
Vt0
EiVt Pt = sin X
U. P. National Engineering Center National Electrification Administration
EiVt Vt 2 Qt = cos X X
Competency Training & Certification Program in Electric Power Distribution System Engineering
138
EiVt = Pt sin X
Observations:
EiVt Vt 2 Qt = cos X X
1. Real Power is injected into the bus (Generator Operation), must be positive (Ei leads Vt) 2. Real Power is drawn from the bus (Motor Operation), must be negative (Ei lags Vt) 3. In actual operation, the numeric value of is small & since the slope of Sine function is maximum for small values, a minute change in can cause a substantial change in Pt
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
139
EiVt = Pt sin X
Observations:
EiVt Vt 2 Qt = cos X X
4. Reactive Power flow depends on relative values of EiCos and Vt 5. Since the slope of Cosine function is minimum for small values of angle, Reactive Power is controlled by varying Ei
Over-excitation (increasing Ei) will deliver Reactive Power into the Bus Under-excitation (decreasing Ei) will absorb Reactive Power from the Bus
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
140
Eq = Vp
X pqQq Vp
X pq Pq Vp
Observations:
1. The Reactive Power Qq causes a voltage drop and thus largely affects the magnitude of Eq 2. A capacitor bank connected to bus q will reduce Qq that will consequently reduce voltage drop
U. P. National Engineering Center National Electrification Administration Competency Training & Certification Program in Electric Power Distribution System Engineering
141
1 y pq a
Observation:
The voltage drop in the transformer is affected by the transformation ratio a
1 a y pq 2 a
a 1 y pq a
Competency Training & Certification Program in Electric Power Distribution System Engineering
142
Sensitivity Analysis with Load Flow Study Analysis of Existing Conditions Analysis for Correcting PQ Problems Expansion Planning Contingency Analysis System Loss Analysis
Competency Training & Certification Program in Electric Power Distribution System Engineering
143
144
Check for branch power flow violations Transformer Overloads Line Overloads Check for system losses Caps on Segregated DSL
Competency Training & Certification Program in Electric Power Distribution System Engineering
145
Competency Training & Certification Program in Electric Power Distribution System Engineering
146
Typical Taps
Tap 1: +5% Tap 2: +2.5% Tap 3: 0% (Rated Voltage) Tap 4: -2.5% Tap 5: -5%
Competency Training & Certification Program in Electric Power Distribution System Engineering
147
Competency Training & Certification Program in Electric Power Distribution System Engineering
148
Competency Training & Certification Program in Electric Power Distribution System Engineering
149
Competency Training & Certification Program in Electric Power Distribution System Engineering
150
Competency Training & Certification Program in Electric Power Distribution System Engineering