ECG 740 GENERATION SCHEDULING (UNIT COMMITMENT)

Similar documents
Contents Economic dispatch of thermal units

Solving the Unit Commitment Problem Using Fuzzy Logic

Sample Problems for Exam 1, EE Note:

Tight and Compact MILP Formulation for the Thermal Unit Commitment Problem

3 Unit Commitment. H, = P P5 MBtu/h 5.1 INTRODUCTION EXAMPLE 5A. Unit 1: Unit 3:

VOL. 2, NO. 11, Dec 2012 ISSN ARPN Journal of Science and Technology All rights reserved.

Economic Operation of Power Systems

Multi-Area Stochastic Unit Commitment for High Wind Penetration in a Transmission Constrained Network

Risk-based Security Constrained Unit Commitment

Power Engineering II. Fundamental terms and definitions

A Simplified Lagrangian Method for the Solution of Non-linear Programming Problem

Systems Operations. PRAMOD JAIN, Ph.D. Consultant, USAID Power the Future. Astana, September, /6/2018

Thermal Unit Commitment Problem

Thermal Unit Commitment

Draft Wholesale Power Price Forecasts

04-Economic Dispatch 2. EE570 Energy Utilization & Conservation Professor Henry Louie

03-Economic Dispatch 1. EE570 Energy Utilization & Conservation Professor Henry Louie

Chapter 2. Planning Criteria. Turaj Amraee. Fall 2012 K.N.Toosi University of Technology

Chapter 2 Deterministic Unit Commitment Models and Algorithms

A Unified Framework for Defining and Measuring Flexibility in Power System

Optimization Basics for Electric Power Markets

Application of Monte Carlo Simulation to Multi-Area Reliability Calculations. The NARP Model

Automatic Generation Control. Meth Bandara and Hassan Oukacha

In the original knapsack problem, the value of the contents of the knapsack is maximized subject to a single capacity constraint, for example weight.

Accelerating the Convergence of MIP-based Unit Commitment Problems

OPTIMAL DISPATCH OF REAL POWER GENERATION USING PARTICLE SWARM OPTIMIZATION: A CASE STUDY OF EGBIN THERMAL STATION

California Independent System Operator (CAISO) Challenges and Solutions

On Optimal Power Flow

Course notes for EE394V Restructured Electricity Markets: Locational Marginal Pricing

Stochastic Dual Dynamic Programming with CVaR Risk Constraints Applied to Hydrothermal Scheduling. ICSP 2013 Bergamo, July 8-12, 2012

Sensitivity Analysis of a Mixed Integer Linear Programming Model For Optimal Hydrothermal Energy Generation For Ghana

International Studies about the Grid Integration of Wind Generation

Multi-Area Stochastic Unit Commitment for High Wind Penetration

Stochastic Unit Commitment with Topology Control Recourse for Renewables Integration

5 Integer Linear Programming (ILP) E. Amaldi Foundations of Operations Research Politecnico di Milano 1

Information Document Calculation of Pool Price and Transmission Constraint Rebalancing Costs During a Constraint Event ID # R

MEDIUM-TERM HYDROPOWER SCHEDULING BY STOCHASTIC DUAL DYNAMIC INTEGER PROGRAMMING: A CASE STUDY

Part 4. Decomposition Algorithms

Controlling variability in power systems

Power System Seminar Presentation Wind Forecasting and Dispatch 7 th July, Wind Power Forecasting tools and methodologies

CE 191: Civil & Environmental Engineering Systems Analysis. LEC 17 : Final Review

APPENDIX 7.4 Capacity Value of Wind Resources

Javier Contreras Sanz- Universidad de Castilla-La Mancha Jesús María López Lezama- Universidad de Antioquia Antonio Padilha-Feltrin- Universidade

Total Market Demand Wed Jan 02 Thu Jan 03 Fri Jan 04 Sat Jan 05 Sun Jan 06 Mon Jan 07 Tue Jan 08

UNIT-I ECONOMIC OPERATION OF POWER SYSTEM-1

Lagrangian Relaxation and Tabu Search Approaches for the Unit Commitment Problem

THE STOCHASTIC UNIT COMMITMENT PROBLEM: A CHANCE CONSTRAINED PROGRAMMING APPROACH CONSIDERING EXTREME MULTIVARIATE TAIL PROBABILITIES

A Progressive Hedging Approach to Multistage Stochastic Generation and Transmission Investment Planning

526 IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 27, NO. 1, FEBRUARY 2012

Introduction to Bin Packing Problems

LINEAR PROGRAMMING APPROACH FOR THE TRANSITION FROM MARKET-GENERATED HOURLY ENERGY PROGRAMS TO FEASIBLE POWER GENERATION SCHEDULES

Benders' Method Paul A Jensen

Value of Forecasts in Unit Commitment Problems

I. Why power? II. Model definition III. Model uses IV. The basics. A. Supply Dispatch B. Demand Bidding. D. Transmission

SUBMITTED TO IEEE TRANSACTIONS ON POWER SYSTEMS, AUGUST

Integer Programming ISE 418. Lecture 8. Dr. Ted Ralphs

Optimal Operation of Large Power System by GA Method

A Robust Optimization Approach to the Self-scheduling Problem Using Semidefinite Programming

On-line supplement to: SMART: A Stochastic Multiscale Model for the Analysis of Energy Resources, Technology

R O B U S T E N E R G Y M AN AG E M E N T S Y S T E M F O R I S O L AT E D M I C R O G R I D S

IE418 Integer Programming

Energy System Modelling Summer Semester 2018, Lecture 7

2001 Dennis L. Bricker Dept. of Industrial Engineering The University of Iowa. Reducing dimensionality of DP page 1

Decomposition Techniques in Mathematical Programming

IE 400 Principles of Engineering Management. Graphical Solution of 2-variable LP Problems

Theory of constraints and linear programming: a reexamination

Discrete lot sizing and scheduling on parallel machines: description of a column generation approach

BRANCH-AND-BOUND SCHEDULING FOR THERMAL GENERATING UNITS

IBM Research Report. Stochasic Unit Committment Problem. Julio Goez Lehigh University. James Luedtke University of Wisconsin

Appendix A Solving Systems of Nonlinear Equations

Influence of knn-based Load Forecasting Errors on Optimal Energy Production

Mathematical Foundations -1- Constrained Optimization. Constrained Optimization. An intuitive approach 2. First Order Conditions (FOC) 7

ECE 422/522 Power System Operations & Planning/ Power Systems Analysis II 4 Active Power and Frequency Control

Mixed Integer Linear Programming Formulation for Chance Constrained Mathematical Programs with Equilibrium Constraints

IV. Violations of Linear Programming Assumptions

Incorporating Demand Response with Load Shifting into Stochastic Unit Commitment

System Planning Lecture 7, F7: Optimization

Analysis of Coupling Dynamics for Power Systems with Iterative Discrete Decision Making Architectures

Integer Linear Programming Modeling

1 Column Generation and the Cutting Stock Problem

A three-level MILP model for generation and transmission expansion planning

20.1 Xenon Production Xe-135 is produced directly in only 0.3% of all U-235 fissions. The following example is typical:

An energy management system for off-grid power systems

Nonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control

SIGNIFICANT increase in amount of private distributed

Analysis of Adaptive Certainty-Equivalent Techniques for the Stochastic Unit Commitment Problem

Bringing Renewables to the Grid. John Dumas Director Wholesale Market Operations ERCOT

Multistage Adaptive Robust Optimization for the Unit Commitment Problem

2. Linear Programming Problem

On Clearing Coupled Day-Ahead Electricity Markets

Valid Inequalities for Optimal Transmission Switching

Tight MIP Formulations of the Power-Based Unit Commitment Problem

Electrical Power and Energy Systems

Research Article A Hybrid Search Algorithm for Midterm Optimal Scheduling of Thermal Power Plants

Congestion Management in a Smart Grid via Shadow Prices

RO: Exercices Mixed Integer Programming

Energy Forecasting Customers: Analysing end users requirements Dec 3rd, 2013 Carlos Alberto Castaño, PhD Head of R&D

Chapter 5. Transmission networks and electricity markets

Perfect and Imperfect Competition in Electricity Markets

Capacity Scarcity Condition Monday, September 3, 2018 Two primary factors led to the implementation of OP 4 event Significant generation outages and r

Transcription:

1 ECG 740 GENERATION SCHEDULING (UNIT COMMITMENT)

2 Unit Commitment Given a load profile, e.g., values of the load for each hour of a day. Given set of units available, When should each unit be started, stopped and how much should it generate to meet the load at minimum cost? Load 1000 500 0 6 12 18 24 Time??? G G G Load Profile

3 A Simple Example Unit 1: P Min = 250 MW, P Max = 600 MW C 1 = 510.0 + 7.9 P 1 + 0.00172 P 1 2 $/h Unit 2: P Min = 200 MW, P Max = 400 MW C 2 = 310.0 + 7.85 P 2 + 0.00194 P 2 2 $/h Unit 3: P Min = 150 MW, P Max = 500 MW C 3 = 78.0 + 9.56 P 3 + 0.00694 P 3 2 $/h What combination of units 1, 2 and 3 will supply a particular load at minimum cost? How much should each unit in that combination generate?

Optimal combination for each hour 4

5 Matching the combinations to the load Load Unit 3 Unit 2 Unit 1 0 6 12 18 24 Time

Issues that must be considered Constraints Unit constraints System constraints Some constraints create a link between periods Start-up costs Cost incurred when we start a generating unit Different units have different start-up costs Curse of dimensionality 6

Unit Constraints Constraints that affect each unit individually: Maximum generating capacity Minimum stable generation Minimum up time Minimum down time Ramp rate 7

8 Notations u(i,t) : u(i,t) = 1: u(i,t) = 0 : x(i,t) : Status of unit i at period t Unit i is on during period t Unit i is off during period t Power produced by unit i during period t

9 Minimum up- and down-time Minimum up time Once a unit is running it may not be shut down immediately: If u(i,t) = 1 and t i up < t i up,min then u(i,t +1) = 1 Minimum down time Once a unit is shut down, it may not be started immediately If u(i,t) = 0 and t i down < t i down,min then u(i,t +1) = 0

Ramp rates Maximum ramp rates To avoid damaging the turbine, the electrical output of a unit cannot change by more than a certain amount over a period of time: Maximum ramp up rate constraint: up,max x( i,t +1) - x( i,t) DP i Maximum ramp down rate constraint: x(i,t) - x(i,t +1) DP i down,max 10

System Constraints Constraints that affect more than one unit Load/generation balance Reserve generation capacity Emission constraints Network constraints 11

12 Load/Generation Balance Constraint N å i=1 u(i,t)x(i,t) = L(t) N : Set of available units

13 Reserve Capacity Constraint Unanticipated loss of a generating unit or an interconnection causes unacceptable frequency drop if not corrected rapidly Need to increase production from other units to keep frequency drop within acceptable limits Rapid increase in production only possible if committed units are not all operating at their maximum capacity N å i=1 u(i,t) P max i ³ L(t) + R(t) R(t): Reserve requirement at time t

Cost of Reserve Reserve has a cost even when it is not called More units scheduled than required Units not operated at their maximum efficiency Extra start up costs Must build units capable of rapid response Cost of reserve proportionally larger in small systems Important driver for the creation of interconnections between systems 14

Environmental constraints Scheduling of generating units may be affected by environmental constraints Constraints on pollutants such SO 2, NO x Various forms: Limit on each plant at each hour Limit on plant over a year Limit on a group of plants over a year Constraints on hydro generation Protection of wildlife Navigation, recreation 15

Network Constraints Transmission network may have an effect on the commitment of units Some units must run to provide voltage support The output of some units may be limited because their output would exceed the transmission capacity of the network A B Cheap generators May be constrained off More expensive generator May be constrained on 16

Start-up Costs Thermal units must be warmed up before they can be brought on-line Warming up a unit costs money Start-up cost depends on time unit has been off α i + β i t SC i (t OFF - i OFF i ) = a i + b i (1 - e t i ) α i t i OFF 17

18 Start-up Costs Need to balance start-up costs and running costs Example: Diesel generator: low start-up cost, high running cost Coal plant: high start-up cost, low running cost Issues: How long should a unit run to recover its start-up cost? Start-up one more large unit or a diesel generator to cover the peak? Shutdown one more unit at night or run several units partloaded?

19 Summary Some constraints link periods together Minimizing the total cost (start-up + running) must be done over the whole period of study Generation scheduling or unit commitment is a more general problem than economic dispatch Economic dispatch is a sub-problem of generation scheduling

20 Solving the Unit Commitment Problem Decision variables: Status of each unit at each period: u(i,t) Î{ 0,1} " i,t Output of each unit at each period: { } " i,t x(i,t) Î 0, P min max é ë i ;P i ù û Combination of integer and continuous variables

21 Optimization with integer variables Continuous variables Can follow the gradients or use LP Any value within the feasible set is OK Discrete variables There is no gradient Can only take a finite number of values Problem is not convex Must try combinations of discrete values

How many combinations are there? 111 110 101 Examples 3 units: 8 possible states N units: 2 N possible states 100 011 010 001 000 22

How many solutions are there? Optimization over a time horizon divided into intervals A solution is a path linking one combination at each interval How many such paths are there? T= 1 2 3 4 5 6 2011 Daniel Kirschen and the University of Washington 23

24 How many solutions are there? ( 2 N )( 2 N ) ( 2 N ) = ( 2 N ) T T= 1 2 3 4 5 6

25 The Curse of Dimensionality Example: 5 units, 24 hours ( 2 N ) T = ( 2 5 ) 24 = 6.2 10 35 combinations Processing 10 9 combinations/second, this would take 1.9 10 19 years to solve There are dozens of units in large power systems... Luckily, many of these combinations do not satisfy the constraints

How do you Beat the Curse? Brute force approach won t work! Need to be smart Try only a small subset of all combinations Can t guarantee optimality of the solution Try to get as close as possible within a reasonable amount of time 26

27 Main Solution Techniques Characteristics of a good technique Solution close to the optimum Reasonable computing time Ability to model constraints Priority list / heuristic approach Dynamic programming Lagrangian relaxation Mixed Integer Programming

Simple Unit Commitment Problem, N = 3 Unit P min (MW) P max (MW) Min up (h) Min down (h) No-load cost ($) Marginal cost ($/MWh) Start-up cost ($) Initial status A 150 250 3 3 0 10 1,000 ON B 50 100 2 1 0 12 600 OFF C 10 50 1 1 0 20 100 OFF 2011 Daniel Kirschen and the University of Washington 28

29 Demand Data Hourly Demand Load 350 300 250 200 150 100 50 0 1 2 3 Hours Reserve requirements are not considered

30 Feasible Unit Combinations (states) Combinations A B C P min P max 1 1 1 210 400 1 1 0 200 350 1 0 1 160 300 1 0 0 150 250 0 1 1 60 150 0 1 0 50 100 0 0 1 10 50 0 0 0 0 0 1 2 3 150 300 200

31 Transitions between feasible combinations A B C 1 1 1 1 1 0 1 0 1 1 0 0 0 1 1 Initial State 1 2 3

32 Infeasible transitions: Minimum down time of unit A A B C 1 1 1 1 1 0 1 0 1 1 0 0 0 1 1 Initial State T D T U A 3 3 B 1 2 C 1 1 1 2 3

33 Infeasible transitions: Minimum up time of unit B A B C 1 1 1 1 1 0 1 0 1 1 0 0 0 1 1 Initial State T D T U A 3 3 B 1 2 C 1 1 1 2 3

34 Feasible transitions A B C 1 1 1 1 1 0 1 0 1 1 0 0 0 1 1 Initial State 1 2 3

35 Operating costs 1 1 1 4 1 1 0 3 7 1 0 1 1 0 0 1 2 6 5

Economic dispatch State Load P A P B P C Cost 1 150 150 0 0 1500 2 300 250 0 50 3500 3 300 250 50 0 3100 4 300 240 50 10 3200 5 200 200 0 0 2000 6 200 190 0 10 2100 7 200 150 50 0 2100 Unit P min P max No-load cost Marginal cost A 150 250 0 10 B 50 100 0 12 C 10 50 0 20 36

37 Operating costs 1 1 1 4 $3200 1 1 0 1 0 1 1 0 0 1 $1500 3 $3100 2 $3500 7 $2100 6 $2100 5 $2000

Start-up costs 1 1 1 1 1 0 1 0 1 $700 $600 $100 4 $3200 3 $3100 2 $3500 $0 $0 $600 $0 $0 7 $2100 6 $2100 1 0 0 $0 1 $1500 Unit Start-up cost 5 $2000 A 1000 B 600 C 100 38

39 Accumulated costs 1 1 1 1 1 0 1 0 1 $1500 1 0 0 $0 1 $1500 $700 $600 $100 $5400 4 $3200 $5200 3 $3100 $5100 2 $3500 $0 $0 $600 $0 $0 $7300 7 $2100 $7200 6 $2100 $7100 5 $2000

40 Total costs 1 1 1 4 1 1 0 1 0 1 1 0 0 1 3 2 $7300 7 $7200 6 $7100 5 Lowest total cost

41 Optimal solution 1 1 1 1 1 0 1 0 1 2 1 0 0 1 $7100 5

42 Notes This example is intended to illustrate the principles of unit commitment Some constraints have been ignored and others artificially tightened to simplify the problem and make it solvable by hand Therefore it does not illustrate the true complexity of the problem The solution method used in this example is based on dynamic programming. This technique is no longer used in industry because it only works for small systems (< 20 units)