Next Article in Journal
Effect of Temperature on Space Charge Distribution of Oil–Paper Insulation under DC and Polarity Reversal Voltage
Next Article in Special Issue
Dynamic Optimization of Combined Cooling, Heating, and Power Systems with Energy Storage Units
Previous Article in Journal
Thermal Conductivity of Korean Compacted Bentonite Buffer Materials for a Nuclear Waste Repository
Previous Article in Special Issue
Adaptive Consensus Algorithm for Distributed Heat-Electricity Energy Management of an Islanded Microgrid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Hybrid DA-PSO Optimization Algorithm for Multiobjective Optimal Power Flow Problems

by
Sirote Khunkitti
1,
Apirat Siritaratiwat
1,
Suttichai Premrudeepreechacharn
2,
Rongrit Chatthaworn
1,* and
Neville R. Watson
3
1
Department of Electrical Engineering, Faculty of Engineering, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Electrical Engineering, Faculty of Engineering, Chiang Mai University, Chiang Mai 50200, Thailand
3
Department of Electrical and Computer Engineering, University of Canterbury, Christchurch 8140, New Zealand
*
Author to whom correspondence should be addressed.
Energies 2018, 11(9), 2270; https://doi.org/10.3390/en11092270
Submission received: 3 August 2018 / Revised: 20 August 2018 / Accepted: 27 August 2018 / Published: 29 August 2018
(This article belongs to the Special Issue Optimization Methods Applied to Power Systems)

Abstract

:
In this paper, a hybrid optimization algorithm is proposed to solve multiobjective optimal power flow problems (MO-OPF) in a power system. The hybrid algorithm, named DA-PSO, combines the frameworks of the dragonfly algorithm (DA) and particle swarm optimization (PSO) to find the optimized solutions for the power system. The hybrid algorithm adopts the exploration and exploitation phases of the DA and PSO algorithms, respectively, and was implemented to solve the MO-OPF problem. The objective functions of the OPF were minimization of fuel cost, emissions, and transmission losses. The standard IEEE 30-bus and 57-bus systems were employed to investigate the performance of the proposed algorithm. The simulation results were compared with those in the literature to show the superiority of the proposed algorithm over several other algorithms; however, the time computation of DA-PSO is slower than DA and PSO due to the sequential computation of DA and PSO.

1. Introduction

For the past few decades, the optimal power flow (OPF) problem has played an essential role in studying the economy terms of power systems [1,2]. The OPF problem is a nonlinear, nonconvex, large-scale, and static programming problem [3] that optimizes selected objective functions while satisfying a set of equality and inequality constraints. The power balance equations are the equality constraints, and the limits of state and control variables are the inequality constraints of the OPF problem. The state variables consist of slack bus active power generation, load bus voltages, reactive power generation, and apparent power flow. The control variables involve active power generation except at slack bus, generator bus voltages, tap ratios of transformers, and reactive powers of shunt compensation capacitors. In recent years, because of the rise in fuel cost, which increases generation cost, fuel cost has become the objective function to be optimized in the OPF problem. Moreover, due to the release of emissions from thermal power plants into the atmosphere, emissions are yet another concern for power system operation and planning [4]. At the same time, because the demand for electricity has outpaced the expansion of transmission capacity, the inadequate reactive power sources of power systems have increased losses in transmission lines. Thus, emissions and transmission losses must also be considered as part of the objective functions of the OPF problem.
To solve the OPF problem, several traditional optimization techniques, such as nonlinear programming [5], quadratic programming [6], and the interior point method [7], have been successfully applied. However, these algorithms’ nonlinear characteristics make them impractical to use in practical systems. The nonlinear characteristics may cause the obtained solutions to be trapped in local optima, and these algorithms require an enormous amount of computational effort and time. Therefore, many optimization methods need to be improved to overcome these shortcomings [8,9]. Recently, several population-based optimization algorithms, including the OPF problem, have been employed to solve a complex constrained optimization problem in the field of power systems. Some of the other proposed techniques include the genetic algorithm (GA) [10], tabu search (TS) [11], differential evolution (DE) [12], evolutionary programming (EP) [13], probabilistic optimal power-flow (P-OPF) [14], preventive security-constrained power flow optimization [15], ant colony optimization (ACO) [16], grey wolf optimizer (GWO) [17], artificial bee colony (ABC) [18], particle swarm optimization (PSO) [19], and the dragonfly algorithm (DA) [20]. Even with the successful optimization of single-objective population-based optimization techniques, minimizing only one objective function is not sufficient in the power system because there are many problems, such as fuel cost, emissions, and transmission losses, which also need to be minimized. Consequently, many objective functions should be considered because this is a multi objective optimization problem. Since there are three independent objective functions in this study (i.e., fuel cost, emissions, and transmission losses), the number of incompatible optimal solutions between the objective functions is infinite, and these optimal solutions are called Pareto optimal solutions [21].
Several optimization algorithms have been proposed and applied to solve the multiobjective OPF (MO-OPF) problem by many researchers. One of these methods was carried out by converting the multiobjective problem into a single-objective problem and then solving the problem by using a single-objective optimizer. However, this method has some drawbacks, such as the limitation of the available choices, the need for weights for each objective, and the requirement of multiple optimizer runs. To overcome these weaknesses, many researchers have proposed multiobjective evolutionary algorithms, such as the improved strength Pareto evolutionary algorithm (ISPEA2) [22], hybrid modified particle swarm optimization-shuffle frog leaping algorithms (HMPSO-SFLA) [23], modified teaching–learning-based optimization (MTLBO) [24], GWO [17], DE [17], multiobjective modified imperialist competitive algorithm (MOMICA) [25], differential search algorithm (DSA) [26], modified shuffle frog leaping algorithm (MSFLA) [27], modified Gaussian bare-bones multiobjective imperialist competitive algorithm (MGBICA) [28], multiobjective harmony search (MOHS) [29], adaptive real coded biogeography-based optimization (ARCBBO) [30], multiobjective differential evolution algorithm (MO-DEA) [31], hybrid modified imperialist competitive algorithm and teaching–learning algorithm (MICA-TLA) [32], etc., to successfully solve the OPF problem. In the past few decades, various well-proposed multiobjective evolutionary algorithms have been successfully applied and improved in many applications; however, most of them have not been extensively investigated in the OPF problem. Moreover, improving the search performance of the multiobjective evolutionary algorithm for solving the OPF problem is also important. In this paper, a hybrid DA-PSO algorithm is proposed to deal with the MO-OPF problem. The concept of the hybrid algorithm is the combination of the exploration and exploitation phases of the DA and PSO algorithms, respectively. The performance of the proposed algorithm was evaluated on the standard IEEE 30-bus and IEEE 57-bus power systems. Three different objective functions—fuel cost, emissions, and transmission losses—were individually and simultaneously considered as parts of the objective function in the OPF problem. The obtained results were compared with other evolutionary algorithms and the traditional DA and PSO.
The rest of the article is classified into five sections as follows. Section 2 introduces the formulation and constraints of the multiobjective optimization. In Section 3, the traditional DA and PSO are explained, and Section 4 depicts the concept of the proposed algorithm. Section 5 presents the optimization results and the comparisons between the solutions from the proposed algorithm and the solution from other algorithms based on IEEE 30-bus and IEEE 57-bus systems. Finally, in Section 6, the conclusions of the simulation results of the proposed algorithm are described.

2. Problem Formulation and Constraints for Multi Objective Optimization for OPF

Multi-objective optimization is a model that optimizes more than one objective function to find optimal control variables while simultaneously satisfying equality and inequality constraints. The compromised solutions, nondominated solutions, which have more than one optimal solution between each objective, are the optimal solutions referred to as the Pareto front. The multiobjective problem is mathematically formulated as follows:
min   f = { f 1 ( x , u ) , f 2 ( x , u ) , , f N o b j ( x , u ) }
subject to
g ( x , u ) = 0
h ( x , u ) 0
where f is a vector of objective functions to be optimized, Nobj is the number of objective functions, g(x,u) are the equality constraints, and h(x,u) are the inequality constraints.
x is a vector of state variables including slack bus active power, load bus voltages, generator reactive powers, and apparent power flows, expressed as follows:
x = [ P g s l a c k , V L 1 , , V L N L , Q g 1 Q g N g e n , S l 1 S l N l ]
where Pgslack is the active power generation at slack bus, VLi is the load voltage at bus i, NL is number of load buses, Qgi is the reactive power generation at bus i, Ngen is the number of total generators, Sli is the apparent power flow at branch i, and Nl is the number of transmission lines.
u is a vector of control variables consisting of active power generations except at slack bus, generator bus voltages, transformer tap ratios, and reactive powers of shunt compensation capacitors, expressed as:
u = [ P g i ; i P V b u s P g N g e n , V g 1 , , V g N g e n , T 1 T N t r a n , Q c 1 Q c N c a p ]
where Pgi is the active power generation at bus i, PVbus is the set of generator buses except at slack bus, Vgi is the generator bus voltage at bus i, Ti is the transformer tap ratio at bus i, Ntran is the number of transformer taps, Qci is the shunt compensation capacitor at bus i, and Ncap is the number of compensation capacitors.

2.1. Objective Functions

In this study, the objective functions of the OPF, consisting of fuel cost, emissions, and transmission line losses, are considered as shown below.

2.1.1. Fuel Cost

The total fuel cost of the generators is considered to be minimized and is given as follows:
f C ( x , u ) = i = 1 N g e n ( a i P g i 2 + b i P g i + c i )
where fC is the total fuel cost of generators function ($/h), and ai, bi and ci are the fuel cost coefficients of the ith generator units.

2.1.2. Emissions

The emissions function can be represented as the sum of all considered emission types, such as sulphur oxides (SOx), nitrogen oxides (NOx), thermal emission, etc. However, in the present study, two important emission types, NOx and SOx, are taken into account, as expressed below:
f E ( x , u ) = i = 1 N g e n ( γ i P g i 2 + β i P g i + α i + ξ i exp ( λ i P g i ) )
where fE is the total emission generations function (ton/h), and γi, βi, αi, ζi and λi are emission coefficients of the ith generator units.

2.1.3. Transmission Line Losses

The system active power loss in the transmission line is formulated as follows:
f L ( x , u ) = k = 1 N l g k ( V i 2 + V j 2 2 V i V j cos ( θ i j ) )
where fL is the total transmission loss function (MW), gk is the conductance of the kth line, Vi is the voltage at bus i, Vj is the voltage at bus j, and θij is the voltage phase angle difference between buses i and j.

2.2. Constraints

2.2.1. Equality Constraints

The OPF equality constraints are the active and reactive power balance constraints, as follows:
P g i P d i = j = 1 N b u s V i V j ( G i j cos ( θ i j ) + B i j sin ( θ i j ) )
Q g i Q d i = j = 1 N b u s V i V j ( G i j sin ( θ i j ) B i j cos ( θ i j ) )
where Pdi is the active power demand at bus i, Nbus is the number of buses, Gij is the transfer conductance between buses i and j, Bij is the transfer susceptance between buses i and j, and Qdi is the reactive power demand at bus i.

2.2.2. Inequality Constraints

P g i min P g i P g i max i = 1 , 2 , , N g e n
Q g i min Q g i Q g i max i = 1 , 2 , , N g e n
V g i min V g i V g i max i = 1 , 2 , , N g e n
| S l i | S l i max
V L i min V L i V L i max i = 1 , 2 , , N L
Q c i min Q c i Q c i max i = 1 , 2 , , N c a p
T i min T i T i max i = 1 , 2 , , N t r a n
where Pgimin and Pgimax are the minimum and maximum active power generations at bus i, respectively, Qgimin and Qgimax are the minimum and maximum reactive power generations at bus i, respectively, Vgimin, Vgimax are the minimum and maximum generator voltage at bus i, respectively, Slimax is the maximum apparent power flow at branch i, VLimin, VLimax are the minimum and maximum load voltage at bus i, respectively, Qcimin and Qcimax are the minimum and maximum shunt compensation capacitor at bus i, respectively, Timin, Timax are the minimum and maximum transformer tap-ratio at bus i, respectively.

2.2.3. Constraints Handling

The inequality of dependent variables, including slack bus active power generation, load bus voltage magnitudes, reactive power generations, and apparent power flows, are integrated into the penalized objective function to maintain these variables within their limits and to refuse infeasible solutions. The penalty function can be expressed as follows [27]:
J ( x , u ) = f ( x , u ) + K P ( P g s l a c k P g s l a c k lim ) 2 + K V i = 1 N l o a d ( V L i V L i lim ) 2 + K Q i = 1 N l i n e ( Q g i Q g i lim ) 2 + K S i = 1 N l i n e ( S L i S L i max ) 2
where J(x,u) is the penalized objective function, Kp, KQ, KV and Ks are the penalty factors, and xlim is the limit value of the dependent variables, determined as follows:
x lim = { x max x x min   i f i f i f   x > x max x min < x < x max x < x min

3. Related Optimization Techniques

3.1. DA

DA is a metaheuristic algorithm which was inspired by the static and dynamic swarming behaviors of dragonflies in nature [33]. Dragonflies swarm for two goals: Hunting (static swarm) and migration (dynamic swarm). In the dynamic swarm, many dragonflies swarm when roaming over long distances and different areas, which is the purpose of the exploration phase. In the static swarm, dragonflies move in larger swarms and along one direction with local movements and sudden changes in the flying path, which is suitable in the exploitation phase.
The behavior of dragonflies can be represented through five principles, which are separation, alignment, cohesion, attraction to a food source, and distraction of an enemy. These five behaviors are described and calculated as follows:
Separation, which is the avoidance of the static crashing of individuals into other individuals in the neighborhood, is calculated by Equation (20).
S i = j = 1 N X X j
where Si is the separation of the ith individual, N is the number of neighboring individuals, X is the position of the current individual, and Xj is the position of jth neighboring individual.
Alignment, which refers to the velocity matching of individuals to the velocity of others in the neighborhood, is computed by Equation (21).
A i = j = 1 N V j N
where Ai is the alignment of the ith individual, and Vj is the velocity of jth neighboring individual.
Cohesion, which is the propensity of individuals towards the center of mass of the neighborhood, is formulated by Equation (22).
C i = j = 1 N X j N X
where Ci is the cohesion of the ith individual
Attraction towards a food source computed by Equation (23), should be the main objective of any swarm to survive.
F i = X + X
where Fi is the food source of the ith individual, and X+ is the position of the food source.
Distraction of an enemy, which is computed by Equation (24), is another survival objective of the swarm.
E i = X + X
where Ei is the position of enemy of the ith individual, and X is the position of the enemy source.
To simulate the movement of artificial dragonflies and update their positions, step vector (ΔX) and position vector (X) are considered. The step vector represents the direction of the movement of the artificial dragonflies and is formulated as follows:
Δ X t + 1 = ( s S i + a A i + c C i + f F i + e E i ) + ω t Δ X t
where ΔXt + 1 is the step vector at iteration t + 1, ΔXt is the step vector at iteration t, s, a, c, f and e are the separation weight, alignment weight, cohesion weight, food factor and enemy factor, respectively, and ωt is the inertia weight factor at iteration t and is calculated by Equation (26).
ω t = ω max ω max ω min I t e r max × I t e r
where ωmax and ωmin are set to 0.9 and 0.4, respectively, Iter is the iteration, and Itermax is the maximum iteration.
The position of the artificial dragonflies can be updated by the following equation:
X t + 1 = X t + Δ X t + 1
where Xt+ 1 is the position at iteration t + 1, and Xt is the position at iteration t.
When the search space does not have a neighboring solution, the artificial dragonflies need to move around the search space by applying random walk (Levy flight) to improve their stochastic behavior. So, in this case, the position of the dragonflies can be calculated by Equation (28).
X t + 1 = X t + L e v y ( d ) × X t
where d is the dimension of the position vectors, and the Levy is the Levy flight which is computed by Equation (29).
L e v y ( d ) = 0.01 × r 1 × σ | r 2 | 1 β
where r1 and r2 are two uniform random values in a range of [0, 1], and σ is calculated by Equation (30).
σ = ( Γ ( 1 + β ) × sin ( π β   2 ) Γ ( 1 + β 2 ) × β × 2 ( β 1 2 ) ) 1 / β
where β is the constant (which is equal to 1.5 in this work), and Γ ( x ) = ( x 1 ) ! .

3.2. PSO

PSO is a population-based stochastic global optimization technique which was first introduced by Eberhart and Kennedy [34]. The idea of PSO came from the flocking behavior of birds or the schooling of fishes in their food hunting. In the PSO system, the population moves around a multidimensional search space where each particle represents a possible solution. Each particle contains the information of control variables and is associated with a fitness value that indicates its performance in the fitness space. Each particle i consists of its position Xi = (xi,1, xi,2, …, xi,Nvar), where Nvar represents the number of control variables, velocity Vi = (vi,1, vi,2, …, vi,Nvar) and personal best experience Xpbesti = (xpbesti,1, xpbesti,2, …, xpbesti,Nvar), and a swarm has a global best experience Xgbest = (xgbest1, xgbest2, …, xgbestNvar). During each iteration, each particle moves in the direction of its own personal best position provided so far as well as in the direction of the global best position obtained so far by particles in the swarm. The particles are operated according to the equations expressed as follows:
V i t + 1 = ω t × V i t + C 1 × r a n d 1 × ( X p b e s t i t X i t ) + C 2 × r a n d 2 × ( X g b e s t t X i t )
X i t + 1 = X i t + V i t + 1
where Vit+ 1 is the velocity of particle i at iteration t + 1, Vit is the velocity of particle i at iteration t, C1 and C2 are two positive acceleration constants, rand1 and rand2 are two uniform random values in a range of [0, 1], Xtpbesti is the personal best position of particle i at iteration t, Xit is the position of particle i at iteration t, Xtgbest is the global best position among all particles at iteration t, and Xit+ 1 is the position of particle i at iteration t + 1.

4. Proposed Hybrid DA-PSO Optimization Algorithm for MO-OPF Problem

Many optimization algorithms have been proposed to overcome the optimization problem of being trapped in the local optima while the algorithms try to find the best solution. PSO has been proven in several works from the literature to find the optimal solution in various problems [35,36,37,38]. Because of its equations in finding the optimal solution by using the best experience of the particles, PSO could quickly converge on the optimal solution, i.e., it is good at exploitation. However, PSO is sometimes still trapped in the local optima because it converges on the optimal solution too quickly. In other words, PSO is poor at exploration, which is an important task of the optimization process. In DA, it applies the Levy flight to improve the randomness and stochastic behavior when there is no neighboring dragonfly. This could significantly improve the exploration process of the algorithm. However, the best experience, which is the personal best, of dragonflies is not applied during the operation. This causes the DA to converge on the optimal solution very slowly and can sometimes cause it to be trapped in the local optima. To overcome these problems, a new algorithm is proposed which combines the prominent points of the DA and PSO algorithms, which are the exploration of DA and the exploitation of PSO. At first, the dragonflies in DA are initialized to explore the search space to find the area of the global solution. Then, the best position of DA is obtained. The obtained best position from DA is then substituted as the global best position in the PSO equation (Equation (31)). After that, the PSO algorithm, which is the exploitation phase, operates by using the global best position from DA, allowing it to provide the expected optimal solution. The velocity and position equations of PSO can be modified as follows:
V i t + 1 = ω t × V i t + C 1 × r a n d 1 × ( X p b e s t i t X i t ) + C 2 × r a n d 2 × ( X D A t + 1 X i t )
X i t + 1 = X i t + V i t + 1
where X D A t + 1 is the best position obtained from DA at iteration t + 1.
The application of the proposed DA-PSO algorithm for solving the MO-OPF problem can be described as follows:
Step 1.
Clarify the system data comprising the fuel cost coefficients of the generators, emission coefficients of the generators, initial values of generator active powers, initial values of generator bus voltages, initial values of transformer tap ratios, initial values of shunt compensation capacitors, upper limit of Sli, lower and upper limits of Pgi, Qgi, Vgi, VLi, Qci, and Ti, the parameters of DA and PSO, the number of dragonflies and particles, the number of iterations, and the archive size.
Step 2.
Generate the initial population of dragonflies and particles.
Step 3.
Convert the constrained multi objective problem to an unconstrained one by using Equation (18).
Step 4.
Perform the power flow and calculate the objective functions for the initial population of dragonflies.
Step 5.
Find the nondominated solutions and save them to the initial archive.
Step 6.
Set the fitness value of the initial population as the food source.
Step 7.
Calculate the parameters of DA (s, a, c, f, and e).
Step 8.
Update the food source and enemy of DA.
Step 9.
Calculate the S, A, C, F, and E by Equations (20)–(24).
Step 10.
Check if a dragonfly has at least one neighboring dragonfly, then update step vector (ΔX) and the position of dragonfly (XDA) by Equations (25) and (27), respectively, and if each dragonfly has no neighboring dragonfly, then update XDA by Equation (28) and set ΔX to be zero.
Step 11.
If any component of each population breaks its limit, then ΔX or XDA of that population is moved into its minimum/maximum limit.
Step 12.
Set the best position obtained from DA as the global best of PSO (Xgbest).
Step 13.
Update the velocity of the particle (V) and the position of the particle (XPSO) by Equations (33) and (34), respectively.
Step 14.
If any component of each population breaks its limit, then V or XPSO of that population is moved into its minimum/maximum limit.
Step 15.
Calculate the objective functions of the new produced population.
Step 16.
Employ the Pareto front method to save the nondominated solutions to the archive and update the archive.
Step 17.
If the maximum number of iterations is reached, the algorithm is stopped; otherwise, go to step 7.
The flowchart of the DA-PSO algorithm for the MO-OPF problem is shown in Figure 1.

5. Simulation Results

To investigate the performance of the proposed algorithm, the IEEE 30-bus and IEEE 57-bus test systems were employed. The proposed algorithm operated for 30 independent runs for each test system. To validate the superiority of the proposed algorithm for solving the economic dispatch optimization problem, the results provided by the proposed algorithm were compared with those of other metaheuristic algorithms from the literature. In order to investigate both the single-objective optimization and multiobjective optimization, the simulation was divided into two cases. Single-objective optimization was evaluated in the first case. In the second case, multiobjective optimization to solve the MO-OPF by using the proposed DA-PSO algorithm was evaluated.

5.1. IEEE 30-Bus Test System

The proposed DA-PSO algorithm was applied to the IEEE 30-bus system to evaluate its performance. The IEEE 30-bus test system was composed of 6 generators at buses 1, 2, 5, 8, 11, and 13, 4 transformers between buses 6 and 9, buses 6 and 10, buses 4 and 12, and buses 27 and 28, and 41 transmission lines. The total system demand was 283.4 MW and 126.6 MVAR. The bus and branch data is given in [39]. The population number and the size of the Pareto archive were set to be 100 and 100, respectively.

5.1.1. Single-Objective OPF

To evaluate the performance of the proposed algorithm for solving the single-objective optimization, three different objective functions consisting of fuel cost, emissions, and transmission loss minimizations were individually considered as part of the objective function. The obtained results by the traditional DA, PSO, and the proposed DA-PSO algorithm for three individual objective functions are shown in Table 1. In Table 2, the best fuel cost of generators provided by the DA-PSO algorithm are compared with other algorithms in the literature, including TS [11], EP [13], ACO [16], SFLA [27], MSFLA [27], improved evolutionary programming (IEP) [40], modified differential evolution optimal power flow (MDE-OPF) [41], stochastic genetic algorithm (SGA) [42], evolutionary-programming-based optimal power flow (EP-OPF) [43], honey bee mating optimization (HBMO) [44], PSO, and DA. The comparison of the best emission values of the DA-PSO algorithm with various algorithms in the literature, including ACO [16], HMPSO-SFLA [23], TLBO [24], MTLBO [24], DSA [26], MSFLA [27], SFLA [27], GA [27], GBICA [28], improved particle swarm optimization (IPSO) [45], PSO, and DA, is shown in Table 3. In Table 4, the best transmission losses provided by DA-PSO algorithm are compared with other algorithms in the literature, including GWO [17], DE [17], MOHS [29], enhanced genetic algorithm with decoupled quadratic load flow (EGA-DQLF) [46], efficient evolutionary algorithm (EEA) [47], enhanced genetic algorithm (EGA) [47], PSO, and DA. It can be seen that the proposed DA-PSO provided better results compared with those of other algorithms for all three objective functions, which can be confirmed by the results in Table 1, Table 2, Table 3 and Table 4. However, the computation time of the proposed DA-PSO is much slower than other algorithms in the literature because the proposed algorithm consumed the sequential computation time of DA and PSO as presented in Table 2, Table 3 and Table 4.

5.1.2. MO-OPF

In this subsection, the proposed algorithm is investigated as a multiobjective optimization problem, while every two and three objective functions are optimized simultaneously. The best two-dimensional Pareto fronts obtained from the DA, PSO, and DA-PSO algorithms for the IEEE 30-bus system are shown in Figure 2, Figure 3 and Figure 4. However, DA could not provide the convergent Pareto front when simultaneously considering the emissions and losses as parts of the objective function. This shows that DA is suitable for some objective functions, but that it is not suitable for every objective function for finding optimal solutions. In Figure 5, the Pareto front provided by the DA-PSO algorithm for the three-dimensional Pareto front is shown. For all figures in this system, most of the nondominated solutions obtained by the DA-PSO algorithm are better than those from the DA and PSO algorithms. For instance, at the same level of the fuel cost, the emissions provided by DA-PSO are less than those of DA and PSO. This shows that the new proposed hybrid DA-PSO algorithm, which adopts the exploration phase of the DA and the exploitation phase of the PSO, could improve the performance of the original DA and PSO algorithms.

5.2. IEEE 57-Bus Test System

The proposed hybrid DA-PSO was also tested on the IEEE 57-bus system to investigate its performance. The system active and reactive power demands were 1250.8 MW and 336.4 MVAR, respectively. It consisted of 7 generators located at buses 1, 2, 3, 6, 8, 9, and 12, 15 transformers, and 80 transmission lines. The detail data were taken from [48]. The population number was 100 and the size of the Pareto archive was 100.

5.2.1. Single-Objective OPF

To verify its performance for solving the single-objective OPF in a larger system, the proposed algorithm was also applied to the IEEE 57-bus test system. Three different objective functions, i.e., fuel cost, emissions, and transmission losses, were individually considered as part of the objective function. The results provided by DA, PSO, and the proposed DA-PSO algorithm for the three individual objectives are shown in Table 5. The best results from DA-PSO are compared with those of: MTLBO [23], DSA [25], GBICA [27], MGBICA [27], ARCBBO [29], MO-DEA [30], MICA-TLA [31], TLBO [48], Levy mutation teaching–learning-based optimization (LTLBO) [49], new particle swarm optimization (NPSO) [50], fuzzy genetic algorithm (Fuzzy-GA) [51], differential evolution pattern search (DE-PS) [52], ABC [53], particle swarm optimization algorithm with linearly decreasing inertia weight (LDI-PSO) [53], evolving ant direction differential evolution (EADDE) [54], gravitational search algorithm (GSA) [55], adaptive particle swarm optimization strategy (APSO) [56], PSO, and DA for the fuel cost objective function; GBICA [27], MGBICA [27], PSO, and DA for the emission objective function; and PSO and DA for the transmission loss objective function—all of which is summarized in Table 5, Table 6, Table 7 and Table 8. From these tables, it is obvious that the proposed algorithm could provide more optimized results than the compared algorithms for all three objective functions.

5.2.2. MO-OPF

This case proposes a multiobjective optimization problem by using the proposed DA-PSO algorithm to evaluate its performance for the IEEE 57-bus test system. The two-dimensional Pareto fronts provided by the PSO and DA-PSO algorithms for this system are shown in Figure 6, Figure 7 and Figure 8, while DA could not provide the convergent Pareto fronts for any multiobjective functions in this system. In Figure 9, the three-dimensional Pareto front obtained from DA-PSO is shown. From all figures for this system, the fronts obtained from DA-PSO algorithm are superior to those from PSO, while the fronts obtained from DA could not converge because the best experience of dragonflies in DA is not applied during the operation and the obtained solutions are trapped in the local optima. From the results, it can be seen that the proposed hybrid DA-PSO performs better than the original DA and PSO algorithms once again.

6. Conclusions

In this paper, a hybrid DA-PSO algorithm is proposed to solve the MO-OPF problem in a power system. As the DA is an algorithm that applies Levy flight to improve its randomness and stochastic behavior, this could significantly develop the exploration phase of the algorithm in an optimization. The PSO could quickly converge on the optimal solution because of its equations for finding optimal solutions by using the best experience of the particles. This makes PSO perform well at the exploitation phase in an optimization. The new hybrid DA-PSO algorithm combines the prominent points of these two algorithms, which are the exploration phase of DA and the exploitation phase of the PSO, to improve its performance for finding the optimal solution of the OPF problem. The proposed algorithm was used to minimize fuel cost, emissions, and transmission losses, which are considered to be parts of the objective function. The standard IEEE 30-bus and 57-bus systems were employed to investigate the performance of the proposed algorithm to find the optimal settings of the control variables. In order to investigate the single-objective and multiobjective optimizations, the simulation was divided into two cases. First, the proposed algorithm was used to solve a single-objective function. The results from the proposed algorithm show its superiority over other optimization algorithms in the literature. For the other case, the DA-PSO was successfully employed to solve the MO-OPF problem because the Pareto fronts generated by DA-PSO are better than those obtained by the original DA and PSO algorithms. All simulation results support the applicability, potential, and effectiveness of the proposed algorithm. However, the computation time of the DA-PSO is much slower than other algorithms in the literature because of the sequential computation of DA and PSO.

Author Contributions

Conceptualization, S.K. and R.C.; Methodology, S.K.; Software, S.K.; Validation, S.K., A.S., R.C., and N.R.W.; Formal Analysis, S.K.; Investigation, S.K.; Resources, A.S.; Data Curation, S.K.; Writing-Original Draft Preparation, S.K.; Writing-Review & Editing, S.K., R.C., N.R.W. and S.P.; Visualization, S.K.; Supervision, A.S.; Project Administration, R.C. and A.S.; Funding Acquisition, A.S.

Funding

This research was funded by the Thailand Research Fund through the Royal Golden Jubilee Ph.D. Program (Grant no. PHD/0192/2557) to Sirote Khunkitti and Apirat Siritaratiwat.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ABCartificial bee colony
ACOant colony optimization
APSOadaptive particle swarm optimization strategy
ARCBBOadaptive real coded biogeography-based optimization
DAdragonfly algorithm
DEdifferential evolution
DE-PSdifferential evolution pattern search
DSAdifferential search algorithm
EADDEevolving ant direction differential evolution
EEAefficient evolutionary algorithm
EGAenhanced genetic algorithm
EGA-DQLFenhanced genetic algorithm with decoupled quadratic load flow
EPevolutionary programming
EP-OPFevolutionary-programming-based optimal power flow
Fuzzy-GAfuzzy genetic algorithm
GAgenetic algorithm
GSAgravitational search algorithm
GWOgrey wolf optimizer
HBMOhoney bee mating optimization
HMPSO-SFLAhybrid modified particle swarm optimization-shuffle frog leaping algorithms
IEPimproved evolutionary programming
IPSOimproved particle swarm optimization
ISPEA2improved strength Pareto evolutionary algorithm
LDI-PSOparticle swarm optimization algorithm with linearly decreasing inertia weight
LTLBOLevy mutation teaching–learning-based optimization
MDE-OPFmodified differential evolution optimal power flow
MGBICAmodified Gaussian bare-bones multiobjective imperialist competitive algorithm
MICA-TLAhybrid modified imperialist competitive algorithm and teaching–learning algorithm
MO-DEAmultiobjective differential evolution algorithm
MOHSmultiobjective harmony search
MOMICAmultiobjective modified imperialist competitive algorithm
MO-OPFmultiobjective optimal power flow
MSFLAmodified shuffle frog leaping algorithm
MTLBOmodified teaching–learning-based optimization
NPSOnew particle swarm optimization
OPFoptimal power flow
P-OPFprobabilistic optimal power flow
PSOparticle swarm optimization
SGAstochastic genetic algorithm
TStabu search

References

  1. Sanseverino, E.R.; Buono, L.; Di Silvestre, M.L.; Zizzo, G.; Ippolito, M.G.; Favuzza, S.; Quynh, T.T.T.; Ninh, N.Q. A distributed minimum losses optimal power flow for islanded microgrids. Electr. Power Syst. Res. 2017, 152, 271–283. [Google Scholar] [CrossRef]
  2. Christakou, K.; Tomozei, D.C.; Le Boudec, J.Y.; Paolone, M. AC OPF in radial distribution networks—Part I: On the limits of the branch flow convexification and the alternating direction method of multipliers. Electr. Power Syst. Res. 2017, 143, 438–450. [Google Scholar] [CrossRef]
  3. Roy, P.K.; Ghoshal, S.P.; Thakur, S.S. Biogeography based optimization for multi-constraint optimal power flow with emission and non-smooth cost function. Expert Syst. Appl. 2010, 37, 8221–8228. [Google Scholar] [CrossRef]
  4. Ma, H.; Hart, J.L. Economic Dispatch in View of the Clean Air Act of 1990. IEEE Trans. Power Syst. 2000, 9, 972–978. [Google Scholar]
  5. Momoh, J.A.; El-Hawary, M.E.; Adapa, R. A review of selected optimal power flow literature to 1993 part I: Nonlinear and quadratic Programming Approaches. IEEE Trans. Power Syst. 1999, 14, 96–103. [Google Scholar] [CrossRef]
  6. Burchett, R.C.; Happ, H.H.; Vierath, D.R. Quadratically Convergent Optimal Power Flow. IEEE Trans. Power Appar. Syst. 1984, PAS-103, 3267–3275. [Google Scholar] [CrossRef]
  7. Yan, X.; Quintana, V.H. Improving an interior-point-based off by dynamic adjustments of step sizes and tolerances. IEEE Trans. Power Syst. 1999, 14, 709–716. [Google Scholar] [CrossRef]
  8. Ghasemi, M.; Ghavidel, S.; Ghanbarian, M.M.; Habibi, A. A new hybrid algorithm for optimal reactive power dispatch problem with discrete and continuous control variables. Appl. Soft Comput. 2014, 22, 126–140. [Google Scholar] [CrossRef]
  9. Ghasemi, M.; Ghanbarian, M.M.; Ghavidel, S.; Rahmani, S.; Mahboubi Moghaddam, E. Modified teaching learning algorithm and double differential evolution algorithm for optimal reactive power dispatch problem: A comparative study. Inf. Sci. (N. Y.) 2014, 278, 231–249. [Google Scholar] [CrossRef]
  10. Lai, L.L.; Ma, J.T.; Yokoyama, R.; Zhao, M. Improved genetic algorithms for optimal power flow under both normal and contingent operation states. Int. J. Electr. Power Energy Syst. 1997, 19, 287–292. [Google Scholar] [CrossRef]
  11. Abido, M.A. Optimal Power Flow Using Tabu Search Algorithm. Electr. Power Compon. Syst. 2002, 30, 469–483. [Google Scholar] [CrossRef]
  12. El Ela, A.A.A.; Abido, M.A.; Spea, S.R. Optimal power flow using differential evolution algorithm. Electr. Power Syst. Res. 2010, 80, 878–885. [Google Scholar] [CrossRef]
  13. Yuryevich, J.; Wong, K.P. Evolutionary Programming Based Optimal Power Flow Algorithm. IEEE Trans. Power Syst. 1999, 14, 1245–1250. [Google Scholar] [CrossRef]
  14. Deng, X.; He, J.; Zhang, P. A novel probabilistic optimal power flow method to handle large fluctuations of stochastic variables. Energies 2017, 10, 1623. [Google Scholar] [CrossRef]
  15. Wu, X.; Zhou, Z.; Liu, G.; Qi, W.; Xie, Z. Preventive security-constrained optimal power flow considering UPFC control modes. Energies 2017, 10, 1199. [Google Scholar] [CrossRef]
  16. Sliman, L.; Bouktir, T. Economic Power Dispatch of Power System with Pollution Control Using Multiobjective Ant Colony Optimization. Int. J. Comput. Intell. Res. 2007, 3, 145–153. [Google Scholar] [CrossRef]
  17. El-Fergany, A.A.; Hasanien, H.M. Single and Multi-objective Optimal Power Flow Using Grey Wolf Optimizer and Differential Evolution Algorithms. Electr. Power Compon. Syst. 2015, 43, 1548–1559. [Google Scholar] [CrossRef]
  18. Bai, W.; Lee, D.; Lee, K. Stochastic Dynamic AC Optimal Power Flow Based on a Multivariate Short-Term Wind Power Scenario Forecasting Model. Energies 2017, 10, 2138. [Google Scholar] [CrossRef]
  19. Abido, M.A. Optimal power flow using particle swarm optimization. Int. J. Electr. Power Energy Syst. 2002, 24, 563–571. [Google Scholar] [CrossRef]
  20. Bashishtha, T.K. Nature Inspired Meta-heuristic dragonfly Algorithms for Solving Optimal Power Flow Problem. Int. J. Electron. Electr. Comput. Syst. 2016, 5, 111–120. [Google Scholar]
  21. Shang, R.; Wang, Y.; Wang, J.; Jiao, L.; Wang, S.; Qi, L. A multi-population cooperative coevolutionary algorithm for multi objective capacitated arc routing problem. Inf. Sci. (N. Y.) 2014, 277, 609–642. [Google Scholar] [CrossRef]
  22. Yuan, X.; Zhang, B.; Wang, P.; Liang, J.; Yuan, Y.; Huang, Y.; Lei, X. Multi-objective optimal power flow based on improved strength Pareto evolutionary algorithm. Energy 2017, 122, 70–82. [Google Scholar] [CrossRef]
  23. Narimani, M.R.; Azizipanah-Abarghooee, R.; Zoghdar-Moghadam-Shahrekohne, B.; Gholami, K. A novel approach to multi objective optimal power flow by a new hybrid optimization algorithm considering generator constraints and multi-fuel type. Energy 2013, 49, 119–136. [Google Scholar] [CrossRef]
  24. Shabanpour-Haghighi, A.; Seifi, A.; Niknam, T. A modified teaching-learning based optimization for multi objective optimal power flow problem. Energy Convers. Manag. 2014, 77, 597–607. [Google Scholar] [CrossRef]
  25. Ghasemi, M.; Ghavidel, S.; Ghanbarian, M.M.; Gharibzadeh, M.; Azizi Vahed, A. Multi-objective optimal power flow considering the cost, emission, voltage deviation and power losses using multi objective modified imperialist competitive algorithm. Energy 2014, 78, 276–289. [Google Scholar] [CrossRef]
  26. Abaci, K.; Yamacli, V. Differential search algorithm for solving multi objective optimal power flow problem. Int. J. Electr. Power Energy Syst. 2016, 79, 1–10. [Google Scholar] [CrossRef]
  27. Niknam, T.; rasoul Narimani, M.; Jabbari, M.; Malekpour, A.R. A modified shuffle frog leaping algorithm for multi objective optimal power flow. Energy 2011, 36, 6420–6432. [Google Scholar] [CrossRef]
  28. Ghasemi, M.; Ghavidel, S.; Ghanbarian, M.M.; Gitizadeh, M. Multi-objective optimal electric power planning in the power system using Gaussian bare-bones imperialist competitive algorithm. Inf. Sci. (N. Y.) 2015, 294, 286–304. [Google Scholar] [CrossRef]
  29. Sivasubramani, S.; Swarup, K.S. Multi-objective harmony search algorithm for optimal power flow problem. Int. J. Electr. Power Energy Syst. 2011, 33, 745–752. [Google Scholar] [CrossRef]
  30. Ramesh Kumar, A.; Premalatha, L. Optimal power flow for a deregulated power system using adaptive real coded biogeography-based optimization. Int. J. Electr. Power Energy Syst. 2015, 73, 393–399. [Google Scholar] [CrossRef]
  31. El-Sehiemy, R.A.; Shaheen, A.M.; Farrag, S.M. Solving multi objective optimal power flow problem via forced initialised differential evolution algorithm. IET Gener. Transm. Distrib. 2016, 10, 1634–1647. [Google Scholar] [CrossRef]
  32. Ghasemi, M.; Ghavidel, S.; Rahmani, S.; Roosta, A.; Falah, H. A novel hybrid algorithm of imperialist competitive algorithm and teaching learning algorithm for optimal power flow problem with non-smooth cost functions. Eng. Appl. Artif. Intell. 2014, 29, 54–69. [Google Scholar] [CrossRef]
  33. Mirjalili, S. Dragonfly algorithm: A new meta-heuristic optimization technique for solving single-objective, discrete, and multi objective problems. Neural Comput. Appl. 2016, 27, 1053–1073. [Google Scholar] [CrossRef]
  34. Eberhart, R.; Kennedy, J. A New Optimizer Using Particle Swarm Theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science (MHS’95), Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar]
  35. Naderi, E.; Narimani, H.; Fathi, M.; Narimani, M.R. A novel fuzzy adaptive configuration of particle swarm optimization to solve large-scale optimal reactive power dispatch. Appl. Soft Comput. 2017, 53, 441–456. [Google Scholar] [CrossRef]
  36. Zhou, Y.; Wu, J.; Ji, L.; Yu, Z.; Lin, K.; Hao, L. Transient stability preventive control of power systems using chaotic particle swarm optimization combined with two-stage support vector machine. Electr. Power Syst. Res. 2018, 155, 111–120. [Google Scholar] [CrossRef]
  37. Mehdinejad, M.; Mohammadi-Ivatloo, B.; Dadashzadeh-Bonab, R.; Zare, K. Solution of optimal reactive power dispatch of power systems using hybrid particle swarm optimization and imperialist competitive algorithms. Int. J. Electr. Power Energy Syst. 2016, 83, 104–116. [Google Scholar] [CrossRef]
  38. Shihabudheen, K.V.; Mahesh, M.; Pillai, G.N. Particle swarm optimization based extreme learning neuro-fuzzy system for regression and classification. Expert Syst. Appl. 2018, 92, 474–484. [Google Scholar] [CrossRef]
  39. The University of Washington Electrical Engineering. Power System Test Case Archive, the IEEE 30-Bus Test System Data. Available online: https://www2.ee.washington.edu/research/pstca/pf30/pg_tca30bus.htm (accessed on 9 November 2017).
  40. Ongsakul, W.; Tantimaporn, T. Optimal Power Flow by Improved Evolutionary Programming. Electr. Power Compon. Syst. 2006, 34, 79–95. [Google Scholar] [CrossRef]
  41. Sayah, S.; Zehar, K. Modified differential evolution algorithm for optimal power flow with non-smooth cost functions. Energy Convers. Manag. 2008, 49, 3036–3042. [Google Scholar] [CrossRef]
  42. Bouktir, T.; Slimani, L.; Mahdad, B. Optimal Power Dispatch for Large Scale Power System Using Stochastic Search Algorithms. Int. J. Power Energy Syst. 2008, 28. [Google Scholar] [CrossRef]
  43. Sood, Y. Evolutionary programming based optimal power flow and its validation for deregulated power system analysis. Int. J. Electr. Power Energy Syst. 2007, 29, 65–75. [Google Scholar] [CrossRef]
  44. Niknam, T.; Narimani, M.R.; Aghaei, J.; Tabatabaei, S.; Nayeripour, M. Modified Honey Bee Mating Optimisation to solve dynamic optimal power flow considering generator constraints. IET Gener. Transm. Distrib. 2011, 5, 989. [Google Scholar] [CrossRef]
  45. Niknam, T.; Narimani, M.R.; Aghaei, J.; Azizipanah-Abarghooee, R. Improved particle swarm optimisation for multi objective optimal power flow considering the cost, loss, emission and voltage stability index. IET Gener. Transm. Distrib. 2012, 6, 515–527. [Google Scholar] [CrossRef]
  46. Kumari, M.S.; Maheswarapu, S. Enhanced Genetic Algorithm based computation technique for multi objective Optimal Power Flow solution. Int. J. Electr. Power Energy Syst. 2010, 32, 736–742. [Google Scholar] [CrossRef]
  47. Surender Reddy, S.; Bijwe, P.R.; Abhyankar, A.R. Faster evolutionary algorithm based optimal power flow using incremental variables. Int. J. Electr. Power Energy Syst. 2014, 54, 198–210. [Google Scholar] [CrossRef]
  48. The University of Washington Electrical Engineering. Power System Test Case Archive, the IEEE 57-Bus Test System Data. Available online: https://www2.ee.washington.edu/research/pstca/pf57/pg_tca57bus.htm (accessed on 9 November 2017).
  49. Ghasemi, M.; Ghavidel, S.; Gitizadeh, M.; Akbari, E. An improved teaching-learning-based optimization algorithm using Lévy mutation strategy for non-smooth optimal power flow. Int. J. Electr. Power Energy Syst. 2015, 65, 375–384. [Google Scholar] [CrossRef]
  50. Niknam, T. A new particle swarm optimization for non-convex economic dispatch. Trans. Electr. 2011, 21, 656–679. [Google Scholar] [CrossRef]
  51. Hsiao, Y.T.; Chen, C.H.; Chien, C.C. Optimal capacitor placement in distribution systems using a combination fuzzy-GA method. Int. J. Electr. Power Energy Syst. 2004, 26, 501–508. [Google Scholar] [CrossRef]
  52. Gitizadeh, M.; Ghavidel, S.; Aghaei, J. Using SVC to Economically Improve Transient Stability in Long Transmission Lines. IETE J. Res. 2014, 60, 319–327. [Google Scholar] [CrossRef]
  53. Rezaei Adaryani, M.; Karami, A. Artificial bee colony algorithm for solving multi objective optimal power flow problem. Int. J. Electr. Power Energy Syst. 2013, 53, 219–230. [Google Scholar] [CrossRef]
  54. Vaisakh, K.; Srinivas, L.R. Evolving ant direction differential evolution for OPF with non-smooth cost functions. Eng. Appl. Artif. Intell. 2011, 24, 426–436. [Google Scholar] [CrossRef]
  55. Duman, S.; Güvenç, U.; Sönmez, Y.; Yörükeren, N. Optimal power flow using gravitational search algorithm. Energy Convers. Manag. 2012, 59, 86–95. [Google Scholar] [CrossRef]
  56. Mahdad, B.; Srairi, K. Hierarchical adaptive PSO for multi objective OPF considering emissions based shunt FACTS. In Proceedings of the 38th Annual Conference on IEEE Industrial Electronics Society IECON 2012), Montreal, QC, Canada, 25–28 October 2012; pp. 1337–1343. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the dragonfly algorithm and particle swarm optimization (DA-PSO) algorithm for solving the multiobjective optimal power flow (MO-OPF) problem.
Figure 1. Flowchart of the dragonfly algorithm and particle swarm optimization (DA-PSO) algorithm for solving the multiobjective optimal power flow (MO-OPF) problem.
Energies 11 02270 g001
Figure 2. Two-dimensional Pareto fronts when considering fuel cost and emissions as part of the objective function for the IEEE 30-bus system.
Figure 2. Two-dimensional Pareto fronts when considering fuel cost and emissions as part of the objective function for the IEEE 30-bus system.
Energies 11 02270 g002
Figure 3. Two-dimensional Pareto fronts when considering fuel cost and transmission losses as part of the objective function for the IEEE 30-bus system.
Figure 3. Two-dimensional Pareto fronts when considering fuel cost and transmission losses as part of the objective function for the IEEE 30-bus system.
Energies 11 02270 g003
Figure 4. Two-dimensional Pareto fronts when considering emissions and transmission losses as part of the objective function for the IEEE 30-bus system.
Figure 4. Two-dimensional Pareto fronts when considering emissions and transmission losses as part of the objective function for the IEEE 30-bus system.
Energies 11 02270 g004
Figure 5. Three-dimensional Pareto fronts when considering fuel cost, emissions, and transmission losses as part of the objective function for the IEEE 30-bus system shown in the different views: (a) front view; (b) side view; (c) back view; (d) side view.
Figure 5. Three-dimensional Pareto fronts when considering fuel cost, emissions, and transmission losses as part of the objective function for the IEEE 30-bus system shown in the different views: (a) front view; (b) side view; (c) back view; (d) side view.
Energies 11 02270 g005
Figure 6. Two-dimensional Pareto fronts when considering fuel cost and emissions as part of the objective function for the IEEE 57-bus system.
Figure 6. Two-dimensional Pareto fronts when considering fuel cost and emissions as part of the objective function for the IEEE 57-bus system.
Energies 11 02270 g006
Figure 7. Two-dimensional Pareto fronts when considering fuel cost and transmission losses as part of the objective function for the IEEE 57-bus system.
Figure 7. Two-dimensional Pareto fronts when considering fuel cost and transmission losses as part of the objective function for the IEEE 57-bus system.
Energies 11 02270 g007
Figure 8. Two-dimensional Pareto fronts when considering emissions and transmission losses as part of the objective function for the IEEE 57-bus system.
Figure 8. Two-dimensional Pareto fronts when considering emissions and transmission losses as part of the objective function for the IEEE 57-bus system.
Energies 11 02270 g008
Figure 9. Three-dimensional Pareto fronts when considering fuel cost, emissions, and transmission losses as part of the objective function for the IEEE 57-bus system.
Figure 9. Three-dimensional Pareto fronts when considering fuel cost, emissions, and transmission losses as part of the objective function for the IEEE 57-bus system.
Energies 11 02270 g009
Table 1. Comparison of the simulation results from particle swarm optimization (PSO), dragonfly algorithm (DA), and DA-PSO for IEEE 30-bus system.
Table 1. Comparison of the simulation results from particle swarm optimization (PSO), dragonfly algorithm (DA), and DA-PSO for IEEE 30-bus system.
VariablesBest Fuel CostBest EmissionBest PLoss
PSODADA-PSOPSODADA-PSOPSODADA-PSO
Pg1 (MW)176.2376176.5128176.186164.167864.340764.099751.697451.598751.5893
Pg2 (MW)48.843248.695548.831867.669267.538367.629580.000080.000080.0000
Pg3 (MW)21.518421.443121.511950.000050.000050.000050.000050.000050.0000
Pg4 (MW)22.125722.099522.073735.000035.000035.000035.000035.000035.0000
Pg5 (MW)12.200012.067312.200530.000030.000030.000030.000040.000030.0000
Pg6 (MW)12.000012.009112.000040.000040.000040.000040.000040.000040.0000
Vg1 (p.u.)1.05001.05001.05001.05001.05001.05001.05001.05001.0500
Vg2 (p.u.)1.03811.03791.03791.04591.04721.04591.04771.04761.0476
Vg3 (p.u.)1.01101.01171.01091.02741.03091.02771.02921.02831.0292
Vg4 (p.u.)1.01941.01971.01871.03531.03771.03501.03661.03421.0363
Vg5 (p.u.)1.10001.10001.10001.10001.10001.10001.10001.03871.1000
Vg6 (p.u.)1.09991.08421.08281.08521.01401.07131.08501.06061.0712
T6-9 (p.u.)0.99731.03181.01661.01361.00171.04901.01531.10001.0482
T6-10 (p.u.)0.90000.90040.92100.90001.10000.90000.90000.90000.9000
T4-12 (p.u.)1.01570.99950.99801.00971.00030.99541.01050.97400.9962
T27-28 (p.u.)0.94030.95010.94780.95181.01360.96090.95290.96470.9618
Qc10 (MVar)28.64307.021910.05210.00000.00305.71747.075330.00005.2416
Qc24 (MVar)0.000011.097410.643330.000016.660510.633317.00859.953410.6499
Fuel Cost ($/h)802.5449802.1299802.1241945.0484944.9387944.7159968.1335967.8979967.8756
Emission (ton/h)0.3636190.3644110.3634940.2048860.2048610.2048530.2072940.2072800.207279
PLoss (MW)9.52499.42729.40413.43703.47903.32923.29743.19873.1893
Table 2. Comparison of the results from DA-PSO and other algorithms when considering only fuel cost as part of the objective function for IEEE 30-bus system.
Table 2. Comparison of the results from DA-PSO and other algorithms when considering only fuel cost as part of the objective function for IEEE 30-bus system.
AlgorithmsPg1 (MW)Pg2 (MW)Pg3 (MW)Pg4 (MW)Pg5 (MW)Pg6 (MW)Emission (ton/h)Loss (MW)Cost ($/h)Time (s)
TS [11]176.040048.760021.560022.050012.440012.00000.3630049.4500802.2900-
EP [13]173.848049.998021.386022.630012.928012.00000.3572179.3900802.620051.40
ACO [16]181.945047.001021.459621.446013.207012.01340.3820009.8520802.5780-
SFLA [27]179.033749.258020.318321.326911.542011.66550.3720009.7444802.5092-
MSFLA [27]179.192948.980420.451720.926411.589711.95790.3723009.6991802.2870-
IEP [40]176.235849.009321.502321.811512.338712.01290.36361010.8700802.465099.01
MDE-OPF [41]175.974048.884021.510022.240012.251012.00000.3629009.4590802.376023.25
SGA [42]179.367044.240024.610019.900010.710014.09000.3711299.5177803.6990-
EP-OPF [43]175.029748.952221.420022.702012.904012.10350.3601259.7114803.5710-
HBMO [44]178.464646.274021.459621.446013.207012.01340.3692129.4662802.211028.56
PSO176.237648.843221.518422.125712.200012.00000.3636199.5249802.544992.18
DA176.512848.695521.443122.099512.067312.00910.3644119.4272802.1299103.06
DA-PSO176.186148.831821.511922.073712.200512.00000.3634949.4041802.1241287.13
Table 3. Comparison of the results from DA-PSO and other algorithms when considering only emissions as part of the objective function for IEEE 30-bus system.
Table 3. Comparison of the results from DA-PSO and other algorithms when considering only emissions as part of the objective function for IEEE 30-bus system.
AlgorithmsPg1 (MW)Pg2 (MW)Pg3 (MW)Pg4 (MW)Pg5 (MW)Pg6 (MW)Cost ($/h)Loss (MW)Emission (ton/h)Time (s)
ACO [16]64.372072.160449.543832.909928.611339.7390945.58703.93680.221000-
HMPSO-SFLA [23]64.814868.069250.000034.999930.000040.0000948.30524.48390.205200-
TLBO [24]63.522168.734549.993134.989429.982439.9801947.43923.80160.205030-
MTLBO [24]64.292467.625050.000035.000030.000040.0000945.19653.51740.204930-
DSA [26]64.072567.571150.000035.000030.000040.0000944.40863.24370.205826-
MSFLA [27]65.779868.268850.000034.999929.998239.9970951.51065.64370.205600-
SFLA [27]64.484071.380749.857335.000030.000039.9729960.19117.29490.206300-
GA [27]78.288568.160246.784833.490930.000036.3713936.61529.69570.211700-
GBICA [28]64.312567.493850.000035.000029.992440.0000944.65163.39870.204900-
IPSO [45]67.040068.140050.000035.000030.000040.0000954.24805.36200.205800-
PSO64.167867.669250.000035.000030.000040.0000945.04843.43700.20488691.84
DA64.066767.689750.000035.000030.000040.0000944.88193.35640.204861103.20
DA-PSO64.099767.629550.000035.000030.000040.0000944.71593.32920.204853290.01
Table 4. Comparison of the results from DA-PSO and other algorithms when considering only losses as part of the objective function for IEEE 30-bus system.
Table 4. Comparison of the results from DA-PSO and other algorithms when considering only losses as part of the objective function for IEEE 30-bus system.
AlgorithmsPg1 (MW)Pg2 (MW)Pg3 (MW)Pg4 (MW)Pg5 (MW)Pg6 (MW)Cost ($/h)Emission (ton/h)Loss (MW)Time (s)
GWO [17]51.810080.000050.000035.000030.000040.0000968.38000.2073103.410015.90
DE [17]51.820079.990049.990035.000029.980040.0000968.23000.2073113.380016.50
MOHS [29]66.275979.641346.883534.888029.121330.0558928.50990.2128903.5165-
EGA-DQLF [46]51.600880.000050.000035.000030.000040.0000967.86000.2072813.2008-
EEA [47]59.321674.813249.854734.908428.109939.7538952.37850.2067353.28235.72
EGA [47]51.674079.970050.000035.000030.000040.0000967.93000.2072753.244029.71
PSO51.697480.000050.000035.000030.000040.0000968.13350.2072943.297493.36
DA51.594180.000050.000035.000040.000040.0000967.88690.2072803.1941102.81
DA-PSO51.589380.000050.000035.000030.000040.0000967.87560.2072793.1893292.33
Table 5. Comparison of the simulation results from PSO, DA, and DA-PSO for IEEE 57-bus system.
Table 5. Comparison of the simulation results from PSO, DA, and DA-PSO for IEEE 57-bus system.
VariablesBest Fuel CostBest EmissionBest PLoss
PSODADA-PSOPSODADA-PSOPSODADA-PSO
Pg1 (MW)142.7472154.8513141.4617236.4846246.6610236.4531193.1342269.9574202.6688
Pg2 (MW)88.842776.622787.7806100.000044.6053100.00008.85810.20470.0000
Pg3 (MW)44.902549.444044.6638139.9999140.0000140.0000139.973160.2481140.0000
Pg4 (MW)70.8490100.000073.6254100.000078.0610100.0000100.000055.3529100.0000
Pg5 (MW)458.6003438.7375458.9904292.5686329.7090292.1457309.5411377.9311308.2507
Pg6 (MW)100.0000100.000097.4933100.000082.5653100.0000100.000090.2309100.0000
Pg7 (MW)360.3487347.6947361.7228298.6306344.7194298.4568410.0000410.0000410.0000
Vg1 (p.u.)1.05001.05001.05001.05001.05001.05001.05001.05001.0500
Vg2 (p.u.)1.04941.04531.04881.05131.04861.05061.04581.03511.0450
Vg3 (p.u.)1.04791.04751.04551.05321.05501.05051.05281.03421.0520
Vg4 (p.u.)1.06281.07551.05811.05181.03031.04931.05371.04541.0525
Vg5 (p.u.)1.07921.08021.07451.05511.01421.05061.06031.05631.0566
Vg6 (p.u.)1.04551.05681.04421.02661.01331.02671.03491.02401.0348
Vg7 (p.u.)1.04101.06291.03941.02511.05561.02621.03921.01191.0384
T4–8 (p.u.)0.94291.10001.02210.96291.10000.96910.96250.97630.9730
T4–18 (p.u.)0.99161.01400.99530.97641.06820.98700.98651.02771.0275
T21–20 (p.u.)1.01511.10001.01961.02331.01131.02281.02261.02861.0442
T24–25 (p.u.)0.90000.98571.02120.90821.10001.10000.91111.04301.0140
T24–25 (p.u.)0.93780.98120.98960.90000.96450.96090.91180.99481.0145
T24–26 (p.u.)1.02191.10001.01751.01150.99811.00861.01301.03121.0111
T7–29 (p.u.)0.99011.02530.99830.97910.96580.99040.98260.98160.9945
T34–32 (p.u.)0.92771.07310.95820.92851.00690.96820.92170.98150.9577
T11–41 (p.u.)0.90000.98790.90630.90000.90620.90000.90001.10000.9036
T15–45 (p.u.)0.96670.97700.97140.97501.01780.97860.97790.97610.9807
T14–46 (p.u.)0.95780.98070.96160.96360.97430.95690.95940.93730.9616
T10–51 (p.u.)0.97480.98990.97660.96420.98240.96740.97100.93950.9707
T13–49 (p.u.)0.93001.01530.93010.92570.98540.92740.92920.98500.9375
T11–43 (p.u.)0.97850.97380.97560.96120.97070.96470.97040.93730.9767
T40–56 (p.u.)0.99621.10001.01050.97150.94600.97100.99691.04570.9972
T39–57 (p.u.)0.96921.10000.96210.97281.07990.97420.96290.93730.9645
T9–55 (p.u.)0.98561.07320.99880.96760.99370.98100.97560.97470.9842
Qc18 (MVar)18.74506.675113.28040.000015.69772.449310.924712.97344.5146
Qc25 (MVar)13.86148.722012.63077.304222.116416.816928.743012.994914.5906
Qc53 (MVar)12.068622.201513.97250.02499.341012.555119.743210.714312.9250
Fuel Cost ($/h)41,698.3741,828.4541,674.6245,671.2245,449.1345,648.6744,951.8043,464.1745,039.05
Emission (ton/h)1.90271.68831.90871.08141.30971.07991.38211.75621.4014
PLoss (MW)15.490316.550214.938016.883715.521016.255610.707613.643010.1212
Table 6. Comparison of the results from DA-PSO and other algorithms when considering only fuel cost as part of the objective function for the IEEE 57-bus system.
Table 6. Comparison of the results from DA-PSO and other algorithms when considering only fuel cost as part of the objective function for the IEEE 57-bus system.
AlgorithmsCost ($/h)
MTLBO [23]41,638.3822
DSA [25]41,686.8200
GBICA [27]41,740.2884
MGBICA [27]41,715.7101
ARCBBO [29]41,686.0000
MO-DEA [30]41,683.0000
MICA-TLA [31]41,675.0545
TLBO [48]41,695.6629
LTLBO [49]41,679.5451
NPSO [50]41,699.5163
Fuzzy-GA [51]41,716.2808
DE-PS [52]41,685.2950
ABC [53]41,693.9589
LDI-PSO [53]41,815.5035
EADDE [54]41,713.6200
GSA [55]41,695.8717
APSO [56]41,713.8868
PSO41,698.3672
DA41,828.4473
DA-PSO41,674.6209
Table 7. Comparison of the results from DA-PSO and other algorithms when considering the only emissions as part of the objective function for the IEEE 57-bus system.
Table 7. Comparison of the results from DA-PSO and other algorithms when considering the only emissions as part of the objective function for the IEEE 57-bus system.
AlgorithmsEmission (ton/h)
GBICA [27]1.1881
MGBICA [27]1.1724
PSO1.0814
DA1.3097
DA-PSO1.0799
Table 8. Comparison of the results from DA-PSO and its traditional algorithms when considering only transmission losses as part of the objective function for the IEEE 57-bus system.
Table 8. Comparison of the results from DA-PSO and its traditional algorithms when considering only transmission losses as part of the objective function for the IEEE 57-bus system.
AlgorithmsLoss (MW)
PSO10.7076
DA13.6430
DA-PSO10.1212

Share and Cite

MDPI and ACS Style

Khunkitti, S.; Siritaratiwat, A.; Premrudeepreechacharn, S.; Chatthaworn, R.; Watson, N.R. A Hybrid DA-PSO Optimization Algorithm for Multiobjective Optimal Power Flow Problems. Energies 2018, 11, 2270. https://doi.org/10.3390/en11092270

AMA Style

Khunkitti S, Siritaratiwat A, Premrudeepreechacharn S, Chatthaworn R, Watson NR. A Hybrid DA-PSO Optimization Algorithm for Multiobjective Optimal Power Flow Problems. Energies. 2018; 11(9):2270. https://doi.org/10.3390/en11092270

Chicago/Turabian Style

Khunkitti, Sirote, Apirat Siritaratiwat, Suttichai Premrudeepreechacharn, Rongrit Chatthaworn, and Neville R. Watson. 2018. "A Hybrid DA-PSO Optimization Algorithm for Multiobjective Optimal Power Flow Problems" Energies 11, no. 9: 2270. https://doi.org/10.3390/en11092270

APA Style

Khunkitti, S., Siritaratiwat, A., Premrudeepreechacharn, S., Chatthaworn, R., & Watson, N. R. (2018). A Hybrid DA-PSO Optimization Algorithm for Multiobjective Optimal Power Flow Problems. Energies, 11(9), 2270. https://doi.org/10.3390/en11092270

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop