CHAPTER - 7. Firefly Algorithm based Strategic Bidding to Maximize Profit of IPPs in Competitive Electricity Market

Similar documents
International Journal of Pure and Applied Sciences and Technology

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

6 Roots of Equations: Open Methods

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

Introduction to Numerical Integration Part II

Applied Statistics Qualifier Examination

4. Eccentric axial loading, cross-section core

CALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVEYS

Dynamic Power Management in a Mobile Multimedia System with Guaranteed Quality-of-Service

Multiple view geometry

Lecture 4: Piecewise Cubic Interpolation

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

Zbus 1.0 Introduction The Zbus is the inverse of the Ybus, i.e., (1) Since we know that

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

GAUSS ELIMINATION. Consider the following system of algebraic linear equations

CIS587 - Artificial Intelligence. Uncertainty CIS587 - AI. KB for medical diagnosis. Example.

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus:

Online Appendix to. Mandating Behavioral Conformity in Social Groups with Conformist Members

APPLICATION OF MULTI OBJECTIVE FUZZY LINEAR PROGRAMMING IN SUPPLY PRODUCTION PLANNING PROBLEM

Numerical Solution of Fredholm Integral Equations of the Second Kind by using 2-Point Explicit Group Successive Over-Relaxation Iterative Method

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

Linear and Nonlinear Optimization

Remember: Project Proposals are due April 11.

The Schur-Cohn Algorithm

8. INVERSE Z-TRANSFORM

Lecture 36. Finite Element Methods

6. Chemical Potential and the Grand Partition Function

Quiz: Experimental Physics Lab-I

523 P a g e. is measured through p. should be slower for lesser values of p and faster for greater values of p. If we set p*

Modeling Labor Supply through Duality and the Slutsky Equation

Bi-level models for OD matrix estimation

A New Algorithm Linear Programming

NUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION

Genetic Programming. Outline. Evolutionary Strategies. Evolutionary strategies Genetic programming Summary

ICS 252 Introduction to Computer Design

Using Predictions in Online Optimization: Looking Forward with an Eye on the Past

7.2 Volume. A cross section is the shape we get when cutting straight through an object.

Abhilasha Classes Class- XII Date: SOLUTION (Chap - 9,10,12) MM 50 Mob no

Sequences of Intuitionistic Fuzzy Soft G-Modules

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

Demand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert

Advanced Machine Learning. An Ising model on 2-D image

Variable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia

Frequency scaling simulation of Chua s circuit by automatic determination and control of step-size

Study of Trapezoidal Fuzzy Linear System of Equations S. M. Bargir 1, *, M. S. Bapat 2, J. D. Yadav 3 1

Course Review Introduction to Computer Methods

The Study of Lawson Criterion in Fusion Systems for the

Altitude Estimation for 3-D Tracking with Two 2-D Radars

An Introduction to Support Vector Machines

Jean Fernand Nguema LAMETA UFR Sciences Economiques Montpellier. Abstract

LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB

3/6/00. Reading Assignments. Outline. Hidden Markov Models: Explanation and Model Learning

Online Learning Algorithms for Stochastic Water-Filling

Support vector machines for regression

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

Polynomial Regression Models

Haddow s Experiment:

Effects of polarization on the reflected wave

Least squares. Václav Hlaváč. Czech Technical University in Prague

Systematic Construction of examples for cycling in the simplex method

Wars of attrition and all-pay auctions with stochastic competition

Coordinating Multi-Attribute Procurement Auctions in Finite Capacity Assembly Environments

A Family of Multivariate Abel Series Distributions. of Order k

18.7 Artificial Neural Networks

4.1. Probability Density Functions

ANALOG CIRCUIT SIMULATION BY STATE VARIABLE METHOD

The Minimum Label Spanning Tree Problem: Illustrating the Utility of Genetic Algorithms

Multi-Stage Power Distribution Planning to Accommodate High Wind Generation Capacity

Dynamic Power Management in a Mobile Multimedia System with Guaranteed Quality-of-Service

Mixed Type Duality for Multiobjective Variational Problems

Probabilistic Graphical Models

Exploiting Structure in Probability Distributions Irit Gat-Viks

SAM UPDATING USING MULTIOBJECTIVE OPTIMIZATION TECHNIQUES ABSTRACT

Acceptance Sampling by Attributes

Name: SID: Discussion Session:

Strong Gravity and the BKL Conjecture

Electrochemical Thermodynamics. Interfaces and Energy Conversion

1B40 Practical Skills

1 Online Learning and Regret Minimization

INTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

An Ising model on 2-D image

Statistical Timing Analysis for Intra-Die Process Variations with Spatial Correlations

Trade-offs in Optimization of GMDH-Type Neural Networks for Modelling of A Complex Process

The Study of Teaching-learning-based Optimization Algorithm

Chapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses

5.1 How do we Measure Distance Traveled given Velocity? Student Notes

Credit Card Pricing and Impact of Adverse Selection

Section 6: Area, Volume, and Average Value

Chemical Reaction Engineering

DESIGN OF MULTILOOP CONTROLLER FOR THREE TANK PROCESS USING CDM TECHNIQUES

6.6 The Marquardt Algorithm

Forecast of Next Day Clearing Price in Deregulated Electricity Market

Solution of Tutorial 5 Drive dynamics & control

A New Markov Chain Based Acceptance Sampling Policy via the Minimum Angle Method

Chemical Reaction Engineering

Transcription:

CHAPTER - 7 Frefly Algorthm sed Strtegc Bddng to Mxmze Proft of IPPs n Compettve Electrcty Mrket 7. Introducton The renovton of electrc power systems plys mjor role on economc nd relle operton of power system. The generton compnes nd tl end customers ppers to undergo mjor multple tsk of desgnng of proper opertng methodologes. A sgnfcnt mount of reserch work s n vogue n the pst pertnng to the mrket structure nd to the development of optml ddng strteges. The tools supportng the ddng process seem to wy from glol soluton n vew of the complexty nd szes of the prctcl prolems. Therefore n exhustve formulton of optml ddng strtegy forys contenton of the generton compnes nd the end consumers. An nnovtve pproch for defnng optml ddng strtegy for Independent Power Producers sngle sde cton s presented s stochstc optmzton prolem nd solved y Frefly lgorthm.fa. The Frefly Algorthm s Met heurstc, nture nspred, optmzton lgorthm whch s sed on the socl flshng ehvor of frefles nd hs een ntroduced for the ddng prolem to otn the glol optml soluton. It effectvely mxmzes the proft of power supplers. The numercl exmples wth Sx genertors power supplers test system nd IEEE 30 us system re consdered to revel the essentl fetures of the proposed method nd test results tulted. The smulton result shows tht ths pproch effectvely 07

mxmzes the Proft of Power supplers, converge much fster nd more relle when compred wth the exstng methods [3]. 7. Prolem formulton 7.. Mthemtcl Model Consder system consst of m power supplers prtcptng n pool-sed sngle-uyer electrcty mrket n whch the seled ucton wth unform Mrket Clerng Prce MCP s employed nd pctorlly represented n Fg, 7.. Assume tht ech suppler s requred to d lner supply functon to the pool. Fg. 7.. A Typcl Model of Electrcty Mrket for Sngle sde cton Fg. 7.. Mrket Equlrum Pont for MCP 08

The th suppler d wth lner supply curve s denoted y G P = + P where =,..m. The P s the ctve power, nd re non-negtve ddng coeffcents of the th suppler. The Independent System Opertor ISO receves d from ll mrket prtcpnts nd usng predcted ggregte lod from the smll users, the ISO determnes MCP tht ttempts to lnce the energy demnd nd Supply. The process s grphclly expressed n Fg. 7.. The ojectve of IPPs s to mxmze ther own proft. Suppose the power producers expressed y cost functon n equton 7. P = ep fp C + 7. The ojectve functon of power producer cn e defned s n equton 7. Mx : F, = RPC P 7. Mrket Clerng Prce R s represented y the equton 7.3 R Q + 0 = = m K+ m = 7.3 The ggregted lod demnd cn e formulted s n equton 7.4 Q R = Q KR 7.4 o Constrnts. Power lnce constrnts: m = P = Q R 7.5 p R = =,... m 7.6 09

0. Power generton lmt constrnts: mx mn p p p =,... m 7.7 7... Development of ddng strtegy The GENCOs my not hve the posslty to know the exct nformton of ther compettor. Hence t s mndtory for GENCO to understnd the opponents unknown nformton. It s ssumed tht the prevous ddng coeffcents re vlle n pulc domn. The th GENCO cn therefore estmte the ddng prmeters usng prolty densty functon pdf otned through sttstcl nlyss of hstorcl ddng dt. Usully, pdf s used to represent the dstruton of rvls ddng prmeters whch cn e expressed s n equton 7.8 = exp x x pdf µ π 7.8 Where, - Stndrd devton µ - Men vlues The ddng coeffcents, re dependent of ech other. Hence one of the co-effcent s kept constnt nd the other s rtrrly chnced usngpdf. The pdf of rndom vrle s formulton whch cn e ntegrted so s to get the prolty tht the vrle choose vlue s specfed ntervl, + Π = exp, pdf µ µ µ ρ µ ρ ρ 7.9

Where ρ s the correlton co-effcent etween nd.the men µ, µ nd stndrd devton, re the prmeters of the jont dstruton. The mrgnl dstruton of, re norml wth men vlues µ, µ nd stndrd devtons, respectvely. Bsed on hstorcl ddng dt these dstrutons cn e determned. The prolty densty functon equton 7.9 represents the jont dstrutons etween nd, the tsk of optmlly coordntng the ddng strteges for suppler wth ojectve functon 7. nd constrnts 7.5 to 7.7, ecomes stochstc optmzton prolem. The Frefly lgorthm s ppled to solve the ove stochstc optmzton prolem. 7.3 Implementton of Frefly lgorthm to solve Bddng prolem Wth vew to sell electrcty t optml prces nd to mxmze proft, the power producers nd consumers need exclusve ddng strteges tht must consder constrnts such s Power lnce, Genertor lmts nd Lod consumpton lmts of mrket prtcpnts. The Frefly Algorthm cn drectly solve optml ddng prolem Mxmze proft ecuse of ts mxmzton chrcterstcs nd the flow chrt of the method s shown n Fg. 7.3. The Frefly Algorthm ncludes four essentl prmeters, Populton sze n, Attrctveness β, rndomzton prmeter α nd Asorpton coeffcent γ.the fesle prmeters otned y tertve processes re s follows. α = 0. 0.9, β = 0..0, γ = 0. 0 nd n = 5 50. Therefore, the followng prmeters of the proposed FA re consdered to solve the optml ddng prolem of sx ndependent power producers nd two lrge consumers. Where n = 30, β =0.0, = 0.5, α γ = nd the mxmum numer of tertons = 5000. Owng to the rndom nture of the FA, ther performnce cnnot e judged y the result of sngle run. Mny trls wth ndependent populton ntlztons re

necesstted to e mde to otn useful concluson of the performnce of the pproch. Strt Red system dt, Genertor dt Cost coeffcents, Genertor lmts, Aggregted lod nd Prce Elstcty Intlze the FA prmeters: Populton sze n. Attrctveness β, rndomzton prmeter α, Asorpton coeffcent γ nd numer of tertons Crete the ntl rndom populton of ddng co effcent Usng ddng coeffcents clculte Mrket Clerng Prce R m Q0 + = R = m K + = From MCP clculte ftness of ech populton usng the equton Mxmze: F, = RP C P Apply FA prmeters to otn the optml soluton Mx proft No Whether optml soluton s reched Yes Prnt the profts of power supplers Stop Fg.7.3. Flow chrt for Proft Mxmzton of IPPs y Proposed method

The supremcy of the proposed FA, s rought out through the test results nd vldted those reported n the recently pulshed methods such s PSO, GA nd Conventonl method for solvng the ddng prolem. The scenros re progrmmed n MATLAB 9.0 nd smulton crred on computer wth Pentum IV, Intel Dul core. GHz, GB RAM. 7.4 Cse study nd Results The GENCOs uld optml ddng strteges to mxmze the proft of power supplers nd mplements the modelng. The smulton s crred out on two cses of test system. The frst system conssts of Sx power supplers s Cse- I nd second test system of IEEE 30 us s Cse-. Test cse: Sx power supplers Test System The proposed Frefly pproch s ppled to test system gven n [] whch conssts of sx power supplers. The cost coeffcents of power generton nd mxmum/ mnmum lmts of sx power supplers re mentoned n ppendx A.9 Tle A.9.. Tle 7. Smulton Results for Sx Power supplers GENCOs Bddng Strtegy $/MW Bddng PowerMW Proft 0.038 58.664 38.935 0.0490 50.0000 78.7500 3 0.0837 34.30 76.6493 4 0.047 7.3485 66.466 5 0.0463 5.8886 6.0396 6 0.0463 5.8886 6.0376 Totl proft $ 5.866 Mrket clerng prce $/h 4.0 Computtonl tme sec.98 $ 3

The fuel cost functon of ech genertor s expressed s qudrtc equton. The prmeters ssocted wth the lod chrcterstcs re consdered from the sme reference where n the ggregted lodq 0 s equl to 450 MW nd the prce elstcty K equls to 0. Tle 7. Comprson of MCP nd Generted power of Proposed wth Conventonl method GENCOs FA Proposed Conventonl Method [] MCP PMW MCP PMW 58.664 56.0 50.0000 6.78 3 4.000 34.30 4.0386 37.95 4 7.3485 64.84 5 5.8886 4.47 6 5.8886 4.47 The smulton results of power supplers re presented n Tle 7.. It ncludes ddng strtegy, ddng power, MCP nd proft of power supplers. The comprtve studes wth conventonl method [] re dsplyed n Tle 7. to nlyze the MCP nd ddng power of power supplers. It cn e seen tht the proposed method provde mxmzed profts whch s etter thn the conventonl method. Besdes, t converges much fster nd more relle thn the other vlle methods. Test cse: IEEE 30-us system The IEEE 30 us system conssts of Sx Power supplers who supply electrcty to n ggregte lod. The genertor dt tken from reference [] s shown n ppendx A.0 Tle A.0.. The fuel cost functon of ech genertor s vlle s qudrtc equton. The prmeters ssocted wth the lod chrcterstcs re consdered from the sme reference where n the ggregted lodq 0 s equl to 500 MW nd prce elstctykequls to 0. 4

Tle 7.3 Smulton Results of Sx Power supplers for IEEE 30 us system GENCOs Bddng Strtegy $/MW Bddng Power MW Revenue $ Fuel Cost $ Proft 0.0353 60.00 4.00 46.00 808.00 0.07334 80.44 65.37 54.00 36.36 3 0.883 48.39 370.8 94.74 75.44 4 0.04435 99. 758.96 404.5 354.44 5 0.06700 55.97 48.7 46.3 8.94 6 0.06700 55.97 48.7 46.3 8.94 Totl Proft $ 063. Mrket clerng prce $/MWh 7.65 $ The smulton results of Sx power supplers for IEEE 30 us system s presented n Tle 7.3. It ncludes ddng strtegy, ddng power, MCP, revenue, fuel cost nd proft of the sx power supplers. The totl proft of sx power suppler s equls to 063. $ nd computtonl tme.98 sec. It s due to fct tht the Frefly lgorthm plys vtl role n serch of the glol optml soluton. Tle 7.4 Comprson of Bddng Strteges of Sx Power supplers for IEEE 30 us system GENCOs FA Proposed PSO [7] GA [4] Trdtonl GSS method [04] 0.0353 0.0009 0.00045 0.5800 0.07334 0.050953 0.048786 0.04745 3 0.883 0.8976 0.7454 0.3099 4 0.04435 0.0483 0.0350 0.0458 5 0.06700 0.0779 0.069694 0.0564 6 0.06700 0.0779 0.069694 0.0564 5

GENCOs Tle 7.5 Comprson of Bddng Power nd Proft of Sx Power supplers for IEEE 30 us system FA Proposed PSO [7] GA [4] Trdtonl GSS method [04] PMW Proft$ PMW Proft$ PMW Proft$ PMW Proft$ 60.00 808.00 60.00 77.4 60.00 74.45 60.00 557.00 80.44 36.36 00.83 340.0 0.0 3.3 9.30 49.00 3 48.39 75.44 3.35 5.06 3.68 9.33 38.80 03.00 4 99. 354.44 00.00 80.36 00.00 6.0 00.00 00.00 5 55.97 8.94 53.40 36.3 53.00 5.56 54.90 94.00 6 55.97 8.94 53.40 36.3 53.00 5.56 54.90 94.00 Tle 7.6 Comprson of MCP, Totl proft nd Computtonl tme of Sx Power supplers for IEEE 30 us system FA Proposed PSO [7] GA [4] Trdtonl GSS method [04] MCP $/hr 7.65 6.88 6.69 6.08 Totl Profts $ Computtonl tme sec 063. 790.57 694.3 97.00.98.06 6.4 -- Tle 7.7 Performnce Comprson of Proposed method wth Exstng methods for IEEE 30 us system Totl proft $ Methods Best Averge Worst FA Proposed 063. 047.5 0.67 PSO [7] 75.60 689.97 67.34 GA [4] 83.89 87.97 09.06 6

The comprtve studes wth Prtcle Swrm Optmzton [7], Genetc Algorthm [4] nd Trdtonl GSS method [04] re mde to nlyze the ddng coeffcents of power supplers nd dsplyed n Tle 7.4. The Tle 7.5 elortes the Bddng power nd Proft of dfferent methods. The comprson of mrket clerng prce MCP, totl profts nd computtonl tme of power supplers for dfferent methods re presented n Tle 7.6. The performnce of totl profts of power supplers re compred wth proposed nd exstng methods n Tle 7.7. It s evdent from Tle 7.6 the totl proft of the proposed method s mproved wth less computtonl tme thn the other vlle methods. 7.5 Summry The Frefly lgorthm hs een ppled to solve ddng strtegy n order to mprove the proft of GENCOs Power supplers n n open electrcty mrket. The numercl exmples wth Sx genertors power supplers Test system nd IEEE 30 us system hve een consdered to llustrte the essentl fetures of the proposed method. The Frefly lgorthm hs een used to determne the optml ddng strtegy n dfferent mrket rule, dfferent fxed lod, dfferent cpcty of uyers nd sellers. The results hve een projected to rng out the promsng nture of technque for solvng complcted power system optmzton prolem under deregulted envronment. 7