Fuzzy Multiobjective Optimal Power Flow Based on Modified Artificial Bee Colony Algorithm

This paper presents amodified artificial bee colony (MABC) algorithm to solve optimal power flow (OPF) problem. In the proposed MABC algorithm, the searching operation for new food source of artificial bee colony (ABC) algorithm is replaced with mutation and crossover operation of differential evolution (DE) algorithm to improve exploitation capacity. The OPF objective functions involve minimization of total fuel cost of generating units, minimization of emission of atmospheric pollutants, minimization of active power losses, and minimization of voltage deviations. The fuzzy satisfaction-maximizing method is utilized to convert the multiobjectives problem into single objective problem. The proposed approach is applied to the OPF problem on IEEE 30-bus test system. And the results are compared with those obtained by other heuristic algorithms, which demonstrate that the MABC algorithm not only has a better exploration capacity but also possesses stronger exploitation capacity and can effectively solve the OPF problem.


Introduction
Optimal power flow (OPF) was first proposed by French scholar Carpentier in the 1960s.Its definition can be described as the optimal power flow distribution of power system with fixed structure parameters and loads, which optimizes the objective functions by optimal settings of the OPF control variables, while at the same time satisfying various constraints [1,2].OPF problem has received much attention by the academics due to the fact that it is an effective analysis tool for the safe operation and economic dispatch of power system.
On the aspect of choosing OPF objective function, Paranjothi and Anburaja in [3] established the single objective OPF model which considered the minimization of total fuel cost as objective function.The results illustrated that the economics of power system could be further improved in some extent.Because the active power loss caused by transmission and distribution and the emission cost brought by the some generators affect the economics of the power system [4], the OPF objective functions which are established from the aspect of system economics also need to consider these two factors.In order to further improve the system stability, Pouyan et al. in [5] proposed the OPF objective function which furthermore considered voltage stability index, and the obtained optimization scheme considers both economics and stability of power system.The multiobjective OPF model can comprehensively reflect the operation performance of the optimized power system, which improves not only the economy of power system but also the stability of operation.
The weighting method [6] and fuzzy mathematics method [7] et al. are usually used to handle each objective weight of multiobjective optimization problem.The weighting method adopts weighting coefficients to treat each objective function to form a single objective.This method can solve multiobjective OPF effectively, but the weighting coefficients are determined either by subjective settings or by several tests for optimal scheme to get a more satisfactory function value.The fuzzy mathematics method uses membership function to make objective function fuzzification to form the fuzzy multiobjective OPF.There is no need to use weighting 2 Mathematical Problems in Engineering coefficients to solve multiobjective problem in fuzzy mathematics method, and the solution is obtained objectively.
The OPF problem is considered as complex multiconstraints, nonlinear, and noncontinuous optimization problem regardless of OPF model with single objective or multiobjective function.Many classical algorithms such as linear programming [8], quadratic programming [9], simplified gradient method [10], Newton method [11], and interior point method [12] have been widely applied to solve the OPF problem.These algorithms are taking an order or second order gradient of objective functions as the main information to search for the optimal solution and have relatively rapid calculation.They are suitable for online calculation but not suitable for optimization problem with discrete variables.The heuristic algorithms are novel algorithms for solving the optimization problem.And the typical heuristic methods include genetic algorithm (GA), particle swarm optimization (PSO), differential evolution (DE), artificial immune algorithm, and artificial bee colony (ABC) algorithm [13][14][15][16][17].These algorithms are based on multipoint stochastic searching, and they can effectively solve optimization problem with discrete variables.Their global convergence capacity is better than classical algorithms.
The ABC algorithm possesses superior performance on many benchmark functions compared to other heuristic algorithms [18].It was proposed by Karaboga in 2005 and it simulated the intelligent foraging behavior of honeybee swarms [19].Its advantages include easy implementation and better global search ability.The ABC algorithm was applied to the power system network reconfiguration problem which considered the power loss as the optimization objective function for the first time in the literature [20] and was simulated on test system.The results show that the power loss of reconfiguration system by the ABC algorithm is less and running time is shorter than other algorithms.In [21], the authors utilized the ABC algorithm to solve the economic power dispatch problem, and the proposed method was applied to the test system.The results indicate that the exploration capacity of ABC algorithm was better than other heuristic algorithms.
The ABC algorithm is also easily to get trapped in local optima as a heuristic algorithm [22].In order to overcome the shortage of ABC algorithm, the scholars have proposed some improved ABC algorithms [23][24][25].Gao and Liu in [26] proposed the new solution search equation of ABC algorithm which was inspired by the DE/best/1 algorithm and adopted the chaotic systems and opposition-based learning method to generate the initial population.In [27] some modified search equations were proposed and the orthogonal learning strategy was applied to discover the useful information from the search experiences to improve the optimization performance of ABC algorithm.However, the exploration capacity and exploitation capacity are not very well balanced in some improved algorithms.For comprehensively improving the exploration and exploitation capacity, the paper proposes a modified ABC algorithm named MABC algorithm whose searching operation for new food source of ABC algorithm is replaced with mutation and crossover operation of DE algorithm and uses the MABC algorithm to solve OPF problem.For comprehensively considering the economics and stability of the power system, the OPF objective functions involve minimization of total fuel cost, minimization of emission, minimization of active power loss, and minimization of voltage deviation.The fuzzy satisfaction-maximizing method is utilized to convert the multiobjective problem into single objective problem to form fuzzy multiobjective OPF.Finally, the fuzzy multiobjective OPF based MABC algorithm is tested on the standard IEEE 30-bus test systems.The simulation results are compared with those obtained by other heuristic methods, and they demonstrate that MABC algorithm possesses not only better exploration capacity but also stronger exploitation capacity.

OPF Mathematical Model
The main purpose of the OPF is to determine the optimal control variables values for minimizing one or more objective functions while satisfying the several equality and inequality constraints.Its mathematical model [28] can be established as min  (, ) where  is the objective function to be optimized,  is the equality constraints representing nonlinear power flow equations, and ℎ is the system operating constraints;  is the vector of dependent state variables including slack bus generated active power, generator reactive power output, load () bus voltage, and transmission line flow;  is the vector of independent control variables including generator active power output, generator bus voltage, transformer tap settings, and shunt VAR compensation.The paper considers four optimization objectives of OPF.Their mathematical models and OPF constraints are described as follows.

Minimization of Total Fuel Cost.
The fuel cost of generating units is one of the important factors affecting the system economic.And how to allocate the active power output at each generating unit decides the total fuel cost under the conditions of fixed demand active power.Therefore, the total fuel cost is considered as one of the OPF objective functions, and it can be expressed as where FC is the total fuel cost;   ,   , and   are the fuel cost coefficients of the th generator;   is the active power output of thermal unit ;   is the number of generators.

Minimization of Total Emission.
The total emission of atmospheric pollutants is caused by the operation of fossil-fueled thermal generation.And the treatment cost of pollutants will be produced, which have an influence on the economical operation of power system.Therefore, the total emission is also considered as one of objective functions, and it can be described as where EM is the total emission;   ,   , and   are the emission coefficients of the th generator.

Minimization of Total Power Loss.
Active power loss is produced by the transmission and distribution in power system.It is an important measurable index of the economic operation.Smaller active power loss will generate the better economic benefits; therefore, the total power loss is considered as one of objective functions and it can be expressed as where  loss is active power loss;   ,   are the voltage magnitudes at bus  and bus , respectively;   ,   are the branch conductance and phase angles difference between bus  and bus , respectively;  is the total number of nodes.

Minimization of Voltage
Deviation.The stability of power system is considered as one of the optimization objectives besides the economics of system.The bus voltage deviation is an important index of voltage stability, and it can measure size of the deviations from 1.0 per unit and represent the overall level of voltage stability.This objective function can be described as where   is the number of  nodes and 1.0 is per unit.

OPF Constraints.
It can be seen that the control variables and state variables must be satisfied with some constraints to make the power system on the stability operation from the definition of OPF.The OPF constraints of power system include equality constraints and inequality constraints.The equality constraints are the power flow equations, and they can be described as where   ,   are active power and reactive power outputs of generating units at bus , respectively;   ,   are the demanded active power and reactive power of loads at bus , respectively;   ,   , and   are transfer conductance, susceptance, and voltage angle difference between bus  and bus , respectively;   ,   are voltage magnitude at bus  and bus , respectively;  is total number of nodes.
The inequality constraints include control variables constraints and state variables constraints in the OPF model.And the control variables constraints can be described as where   ,   ,   , and   are active power output of the generator, terminal voltage of the generator, transformer tap settings, and reactive power compensation capacity, respectively.They are automatically satisfied by setting the constraint boundary.
The state variables constraints can be described as where   ,   , and   are the load bus voltage magnitudes, reactive power output of generator, and transmission line flow, respectively.They are satisfied by penalty functions that are brought in the objective functions.

Fuzzy Multiobjective OPF
The multiobjective OPF model is established, which involves the minimization of total fuel cost, the minimization of total emission, the minimization of power loss, and the minimization of voltage deviation.For objectively considering the weight of each optimization objective, the linear membership function is used to make the objective functions fuzzification to form the fuzzy multiobjective function OPF.The linear membership function can be described as where   () is th objective function of fuzzy multiobjective OPF;   () is the membership function of   ();  is the number of objective functions;   max ,   min are upper and lower limit values of   (), respectively,   min is the optimal value obtained by single objective OPF for each objective function, and   max is initial value of each objective function.
The curves of each membership function by ( 9) are shown in Figure 1.
Abscissa [4] represents each objective function value of fuzzy multiobjective OPF, and the ordinate represents the membership function value of each optimization objective at (a), (b), (c), and (d) in Figure 1.From (9) or Figure 1, it can be seen that the larger the membership function value, the smaller the corresponding objective function.The objective function   () of fuzzy multiobjective OPF is chosen from membership function   () by the fuzzy satisfactionmaximizing method [29], and it can be expressed as where   () is membership of th objective function;   () is objective function of fuzzy multiobjective OPF;  ∈ {1, 2, . . ., },  is the number of objective functions.The model of fuzzy multiobjective OPF converted by ( 10) is described as min   () The optimal scheme of fuzzy multiobjective OPF described by (11) is the optimal solution of multiobjective OPF, and it can make each objective function achieve satisfactory values at the same time.

Modified Artificial Bee Colony Algorithm
4.1.Artificial Bee Colony Algorithm.Artificial bee colony algorithm is a swarm intelligence optimization algorithm which simulates the behavior of honeybee swarms foraging for maximum nectar amount [18].The honeybee swarms, food sources, and nectar amount are important parts of the ABC algorithm.The food sources foraged by honeybees represent a feasible solution of the optimization problem, and the nectar amount of food sources corresponds to the fitness value of associated feasible solution.The honeybee swarms include employed bees and unemployed bees, and the unemployed bees are divided into two groups: onlooker bees and scout bees.The number of employed bees   is equal to the number of onlooker bees   and is equal to half of the number of honeybee swarms   [19].The whole optimization process of ABC algorithm mainly includes initialization, employed bee phase, onlooker bee phase, and scout bee phase.And they are described by the following subsections.

Initialization.
The ABC algorithm randomly generates   numbers of initial swarms which are described as  = {  |  = 1, 2, . . .,   ;  = 1, 2, . . ., };  is the number of parameter dimensions and is equal to the number of control variables.The initial swarms are generated as where   is the th dimension parameter of food source ;   ,   are the upper and lower bounds for the dimension , respectively; rand is random number between 0 and 1.
The fitness values of the initial swarms are calculated by the following expression: where   is the objective function value of employed bee , and fit  is the fitness value of employed bee .
The higher fitness value indicates the smaller value of objective function and the better food source.The values of initial swarms are sorted from the largest to the smallest.And the corresponding food sources of the previous   numbers fitness values are chosen to be considered as food sources of employed bees.

Employed Bee Phase.
After initialization, the employed bees search for new candidate food sources V  within the neighborhood of the associated food sources   .The searching equation is described as where  ∈ {1, 2, . . ., } and  ∈ {1, 2, . . .,   } are randomly generated and  ̸ = ; V  is a new candidate food source;   is a randomly chosen food source different from   ;   is random number between −1 and 1, if a parameter value exceeds its limits, it will be fixed on its limit value.
The greedy selection mechanism is utilized to select the better solution between V  and   .

Onlooker Bee
Phase.After all employed bees complete the searching process, the probability value that a food source will be chosen by the onlooker bee is calculated by the following expression: where   is the probability value.
Each onlooker bee chooses a food source according to the probability   to search for a new food source by (14) and records the better solution by the greedy selection mechanism.By increasing the fitness value of that food source, the probability of a food source chosen by the onlooker bees increases.

Scout Bee
Phase.After all the employed bees and onlooker bees complete their operation phase, the scout bee checks the numbers of fitness unimproved of a food source.If a food source cannot be improved further over the predetermined number of cycles limit, which indicates that the fitness value of the food source is local optima, the food source will be replaced with a new food source produced by the scout bee using (12).

Modified Artificial Bee Colony
Algorithm.The ABC algorithm possesses the better global searching ability due to the employed bees and scout bee phase.The only one dimension parameter will be changed by (14) for new food sources searching, and it may find the global optimal solution in low dimension optimization problems.But the optimization efficiency will be decreased by (14) in the high dimension optimization problems, and the algorithm gets trapped in the local optimal solution easily.Therefore, the local searching ability of ABC algorithm is poor in the high dimension optimization problems.
The DE algorithm possesses the strong exploitation ability due to one or more dimension parameters which will be changed at the mutation operation process [30].In order to overcome the shortage of poor exploitation capacity of ABC algorithm, the paper proposes the modified ABC algorithm which brings the mutation and crossover operation of DE algorithm into the employed bee phase and onlooker bee phase instead of the searching operation of ( 14).The modified ABC algorithm is named MABC algorithm which combines the advantages of both ABC and DE algorithms.
On the aspect of choosing the mutation operation, the paper considers the influence of current optimal food source  best on searching for new food source   at the th iteration.The convergence speed is improved by the effect of current optimal food source  best .The new food source   of mutation operation is generated by where   is a new food source generated by the mutation operation;  best is the current best food source having the highest fitness value;  1 and  2 are different food sources generated randomly, and 1 ̸ = 2 ̸ = ;  and  are scaling factor of mutation operation.
The new food source V  of crossover operation is generated by where V  is a new food source generated by the crossover operation;   is random parameter of dimension  in the range [0, 1], generated anew for each value of ; CR is the constant parameter of crossover operation in the range [0, 1];  is random parameter in the range [0, ] which ensures that V  gets at least one parameter from   .The flowchart of MABC algorithm for solving the optimization problem is shown in Figure 2.

Start
Generate the initial swarms using (12) calculate the fitness value using (13), and record the better solution by greedy selection mechanism

Memorize the global optimal solution
Print the optimal solution Yes No End Calculate the fitness value using (13) and record the better solution by greedy Whether thel local optimal solution exists Generate a new solution using (12) to replace the local optimal one Yes No Calculate the fitness values using (13), sort them from the largest to the smallest, and choose previous N e numbers of bees as the employed bees Search for new solution i using ( 16) and ( 17), Calculate the probability value P i using ( 15), choose a solution x i by P i to search for new one i selection mechanism between x i and i k > MCN

Fuzzy Multiobjective OPF Base on MABC Algorithm
The proposed MABC algorithm is applied to fuzzy multiobjective OPF in power system.The capacity of optimization to multiobjective functions simultaneously can be reflected by the fuzzy fitness value which is calculated by where   is the objective function value of th food source; fit  is the fuzzy fitness value of th food source;   is the number of food sources.The higher fuzzy fitness value represents the closer to optimal value of each objective function simultaneously and the better optimization capacity.The detailed calculation processes of fuzzy multiobjective OPF based on MABC algorithm are described as follows.
(3) Calculate fuzzy fitness value.Calculate the value of objective functions using (2)∼( 5) based on the results of Newton-Raphson power flow calculation.Utilize (9) to convert each objective function value into corresponding fuzzy membership function values.Choose a fuzzy optimization objective function of multiobjective OPF using (10) according to each membership function value.Finally, calculate the fuzzy fitness value by (18).
(4) Select the employed bees.Sort the   numbers of fuzzy fitness values obtained by the step (3) from the largest to the smallest and select food sources of the previous   numbers of fitness values to be considered as employed bees.
(5) Employed bees phase: each employed bee   adopts ( 16) and ( 17) to search for new food source V  .
(6) Greedy selection mechanism: the fuzzy fitness value fit(V  ) of the new food source V  is calculated by step (3).It is compared with fit(  ) of corresponding food source   ; if the fit(V  ) is better than fit(  ),   will be replaced by the new food source V  ; otherwise   is retained.
(7) Calculate probability value   .The probability value   of each employed bee   is calculated by (15).( 8) Onlooker bees phase: a parameter is randomly generated by th onlooker bee between 0 and 1, and it is compared with   ; if it is smaller than   , then the th onlooker bee searches for a new food source V  within the neighborhood of food source   by ( 16) and (17).
Greedy selection mechanism is used to choose the better feasible solution between V  and   , and  =  + 1.Otherwise,  =  + 1, and the th onlooker bee regenerates a random parameter until it is assigned.
(9) Judge whether all the onlooker bees are assigned.If all the onlooker bees are assigned ( >   ), then proceed to step (10).Otherwise, return to step (8).
(10) Scout bee phase: the scout bee checks the unimproved fitness number; if a food source cannot be improved further over the number limit, the food source will be abandoned and a new food source will be generated by (12) to replace the associated food source.
(11) Record the current optimal food source, current optimal value of each objective function, and current optimal fuzzy fitness value of MABC algorithm, and the cycle parameter  =  + 1.

Numerical Examples
Case 1 (minimization of total fuel cost for single objective OPF).The minimization of total fuel cost described by ( 2) is selected as objective function for single objective OPF.
The minimum total fuel cost obtained by the proposed MABC approach over 20 independent runs was 800.3981 $/h with average of 800.4043 $/h and maximum of 800.4446 $/h, which are compared to those obtained by ABC and also compared to gravitational search algorithm (GSA) [6], enhanced genetic algorithm (EGA) [33], improved evolutionary programming (IEP) [34], and DE [35].The results are given in Table 1.
From Table 1, it can be seen that the average total fuel cost obtained by the proposed MABC approach was 800.4043 $/h which is less than the 800.5346 $/h obtained by ABC algorithm and is smaller than 812.1935 $/h obtained by the GSA [6] obviously.The maximum value was 800.4446 $/h obtained by MABC algorithm, and it is better than maximum values obtained by other algorithms in Table 1.The maximum value, average value, minimum value, and standard deviation value obtained by MABC algorithm are better than those obtained by ABC algorithm.The results demonstrate that the MABC can solve the OPF problem effectively, and the optimal value obtained by MABC algorithm is better than other algorithms.The convergence characteristics of MABC and ABC algorithms for Case 1 are shown in Figure 4.  [36] 0.2096 GA [36] 0.2117 SFLA [36] 0.2063 MSFLA [36] 0.2056 In order to compare the results by MABC with those by ABC obviously, Figure 4 has only described the data of previous 200 iterations.From Figure 4 and simulation data, it can be seen that when the iteration number reaches up to 20, the fuel cost values obtained by MABC and ABC are 800.5810$/h and 815.4091 $/h, respectively.The results demonstrate that the convergent rate of the MABC algorithm is better than that of ABC algorithm.
Case 2 (minimization of total emission for single objective OPF).The minimization of total emission described by ( 3) is selected as objective function for single objective OPF.The minimum total emission obtained by MABC is 0.1943 ton/h, which are compared to those by ABC and also compared to those obtained by particle swarm optimization (PSO) [36], genetic algorithm (GA) [36], shuffled frog leaping algorithm (SFLA) [36], and modified SFLA (MSFLA) [36].The results are given in Table 2.
From Table 2, it can be seen that the optimal value obtained by MABC is 0.1943 ton/h which reduces by 8.22% compared with 0.2117 ton/h obtained by GA [36] and is smaller than the values obtained by other algorithms obviously.The results demonstrate that the optimal value obtained by MABC algorithm is better than other algorithms.The convergence characteristics of the MABC and ABC for Case 2 are shown in Figure 5.
From Figure 5 and simulation data, it can be seen that when the iteration number reaches up to 20, the minimum value obtained by MABC is already equal to global optimal value 0.1943 ton/h, and at the same time the objective function value obtained by ABC algorithm got tripped in local optimal value, and the ABC algorithm converges to global optimal value in about 160th iteration.The results demonstrate that the convergent rate of the MABC algorithm is better than that of ABC obviously.
Case 3 (minimization of active power loss for single objective OPF).The minimization of active power loss described by (4) is selected as objective function for single objective OPF.The minimum active power loss obtained by the MABC is 3.0819 MW, which is compared to the ABC and also compared to those obtained by enhanced genetic algorithm (EGA) [37], PSO [37], and DE [38].The results are given in Table 3.

Method
Active power loss (MW) MABC 3.0819 ABC 3.0938 EGA [37] 3 .2 0 0 8 PSO [37] 3 .6 2 9 4 DE [38] 3 .2 4 0 0 From Table 3, it can be seen that the optimal value obtained by MABC algorithm is 3.0819 MW which reduces by about 15.09% compared with the 3.6294 MW obtained by PSO [37] and reduces by about 4.88% compared with the 3.2400 MW obtained by DE [38] and is smaller than the values obtained by other algorithms.The results demonstrate that the MABC leads to better results than the other algorithms.The convergence characteristics of the MABC and ABC algorithms for Case 3 are shown in Figure 6.
From Figure 6 and simulation data, it can be seen that the MABC algorithm converges to the highest quality solution among the two algorithms in less iterations.
Case 4 (minimization of voltage deviation for single objective OPF).The minimization of voltage deviation described by ( 5) is selected as objective function for single objective OPF.The minimum voltage deviation obtained by MABC algorithm is 0.0841 p.u., which is compared with those obtained by ABC algorithm and also compared to those obtained by DE [15] and biogeography-based optimization algorithm (BBO) [39].The results are given in Table 4.
From Table 4, it can be seen that the optimal value obtained by MABC is 0.0841 p.u. which reduces by 19.13% compared with optimal value 0.1040 p.u. obtained by ABC   [15] 0.1357 BBO [39] 0.1020 algorithm and reduces by 38.03% compared with optimal value 0.1357 p.u. obtained by DE algorithm [15] and reduces by 17.55% compared with optimal value 0.1020 p.u. obtained by BBO algorithm [39].The results indicate that the optimal value found by MABC is better than other heuristic algorithms obviously.The convergence characteristics of the MABC and ABC algorithms for Case 4 are shown in Figure 7.
From Figure 7 and simulation data, it can be seen that the convergent capacity is better than other heuristic algorithms.
From the results of single objective OPF in Case 1∼Case 4, it can be seen that the convergent capacity and the optimal value of MABC algorithm are better than ABC algorithm and other heuristic algorithms, which can demonstrate that the proposed MABC algorithm possesses not only the stronger exploration capacity but also the better exploitation capacity.
Case 5 (fuzzy multiobjective OPF).The optimal values obtained by Case 1∼Case 4 are considered as the reference values of fuzzy membership function, and ( 9) and ( 10) are used to convert OPF with four objective functions into fuzzy multiobjective OPF which is solved by MABC algorithm.The minimum values of each objective function obtained by the proposed MABC approach and ABC algorithm over 20 independent runs for Case 5 are given in Table 5.The maximum value, average value, minimum value, and standard deviation value of the fuzzy fitness obtained by MABC and ABC algorithms over 20 independent runs are given in Table 6.
From Table 6, it can be seen that the maximum value obtained by MABC is 0.7716 which is higher than the value 0.7645 obtained by ABC algorithm.The minimum value and average value obtained by MABC are also higher than those obtained by ABC.The results indicate that the optimal solution obtained by MABC algorithm is better than that by ABC algorithm.In other words, the control variables setting values obtained by MABC are able to make the system more efficient and reliable than those by ABC algorithm.
From Table 6, it also can be seen that the standard deviation value by MABC is 0.0003 which is smaller than the value 0.0044 by ABC obviously.The results demonstrate that the whole optimization capacity of MABC is obviously better than that of ABC algorithm.
The convergence characteristics of MABC and ABC algorithms for Case 5 are shown in Figure 8.
From Figure 8, it can be seen that the MABC algorithm converges to the highest fuzzy fitness value among the two algorithms in less iterations.In addition, the convergent characteristic curve of MABC algorithm is smoother, and the ABC algorithm is easily to get trapped in local optima to solve the multiobjective OPF.The results illustrate that the convergence capacity of proposed MABC algorithm is better than ABC algorithm.
In Case 1∼Case 5, the optimal control variables and the corresponding objective function values obtained by the proposed MABC and standard ABC algorithms are given in Table 7.
From Table 7, it can be seen that the optimal solutions obtained by the proposed MABC algorithm can make the results better than those obtained by the ABC algorithm in different optimization objectives.

Conclusions
The paper proposes the MABC algorithm which adopts the mutation and crossover operation of DE algorithm instead of the searching operation of ABC algorithm.And the MABC algorithm is used to solve the multiobjective OPF problem which considers the total fuel cost of generating units, the emission of atmospheric pollutants, the active power losses, and the voltage deviations as the optimization objective.The simulation results show that the convergent capacity and the optimal value of MABC algorithm are better than ABC algorithm and other heuristic algorithms, which demonstrate that the proposed MABC algorithm possesses not only the stronger exploration capacity but also the better exploitation capacity and indicate the efficiency and superiority of MABC algorithm for solving the nonlinear optimization problem.However, the influences of algorithm parameters settings on the optimization performance of MABC algorithm will be further discussed in the future researches.

Figure 4 :
Figure 4: Comparisons of convergent characteristic obtained by different algorithms for Case 1.

Figure 5 :
Figure 5: Comparisons of convergent characteristic obtained by different algorithms for Case 2.

Figure 6 :
Figure 6: Comparisons of convergent characteristic obtained by different algorithms for Case 3.

Figure 7 :
Figure 7: Comparisons of convergent characteristic obtained by different algorithms for Case 4.

Figure 8 :
Figure 8: Comparisons of convergent characteristic obtained by different algorithms for Case 5.

Table 1 :
Comparisons of the results obtained by different algorithms for Case 1.
[14]IEEE 30-Bus TestSystem.The IEEE 30-bus test system[14]has six generators, four transformers, and nine shunt VAR compensation devices and has a total of twenty-four

Table 2 :
Comparisons of the results obtained by different algorithms for Case 2.

Table 3 :
Comparisons of the results obtained by different algorithms for Case 3.

Table 4 :
Comparisons of the results obtained by different algorithms for Case 4.

Table 5 :
Minimum values obtained by different algorithms over 20 independent runs.

Table 6 :
Fuzzy fitness values obtained by different algorithms over 20 independent runs.

Table 5 ,
it can be seen that the minimum values of each objective function obtained by MABC algorithm are 850.5996$/h, 0.2224 ton/h, 4.8331 MW, and 0.3253 p.u., respectively.And the minimum values obtained by ABC algorithm are 852.5226$/h, 0.2224 ton/h, 4.9129 MW, and 0.3474 p.u., respectively.The results demonstrate that the minimum values of each objective function obtained by MABC algorithm are better than those obtained by ABC algorithm for fuzzy multiobjective OPF.

Table 7 :
Comparison of the OPF control variables for different cases.