Constrained Evolutionary Programming Approaches to Power System Economic Dispatch

Similar documents
EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

The Study of Teaching-learning-based Optimization Algorithm

Evolutionary Computational Techniques to Solve Economic Load Dispatch Problem Considering Generator Operating Constraints

Single-Facility Scheduling over Long Time Horizons by Logic-based Benders Decomposition

Self-adaptive Differential Evolution Algorithm for Constrained Real-Parameter Optimization

Energy Conversion and Management

Solving Nonlinear Differential Equations by a Neural Network Method

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

Research Article Multiobjective Economic Load Dispatch Problem Solved by New PSO

Economic dispatch solution using efficient heuristic search approach

Global Optimization Using Hybrid Approach

Combined Economic Emission Dispatch Solution using Simulated Annealing Algorithm

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Multi-Objective Evolutionary Programming for Economic Emission Dispatch Problem

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach

An Admission Control Algorithm in Cloud Computing Systems

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Differential Evolution Algorithm with a Modified Archiving-based Adaptive Tradeoff Model for Optimal Power Flow

VQ widely used in coding speech, image, and video

Lecture Notes on Linear Regression

Chapter 2 Real-Coded Adaptive Range Genetic Algorithm

Population Variance Harmony Search Algorithm to Solve Optimal Power Flow with Non-Smooth Cost Function

Solving of Single-objective Problems based on a Modified Multiple-crossover Genetic Algorithm: Test Function Study

CHAPTER 7 STOCHASTIC ECONOMIC EMISSION DISPATCH-MODELED USING WEIGHTING METHOD

HYBRID FUZZY MULTI-OBJECTIVE EVOLUTIONARY ALGORITHM: A NOVEL PARETO-OPTIMIZATION TECHNIQUE

Some modelling aspects for the Matlab implementation of MMA

Optimal choice and allocation of distributed generations using evolutionary programming

Lecture 14: Bandits with Budget Constraints

A SEPARABLE APPROXIMATION DYNAMIC PROGRAMMING ALGORITHM FOR ECONOMIC DISPATCH WITH TRANSMISSION LOSSES. Pierre HANSEN, Nenad MLADENOVI]

Markov Chain Monte Carlo Lecture 6

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

The Minimum Universal Cost Flow in an Infeasible Flow Network

Using Immune Genetic Algorithm to Optimize BP Neural Network and Its Application Peng-fei LIU1,Qun-tai SHEN1 and Jun ZHI2,*

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

The Convergence Speed of Single- And Multi-Objective Immune Algorithm Based Optimization Problems

ECE559VV Project Report

Determining Transmission Losses Penalty Factor Using Adaptive Neuro Fuzzy Inference System (ANFIS) For Economic Dispatch Application

Problem Set 9 Solutions

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

Computing Correlated Equilibria in Multi-Player Games

Chapter Newton s Method

C-GRASP APPLICATION TO THE ECONOMIC DISPATCH PROBLEM

Particle Swarm Optimization with Adaptive Mutation in Local Best of Particles

STRATEGIES TO IMPROVE THE DIRECT SEARCH METHOD FOR NON-CONVEX WIND-THERMAL COORDINATION DISPATCH PROBLEM CONSIDERING TRANSMISSION CAPACITY

Applied Energy 86 (2009) Contents lists available at ScienceDirect. Applied Energy. journal homepage:

Kernel Methods and SVMs Extension

GHHAGA for Environmental Systems Optimization

A Novel Evolutionary Algorithm for Capacitor Placement in Distribution Systems

Chapter 2 A Class of Robust Solution for Linear Bilevel Programming

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business

Optimal Solution to the Problem of Balanced Academic Curriculum Problem Using Tabu Search

Markov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

E O C NO N MIC C D I D SP S A P T A C T H C H A N A D N D UN U I N T T CO C MMITM T EN E T

Lecture 20: November 7

Real-Time Systems. Multiprocessor scheduling. Multiprocessor scheduling. Multiprocessor scheduling

Chapter - 2. Distribution System Power Flow Analysis

ECONOMIC POWER DISPATCH USING THE COMBINATION OF TWO GENETIC ALGORITHMS

Tornado and Luby Transform Codes. Ashish Khisti Presentation October 22, 2003

Evolutionary Multi-Objective Environmental/Economic Dispatch: Stochastic vs. Deterministic Approaches

An Integrated Asset Allocation and Path Planning Method to to Search for a Moving Target in in a Dynamic Environment

Generalized Linear Methods

Static security analysis of power system networks using soft computing techniques

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

On the Multicriteria Integer Network Flow Problem

A Hybrid Variational Iteration Method for Blasius Equation

EP BASED PSO METHOD FOR SOLVING PROFIT BASED MULTI AREA UNIT COMMITMENT PROBLEM

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique

ECONOMIC DISPATCH SOLUTION USING A REAL-CODED GENETIC ALGORITHM

CHAPTER 2 MULTI-OBJECTIVE GENETIC ALGORITHM (MOGA) FOR OPTIMAL POWER FLOW PROBLEM INCLUDING VOLTAGE STABILITY

An Improved Shuffled Frog Leaping Algorithm Approach for Unit Commitment Problem

An improved multi-objective evolutionary algorithm based on point of reference

Transient Stability Constrained Optimal Power Flow Using Improved Particle Swarm Optimization

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM

Wavelet chaotic neural networks and their application to continuous function optimization

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

RECENTLY, the reliable supply of electric power has been

arxiv: v1 [math.oc] 3 Aug 2010

An Improved Clustering Based Genetic Algorithm for Solving Complex NP Problems

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

UNIT COMMITMENT SOLUTION USING GENETIC ALGORITHM BASED ON PRIORITY LIST APPROACH

Thin-Walled Structures Group

Capacitor Placement In Distribution Systems Using Genetic Algorithms and Tabu Search

A HYBRID DIFFERENTIAL EVOLUTION -ITERATIVE GREEDY SEARCH ALGORITHM FOR CAPACITATED VEHICLE ROUTING PROBLEM

Stochastic Weight Trade-Off Particle Swarm Optimization for Optimal Power Flow

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Comparative Analysis of SPSO and PSO to Optimal Power Flow Solutions

A New Evolutionary Computation Based Approach for Learning Bayesian Network

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL

Economic Load Dispatch with Nonsmooth Cost Functions Using Evolutionary Particle Swarm Optimization

Semi-supervised Classification with Active Query Selection

An Extended Hybrid Genetic Algorithm for Exploring a Large Search Space

Lecture 10 Support Vector Machines II

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances

Global Optimization of Truss. Structure Design INFORMS J. N. Hooker. Tallys Yunes. Slide 1

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Transcription:

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) Constraned Evolutonary Programmng Approaches to Power System Economc Dspatch K. Shant Swarup Abstract - Ths paper proposes a novel methodology of Constrant evolutonary programmng for solvng dynamc economc dspatch. Dynamc Economc Dspatch s one of the man functons of power generaton operaton and control. It determnes the optmal settngs of generator unts wth predcted load demand over a certan perod of tme. The objectve s to operate an electrc power system most economcally whle the system s operatng wthn ts securty lmts. Ten unts test system wth smooth and nonsmooth fuel cost functons are consdered to llustrate the sutablty and effectveness of the proposed method. Keywords: Constraned Optmzaton, Evolutonary Programmng, Stochastc Optmzaton technques, Power system operaton and control, Dynamc Economc Dspatch. I INTRODUCTION Dynamc Economc Dspatch s a method to schedule the onlne generator outputs wth the predcted load demands over a certan perod of tme so as to operate an electrc power system most economcally [1-4]. It s a dynamc optmzaton problem takng nto account the constrants mposed on system operaton by generatng rampng rate lmts. The dynamc economc dspatch s not only the most accurate formulaton of the economc dspatch problem but also the most dffcult to solve because of ts large dmensonalty. Normally t s solved by dvdng the entre dspatch perod nto a number of small tme ntervals, and then a statc economc dspatch has been employed to solve the problem n each nterval. However all of those methods may not be able to provde an optmal soluton and usually gettng stuck at a local optmal. Recently, stochastc optmzaton technques [5-15] such as Smulated Annealng (SA), Genetc Algorthm (GA) and Evolutonary Programmng (EP) have been gven much attenton by researchers due to ther ablty to seek for the near global optmal soluton. Approprate settng of the control parameters of the SA algorthm s a dffcult task and the speed of the algorthm s slow when appled to a real power system. Genetc Algorthm, nvented by Holland [15] n the early 1970s, s a stochastc global search method that mmcs the metaphor of natural bologcal evaluaton. The GA uses only the objectve functon nformaton, not dervatve or other auxlary knowledge. Therefore GA can deal wth the non-smooth, non-contnuous and non dfferentable functons, whch actually exst n a practcal optmzaton problem. Evolutonary Programmng (EP) ntroduced by Lawrence J. Fogel n the 1960s, on the other hand, s also a global search method startng from a populaton of canddate solutons, and fnds soluton n parallel usng evaluaton process. Both GA and EP can provde a near global soluton. However the encodng and decodng schemes essental n the GA approach are not needed n dynamc economc dspatch problem. Therefore EP s faster n speed than GA n ths case. Any evolutonary computaton technque appled to a partcular problem should address the ssue of handlng nfeasble ndvduals. The presence of constrants sgnfcantly affects the performance of optmzaton algorthm, ncludng evolutonary search methods. The general way of handlng constrants s by penalzng the nfeasble ponts. Constrant handlng n evolutonary computaton s more or less problem dependent [21-23]. In ths paper, the applcaton of constraned optmzaton usng Evolutonary programmng technque for solvng dynamc economc dspatch has been proposed. The proposed method starts wth only feasble canddate solutons generated by an teratve procedure of correctng volatons. In order to show the effectveness of the proposed method, 10 unts test system wth smooth and non-smooth fuel cost functons are consdered n ths paper. The results of the proposed method are compared wth those reported n lterature [18-20]. II PROBLEM FORMULATION The dynamc economc dspatch problem can be formulated as follows. The objectve functon: Where, T Mn F = ( ) N t= 1 = 1 F t P t (1)

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) F s the total operatng cost over the whole dspatch perods, T s the number of hours n the tme horzon, N s the number of dspatchable unts, F P ) s the ndvdual generaton producton t ( t cost n terms of real power output P at any tme t. The fuel cost functon of each unt wth smooth cost functon can be expressed as 2 F ( P ) = a P + b P + c (2) The fuel cost functon of each unt wth valve pont loadng effects can be expressed as 2 F P a P + b P + c + e sn f P ( ) ( ( P )) = mn Where a, b, c, e and f are constants of fuel cost functon of the unt. Subject to: 1) Real power balance constrant N = 1 P j P Dj = 0 (3) Where PDj s the total assumed load demand durng j th nterval. 2) Real power operatng lmts P P P, = 1,..., N (4) mn max Where P mn and P max are the mnmum and the maxmum real power outputs of the th generator respectvely. 3) Generatng unt ramp rate lmts Pt P( t 1) UR P( t 1) Pt DR = 1,..., N t = 1,..., T (5) Where UR and DR are the ramp-up and rampdown lmts of the th generator respectvely. III METHODOLOGY Evolutonary Programmng s a near global stochastc optmzaton method whch places emphass on the behavoral lnkage between parents and ther offsprng, rather than seekng to emulate specfc genetc operators as observed n nature to fnd a soluton. In order to obtan hgh qualty soluton, constraned optmzaton procedure has been adopted n ths paper. Constraned Optmzaton usng Evolutonary Programmng method starts from multple feasble ponts. 1) Representaton For T ntervals n the generaton schedulng horzon, there are T dspatched for the N generators. An array of control varable vectors can be shown as: P11P12... P1 T P21P22... P2 T (6)... PN 1... PNT Where Pt s the real power output of generator at ntervalt. 2) Intalzaton In order to avod usng penaltes and unnecessarly searchng n nfeasble regon, the ntal populaton s flled wth only feasble canddate solutons generated by an teratve procedure of correctng constrant volatons. 3) Ftness Evaluaton Snce the constrants are always satsfed by every member at any pont of evoluton only the objectve functon s used to calculate the cost of each canddate soluton, whch s taken as ts ftness value. 4) Creaton of offsprng A new populaton of soluton (Offsprng) s produced from the exstng populaton by selectng two dfferent generators randomly and at a randomly selected tme and transferrng some of the power from one generator to the other n such a way that the constrants are not volated. The amount to be transferred s chosen by the followng procedure: a) Calculate how much power can be removed from frst randomly selected generator so that t wll not volate ts constrants (both real power operatng lmts and generatng unt ramp rate lmts) and take a random fracton of t. b) Calculate how much power can be added to the second randomly selected generator so that t wll not exceed ts constrants and take a random fracton of t. c) Select the smaller of the two fractons and transfer t from the frst to the second generator. Snce both are n the same tme nterval power balance constrant wll not be volated. 5) Selecton and Competton The selecton process used n ths paper s the stochastc tournament method. The 2p ndvduals compete wth each other for selecton. p potental canddates are selected for the next generaton. IV COMPUTATIONAL PROCEDURE:

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) The pseudo-code for the proposed method can be descrbed as follows: Step 1: Read the system data from the nput fle. Step 2: Intalze the populaton Step 3: Check for feasblty and apply repar operator f populaton s not feasble. Step 3: Evaluate ftness functon. Step 4: Creaton of offsprng by mutaton. Step 5: Competton and selecton. Step 6: If termnaton not reached loop to Step 3. Step 7: Get the fnal result and qut the program. The above strateges are llustrated clearly n Fg 1. V SIMULATION RESULTS The 10 unts system wth smooth and non-smooth Fg 1. Flow Chart of Constraned Optmzaton. fuel cost functon s used n ths paper to demonstrate the performance of the proposed method. Unt data can be found n Appendx. The smulatons were carred out on a Pentum III 550MHz processor. Case 1: In ths case, 10 unts test system wth smooth fuel cost functon was used (Tables A.1 and A.2). The demand of the system was dvded nto 12 ntervals. In ths smulaton, populaton sze s 20 and the algorthm wll stop when 10000 generatons are reached. The generaton schedule (best soluton) for the 10 unts test system wth smooth fuel cost functon s gven n Table 1. The smulaton result of the proposed method compared to those obtaned from (Lnear Programmng) LP[19] and (Evolutonary Programmng) EP (Drect) [19] s gven n Table 2. Hgh qualty soluton can be obtaned at the expense of computng tme by ncreasng the populaton sze. The performance of the proposed method for 10-unts 12-hour demand system s gven n Fg 2. The contrbuton of each unt towards the load s gven n Fg 3. Case II: In ths case, 10 unts test system wth non-smooth fuel cost functon was used 9Tables A.3 and A.4). The demand of the system was dvded nto 24 ntervals. In ths smulaton, a populaton sze of 20 was chosen and the algorthm wll stop after 30000 teratons. The generaton schedule (best soluton) for the 10 unts test system wth non-smooth fuel cost functon s gven n Table 3. The smulaton result of the proposed method compared to those obtaned from (Sequental Quadratc Programmng) SQP and (Evolutonary Programmng) EP [20] s gven n Table 4. The performance of the proposed method for 10-unts 24-hour demand system s gven n Fg 4. The contrbuton of each unt towards the load s gven n Fg 5. VI DISCUSSION In the past, to solve the economc dspatch problem, most algorthms requre the ncremental cost curves to be of monotoncally smooth ncreasng nature and contnuous. The conventonal methods gnore or flatten out the portons of the ncremental cost curve that are monotoncally ncreasng or not contnuous. Approaches wthout restrcton on the shape of fuel cost functons are needed to gve accurate results. The constraned optmzaton usng Evolutonary Programmng approach has proved to be a very effcent method to solve the Dynamc Economc Dspatch problem wth smooth and non-smooth fuel cost functons. EP seeks a soluton wth exponental convergence rate,.e., the speed of convergence wll be very fast at the begnnng of the smulaton and wll slow down when the generaton ncreases. When coupled wth Constraned Optmzaton usng effcent constrant handlng technques, the algorthm gves an optmal soluton at a faster rate. The advantage of the proposed method s that, fuel cost functon wth any arbtrary shape can be dealt wth. Therefore, the proposed methodology can be appled to the system wth non-smooth, nondfferentable and non-contnuous functon, whch exsts n any real power system. VII CONCLUSION Ths paper presents a new methodology whch uses effcent constrant handlng technques wthout adoptng penalty methods for obtanng

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) the global soluton of the Dynamc Economc Dspatch problem. Smulaton results demonstrate that the proposed method gve a cheaper total producton cost than those obtaned by other classcal technques and Evolutonary Paradgms. Ten unts system wth smooth and non-smooth fuel cost functons were used n ths paper. Increase n the complexty of the objectve functon won t affect the complexty of the proposed method. VIII REFERENCES [1] A.J. Wood, B. F. Wollenberg, Power Generaton Operaton and Control, 2 nd ed., New York: John Wley & Sons, Inc., 1996. [2] F. LI, R. Morgan and D. Wllams, Hybrd Genetc Approaches to rampng rate constraned Dynamc Economc Dspatch, Electrc Power Systems Research, vol. 43, no. 2, pp. 97-103, Nov.1997. [3] Q. Xa, Y. H. Song, BM. Zhang, CQ. Kang and N. Xang, Dynamc queung approach to power system short term economc and securty dspatch, IEEE Trans. Power Systems, vol. 13, no. 2, pp. 280-285, May. 1998. [4] X. S. Han, H. B. Goo and D. S Krschen, Dynamc Economc Dspatch: Feasble and Optmal Solutons, IEEE Trans. Power Systems, vol. 16, no. 1, pp. 22-28, Feb. 2001. [5] D. B. Fogel, Evolutonary Computaton, Toward a New Phlosophy of Machne Intellgence, 2nd ed., New York: IEEE Press, 2000. [6] Z. Mchalewcz, Genetc Algorthm + Data Structures = Evoluton Programs, Sprnger-Verlag, Berln Hedelberg, 1992. [7] D.E. Goldberg, Genetc Algorthm n Search, Optmzaton and Machne Learnng, Addson Wesley, Readng, 1989. [8] T. Back, Evolutonary Algorthm n Theory and Practce, New York: Oxford Unv. Press, 1996. [9] V. Mranda, D. Srnvasan and L.M. Proenca, "Evolutonary computaton n power systems," Electrcal Power & Energy Systems, vol. 20, No. 2, pp. 89-98, 1998. [10] L.J. Fogel, A.J. Owens and M.J. Walsh, Artfcal Intellgence through Smulated Evoluton, New York: John Wley & Sons, Inc., 1966. [11] K.P. Wong and C.C. Fung, "Smulated annealng based economc dspatch algorthm," IEE Proc. C, vol. 140, no. 6, pp. 509-515, Nov. 1993. [12] H.T. Yang, P.C. Yang and C.L. Huang, "Evolutonary Programmng based economc dspatch for unts wth non-smooth fuel cost functons," IEEE Trans. Power Systems, vol. 11, no. 1, pp. 112-118, Feb. 1996. [13] D. Walters, C. Sheble and B. Gerald, "Genetc algorthm soluton of economc dspatch wth valve pont loadng," IEEE Trans. Power Systems, vol. 8, no. 3, pp. 1325-1332, Aug. 1993. [14] K.P. Wong and J. Yuryevch, "Evolutonary- Programmng-Based Algorthm for Envronmentally-Constraned Economc Dspatch," IEEE Trans. Power Systems, vol. 13, no. 2, pp. 301-306, May. 1998. [15] J.H. Holland, Adaptaton n Natural and Artfcal Systems, The Unversty of Mchgan Press, Ann Arbor, USA, 1975. [16] H. Duo, H. Sasak, T. Nagata and H. Fujta, "A soluton for unt commtment usng Lagrangan relaxaton combned wth evolutonary programmng," Electrc Power Systems Research, vol. 51, pp. 71-77, 1999. [17] F. L, R. Morgan and D. Wllams, "Hybrd genetc approaches to rampng rate constraned dynamc economc dspatch," Electrc Power Systems Research, vol. 43, pp. 97-103, 1997. [18] P. Attavryanupap, H. Kta, E. Tanaka and J. Hasegawa, "A hybrd method of EP and SQP for Dynamc Economc Dspatch," n 2001 Natonal Conventon Records (6), IEE Japan, pp. 2275-2276. [19] P. Attavryanupap, H. Kta, E. Tanaka and J. Hasegawa, "Dynamc economc dspatch usng evolutonary programmng combne wth sequental quadratc programmng," n The Intellgent Systems Applcaton to Power Systems Conference 2001, pp. 137-141. [20] P. Attavryanupap, H. Kta, E. Tanaka and J. Hasegawa, A Hybrd EP and SQP for Dynamc Economc Dspatch wth Non- Smooth Fuel Cost Functon, IEEE Transactons on Power Systems, vol. 17, no. 2, pp. 411-416, May 2002. [21] Zbgnew Mchalewcz. A Survey of Constrant Handlng Technques n Evolutonary Computaton Methods. In J. R. McDonnell, R. G. Reynolds and D. B. Fogel, edtors, Proceedngs of the 4th Annual Conference on Evolutonary Programmng, pages 135-155. The MIT Press, Cambrdge, Massachusetts, 1995. [22] Zbgnew Mchalewcz and Cezary Z. Jankow. Handlng Constrants n Genetc Algorthms. In R. K. Belew and L. B. Booker, edtors, Proceedngs of the Fourth Internatonal Conference on Genetc Algorthms, pages 151-157, San Mateo, Calforna, 1991. Morgan Kaufmann Pub. [23] Zbgnew Mchalewcz and Marc Schoenauer. Evolutonary Algorthms for Constraned Parameter Optmzaton Problems. Evolutonary Computaton, 4(1):1-32, 1996.

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) Anand N receved B.Tech. Degree n Electrcal Engneerng n 2003 from IIT Madras. Hs research nterests are Evolutonary Computaton, Optmzaton technques and ther applcaton to Power Systems,. K.S. Swarup s wth the Department of Electrcal Engneerng, Indan Insttute of Technology, Madras, Inda. Pror to jonng the department as vstng faculty, he held postons at the Mtsubsh Electrc Corporaton, Osaka, Japan, and Ktam Insttute of Technology, Hokkado, Japan, as vstng research scentst and vstng professor, respectvely durng 1992 to 1999. Hs areas of research are AI, knowledge-based systems, computatonal ntellgence, soft computng and Optmzaton Technques for electrc power systems. 2.32 x 106 Contrbuton of each unt towards total load 2.3 2.28 2.26 2.24 2.22 2.2 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 Fg 2. Performance of CEP for 10-unts 12-hour demand test system Generaton 7000 6000 5000 4000 3000 2000 1000 0 1 2 3 4 5 6 7 8 9 10 11 12 Hours T/U P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 Fg 3. Load dstrbuton among unts for 12-hour demand 1.07 x 106 1.065 1.06 1.055 1.05 1.045 1.04 1.035 1.03 Generaton Contrbuton of each unt towards total load 2500 2000 1500 1000 500 0 1 3 5 7 9 11 13 15 17 19 21 23 Hours 1.025 0 0.5 1 1.5 2 2.5 3 x 10 4 Unts/Hours P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 Fg 4. Performance of CEP for 10-unts 24-hours Fg 5. Load dstrbuton among unts for 24-hours TABLE 1. GENERATION SCHEDULE FOR 10-GENERATING UNIT 12-HOUR SYSTEM GENERATING UNITS Hours 1 2 3 4 5 6 7 8 9 10 Load (MW) 1 238.5 354.9 411.9 475.6 415.6 546.5 620.0 643.0 908.0 946.1 5560 2 243.5 361.1 417.4 481.9 421.3 555.7 620.0 643.0 918.9 957.3 5620 3 261.7 381.1 439.2 505.3 442.3 588.5 620.0 643.0 920.0 998.9 5800 4 262.6 382.7 440.3 506.6 443.4 590.2 620.0 643.0 920.0 1001.2 5810 5 280.9 402.7 462.2 530.3 464.5 623.1 620.0 643.0 920.0 1043.2 5990 6 286.6 409.9 469.1 537.7 471.1 633.5 620.0 643.0 920.0 1050.0 6041 7 282.0 404.5 463.5 531.7 465.6 625.1 620.0 643.0 920.0 1045.6 6001 8 260.6 380.1 437.9 504.1 441.1 586.6 620.0 643.0 920.0 996.6 5790 9 249.5 367.4 424.6 489.6 428.3 566.5 620.0 643.0 920.0 971.1 5680 10 236.9 353.1 409.7 473.4 413.8 543.9 620.0 643.0 904.1 942.1 5540 11 250.5 368.6 425.8 491.0 429.4 568.3 620.0 643.0 920.0 973.4 5690 12 256.6 375.5 433.0 498.7 436.4 579.4 620.0 643.0 920.0 987.4 5750

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) TABLE 2. COMPARISON OF TOTALPRODUCTION COST OBTAINED BY CONSTRAINED EVOLUTIONARY PROGRAM Methods Lnear Programmng [19] Evol. Prog (Drect) [19] Constraned EP Total Cost ($) 2196939 2196608 2196500 TABLE 3. GENERATION SCHEDULE FOR 10-UNIT 24-HOUR SYSTEM GENERATING UNITS Tme LOAD 1 2 3 4 5 6 7 8 9 10 (hrs) (MW) 1 150.00 135.00 73.00 60.00 172.60 123.29 129.75 85.31 52.04 55.00 1036 2 226.61 135.58 73.00 108.29 122.87 122.35 129.59 85.27 51.43 55.00 1110 3 226.63 215.58 108.58 120.54 172.87 122.45 129.59 85.32 21.43 55.00 1258 4 226.62 222.25 185.16 170.54 172.67 138.79 129.66 85.32 20.00 55.00 1406 5 226.65 222.27 213.59 180.97 222.67 122.63 129.63 85.31 21.28 55.00 1480 6 303.25 222.28 293.59 130.97 222.60 134.13 129.60 85.31 51.28 55.00 1628 7 379.87 302.28 293.38 120.31 222.58 122.41 129.59 55.31 21.28 55.00 1702 8 456.50 316.80 285.84 120.28 222.59 122.40 129.59 47.00 20.00 55.00 1776 9 456.50 396.80 313.64 130.53 222.49 122.44 129.60 47.00 50.00 55.00 1924 10 456.53 460.00 304.39 180.53 222.61 160.00 129.62 51.26 52.06 55.00 2072 11 456.50 460.00 298.47 230.53 222.61 160.00 129.60 81.26 52.05 55.00 2146 12 456.50 460.00 340.00 241.47 240.28 160.00 129.59 85.31 52.04 55.00 2220 13 456.77 396.80 310.11 241.35 222.43 122.51 129.67 85.31 52.04 55.00 2072 14 379.80 316.80 322.18 241.26 219.70 122.03 129.86 85.32 52.06 55.00 1924 15 303.07 302.75 297.33 240.83 222.60 123.97 123.09 85.30 22.06 55.00 1776 16 227.02 222.75 313.81 241.44 172.73 122.84 93.09 85.31 20.00 55.00 1554 17 150.00 228.95 291.48 241.12 222.61 122.46 63.09 85.28 20.00 55.00 1480 18 226.62 308.95 289.47 240.97 222.64 122.45 56.58 85.30 20.00 55.00 1628 19 303.24 388.95 298.32 241.24 172.70 122.52 86.58 85.39 22.05 55.00 1776 20 379.85 460.00 310.58 291.24 222.61 122.30 93.04 85.31 52.06 55.00 2072 21 303.39 382.53 274.17 299.97 222.69 125.89 123.04 85.31 52.01 55.00 1924 22 226.62 302.53 194.97 249.98 172.73 160.00 129.59 85.31 52.06 55.00 1628 23 150.0 222.53 185.52 206.67 122.87 122.45 129.59 85.31 52.06 55.00 1332 24 150.00 142.53 191.64 180.84 73.00 123.96 129.65 85.31 52.06 55.00 1184 TABLE 4. COMPARISON OF TOTAL PRODUCTION COST AND SIMULATION TIME OBTAINED BY CONSTRAINED EP Methods Seq. Quad. Prog [20] Evol. Prog [20] Constraned EP Total Cost ($) 1051163 1048640 1026500 CPU Tme 1.19 48 45.15

Proceedngs of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lsbon, Portugal, June 16-18, 2005 (pp160-166) VII APPENDIX TABLE A.1 LOAD PATTERN FOR 12-HOUR SYSTEM [19] Hours 1 2 3 4 5 6 7 8 9 10 11 12 Load(MW) 5560 5620 5800 5810 5990 6041 6001 5790 5680 5540 5690 5750 TABLE A.2 DATA FOR 10-UNIT SYSTEM (SMOOTH FUEL COST FUNCTION) [19] UNITS P mn P max a b a UR DR (MW) (MW) ($/MW 2 h) ($/MWh) ($/h) (MW/h) (MW/h) 1 155 360 0.03720 26.4408 180 20 25 2 320 680 0.03256 21.0771 275 20 25 3 323 718 0.03102 18.6626 352 50 50 4 275 680 0.02871 16.8894 792 50 50 5 230 600 0.03223 17.3998 440 50 50 6 350 748 0.02064 21.6180 348 50 50 7 220 620 0.02268 15.1716 588 100 100 8 225 643 0.01776 14.5632 984 100 150 9 350 920 0.01644 14.3448 1260 100 150 10 450 1050 0.01620 13.5420 1260 100 150 TABLE A.3 LOAD PATTERN FOR 24-HOUR SYSTEM [20] Hour 1 2 3 4 5 6 7 8 9 10 11 12 Load(MW) 1036 1110 1258 1406 1480 1628 1702 1776 1924 2072 2146 2220 Hour 13 14 15 16 17 18 19 20 21 22 23 24 Load(MW) 2072 1924 1776 1554 1480 1628 1776 2072 1924 1628 1332 1184 TABLE A.4 DATA FOR 10-UNIT SYSTEM (NON-SMOOTH FUEL COST FUNCTION) [20] UNITS P mn P max a b c e F UR DR (MW) (MW) ($/MW 2 h) ($/MWh) ($/h) ($/h) (rad/mw) (MW/h) (MW/h) 1 150 470 0.00043 21.60 958.20 450 0.041 80 80 2 135 460 0.00063 21.05 1313.6 600 0.036 80 80 3 73 340 0.00039 20.81 604.97 320 0.028 80 80 4 60 300 0.00070 23.90 471.60 260 0.052 50 50 5 73 243 0.00079 21.62 480.29 280 0.063 50 50 6 57 160 0.00056 17.87 601.75 310 0.048 50 50 7 20 130 0.00211 16.51 502.70 300 0.086 30 30 8 47 120 0.00480 23.23 639.40 340 0.082 30 30 9 20 80 0.10908 19.58 455.60 270 0.098 30 30 10 55 55 0.00951 22.54 692.40 380 0.094 30 30 Notatons Used n Tables A.2 and A.4: UR= Up Ramp Rate; DR= Down Ramp Rate.