# A Fuzzy Based Evolutionary Algorithm for Solving Multiobjective Optimal Power Flow with FACTS Devices

Authors: R. Vanitha; J. Baskaran
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

A new Fuzzy Differential Evolution (FDE) algorithm is proposed for solving multiobjective optimal power flow with FACTS devices. This new optimization technique combines the advantages of Weighted Additive Fuzzy Goal Programming (WAFGP) and Differential Evolution (DE) in enhancing the capacity, stability, and security of the power system. As the weights used in WAFGP would have a significant impact on the operational and economical enhancements achieved in the optimization, they are optimized using evolutionary DE algorithm. This provides a way for exploring a balanced solution for a multiobjective problem without sacrificing any individual objective’s uniqueness and priority. The multiple objectives considered are maximizing the loadability condition of the power system with minimum system real power loss and minimum installation cost of the FACTS devices. Indian utility Neyveli Thermal Power Station (NTPS) 23 bus system is used to test the proposed algorithm using multiple FACTS devices. The results compared with that of DE based fuzzy goal programming (FGP) demonstrates that DE based WAFGP algorithm not only provides a balanced optimal solution for all objectives but also provides the best economical solution.


---

## Body

## 1. Introduction

The ever-increasing demand of electric power and deregulation norms followed in certain countries create a heavy traffic, and thereby straining an existing power network. Maintaining the stability and security of power network becomes the foremost challenge for the power suppliers.FACTS devices are power electronic devices that can be inserted in existing power network at desired location without altering the existing network. They provide required transmission line impedance, real power, and reactive power, maintain the voltage levels at buses which in turn increases the load carrying capacity of transmission network, reduces real and reactive power losses, and maintains voltage stability in the overall power network [1].The optimal location and sizing of FACTS controllers in a power network are done for better utilization of these devices in an economical way. Several conventional and intelligent natural inspired algorithms are available in obtaining a global optimal solution for a single objective optimization [2–6]. When more than one objective is considered for optimization, the problem turns out into a multiobjective optimization yielding Pareto optimal solutions. The fuzzy concept is one of the most popular methodology that can be combined with any conventional/intelligent techniques to bring out the best global optimal solution from pareto optimal solutions.Zadeh was the first to introduce the concept of fuzzy set theory in the year 1965 [7]. Narasimhan has illustrated the application of fuzzy subsets concepts to fuzzy goal programming in 1980 [8]. Tiwari et al. have formulated an additive model to solve fuzzy goal programming in 1987 [9]. This model is mostly suited for multiobjective problems where objectives have different levels of importance. Many researchers have used fuzzy technique with other optimization methods to solve multiobjective optimal power flow problems in an effective manner.Rosehart et al. have combined interior point methods with goal programming in optimizing active and reactive power dispatch while maximizing voltage security in power systems in 2003. In 2011, Hazra and Sinha have presented a multiobjective particle swarm optimization method for minimizing generation cost and environmental pollution simultaneously. Liang et al. (2011) have presented a fuzzy based hybrid particle swarm optimization (PSO) approach in minimizing fuel cost, total emission, and total real power loss [10–12].In 2011, Sivasubramani and Swarup have proposed multiobjective harmony search (MOHS) algorithm and He et al., in 2013, have used Artificial Bee Colony (ABC) algorithm to solve multiobjective OPF. Azizipanah-Abarghooee et al. have presented modified Shuffled frog-leaping algorithm (SFLA) in solving multiobjectives like reducing generation cost, decreasing transmission loss, improving voltage stability index, and power system security in 2014 [13–15]. In all the above papers, fuzzy membership is used to choose a comprise solution from the set of Pareto optimal solutions. In this paper, maximization and minimization goals for each objective in multiobjective OPF are determined using DE algorithm and are converted to fuzzy goals using fuzzy membership values. Each fuzzy goal is assigned with a weight and all fuzzy goals are combined together to frame the Weighted Additive Fuzzy Goal model. Taking the importance of each individual objective into consideration, the weights for each objective in multiobjective OPF is optimized using DE algorithm.This paper is organized into seven sections. Section1 provides an introduction. Section 2 presents the implementation of mathematical modelling of FACTS devices in power system network. Section 3 formulates this study’s problem. Section 4 describes the proposed FDE approach in detail for optimal location of FACTS in a multiobjective optimization. Section 5 presents and discusses the results, and Section 6 concludes the benefits of this study.


## 2. Implementation of Facts Devices in Power System Network

The generalized real and reactive power flow equations for a bus “i” without any FACTS devices connected in a power system can be written as:(1)
P
i
=
∑
j
=
1
nbr
V
i
V
j
Y
i
j
cos
⁡
δ
i
-
δ
j
-
θ
i
j
,
Q
i
=
∑
j
=
1
nbr
V
i
V
j
Y
i
j
sin
⁡
δ
i
-
δ
j
-
θ
i
j
,where nbr is the total number of transmission line branches in a power system; θ and Y are the angle and magnitude of bus admittance elements; |
V
|  and  δ are magnitude and phase angle of bus voltages; P and Q are real and reactive powers of the buses.The introduction of FACTS devices in the power systems makes a significant change in the power flow equations that depends on the mathematical model of FACTS devices.


### 2.1. Thyristor Controlled Series Capacitor

TCSC is a series compensated thyristor controlled device that is comprised of a capacitor in parallel with a TCR. Figure1 depicts the equivalent circuit of a TCSC. In this model, the desired transmission line impedance can be achieved for smoother reactive power control by controlling the TCSC both in inductive and capacitive modes.Figure 1
Equivalent circuit of TCSC.The transmission line reactance to which the TCSC is connected is only altered as shown in (2) that in turn makes the significant change in the power flows [16]:(2)
X
line
=
X
line
+
X
c
,where (3)
r
c
min
⁡
∗
X
line
≤
X
c
≤
r
c
max
⁡
∗
X
line
;
r
c is the amount of compensation that can be provided by the TCSC in inductive and capacitive mode.


### 2.2. STATCOM

STATCOM is a shunt-connected device that is comprised of a voltage source converter fed from a capacitor or a constant DC source. It can independently control real and reactive power with the power system. Reactive power can be supplied or absorbed by varying the magnitude of voltage and real power can be controlled by varying the phase angle of the voltage with respect to the bus voltage to which it is connected. The power injection model of STATCOM is shown in Figure2 [17]. In this model, the STATCOM is considered as a voltage source in series with an impedance Z
s.Figure 2
Power injection model of STATCOM.The real power and reactive power at busp to which the STATCOM is connected get altered with the addition of injected real and reactive powers of STATCOM as(4)
P
p
=
P
s
+
∑
q
=
1
nb
V
p
V
q
Y
p
q
cos
⁡
δ
p
-
δ
q
-
θ
p
q
,
Q
p
=
Q
s
+
∑
q
=
1
nb
V
p
V
q
Y
p
q
sin
⁡
δ
p
-
δ
q
-
θ
p
q
,
P
s
=
G
s
V
p
2
-
V
p
E
s
Y
s
cos
⁡
δ
p
-
δ
s
-
θ
s
,
Q
s
=
-
B
s
V
p
2
-
V
p
E
s
Y
s
sin
⁡
δ
p
-
δ
s
-
θ
s
,where |
E
s
| and δ
s are the magnitude and phase angle of STATCOM voltage source, respectively, and  |
Y
s
| and θ
s are the magnitude and phase angle of admittance in series with the STATCOM voltage source.In this model the Jacobian matrix in Newton Raphson load flow method is extended by adding a new row and column to determine the two new variables|
E
s
| and δ
s as in [17, 18]. The initial value of STATCOM voltage is considered as 1.0 p.u and the control parameter of STATCOM in this model is X
s and in turn varies with the real and reactive power injection/absorption of STATCOM with respect to the bus p connected.


### 2.3. Unified Power Flow Controller

UPFC is a combined FACTS controller which has two Voltage Source Converters (VSC) connected back to back through a DC link. It is the most sophisticated FACTS device that works like a series and/or shunt compensator, voltage regulator, and phase shifter. In this paper, the power injection model of a UPFC device proposed by Noroozian in the year 1997 is used.In this model, the series source converter is modeled as a voltage source in series with impedanceX
s as shown in Figure 3:(5)
V
p
′
¯
=
V
s
¯
+
V
p
¯
V
s
¯
=
r
V
p
¯
e
j
γ
,where r is the magnitude and γ is the phase angle of voltage injected by series converter of UPFC, and their value ranges are 0
<
r
<
r
max
⁡ and 0
<
γ
<
2
π. The power injection model of UPFC including the shunt converter is depicted in Figure 4.Figure 3
Equivalent circuit of series connected VSC.Figure 4
Power injection model of UPFC.The injected real and reactive powers of UPFC atpth and qth bus as in [19] are as follows:(6)
P
s
p
=
r
b
s
V
p
V
q
sin
⁡
θ
p
q
+
γ
,
Q
s
p
=
r
b
s
V
p
2
cos
⁡
γ
,
P
s
q
=
-
r
b
s
V
p
V
q
sin
⁡
θ
q
p
+
γ
,
Q
s
q
=
-
r
b
s
V
p
V
q
cos
⁡
θ
p
q
+
γ
,where b
s
=
1
/
X
s.This model can be easily implemented in Newton Raphson load flow method. If a UPFC is located between busp and bus q in a power system, the bus admittance matrix of the system is modified by adding a reactance equivalent to X
s between bus p and bus q. The Jacobian matrix is modified by addition of appropriate powers as in [19]. The control parameters in this model are r and γ.


## 2.1. Thyristor Controlled Series Capacitor

TCSC is a series compensated thyristor controlled device that is comprised of a capacitor in parallel with a TCR. Figure1 depicts the equivalent circuit of a TCSC. In this model, the desired transmission line impedance can be achieved for smoother reactive power control by controlling the TCSC both in inductive and capacitive modes.Figure 1
Equivalent circuit of TCSC.The transmission line reactance to which the TCSC is connected is only altered as shown in (2) that in turn makes the significant change in the power flows [16]:(2)
X
line
=
X
line
+
X
c
,where (3)
r
c
min
⁡
∗
X
line
≤
X
c
≤
r
c
max
⁡
∗
X
line
;
r
c is the amount of compensation that can be provided by the TCSC in inductive and capacitive mode.


## 2.2. STATCOM

STATCOM is a shunt-connected device that is comprised of a voltage source converter fed from a capacitor or a constant DC source. It can independently control real and reactive power with the power system. Reactive power can be supplied or absorbed by varying the magnitude of voltage and real power can be controlled by varying the phase angle of the voltage with respect to the bus voltage to which it is connected. The power injection model of STATCOM is shown in Figure2 [17]. In this model, the STATCOM is considered as a voltage source in series with an impedance Z
s.Figure 2
Power injection model of STATCOM.The real power and reactive power at busp to which the STATCOM is connected get altered with the addition of injected real and reactive powers of STATCOM as(4)
P
p
=
P
s
+
∑
q
=
1
nb
V
p
V
q
Y
p
q
cos
⁡
δ
p
-
δ
q
-
θ
p
q
,
Q
p
=
Q
s
+
∑
q
=
1
nb
V
p
V
q
Y
p
q
sin
⁡
δ
p
-
δ
q
-
θ
p
q
,
P
s
=
G
s
V
p
2
-
V
p
E
s
Y
s
cos
⁡
δ
p
-
δ
s
-
θ
s
,
Q
s
=
-
B
s
V
p
2
-
V
p
E
s
Y
s
sin
⁡
δ
p
-
δ
s
-
θ
s
,where |
E
s
| and δ
s are the magnitude and phase angle of STATCOM voltage source, respectively, and  |
Y
s
| and θ
s are the magnitude and phase angle of admittance in series with the STATCOM voltage source.In this model the Jacobian matrix in Newton Raphson load flow method is extended by adding a new row and column to determine the two new variables|
E
s
| and δ
s as in [17, 18]. The initial value of STATCOM voltage is considered as 1.0 p.u and the control parameter of STATCOM in this model is X
s and in turn varies with the real and reactive power injection/absorption of STATCOM with respect to the bus p connected.


## 2.3. Unified Power Flow Controller

UPFC is a combined FACTS controller which has two Voltage Source Converters (VSC) connected back to back through a DC link. It is the most sophisticated FACTS device that works like a series and/or shunt compensator, voltage regulator, and phase shifter. In this paper, the power injection model of a UPFC device proposed by Noroozian in the year 1997 is used.In this model, the series source converter is modeled as a voltage source in series with impedanceX
s as shown in Figure 3:(5)
V
p
′
¯
=
V
s
¯
+
V
p
¯
V
s
¯
=
r
V
p
¯
e
j
γ
,where r is the magnitude and γ is the phase angle of voltage injected by series converter of UPFC, and their value ranges are 0
<
r
<
r
max
⁡ and 0
<
γ
<
2
π. The power injection model of UPFC including the shunt converter is depicted in Figure 4.Figure 3
Equivalent circuit of series connected VSC.Figure 4
Power injection model of UPFC.The injected real and reactive powers of UPFC atpth and qth bus as in [19] are as follows:(6)
P
s
p
=
r
b
s
V
p
V
q
sin
⁡
θ
p
q
+
γ
,
Q
s
p
=
r
b
s
V
p
2
cos
⁡
γ
,
P
s
q
=
-
r
b
s
V
p
V
q
sin
⁡
θ
q
p
+
γ
,
Q
s
q
=
-
r
b
s
V
p
V
q
cos
⁡
θ
p
q
+
γ
,where b
s
=
1
/
X
s.This model can be easily implemented in Newton Raphson load flow method. If a UPFC is located between busp and bus q in a power system, the bus admittance matrix of the system is modified by adding a reactance equivalent to X
s between bus p and bus q. The Jacobian matrix is modified by addition of appropriate powers as in [19]. The control parameters in this model are r and γ.


## 3. Problem Formulation

The scope of this paper is to determine the optimal location of FACTS devices to obtain a maximum loadability in a power system with minimum total system real power loss and minimum installation cost of FACTS devices.The objective function is(7)
F
=
max
⁡
f
1
+
min
⁡
f
2
+
min
⁡
f
3
,where(8)
f
1
:  Loadability of the system
λ
,
f
2
:
P
Loss
=
∑
k
=
1
nl
g
i
v
s
2
+
v
r
2
-
2
v
s
v
r
cos
⁡
δ
s
-
δ
r
,and nl is the number of transmission lines; g is the conductance of the transmission line branch:(9)
f
3
:
I
C
=
C
∗
Q
∗
1000
US$
,
C
=
C
TCSC
+
C
STATCOM
+
C
UPFC
.Based on Siemens AG Database, the cost functions for TCSC, STATCOM, and UPFC are(10)
C
UPFC
=
0.0003
Q
2
-
0.2691
Q
+
188.22
,
C
TCSC
=
0.0015
Q
2
-
0.7130
Q
+
153.75
,
Q
=
Q
2
-
Q
1
,
C
STATCOM
=
0.000375
Q
2
-
0.3041
Q
+
162.4
,where Q is the operating reactive range of FACTS controllers in MVAR; IC is the installation cost of FACTS devices in US$; Q
1 and Q
2 are the reactive power flows in a line branch before and after placing FACTS device in MVAR. The above mentioned objectives are subject to the following constraints.Equality Constraint: Power Flow Equation(11)
P
G
p
-
P
D
p
-
∑
r
=
1
nb
V
p
V
q
Y
p
q
cos
⁡
θ
p
q
-
δ
p
+
δ
q
=
0
,
Q
G
p
-
Q
D
p
-
∑
r
=
1
nb
V
p
V
q
Y
p
q
sin
⁡
θ
p
q
-
δ
p
+
δ
q
=
0
,where nb is the total number of buses; P
G and Q
G are real and reactive power generations; P
D and Q
D are real and reactive power demands.Inequality Constraints(12)
P
G
p
min
⁡
≤
P
G
p
≤
P
G
p
max
⁡
,
Q
G
p
min
⁡
≤
Q
G
p
≤
Q
G
p
max
⁡
,
0.95
≤
V
p
≤
1.05
.FACTS Devices Constraints(13)
-
0.7
∗
X
line
≤
X
c
≤
0.2
∗
X
line
,
0.5
≤
X
S
≤
3.0
,
0
≤
r
≤
0.4
,
0
∘
≤
γ
≤
36
0
∘
.


## 4. Proposed Fuzzy DE Algorithm

### 4.1. Differential Evolution Algorithm

DE is a population-based stochastic optimization algorithm proposed by Storn and Price in the year 1996. It optimizes the real valued functions which are multidimensional, nonlinear, and discontinuous. The DE optimization is a multistage process with stages of initialization, mutation, recombination, and selection [20]. The implementations of various stages of DE in OPF are explained as follows.(1) Set NG as Number of Generations, NP as Number of Populations,n as Number of Vector Particles (control variables), minimum and maximum limits of each control variable, S as Scaling Factor [
0,1.2
], and CR as crossover constant [
0,1
]. The control parameters chosen in this work are generator’s real power, generator’s reactive power, generator’s voltages, location of FACTS devices in the system, FACTS devices control parameters (
r
c
,
X
s
,
r
,
γ
), and load factor (λ) of the system.(2) In this initialization step, a population of size (NP
×
n) is randomly generated comprised of NP population vectors. Each population vector (
X
j
;
j
=
1 : NP
) is comprised of i control variables (
Y
i
;
i
=
1 :
n
) and each vector particle (Control variable) is randomly generated within its maximum and minimum limits using (14) and the number of population vectors are maintained throughout all generations. X
j
G vector is otherwise called as parent vector or target vector:(14)
Y
i
G
=
Y
i
min
⁡
+
Y
i
max
⁡
-
Y
i
min
⁡
∗
rand
⁡
0,1
.(3) Determine the fitness function (objective function) for each population vector (X
j). In this paper the fitness functions are minimizing and maximizing the loadability, system real power loss, and installation cost of FACTS devices.(4) Randomly choose any three population vectors from the population pool (sayX
a, X
b, and X
c) and perform mutation operation using (15) to determine the mutated vector (X
M). This operation helps in introducing new vectors in existing population. S value is chosen in a way to have high convergence rate without trapping into local minima:(15)
X
M
G
=
X
a
G
+
S
X
b
G
~
X
c
G
,
(16)
X
a
≠
X
b
≠
X
c
.(5) A random number (η) is generated using binominal recombination and using this random number cross over operation is done to determine the trial vector (X
T) as in(17)
X
T
G
=
X
M
G
;
if
η
≤
CR,
X
j
G
;
if
η
>
CR.(6) New offsprings for next generation are produced by evaluating the fitness of trial vectors (X
T) and original population vectors or target vectors (X
j) using (18)
X
j
G
+
1
=
X
T
G
;
iff
X
T
G
≤
f
X
j
G
,
X
j
G
;
iff
X
j
G
≤
f
X
T
G
.(7) Repeat steps (4) to (6) until the stopping criteria are satisfied. In this paper Max Generation = NG is chosen as the stopping criteria.(8) Repeat steps (1) to (7) to determine the optimal solution (goal) for individual objectives. The goals for minimizing and maximizing individual objective are specified asZ
m
i and Z
n
i.


### 4.2. Weighted Additive Fuzzy Goal Programming

Fuzzy linear programming was used by Zimmermann in solving multiobjective optimization problems in the year 1978 [21]. Tiwari in the year 1987 modified this approach as Weighted Additive Fuzzy Goal Programming. In this method, each objective function (Z
i) of a multiobjective optimization problem is resolved to obtain its minimum (Z
m) and maximum (Z
n) goals individually satisfying all the constraints. All individual objectives are converted into fuzzy goals using their minimum and maximum goals as in (19) and (20), respectively [22]:(19)
μ
Z
m
y
=
1
;
Z
m
≤
Z
m
min
⁡
Z
m
max
⁡
-
Z
m
y
Z
m
max
⁡
-
Z
m
min
⁡
;
Z
m
min
⁡
≤
Z
m
y
≤
Z
m
max
⁡
0
;
Z
m
≥
Z
m
max
⁡
m
=
1,2
,
…
,
p
,
(20)
μ
Z
n
y
=
1
;
Z
n
≥
Z
n
max
⁡
Z
n
y
-
Z
n
min
⁡
Z
n
max
⁡
-
Z
n
min
⁡
;
Z
n
min
⁡
≤
Z
n
y
≤
Z
n
max
⁡
0
;
Z
n
≤
Z
n
min
⁡
n
=
p
+
1
,
p
+
2
,
…
,
q
.Fuzzy goals (FG) are combined together and WAFGP model is formulated as shown in (21) subjected to constraints as in (22) and solved to obtain optimal solutions. The optimal solutions depend on the weights (W) assigned to each fuzzy goal:(21)
Max
⁡
⁡
∑
i
=
1
q
W
i
F
G
isubject  to(22)
F
G
i
≤
μ
z
i
y
;
0
≤
μ
Z
i
i
≤
1
;
0
≤
W
i
≤
1
.


### 4.3. Implementation of FDE (WAFGP with DE) Algorithm in OPF

FDE method consists of following steps.(1) Set NG as Number of Generations, NP as Number of Populations,n as Number of Vector particles (control variables), minimum and maximum limits of each control variable, S as Scaling Factor [
0,1.2
], and CR as crossover constant [
0,1
].(2) Determine the minimum (Z
m) and maximium (Z
n) goal for individual objective using DE algorithm. In this work, minimum loadability (Z
m
1), maximum loadability (Z
n
1), minimum system real loss (Z
m
2), maximum system real loss (Z
n
2), minimum installation cost of FACTS devices (Z
m
3), and maximum installation cost of FACTS devices (Z
n
3) are obtained using DE algorithm.(3) Determine the fuzzy linear membership values for minimization goals (Z
m
i) and maximization goals (Z
n
i), where i is the number of objectives (
i
=
3
) using (19) and (20) for NP population vectors.(4) Formulate the WAFGP model as in (21):(23)
Max
⁡
W
1
F
G
1
+
W
2
F
G
2
+
W
3
F
G
3
.(5) Determine the optimal weights (W
1
,
W
2
,
W
3) using DE Algorithm for NP population vectors and find the optimal solution.(6) Repeat steps (3) to (5) for NG number of generations.(7) Obtain the best optimal maximum loadability of the system with minimum system real power loss and minimum installation cost of FACTS devices.


## 4.1. Differential Evolution Algorithm

DE is a population-based stochastic optimization algorithm proposed by Storn and Price in the year 1996. It optimizes the real valued functions which are multidimensional, nonlinear, and discontinuous. The DE optimization is a multistage process with stages of initialization, mutation, recombination, and selection [20]. The implementations of various stages of DE in OPF are explained as follows.(1) Set NG as Number of Generations, NP as Number of Populations,n as Number of Vector Particles (control variables), minimum and maximum limits of each control variable, S as Scaling Factor [
0,1.2
], and CR as crossover constant [
0,1
]. The control parameters chosen in this work are generator’s real power, generator’s reactive power, generator’s voltages, location of FACTS devices in the system, FACTS devices control parameters (
r
c
,
X
s
,
r
,
γ
), and load factor (λ) of the system.(2) In this initialization step, a population of size (NP
×
n) is randomly generated comprised of NP population vectors. Each population vector (
X
j
;
j
=
1 : NP
) is comprised of i control variables (
Y
i
;
i
=
1 :
n
) and each vector particle (Control variable) is randomly generated within its maximum and minimum limits using (14) and the number of population vectors are maintained throughout all generations. X
j
G vector is otherwise called as parent vector or target vector:(14)
Y
i
G
=
Y
i
min
⁡
+
Y
i
max
⁡
-
Y
i
min
⁡
∗
rand
⁡
0,1
.(3) Determine the fitness function (objective function) for each population vector (X
j). In this paper the fitness functions are minimizing and maximizing the loadability, system real power loss, and installation cost of FACTS devices.(4) Randomly choose any three population vectors from the population pool (sayX
a, X
b, and X
c) and perform mutation operation using (15) to determine the mutated vector (X
M). This operation helps in introducing new vectors in existing population. S value is chosen in a way to have high convergence rate without trapping into local minima:(15)
X
M
G
=
X
a
G
+
S
X
b
G
~
X
c
G
,
(16)
X
a
≠
X
b
≠
X
c
.(5) A random number (η) is generated using binominal recombination and using this random number cross over operation is done to determine the trial vector (X
T) as in(17)
X
T
G
=
X
M
G
;
if
η
≤
CR,
X
j
G
;
if
η
>
CR.(6) New offsprings for next generation are produced by evaluating the fitness of trial vectors (X
T) and original population vectors or target vectors (X
j) using (18)
X
j
G
+
1
=
X
T
G
;
iff
X
T
G
≤
f
X
j
G
,
X
j
G
;
iff
X
j
G
≤
f
X
T
G
.(7) Repeat steps (4) to (6) until the stopping criteria are satisfied. In this paper Max Generation = NG is chosen as the stopping criteria.(8) Repeat steps (1) to (7) to determine the optimal solution (goal) for individual objectives. The goals for minimizing and maximizing individual objective are specified asZ
m
i and Z
n
i.


## 4.2. Weighted Additive Fuzzy Goal Programming

Fuzzy linear programming was used by Zimmermann in solving multiobjective optimization problems in the year 1978 [21]. Tiwari in the year 1987 modified this approach as Weighted Additive Fuzzy Goal Programming. In this method, each objective function (Z
i) of a multiobjective optimization problem is resolved to obtain its minimum (Z
m) and maximum (Z
n) goals individually satisfying all the constraints. All individual objectives are converted into fuzzy goals using their minimum and maximum goals as in (19) and (20), respectively [22]:(19)
μ
Z
m
y
=
1
;
Z
m
≤
Z
m
min
⁡
Z
m
max
⁡
-
Z
m
y
Z
m
max
⁡
-
Z
m
min
⁡
;
Z
m
min
⁡
≤
Z
m
y
≤
Z
m
max
⁡
0
;
Z
m
≥
Z
m
max
⁡
m
=
1,2
,
…
,
p
,
(20)
μ
Z
n
y
=
1
;
Z
n
≥
Z
n
max
⁡
Z
n
y
-
Z
n
min
⁡
Z
n
max
⁡
-
Z
n
min
⁡
;
Z
n
min
⁡
≤
Z
n
y
≤
Z
n
max
⁡
0
;
Z
n
≤
Z
n
min
⁡
n
=
p
+
1
,
p
+
2
,
…
,
q
.Fuzzy goals (FG) are combined together and WAFGP model is formulated as shown in (21) subjected to constraints as in (22) and solved to obtain optimal solutions. The optimal solutions depend on the weights (W) assigned to each fuzzy goal:(21)
Max
⁡
⁡
∑
i
=
1
q
W
i
F
G
isubject  to(22)
F
G
i
≤
μ
z
i
y
;
0
≤
μ
Z
i
i
≤
1
;
0
≤
W
i
≤
1
.


## 4.3. Implementation of FDE (WAFGP with DE) Algorithm in OPF

FDE method consists of following steps.(1) Set NG as Number of Generations, NP as Number of Populations,n as Number of Vector particles (control variables), minimum and maximum limits of each control variable, S as Scaling Factor [
0,1.2
], and CR as crossover constant [
0,1
].(2) Determine the minimum (Z
m) and maximium (Z
n) goal for individual objective using DE algorithm. In this work, minimum loadability (Z
m
1), maximum loadability (Z
n
1), minimum system real loss (Z
m
2), maximum system real loss (Z
n
2), minimum installation cost of FACTS devices (Z
m
3), and maximum installation cost of FACTS devices (Z
n
3) are obtained using DE algorithm.(3) Determine the fuzzy linear membership values for minimization goals (Z
m
i) and maximization goals (Z
n
i), where i is the number of objectives (
i
=
3
) using (19) and (20) for NP population vectors.(4) Formulate the WAFGP model as in (21):(23)
Max
⁡
W
1
F
G
1
+
W
2
F
G
2
+
W
3
F
G
3
.(5) Determine the optimal weights (W
1
,
W
2
,
W
3) using DE Algorithm for NP population vectors and find the optimal solution.(6) Repeat steps (3) to (5) for NG number of generations.(7) Obtain the best optimal maximum loadability of the system with minimum system real power loss and minimum installation cost of FACTS devices.


## 5. Case Studies and Discussions

The proposed combined approach of WAFGP and DE algorithm for solving multiobjective problem is tested on Indian utility Neyveli Thermal Power Station (NTPS) 23 bus system [23] using Matlab 7.5. This practical thermal utility system has one slack bus, 3 generator buses, 19 load buses, and 22 transmission lines. The load factor is varied from 1.0 to 1.5 allowing a 50% increase in real power loading. The values chosen for NP and NG are 100 and 200, respectively, in DE and in RCGA. In DE algorithm the value chosen for SF is 0.6 and for CF is 0.7.The optimal locations of TCSC, STATCOM, and UPFC with their rated values are obtained to achieve the best optimal results for the objective functionsf
1, f
2, and f
3 individually using DE and RCGA. f
1 is the objective dealing with maximization of the load factor of the system, f
2 deals with the minimization of system real loss, and f
3 deals with the minimum installation cost of the FACTS devices. The optimal rated control values of FACTS devices along with their location using DE algorithm are depicted in Table 1 and the comparative results of best optimal solutions for various individual objectives with DE and RCGA are shown in Figures 5, 6, and 7.Table 1
Optimal location and rating of facts devices in NTPS-23 bus system.
FACTS devices
Best
f
1
Best
f
2
Best
f
3
TCSC
Location
# 1-2
# 1-5
# 1-2
r
c
−0.4828
−0.1910
−0.2451
STATCOM
Location
Bus 7
Bus 8
Bus 7
X
s
2.6782
2.8615
2.2981
UPFC
Location
# 3-13
# 3-13
# 3-13
r
0.2410
0.0238
0.1759
γ
297.1062
293.4171
333.3421
# denotes the transmission line branch of NTPS 23 bus system.Figure 5
Comparison results of maximum loadability for NTPS23 bus system.Figure 6
Comparison results of minimum real power loss for NTPS23 bus system.Figure 7
Comparison results of minimum installation cost of FACTS devices for NTPS23 bus system.The comparison results in Figures5, 6, and 7 demonstrate that DE algorithm performs better than GA and shows the potential of individual and combined FACTS devices in achieving the individual objectives. In the next level TCSC, STATCOM, and UPFC are used simultaneously to achieve the minimum and maximum optimal goals individually and the results are displayed in Table 2.Table 2
Best pareto optimal solutions for combined facts devices in NTPS-23 bus system using DE.
f
1 (%)
f
2 (×100 Mw)
f
3 (×106 US$)
Min (Z
m
1)
Max (Z
n
1)
Min (Z
m
2)
Max (Z
n
2)
Min (Z
m
3)
Max (Z
n
3)
Location of TCSC
# 3-1
# 1-2
# 3-1
# 1-2
# 1-2
# 1-5
Location of STATCOM
Bus 8
Bus 9
Bus 6
Bus 7
Bus 6
Bus 9
Location of UPFC
# 2-4
# 3-13
# 3-13
# 2-4
# 2-4
# 3-1
r
c
-
0.165
-
0.565
-
0.565
-
0.691
-
0.024
-
0.061
X
s
2.573
1.423
2.249
0.925
0.723
2.672
r
γ
0.287300.553
0.009345.110
0.299280.607
0.003342.548
0.039274.248
0.323291.103
f
1 (%)
100
141
100
139
100
140
f
2 (×100 Mw)
0.150
0.320
0.137
0.280
0.145
0.296
f
3 (×106 US$)
1.4123
22.172
1.2972
20.272
1.2859
19.017Results in Table2 show that when increasing the loadability of the system is considered as the only objective, the real power loss in the system is very high. When minimizing real power loss or installation cost of FACTS devices is considered as the only objective, the loadability of the system is low. This result highlights the need for balancing the multiple objectives rather than focusing on a single objective.To achieve this balance, DE based WAFGP multiobjective technique is used. This multiobjective technique aims to achieve the optimal results for maximum loadability, minimum real power loss, and minimum installation cost of FACTS devices by obtaining the optimal location and control parameter settings for TCSC, STATCOM, and UPFC combination.The minimum and maximum goals of individual objectives are converted into fuzzy numbers using fuzzy membership functions as in (19) and (20). WAFGP model is framed as in (21) and DE algorithm is used to optimize the weights. The final optimal results are depicted in Table 3.Table 3
Best optimal solutions using FDE for NTPS23 bus system.
Location of TCSC
# 3-13
r
c
-
0.5101
Location of STATCOM
Bus 7
X
s
2.8364
Location of UPFC
# 2-4
r
0.1880
γ
279.0781
WeightW
1
0.5011
f
1 (%)
139
WeightW
2
0.1465
f
2 (×100 Mw)
0.204
WeightW
3
0.3524
f
3 (×106 US$)
6.5779The results in Table3 show that different objectives are achieved in a balanced manner than achieving one objective at the cost of other objectives as in Table 2. To show the merits of optimal weights, the installation cost of combined FACTS devices computed using WAFGP is compared against FGP approach and the results are shown in Table 4. The results show that, for the same loadability condition, installation cost of FACTS devices obtained by WAFGP approach is significantly lower than the installation cost of FACTS devices obtained by FGP.Table 4
Comparison results for installation cost of facts devices for NTPS23 bus system.
Optimalλ  (%)
Optimal IC for FACTS devices (×106 US$)
With WAFGP optimal weights (FDE)
Without weights
139
6.5779
6.9982From the results of Tables3 and 4, it is concluded that DE based WAFGP approach provides an optimal solution for a multiobjective power flow problems and can be implemented for any higher order systems with any number of objectives.


## 6. Conclusion

In this paper, analysis performed on the Indian utility NTPS 23 bus system demonstrates that the best optimal solution for single objectives are obtained by optimal placing and configuration of TCSC, STATCOM, and UPFC devices using DE algorithm when compared to that of RCGA. The results also indicate that when multiple objectives have to be optimized, the focus on a single objective has severe negative impact on other objectives. To consider multiple objectives simultaneously, WAFGP technique is used. Since the weights in WAFGP have a significant importance in deciding the optimal solution, DE algorithm is used for optimizing the weights. The single optimal solution for multiobjective OPF by the DE based WAFGP algorithm clearly indicates that each objective obtains its optimal solution without compromising other objectives. Thus, this algorithm provides a balance in optimizing multiple objectives and thereby preserves the uniqueness of each individual objective. Further, comparison of results with that of DE based FGP shows that the use of WAFGP with optimized weights results in minimum installation cost of FACTS devices and thereby better economical solution.



---

*Source: 275129-2015-08-04.xml*