Open access peer-reviewed chapter

Solution of Combined Economic Emission Dispatch Problem with Valve-Point Effect Using Hybrid NSGA II-MOPSO

Written By

Arunachalam Sundaram

Submitted: 27 May 2017 Reviewed: 29 November 2017 Published: 20 December 2017

DOI: 10.5772/intechopen.72807

From the Edited Volume

Particle Swarm Optimization with Applications

Edited by Pakize Erdoğmuş

Chapter metrics overview

1,463 Chapter Downloads

View Full Metrics

Abstract

This chapter formulates a multi-objective optimization problem to simultaneously minimize the objectives of fuel cost and emissions from the power plants to meet the power demand subject to linear and nonlinear system constraints. These conflicting objectives are formulated as a combined economic emission dispatch (CEED) problem. Various meta-heuristic optimization algorithms have been developed and successfully implemented to solve this complex, highly nonlinear, non-convex problem. To overcome the shortcomings of the evolutionary multi-objective algorithms like slow convergence to Pareto-optimal front, premature convergence, local trapping, it is very natural to think of integrating various algorithms to overcome the shortcomings. This chapter proposes a hybrid evolutionary multi-objective optimization framework using Non-Dominated Sorting Genetic Algorithm II and Multi-Objective Particle Swarm Optimization to solve the CEED problem. The hybrid method along with the proposed constraint handling mechanism is able to balance the exploration and exploitation tasks. This hybrid method is tested on IEEE 30 bus system with quadratic cost function considering transmission loss and valve point effect. The Pareto front obtained using hybrid approach demonstrates that the approach converges to the true Pareto front, finds the diverse set of solutions along the Pareto front and confirms its potential to solve the CEED problem.

Keywords

  • multi-objective optimization
  • economic emission dispatch
  • Pareto optimality
  • NSGAII
  • MOPSO
  • B-loss coefficients

1. Introduction

In order to operate the power system economically and also to protect the environment from pollution the power system operator has to carry out optimal scheduling of active power to simultaneously minimize the fuel cost and the emissions from the fossil fuel-fired power plants. These objectives are desirable to obtain great economic benefit [1] and to reduce the nitrogen oxide (NOx), sulfur oxide (SOx) and carbon dioxide (CO2) pollutants which cause harmful effect on human beings [2]. These conflicting objectives can be formulated as a multi-objective combined economic emission dispatch (CEED) problem. This CEED problem can be solved using traditional mathematical programming techniques such as lambda iteration, gradient search [1] and can also be solved using modern heuristics optimization techniques. The numerous advantages of solving the CEED problem using heuristic optimization methods compared to the traditional mathematical programming techniques are they are population-based, do not require any derivative information, do not use gradient information in search process, use stochastic operators in search process, they are simple to implement and flexible, have inbuilt parallel architecture and they are scalable and are also computationally quick.

A single optimal solution cannot be obtained for a multi-objective CEED problem which simultaneously minimizes the conflicting objectives of fuel cost and emission. Thus the simultaneous minimization of conflicting objectives in a multi-objective optimization problem (MOP) gives rise to a set of tradeoff solution called as Pareto-optimal (PO) solutions [3] which needs further processing to arrive at a single preferred solution. In literature domination based framework using multi-objective evolutionary algorithms (MOEA) which simultaneously minimizes the fuel cost and emission have been employed to solve the CEED problem. These population-based approaches can obtain the multiple non dominated solutions in a single simulation run. These non-dominated solutions portray the tradeoff between fuel cost and emission objectives of CEED problem. Modern meta-heuristic optimization algorithms like Genetic Algorithm [4, 5], Biogeography Based Optimization [6], Particle Swarm Optimization [7], Bacterial Foraging Algorithm [8], Scatter Search [9], Teaching Learning Based Optimization [10], Differential Evolution [11] and Harmony Search Algorithm [12] have been developed and successfully implemented to solve this complex, highly nonlinear, non-convex CEED problem.

The multiple objective CEED problem can also be transformed into a single objective problem using a weighted sum approach and h parameter values. The h parameters are used to overcome the dimensionality problem when combining multi-objectives and the converted single objective problem is then solved using evolutionary algorithms [13, 14, 15]. Another technique to solve CEED problem without the h parameter is to normalize the fuel cost and emission components [6] and solve the single objective function using evolutionary algorithms (EA). In these approaches for the chosen value of weights will give one particular PO solution at a time. However, the disadvantage of these methods is that it requires multiple runs to find the set of PO solutions.

Each evolutionary algorithm has its own characteristics and merits; therefore it is natural to think of integrating these different algorithms to handle a complex problem like CEED. In the research field of Evolutionary Algorithms merging of two or more optimization algorithms into a single framework is called hybridization. In [16, 17, 18, 19, 20, 21] hybrid multi-objective optimization algorithms have been successfully applied to solve CEED, various complex engineering problems, and standard test functions. The results indicate that the hybrid algorithms are effective, can exchange elite knowledge within the hybrid framework, can do parallel processing, can improve the exploration and exploitation capabilities and can yield more favorable performance than any single algorithm.

In order to obtain a globally optimal solution without being trapped in local optima requires a tradeoff between exploration and exploitation task in the search process. Exploration phase in any algorithm is important to search every part of the solution domain to provide an estimate of the global optimal solution. On the other hand exploitation phase in any algorithm is important to improve the best solutions found so far by searching in their neighborhood. In this chapter, a hybrid framework using Non-Dominated Sorting Genetic Algorithm II (NSGA II) [22] and Multi-objective Particle Swarm Optimization (MOPSO) [23] is used to solve the CEED problem. This hybrid framework integrates the desirable features of the NSGA II and MOPSO while curbing their individual flaws. These population-based approaches use different techniques for exploring the search space and when they are combined will improve the tradeoff between the exploration and exploitation tasks to converge around the best possible solutions. The main purpose of this hybridization technique is to obtain a well-spread and well-diverse PO solution. When the proposed hybrid algorithm is used to solve the highly complex CEED problem the PO solution is obtained in less number of iteration and is also computationally fast when compared to MOPSO.

The rest of the chapter is organized as follows. The next section formulates the CEED problem. In Section 3, the transmission loss handling procedure and the constraint handling procedure is explained. In Section 4 the short review of NSGA II and MOPSO is provided. Section 5 is devoted to explaining the hybrid algorithm. In Section 6 the hybrid algorithm is applied on standard IEEE 30 bus systems and it also discusses the simulation results. Finally, the conclusion is drawn in Section 7.

Advertisement

2. Formulation of combined economic emission dispatch (CEED) problem

The combined economic emission dispatch problem has two conflicting objectives. The first objective can be stated as determining the optimal power generation schedule from a set of online generating units to satisfy the load demand subject to several physical and operational constraints to minimize the fuel cost. The second objective can be stated as determining the optimal power generation schedule from a set of online generating units to satisfy the load demand to minimize the pollutant emissions produced by the generating units. Both the conflicting objectives have to be minimized at the same time because operating the system with minimum cost will result in higher emission and considering only the minimum environmental impact is not practical which results in high production cost of the system. This section formulates the objective functions of the CEED problem along with equality and inequality constraints to maintain rigorous standards to meet the practical requirements of the power system. The goal of this chapter is to find the Pareto-optimal solutions of the CEED problem which minimize both these objectives subject to constraints. The mathematical formulation is as follows.

2.1. Objective functions of CEED problem

The general formulation for a multi-objective optimization problem (MOOP) is to minimize the number of objective functions simultaneously. A general mathematical model is represented as follows [21]:

Minimizefx=f1xf2xfmx,xDE1

where fx represents the vector of objectives and fix,i=1,2,,m is a scalar decision variable which maps decision variable x into objective space fi=RnR. The n-dimensional variable x is restricted to lie in a feasible region D which is constrained by jin-equality constraint and kequality constraint, i.e.

D=x:gjx0hkx=0j=12Jk=12KE2

The decision variable x can be written more suitably as

x=x1x2x3xnTE3

where T is the transposition of the column vector to the row vector. The decision variables are restricted to take a value within a lower ximin and upper ximax bounds. These bounds are called the decision space [3].

In MO CEED problem the number of objectives m=2.The mathematical model of CEED is represented as follows:

Minimizefx=f1xf2x,xDE4

subject to power balance equality constraints hxand bounds. The function f1x represents the minimization of total fuel cost function and the function f2x represents the minimization of the emissions from the fossil fuel fired plants. The decision variable x consists of the real power generation of the ngenerating units and can be written as

x=Pg1Pg2Pg3PgnTE5

where Pgi is the real power output of the ith generator.

Power plants commonly have multiple valves that are used to control the power output of the units. In a practical generating unit, when steam admission valves in thermal units are first opened, a sudden increase in losses is registered which results in ripples in the cost function. In order to model these ripples accurately, sinusoidal functions are added to the quadratic cost function [24]. The resulting cost function contains higher order nonlinearity and makes the problem non-differentiable and non-convex. Hence there are two versions of the fuel cost function, the quadratic function represented by f1xand the combination of quadratic and a sinusoidal (valve-point) function represented byf1,Vx . The two versions of the fuel cost functions are given below

f1Pg=i=1naiPgi2+biPgi+ciE6
f1,VPg=i=1naiPgi2+biPgi+ci+eisinfiPgiminPgiE7

where ai,bi,ci represent the cost coefficients of the generator i. ei and fi are coefficients to model the effect of valve point of the generator i.

The second objective function f2xis an emission function which takes into account the major pollutants caused by the fossil fuel fired power plants. The main pollutants from the power plants are the sulfur oxides and nitrogen oxides. The sulfur oxide emissions are proportional to the fuel consumed by the power plants and have the same form as that of the fuel cost function given by (6). The sulfur oxide emission function can be stated as follows [7].

f2,soPg=i=1nSi,1+Si,2Pgi+Si,3Pgi2E8

The nitrogen oxides emissions are difficult to evaluate as the nitrogen is available in air and also in the fuel. The production of nitrogen gas is related to boiler temperature and air content. The modeling of the nitrogen oxides consists of straight lines and exponential terms. The nitrogen oxides emission function can be stated as follows

f2,NoPg=i=1nNi,1+Ni,2Pgi+Ni,3eNi,4PgiE9

The total emission function is obtained by adding the coefficients of (8) and (9) which gives the combination of the mixture of sulfur oxides and nitrogen oxides pollutants [7]. The total emission function can be stated as follows

f2Pg=i=1n102αi+βiPgi+γiPgi2+ηieδiPgiE10

The total emission function given by (10) has a quadratic term and an exponential term which makes the function highly nonlinear. In (10) αi,βi,γi,ηi,δi are the emission coefficients of the generatori. The modeling of the emission function is very important because according to the Amendments of the Clean Air Act regulatory agencies might decide to limit power plant emission in the areas where there are high concentrations of harmful contaminants.

2.2. Active power balance equality constraint and bounds

In order ensure that the total real power generation exactly match with the total load demandPd and transmission lossPl in the system a power balance equality constraint given in (11) should be satisfied.

hx=i=1nPgiPdPl=0E11

The transmission losses in the power network are function of Pg and can be represented using B-matrix coefficients (Kron’s loss formula [1]) as follows

PlPg=i=1nj=1nPgiBijPgj+i=1nB0iPgi+B00E12

where Bij,B0i,B00 are transmission loss coefficients. There are instances in literature where the power losses in the system is neglected and the power balance equation given by (11) is curtailed as follows

hx=i=1nPgiPd=0E13

The above equations given by (11) and (13) are most common form of power balance equation found in the literature.

The power output of each generator ishould lie within its minimum limit (Pgimin) and maximum limit (Pgimax) given by

PgiminPgiPgimax;i=1,2,3,nE14

2.3. Combined economic emission dispatch

The purpose of the CEED problem is to determine the Pareto-optimal real power generation vector x=Pg1Pg2Pg3PgnT that minimize the two conflicting objective given by (7) and (10) while satisfying the real power equality constraint given by (11) and the bounds given by (14). The bi-objective CEED problem can be formulated as

Minimizef=f1,VPgf2PgE15

In MO CEED problem, the economic and emission objectives will conflict with each other and is not possible to satisfy them simultaneously. There is no way of improving these objectives without degrading at least one of these objectives and the resulting set of non-dominated solutions thus obtained are called Pareto-optimal set. The objective function values of all elements in the PO set in the objective space constitute the Pareto front. When the sufficient number of PO solutions is available for the CEED problem then it is possible to find a convex curve containing these solutions to produce the Pareto front. The two main goals of MO CEED problem:

  1. Find a set of non-dominated solutions which lie on the Pareto-optimal front

  2. Find a wide spread of non-dominated solutions to represent the entire range of the Pareto-optimal front.

Advertisement

3. Constraint handling mechanism

At any stage of the algorithm whenever a new population is being generated it is very important to make sure that the population lies within the decision space. While solving the CEED problem this implies that the population should satisfy the equality constraints and bounds. If the transmission losses are neglected than the kth variable of the candidate solution Pgk can be calculated by subtracting the sum of the power generations (excluding the kthvariable) i=1n1Pgifrom the power demandPd. If the power transmission losses are considered, to determine Pgk and to maintain the equality constraint becomes hard. It is done using the following steps.

  1. Step 1. Update the variables belonging to the set αn by normal optimization process of an evolutionary algorithm.

    Pgi=PgiminrandPgiminPgimax;iαnE16

    Here randis a uniformly distributed random number in the range of01. The set αn contains all the integers in the range 1n exceptk, where kis a randomly generated integer which lies in the range of1n

  2. Step 2. If updating of the variables is carried out using any other technique then regulate the updated variables which violate the lower bounds asPgi=Pgimin;iαn. Regulate the updated variables which violate the upper bounds asPgi=Pgimax;iαn.

  3. Step 3. Obtain the value of the kth variable of the candidate solution Pgk by solving the following quadratic equation (17) whose coefficients are associated with the variables belonging to the set αn and the transmission loss coefficients [7]. To improve the potential candidate solution and also to improve the flexibility and diversity of the optimization algorithm the value of k is randomly generated integer between 1 andn.

    BkkPgk2+2iαnBkiPgi+B0k1Pgk+Pd+iαnjαnPgiBijPgj+iαnB0iPgiiαnPgi+B00=0E17

    Out of the two roots of the quadratic equation (17), one root will be selected as the value of the variable Pgkusing the following procedure. If both the roots of the quadratic equation lie within the bounds then the root which has the minimum value is selected. If only one root lies within the bounds, this root is selected as the value of Pgk and the other root which lies outside the bounds is neglected. If both the roots lay outside the bounds the value of Pgkis set equal toPgkmin.

  4. Step 4. Calculate the residue PRD by subtracting the total system demand Pdand the total system transmission loss Plfrom the sum of the total power generationi=1nPgi. IfPRD<tol, then go to step 7; otherwise go to step 5. Here, tolis the demand tolerance usually set as 0.001p.u.

  5. Step 5. Recalculate Pgiusing Eq. (16).

  6. Step 6. Repeat step 3, step 4 and step 5 untilPRD<tol. This step will ensure that the candidate solution will always lie within the decision space.

  7. Step 7. Stop the constraint handling procedure.

The main purpose of this constraint handling mechanism is to increase the flexibility and diversity of the algorithm and to make sure that the candidate solution generated at any point of the algorithm always lies within the decision space.

Advertisement

4. NSGA II and MOPSO algorithms for solving CEED problem

Several Evolutionary Multi-objective (EMO) algorithms like NSGAII, MOPSO, SPEA 2 (Strength Pareto Evolutionary Algorithm), GDE 3 (Generalized Differential Equation) have been designed and used in solving numerous complex real word problems involving two or more objectives. All these algorithms can find the multiple Pareto-optimal solutions in a single run. Out of all these available algorithms, two of the widely used reliable methods for solving bi-objective optimization problems are the NSGA II and MOPSO. This section provides the review of these two EMO algorithms.

NSGA II was proposed in [22] as an improvement of the NSGA proposed in [25]. This NSGA II algorithm was the revised version of NSGA to overcome the following criticisms:

  • Computational complexity associated with non-dominated sorting.

  • Lack of elite-preserving strategy.

  • Lack of maintaining diversity among obtained solutions.

The NSGA II algorithm is very efficient for solving multi-objective optimization problems since it incorporates an efficient elitism preserving technique using non-domination sorting. The population is ranked based on non-domination sorting before the selection is performed. All non-dominated individuals are classified into one category. Another layer of non-dominated individuals are considered after the group of classified individuals are ignored. This process is continued until all individuals in the population are classified. NSGA II also uses a mechanism for preserving the diversity and spread of the solutions without specifying any additional parameters (NSGA uses fitness sharing). This crowding distance operator guides the selection process towards a uniformly spread out Pareto-optimal front. The NSGA II algorithm for solving the CEED problem is stated below:

  • Specify the parameters for the CEED problem

    • The total demand of the power system Pd

    • Fuel cost and emission coefficients for each generating unit

    • B matrix coefficients for transmission loss calculations

    • Number of decision variables nVar

    • Lower bounds of the decision variables VarMin

    • Upper bounds of the decision variables VarMax

  • Specify the parameters for NSGA II Algorithm

    • Population Size nPop

    • Maximum number of iteration MaxIt

    • Crossover Percentage pCrossover

    • Mutation Percentage pMutation

    • Mutation rate mu

    • Mutation step size sigma

  • Initialize Population

    • Generate a random nPop size population

    • Once the random population is initialized the Constraint Handling Mechanism proposed in Section 3 is carried out.

  • Evaluate the objective functions

    • Evaluate the fuel cost objective function E and emission objective function F

  • Perform Non Domination Sorting

  • Calculate Crowding Distance and rank the population based on Non Dominated fronts

  • For each generation do

    • Create offspring population

      • Selection, Crossover and Mutation

    • Apply Constraint Handling Mechanism

    • Evaluate the fuel cost objective function E and emission objective function F

    • Merge the parent and offspring population

    • Perform non domination sorting

    • Calculate crowding distance and rank based on non-domination fronts

    • Select solutions

      • Each front is filled in ascending order

      • Last front-descending order of crowding distance

    • Store the non-dominated solutions in list Ϝ1

    • Plot the non-dominated solutions in list Ϝ1

    • Increment generation count

  • End for

In order to handle multiple objectives Pareto dominance is incorporated into PSO algorithm and the MOPSO algorithm is proposed in [23]. The algorithm proposed in [23] uses an external repository of particles to keep a record of the non-dominated vectors found along the search process. At each generation, for each particle in the swarm, by using Roulette wheel selection, a leader is selected from the external repository. This leader then guides other particles towards better regions of the search space by modifying the flight of the particles. A special mutation operator is applied to the particles of the swarm and also to the range of each design variable of the problem to be solved to improve the explorative behavior of the algorithm. The value of the mutation operator is decreased during the iteration. To produce well spread Pareto fronts the MOPSO algorithm in [23] uses an adaptive grid. The MOPSO algorithm for solving the CEED problem is stated below:

  • Specify the parameters for the CEED problem

    • The total demand of the power system Pd

    • Fuel cost and emission coefficients for each generating unit

    • B matrix coefficients for transmission loss calculations

    • Number of decision variables nVar

    • Lower bounds of the decision variables VarMin

    • Upper bounds of the decision variables VarMax

  • Specify the parameters for MOPSO Algorithm

    • Maximum number of iteration MaxIt

    • Population Size nPop

    • Repository size nRep

    • Inertia Weight w and Inertia Weight damping rate wdamp

    • Personal learning coefficient c1 and Global learning coefficient c2

    • Number of grids per dimension nGrid

    • Inflation Rate alpha, leader selection pressure beta, Deletion selection pressure gamma

    • Mutation rate mu

  • Initialize Swarm Population

    • Generate a random swarm particles

    • Once the random particles are initialized the Constraint Handling Mechanism (Section 3) is carried out.

    • Store the values of the particles as their personal best pBest

  • Determine Domination

    • Initialize external repository rep

    • Create grid and find grid index

  • For each generation do

  • For each particle do

    • Select leader from external repository

    • Update particle position and velocity

    • Apply Constraint Handling Mechanism

    • Evaluate the fuel cost objective function E and emission objective function F

    • Apply Mutation and calculate new solutions

    • Apply Constraint Handling Mechanism

    • Determine Domination

    • Update pBest

  • End for

    • Add non dominated particles to the repository

    • Determine domination of new repository members

    • Keep only the non-dominated members in the repository

    • Update grid and grid index

    • If repository is full delete members

    • Plot the members in the external repository

    • Modify inertia weight

  • End for

Advertisement

5. Hybrid NSGA II and MOPSO algorithm for solving CEED problem

The mechanism of the proposed hybrid approach for solving the CEED problem is to integrate the desirable features of NSGA II (retaining the elitism feature) and MOPSO (exploitation capability) while curbing the individual flaws (NSGAII––does not have an efficient feedback mechanism, PSO overutilization of resources). The mechanism to explore the search space differs in both the algorithms. GA uses mutation and crossover operators which will enhance the exploration task of the hybrid algorithm. The particles in PSO are influenced by their own knowledge and information shared among swarm members. PSO enhances the exploitation task of the hybrid algorithm by finding better solutions from the good ones by searching the neighborhood of good solutions. In this hybrid algorithm at every generation, the Pareto dominance of the population is computed and based on these values non dominated sorting is performed [19]. In order to avoid premature convergence, the elite upper half of the population are enhanced by NSGA II algorithm while the lower half of the population are considered as swarm particles and are optimized by MOPSO to make them converge around the best possible solutions. The hybrid NSGA II-MOPSO algorithm for solving the CEED problem is stated below:

  • Specify the parameters for the CEED problem

  • Specify the parameters for NSGA II Algorithm

    • Population Size nPop

    • Maximum number of iteration MaxIt

    • Crossover Percentage pCrossover

    • Mutation Percentage pMutation

    • Mutation rate mu

    • Mutation step size sigma

  • Specify the parameters for MOPSO Algorithm

    • Repository size nRep

    • Inertia Weight w and Inertia Weight damping rate wdamp

    • Personal learning coefficient c1 and Global learning coefficient c2

    • Number of grids per dimension nGrid

    • Inflation Rate alpha, leader selection pressure beta, Deletion selection pressure gamma

    • Mutation rate mu

  • Initialize Population

    • Generate a random nPop size population

    • Once the random population is initialized the Constraint Handling Mechanism proposed in Section 3 is carried out.

  • Evaluate the objective functions

    • Evaluate the fuel cost objective function E and emission objective function F

  • For each generation do

    • Perform Non Domination Sorting

    • Calculate Crowding Distance and rank the population based on Non Dominated fronts

    • Truncate and divide the population into two halves.

    • Using the upper half of the population create offspring population

      • Selection, Crossover and Mutation

    • Perform Constraint Handling Mechanism

    • Evaluate the fuel cost objective function E and emission objective function F

    • Merge the parent and offspring population

    • Perform non domination sorting

    • Calculate crowding distance and rank based on non-domination fronts

    • Select solutions

      • Each front is filled in ascending order

      • Last front- descending order of crowding distance

    • Store the non-dominated solutions in list Ϝ1

    • Plot the non-dominated solutions in list Ϝ1

    • Position and cost of the particle are initialized from the lower half of the population

    • Store the values of the particles as their personal best pBest

    • Determine Domination

      • Initialize external repository rep

      • Create grid and find grid index

    • For each particle do

      • Select leader from external repository

      • Update particle position and velocity

      • Constraint Handling Mechanism

      • Evaluate the fuel cost objective function E and emission objective function F

      • Apply Mutation and calculate new solutions

      • Apply Constraint Handling Mechanism

      • Determine Domination pBest

    • End for

    • Add non dominated particles to the repository

    • Determine domination of new repository members

    • Keep only the non-dominated members in the repository

    • Update grid and grid index

    • Modify inertia weight

    • Create a new set of particles half the size nPop and fill it with the non-dominated solutions in the repository followed by the pBest

    • Combine the populations of NSGA II and the new set of particles of the MOPSO

    • Increment generation count

  • End for

Advertisement

6. Numerical tests

In order to validate the proposed hybrid algorithm, the CEED problem was solved for IEEE 30-bus system and the results are presented in this section. The fuel cost coefficients with valve-point loading, emission coefficients, and generator limits are adapted from [26] and is given in Table 1. The transmission loss B-matrix coefficients are obtained by running a load flow program and is in [26] is adapted here and given in Table 2. The total power demand in the system is 2.834p.u. to the base of100MVA. Program in MATLAB was developed for the Hybrid Algorithm to perform CEED and executed on1.60GHz, Intel T2050 processor, 1.5GB RAM HP Pavilion Laptop with WINDOWS 7 operating system. Various test cases are considered to compute the Pareto front of the multi-objective CEED problem. The Pareto-optimal front is obtained using the NSGA II algorithm and also using the MOPSO algorithm given in Section 4. The Pareto front obtained from the hybrid approach given in Section 5 is then compared with the Pareto front obtained using NSGAII and MOPSO algorithm.

Unit iGeneration LimitsFuel Cost Coefficients with valve point loadingEmission Coefficients
iminimaxaibicieifiαiβiγiηiδi
10.050.510200100156.2834.091−5.5546.4902e−42.857
20.050.6010150120108.9762.543−6.0475.6385e−43.333
30.051.0020180401014.7844.258−5.0944.5861e−68.000
40.051.201010060520.9445.326−3.5503.3802e−32.000
50.051.002018040525.1334.258−5.0944.5861e−68.000
60.050.6010150100518.4806.131−5.5555.1511e−56.667

Table 1.

Fuel costs Coefficients with valve point loading, Emission Coefficients, Generator limits of IEEE 30 bus system.

B0.021800.01070−0.00036−0.001100.000550.00330
0.010700.01704−0.00010−0.001790.000260.00280
−0.00040−0.000200.02459−0.01328−0.01180−0.00790
−0.00110−0.00179−0.013280.026500.009800.00450
0.000550.00026−0.011800.009800.02160−0.00010
0.003300.00280−0.007920.00450−0.000120.02978
B01.0731e−050.0017704−0.00406450.00384530.00138320.0055503
B000.0014

Table 2.

B−Loss Coefficients for IEEE 30 bus test system.

In case 1 the fuel cost function is modeled as a quadratic function with sine term to incorporate the valve-point effect. The transmission losses are also considered in this case. The Pareto front obtained using NSGA II, MOPSO, and Hybrid NSGAII-MOPSO is shown in Figures 1, 2 and 3 respectively. In all these figures there is a discontinuity in the Pareto front due to modeling of the valve point loading effect of generators.

Figure 1.

Pareto-optimal curve for IEEE 30 bus system obtained using NSGA II.

Figure 2.

Pareto-optimal curve for IEEE 30 bus system obtained using MOPSO.

Figure 3.

Pareto-optimal curve for IEEE 30 bus system obtained using Hybrid NSGAII and MOPSO Algorithm.

The parameter settings for NSGA II are obtained using trial and error is as follows: M=2; Population SizenPop=100; Maximum number of iterationMaxIt=100; Crossover Percentage pCrossover=0.7; Mutation PercentagepMutation=0.4; Mutation ratemu=0.02. The extreme points of the Pareto front and time for execution of NSGAII algorithm are provided in Table 3.

Method123456PlFuel Cost ($/h)Emission (Tons/h)Time Taken (s)
NSGA II0.06490.38660.68510.79990.53990.38860.03126616.4260.2121367
0.40700.45280.54160.41980.53650.50870.03279677.9410.1942
MOPSO0.06260.41060.68850.79940.54720.35640.03090618.2110.21251507
0.44120.45740.55010.38210.55230.48320.03242678.7020.1943
Hybrid NSGAII- MOPSO0.05000.38930.68610.80010.54900.39110.03178613.850.2127662
0.41090.45630.54290.40020.54350.51280.03279678.300.1942

Table 3.

Comparison of extreme points (shown in bold) and time taken for convergence using NSGAII, MOPSO and Hybrid NSGA II-MOPSO for IEEE30 bus system with valve point loading.

The parameter settings for MOPSO is obtained using trial and error is as follows: M=2; Maximum number of iterationMaxIt=500; Population SizenPop=250; Repository size nRep=100; Inertia Weightw=0.5; Inertia Weight damping ratewdamp=0.99; Personal learning coefficientc1=1; Global learning coefficient c2=2; Number of grids per dimension nGrid=10; Inflation Rate alpha=0.1, leader selection pressure beta=2, Deletion selection pressure gamma=2; Mutation rate mu=0.1. The extreme points of the Pareto front and time for execution of MOPSO algorithm are provided in Table 3. We can observe from Figure 2 and Table 3 that there are difficulties in MOPSO algorithm in obtaining well spread Pareto front and also very slow convergence to the Pareto front when compared to NSGA II. This can be improved if the proposed hybrid approach is used to solve the CEED problem.

The Parameter setting for the hybrid algorithm is same as those given above expect for the settings provided here Population SizenPop=200; Maximum number of iterationMaxIt=50; Repository sizenRep=20. The extreme points of the Pareto front and time for execution of the proposed NSGAII-MOPSO hybrid algorithm are provided in Table 3. From Table 3 it is clear that the extreme points found by the hybrid algorithm are better than NSGA II and MOPSO algorithm. Even though the time of execution of the Hybrid algorithm is slower than NSGA II it is able to find well spread Pareto front compared to NSGA II. The hybrid algorithm is far superior to MOPSO in terms of converge speed and also in finding well spread Pareto-optimal front.

In case II the valve point effect is neglected from the fuel cost curve and is solved using the proposed hybrid approach using the same parameters. The Pareto front obtained is shown in Figure 4 and is a continuous curve when compared to the Pareto front shown in Figure 3. In Figure 3 the Pareto front is discontinuous due to the effect of the Valve point loading in the cost curve. Both these case studies indicate that the hybrid approach is effective to solve the CEED problem.

Figure 4.

Pareto-optimal curve for IEEE 30 bus system without valve point effect obtained using hybrid NSGAII and MOPSO algorithm.

Advertisement

7. Conclusion

In this chapter, a hybrid multi-objective optimization algorithm based on NSGA II and MOPSO have been proposed to solve the highly nonlinear, highly constrained combined economic emission dispatch problem. At any stage of the algorithm, only feasible solution is created because of the incorporation of the proposed constraint handling mechanism. During every iteration of the hybrid algorithm new population is created and NSGA II is applied on best performing individuals whereas MOPSO is applied on the lower ranked individuals to strengthen the exploration and exploitation capability of the algorithm. This hybrid approach is tested on an IEEE 30 bus system. The results obtained shows that the hybrid approach is efficient for solving CEED problem and is also able to quickly converge to a better Pareto-optimal front when compared to MOPSO algorithm. The result obtained by the hybrid approach also demonstrates it is able to yield a wide spread of solutions and convergence to true Pareto-optimal fronts.

Advertisement

Acknowledgments

The author gratefully acknowledges the authorities of Royal Commission for Jubail and Yanbu and the authorities of Jubail Industrial College for the facilities offered to carry out the work. This work is part of the JRICH (Jubail Research and Innovation Cluster Hub) initiative of the Jubail Industrial College.

References

  1. 1. Wood AJ, Woolenberg BF. Power Generation, Operation and Control. 2nd ed. New York: John Wiley and Sons; 1996
  2. 2. Le KD, Golden JL, Stansberry CJ, Vice RL, Wood JT, Ballance J, Cauley GW. Potential impacts of clean air regulations on system operations. IEEE Transactions on Power Systems. 1995;10(2):647-656. DOI: 10.1109/59.387899
  3. 3. Deb K. Multi-Objective Optimization Using Evolutionary Algorithms. New York: Wiley; 2001. p. 518. ISBN: 978-0-471-87339-6
  4. 4. Abido M. A niched Pareto genetic algorithm for multi-objective environmental/economic dispatch. Electric Power and Energy System. 2003;25(2):97-105. DOI: 10.1016/S0142-0615(02)00027-3
  5. 5. Basu M. Dynamic economic emission dispatch using non dominated sorting genetic algorithm-II. International Journal of Electric Power and Energy Systems. 2008;30(2):140-149. DOI: 10.1016/j.ijepes.2007.06.009
  6. 6. Rajasomashekar S, Aravindhababu P. Biogeography based optimization technique for best compromise solution of economic emission dispatch. Swarm and Evolutionary Computation. 2012;7:47-57. DOI: 10.1016/j.swevo.2012.06.001
  7. 7. Zou D, Li S, Li Z, Kong X. A new global Particle Swarm Optimization for the economic emission dispatch with or without transmission losses. Energy Conversion and Management. 2017;139:45-70. DOI: 10.1016/j.swevo.2012.06.001
  8. 8. Panigrahi BK, Ravikumar Pandi V, Sanjay D, Das S. Multiobjective fuzzy dominance based bacterial foraging algorithm to solve economic emission dispatch problem. Energy. 2010;35(12):4761-4770. DOI: 10.1016/j.energy.2010.09.014
  9. 9. Silva MAC, Klein CE, Mariani VC, Coelho LS. Multiobjective scatter scearch approach with new combination scheme aplied to solve environmental/economic dispatch problem. Energy. 2013;53(5):14-21. DOI: 10.1016/j.energy.2013.02.045
  10. 10. Roy PK, Bhui S. Multi-objective quasi oppositional teaching learning based optimization for economic emission load dispatch problem. Electrical Power and Energy system. 2013;53(4):937-948. DOI: 10.1016/j.ijepes.2013.06.015
  11. 11. Abou El Ela AA, Abido MA, Spea SR. Differntial evolution algorithm for emission constrained economic power dispatch problem. Electric Power Systems Research. 2010;80(10):1286-1292. DOI: 10.1016/j.epsr.2010.04.011
  12. 12. Sivasubramani S, Swarup KS. Environemntal/economic dispatch using multi-objective harmony search algorithm. Electric Power Systems Research. 2011;81(9):1778-1785. DOI: 10.1016/j.epsr.2011.04.007
  13. 13. Arunachalam S, Saranya R, Sangeetha N. Hybrid artificial bee colony algorithm and simulated annealing algorithm for combined economic and emission dispatch including valve point effect. LNCS. 2013;8297(Part I):354-365. DOI: 10.1007/978-3-319-03753-0_32
  14. 14. Arunachalam S, Agnes Bhomila T, Ramesh Babu M. Hybrid Particle Swarm Optimization algorithm and firefly algorithm based combined economic and emission dispatch including valve point effect. LNCS. 2015;8947:647-660. DOI: 10.1007/978-3-319-20294-5_56
  15. 15. Palanichamy C, Babu NS. Analytical solution for combined economic and emission dispatch. Electric Power Systems Research. 2008;78(7):1129-1139. DOI: 10.1016/j.epsr.2007.09.005
  16. 16. Zhiyong L, Weiyou W, Yanyan Y, Li Z. PS-ABC: A hybrid algorithm based on particle swarm and artificial bee colony for high-dimensional optimization problems. Expert System with Applications. 2015;42:8881-8895. DOI: 10.1016/j.eswa.2015.07.043
  17. 17. Gherbi YA, Bouzeboudja H, Gherbi FZ. The combined economic environmental dispatch using new hybrid meta heuristics. Energy. 2016;115:468-477. DOI: 10.1016/j.energy.2016.08.079
  18. 18. Gong D-W, Zhang Y, Qi C-L. Environmental/economic power dispatch using hybrid multi-objective optimization algorithm. Electric Power and Energy System. 2010;32:607-614. DOI: 10.1016/j.ijepes.2009.11.017
  19. 19. Agarwal A, Nanavati N. Association rule mining using hybrid GA-PSO for multi-objective optimisation. In: IEEE International Conference on Computational Intelligence and Computing Research; Dec 2016; Agni College of Technology. Chennai: IEEE; 2017. DOI: 10.1109/ICCIC.2016.7919571
  20. 20. Gandelli A, Grimaccia F, Mussetta M, Pirinoli P, Zich RE. Development and validation of different hybridization strategies between GA and PSO. In: IEEE Congress on Evolutionary Computation; 25–28 September 2007; Singapore. IEEE; 2008. DOI: 10.1109/CEC.2007.4424823
  21. 21. Janga Reddy M, Nagesh Kumar D. An efficient multi-objective optimization algorithm based on swarm intelligence for engineering design. Engineering Optimization. 2007;39(1):49-68. DOI: 10.1080/03052150600930493
  22. 22. Kalyanmoy D, Pratap A, Agarwal S, Meyarivan T. A fast and elitist multi-objective genetic algorithm: NSGA II. IEEE Transaction on Evolutionary Computations. 2002;6(2):182-197. DOI: 10.1109/4235.996017
  23. 23. Coello Coello CA, Lechuga MS. MOPSO: A proposal for multiple objective Particle Swarm Optimization. In: Proceedings of the 2002 Congress on Evolutionary Computation; 12–17 May 2002; 07 August 2002. Honolulu: IEEE; 2002. DOI: 10.1109/CEC.2002.1004388
  24. 24. Chaing C-L. Improved genetic algorithm for power economic dispatch of units with valve-point effects and multiple fuels. IEEE Transactions on Power System. 2005;20(4):1690-1699. DOI: 10.1109/TPWRS.2005.857924
  25. 25. Srinivas N, Ded K. Multi-objective optimization using non-dominated sorting in genetic algorithm. Evolutionary Computation. 1994;2(3):221-248. DOI: 10.1162/evco.1994.2.3.221
  26. 26. Hemamalini S, Sishaj PS. Emission constrained economic dispatch with valve point effect using Particle Swarm Optimization. In: TENCON; 18–21 November 2008; Hyderabad. IEEE; 2009. DOI: 10.1109/TENCON.2008.4766473

Written By

Arunachalam Sundaram

Submitted: 27 May 2017 Reviewed: 29 November 2017 Published: 20 December 2017