Economic and Emission Load Dispatch Using Cuckoo Search Algorithm
Economic and Emission Load Dispatch Using Cuckoo Search Algorithm
Abstract—To examine the optimal power output for all the Raphson, Gradient based method, Lambda iteration and
generating units the economic load dispatch problem in a power Dynamic programming which doesn’t provide fast and global
system is used by which will decreases the total fuel cost whereas solution and the problem complexity increases[10]. Thus the
fulfilling all practical constraints. Since optimal economic power system optimization problem need algorithm with faster
dispatch is not globally the excellent solution so such it was rate of convergence, global solution provision and ability to
required to reduce the toxic waste or emission dispatch problem.
handle very large complexity. Various artificial intelligent
The Combined Economic Emission Dispatch problem is
formulated as bi-objective problem by considering both fuel cost techniques such as the Pattern Search (PS), Artificial bee
and emission concurrently. By using price penalty factor the bi- colony, Gravitational search algorithm, Genetic Algorithm
objective Economic and Emission Dispatch problem is converted (GA)[5,15,18], Evolution Programming (EP)[11], Particle
into single objective problem. Here, to resolve the combined Swarm Optimization[13], Firefly algorithm, Tabu search, have
economic-emission load dispatch problem an effort has been been effectively used to resolve power optimization
performed using Cuckoo search algorithm, where CSA have problem[17]. In this paper an appropriate and robust technique
number of application and it is a useful tool developed for is used to resolve economic dispatch problem i.e. Cuckoo
optimization in engineering. On the basis of the obligate brood Search Algorithm (CSA). It is a metaheuristic search
parasitic behavior of some cuckoo species this algorithm is
algorithm, developed by Yang and Deb in 2009[3,9].
explained. In this paper, Cuckoo Search Algorithm (CSA)
method is used to obtain the best optimal solution considering for II. FORMULATION OF ECONOMIC AND EMISSION DISPATCH
Combined Economic and Emission dispatch problem for simple
PROBLEM
six unit test system.
The main objective of the power system is to obtain the
Keywords—Combined Economic Emission Dispatch problem, optimal dispatch of power and minimization of the two
Cuckoo Search Algorithm, Price Penalty Factor, Best optimal objectives, i.e. fuel cost and emission content subject to
solution. equality and inequality constraints simultaneous. First the two
I.INTRODUCTION objectives are optimized separately by single objective
technique which is defined in below. Then the two objectives
As Unit Commitment & Load Forecasting, the Economic are concurrently optimized by using bi-objective optimization
Load Dispatch is also a crucial part of modern power technique which is also termed in below as formulation. The
system[1]. The ELD where purposely used to find the formulation of Combined Economic and Emission Load
optimum generation among the existing units, such that the Dispatch problem are given as
total cost of generating units is decreased while satisfying ൌ ሺσ୲ୀଵ σ୬୨ୀଵ ୨ ൫୨ ൯ ǡ୨ ሻ (1)
equality and inequality constraints in the system
concurrently[7]. For any load demand the ELD problem
Where,
provide power output of each unit in MW and also calculate
Pj = Real power output of the ith generation
the overall fuel cost of generating units[6]. The thermal plant
Fj(Pj) = Generation cost function
releases various pollutants in the atmosphere during
ET,j = Environmental cost function.
generating electricity[16]. Thus to reduce the fuel cost is not
Environmental cost is an additional cost added on the
the only single objective therefore in this paper we also
thermal generators by considering thermal units emission
consider emission constraints in order to reduce the
dispatch in the electricity market typical based on the cost
emissions[2,8,12]. The fuel cost curve is a smooth convex
function.
curve but due to insertion of valve point loading this fuel cost
A system consisting of n thermal generating units
curve become non-linear and the problem become more
connected to a transmission network with the electrical load
complex. Valve point loading is also considered in this
PD [MW].
paper[4,14]. Economic dispatch analysis has been studied by
many researchers using diverse methods. Previous methods
include various mathematical procedures such as Newton ୨ ൫୨ ൯ ൌ σ୬୨ୀଵሺ ୨ ୨ଶ ୨ ୨
୨ ሻ (2)
2
International Conference on Innovations in Power and Advanced Computing Technologies [i-PACT2017]
A. Levy flights 9: if fi < fj then
This algorithm is based by means of obligate brood 10: xj = xi; fj = fi; // replace j-th solution
parasitism regarding partial cuckoo species which lays their 11.end if
eggs among according to the host bird’s nest. Those female 12: if rand(0, 1) < pa then
parasitic cuckoo species can mimic the shade or pattern 13: init_ nest(x worst);
concerning the eggs of the host species. The initial solutions 14: end if
are represented as the existing egg in the host nest. Thus the 15: if fi < f min then
new solution generated in this algorithm represents by the egg 16: x best = xi; f min = fi; // save the local best sol.
laid by the cuckoo species. On the basis of the following three 17: end if
assumptions this algorithm works: 18: end for
• The egg laid by the cuckoo selects the host bird nest 19: end while
randomly and only one egg is laid by the cuckoo at a
time, C. Flowchart for Cuckoo Search Algorithm
• The best nest which has the high superiority egg
(solution) will carry over to the next generations. Start
• The host bird can conclude a cuckoo’s egg with a
probability and the host nests are fixed with the total
number of accessible nest, Initialization of the random population of n host
ܲܽ א0, 1.By switching probability (Pa) the third assumption
can be estimated and the total n nests get replaced by the new By Levy flights cuckoo egg are get randomly, i
nest which contains a new random solution. When generating
new solutionࢄ ࢚+ from levy flight is achieved as
ࢄ࢚+ = ࢄ࢚ + ࢻ ْ ࡸe࢜࢟ ࣅ Fitness is evaluated, Fi
where Į > 0 is the step size which are related to the problem of
interests with scales. In most cases, Į = 1 is used.
For random walk the stochastic equation is formed. randomly Select a nest among ,n
Generally, Markov chain is a random walk where depending
on the present location and the next location is obtained by the
transition probability. The product ْ represents the entry wise Yes
multiplications. A random walk was provided by the Levy Fi Fj
flight where the random step interval is drained from a Levy
distribution. No
Levy~ ݐ = ݑí, (1 < ߣ 3) which has an infinite mean and
Let j as a solution
variance. Here the random walk process formed by the steps
with a power law step size distribution with a heavy tail. The
algorithm which are stretched to an additional complicated Replace j by the new solution
case by which that the each nest contains multiple eggs
(solution)
CS can be characterized by the following three features:
• It satisfies the global convergence requirements. Abandon a fraction, Pa of worst nests and build new
• It supports local and global search capabilities. ones at new locations via Levy flights
• It uses Levy flights as a global search strategy.
B. Pseudo Code for Cuckoo Search Algorithm The current best is fixed
The pseudo code obtained by the summarization of the
algorithm
Input: Population of nests xi = (xi1, …, xi D)T for i =1 … n,
MAX FE. No
Output: The best solution x and its equivalent value f min = j max
min(f(x)). iterations
1: generate initial host nest locations();
2: eval = 0;
Yes
3: while termination condition not meet do
4: for i = 1 to n do Find the best objective (the best nest)
5: xi = generate new solution(xi);
6: fi = evaluate the new solution(xi);
7: eval = eval + 1; End
8: j = ۏrand(0, 1) כn + 1;ے
3
International Conference on Innovations in Power and Advanced Computing Technologies [i-PACT2017]
IV. RESULT AND DISCUSSION
Thus the proposed approach for resolving the bi-objective 7420
ECONOMIC LOAD DISPATCH FOR DEMAND 600MW
along single load require and for load varied for every hour on 7380
7360
/hr)
7340
($
units system for fuel cost ($/h) only using cuckoo search
OST
7320
C
algorithm together with genetic algorithm approach which are 7300
7240
concerning for simple six- generator test system which are Fig. 1. Convergence characteristics of fuel cost using cuckoo search algorithm
(Demand=600MW)
explained on case2 results. In this paper, the maximum price
penalty factor has been selected because combining fuel cost in ECONOMIC LOAD DISPATCH FOR DEMAND 800MW
/hr)
applied in accordance with six-generator systems for load
9500
O ST($
varied for every hour on one day and also consider the cost
C
9498
9494
CASE 1 9492
0 10 20 30 40
GENERTION
50 60 70
1.1933
r)
1.1932
O
C S
1.1931
1.193
Load Demand
Economic dispatch Economic dispatch
/hr)
1.5355
O ST($
($/h) ($/h)
C
1.535
1.5345
1000 MW 11943 11938 Fig. 4. Convergence characteristics of fuel cost using cuckoo search algorithm
(Demand=1263MW)
1263 MW 15346 15338
4
x 10 ECONOMIC LOAD DISPATCH FOR DEMAND 1450MW
1.791
1.79
From Table 2 Shows that the fuel cost attained from the
/hr)
1.7895
($
1.789
4
International Conference on Innovations in Power and Advanced Computing Technologies [i-PACT2017]
CASE 2
Table 3 Fuel cost, emission coefficients and generation limits of six-generator system
Gen. no
aj bj cj ߙj ߚj ߛ݆ min load max load
($/MW2h) ($/MW h) ($/h) (Kg/MW2h) (Kg/MW h) (Kg/h) (MW) (MW)
1 0.1525 38.540 756.800 0.0042 0.33000 13.86 10 125
2 0.1060 46.160 451.325 0.0042 0.33000 13.86 10 150
3 0.0280 40.400 1050.00 0.0068 -0.5455 40.267 35 225
4 0.0355 38.310 1243.53 0.0068 -0.5455 40.267 35 210
5 0.0211 36.328 1658.57 0.0046 -0.5112 42.90 130 325
6 0.0180 38.270 1356.66 0.0046 -0.5112 42.90 125 315
Table 4 Comparison between the proposed method for different PPF types at
700MW
Total
PPF Emission Emission Economic
Method cost
($/kg) (kg/h) cost ($/h) cost ($/h)
($/h)
fi max Cuckoo 435 20,929 36,389 57,317
fi min Search 424 37,524 37,828 75,352
fi avg Algorithm 421 31,046 37,700 68,746
Method
fi com 432 30,601 37,080 67,682
CASE 3
Table 5 Fuel cost characteristics coefficient for the 6 conventional thermal
generators
Gen no ߙj ߚj ߛj ߜj ߠj
1 4.073 -5.112 5.534 0.0021000 3.876
2 3.216 -5.542 6.123 0.0005200 4.223
3 3.972 -4.521 5.785 0.0065000 4.756
4 4.333 -6.281 4.432 0.0000370 5.789
5 4.112 -5.776 6.395 0.0043000 6.773 Fig. 8. Convergence characteristics of emission cost obtained for 6
6 4.112 -6.043 5.295 0.0005600 3.156 unit system considering price penalty factor (hi max)
5
International Conference on Innovations in Power and Advanced Computing Technologies [i-PACT2017]
Table 7 Total generation cost on Dth day using the Cuckoo Search Algorithm
Time(hour) Load(MW) P1(MW) P2(MW) P3(MW) P4(MW) P5(MW) P6(MW) COST ($/h)
1 308 46.30 45.80 98.80 36.0 45.5 35.7 4,662.1
2 297 68.30 48.60 50.00 39.7 55.4 35.0 4,381.1
3 295 29.90 78.90 46.00 58.2 38.2 43.8 4,533.9
4 302 40.60 69.80 36.20 75.5 44.8 35.0 4,433.3
5 308 72.00 35.00 94.40 35.0 35.0 36.6 4,497.7
6 310 38.00 38.70 113.9 35.0 49.4 35.0 4,909.7
7 350 131.6 57.00 46.00 44.6 35.8 35.0 4,874.7
8 430 92.30 89.80 48.90 80.1 40.2 78.7 6,372.3
9 460 137.0 100.1 35.00 69.1 83.9 35.0 6,598.6
10 465 139.8 102.5 57.60 82.2 41.9 41.1 6,677.5
11 562 95.80 108.3 101.5 75.9 39.8 40.7 6,659.4
12 455 156.4 66.90 80.30 81.3 35.0 35.0 6,549.7
13 430 95.20 104.0 35.50 79.3 73.6 42.5 6,172.6
14 432 139.4 99.00 40.10 35.7 82.7 35.0 6,079.1
15 456 151.7 94.40 77.80 35.1 60.3 36.7 6,566.8
16 452 153.6 125.5 41.30 47.5 47.7 36.3 6,444.7
17 446 125.5 104.6 105.3 37.7 37.8 35.0 6,297.0
18 435 120.8 116.7 45.20 74.5 42.5 35.2 6,141.8
19 426 153.9 108.8 35.00 35.0 35.0 58.3 5,947.7
20 407 104.1 110.1 43.00 79.4 35.0 35.4 5,584.2
21 370 96.80 64.60 66.80 35.0 65.4 41.4 5,412.2
22 350 88.30 94.70 45.20 35.7 45.5 40.6 4,878.4
23 320 53.50 108.7 38.90 35.0 46.0 37.9 4,512.3
24 300 47.50 54.70 57.40 46.6 55.8 38.0 4,615.7
TOTAL COST (FUEL +EMISSION) 1,33,800
6
International Conference on Innovations in Power and Advanced Computing Technologies [i-PACT2017]
[12] Ramarathnam,R., “Emission constrained economic dispatch,” IEEE
Transactions on Power Systems, Vol. 9, No. 4, pp. 1994–2000, Nov1994
[13] Revathy K, Nithiyanandham, “Economic dispatch using particle swarm
optimization,” International Journal of Advanced Research in Electrical,
Electronics and Instrumentation Engineering, vol. 3, Special Issue 4, pp:59-
66, May 2014.
[14] Chandrakala G, Ramakrishna V, Shaik J, “Solving Economic Dispatch
with Valve Point Loading Effects by using Optimization Techniques”,
International Journal of Engineering Science Invention, Vol. 3, Issue 9, pp:01-
11,Sep. 2014.
[15] Chiragkumar K. Patel, Mihir B. Chaudhari, “Economic load dispatch
problem using genetic algorithm,” Indian Journal of Applied Research, vol.4,
Issue 11,pp:173-175, November 2014.
[16] Lalita Rani, Ms. Manju Mam, Sanjeev Kumar,” Economic Load Dispatch
In Thermal Power Plant Taking Real Time Efficiency As An Additional
Constraints” International Journal of Engineering Research & Technology
(IJERT) Vol. 2 Issue 7, July 2013
[17] Balaji M, Naresh G, ''Solution to CEED Problem using Classical
Approach of Lagrange’s and PSO Algorithms'', International Journal of
Innovations in Engineering and Technology, Volume 4, Issue 4, pp:13-19,
Dec. 2014.
[18] Koridak L, Rahli M, Younes M,''Hybrid Optimization of the Emission
and Economic Dispatch by the Genetic Algorithm'', Leonardo Journal of
Sciences, Issue 14, January, pp: 193-203,June 2008.