An Optimal Power Flow Control Method of Power System PDF
An Optimal Power Flow Control Method of Power System PDF
An Optimal Power Flow Control Method of Power System PDF
S. Teerathana and A. Yokoyama are with the Electrical Engineering Static synchronous series
Department, the University of Tokyo, 7-3-1 Hongo Bunkyo-ku Tokyo 113- compensator (SSSC)
employed to increase the transferable active power on a given 360 degree range, the end of voltage V pq moves along a circle
•
line and to balance the loading of a transmission network. with its center located at the end of voltage Vs . The area within
•
In addition, active power can be exchanged through these this circle defines the operating range of voltage V pq .
•
two series converters via the common DC link in IPFC. It is According to the equivalent circuit of the IPFC shown in
noted that the sum of the active powers outputted from VSCs Fig. 3, the power flow equations can be obtained as follows
to transmission lines should be zero when the losses of the [5]:
converter circuits can be ignored. A combination of the series
n
connected VSC can inject a voltage with controllable Pi = Vi2 g ii − ∑ Vi V j ( g ij cos( θ j − θ i ) + bij sin( θ j − θ i ))
magnitude and phase angle at the fundamental frequency j =1, j ≠ i
(1)
while DC link voltage can be maintained at a desired level. n
− ∑ ViVseij ( g ij cos( θ i − θseij ) + bij sin( θ i − θseij ))
The common dc link is represented by a bidirectional link for j =1, j ≠ i
active power exchange between voltage sources.
As is illustrated in Fig.2, a phasor diagram of voltage, for n
instance, controlled by one converter of IPFC defines the Qi = −Vi2 bii − ∑ ViV j ( g ij sin( θ j − θ i ) − bij cos( θ j − θ i ))
j =1, j ≠ i
relationship between sending end-voltage ( Vs ), receiving end-
•
(2)
n
voltage ( Vr ), the voltage across line impedance X ( V x ) and the
• •
− ∑ ViVseij ( g ij sin( θ i − θseij ) − bij cos( θ i − θseij ))
• j =1, j ≠i
inserted voltage V pq , with controllable magnitude and phase
•
angle. When V pq is added to the sending-end voltage, an
effective sending-end voltage would be received as Vseff = Vs
• • P ji = V j2 g jj − Vi V j ( g ij cos( θ j − θ i ) + bij sin( θ j − θ i ))
(3)
• • •
+ V pq .So the difference Vseff – Vr , sets the compensated + V j Vseij ( g ij cos( θ j − θseij ) + bij sin( θ i − θse ij ))
voltage or V x across reactance X. As ρ is varied over its full
•
voltage source
leakage V min ≤ ≤ max
seij V seij V seij
(7)
converter
reactance jxline
θ min ≤ ≤ max
seij θ seij θ seij
(8)
Fig. 2. Phasor diagram of voltage control
where
Z seij V se ij g in + jbin = 1 / Zsein , g nn + jbnn = 1 / Zse nn (9)
V i V
+ − j
g ii = ∑ g in , bii = ∑ bin (10)
I ji n n
P ji + jQ ji and
P i + jQi Pji : active power flow from node j to node i
Re Vseij I ji + Vseik I ki = 0 θse : angle of injected voltage
Vse : magnitude of injected voltage of FACTS controllers
Vk θ : bus angle
+ − V : bus voltage magnitude
I ki
Z seik V seik Pki + jQki
III. NONLINEAR OPTIMAL POWER FLOW CONTROL IV. NUMERICAL EXAMPLES IN MULTI-MACHINE
The optimal power flow (OPF) is a very important POWER SYSTEM
mathematical programming problem in the power engineering In this paper, the proposed Optimal Power Flow (OPF)
field. In this paper, we propose the concept of OPF control control method is applied to 6-machine 22-bus test system as
method for coping with an overload constraint in a shown in Fig.4 for coping with overloaded problems. For
transmission network under IPFC operation. The entire set of simplicity, IPFC with two converters applied to the test system
power flow equations are solved simultaneously with are focused and its minimum capacity at some critical
minimization of the total capacity of the installed IPFC branches is discussed below.
required for solving the overloaded of the transmission lines. In the simulation, leakage reactance of both two converters
We discuss a method for solving a nonlinear OPF problem is assigned at 0.01 in the overloaded power system. As a high
that can be formulated as: demand load, L5 is examined in different three conditions at
6.37[GW], 6.41[GW] and 6.46[GW]. Consequently, G5 is
Objective function: needed to supply active power of 7.88[GW], 7.92[GW] and
7.97[GW] respectively. Therefore, this leads to the overload
minimize f ( x ) = PQ1 2 + PQ2 2 + (11) problem at branch 21 between the heavy load bus 21 serving
L5 and bus 20 that the active power flow through the branch
s.t. becomes more than its thermal constraint 2.784[GW] per line
hi ( x ) = 0 (12) when its one circuit is tripped.
gi ( x ) ≤ 0 (13) When the test system is forced that the active power flow
over branch 21 is more than its thermal constraint 2.784 [GW]
where per line by tripping one circuit simultaneously increasing load
PQ : capacity of each VSCs of IPFC demand at bus 5, the overload occurs.
h(x) : equality constraints such as bus power flow equations Case 1: Overload occurs at branch no.21 due to its one circuit
under the IPFC operating and control. trip. Installation of IPFC for a satisfactory solution at the
g(x) : inequality constraints such as thermal line constraints minimum capacity of two converters (PQ1 & PQ2) depends on
and voltage limitation. load demand and line current at each branch. Table1 showed
x : (e, f) the simulation results by OPF method for the minimum
e : the real part of bus voltage capacity of IPFC installed at branch no. 15&16, 16&17,
f : the imaginary part of bus voltage 21&22, 17&21 and 17&22 respectively.
Z ij Z ik [15]
[19] [17]
G1
(14)
[20] [21] [22] [23]
0.27[GW] 0.27[GW]
16 16
0.34[GW] → -0.35[GW] 2.80[GW] → 2.82[GW] 0.34[GW] → 0.10[GW]
2.80[GW] → 2.83[GW]
L5 =6.46[GW]
L5 =6.46[GW]
Fig. 9. Active Power flow controlled by IPFC at branch no.22 & 17 for
Fig. 5. Active Power flow controlled by IPFC at branch no.15 & 16 to
overload at branch no.21
overload at branch no.21
0.27[GW]
0.27[GW]
2.79[GW] → 2.80[GW] 16 0.40[GW] → 0.20[GW]
16
2.80[GW] → 2.83[GW] 0.34[GW] → 0.26[GW]
[16] [15]
[16] [15]
3.41[GW] → 3.27[GW]
Pex=0.01081[GW] 3.41[GW] → 3.36[GW] [17]
[17]
2.81[GW] → 2.784[GW] Pex=0.0089[GW] 0.19[GW] → 0.41[GW]
2.805[GW] → 2.784[GW] 0.24[GW] → 0.32[GW]
[21] [22]
[21] [22] 21
21
L5 =6.46[GW]
L5 =6.46[GW]
Fig. 10. Active Power flow controlled by IPFC at branch no.21 & 22 when
Fig. 6. Active Power flow controlled by IPFC at branch no.16 & 17 for one circuit trip of branch 22 during overload at branch 21
overload at branch no.21
0.27[GW]
0.27[GW]
16
16 2.79[GW] → 2.81[GW] 0.40[GW] → 0.27[GW]
2.79[GW] → 2.82[GW] 0.34[GW] → 0.099[GW]
[16] [15]
[16] [15]
Pex=0.00625[GW] 3.41[GW] → 3.35[GW]
3.41[GW] → 3.18[GW]
[17]
[17] 2.805[GW] → 2.784[GW] 0.19[GW] → 0.33[GW]
2.805[GW] → 2.784[GW] Pex=0.00705[GW] 0.24[GW] → 0.50[GW]
[21] [22]
[21] [22] 21
21
L5 =6.46[GW]
L5 =6.46[GW]
Fig. 11. Active Power flow controlled by IPFC at branch no.21 & 17 when
Fig. 7. Active Power flow controlled by IPFC at branch no.21 & 22 for one circuit trip of branch 22 during overload at branch 21
overload at branch no.21
As is shown in Table2, to solve the problem of overload at
branch 21 by using the minimum capacity of IPFC, it would
0.27[GW]
be a great way to install one converter of IPFC at the branch
16
2.80[GW] → 2.82[GW] 0.34[GW] → 0.17[GW] 21, the overloaded line. As for the other converter installation,
[15] the total required capacity of IPFC does not make a big
[16]
3.41[GW] → 3.25[GW]
difference between installations of branches from the same
Pex=0.00473[GW] [17] bus 21, branches 17 and 22.Figs.5 to 9 illustrates only active
2.805[GW] → 2.784[GW] 0.24[GW] → 0.43[GW] power flows when IPFC is employed to solve the overload
[21] [22] problem at branch no.21. It can be seen from Figs.5 to 9 that
21 each of converters is adjusted at a proper magnitude and phase
angle of the inserted voltage so that the sum of active powers
L5 =6.46[GW] exchanged between two series converters via the common
DC link becomes zero if we can neglect the losses of the
Fig. 8. Active Power flow controlled by IPFC at branch no.21 & 17 for converter circuits.
overload at branch no.21
5
As is shown in Fig. 5, to solve the overload problem so that overloaded line should be the first step to solve the problem
active power flow in branch 21 become less than 2.784[GW], before installation of the other converter at another lines from
the IPFC with two converters installed at branches no.15 and the same bus. Also, Table 2 shows that one converter should
16 respectively provides a capability to change the direction of be installed at branch 21 and then the other one may be at
power flow over branch 15 into the inverse direction. branch 17or branch22 both of which are underloaded.
Case2: The same system model as in case 1 is still used here. Although the total required capacity in case of two converters
One circuit of branch 22 in the other loop is forced to be installed at branch 17 and branch 21 is smaller than at branch
tripped one circuit during overload line at branch 21 caused by 21 and branch 22, the exchanged power between two series
its one circuit trip. By the proposed OPF method, the converters in branch 21 and branch 22 becomes larger at every
minimum capacities of IPFC inserted at various branches are load demand at load 5. That the power flow on branch 22 in
obtained as shown in Table3. normal condition is less than that on branch 17 offers more
Figures 10 and 11 show the active power flow under IPFC transferable power flow on branch 22. In addition to the
operation at branch 21&22 and branch 17&21 and the condition that the power flow goes worse, for instance, when
transferred power between two converters. one circuit trips, the total required capacity becomes larger as
It can be seen that to locate one converter of IPFC at illustrated in Tables 2 and 3.
Table II: Capacity of IPFC in case of overload at branch 21 Table III: Capacity of IPFC in case of overload at both branches 21 & 22