Multi Objective Dynamic Economic Dispatch with Cubic Cost Functions
Moses Peter Musau, Nicodemus Odero Abungu, Cyrus Wabuge Wekesa
Department of Electrical and Information Engineering, School of Engineering, The University of Nairobi, Nairobi, Kenya
To cite this article:
Moses Peter Musau, Nicodemus Odero Abungu, Cyrus Wabuge Wekesa. Multi Objective Dynamic Economic Dispatch with Cubic Cost Functions. International Journal of Energy and Power Engineering. Vol. 4, No. 3, 2015, pp. 153-167. doi: 10.11648/j.ijepe.20150403.14
Abstract: The formulation and solution of the Dynamic Economic Dispatch (DED) problem is one of the key disciplines in modern power system operation, planning, operation and control. Past researches have considered DED on Quadratic Cost Functions (QCF), with only few works considering higher order cost functions which are more accurate. The Static Economic Dispatch (SED) has been widely tackled in past researches, however, it is the DED problem that represents a real life power system. There is need to review this problem and establish a more practical formulation of the same taking into consideration all the objectives and constraints possible. The methods used in the solution of DED problem have evolved from the traditional deterministic ones, to the pure heuristic, and finally to the state of the heart hybrids. The hybrids methods have been developed to exalt the strengths and improve the weaknesses of the base method. Such optimisation methods need to be reviewed and classified. This paper will do an in depth review of the DED problem on both quadratic and cubic cost functions. Further it will provide a detailed classification of the methods used to solve the problem as its complexity increases. Consequently, three method Hybrids is the way go as far as the solution of the cubic Multi Objective DED (MODED) with five objectives is concerned.
Keywords: Cubic Cost Functions (CCF), Dynamic Economic Dispatch (ded), Hybrid Methods, Multi Objective DED (Moded)
The fuel cost in in existing textbooks, for example Optimization of Power System Operation  is generally assumed to be a smooth Static ED (SED) modeled as a Quadratic Cost Function(QCF).According to Jizhong Zhu, pp85-88, (2006) , this definition however makes many assumptions which are impractical to real systems. These assumptions include ignoring emissions, uncertainties, reactive power dispatch, ramp rates, valve points and integration of Renewable Energy (RE) generators. These assumptions are impractical in the real time power system and can no longer be ignored. Consideration of such assumptions in ED formulation leads to the DED formulation and solution.
DED optimization problem is one of the most important issues which must be taken into consideration in power systems planning and operation. DED is aimed at planning the power output for each devoted generator unit in such a way that the operating cost is minimized and simultaneously, matching real and reactive load demand, power operating limits and above all maintaining the system stability. Further, the security challenges and power trade issues need to be addressed in such a formulation. Based on convention, electrical power systems are operated based on minimizing operational cost while maintaining all the system constraints in place.
The integration of renewable energy into the grid and the need to account for the emission and transmission (real and reactive) losses in a more accurate sense has greatly affected the way in which the DED problem is formulated and solved. This is due to the use of different generators which utilize various fuels and the fact that their real and reactive power outputs are not deterministic but stochastic in such a case. Thus, there has been a desperate need to review the formulation and solution of the existing DED problem.
The Thermal, Renewable Energy (RE), Emissions and Transmission Loss cost functions when considered individually, results in a Single Objective DED (SODED). The SODED has been considered on a Quadratic Cost Function which is less accurate. Where more than two objectives are taken into consideration, a MODED problem results. In most researches, the thermal cost functions has been taken as the base cost function then the other functions and the corresponding constraints have been added depending on the interest of the researcher. However there is need to consider all these objectives and constraints simultaneously since they do affect the modern power system in such a practical manner. The Multi Objective DED(MODED) which is well suited for the modern smart grid need to be revisited to include wind ,solar, emission and transmission line losses in a more accurate scenario in which the five main objectives are handled simultaneously in a fully constrained environment.
Dealing with multiple and conflicting objectives in optimization engineering has always been a challenging undertaking. Many practical applications in optimization often involve many goals to be satisfied simultaneously. In power system planning, operation and control, MODED is a good example of such problems. As the MODED gets more complex, improved and higher order cost functions will continue to replace the Quadratic ones. This is because the higher order cost functions have proved to be more accurate and realistic than the traditional Quadratic ones. In addition, Hybrid methods which do accommodate the uncertain RE will replace the static ones which are restricted to one or two objectives. Also, the existing hybrids must also be tested on higher order cost functions for the MODED with more objectives. Further, three and four method hybrids will be embraced for more accurate realistic results.
Most of the existing security, reactive power and multi area problems of the economic dispatch has been considered in a static environment and using the classical methods. There is need to reconsider these problems in the MODED with the stochastic RE. This will lead to improved ED, increased system security, advanced voltage profile and better trade on electrical power.
This research paper aims at addressing the MODED problem by providing a detailed review of the existing DED works in terms of the problem formulation and the solution methods used. The paper has been divided into six main sections. There is an introduction to power system operation and planning, DED overview, MODED review, formulation of MODED and the existing sub problems, a classification of DED optimization methods and the review conclusion.
2. Modern Power System Operation and Planning
Power system operation and control problem include the Unit Commitment (UC), Optimal Power Flow (OPF), SED, Hydro-Thermal Scheduling (HTS), and DED which are integrated to each other. According to a review of the Science Direct databases, there is not as many studies of the DED problem and it has not been as thoroughly investigated as other electric power system optimization areas of study .This is evident in Table 1. Similar to most of the real-world complex engineering optimization problems, the nonlinear and non-convex characteristics of the cubic cost functions(CCF) are more prevalent in the DED problem more so with the increasing number of objectives and constraints.
Survey of published works dealing with these optimization problems include N.P. Padhy,2004  for the UC, J.A. Momoh et al ,1999  and M.R. Al Rashid and M.E. El-Hawary,2009  for the OPF and A. Mahor et al ,2009  ,H. Altun and T. Yalcinoz ,2008  for the SED.
|Optimization Area||% Publications|
3. Dynamic Economic Dispatch (DED)
Dynamic Economic Dispatch (DED) involves the formulation and solution of the thermal ED such that the ramp rates, valve points and all the possible constraints are accounted for. As evident in Table 1, the DED occupies a prominent place in a power system’s planning, operation and control. The goal of DED is to determine the optimal power outputs of online generating units in order to meet the load demand subject to satisfying various operational equality and inequality constraints over finite dispatch periods. In practice, there are SED and DED problems where the latter considers additional practical constraints. Early research works addressing this aspect were published by T.E.Bechert and H.G. Kwatny, 1972  and T.E. Bechert and Nanming Chen, 1977 .The DED problem has not been addressed and researched as thoroughly as other power optimization problems. Xia et al, 2010,  offered a review of DED highlighting the problem constraints and some available solution methods. The work presented in this research paper is different in that it categorizes earlier work based on the optimization algorithms, number of objective functions involved, number and type of constraints, type of coast function and published works with RE cost functions. Figure 1 illustrates annual published research on DED according to a review of the Science Direct databases .As shown in the figure, the number of publications on DED is on increase. This signifies the undeniable importance of investigating the DED problem.
From the published works sampled, the number of works integrating the Renewable Energy (RE) into the DED problem have been increasing significantly. This means that the inclusion of the RE into the ED problem can no longer be ignored since the security, reactive power and multi area aspects of the power system are affected . This is as shown in the Figure 2.
4. MODED Review
The history of the MODED problem dates back to early 2000 when transmission losses were considered as a cost function other than a constraint leading to the formulation of the first 2-objective MODED
Table 2.2 summarizes the MODED problem so far.The summary includes the number and types of objectives, the type of method used and the number of constraints. From the table, it is clear that there is need to solve the MODED with all the five major objectives considered simultaneously on cubic cost functions. Further the problem need to be constrained further for more realistic results.
|Faisal .A. Mohamed ,2009||4||T,W,S,E||Non Linear Constrained Method (NLCM )||2||Q|
|M.A Abibo, 2003 ||2||T,TL||Strength Pareto Evolutionary Algorithm (SPEA)||3||Q|
|Zhao Bo and Cao Yi-Jia,2005||3||T,E,TL||Multi Objective Particle Swarm Optimization (MOPSO)||5||Q|
|M.A Abibo, 2006||2||T,TL||SPEA and Strength Pareto Genetic Algorithm (SPGA)||4||Q|
|M. Basu, 2006 , Amitah Mahor ,2009,||2||T E||Particle Swarm Optimization (PSO)||3||Q|
|M .Basu, 2007 ||2||T,E||Evolutionary Programming (EP) and Fuzzy Logic (FL)-(EP-FL)||4||Q|
|C .Chen, 2007 ||2||T,W||Simulated Annealing (SA) and Direct Search Method (DSM)-(SA-DSM)||4||Q|
|M.Basu, 2008 ||2||T,E||Elitist Non-Dominated Sorting Genetic Algorithm-II (ENSGA-II)||3||Q|
|R.P Brent,2010 ||2||T,E||Preference-based Non-Dominated Sorting Genetic Algorithm (PNDSGA)||5||Q|
|Zwe-Lee Gaing et al ,2009 ||T,E||Bacterial Foraging(BF)-PSO-Differential Evolution-(DE) (BF-PSO-DE)||4||Q|
|Alsumait et al, 2010 ||2||T,E||Pattern Search (PS)||3||Q|
|S.Subramanial et al, 2010 ||3||T,E,TL||Sequential Approach with Matrix Framework (SAMF)||6||Q|
|Chunghun et al,2013 ||3||T,W,S||Model Predictive Control (MPC)||5||Q|
|Haiwang Zhong et al,2013||2||T,TL||Penalty-Based Algorithm (PBA)||4||Q|
|Azza A. ElDesouky, 2013 ||4||T,W,S,E||Weighted Aggregation (WA) and PSO(WAPSO)||7||Q|
|Bakirtzis et al ,2014 ||2||T,W||Variable Time Resolution (VTR) and Scheduling Horizon (SH)-(VTRSH)||4||Q|
|Soubache ID and Sudhakara Reddy,2014||2||T,E||Multi Objective Computer Programming Based Method(MOCPBM)||5||Q|
n- Number of objective
T-Thermal cost objective function
W-Wind, cost objective function
S-Solar cost objective function
E-Emissions cost objective function
Con-Number of constraints
ED with cubic cost functions has been studied in the recent past . According Z.X Liang and J.D Glover ,1991,a very crucial issue in DED studies is to determine the order and approximate the coefficients of the polynomial used to model the cost function. This helps in reducing the error between the approximated polynomial along with its coefficients and the actual operating cost. According to Z.X Liang and J.D Glover,1992  and A.Jiang and S.Ertem,1995  to obtain accurate ED results, a third order polynomial is realistic in modelling the operating cost for a non-monotonically increasing cost curve. DED works using cubic cost functions include Bharathkumar S et al, 2013, Hari M.D et al, 2014, Deepak Mishra et al, 2006, and N.A.Amoli et al, 20 12. Krishnamurthy, 2012  used the static cubic function of the emissions dispatch in the Multi Objective Static ED(MOSED) using the Lagrange Method(LM).This provided better results as compared to the quadratic functions. In all these studies, the cubic cost function provided more accurate and practical results as compared to lower order cost functions. A summary of ED works using cubic cost functions is provided in Table 3.From the table, it is clear that cubic renewable and transmission losses cost functions have not been considered. Further, there is need to use more advanced hybrid methods for better results in these complex problems
|Reference||Objectives||Nature of objective functions||Number of Constraints(n)||Method|
|Z.X Liang et al,1991||T||Static||-||Gram-Schmidt, Least Squares|
|Z.X Liang et al,1992 ||T||Static||3||Dynamic programming(DP)|
|A.Jiang and S.Ertem,1995 ||T||Static||2||Newton Method(NM)|
|B.S et al ,2013||T,E||DED with ramp rates and valve points||4||Fuzzy Logic (FL), BF and Nelder-Mead(NM) (FL-BF-NM)|
|Hari Mohan D.et al,2014||T,E||Static||5||PSO-General Search Algorithm (PSO-GSA)|
|Yusuf Somez,2013||T||General static||2||ABC|
|Deepak Mishra et al,2006 ||T||General static||2||OR-Hopefield Neural Network|
|N.A Amoli et al,2012 ||T||Static||2||Firefly Algorithm|
|Krishnamurthy .S et al ,2012||T,E||Static||2||Langrange Method|
|T.Adhinarayananand M.Sydulu,2006||T||Static||2||Lambda-logic based|
|T.Adhinarayananand M.Sydulu,2010 ||T||Static||2||Lambda-logic based|
|E.B Elanchezhian et al,2014||T||Static||8||Teaching learning based optimization (TLBO)|
n- Number of objective
T-Thermal cost objective function
W-Wind, cost objective function
S-Solar cost objective function
E-Emissions cost objective function
Con-Number of constraints
5. MODED Formulation
DED considers change-related costs. The DED takes the ramp rate limits, valve points and prohibited operating zone of the generating units into consideration. The general form of DED was formulated by Yusuf Somez, 2013 as is given by
Where are the cost coefficients of the unit, is the lower generation bound for it unit and is the error associated with the ith equation.When L=3, the cubic form of the SODED results. This can be expressed as
Also the SODED can be written as
where is the output power of the unit at time m, N is the number of generation units and M is the number of hours in the time horizon
To maximize the outcome of the Renewable energy(RE) systems, the implementation of ED combined Thermal and RE(TRE) systems has become of paramount importance. The formulations adopted in this paper are based on the work of Azza A. El Desouky,2013 for the solar and John Hetzer et al ,2008  for the wind cost functions. The operational cost objective function for wind power generation is formulated as
In this case, is the scheduled output of the wind generator in the hour
is the weighted cost function representing the cost based on wind speed profile , is the penalty cost for not using all the available wind power and is the penalty reserve requirement cost which is due to the fact that that actual or available power is less than the scheduled wind power.
Similarly, the operational cost objective function for the PV power generation plant is formulated as 
is the weighted cost function representing cost based on solar irradiance profile is the penalty cost for not using all the available solar power, and is the penalty reserve requirement cost due to the fact that the actual or available power is less than the scheduled PV power.
While attempting to find an optimum DED of generation, transmission loss is one important constraint since the generating centers and the connected load exist in geographically distributed scenario. Since the power stations are usually spread out geographically, transmission network losses must be taken into account to achieve true economic dispatch. Network loss is a function of power injection at each node . Where the real power system transmission losses, PL, is expressed using B- coefficients by the relation
Where is the number of generators , is the number of buses in the system, is the element of the loss coefficient Square matrix, is the element of the loss coefficient matrix and is the constant loss coefficient .
When the loss expression is expanded for a particular number of generators and buses, the equation becomes
This can be expressed as a cubic function
The cost of transmission line losses between plants are accounted with the actual fuel cost function by using a price factor .This factor is defined as the ratio between the fuel cost at its maximum power output to the maximum power output .That is for this multi objective case
Thus, the cost function for the losses at a particular time becomes
Several works on combined economic and emissions dispatch (CEED) have been considered in the past research. Some of these can be found in ,,,,,,,,,  and ;just to list a few. All these works have considered emissions using the quadratic function which is less accurate. Consequently, the more accurate equation for estimating the emissions cost is the third (cubic ) form and it is the function which will be used in the proposed research. It is given by
From equations (2), (4), (5), (10) and (12), the MODED can be formulated as
, are non-negative weights used to make tradeoff (relative importance )between emission security and total fuel cost considering the three fuels such that + =1.N, M and S are the number of thermal wind farm and PV power plants. Is the algebraic sum of the individual weight of the four objectives?
6. Applications of 5-Objective MODED
The MODED formulated can be used to investigate Multi Area Multi Objective Dynamic Economic Dispatch (MAMODED, Dynamic Reactive Power MODED (DRPMODED) and the Security Constrained MODED (SCMODED).These are discussed in this section.
The multi area single objective economic dispatch (MASOED) problem formulation is given by
Where Is the power output of generator n in area m, are the fuel cost coefficients, M is the number of are .
According to Jizhong Zhu, pp. 211-248, (2006) , many approaches have been considered for static MAED. These include NFP with multi area wheeling (NFPMAW), ONN, AHP, DP, Spatial DP(SDP), EP and NLCNFP. Other methods of solving the MAED problems are explained in -. In these methods however, only static MAED is considered. When a MAED problem is solved with spinning reserve constraints and RE cost functions the problem becomes further complicated. The power allocation to each unit is done in such a manner that after supplying the total load, some specified reserve is left for security reasons. This is the dynamic MAED(MAMODED) problem In the event that the power in these areas is changing significantly, the classical methods of MAED solution are no longer applicable in such stochastic environments .Thus more advanced hybrids must be applied to solved the MADED problem which has gained a lot of interest with the integration of renewable energy into the grid .
Thus, the multi area single objective DED (MASODED) with a cubic objective function that is of interest in the proposed research is defined by the relation
The Multi Area Multi Objective DED (MAMODED) which incorporates the Renewable Energy, Transmission line and Emissions Cost functions can be represented by
Subject to Area power balance constraints (APBC)
Generator capacity limit constraint
Tie-line power flow limit
Generation constraint (GRC)
In this formulation, is the MAMODED problem to be solved in the proposed research, , and is the power output of generator n in area mare the fuel cost coefficients for area m and unit n,, is the active load at node n in the area m, is the total real power loss for multi area system, is the number of loads in area m, is the number of tie lines and is the active power flow in the tie line t
The objectives of reactive power (VAR) optimization, which include RPED, are to improve the voltage proﬁle, to minimize system active power losses, and to determine optimal VAR compensation placement under various operating conditions. To achieve these objectives, power system operators utilize control options such as adjusting generator excitation, transformer tap changing, shunt capacitors, and SVC. However, the size of power systems and prevailing constraints produce strenuous circumstances for system operators to correct voltage problems at any given time. In such cases, there is certainly a need for decision - making tools in predominantly ﬂuctuating and uncertain computational environments. There has been a growing interest in VAR optimization problems over the last decade. Methods applied so far to solve the SORPED are as found in -.Solving ORPED is gaining more importance due to their effectiveness in handling the inequality constraints and discrete values compared to that of conventional gradient-based methods. EAs generally perform unconstrained searches, and they require some additional mechanism to handle constraints. In the literature, various constraint handling techniques have been proposed. However, to solve ORPD the penalty function approach has been commonly used, while the other constraint handling methods remain untested. However, the excessive time consumption of EP and GAs will limit their applications in power systems, especially during real - time operation. Thus better method are needed to handle the more complex problems where stochastic reactive power from wind and solar generators are involved.
The fuel cost in terms of reactive power output can be expressed as
where and are the reactive power cost coefficients calculated using a curve fitting method is the reactive power generated by generator and is the order of the fuel cost function
Power balance constraints
Continuous control variable (Generator Bus Voltage)
Discrete control variable (Transformer Tap Settings)
where is the tap setting of transformer at branch k
Reactive power balance
In this formulation is the reactive power generated by the ith capacitor bank, is the reactive power generated at bus i, is the apparent power flow through the ith branch, is the total number of buses,the number of tap setting transformer branches, is the number of capacitor banks and is the number of generator buses, Reactive power generated by generator i, Reactive power generated and absorbed by VAR compensation device j such as capacitors,SVC,Wind Based DFIGs, and PV generators,Reactive power load at load bus k and Power system reactive power power loss and absorption .Further , Where is the voltage magnitude at bus i , is the voltage magnitude at bus j, is the real and reactive powers injected at bus i, is the mutual conductance and suspectance between bus I and j , is the total number of buses excluding the slack bus, is the number of PQ buses and is the voltage angle difference between bus I and bus j
The security - constrained economic dispatch (SCED) is one of the simpliﬁed optimal power ﬂow (OPF) problems. It is widely used in power industry. The main objective of electric power dispatch is to provide electricity to the customers at low cost and high reliability. Transmission line failures constitute a great threat to the electric power system security.
The static single objective SCED (SSOSCED) problem is formulated as
The SOSCDED in the proposed research uses the more accurate cubic function and incorporates the error ,valve points and ram rate terms. The multi objective SCED(MOSCDED) which investigates the effects of renewable energy ,power loss and emissions functions to the energy security of the power system is formulated as
7. MODED Optimization Methods
The methods that have been to solve the DED Problem so far can be classified into three groups. These are, Deterministic, Heuristic and Hybrid methods. These methods have been discussed in the following sub sections.
7.1. Deterministic (D)
Unconstrained Methods (UCM) convert constrained problems into unconstrained form. This forms the basics of the formulation of the constrained optimization Algorithms. These methods include Gradient Search (GS),Line Search(LS),Lagrange Multiplier(LM),Newton Raphson(NR) ,Trust Region (TR) ,Quasi-Newton (QN),Double-dogleg(DD) and Conjugate Gradient(CG) Methods. In DED Only LM method has been used by W.G. Wood, 1982  and W.R. Barcelo and P. Rastgoufard,1997  and Maclaurin Series-Based Lagrangian Method (MSBLM) was later applied to the DED problem by Hemamalini et al, 2010 .
Linear Programming (LP) linearizes the nonlinear power system DED problem so that the objective function and the constraints have the linear form. It was Ahmed Farag et al, 1995  and Y.H. Song and I. Yu, 1997  and it proved reliable in terms of convergence .The method was found Quick in identifying infeasibility and thus able to accommodate large variety of power system operating limits and contingency constraints. However the method was found inaccurate in evaluating the power system losses and therefore had insufficient ability to find an exact solution as compared with more accurate Non Linear models. However the DED solutions generally met the requirements of Engineering precision when applied to the SCED and RPED optimization .
Non Linear Programming (NLP) is a first order method .Since DED problems are nonlinear, the method able to handle the nonlinear objectives and constraints .A search direction is chosen and the search is done using reduced gradient method. As applied by P.P.J. van den Bosch, 1985  and an improved version of NLP called Non Linear Constrained Method (NLCM) developed by Faisal .A. Mohamed, 2009 ,the method has several advantages as compared to the LP Method .These include higher accuracy and global convergence ,that is, convergence is guaranteed independent of the starting point. However the method has a slow rate of convergence due to a zigzagging problem in the search direction. The method has also been used in SCED, MAED and RPED problems .
Other deterministic methods that have been used to handle DED problems include Economic Load Allocation and Supplementary Control Action (ELASCA), Bechert and Kwatny, 1972  ,Muller Algorithm (MA), Chandram et al, 2006 ,Brent Method (BM), K. Chandram et al, 2008  and Model Predictive Control (MPC), X. Xia et al, 1999 , X. Xia et al,2011  just to mention a few.
7.2. Non Quantity Approaches (N)
These methods are applied to problems with uncertainties in their parameter variations. Uncertainty due to insufficient information generate an uncertain region of decisions and therefore results from average values cannot be used to represent the uncertainty level. NQA methods are meant to account for the uncertainties in information and goals related to multiple and usually conflicting objectives in power system optimization .These methods find application in ED,HTS,UC,RPED and State Estimation problems with uncertainties. Compared to the deterministic methods, the NQA Methods are effective in handling uncertainties, that is, they compute the unavailable or uncertain data so that the MODED problem can be solved even when some data for the wind and solar objective cost functions are not available. These methods include Optimization Neural Networks (ONN), Artificial Neural Network(ANN), Probabilistic /Characterization Theory (PT),Fuzzy Set Theory (FST),Analytical Hierarchical Process (AHP),Risk Management Tools(RMT) and Cost Benefit Analysis(CBA).These Methods have been used in the formulation of hybrids used in the solution of DED problems with uncertainties. Examples of such include,ANN ,Heuristic Neural Network(HNN)[83-84],Trust Region (TR) and Goal Programming (GP)  , Chaotic Fuzzy(CF) and Variable Step Size (VSS) .
7.3. Heuristic (H)
Evolutionary Algorithms (EAs) are based on natural evolution .They are population based optimization processes. There is no need to differentiate the cost functions and the constraints in the ED formulation. The method can be categorized as GA, EP and DE. The three however use the same mutation, recombination, reproduction, cross over and selection operators . DED using Genetic Algorithms (GA) and the derived methods include Elitist Genetic Algorithm (EGA),Fly et al, 1997 , Parallel Micro GA (PMGA), Ongsakul et al ,2002  and Elitist non-dominated sorting genetic algorithm-II (ENSGA-II), Robert.T.F et al,2004.Evolutionary Programming (EP) include Constrained Evolution Programming (CEP),Shailti et al, 2005 . EP, Joned et al, 2006  and Quantum Evolutionary Algorithm (QEA), Babu et al, 2008 . Differential Evolution (DE) is the most recent EA and the DED works done include Modified DE (MDE),Yuang et al ,2008 , K.Deb,2000 ,Variable Scaling Hybrid Differential Evolution (VSHDE), J.Chiou,2009  and Cellular Differential Evolution (CDE) Noman and Iba, 2011 .EAs methods have been used to solve simplified CED, DED SCED RPED and MAED problems. However these methods require all information to be included in the fitness function therefore it is difficult to consider all the MODED problem objectives and constraints in the objective function.
Particle Swarm Optimization (PSO) includes all the swarm intelligence algorithms inspired by the social dynamics and an emergent behavior that arises in socially organized colonies. These methods have been used in all hybrid algorithms where an accurate search is needed in DED problems. They include Particle Swarm Optimization (PSO) Z.L. Gaing, 2004 , Hardiansyah et al 2012, Zhao et al, 2004  and Modified Adaptive PSO (MAPSO), Niknam et al, 2011.More specific PSO methods used in DED problems are BAT Algorithm, Anti Colony Optimization(ACO),Bee Colony Optimization(BCO) ,Artificial Bee Colony (ABC),Bacterial Foraging (BF),BAT Algorithm, Teacher Learning Algorithm(TLA) , Bandi Ramesh et al, 2013,Biogeography-Based Optimization (BBO), Divya Mathur, 2013, .
Other Methods that have been used in DED include Simulated Annealing (SA), Panigrahi et al, 2006 ,Pattern Search (PS),Alsumait et al, 2010 ,Adaptive Look-Ahead (ALA), Han et al, 2001 ,Feasibility-Based Selection Comparison (FBSC),K.Deb,2000 ,Ongsakul et al, 2002 , X.Yuan et al,2009,Artificial Immune System (AIS), M.Basu, 2011 Hemmalini et al,2011 ,Modified Teaching-Learning Algorithm (MTLA) Niknam, T et al ,2013  and Optimality Condition Decomposition (OCD), Rabiee, A et al , 2014 
The reason for increased use of hybrids is because they exalt the strengths and improve the weaknesses of the methods concerned. Deterministic,Heuristic and Non-Quantity Approaches are used in hybrid formation .There are two types of hybrids depending on the number and Type of methods used.
7.5. Two Method Hybrids
There are eight possible types of two-method hybrids. These include N-N, N-D, and N-H, D-D, D-H, D-N, H-D, H-H and H-H approaches.
Deterministic-Deterministic (D-D): These method involve two deterministic methods .They include Quadratic Programming (QP) and Linear Programming (LP) – (QP-LP), Somuah and Khunaizi, 1990 , QP-LP, Han et al, 2001  and Interior Point (IP) and Quadratic Programming (QP) - (IP-QP), Lin and Chen, 2002
Deterministic-Heuristic (D-H): These involve a deterministic method with a Heuristic method to fine tune it. These include SA and GA – (SA-GA), Ongsakul et al, 2001 , SQP-EP, Deterministically Guided (DG) PSO (DG-PSO), Victoire et al, 2005 [113-115] and SA and Direct Search Method (DSM)-(SA-DSM), C .Chen, 2007 
Deterministic- Non Quantity Approaches (D-N) or Non Quantity Approaches-Deterministic (N-D): This involves a base deterministic method with a non-quantity method to deal with the uncertainities.The vice versa is also true. Examples in this context include Lambda-Iteration Technique (LIT) and Artificial Neural Networks (ANN)-(LIT-ANN), R.H Liang, 1999  and Heuristic Neural Networks (HNN) and QP (HNN-QP), Abdul-Aziz et al, 2008 [83-84].
Heuristic -Deterministic (H-D): In this case a hybrid method is used as the base algorithm then its weaknesses are strengthened using a deterministic method. In DED ,examples include Genetic Algorithm (GA) and Gradient Search (GS)-(GA-GS), F.Li et al, 1997 ,Relaxed GA (RGA) and Gradient Search (GS)-(RGA-GS) ,Li and Aggarwal,2000 ,EP and Sequential QP (EP-SQP), P.A et al, 2002 ,PSO-SQP, Victoire et al,2005 [113-115],Swarm Direction Technique (SDT) and Fast Evolutionary Programming (FEP) -(SDT-FEP), Zwe-Lee et al,2009 ,Evolution (DE) and Local Random Search (LRS)-(DE-LRS), Y. Lu et al, 2009 ,Improved Differential Evolutionary (IDE) and Shor’s Algorithm (SHA)-(IDE-SHA), Yuan et al,2009 , F.Kappel and A.V Kuntsevich,2000 ,PSO and Harmony Search (HS)-(PSO-HS), Ravikumar et al,2011  ,Seeker Optimization Algorithm (SOA) and SQP- (SOA-SQP), S. Sivasurbramani et al,2010 ,Hybrid Differential Evolution (HDE) and Sequential Quadratic Programming (SQP)-(HDE-SQP), A. M. Elaiw,2012  ,Real Coded GA (RCGA) with Quasi-Simplex (QS)-(RCGA-QS), Zhang et al, 2006 
Heuristic- Heuristic (H-H): In this type two heuristic methods are involved, on serves as the base and the other is used to improve the weaknesses of the first for accurate results. Examples include PSO and DE (PSO-DE), Dun –Wei Gong et al, 2010,Firefly Algorithm (FA) with GA Mutation (GAM)-(FA-GAM) , Niknam, T,2012 ,Multi Objective PSO (MOPSO) and Brent method(BM)-(MOPSO-BM), H. Shayeghi and A. Ghasemi ,2012,Modified Artificial Bee Colony (MABC) with Differential Evolution (DE)-(MABC-DE), Hardiansyah ,2013 and PSO and Simulated Annealing (SA) –(PSO-SA),V.Karthikeyan et al,2013
Non Quantity Approaches-Heuristic (N-H) 0r Heuristic- Non Quantity Approaches (H-N): In this case a non –quantity method is used to address an uncertainty before a heuristic method is applied. The vice versa is also true Examples include Fuzzy Optimization Technique (FOT) and Goal Satisfaction Concept (GSC)-(FOT-GSC) P.A et al, 2004  and Chaotic Differential Evolution (CDE), Y.Lu et al,2011 .
Non Quantity Approaches-Non Quantity Approaches (N-N): Since non quantity methods are only used to address uncertainities,then a combination of two or more of such methods in the formulation of an hybrid is not feasible.
7.6. Three Method Hybrids
Recent trends in Hybrid formation involve the use of three methods. There are twenty seven possible three-method hybrids that can be used. In this section we discuss such methods that have been used in DED solution. These are summarized in Table 4
Heuristic algorithms seem to have shared the same dominance as deterministic algorithms. This is because heuristic algorithms, unlike deterministic ones, are derivative-free, and capable of solving optimization problems without requiring convexity.They are also independent of the initial solution, and have the ability to avoid being trapped in local optima.
On the other hand, heuristic algorithms have drawbacks such as being problem dependent, requiring parameter tuning, and unable to guarantee global solution attainment. A combination of a both methods into a hybrid has proved to solve this problem. This explains why publishing research work using hybrid methods has been becoming increasingly popular. This is as shown in Figure 3.
Also in the last decade integration of RE cost function has been of increasing interest. This is because RE is well utilized once in the national grid. Therefore, better methods for handling DED with wind and solar cost functions need to be developed.
|Type||Examples in DED|
|H-H-D||EP, PSO, and SQP (EP-PSO-SQP), S. Titus and A.E Jeyakumar,2008 |
|H-H-H||BF-PSO-DE, Praveena et al ,2010 |
|H-D-H||GA and Local Search(LS) -Goal Programming(GP) -(GA-LS-GP), A.A.Mousa et al,2011 |
|D-H-H||Mixed Integer Quadratic Programming (MIQP), the Warm Start Technique (WST) and the Range Restriction Scheme (RRS)-(MIQP-WST-RRS), Wang et al, 2014 |
|N-N-H||Variable Step Size(VSS) Chaotic Fuzzy(CF) Quantum Genetic Algorithm, (VSS-CF-QGA), Wenxia Liu, et al,2014 Trust Region (TR),Goal Programming(GP) and PSO (TR-GP-PSO), Amhed et al,2013|
This paper has psesented a detailed review of the MODED problem in terms of the formulation and the methods applied in the solution.It is apparent that cubic cost functions are found to be more accurate than quadratic ones. Further, fully constrained 5-Objective MODED with cubic functions to give more realistic results. However fourth order functions are predicted to be even better than the third order systems.
Dynamic trade in power is supposed be cheaper than the static one as the Improved Wheel Method (IWM) will be more realistic in modelling the areas and the tie lines. Dynamic reactive power (two-way) from the renewable sources will provide a far much improved voltage profile as compared to the static reactive power. In addition, there is foreseen better security due to the full constraining (with, uncertainty) and two way reactive power.
Therefore ,there will be a more accurate MODED modelling ,new and better method for MODED Solution, Improved power system Security ,dynamic MAED and a new way for Reactive Power Dispatch.
A more recent trend for solving MODED is the three-method hybrids formulation in which all the weaknesses of the base methods are suppressed and the strengths exalted. This leads to increased accuracy and speed in handling higher order cost functions with more objectives.