Engineering, Technology & Applied Science Research Vol. 8, No. 5, 2018, 3321-3328 3321 www.etasr.com Marouani et al.: A Modified Artificial Bee Colony for the Non-Smooth Dynamic … A Modified Artificial Bee Colony for the Non- Smooth Dynamic Economic/Environmental Dispatch I. Marouani National School of Engineering of Sfax (ENIS), University of Sfax, Sfax, Tunisia ismailmarouani@yahoo.fr A. Boudjemline University of Hail, Hail, Saudi Arabia a.boudjemline@uoh.edu.sa T. Guesmi College of Engineering, University of Hail, Hail, Saudi Arabia tawfiq.guesmi@gmail.com H. H. Abdallah National School of Engineering of Sfax (ENIS), University of Sfax, Sfax, Tunisia Hsan.Hadj@enis.rnu.tn Abstract—This paper presents an improved artificial bee colony (ABC) technique for solving the dynamic economic emission dispatch (DEED) problem. Ramp rate limits, valve-point loading effects and prohibited operating zones (POZs) have been considered. The proposed technique integrates the grenade explosion method and Cauchy operator in the original ABC algorithm, to avoid random search mechanism. However, the DEED is a multi-objective optimization problem with two conflicting criteria which need to be minimized simultaneously. Thus, it is recommended to provide the best solution for the decision-makers. Shannon’s entropy-based method is used for the first time within the context of the on-line planning of generator outputs to extract the best compromise solution among the Pareto set. The robustness of the proposed technique is verified on six- unit and ten-unit system tests. Results proved that the proposed algorithm gives better optimum solutions in comparison with more than ten metaheuristic techniques. Keywords-evolutionary computation; power generation dispatch; optimal scheduling; decision making; cost function I. INTRODUCTION Emission dispatch aims at minimizing emission of harmful gases, caused by fossil-fueled thermal units, such as CO, CO2, NOx and SO2 [1-2]. The combination of the above problems is called the economic emission dispatch (EED) problem. However, due to the dynamic nature of today’s network loads, it is required to schedule the thermal unit outputs in real time according to the variation of power demands during a certain time period [3]. To solve this modified EED problem, known as dynamic economic emission dispatch (DEED), several mathematical formulations have been suggested [3-6]. Usually the DEED problem is considered as a dynamic optimization problem having the same objectives as the EED over a time- period of one day, subdivided into definite time intervals of one hour with respect to the constraints imposed by the generator ramp-rate limits (RRL) [3]. Therefore, the operational decision at an hour may be influenced by the one taken at a previous hour. Other constraints such as prohibited operating zones (POZ) and valve-point loading effects (VPLE) have been considered [7-8]. However, incorporating VPLE in the fuel cost function introduces ripples in the latter and the problem will be with multiple minima. On the other hand, POZ constraints due to physical operation limitation like vibrations in the shaft bearing [9] create discontinuities in the objective functions. Therefore, the DEED becomes a highly nonlinear problem with non-convex and discontinuous fitness functions. Classical methods like, dynamic programming [10] and linear programming [11], have been used to solve the static EED. However, these techniques are iterative and require an initialization step, which can cause the convergence of the search process into local optima. Moreover, they may fail to solve the dynamic case including the above constraints. Recently, metaheuristic search algorithms have demonstrated good performance and high efficiency when applied to complex optimization problems. These optimization procedures are classified into various groups in terms of the optimization methodology. Swarm intelligence-based evolutionary algorithms are the most used algorithms. Among metaheuristic-based optimization techniques, genetic algorithm [12], particle swarm optimization [13], simulated annealing [14], artificial bee colony (ABC) [7], tabu search [15], differential evolution [4] and bacterial foraging [5] have been suggested for solving the EED problem. Despite the fact that these techniques have been proven to have a clear edge over traditional methods, they have been criticized [16] because their efficiencies are sensitive to the form of problem constraints and number of units. Most of the above-mentioned works have concentrated only on the static EED problem. Only a few considered the multi-objective DEED problem. In addition, RRL and POZ constraints were not considered during the transition from the last hour of the current day to the first hour of the next day. ABC algorithms attracted much attention for EED problems [9]. ABC algorithm [17] simulates the foraging behavior of a real bee colony for maximizing the nectar amount stocked in the hive. Compared to several population-based techniques like, PSO and GA, the ABC algorithm is simple in concept with a few setting parameters, easy for combination with other optimization approaches and more effective. Unfortunately, like other evolutionary algorithms, the ABC method has also been criticized for its poor convergence rate and premature convergence due to the unbalanced exploration-exploitation processes [16]. Exploration corresponds to the capability to avoid convergence toward local optima by expanding the search into new areas, while exploitation is the capability to Engineering, Technology & Applied Science Research Vol. 8, No. 5, 2018, 3321-3328 3322 www.etasr.com Marouani et al.: A Modified Artificial Bee Colony for the Non-Smooth Dynamic … search the space to find interesting new solutions. Thus, many improved techniques have been proposed to further enhance its performance [16-17]. In order to increase ABC's exploitation ability, the classical employed bee and onlooker bees’ phases have been changed by incorporating the so-called grenade explosion method (GEM) [18]. The GEM is an optimization method proposed in [19] and, as the name suggests, it imitates the mechanism of a grenade explosion. The effectiveness of this modified version of ABC, symbolized by GABC, was verified on a set of standard reference functions. In [16], the Cauchy operator is embedded into the scout bees' phase of the ABC in order to increase the exploration ability by generating a larger set of solutions instead of one random solution for each scout bee. A new method exploiting the advantages of GEM and Cauchy operator is proposed in the present work in order to solve the DEED problem with respect to all the above- mentioned constraints. This optimization, symbolized by GCABC, integrates the GEM and Cauchy operator into the ABC technique. In addition, a new decision-making method based on Shannon’s entropy, called extended entropy-weighted reference (EEWR) approach, is developed and incorporated into the GCABC algorithm to select the most suitable solution among all non-dominated solutions provided by the optimization algorithm. Unlike other techniques like those based on graph theory [20] and Z-transformation [21], the EEWR is characterized by uncomplicated mathematics [22]. The main contributions of the current work are:  A new optimization technique, called GCABC, for scheduling power production of thermal units according to the expected load variations is proposed. To the best of our knowledge, the present work is the first attempt to solve the EED problem by the use of GCABC algorithm. In addition, a modified EWR-based technique, called extended entropy- weighted reference (EEWR), is proposed for decision making. This technique has not been used in any field of power systems.  All aforementioned constraints are considered simultaneously in the DEED problem.  The RRL constraints are taken into account during the transition from the last hour of the one day to the first hour of the next day. II. PROBLEM FORMULATION The DEED problem is considered as a multi-objective optimization problem (MOP). It aims to minimize simultaneously total fuel cost and total emission by finding the power production of thermal power plants according to the predicted load demands. The resolution of the DEED problem can be accomplished by solving the static EED (SEED) problem over a certain period of time subdivided into smaller time intervals. DEED problem objectives and constraints are described below. A. Objective Functions Thermal units with multi-steam admission valves that work sequentially to cover the demand, involves higher order nonlinearity to total fuel cost due to VPLE, as illustrated in Figure 1. Unfortunately, neglecting the VPLE, which is required when using classical methods, causes some inaccuracy in the solution of the DEED problem. Taking into account VPLE constraints, a sinusoidal form is included in the total non-smooth cost function expressed in $/h, as given in (1). The second objective corresponding to the total emission in ton/h is described by (2):      2 min 1 1 sin        T N t t t T i i i i i i i i i t i C a b P c P d e P P (1)    2 1 1 exp       T N t t t T i i i i i i i i t i E P P P     (2) where, ia , ib , ic , id and ie are the cost coefficients of the i-th unit. While, i , i , i , i and i are the emission coefficients. tiP is the output power in MW at the t-th interval. T is the number of hours. In this study, 24T  . The bi-objective DEED problem is converted into a mono- objective optimization problem in [23]. In this study, the price penalty factor (PPF)-based method is adopted. Thus the combined economic-emission objective function FT can be described by (3):  1  T T TF C E   (3) where,  0,1rand  . For each generated value of , the function FT is minimized to obtain the optimum solution that can be a candidate solution to be in the Pareto front. The parameter  is the average of the PPF of all thermal units. As shown in (4), the PPF of the i-th unit is the ratio between its fuel cost, maxi C , and its emission, maxi E , for maximum generation capacity. max max  i i i C PPF E (4) Fig. 1. Fuel cost function with five valves (A, B, C, D, E) B. Problem Constraints The DEED problem is solved by minimizing the function FT defined in (3) with respect to the following constraints:  Generation capacity The real power output of each unit i should be within its minimum and maximum limits miniP and max iP : Generation (MW) F u el c o st ( $ /h ) Without VPLE With VPLE E C A B D Engineering, Technology & Applied Science Research Vol. 8, No. 5, 2018, 3321-3328 3323 www.etasr.com Marouani et al.: A Modified Artificial Bee Colony for the Non-Smooth Dynamic … min max , 1, ,   ti i iP P P i N (5)  Power balance constraints At each time period t, the total power generation must cover the total demand power tDP plus the total transmission losses t LP . Thus, the power balance constraints can be described by: 1 0, 1,...,      N t t t i D L i P P P t T (6) where tLP can be calculated using the constant-loss formula [3] given by (7): 1 1 1       N N N t t t t L i ij j oi i oo i j i P P B P B P B (7) where, ijB , oiB , ooB are the loss parameters also called B- coefficients.  Generating unit RRL In practice, power generation of each unit i during two consecutive time periods is limited by its RRLs defined by (8) and (9): 1  t t downi i iP P R (8) 1  upt ti i iP P R (9) where, 1tiP  is the previous output real power of the i-th machine. downiR and up iR are the down-ramp and up-ramp limits of the i-th unit in (MW/time period). As one of the contributions of the present work, the RRL constraints are taken into account during the transition from the last hour of the current day to the first hour of the next day. Two constraints are embedded in the problem formulation and they are described by (10) and (11): 24 1  downi i iP P R (10) 1 24  upi i iP P R (11)  POZ constraints The POZ constraints are described as: min ,1 , 1 , max , , 2,...,              i t down i i i upt t down i i ii k i k up t i ii z P P P P P P P k z P P P (12) where, , down i kP and , up i kP are down and up bounds of POZ number k, iz is the number of POZ for the i-th unit due to the vibrations in the shaft or other machine faults. Therefore, the machine has discontinuous input–output characteristics [19]. Figure 2 shows the fuel cost function for a typical thermal unit with POZ constraints. Fig. 2. Cost function for a thermal unit with POZ constraints. By considering the generation capacity, RRL and POZ constraints, the minimum and maximum limits of the power generation tiP of the i-th unit for the period t are modified as:  min 1max , max 1 min , , ,1 min 1max , , , 1 max 1 min , , , min 1max , , , mi                               i t down tP P R Pi i i i upt downP P R Pi i i i upt down tP P R P Pi i i ii k tPi upt downP P R Pi i i i k upt down tP P R P Pi i i ii z max 1n ,                      uptP P Ri i i (13) where 2,..., ik z III. ORIGINAL ABC ALGORITHM OVERVIEW ABC algorithm [17] is an efficient and robust technique for several optimization problems. As all swarm-based techniques, ABC algorithm starts by generating randomly an initial population of SN solutions. Each solution is considered as a food source and it corresponds to an employed bee. SN is half of the entire population size (PN). The onlooker bees constitute the second half. If D is the number of the decision variables, an i-th solution iX will be represented by 1 2 ...    i i i i DX x x x . The fitness function evaluated at the solution iX , signifies the nectar quantity of the corresponding food source estimated by an employed bee:           1 , 0 1 1 , 0        i i i i i f X f Xfit X f X f X (14) where,  if X is the objective function estimated at iX . The probability ip to choose the candidate solution iX by an onlooker bee is expressed as follows. Power output (MW) F u el c o st ( $ /h ) P i min P i,1 down P i,1 up P i,2 down P i,3 down P i,3 upPi,2 up P i max POZ 3 POZ 2 POZ 1 Engineering, Technology & Applied Science Research Vol. 8, No. 5, 2018, 3321-3328 3324 www.etasr.com Marouani et al.: A Modified Artificial Bee Colony for the Non-Smooth Dynamic …     1   i i SN n n fit X p fit X (15) The onlooker bees update the selected food source iX to discover a new one. The new solution iV is generated by modifying only one parameter ijx of iX as:    i i i i kj j j j jv x x x (16) where, indices k and j are chosen randomly from  1, 2,..., SN and  1, 2,..., D , respectively. The index k must be different from i. The parameter ij is a real number   0,1 generated from the uniform distribution. During the onlooker phase, a greedy selection between food sources iV and iX will be done. If an employed bee’s food source cannot be improved through a pre-specified triggering threshold, called LIMIT, it becomes a scout and its solution will be abandoned. If iX is an abandoned solution, the converted scout-bee starts to search for a new solution randomly according to (17):   min max min0,1  ij j j jx X rand X X (17) where, max jX and min jX are bounds of the food source in dimension j. IV. PROPOSED OPTIMIZATION ALGORITHM An enhanced version of the classical ABC technique is presented to enhance its exploitation and exploration abilities related to the onlooker and scout bees respectively. In [18-19], the ABC method was criticized for the random selection of the j-th dimension because it may slow the convergence of the algorithm and increase the risk of convergence of the search phase to a local optima. To ovoid these limitations, authors in [19] incorporated a new method, called grenade explosion based method (GEM), into the employed bee and onlooker bee phases. The idea behind this technique is to mimic the grenade explosion principle where, the fitness function is the overall damage of the explosion. In each cycle of the proposed method, it is assumed that there is only one grenade with one piece of shrapnel for each decision parameter. A procedure guideline for tuning the parameters of GEM is given in [19]. The shrapnel pieces are thrown in all dimensions to collect information about the area of the explosion which is considered as the old food source as given in (18). Then, a set of new food sources is proposed by the onlooker bees based on the damage- per-shrapnel degree. This allows reaching the global solution more quickly. As given in (19), the optimal search dimension (OSD) of the new candidate iOSDV corresponds to the maximum damage in all directions.    i i i i kt t t t tv x x x (18) where,  1, 2,...,k SN is a randomly chosen index and k i and  1, 2,...,t D .  0,1it  is a random number.     max | 1, 2,..., i iOSD tfit V fit V t D (19) As explained above, a greedy selection between the new solution iV and the old one iX is applied. In order to improve the global and local exploration abilities of the optimization algorithm and ensure the convergence into the global solution within a short calculation time, the Cauchy operator is embedded in the scout bee phase. The incorporation of a Cauchy operator in optimization techniques has been employed in some algorithms to enhance the global search ability [16]. This enhancement is due to the long tail of Cauchy distribution compared to other operators such as Gaussian distribution. As given in Figure 3, the standard Gaussian function has a large probability within the interval [-3,3]. Nevertheless, it is possible to make larger jumps in the search space using Cauchy operator. Fig. 3. Standard Cauchy and Gaussian distributions In this study, the origin-centered Cauchy distribution with unit scale parameter is used. Thus, the new solution provided from an abandoned solution iX will be obtained using (20) instead of (17).  0,1i ij jx x CAUCHY (20) where,    2 1 0,1 1 ( )   ij CAUCHY x (21) Based on the above explanation, the main steps of the modified ABC algorithm are listed below.  Step 1: Initialization Preset population size PN (SN ) Initialize iteration: Iter = 0 Preset Maximum Cycle Number (MCN) Fix triggering threshold (LIMIT) Initialize food sources using (4)  Step 2: Population evaluation Estimate fitness value of each source using (1) Initialize failure counter of each source Xi,   0FC i  Step 3: Cycle increment  Step 4: Employed bees’ phase for i= 1 to SN do Produce new source iOSDV from iX using (19); Estimate the fitness value of iOSDV using (19); -4 -2 0 2 4 0 0.1 0.2 0.3 0.4 x P ro b ab il it y d is tr ib u ti o n Gaussian (0,1) Cauchy (0,1) Engineering, Technology & Applied Science Research Vol. 8, No. 5, 2018, 3321-3328 3325 www.etasr.com Marouani et al.: A Modified Artificial Bee Colony for the Non-Smooth Dynamic … Apply greedy selection between iOSDV and iX ; end for  Step 5: Onlooker bees’ phase t = 0;i = 1; while (t