Artificial Bee Colony Algorithm Applied To Optimal Power Flow Solution Incorporating Stochastic Wind Power
Artificial Bee Colony Algorithm Applied To Optimal Power Flow Solution Incorporating Stochastic Wind Power
Corresponding Author:
Yasir G. Rashid
Department of Electronic Engineering
College of Engineering, University of Diyala
Baqubah, Diyala, Iraq
Email: [email protected]
1. INTRODUCTION
The majority of the world's fossil-fuel power generation operations use coal and natural gas to
generate electricity, which is one of the most expensive commodities used to generate electric power.
Polluting emissions from electricity generation based on the combustion of fossil fuels account for a sizable
portion of global greenhouse gas emissions [1], [2]. As a result of economic and environmental reasons,
workers in the field of electric energy were encouraged to increase and develop renewable energy. The
electrical power control are experiencing noteworthy changes due to an increase in wind energy penetration
level, causing unused challenges to system operation and planning [3], [4]. Therefore, the operators of power
systems both in the planning and operating stage are very interested in optimal power flow (OPF) [5]. The
main objective of an optimal power flow methodology is to find the ideal working of a power system by
optimizing a specific objective whereas fulfilling certain indicated physical and security limitations [6], [7].
In recent years, the rapid development of computational intelligence have motivated researchers in
the field of optimization algorithms to resolve various complex optimization cases such as particle swarm
optimization algorithm (PSO) [8], [9], improved colliding bodies optimization method [10], imperialist
competitive method [11], black-hole-based optimization technique [12], differential evolutionary technique
[13], hybrid algorithm of PSO and GSA algorithms [14], gravitational search method (GSM) [15], [16],
improved PSO algorithm [17], biogeography-based optimization technique [18], chaotic self-adaptive
differential harmony search method [19], grey wolf optimizer [20], fuzzy-based hybrid PSO algorithm [21],
differential search technique [22], multiphase search optimization technique [23], harmony search technique
[24], Jaya optimization technique [25], artificial bee colony (ABC) algorithm [26], differential evolution
(DE) [27], biogeography-based optimization (BBO) [28], teaching-learning-optimization algorithm [29], and
the firefly algorithm (FA) [30]. This paper was motivated by two factors. First, the application of the artificial
bee colony algorithm to solve the optimum power flow problem has been studied. Second, solving OPF
considering wind power penetration of different sites (single & multiple) and studying the impact of the wind
power penetration on the slack bus generation, the total production cost, active power losses and voltage
deviation.
(1)
Subject to
(2)
(3)
Where F refers to the target (objective) function to be minimized, x and u are state and control variables
respectively. The state vector x including; i) PG , generating power at swing (slack) bus, ii) QG, reactive
1
generating power outputs, and iii) VL, load bus voltage. x can be written as:
(4)
Where NG, NL, NTL and SL are the number of generator buses, number of load buses, transmission lines
and number of transmission line loading, respectively. The control vector u including; i) PG, generator active
power outputs, ii) VG, generator voltages, iii) QC, shunt VAR compensations, and iv) T transformer tap
settings. u can be written as:
(5)
Where NC and NT are the shunt VAR compensators output and the transformers regulated number,
respectively [31].
∑ (6)
Where , , cost coefficients of fuel generators , N number of generation units, active power
generation of generators .
∑ (7)
Artificial bee colony algorithm applied to optimal power flow solution … (Vian H. Ahgajan)
1892 ISSN: 2088-8694
∑ (8)
Where wj is active power generated by jth wind turbine and qj is direct cost coefficient.
∑ (9)
= + + ( + wj )
(10)
+( ) (1+ ,( ) (1+ ,( )
and
∑ (11)
= ( + wj )
(12)
+( ) (1+ ,( ) (1+ ,( )
Where Cpwj and Crwj are the overestimation and underestimation cost coefficient of jth wind
generator in $/MW h respectively. ( ) and are the overestimation and underestimation
anticipated value of wind power for jth wind turbine. kj and cj are a shape factor and a scale of the jth wind
generator respectively estimating of wind speed in the Weibull probability density function (pdf). v inj, vout,I,
vr,j are cut-in, cut-out and rated wind speed respectively. v1 = vin + (vr − vin) w1/wr is an intermediary
parameter in [6]. Minimize the total production cost in wind-thermal power system can be expressed as [33]:
(13)
∑ ( ) (14)
Where m is the total number of lines in the system, G k is the conductance of the kth line, Vj and Vi are the
voltage magnitude at bus j and bus i respectively, δj and δi are the voltage phase angle at bus j and i
respectively [34]. Second, minimize the voltage deviation (VD) of all load buses to improve the voltage
profile on load buses. The voltage deviation given by (15) [35]:
∑ (15)
Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1890 – 1899
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1893
members. The colony of artificial bees consists of three main groups, which are the employed bees, onlooker
bees and scout bees. This breed of bee features a distinctive part within the optimization preparation. The
employed bee can remember the location of the extra nectar as well as it chooses the best of the others to
drink from, while the onlooker bees use what the employed bees have collected to come up with a solution
for what nectar they can't remember. For an optimization problem, an algorithm consists of three steps is as
follows: In the first step, the employed bees are dispatched to find all the resources needed, and then the
nectar amount is calculated. Step two, the onlooker bees choose an asset that matches the information from
the already-discovered honeydew assets. The employed bumblebee was sent out to the fields to select new
locations in order to identify potential food sources. "Looking" bees would be further broken into two
categories: the "used" bees and the "observing" bees. The algorithm works on the basis that the number of
employable bees equals the number of available sources of nectar. When we understand where the issues
likely lie, we'll be better equipped to deal with them [36].
ABC algorithm:
a) Initialization phase
In the first step, variables ( = 1, 2, 3, … ) that have not been measured yet are selected at random,
using some sort of random methodology.
b) Employed bee phase
The new sources are identified by each employed bee whose amounts are equal to the half of the total
sources. a new source can be found by:
(16)
Where j is a randomly selected parameter index, is a random number between [0, 1] and it has to be
different from , is a random number within the range [-1, 1], is the current position of food source
which comparing two food postion visually by bee from this parameter the production of the neighbor
food source can be controlled. The new food source postion is produced and evaluated by the artificial
bee,by comparing the current food source with previous source taking its performance in the consider.
From the information that obtained if the new source has equal or better amount of food or nectar than the
old source,it used to replace the old source in the memory. Otherwise, the old source would be retained in
memory.
c) Onlooker bee phase
In this phase ,the onlooker bees are work on the principle of probability by selecting the food source with
probability can be written as:
= (17)
∑
Where and are the fitness value and probability associated with solution respectively. In each
colony, great responsibility for random research is scout bees’ bear.
d) Scout bee phase
In this stage, the scout bee randomly investigates food sources without direction from the queen. Every
scout in the swarm thinks that he or she is an explorer. If the supply of food decreases below the gainful
level or as a result of applying a given level of the food application of the nectar, the bees associated with
it cease feeding. When you have new information, a new understanding, or a new insight, the limit on the
number of bees tells you how many from the source and how many to the destination.
(18)
Where and are the maximum and minimum limits for optimization parameter, rand (0, 1) is
a random number within the range [0, 1]. The number of iterations in ABC algorithm considered as the
important criterion for stopping an ABC algorithm.
An optimization algorithm might therefore determine that the stopping criteria to be:
1. Number of maximum iterations
2. Maximum error between two consecutive iterations
Figure 1 shown the flowchart of the ABC algorithm based OPF problem.
Artificial bee colony algorithm applied to optimal power flow solution … (Vian H. Ahgajan)
1894 ISSN: 2088-8694
4. CASE STUDY
In this paper, two wind farms connecting to bus 10 and bus 24 are suggested. Figure 2 shown the
standard IEEE 30 system with two wind farms. The wind power penetration level is defined as the ratio of
the installed wind power capacity to the total-installed system generation capacity of 10%. The total power
generation of six thermal generating in system are around 400MW, therefore the installed wind power
capacity is 40 MW. Two wind farms included 10 wind turbines each one has rating 2 MW (Vestas V90, 2
MW) and connected at bus 10 and bus 24 (20 MW in each bus) is used to analyse the impact of incorporating
wind farm on different performance analysis of system. Several scenarios with dispersed wind penetration
levels from 0% to 100% have been investigated.
Figure 2. Single-line diagram of the IEEE 30-bus system including two wind farms at bus10 & bus 24
Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1890 – 1899
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1895
Table 2. Comparison of ABC with LP when wind farm with different wind power penetration levels
connected to system
Wind Bus no. ABC LP [33]
penetration
slack bus Cost($hr) Losses VD slack bus Cost($hr) Losses VD
(MW) (p.u.) (MW) (p.u.)
0% 10 177.05 800.638 8.9246 0.8977 176 802.46 10.31 0.8513
25% 171.38 783.142 9.152 0.604 175 791.03 9.12 0.594
50% 166.01 765.221 8.246 0.920 171 776.83 9.05 0.900
75% 160.65 747.941 7.888 0.931 166 764.63 8.79 0.912
100% 155.31 730.943 7.550 0.941 160 753.42 8.27 0.930
0% 24 177.05 800.638 8.9246 0.8977 176 802.46 10.31 0.8513
25% 171.28 782.447 8.522 0.924 177 790.89 9.03 0.900
50% 165.83 764.657 8.072 0.948 172 776 8.85 0.915
75% 160.43 746.814 7.674 0.972 168 764.63 8.34 0.960
100% 155.09 730.234 7.326 0.993 163 753.42 7.96 0.984
1.08
without WF bus 10
1.06
1.04
1.02
1
0.98
0.96
0.94
2 3 4 6 7 9 10 12 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Bus No.
Figure 3. Load bus voltage profile 30-bus IEEE system
Artificial bee colony algorithm applied to optimal power flow solution … (Vian H. Ahgajan)
1896 ISSN: 2088-8694
Iterations Iterations
Iterations Iterations
Figure 4. Convergence characteristics of the ABC for penetration levels of wind power at bus 10 only
Table 3. Comparison of ABC with LP when wind farm with different wind power penetration levels
connected to system
Wind Bus no. ABC LP [ 33]
penetration
slack Cost($hr) Losses VD slack bus Cost($hr) Losses VD
bus (MW) (p.u.) (MW) (p.u.)
0% 10 177.05 800.6380 8.9246 0.8977 176 802.46 10.31 0.8513
25% & 165.90 764.670 8.614 0.639 170 758.89 9.03 0.602
50% 24 155.12 730.526 8.502 0.643 159 736.91 8.85 0.609
75% 144.44 694.795 7.599 0.905 148 701.78 8.34 0.885
100% 133.86 661.621 7.565 0.813 138 680.59 7.96 0.803
Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1890 – 1899
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1897
1.08
without WF bus 10& 24
1.06
1.04
1.02
1
0.98
0.96
0.94
2 3 4 6 7 9 10 12 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Bus No.
Iterations Iterations
Total production cost ($/h)
Total production cost ($/h)
Iterations Iterations
Figure 6. Convergence characteristics of the ABC for penetration levels of wind power on bus 10 &
24 together
5. CONCLUSION
This paper proposes the application of artificial bee colony (ABC) algorithm optimal power flow for
a system that incorporates thermal units and wind farms during normal operation. The performance of the
ABC was applied to standard IEEE-30 bus system with and without incorporating wind farm to show its
impact on the the slack bus generation, the total production cost, active power losses and voltage deviation,
and compared its simulation results with another method. Based on technical results obtained are it can be
noticed that the ABC high performance than the rest methods, and concluded that an optimal integration and
location of wind farms give significant to system, such as reducing in the total production cost, active power
losses and improvement in the load bus voltage profile, while high performance can be noticed when a wind
farm site on bus 24 rather than its site on bus 10. Finally, the results are exceptionally much promising.
Artificial bee colony algorithm applied to optimal power flow solution … (Vian H. Ahgajan)
1898 ISSN: 2088-8694
REFERENCES
[1] M. Čepin, “Evaluation of the power system reliability if a nuclear power plant is replaced with wind power plants,”
Reliability Engineering & System Safety, vol. 185, pp. 455-464, 2019, doi: 10.1016/j.ress.2019.01.010.
[2] Y. G. Rashid, F. M. Tuaimah, and K. Mohammed, “Assessment of integrating wind energy system on iraqi power
grid capability limit,” Journal of Engineering and Applied Sciences, vol. 14, no.17, pp. 5942-5954, 2019, doi:
10.36478/jeasci.2019.5942.5954.
[3] D. Wanget al., “Optimal scheduling strategy of district integrated heat and power system with wind power and
multiple energy stations considering thermal inertia of buildings under different heating regulation modes,” Applied
Energy, vol. 240, pp. 341-358, 2019, doi: 10.1016/j.apenergy.2019.01.199.
[4] J. Liu, Q. Yao, and Y. Hu, “Model predictive control for load frequency of hybrid power system with wind power
and thermal power,” Energy, vol. 172, pp. 555-565, 2019, doi: 10.1016/j.energy.2019.01.071.
[5] P. P. Biswas, P. N. Suganthan, and G. A. J. Amaratunga, “Optimal power flow solutions incorporating stochastic
wind and solar power,” Energy Conversion and Management, vol. 148, pp. 1194-1207, 2017, doi:
10.1016/j.enconman.2017.06.071.
[6] R. Roy, and H. T. Jadhav, “Optimal power flow solution of power system incorporating stochastic wind power
using Gbest guided artificial bee colony algorithm,” International Journal of Electrical Power & Energy Systems,
vol. 64, pp. 562-578, 2015, doi: 10.1016/j.ijepes.2014.07.010.
[7] J. Luo, L. Shi, and Y. Ni, “A solution of optimal power flow incorporating wind generation and power grid
uncertainties,” IEEE Access, vol. 6, pp. 19681-19690, 2018, doi: 10.1109/ACCESS.2018.2823982.
[8] C. Shilaja, and T. Arunprasath, “Optimal power flow using moth swarm algorithm with gravitational search
algorithm considering wind power,” Future Generation Computer Systems, vol. 98, pp. 708-715, 2019, doi:
10.1016/j.future.2018.12.046.
[9] Y. Merzoug, B. Abdelkrim, and B. Larbi, “Optimal placement of wind turbine in a radial distribution network using
PSO method,” International Journal of Power Electronics and Drive Systems (IJPEDS), vol. 11, no. 2, pp. 1074-
1081, 2020, doi: 10.11591/ijpeds.v11.i2.pp1074-1081.
[10] H. R. E. H. Bouchekara, A. E. Chaib, M. A. Abido, and R. A. El-Sehiemy, “Optimal power flow using an improved
colliding bodies optimization algorithm,” Applied Soft Computing, vol. 42, pp. 119-131, 2016, doi:
10.1016/j.asoc.2016.01.041.
[11] A. J. Ghanizadeh, G. Mokhtari1, M. Abedi, and G. B. Gharehpetian, “Optimal power flow based on imperialist
competitive algorithm,” International Review of Electrical Engineering, vol. 6, no. 4, pp. 1847-1852, 2011.
[12] H. R. E. H. Bouchekara, “Optimal power flow using black-hole-based optimization approach,” Applied Soft
Computing, vol. 24, pp. 879-888, 2014, doi: 10.1016/j.asoc.2014.08.056.
[13] A. M. Shaheen, R. A. El-Sehiemy, and S. M. Farrag, “Solving multi-objective optimal power flow problem via
forced initialised differential evolution algorithm,” IET Generation, Transmission & Distribution,
vol. 10, no. 7, pp. 1634-1647, 2016, doi: 10.1049/iet-gtd.2015.0892.
[14] J. Radosavljević, D. Klimenta, M. Jevtić, and N. Arsić, “Optimal power flow using a hybrid optimization algorithm
of particle swarm optimization and gravitational search algorithm,” Electric Power Components and Systems, vol.
43, no. 17, pp. 1958-1970, 2015, doi: 10.1080/15325008.2015.1061620.
[15] S. Duman, U. Güvenç, Y. Sönmez, and N. Yörükeren, “Optimal power flow using gravitational search algorithm,”
Energy conversion and management, vol. 59, pp. 86-95, 2012, doi: 10.1016/j.enconman.2012.02.024.
[16] M. S. Jahan, and N. Amjady, “Solution of large-scale security constrained optimal power flow by a new bi-level
optimisation approach based on enhanced gravitational search algorithm,” IET Generation, Transmission &
Distribution, vol. 7, no. 12, pp. 1481-1491, 2013, doi: 10.1049/iet-gtd.2012.0697.
[17] Dieu Ngoc V, and Peter Schegner, “An improved particle swarm optimization for optimal power flow,” in Meta-
heuristics Optimization Algorithms in Engineering, Business, Economics, and Finance, IGI Global, 2013, pp. 1-40,
doi: 10.4018/978-1-4666-2086-5.ch001.
[18] P. Pravina, M. R. Babu, and A. R. Kumar, “Solving optimal power flow problems using adaptive quasi-
oppositional differential migrated biogeography-based optimization,” Journal of Electrical Engineering &
Technology, vol. 16, pp. 1891-1903, 2021, doi: 10.1007/s42835-021-00739-z.
[19] R. Arul, G. Ravi, and S. Velusami, “Solving optimal power flow problems using chaotic self-adaptive differential
harmony search algorithm,” Electric Power Components and Systems, vol. 41, no. 8, pp. 782-805, 2013, doi:
10.1080/15325008.2013.769033.
[20] A. Meng et al., “A high-performance crisscross search based grey wolf optimizer for solving optimal power flow
problem,” Energy, vol. 225, p. 120211, 2021, doi: 10.1016/j.energy.2021.120211.
[21] R.-H. Liang, S.-R. Tsai, Y.-T. Chen, and W.-T. Tseng, “Optimal power flow by a fuzzy based hybrid particle
swarm optimization approach,” Electric Power Systems Research, vol. 81, no. 7, pp. 1466-1474, 2011, doi:
10.1016/j.epsr.2011.02.011.
[22] H. R. E.-H. Bouchekara, and M. A. Abido, “Optimal power flow using differential search algorithm,” Electric
Power Components and Systems, vol. 42, no. 15, pp. 1683-1699, 2014, doi: 10.1080/15325008.2014.949912.
[23] H. Pulluri, R. Naresh, and V. Sharma, “A solution network based on stud krill herd algorithm for optimal power
flow problems,” Soft Computing, vol. 22, no. 1, pp. 159-176, 2018, doi: 10.1007/s00500-016-2319-3.
[24] S. S. Reddy, “Optimal power flow using hybrid differential evolution and harmony search algorithm,” International
Journal of Machine Learning and Cybernetics, vol. 10, no. 5, pp. 1077-1091, 2019, doi: 10.1007/s13042-018-0786-
9.
Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1890 – 1899
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1899
[25] W. Warid, H. Hizam, N. Mariun, and N. I. Abdul-Wahab, “Optimal power flow using the Jaya algorithm,”
Energies, vol. 9, no. 9, p. 678, 2016, doi: 10.3390/en9090678.
[26] X. He, W. Wang, J. Jiang, and L. Xu, “An improved artificial bee colony algorithm and its application to multi-
objective optimal power flow,” Energies, vol. 8, no. 4, pp. 2412-2437, 2015, doi: 10.3390/en8042412.
[27] P. P. Biswas, P. N. Suganthan, R. Mallipeddi, and G. A. J. Amaratunga, “Optimal power flow solutions using
differential evolution algorithm integrated with effective constraint handling techniques,” Engineering Applications
of Artificial Intelligence, vol. 68, pp. 81-100, 2018, doi: 10.1016/j.engappai.2017.10.019.
[28] A. A. Christy, and P. A. D. V. Raj, “Adaptive biogeography based predator–prey optimization technique for
optimal power flow,” International Journal of Electrical Power & Energy Systems, vol. 62, pp. 344-352, 2014, doi:
10.1016/j.ijepes.2014.04.054.
[29] S. S. Reddy, and B. K. Panigrahi, “Optimal power flow using clustered adaptive teaching learning-based
optimisation,” International Journal of Bio-Inspired Computation, vol. 9, no. 4, pp. 226-234, 2017, doi:
10.1504/IJBIC.2017.084316.
[30] P. Balachennaiah, M. Suryakalavathi, and P. Nagendra, “Firefly algorithm based solution to minimize the real
power loss in a power system,” Ain Shams Engineering Journal, vol. 9, no. 1, pp. 89-100, 2018, doi:
10.1016/j.asej.2015.10.005.
[31] Z. Ullah, S. Wang, J. Radosavljević, and J. Lai, “A Solution to the Optimal Power Flow Problem Considering WT
and PV Generation,” IEEE Access, vol. 7, pp. 46763-46772, 2019, doi: 10.1109/ACCESS.2019.2909561.
[32] V. H. Ahgajan, and F. M. Tuaimah, “Optimal power flow for a power system under particle swarm optimization
(PSO) based,” International Journal of Computer Applications, vol. 177, no. 33, pp. 0975-8887, 2020, doi:
10.5120/ijca2020919757.
[33] U. Güvenç, and E. Kaymaz, “Economic dispatch integrated wind power using coyote optimization algorithm,” in
7th International Istanbul Smart Grids and Cities Congress and Fair (ICSG), 2019, pp. 179-183, doi:
10.1109/SGCF.2019.8782354.
[34] Y. G. Rashid, and F. M. Tuaimah, “Impact of integrating wind energy into power grid system,” International
Journal of Computer Applications, vol. 177, no. 14, pp. 975-8887, 2019, doi: 10.5120/ijca2019919505.
[35] H. Xiao, Z. Dong, L. Kong, W. Pei, and Z. Zhao, “Optimal power flow using a novel metamodel based global
optimization method,” Energy Procedia, vol. 145, pp. 301-306, 2018, doi: 10.1016/j.egypro.2018.04.055.
[36] M. Ettappan, V. Vimala, S. Ramesh, and V. T. Kesavan, “Optimal reactive power dispatch for real power loss
minimization and voltage stability enhancement using Artificial Bee Colony Algorithm,” Microprocessors and
Microsystems, vol. 76, p. 103085, 2020, doi: 10.1016/j.micpro.2020.103085.
[37] V. H. Ahgajan, “Optimal power flow for grid using artificial bees algorithm (ABC),” M.Sc thesis, Razi University,
Iran, 2018.
Artificial bee colony algorithm applied to optimal power flow solution … (Vian H. Ahgajan)