SlideShare a Scribd company logo
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 34
Paper Publications
Dwindling of real power loss by using Improved
Bees Algorithm
K. Lenin, B. Ravindranath Reddy, and M. Surya Kalavathi
Jawaharlal Nehru Technological University Kukatpally, Hyderabad 500 085, India
Abstract: In this paper, a new Improved Bees Algorithm (IBA) is proposed for solving reactive power dispatch
problem. The aim of this paper is to utilize an optimization algorithm called the improved Bees Algorithm,
inspired from the natural foraging behaviour of honey bees, to solve the reactive power dispatch problem. The
IBA algorithm executes both an exploitative neighbourhood search combined with arbitrary explorative search.
The proposed Improved Imperialist Competitive Algorithm (IBA) algorithm has been tested on standard IEEE 57
bus test system and simulation results show clearly the high-quality performance of the projected algorithm in
reducing the real power loss.
Keywords: Optimal Reactive Power, Transmission loss, honey bee, foraging behaviour, waggle dance, bee’s
algorithm, swarm intelligence, swarm-based optimization, adaptive neighbourhood search, site abandonment,
random search
I. INTRODUCTION
Optimal reactive power dispatch (ORPD) problem is to minimize the real power loss and bus voltage deviation. Various
mathematical techniques like the gradient method [1-2], Newton method [3] and linear programming [4-7] have been
adopted to solve the optimal reactive power dispatch problem. Both the gradient and Newton methods have the
complexity in managing inequality constraints. If linear programming is applied then the input- output function has to be
uttered as a set of linear functions which mostly lead to loss of accuracy. The problem of voltage stability and collapse
play a major role in power system planning and operation [8]. Global optimization has received extensive research
awareness, and a great number of methods have been applied to solve this problem. Evolutionary algorithms such as
genetic algorithm have been already proposed to solve the reactive power flow problem [9, 10]. Evolutionary algorithm
is a heuristic approach used for minimization problems by utilizing nonlinear and non-differentiable continuous space
functions. In [11], Genetic algorithm has been used to solve optimal reactive power flow problem. In [12], Hybrid
differential evolution algorithm is proposed to improve the voltage stability index. In [13] Biogeography Based algorithm
is projected to solve the reactive power dispatch problem. In [14], a fuzzy based method is used to solve the optimal
reactive power scheduling method. In [15], an improved evolutionary programming is used to solve the optimal reactive
power dispatch problem. In [16], the optimal reactive power flow problem is solved by integrating a genetic algorithm
with a nonlinear interior point method. In [17], a pattern algorithm is used to solve ac-dc optimal reactive power flow
model with the generator capability limits. In [18], F. Capitanescu proposes a two-step approach to evaluate Reactive
power reserves with respect to operating constraints and voltage stability. In [19], a programming based approach is used
to solve the optimal reactive power dispatch problem. In [20], A. Kargarian et al present a probabilistic algorithm for
optimal reactive power provision in hybrid electricity markets with uncertain loads. This paper proposes a new Improved
Bees Algorithm (IBA) to solve the optimal reactive power dispatch problem. The aim of this paper is to solve optimal
reactive power problem by utilizing Bees Algorithm, introduced by Pham [21], inspired from the natural foraging
behaviour of honey bees. The IBA algorithm performs both an exploitative neighbourhood search combined with
arbitrary explorative search. The proposed algorithm IBA has been evaluated in standard IEEE 57 bus test system and
the simulation results show that our proposed approach outperforms all the entitled reported algorithms in minimization
of real power loss.
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 35
Paper Publications
II. PROBLEM FORMULATION
The optimal power flow problem is treated as a general minimization problem with constraints, and can be
mathematically written in the following form:
Minimize f(x, u) (1)
subject to g(x,u)=0 (2)
and
h(x, u) ≤ 0 (3)
where f(x,u) is the objective function. g(x.u) and h(x,u) are respectively the set of equality and inequality constraints. x is
the vector of state variables, and u is the vector of control variables.
The state variables are the load buses (PQ buses) voltages, angles, the generator reactive powers and the slack active
generator power:
x = Pg1, θ2, . . , θN, VL1, . , VLNL , Qg1, . . , Qgng
T
(4)
The control variables are the generator bus voltages, the shunt capacitors/reactors and the transformers tap-settings:
u = Vg, T, Qc
T
(5)
or
u = Vg1, … , Vgng , T1, . . , TNt , Qc1, . . , QcNc
T
(6)
Where ng, nt and nc are the number of generators, number of tap transformers and the number of shunt compensators
respectively.
III. OBJECTIVE FUNCTION
A. Active power loss
The objective of the reactive power dispatch is to minimize the active power loss in the transmission network, which can
be described as follows:
𝐹 = 𝑃𝐿 = 𝑔 𝑘𝑘∈𝑁𝑏𝑟 𝑉𝑖
2
+ 𝑉𝑗
2
− 2𝑉𝑖 𝑉𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 (7)
Or
𝐹 = 𝑃𝐿 = 𝑃𝑔𝑖 − 𝑃𝑑 = 𝑃𝑔𝑠𝑙𝑎𝑐𝑘 + 𝑃𝑔𝑖 − 𝑃𝑑
𝑁𝑔
𝑖≠𝑠𝑙𝑎𝑐𝑘𝑖∈𝑁𝑔 (8)
where gk : is the conductance of branch between nodes i and j, Nbr: is the total number of transmission lines in power
systems. Pd: is the total active power demand, Pgi: is the generator active power of unit i, and Pgsalck: is the generator
active power of slack bus.
B. Voltage profile improvement
For minimizing the voltage deviation in PQ buses, the objective function becomes:
𝐹 = 𝑃𝐿 + 𝜔 𝑣 × 𝑉𝐷 (9)
where ωv: is a weighting factor of voltage deviation.
VD is the voltage deviation given by:
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 36
Paper Publications
𝑉𝐷 = 𝑉𝑖 − 1
𝑁𝑝𝑞
𝑖=1 (10)
C. Equality Constraint
The equality constraint g(x,u) of the ORPD problem is represented by the power balance equation, where the total power
generation must cover the total power demand and the power losses:
𝑃𝐺 = 𝑃𝐷 + 𝑃𝐿 (11)
This equation is solved by running Newton Raphson load flow method, by calculating the active power of slack bus to
determine active power loss.
D. Inequality Constraints
The inequality constraints h(x,u) reflect the limits on components in the power system as well as the limits created to
ensure system security. Upper and lower bounds on the active power of slack bus, and reactive power of generators:
𝑃𝑔𝑠𝑙𝑎𝑐𝑘
𝑚𝑖𝑛
≤ 𝑃𝑔𝑠𝑙𝑎𝑐𝑘 ≤ 𝑃𝑔𝑠𝑙𝑎𝑐𝑘
𝑚𝑎𝑥
(12)
𝑄 𝑔𝑖
𝑚𝑖𝑛
≤ 𝑄 𝑔𝑖 ≤ 𝑄 𝑔𝑖
𝑚𝑎𝑥
, 𝑖 ∈ 𝑁𝑔 (13)
Upper and lower bounds on the bus voltage magnitudes:
𝑉𝑖
𝑚𝑖𝑛
≤ 𝑉𝑖 ≤ 𝑉𝑖
𝑚𝑎𝑥
, 𝑖 ∈ 𝑁 (14)
Upper and lower bounds on the transformers tap ratios:
𝑇𝑖
𝑚𝑖𝑛
≤ 𝑇𝑖 ≤ 𝑇𝑖
𝑚𝑎𝑥
, 𝑖 ∈ 𝑁 𝑇 (15)
Upper and lower bounds on the compensators reactive powers:
𝑄𝑐
𝑚𝑖𝑛
≤ 𝑄𝑐 ≤ 𝑄 𝐶
𝑚𝑎𝑥
, 𝑖 ∈ 𝑁𝐶 (16)
Where N is the total number of buses, NT is the total number of Transformers; Nc is the total number of shunt reactive
compensators.
IV. BEHAVIOUR OF HONEY BEES
A colony of honey bees can exploit a huge number of food sources in big fields and they can fly up to 12 km to exploit
food sources [22, 23]. The colony utilize about one-quarter of its members as searcher bees. The foraging process begins
with searching out hopeful flower patches by scout bees. The colony keeps a proportion of the scout bees during the
harvesting season. When the scout bees have found a flower area, they will look further in hope of finding an even
superior one [23]. The scout bees search for the better patches randomly [24]. The scout bees notify their peers waiting in
the hive about the eminence of the food source, based amongst other things, on sugar levels. The scout bees dump their
nectar and go to the dance floor in front of the hive to converse to the other bees by performing their dance, known as the
waggle dance [22]. The waggle dance is named based on the wagging run, which is used by the scout bees to
communicate information about the food source to the rest of the colony. The scout bees present the following
information by means of the waggle dance: the quality of the food source, the distance of the source from the hive and
the direction of the source [23- 25]. Figure 1a,b [25]. The scout then circles back, alternating a left and a right return path
. The speed/duration of the dance indicates the distance to the food source; the frequency of the waggles in the dance and
buzzing convey the quality of the source; see Figure 1c [25]. This information will influence the number of follower
bees.
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 37
Paper Publications
Fig1. (a) Orientation of waggle dance with respect to the sun; (b) Orientation of waggle dance with respect to the food source, hive
and sun; (c) The Waggle Dance and followers.
Fundamental parameters of the Bees Algorithm:
Quantity of scout bees in the selected patches - n
Quantity of best patches in the selected patches - m
Quantity of elite patches in the selected best patches- e
Quantity of recruited bees in the elite patches -nep
Quantity of recruited bees in the non-elite best patches- nsp
The size of neighbourhood for each patch - ngh
Quantity of iterations- Maxiter
Variation between value of the first and last iterations- diff
Bees Algorithm:
Create the initial population size as n, m, e, nep, set nsp, ngh, MaxIter, and set the error limit as Error.
i = 0
Generate preliminary population.
Calculate Fitness Value of initial population.
Arrange the initial population based on the fitness result.
While 𝑖 ≤ 𝑚𝑎𝑥𝐼𝑡𝑒𝑟 𝑜𝑟 𝑓𝑖𝑡𝑛𝑒𝑠𝑠 𝑣𝑎𝑙𝑢𝑒𝑖 − 𝑓𝑖𝑡𝑛𝑒𝑠𝑠𝑣𝑎𝑙𝑢𝑒𝑖−1 ≤ 𝐸𝑟𝑟𝑜𝑟
i. i = i + l;
ii. Choose the elite patches and non-elite best patches for neighbourhood search.
iii. Engage the forager bees to the elite patches and non-elite best patches.
iv. Calculate the fitness value of each patch.
v. Arrange the results based on their fitness.
vi. Distribute the rest of the bees for global search to the non-best locations.
vii. Calculate the fitness value of non-best patches.
viii. Arrange the overall results based on their fitness.
x. Run the algorithm until stop criteria met.
End
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 38
Paper Publications
V. IMPROVED BEES ALGORITHM BY ADAPTIVE NEIGHBOURHOOD SEARCH AND
SITE ABANDONMENT STRATEGY
This segment explains the proposed improvements to the bee‟s algorithm (BA) by applying adaptive transform to the
neighbourhood size and site abandonment approach simultaneously. Collective neighbourhood size change and site
abandonment (NSSA) approach has been attempted on the BA by Koc [26] who found that the convergence rate of a
NSSA-based BA can be sluggish, when the promising locations are far from the current best sites. Here an adaptive
neighbourhood size change and site abandonment (ANSSA) approach is proposed which will keep away from local
minima by changing the neighbourhood size adaptively. The ANSSA-based BA possesses both shrinking and
augmentation strategies according to the fitness evaluation. The primary move is to implement the shrinking approach.
This approach works on a best site after a definite number of repetitions. The approach works until the repetition stops.
If, in spite of the shrinking approach, the number of repetitions still increases for a definite number of iterations, then an
augmentation approach is utilized. Finally, if the number of repetitions still increases for a number of iterations after the
use of the augmentation approach, then that site is abandoned and a new site will be generated. Koc [26] utilized the
following parameter for shrinking the neighbourhood size and site abandonment approach: neighbourhood size = ngh,
the shrinking constant = sc, the abandoned sites = aband_site. In this study four more parameters are introduced. The first
is the number of repetitions for each site, denoted as keep_point. The keep_point records the number of repetitions for all
the repetitive results for best sites. The second parameter is called the “Repetition Number for the shrinking” is denoted
as rep_nshr; the number of shrinking is the number of repetitions necessary to start the shrinking strategy, as given in
Equations (17) and (18). The parameter is the “Repetition Number for the enhancement” is denoted as rep_nenh. This
parameter defines the number of repetitions until the end of the shrinking process, and the beginning of the enhancement
process as shown in Equations (17) and (19) [27,28]. The enhancement process works until the number of the repetitions
is equal to the rep_naban, which denotes the “Repetition Number for the abandonment process”. Hence a non-productive
site is abandoned and it is stored in aband_site list. If there is no better solution than the abandoned site at the end of the
searching process, this is the final solution.
𝑛𝑒𝑤 𝑛𝑔 𝑕 =
𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 ≤ 𝑟𝑒𝑝 𝑛𝑠 𝑕𝑟 𝑛𝑔𝑕
𝑟𝑒𝑝 𝑛𝑠 𝑕𝑟 < 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 ≤ 𝑟𝑒𝑝 𝑛𝑒𝑛 𝑕 𝑅1
𝑟𝑒𝑝 𝑛𝑒𝑛 𝑕 < 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 ≤ 𝑟𝑒𝑝 𝑛𝑎𝑏𝑎𝑛 𝑅2
𝑟𝑒𝑝 𝑛𝑎𝑏 𝑎𝑛 < 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 𝑛𝑔𝑕
(17)
𝑅1 = 𝑛𝑔𝑕 − 𝑛𝑔𝑕 ∗
𝑘𝑒𝑒 𝑝 𝑝𝑜𝑖𝑛𝑡 −𝑟𝑒𝑝 _𝑛𝑠 𝑕𝑟
100
∗ 𝑠𝑐 (18)
𝑅2 = 𝑛𝑔𝑕 + 𝑛𝑔𝑕 ∗
𝑘𝑒𝑒 𝑝 𝑝𝑜𝑖𝑛𝑡 −𝑟𝑒𝑝 _𝑛𝑒𝑛 𝑕
100
∗ 𝑠𝑐 (19)
VI. SIMULATION RESULTS
The proposed Improved Bees Algorithm (IBA) algorithm for solving ORPD problem is tested for standard IEEE-57 bus
power system. The IEEE 57-bus system data consists of 80 branches, seven generator-buses and 17 branches under load
tap setting transformer branches. The possible reactive power compensation buses are 18, 25 and 53. Bus 2, 3, 6, 8, 9 and
12 are PV buses and bus 1 is selected as slack-bus. In this case, the search space has 27 dimensions, i.e., the seven
generator voltages, 17 transformer taps, and three capacitor banks. The system variable limits are given in Table I. The
initial conditions for the IEEE-57 bus power system are given as follows:
Pload = 12.310 p.u. Qload = 3.322 p.u.
The total initial generations and power losses are obtained as follows:
𝑃𝐺 = 12.7634 p.u. 𝑄 𝐺 = 3.4468 p.u.
Ploss = 0.27351 p.u. Qloss = -1.2248 p.u.
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 39
Paper Publications
Table II shows the various system control variables i.e. generator bus voltages, shunt capacitances and transformer tap
settings obtained after IBA based optimization which are within their acceptable limits. In Table III, a comparison of
optimum results obtained from proposed IBA with other optimization techniques for ORPD mentioned in literature for
IEEE-57 bus power system is given. These results indicate the robustness of proposed IBA approach for providing better
optimal solution in case of IEEE-57 bus system.
TABLE I: VARIABLES LIMITS FOR IEEE-57 BUS POWER SYSTEM (P.U.)
REACTIVE POWER GENERATION LIMITS
BUS NO 1 2 3 6 8 9 12
QGMIN -1.2 -.014 -.02 -0.06 -1.2 -0.03 -0.3
QGMAX 2 0.4 0.5 0.24 2 0.08 1.54
VOLTAGE AND TAP SETTING LIMITS
VGMIN VGMAX VPQMIN VPQMAX TKMIN TKMAX
0.7 1.3 0.95 1.06 0.7 1.3
SHUNT CAPACITOR LIMITS
BUS NO 18 25 53
QCMIN 0 0 0
QCMAX 10 5.3 6.5
TABLE II: CONTROL VARIABLES OBTAINED AFTER OPTIMIZATION BY IBA METHOD FOR IEEE-57 BUS SYSTEM (P.U.)
Control
Variables
IBA
V1 1.2
V2 1.084
V3 1.073
V6 1.051
V8 1.074
V9 1.052
V12 1.061
Qc18 0.0843
Qc25 0.333
Qc53 0.0628
T4-18 1.016
T21-20 1.072
T24-25 0.973
T24-26 0.945
T7-29 1.092
T34-32 0.957
T11-41 1.015
T15-45 1.074
T14-46 0.943
T10-51 1.055
T13-49 1.075
T11-43 0.921
T40-56 0.911
T39-57 0.973
T9-55 0.985
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 40
Paper Publications
TABLE III: COMPARATIVE OPTIMIZATION RESULTS FOR IEEE-57 BUS POWER SYSTEM (P.U.)
S.No. Optimization
Algorithm
Best Solution Worst Solution Average
Solution
1 NLP [29] 0.25902 0.30854 0.27858
2 CGA [29] 0.25244 0.27507 0.26293
3 AGA [29] 0.24564 0.26671 0.25127
4 PSO-w [29] 0.24270 0.26152 0.24725
5 PSO-cf [29] 0.24280 0.26032 0.24698
6 CLPSO [29] 0.24515 0.24780 0.24673
7 SPSO-07 [29] 0.24430 0.25457 0.24752
8 L-DE [29] 0.27812 0.41909 0.33177
9 L-SACP-DE [29] 0.27915 0.36978 0.31032
10 L-SaDE [29] 0.24267 0.24391 0.24311
11 SOA [29] 0.24265 0.24280 0.24270
12 LM [30] 0.2484 0.2922 0.2641
13 MBEP1 [30] 0.2474 0.2848 0.2643
14 MBEP2 [30] 0.2482 0.283 0.2592
15 BES100 [30] 0.2438 0.263 0.2541
16 BES200 [30] 0.3417 0.2486 0.2443
17 Proposed IBA 0.22359 0.23492 0.23121
VII. CONCLUSION
IBA has been fruitfully applied for ORPD problem. The IBA based ORPD is tested in standard IEEE-57 bus system.
Performance comparisons with well-known population-based algorithms give cheering results. IBA emerges to find good
solutions when compared to that of other algorithms. The simulation results presented in previous section prove the
ability of IBA approach to arrive at near global optimal solution.
REFERENCES
[1] O.Alsac,and B. Scott, “Optimal load flow with steady state security”,IEEE Transaction. PAS -1973, pp. 745-751.
[2] Lee K Y ,Paru Y M , Oritz J L –A united approach to optimal real and reactive power dispatch , IEEE Transactions
on power Apparatus and systems 1985: PAS-104 : 1147-1153
[3] A.Monticelli , M .V.F Pereira ,and S. Granville , “Security constrained optimal power flow with post contingency
corrective rescheduling” , IEEE Transactions on Power Systems :PWRS-2, No. 1, pp.175-182.,1987.
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 41
Paper Publications
[4] Deeb N ,Shahidehpur S.M ,Linear reactive power optimization in a large power network using the decomposition
approach. IEEE Transactions on power system 1990: 5(2) : 428-435
[5] E. Hobson ,‟Network consrained reactive power control using linear programming, „ IEEE Transactions on power
systems PAS -99 (4) ,pp 868=877, 1980
[6] K.Y Lee ,Y.M Park , and J.L Oritz, “Fuel –cost optimization for both real and reactive power dispatches” , IEE Proc;
131C,(3), pp.85-93.
[7] M.K. Mangoli, and K.Y. Lee, “Optimal real and reactive power control using linear programming” , Electr.Power
Syst.Res, Vol.26, pp.1-10,1993.
[8] C.A. Canizares , A.C.Z.de Souza and V.H. Quintana , “ Comparison of performance indices for detection of
proximity to voltage collapse ,‟‟ vol. 11. no.3 , pp.1441-1450, Aug 1996 .
[9] S.R.Paranjothi ,and K.Anburaja, “Optimal power flow using refined genetic algorithm”, Electr.Power Compon.Syst ,
Vol. 30, 1055-1063,2002.
[10] D. Devaraj, and B. Yeganarayana, “Genetic algorithm based optimal power flow for security enhancement”, IEE
proc-Generation.Transmission and. Distribution; 152, 6 November 2005.
[11] A. Berizzi, C. Bovo, M. Merlo, and M. Delfanti, “A ga approach to compare orpf objective functions including
secondary voltage regulation,” Electric Power Systems Research, vol. 84, no. 1, pp. 187 – 194, 2012.
[12] C.-F. Yang, G. G. Lai, C.-H. Lee, C.-T. Su, and G. W. Chang, “Optimal setting of reactive compensation devices
with an improved voltage stability index for voltage stability enhancement,” International Journal of Electrical Power and
Energy Systems, vol. 37, no. 1, pp. 50 – 57, 2012.
[13] P. Roy, S. Ghoshal, and S. Thakur, “Optimal var control for improvements in voltage profiles and for real power
loss minimization using biogeography based optimization,” International Journal of Electrical Power and Energy
Systems, vol. 43, no. 1, pp. 830 – 838, 2012.
[14] B. Venkatesh, G. Sadasivam, and M. Khan, “A new optimal reactive power scheduling method for loss
minimization and voltage stability margin maximization using successive multi-objective fuzzy lp technique,” IEEE
Transactions on Power Systems, vol. 15, no. 2, pp. 844 – 851, may 2000.
[15] W. Yan, S. Lu, and D. Yu, “A novel optimal reactive power dispatch method based on an improved hybrid
evolutionary programming technique,” IEEE Transactions on Power Systems, vol. 19, no. 2, pp. 913 – 918, may 2004.
[16] W. Yan, F. Liu, C. Chung, and K. Wong, “A hybrid genetic algorithminterior point method for optimal reactive
power flow,” IEEE Transactions on Power Systems, vol. 21, no. 3, pp. 1163 –1169, aug. 2006.
[17] J. Yu, W. Yan, W. Li, C. Chung, and K. Wong, “An unfixed piecewiseoptimal reactive power-flow model and its
algorithm for ac-dc systems,” IEEE Transactions on Power Systems, vol. 23, no. 1, pp. 170 –176, feb. 2008.
[18] F. Capitanescu, “Assessing reactive power reserves with respect to operating constraints and voltage stability,” IEEE
Transactions on Power Systems, vol. 26, no. 4, pp. 2224–2234, nov. 2011.
[19] Z. Hu, X. Wang, and G. Taylor, “Stochastic optimal reactive power dispatch: Formulation and solution method,”
International Journal of Electrical Power and Energy Systems, vol. 32, no. 6, pp. 615 – 621, 2010.
[20] A. Kargarian, M. Raoofat, and M. Mohammadi, “Probabilistic reactive power procurement in hybrid electricity
markets with uncertain loads,” Electric Power Systems Research, vol. 82, no. 1, pp. 68 – 80, 2012.
[21] Pham, D.T.; Ghanbarzadeh, A.; Koc, E.; Otri, S.; Rahim, S.; Zaidi, M. The Bees Algorithm, Technical Note;
Manufacturing Engineering Center, Cardiff University: Cardiff, UK, 2005.
[22] Seeley, T.D. The Wisdom of the Hive: The Social Physiology of Honey Bee Colonies; Harvard University Press:
Cambridge, MA, USA, 2009.
[23] Gould, J.L.; Gould, C.G. The Honey Bee; Scientific American Library: New York, NY, USA, 1988.
[24] Von Frisch, K. Bees: Their Vision, Chemical Senses, and Language; Cornell University Press: Ithaca, NY, USA,
1950.
[25] Huang, Z. Behavioral Communications: The Waggle Dance. Available online: http://photo.bees.net/
biology/ch6/dance2.html (accessed on 29 June 2013).
[26] Koc, E. The Bees Algorithm Theory, Improvements and Applications. Ph.D Thesis, Cardiff University, Cardiff, UK,
2010.
International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE)
Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org
Page | 42
Paper Publications
[27] Yuce, B. Novel Computational Technique for Determining Depth Using the Bees Algorithm and Blind Image
Deconvolution. Ph.D Thesis, Cardiff University, Cardiff, UK, 2012.
[28] Baris Yuce, Michael S. Packianather , Ernesto Mastrocinque , Duc Truong Pham and Alfredo Lambiase “Honey
Bees Inspired Optimization Method” : The Bees Algorithm Insects 2013, 4, 646-662; doi:10.3390/insects4040646
[29] Chaohua Dai, Weirong Chen, Yunfang Zhu, and Xuexia Zhang, “Seeker optimization algorithm for optimal
reactive power dispatch,” IEEE Trans. Power Systems, Vol. 24, No. 3, August 2009, pp. 1218-1231.
[30] J. R. Gomes and 0. R. Saavedra, “Optimal reactive power dispatch using evolutionary computation: Extended
algorithms,” IEE Proc.-Gener. Transm. Distrib.. Vol. 146, No. 6. Nov. 1999.
Author Biography:
K. Lenin has received his B.E., Degree, electrical and electronics engineering in 1999 from
university of madras, Chennai, India and M.E., Degree in power systems in 2000 from Annamalai
University, Tamil Nadu, India. At present pursuing Ph.D., degree at JNTU, Hyderabad, India.
Bhumanapally. RavindhranathReddy, Born on 3rd September,1969. Got his B.Tech in Electrical &
Electronics Engineering from the J.N.T.U. College of Engg., Anantapur in the year 1991. Completed
his M.Tech in Energy Systems in IPGSR of J.N.T.University Hyderabad in the year 1997. Obtained
his doctoral degree from JNTUA,Anantapur University in the field of Electrical Power Systems.
Published 12 Research Papers and presently guiding 6 Ph.D. Scholars. He was specialized in Power
Systems, High Voltage Engineering and Control Systems. His research interests include Simulation
studies on Transients of different power system equipment.
M. Surya Kalavathi has received her B.Tech. Electrical and Electronics Engineering from SVU,
Andhra Pradesh, India and M.Tech, power system operation and control from SVU, Andhra Pradesh,
India. she received her Phd. Degree from JNTU, hyderabad and Post doc. From CMU – USA.
Currently she is Professor and Head of the electrical and electronics engineering department in
JNTU, Hyderabad, India and she has Published 16 Research Papers and presently guiding 5 Ph.D.
Scholars. She has specialised in Power Systems, High Voltage Engineering and Control Systems.
Her research interests include Simulation studies on Transients of different power system equipment. She has 18 years
of experience. She has invited for various lectures in institutes.

More Related Content

PDF
Dwindling of real power loss by using Improved Bees Algorithm
PDF
5. 9375 11036-1-sm-1 20 apr 18 mar 16oct2017 ed iqbal qc
PDF
Enriched Firefly Algorithm for Solving Reactive Power Problem
PDF
Voltage Profile Enhancement and Reduction of Real Power loss by Hybrid Biogeo...
PDF
Al36228233
PDF
Optimal Power Flow with Reactive Power Compensation for Cost And Loss Minimiz...
PDF
Cuckoo Search Algorithm for Congestion Alleviation with Incorporation of Wind...
PDF
Economic/Emission Load Dispatch Using Artificial Bee Colony Algorithm
Dwindling of real power loss by using Improved Bees Algorithm
5. 9375 11036-1-sm-1 20 apr 18 mar 16oct2017 ed iqbal qc
Enriched Firefly Algorithm for Solving Reactive Power Problem
Voltage Profile Enhancement and Reduction of Real Power loss by Hybrid Biogeo...
Al36228233
Optimal Power Flow with Reactive Power Compensation for Cost And Loss Minimiz...
Cuckoo Search Algorithm for Congestion Alleviation with Incorporation of Wind...
Economic/Emission Load Dispatch Using Artificial Bee Colony Algorithm

What's hot (20)

PDF
Reduction of Active Power Loss byUsing Adaptive Cat Swarm Optimization
PDF
Economic and Emission Dispatch using Whale Optimization Algorithm (WOA)
PDF
Economic Load Dispatch Problem with Valve – Point Effect Using a Binary Bat A...
PDF
NOVEL PSO STRATEGY FOR TRANSMISSION CONGESTION MANAGEMENT
PDF
Optimal Power Generation in Energy-Deficient Scenarios Using Bagging Ensembles
PDF
Gy3312241229
PDF
A Solution to Optimal Power Flow Problem using Artificial Bee Colony Algorith...
PDF
International Journal of Engineering Research and Development
PDF
Optimal Allocation of Capacitor Bank in Radial Distribution System using Anal...
PDF
Capacitor Placement Using Bat Algorithm for Maximum Annual Savings in Radial ...
PDF
01 16286 32182-1-sm multiple (edit)
PDF
Optimal Power System Planning with Renewable DGs with Reactive Power Consider...
PDF
Atmosphere Clouds Model Algorithm for Solving Optimal Reactive Power Dispatch...
PDF
Capacitor Placement and Reconfiguration of Distribution System with hybrid Fu...
PDF
Power losses reduction of power transmission network using optimal location o...
PDF
Hybrid Quantum Genetic Particle Swarm Optimization Algorithm For Solving Opti...
PDF
OPTIMIZATION OF COMBINED ECONOMIC EMISSION DISPATCH PROBLEM USING ARTIFICIAL ...
PDF
Application of-computational-intelligence-techniques-for-economic-load-dispatch
PDF
Multi Objective Directed Bee Colony Optimization for Economic Load Dispatch W...
PDF
HYBRID PARTICLE SWARM OPTIMIZATION FOR SOLVING MULTI-AREA ECONOMIC DISPATCH P...
Reduction of Active Power Loss byUsing Adaptive Cat Swarm Optimization
Economic and Emission Dispatch using Whale Optimization Algorithm (WOA)
Economic Load Dispatch Problem with Valve – Point Effect Using a Binary Bat A...
NOVEL PSO STRATEGY FOR TRANSMISSION CONGESTION MANAGEMENT
Optimal Power Generation in Energy-Deficient Scenarios Using Bagging Ensembles
Gy3312241229
A Solution to Optimal Power Flow Problem using Artificial Bee Colony Algorith...
International Journal of Engineering Research and Development
Optimal Allocation of Capacitor Bank in Radial Distribution System using Anal...
Capacitor Placement Using Bat Algorithm for Maximum Annual Savings in Radial ...
01 16286 32182-1-sm multiple (edit)
Optimal Power System Planning with Renewable DGs with Reactive Power Consider...
Atmosphere Clouds Model Algorithm for Solving Optimal Reactive Power Dispatch...
Capacitor Placement and Reconfiguration of Distribution System with hybrid Fu...
Power losses reduction of power transmission network using optimal location o...
Hybrid Quantum Genetic Particle Swarm Optimization Algorithm For Solving Opti...
OPTIMIZATION OF COMBINED ECONOMIC EMISSION DISPATCH PROBLEM USING ARTIFICIAL ...
Application of-computational-intelligence-techniques-for-economic-load-dispatch
Multi Objective Directed Bee Colony Optimization for Economic Load Dispatch W...
HYBRID PARTICLE SWARM OPTIMIZATION FOR SOLVING MULTI-AREA ECONOMIC DISPATCH P...
Ad

Viewers also liked (20)

PDF
Minimization of Overall Losses of a Distribution System under Contingency Con...
PDF
Effect of Mobile Banking on Customer Satisfaction in Selected Banks in Trans-...
PDF
Organizational Efficiency and Internal Customer Satisfaction: A Case Study of...
PDF
THE IMPACT OF SERVICE QUALITY AND BRAND AWARENESS ON BRAND LOYALTY: (A STUDY ...
PDF
Design and Implementation of Area Efficiency AES Algoritham with FPGA and ASIC
PDF
The Impact of Key Performance Indicators (KPIs) on Talent Development
PDF
Role of Procurement on the Completion of Projects in the Energy Sector in Kenya
PDF
EVOLVING TRENDS FOR ENHANCING THE ACCURACY OF FAULT LOCATION IN POWER DISTRIB...
PDF
Comparing the effect of PI and I controllers on a four tank process
PDF
PS2O, Hybrid Evolutionary-Conventional Algorithm, Genetical Swarm Optimizatio...
PDF
A Theoretical Statistical Measurement Model Analysis on Human Capital Economi...
PDF
DETERMINANTS OF SERVICE DELIVERY IN PUBLIC PRIMARY SCHOOLS IN WEST POKOT COUNTY
PDF
Factors Influencing Implementation of Strategic Plans in Non-Governmental Org...
PDF
Prevalence of SHRM Practices in Operational NGOs in Kerala
PDF
TALENT ACQUISITION PROCEDURE OF HR PROCESS: A CASE OF HCL INFOSYSTEM PVT LTD.
PDF
Combination of Immune Genetic Particle Swarm Optimization algorithm with BP a...
PDF
Hybrid Quantum Genetic Particle Swarm Optimization Algorithm For Solving Opti...
PDF
EFFECTS OF SUPPLIER EFFICIENCY ON PROCUREMENT PERFORMANCE IN KENYA SEED COMPA...
PDF
Performance Analysis of Digital Watermarking Of Video in the Spatial Domain
PDF
MAXIMUM POWER POINT TRACKING TECHNIQUES: A REVIEW
Minimization of Overall Losses of a Distribution System under Contingency Con...
Effect of Mobile Banking on Customer Satisfaction in Selected Banks in Trans-...
Organizational Efficiency and Internal Customer Satisfaction: A Case Study of...
THE IMPACT OF SERVICE QUALITY AND BRAND AWARENESS ON BRAND LOYALTY: (A STUDY ...
Design and Implementation of Area Efficiency AES Algoritham with FPGA and ASIC
The Impact of Key Performance Indicators (KPIs) on Talent Development
Role of Procurement on the Completion of Projects in the Energy Sector in Kenya
EVOLVING TRENDS FOR ENHANCING THE ACCURACY OF FAULT LOCATION IN POWER DISTRIB...
Comparing the effect of PI and I controllers on a four tank process
PS2O, Hybrid Evolutionary-Conventional Algorithm, Genetical Swarm Optimizatio...
A Theoretical Statistical Measurement Model Analysis on Human Capital Economi...
DETERMINANTS OF SERVICE DELIVERY IN PUBLIC PRIMARY SCHOOLS IN WEST POKOT COUNTY
Factors Influencing Implementation of Strategic Plans in Non-Governmental Org...
Prevalence of SHRM Practices in Operational NGOs in Kerala
TALENT ACQUISITION PROCEDURE OF HR PROCESS: A CASE OF HCL INFOSYSTEM PVT LTD.
Combination of Immune Genetic Particle Swarm Optimization algorithm with BP a...
Hybrid Quantum Genetic Particle Swarm Optimization Algorithm For Solving Opti...
EFFECTS OF SUPPLIER EFFICIENCY ON PROCUREMENT PERFORMANCE IN KENYA SEED COMPA...
Performance Analysis of Digital Watermarking Of Video in the Spatial Domain
MAXIMUM POWER POINT TRACKING TECHNIQUES: A REVIEW
Ad

Similar to Dwindling of real power loss by using Improved Bees Algorithm (20)

PDF
Optimal economic environmental power dispatch by using artificial bee colony ...
PDF
FLOWER POLLINATION ALGORITHM FOR SOLVING OPTIMAL REACTIVE POWER DISPATCH PRO...
PDF
FLOWER POLLINATION ALGORITHM FOR SOLVING OPTIMAL REACTIVE POWER DISPATCH PROBLEM
PDF
Passerine swarm optimization algorithm for solving optimal reactive power dis...
PDF
Firefly Algorithm to Opmimal Distribution of Reactive Power Compensation Units
PDF
Artificial bee colony algorithm applied to optimal power flow solution incorp...
PDF
Modified Monkey Optimization Algorithm for Solving Optimal Reactive Power Dis...
PDF
Wolf Search Algorithm for Solving Optimal Reactive Power Dispatch Problem
PDF
Optimal Allocation and Sizing of Distributed Generation using Artificial Bee ...
PDF
Improved Spider Algorithm for Solving Optimal Reactive Power Dispatch Problem
PDF
Optimal Reactive Power Dispatch using Crow Search Algorithm
PDF
Improved Spider Algorithm for Solving Optimal Reactive Power Dispatch Problem
PDF
Design of pss using bees colony intelligence
PDF
IRJET- Modified BEE Swarming Algoritm to Emission Constrained Economic Dispat...
PDF
Combination of Immune Genetic Particle Swarm Optimization algorithm with BP a...
PDF
Optimal Reactive Power Scheduling Using Cuckoo Search Algorithm
PDF
Biogeography Based Optimization Approach for Optimal Power Flow Problem Consi...
PDF
K346670
PDF
Economic dispatch by optimization techniques
PDF
A Particle Swarm Optimization for Reactive Power Optimization
Optimal economic environmental power dispatch by using artificial bee colony ...
FLOWER POLLINATION ALGORITHM FOR SOLVING OPTIMAL REACTIVE POWER DISPATCH PRO...
FLOWER POLLINATION ALGORITHM FOR SOLVING OPTIMAL REACTIVE POWER DISPATCH PROBLEM
Passerine swarm optimization algorithm for solving optimal reactive power dis...
Firefly Algorithm to Opmimal Distribution of Reactive Power Compensation Units
Artificial bee colony algorithm applied to optimal power flow solution incorp...
Modified Monkey Optimization Algorithm for Solving Optimal Reactive Power Dis...
Wolf Search Algorithm for Solving Optimal Reactive Power Dispatch Problem
Optimal Allocation and Sizing of Distributed Generation using Artificial Bee ...
Improved Spider Algorithm for Solving Optimal Reactive Power Dispatch Problem
Optimal Reactive Power Dispatch using Crow Search Algorithm
Improved Spider Algorithm for Solving Optimal Reactive Power Dispatch Problem
Design of pss using bees colony intelligence
IRJET- Modified BEE Swarming Algoritm to Emission Constrained Economic Dispat...
Combination of Immune Genetic Particle Swarm Optimization algorithm with BP a...
Optimal Reactive Power Scheduling Using Cuckoo Search Algorithm
Biogeography Based Optimization Approach for Optimal Power Flow Problem Consi...
K346670
Economic dispatch by optimization techniques
A Particle Swarm Optimization for Reactive Power Optimization

Recently uploaded (20)

PPT
Mechanical Engineering MATERIALS Selection
PDF
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
PPTX
bas. eng. economics group 4 presentation 1.pptx
PPTX
UNIT-1 - COAL BASED THERMAL POWER PLANTS
PPTX
Geodesy 1.pptx...............................................
PDF
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
PDF
composite construction of structures.pdf
PPTX
IOT PPTs Week 10 Lecture Material.pptx of NPTEL Smart Cities contd
PDF
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
PPTX
CYBER-CRIMES AND SECURITY A guide to understanding
PPTX
Welding lecture in detail for understanding
PPT
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PDF
Embodied AI: Ushering in the Next Era of Intelligent Systems
PPTX
Construction Project Organization Group 2.pptx
PDF
Digital Logic Computer Design lecture notes
PPTX
Sustainable Sites - Green Building Construction
PDF
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
PPTX
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
PPTX
CARTOGRAPHY AND GEOINFORMATION VISUALIZATION chapter1 NPTE (2).pptx
Mechanical Engineering MATERIALS Selection
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
bas. eng. economics group 4 presentation 1.pptx
UNIT-1 - COAL BASED THERMAL POWER PLANTS
Geodesy 1.pptx...............................................
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
composite construction of structures.pdf
IOT PPTs Week 10 Lecture Material.pptx of NPTEL Smart Cities contd
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
CYBER-CRIMES AND SECURITY A guide to understanding
Welding lecture in detail for understanding
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
Embodied AI: Ushering in the Next Era of Intelligent Systems
Construction Project Organization Group 2.pptx
Digital Logic Computer Design lecture notes
Sustainable Sites - Green Building Construction
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
CARTOGRAPHY AND GEOINFORMATION VISUALIZATION chapter1 NPTE (2).pptx

Dwindling of real power loss by using Improved Bees Algorithm

  • 1. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 34 Paper Publications Dwindling of real power loss by using Improved Bees Algorithm K. Lenin, B. Ravindranath Reddy, and M. Surya Kalavathi Jawaharlal Nehru Technological University Kukatpally, Hyderabad 500 085, India Abstract: In this paper, a new Improved Bees Algorithm (IBA) is proposed for solving reactive power dispatch problem. The aim of this paper is to utilize an optimization algorithm called the improved Bees Algorithm, inspired from the natural foraging behaviour of honey bees, to solve the reactive power dispatch problem. The IBA algorithm executes both an exploitative neighbourhood search combined with arbitrary explorative search. The proposed Improved Imperialist Competitive Algorithm (IBA) algorithm has been tested on standard IEEE 57 bus test system and simulation results show clearly the high-quality performance of the projected algorithm in reducing the real power loss. Keywords: Optimal Reactive Power, Transmission loss, honey bee, foraging behaviour, waggle dance, bee’s algorithm, swarm intelligence, swarm-based optimization, adaptive neighbourhood search, site abandonment, random search I. INTRODUCTION Optimal reactive power dispatch (ORPD) problem is to minimize the real power loss and bus voltage deviation. Various mathematical techniques like the gradient method [1-2], Newton method [3] and linear programming [4-7] have been adopted to solve the optimal reactive power dispatch problem. Both the gradient and Newton methods have the complexity in managing inequality constraints. If linear programming is applied then the input- output function has to be uttered as a set of linear functions which mostly lead to loss of accuracy. The problem of voltage stability and collapse play a major role in power system planning and operation [8]. Global optimization has received extensive research awareness, and a great number of methods have been applied to solve this problem. Evolutionary algorithms such as genetic algorithm have been already proposed to solve the reactive power flow problem [9, 10]. Evolutionary algorithm is a heuristic approach used for minimization problems by utilizing nonlinear and non-differentiable continuous space functions. In [11], Genetic algorithm has been used to solve optimal reactive power flow problem. In [12], Hybrid differential evolution algorithm is proposed to improve the voltage stability index. In [13] Biogeography Based algorithm is projected to solve the reactive power dispatch problem. In [14], a fuzzy based method is used to solve the optimal reactive power scheduling method. In [15], an improved evolutionary programming is used to solve the optimal reactive power dispatch problem. In [16], the optimal reactive power flow problem is solved by integrating a genetic algorithm with a nonlinear interior point method. In [17], a pattern algorithm is used to solve ac-dc optimal reactive power flow model with the generator capability limits. In [18], F. Capitanescu proposes a two-step approach to evaluate Reactive power reserves with respect to operating constraints and voltage stability. In [19], a programming based approach is used to solve the optimal reactive power dispatch problem. In [20], A. Kargarian et al present a probabilistic algorithm for optimal reactive power provision in hybrid electricity markets with uncertain loads. This paper proposes a new Improved Bees Algorithm (IBA) to solve the optimal reactive power dispatch problem. The aim of this paper is to solve optimal reactive power problem by utilizing Bees Algorithm, introduced by Pham [21], inspired from the natural foraging behaviour of honey bees. The IBA algorithm performs both an exploitative neighbourhood search combined with arbitrary explorative search. The proposed algorithm IBA has been evaluated in standard IEEE 57 bus test system and the simulation results show that our proposed approach outperforms all the entitled reported algorithms in minimization of real power loss.
  • 2. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 35 Paper Publications II. PROBLEM FORMULATION The optimal power flow problem is treated as a general minimization problem with constraints, and can be mathematically written in the following form: Minimize f(x, u) (1) subject to g(x,u)=0 (2) and h(x, u) ≤ 0 (3) where f(x,u) is the objective function. g(x.u) and h(x,u) are respectively the set of equality and inequality constraints. x is the vector of state variables, and u is the vector of control variables. The state variables are the load buses (PQ buses) voltages, angles, the generator reactive powers and the slack active generator power: x = Pg1, θ2, . . , θN, VL1, . , VLNL , Qg1, . . , Qgng T (4) The control variables are the generator bus voltages, the shunt capacitors/reactors and the transformers tap-settings: u = Vg, T, Qc T (5) or u = Vg1, … , Vgng , T1, . . , TNt , Qc1, . . , QcNc T (6) Where ng, nt and nc are the number of generators, number of tap transformers and the number of shunt compensators respectively. III. OBJECTIVE FUNCTION A. Active power loss The objective of the reactive power dispatch is to minimize the active power loss in the transmission network, which can be described as follows: 𝐹 = 𝑃𝐿 = 𝑔 𝑘𝑘∈𝑁𝑏𝑟 𝑉𝑖 2 + 𝑉𝑗 2 − 2𝑉𝑖 𝑉𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 (7) Or 𝐹 = 𝑃𝐿 = 𝑃𝑔𝑖 − 𝑃𝑑 = 𝑃𝑔𝑠𝑙𝑎𝑐𝑘 + 𝑃𝑔𝑖 − 𝑃𝑑 𝑁𝑔 𝑖≠𝑠𝑙𝑎𝑐𝑘𝑖∈𝑁𝑔 (8) where gk : is the conductance of branch between nodes i and j, Nbr: is the total number of transmission lines in power systems. Pd: is the total active power demand, Pgi: is the generator active power of unit i, and Pgsalck: is the generator active power of slack bus. B. Voltage profile improvement For minimizing the voltage deviation in PQ buses, the objective function becomes: 𝐹 = 𝑃𝐿 + 𝜔 𝑣 × 𝑉𝐷 (9) where ωv: is a weighting factor of voltage deviation. VD is the voltage deviation given by:
  • 3. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 36 Paper Publications 𝑉𝐷 = 𝑉𝑖 − 1 𝑁𝑝𝑞 𝑖=1 (10) C. Equality Constraint The equality constraint g(x,u) of the ORPD problem is represented by the power balance equation, where the total power generation must cover the total power demand and the power losses: 𝑃𝐺 = 𝑃𝐷 + 𝑃𝐿 (11) This equation is solved by running Newton Raphson load flow method, by calculating the active power of slack bus to determine active power loss. D. Inequality Constraints The inequality constraints h(x,u) reflect the limits on components in the power system as well as the limits created to ensure system security. Upper and lower bounds on the active power of slack bus, and reactive power of generators: 𝑃𝑔𝑠𝑙𝑎𝑐𝑘 𝑚𝑖𝑛 ≤ 𝑃𝑔𝑠𝑙𝑎𝑐𝑘 ≤ 𝑃𝑔𝑠𝑙𝑎𝑐𝑘 𝑚𝑎𝑥 (12) 𝑄 𝑔𝑖 𝑚𝑖𝑛 ≤ 𝑄 𝑔𝑖 ≤ 𝑄 𝑔𝑖 𝑚𝑎𝑥 , 𝑖 ∈ 𝑁𝑔 (13) Upper and lower bounds on the bus voltage magnitudes: 𝑉𝑖 𝑚𝑖𝑛 ≤ 𝑉𝑖 ≤ 𝑉𝑖 𝑚𝑎𝑥 , 𝑖 ∈ 𝑁 (14) Upper and lower bounds on the transformers tap ratios: 𝑇𝑖 𝑚𝑖𝑛 ≤ 𝑇𝑖 ≤ 𝑇𝑖 𝑚𝑎𝑥 , 𝑖 ∈ 𝑁 𝑇 (15) Upper and lower bounds on the compensators reactive powers: 𝑄𝑐 𝑚𝑖𝑛 ≤ 𝑄𝑐 ≤ 𝑄 𝐶 𝑚𝑎𝑥 , 𝑖 ∈ 𝑁𝐶 (16) Where N is the total number of buses, NT is the total number of Transformers; Nc is the total number of shunt reactive compensators. IV. BEHAVIOUR OF HONEY BEES A colony of honey bees can exploit a huge number of food sources in big fields and they can fly up to 12 km to exploit food sources [22, 23]. The colony utilize about one-quarter of its members as searcher bees. The foraging process begins with searching out hopeful flower patches by scout bees. The colony keeps a proportion of the scout bees during the harvesting season. When the scout bees have found a flower area, they will look further in hope of finding an even superior one [23]. The scout bees search for the better patches randomly [24]. The scout bees notify their peers waiting in the hive about the eminence of the food source, based amongst other things, on sugar levels. The scout bees dump their nectar and go to the dance floor in front of the hive to converse to the other bees by performing their dance, known as the waggle dance [22]. The waggle dance is named based on the wagging run, which is used by the scout bees to communicate information about the food source to the rest of the colony. The scout bees present the following information by means of the waggle dance: the quality of the food source, the distance of the source from the hive and the direction of the source [23- 25]. Figure 1a,b [25]. The scout then circles back, alternating a left and a right return path . The speed/duration of the dance indicates the distance to the food source; the frequency of the waggles in the dance and buzzing convey the quality of the source; see Figure 1c [25]. This information will influence the number of follower bees.
  • 4. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 37 Paper Publications Fig1. (a) Orientation of waggle dance with respect to the sun; (b) Orientation of waggle dance with respect to the food source, hive and sun; (c) The Waggle Dance and followers. Fundamental parameters of the Bees Algorithm: Quantity of scout bees in the selected patches - n Quantity of best patches in the selected patches - m Quantity of elite patches in the selected best patches- e Quantity of recruited bees in the elite patches -nep Quantity of recruited bees in the non-elite best patches- nsp The size of neighbourhood for each patch - ngh Quantity of iterations- Maxiter Variation between value of the first and last iterations- diff Bees Algorithm: Create the initial population size as n, m, e, nep, set nsp, ngh, MaxIter, and set the error limit as Error. i = 0 Generate preliminary population. Calculate Fitness Value of initial population. Arrange the initial population based on the fitness result. While 𝑖 ≤ 𝑚𝑎𝑥𝐼𝑡𝑒𝑟 𝑜𝑟 𝑓𝑖𝑡𝑛𝑒𝑠𝑠 𝑣𝑎𝑙𝑢𝑒𝑖 − 𝑓𝑖𝑡𝑛𝑒𝑠𝑠𝑣𝑎𝑙𝑢𝑒𝑖−1 ≤ 𝐸𝑟𝑟𝑜𝑟 i. i = i + l; ii. Choose the elite patches and non-elite best patches for neighbourhood search. iii. Engage the forager bees to the elite patches and non-elite best patches. iv. Calculate the fitness value of each patch. v. Arrange the results based on their fitness. vi. Distribute the rest of the bees for global search to the non-best locations. vii. Calculate the fitness value of non-best patches. viii. Arrange the overall results based on their fitness. x. Run the algorithm until stop criteria met. End
  • 5. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 38 Paper Publications V. IMPROVED BEES ALGORITHM BY ADAPTIVE NEIGHBOURHOOD SEARCH AND SITE ABANDONMENT STRATEGY This segment explains the proposed improvements to the bee‟s algorithm (BA) by applying adaptive transform to the neighbourhood size and site abandonment approach simultaneously. Collective neighbourhood size change and site abandonment (NSSA) approach has been attempted on the BA by Koc [26] who found that the convergence rate of a NSSA-based BA can be sluggish, when the promising locations are far from the current best sites. Here an adaptive neighbourhood size change and site abandonment (ANSSA) approach is proposed which will keep away from local minima by changing the neighbourhood size adaptively. The ANSSA-based BA possesses both shrinking and augmentation strategies according to the fitness evaluation. The primary move is to implement the shrinking approach. This approach works on a best site after a definite number of repetitions. The approach works until the repetition stops. If, in spite of the shrinking approach, the number of repetitions still increases for a definite number of iterations, then an augmentation approach is utilized. Finally, if the number of repetitions still increases for a number of iterations after the use of the augmentation approach, then that site is abandoned and a new site will be generated. Koc [26] utilized the following parameter for shrinking the neighbourhood size and site abandonment approach: neighbourhood size = ngh, the shrinking constant = sc, the abandoned sites = aband_site. In this study four more parameters are introduced. The first is the number of repetitions for each site, denoted as keep_point. The keep_point records the number of repetitions for all the repetitive results for best sites. The second parameter is called the “Repetition Number for the shrinking” is denoted as rep_nshr; the number of shrinking is the number of repetitions necessary to start the shrinking strategy, as given in Equations (17) and (18). The parameter is the “Repetition Number for the enhancement” is denoted as rep_nenh. This parameter defines the number of repetitions until the end of the shrinking process, and the beginning of the enhancement process as shown in Equations (17) and (19) [27,28]. The enhancement process works until the number of the repetitions is equal to the rep_naban, which denotes the “Repetition Number for the abandonment process”. Hence a non-productive site is abandoned and it is stored in aband_site list. If there is no better solution than the abandoned site at the end of the searching process, this is the final solution. 𝑛𝑒𝑤 𝑛𝑔 𝑕 = 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 ≤ 𝑟𝑒𝑝 𝑛𝑠 𝑕𝑟 𝑛𝑔𝑕 𝑟𝑒𝑝 𝑛𝑠 𝑕𝑟 < 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 ≤ 𝑟𝑒𝑝 𝑛𝑒𝑛 𝑕 𝑅1 𝑟𝑒𝑝 𝑛𝑒𝑛 𝑕 < 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 ≤ 𝑟𝑒𝑝 𝑛𝑎𝑏𝑎𝑛 𝑅2 𝑟𝑒𝑝 𝑛𝑎𝑏 𝑎𝑛 < 𝑘𝑒𝑒𝑝 𝑝𝑜𝑖𝑛𝑡 𝑛𝑔𝑕 (17) 𝑅1 = 𝑛𝑔𝑕 − 𝑛𝑔𝑕 ∗ 𝑘𝑒𝑒 𝑝 𝑝𝑜𝑖𝑛𝑡 −𝑟𝑒𝑝 _𝑛𝑠 𝑕𝑟 100 ∗ 𝑠𝑐 (18) 𝑅2 = 𝑛𝑔𝑕 + 𝑛𝑔𝑕 ∗ 𝑘𝑒𝑒 𝑝 𝑝𝑜𝑖𝑛𝑡 −𝑟𝑒𝑝 _𝑛𝑒𝑛 𝑕 100 ∗ 𝑠𝑐 (19) VI. SIMULATION RESULTS The proposed Improved Bees Algorithm (IBA) algorithm for solving ORPD problem is tested for standard IEEE-57 bus power system. The IEEE 57-bus system data consists of 80 branches, seven generator-buses and 17 branches under load tap setting transformer branches. The possible reactive power compensation buses are 18, 25 and 53. Bus 2, 3, 6, 8, 9 and 12 are PV buses and bus 1 is selected as slack-bus. In this case, the search space has 27 dimensions, i.e., the seven generator voltages, 17 transformer taps, and three capacitor banks. The system variable limits are given in Table I. The initial conditions for the IEEE-57 bus power system are given as follows: Pload = 12.310 p.u. Qload = 3.322 p.u. The total initial generations and power losses are obtained as follows: 𝑃𝐺 = 12.7634 p.u. 𝑄 𝐺 = 3.4468 p.u. Ploss = 0.27351 p.u. Qloss = -1.2248 p.u.
  • 6. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 39 Paper Publications Table II shows the various system control variables i.e. generator bus voltages, shunt capacitances and transformer tap settings obtained after IBA based optimization which are within their acceptable limits. In Table III, a comparison of optimum results obtained from proposed IBA with other optimization techniques for ORPD mentioned in literature for IEEE-57 bus power system is given. These results indicate the robustness of proposed IBA approach for providing better optimal solution in case of IEEE-57 bus system. TABLE I: VARIABLES LIMITS FOR IEEE-57 BUS POWER SYSTEM (P.U.) REACTIVE POWER GENERATION LIMITS BUS NO 1 2 3 6 8 9 12 QGMIN -1.2 -.014 -.02 -0.06 -1.2 -0.03 -0.3 QGMAX 2 0.4 0.5 0.24 2 0.08 1.54 VOLTAGE AND TAP SETTING LIMITS VGMIN VGMAX VPQMIN VPQMAX TKMIN TKMAX 0.7 1.3 0.95 1.06 0.7 1.3 SHUNT CAPACITOR LIMITS BUS NO 18 25 53 QCMIN 0 0 0 QCMAX 10 5.3 6.5 TABLE II: CONTROL VARIABLES OBTAINED AFTER OPTIMIZATION BY IBA METHOD FOR IEEE-57 BUS SYSTEM (P.U.) Control Variables IBA V1 1.2 V2 1.084 V3 1.073 V6 1.051 V8 1.074 V9 1.052 V12 1.061 Qc18 0.0843 Qc25 0.333 Qc53 0.0628 T4-18 1.016 T21-20 1.072 T24-25 0.973 T24-26 0.945 T7-29 1.092 T34-32 0.957 T11-41 1.015 T15-45 1.074 T14-46 0.943 T10-51 1.055 T13-49 1.075 T11-43 0.921 T40-56 0.911 T39-57 0.973 T9-55 0.985
  • 7. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 40 Paper Publications TABLE III: COMPARATIVE OPTIMIZATION RESULTS FOR IEEE-57 BUS POWER SYSTEM (P.U.) S.No. Optimization Algorithm Best Solution Worst Solution Average Solution 1 NLP [29] 0.25902 0.30854 0.27858 2 CGA [29] 0.25244 0.27507 0.26293 3 AGA [29] 0.24564 0.26671 0.25127 4 PSO-w [29] 0.24270 0.26152 0.24725 5 PSO-cf [29] 0.24280 0.26032 0.24698 6 CLPSO [29] 0.24515 0.24780 0.24673 7 SPSO-07 [29] 0.24430 0.25457 0.24752 8 L-DE [29] 0.27812 0.41909 0.33177 9 L-SACP-DE [29] 0.27915 0.36978 0.31032 10 L-SaDE [29] 0.24267 0.24391 0.24311 11 SOA [29] 0.24265 0.24280 0.24270 12 LM [30] 0.2484 0.2922 0.2641 13 MBEP1 [30] 0.2474 0.2848 0.2643 14 MBEP2 [30] 0.2482 0.283 0.2592 15 BES100 [30] 0.2438 0.263 0.2541 16 BES200 [30] 0.3417 0.2486 0.2443 17 Proposed IBA 0.22359 0.23492 0.23121 VII. CONCLUSION IBA has been fruitfully applied for ORPD problem. The IBA based ORPD is tested in standard IEEE-57 bus system. Performance comparisons with well-known population-based algorithms give cheering results. IBA emerges to find good solutions when compared to that of other algorithms. The simulation results presented in previous section prove the ability of IBA approach to arrive at near global optimal solution. REFERENCES [1] O.Alsac,and B. Scott, “Optimal load flow with steady state security”,IEEE Transaction. PAS -1973, pp. 745-751. [2] Lee K Y ,Paru Y M , Oritz J L –A united approach to optimal real and reactive power dispatch , IEEE Transactions on power Apparatus and systems 1985: PAS-104 : 1147-1153 [3] A.Monticelli , M .V.F Pereira ,and S. Granville , “Security constrained optimal power flow with post contingency corrective rescheduling” , IEEE Transactions on Power Systems :PWRS-2, No. 1, pp.175-182.,1987.
  • 8. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 41 Paper Publications [4] Deeb N ,Shahidehpur S.M ,Linear reactive power optimization in a large power network using the decomposition approach. IEEE Transactions on power system 1990: 5(2) : 428-435 [5] E. Hobson ,‟Network consrained reactive power control using linear programming, „ IEEE Transactions on power systems PAS -99 (4) ,pp 868=877, 1980 [6] K.Y Lee ,Y.M Park , and J.L Oritz, “Fuel –cost optimization for both real and reactive power dispatches” , IEE Proc; 131C,(3), pp.85-93. [7] M.K. Mangoli, and K.Y. Lee, “Optimal real and reactive power control using linear programming” , Electr.Power Syst.Res, Vol.26, pp.1-10,1993. [8] C.A. Canizares , A.C.Z.de Souza and V.H. Quintana , “ Comparison of performance indices for detection of proximity to voltage collapse ,‟‟ vol. 11. no.3 , pp.1441-1450, Aug 1996 . [9] S.R.Paranjothi ,and K.Anburaja, “Optimal power flow using refined genetic algorithm”, Electr.Power Compon.Syst , Vol. 30, 1055-1063,2002. [10] D. Devaraj, and B. Yeganarayana, “Genetic algorithm based optimal power flow for security enhancement”, IEE proc-Generation.Transmission and. Distribution; 152, 6 November 2005. [11] A. Berizzi, C. Bovo, M. Merlo, and M. Delfanti, “A ga approach to compare orpf objective functions including secondary voltage regulation,” Electric Power Systems Research, vol. 84, no. 1, pp. 187 – 194, 2012. [12] C.-F. Yang, G. G. Lai, C.-H. Lee, C.-T. Su, and G. W. Chang, “Optimal setting of reactive compensation devices with an improved voltage stability index for voltage stability enhancement,” International Journal of Electrical Power and Energy Systems, vol. 37, no. 1, pp. 50 – 57, 2012. [13] P. Roy, S. Ghoshal, and S. Thakur, “Optimal var control for improvements in voltage profiles and for real power loss minimization using biogeography based optimization,” International Journal of Electrical Power and Energy Systems, vol. 43, no. 1, pp. 830 – 838, 2012. [14] B. Venkatesh, G. Sadasivam, and M. Khan, “A new optimal reactive power scheduling method for loss minimization and voltage stability margin maximization using successive multi-objective fuzzy lp technique,” IEEE Transactions on Power Systems, vol. 15, no. 2, pp. 844 – 851, may 2000. [15] W. Yan, S. Lu, and D. Yu, “A novel optimal reactive power dispatch method based on an improved hybrid evolutionary programming technique,” IEEE Transactions on Power Systems, vol. 19, no. 2, pp. 913 – 918, may 2004. [16] W. Yan, F. Liu, C. Chung, and K. Wong, “A hybrid genetic algorithminterior point method for optimal reactive power flow,” IEEE Transactions on Power Systems, vol. 21, no. 3, pp. 1163 –1169, aug. 2006. [17] J. Yu, W. Yan, W. Li, C. Chung, and K. Wong, “An unfixed piecewiseoptimal reactive power-flow model and its algorithm for ac-dc systems,” IEEE Transactions on Power Systems, vol. 23, no. 1, pp. 170 –176, feb. 2008. [18] F. Capitanescu, “Assessing reactive power reserves with respect to operating constraints and voltage stability,” IEEE Transactions on Power Systems, vol. 26, no. 4, pp. 2224–2234, nov. 2011. [19] Z. Hu, X. Wang, and G. Taylor, “Stochastic optimal reactive power dispatch: Formulation and solution method,” International Journal of Electrical Power and Energy Systems, vol. 32, no. 6, pp. 615 – 621, 2010. [20] A. Kargarian, M. Raoofat, and M. Mohammadi, “Probabilistic reactive power procurement in hybrid electricity markets with uncertain loads,” Electric Power Systems Research, vol. 82, no. 1, pp. 68 – 80, 2012. [21] Pham, D.T.; Ghanbarzadeh, A.; Koc, E.; Otri, S.; Rahim, S.; Zaidi, M. The Bees Algorithm, Technical Note; Manufacturing Engineering Center, Cardiff University: Cardiff, UK, 2005. [22] Seeley, T.D. The Wisdom of the Hive: The Social Physiology of Honey Bee Colonies; Harvard University Press: Cambridge, MA, USA, 2009. [23] Gould, J.L.; Gould, C.G. The Honey Bee; Scientific American Library: New York, NY, USA, 1988. [24] Von Frisch, K. Bees: Their Vision, Chemical Senses, and Language; Cornell University Press: Ithaca, NY, USA, 1950. [25] Huang, Z. Behavioral Communications: The Waggle Dance. Available online: http://photo.bees.net/ biology/ch6/dance2.html (accessed on 29 June 2013). [26] Koc, E. The Bees Algorithm Theory, Improvements and Applications. Ph.D Thesis, Cardiff University, Cardiff, UK, 2010.
  • 9. International Journal of Recent Research in Electrical and Electronics Engineering (IJRREEE) Vol. 1, Issue 1, pp: (34-42), Month: April - June 2014, Available at: www.paperpublications.org Page | 42 Paper Publications [27] Yuce, B. Novel Computational Technique for Determining Depth Using the Bees Algorithm and Blind Image Deconvolution. Ph.D Thesis, Cardiff University, Cardiff, UK, 2012. [28] Baris Yuce, Michael S. Packianather , Ernesto Mastrocinque , Duc Truong Pham and Alfredo Lambiase “Honey Bees Inspired Optimization Method” : The Bees Algorithm Insects 2013, 4, 646-662; doi:10.3390/insects4040646 [29] Chaohua Dai, Weirong Chen, Yunfang Zhu, and Xuexia Zhang, “Seeker optimization algorithm for optimal reactive power dispatch,” IEEE Trans. Power Systems, Vol. 24, No. 3, August 2009, pp. 1218-1231. [30] J. R. Gomes and 0. R. Saavedra, “Optimal reactive power dispatch using evolutionary computation: Extended algorithms,” IEE Proc.-Gener. Transm. Distrib.. Vol. 146, No. 6. Nov. 1999. Author Biography: K. Lenin has received his B.E., Degree, electrical and electronics engineering in 1999 from university of madras, Chennai, India and M.E., Degree in power systems in 2000 from Annamalai University, Tamil Nadu, India. At present pursuing Ph.D., degree at JNTU, Hyderabad, India. Bhumanapally. RavindhranathReddy, Born on 3rd September,1969. Got his B.Tech in Electrical & Electronics Engineering from the J.N.T.U. College of Engg., Anantapur in the year 1991. Completed his M.Tech in Energy Systems in IPGSR of J.N.T.University Hyderabad in the year 1997. Obtained his doctoral degree from JNTUA,Anantapur University in the field of Electrical Power Systems. Published 12 Research Papers and presently guiding 6 Ph.D. Scholars. He was specialized in Power Systems, High Voltage Engineering and Control Systems. His research interests include Simulation studies on Transients of different power system equipment. M. Surya Kalavathi has received her B.Tech. Electrical and Electronics Engineering from SVU, Andhra Pradesh, India and M.Tech, power system operation and control from SVU, Andhra Pradesh, India. she received her Phd. Degree from JNTU, hyderabad and Post doc. From CMU – USA. Currently she is Professor and Head of the electrical and electronics engineering department in JNTU, Hyderabad, India and she has Published 16 Research Papers and presently guiding 5 Ph.D. Scholars. She has specialised in Power Systems, High Voltage Engineering and Control Systems. Her research interests include Simulation studies on Transients of different power system equipment. She has 18 years of experience. She has invited for various lectures in institutes.