Sheth 2015
Sheth 2015
, ,
i j
Fig. 2. Typical one line diagram of transmission line
Utilizing the concept of power flow in the line and analyzing
with ‘π’ model representation, the power flow at the receiving
and sending end can be expressed as
Fig. 1. Mathematical Model of TCSC
| || | | |
(1) (4)
(2) | | | || |
(5)
(3)
From these power equations one can separate real and
Where, XL = Reactance of transmission line. reactive power,
XTCSC = Reactance of TCSC. V V V
P cos θ δ δ cosθ (6)
Z Z
rTCSC = Coefficient which represents the
V V V
compensation of degree of TCSC. Q sin θ δ δ sinθ (7)
Z Z
Putting, δ δ δ into equation 4 and solving it for V V wV C rand P S C rand G S (12)
S S V (13)
.
V V ZQ
V (8) Where,
To obtain real values of Vr in terms of Q, the equation must Sik = present position of agent i at kth iteration
have real roots. Thus the following conditions, which can be Sik+1= present position of agent i at (k+1)th iteration
used as a stability criterion, need to be satisfied: Pbesti = Particle best of agent i
V sin θ δ 4xQ 0 (10) Gbest = Global best of the set
Or
B. PSO-TVAC
Q
L 1.00 (11) PSO-TVAC is particle swarm optimization with time
V
varying acceleration co-efficient algorithm is the algorithm
worked on the principle same as that of PSO with the varying
L is termed the stability index of the line. The stability acceleration co-efficient in each of iteration in contrast with the
criterion is used to find the stability index for each line classical PSO where acceleration coefficients are kept constant.
connected between two bus bars in an interconnected network. PSO-TVAC enhance global search of particles in early parts of
Based on the stability indices of lines, voltage collapse can be optimization and allow the particles to converge towards global
accurately predicted. As long as the stability index L remains optimum solution at the end of the algorithm. To achieve so c1
less than 1, the system is stable and when this index exceeds and c2 are varied with time such that cognitive component is
the value 1, the whole system loses its stability and voltage large and social component is small initially and they are
collapse occurs. Thus the proposed method can be used in decreased and increased respectively. This condition allows the
voltage collapse prediction. particles to move around the search space in spite of moving
IV. PSO AND PSO-TVAC towards the Pbest initially and converge towards global optimal
in later parts.
A. Simple PSO
PSO employs a population of individuals known as particles. Following are the formulae of acceleration co-efficient
change in each of iteration.
These particles fly across the problem work field with some
specified initial velocities. The velocities of the particles are C1= C1i + (C1f – C1i)×(iter/max_iter) (14)
stochastically attuned considering, the past best location of the C2= C2i + (C2f – C2i)×(iter/max_iter) (15)
particles and their vicinity best position, in each iteration.
Here, C1i, C1f, C2i and C2f are initial and final values of
These positions are determined according to some predefined
cognitive and social acceleration co-efficient respectively. Due
fitness function and then, the movement of each particle as you
to such variation of acceleration co-efficient, PSO-TVAC
would expect emerges to an optimal or near-optimal solution.
converges faster than classical PSO [8].
Each particle follows the path of its coordinates in the
problem work space which are related with the best solution V. SIMULATION AND DISCUSSION
(fitness), it has attained until now. The fitness value is also Effectiveness of TCSC in the system is analyzed on a 6 bus
recorded and denoted as Pbest. When a particle takes all the test system (modification in standard IEEE 6-bus system) using
population as its geological neighbors, of which the best value power world simulator. This test system is consist of 6 buses,
is a global best and denotes as Gbest. three generators, three loads and eleven transmission line
connected as shown in Fig.3. Data of this system is given in
The particle revises its velocity and positions, with following
APPENDIX.
equation, after getting the two best values.
2. Calculation of LSI of each line for base case and single
line contingency (12 cases).
3. Implementation of PSO algorithm to find optimal rating
of TCSC at different locations. The objective function of the
algorithm is to minimize (Σ LSI).
Table II shows different PSO parameters values used in
simulation.
TABLE II. PSO PARAMETERS
Parameters Value
Number of Particles 15
Inertia weight Linearly decreased
Individual acceleration co- 2.5
efficient(C1)
Fig. 3. 6-Bus Test System
Social acceleration co-efficient(C2) 4-C1
From the given data, the model of the test system is Number of Iterations(max ) 20
prepared using power world simulator and the power flow is rand1 0.3
checked. With normal operating condition, LSI values for all
rand2 0.2
eleven lines are computed. This normal condition is termed as
Velocity bounds [-2,2]
the base case here. Then after some line outage contingencies
are considered and the power flow is analyzed. Again the LSI Position bounds(Search Space) [-0.7,0.2] (range of
r)
values of each line for these cases are computed. The line
C1 initial and C1 final(for PSO- 0.2 and 2.5
outage providing highest value of LSI is considered as the TVAC) respectively
most severe line outage from the viewpoint of voltage C2 initial and C2 final(for PSO- 2.5 and 0.2
stability. TVAC) respectively
From the Table I, it can be seen that LSI of line no. 8 is
A. Classical PSO applied to 6-Bus Test System
maximum in most of the cases. So compensation should
By applying PSO for different cases, to improve voltage
provide at this line. The TCSC having capacitive reactance
stability, location and optimal rating of TCSC can be found.
0.0091p.u.is placed at line 8. The rating of TCSC is found
Results for all cases are shown in Table III.
using trial and error method such that voltages improve at
each bus. Table shows that after placing TCSC voltage For Base Case, Line 8 is the critical line, the line having
stability is improved i.e. maximum LSI is reduced. maximum LSI value. Similarly for different line outage
contingencies, critical lines are different.
TABLE I. LSI VALUES WITH AND WITHOUT TCSC From these results it is observed that, location of TCSC is
not same as the weakest line for each case. But for all the cases
No Case Critical Max LSI Max LSI
Line with TCSC after placement of TCSC summation of LSI is reduced after
1 Base Case 8 0.3274 0.32061 placement of TCSC. So voltage stability is improved of power
2 Outage of Line 1 8 0.3602 0.3534 system.
3 Outage of Line 2 1 0.3895 0.3892
4 Outage of Line 3 8 0.404 0.3940 B. PSO-TVAC applied to 6-Bus Test System
5 Outage of Line 4 8 0.3471 0.34 With the same parameters as of classical PSO and only keep
6 Outage of Line 5 2 0.4602 0.4591 updating C1 and C2 in each iteration, as the principle of PSO-
7 Outage of Line 6 8 0.3824 0.374
TVAC, same analysis is carried out to find optimal location
From this it is clear that LSI can be used to measure the and rating of TCSC. Application of PSO-TVAC gives the same
voltage stability of the system. results as of PSO (Table III.), both TCSC location and rating.
Analysis is carried out using MATLAB using the following Thus PSO and PSO-TVAC can be used as the substitutes of
steps: each other.
APPENDIX
6 BUS TEST SYSTEM DATA
(BASE MVA=100)
(pu) (pu)
1 Slack 0 0 230 1.05 1.05
2 PV 0 0 230 1.05 1.05
3 PV 0 0 230 1.07 1.07
4 PQ 0.85 0.8 230 1.05 0.95
5 PQ 0.85 0.8 230 1.05 0.95
6 PQ 0.85 0.8 230 1.05 0.95