Progressive Hedging for Stochastic Economic Dispatch with AC Power Flow

Similar documents
ECONOMIC OPERATION OF POWER SYSTEMS

Optimization Methods MIT 2.098/6.255/ Final exam

Study on Coal Consumption Curve Fitting of the Thermal Power Based on Genetic Algorithm

DG Installation in Distribution System for Minimum Loss

Optimal Design of Accelerated Life Tests with Multiple Stresses

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 6 - Small Signal Stability

DESCRIPTION OF THE SYSTEM

Scheduling under Uncertainty using MILP Sensitivity Analysis

AP Calculus AB 2006 Scoring Guidelines Form B

Mixtures of Gaussians and the EM Algorithm

EE260: Digital Design, Spring n Binary Addition. n Complement forms. n Subtraction. n Multiplication. n Inputs: A 0, B 0. n Boolean equations:

Mixed Acceptance Sampling Plans for Multiple Products Indexed by Cost of Inspection

Computational Analysis of IEEE 57 Bus System Using N-R Method

The target reliability and design working life

Notes for Lecture 5. 1 Grover Search. 1.1 The Setting. 1.2 Motivation. Lecture 5 (September 26, 2018)

Wind Energy Explained, 2 nd Edition Errata

Exercises and Problems

AP Calculus BC 2005 Scoring Guidelines

Chapter 2 Feedback Control Theory Continued

An Effective Biogeography Based Optimization Algorithm to Slove Economic Load Dispatch Problem

NUMERICAL METHODS FOR SOLVING EQUATIONS

w (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.

Recent Magnetic Measurement Activities at BNL

Economic Evaluation of Transmission Expansion for Investment Incentives in a Competitive Electricity Market

Time-Domain Representations of LTI Systems

LECTURE 5 PART 2 MOS INVERTERS STATIC DESIGN CMOS. CMOS STATIC PARAMETERS The Inverter Circuit and Operating Regions

Beyond simple iteration of a single function, or even a finite sequence of functions, results

A New Solution Method for the Finite-Horizon Discrete-Time EOQ Problem

Research Article A New Second-Order Iteration Method for Solving Nonlinear Equations

EE692 Applied EM- FDTD Method One-Dimensional Transmission Lines Notes- Lecture 4

OPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

THE EFFECT OF DISPERSED GENERATION UNITS WITH AND WITHOUT CAPACITOR BANKS ON THE SYSTEM LOSSES AND THE VOLTAGE REGULATION.

Complex Analysis Spring 2001 Homework I Solution

Bayesian Methods: Introduction to Multi-parameter Models

Live Line Measuring the Parameters of 220 kv Transmission Lines with Mutual Inductance in Hainan Power Grid

Outline. Linear regression. Regularization functions. Polynomial curve fitting. Stochastic gradient descent for regression. MLE for regression

MATHEMATICS. The assessment objectives of the Compulsory Part are to test the candidates :

Research Article Single-Machine Group Scheduling Problems with Deterioration to Minimize the Sum of Completion Times

Efficiency of Linear Supply Function Bidding in Electricity Markets

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

subject to A 1 x + A 2 y b x j 0, j = 1,,n 1 y j = 0 or 1, j = 1,,n 2

Surrogate Affine Approximation based Co-optimization of Transactive Flexibility, Uncertainty, and Energy

Application of Energy Storage Systems for Frequency Regulation Service

Integer Linear Programming

Math 10A final exam, December 16, 2016

Introduction to Machine Learning DIS10

Power supplies for parallel operation Power supplies

A Nash Equilibrium Need Not Exist in the Locational Marginal Pricing Mechanism

A NEW APPROACH TO SOLVE AN UNBALANCED ASSIGNMENT PROBLEM

- if lump-sum taxes are infeasible, level of public good provision must take account of the method of finance

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Dr. Maddah ENMG 617 EM Statistics 11/26/12. Multiple Regression (2) (Chapter 15, Hines)

Control Charts for Mean for Non-Normally Correlated Data

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

AP Calculus BC 2007 Scoring Guidelines Form B

Chapter 23. The Economics of Resources. Chapter Outline. Chapter Summary

Multilayer perceptrons

Vector Quantization: a Limiting Case of EM

Chapter 9: Numerical Differentiation

CS284A: Representations and Algorithms in Molecular Biology

Data Assimilation. Alan O Neill University of Reading, UK

Robust Resource Allocation in Parallel and Distributed Computing Systems (tentative)

Oblivious Gradient Clock Synchronization

Sensitivity Analysis of Daubechies 4 Wavelet Coefficients for Reduction of Reconstructed Image Error

True Nature of Potential Energy of a Hydrogen Atom

Generating subtour elimination constraints for the Traveling Salesman Problem

ROLL CUTTING PROBLEMS UNDER STOCHASTIC DEMAND

The Firm and the Market

AP Calculus BC 2011 Scoring Guidelines Form B

Product measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.

c. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area?

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Application of Imperialist Competitive Algorithm to Solve Constrained Economic Dispatch

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Analysis of Experimental Measurements

MASSACHUSETTS INSTITUTE OF TECHNOLOGY

ABSTRACT 1. INTRODUCTION

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare

INVESTIGATING THE CORRELATION BETWEEN WIND AND SOLAR POWER FORECAST ERRORS IN THE WESTERN INTERCONNECTION

Linear Support Vector Machines

PRELIM PROBLEM SOLUTIONS

Economics 241B Relation to Method of Moments and Maximum Likelihood OLSE as a Maximum Likelihood Estimator

) +m 0 c2 β K Ψ k (4)

Linear Programming and the Simplex Method

Variable selection in principal components analysis of qualitative data using the accelerated ALS algorithm

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Online Convex Optimization with Ramp Constraints

Short-circuit power as important reliability factor for power system planning

TESTING OF THE FORCES IN CABLE OF SUSPENSION STRUCTURE AND BRIDGES

Problem of Technical Losses in Natural Gas Distribution and Transportation Systems of Georgia

Chapter 2 Motion and Recombination of Electrons and Holes

Optimization Methods: Linear Programming Applications Assignment Problem 1. Module 4 Lecture Notes 3. Assignment Problem

OBJECTIVES. Chapter 1 INTRODUCTION TO INSTRUMENTATION FUNCTION AND ADVANTAGES INTRODUCTION. At the end of this chapter, students should be able to:

Castiel, Supernatural, Season 6, Episode 18

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

Probabilistic analyses in Phase 2 8.0

Notes on the GSW function gsw_geostrophic_velocity (geo_strf,long,lat,p)

Transcription:

hotos placed i horizotal positio with eve amout of white space betwee photos ad header rogressive Hedgig for Stochastic Ecoomic Dispatch with AC ower Flow Joh D. Siirola, Cesar A. Silva-Moroy 13 th Iteratioal Coferece o Stochastic rogramg Bergamo, Italy 1 July, 013 Sadia Natioal Laboratories is a multi-program laboratory maaged ad operated by Sadia Corporatio, a wholly owed subsidiary of Lockheed Marti Corporatio, for the U.S. Departmet of Eergy s Natioal Nuclear Security Adistratio uder cotract DE-AC04-94AL85000. SAND NO. 011-XXXX

Outlie Operatig a ower Grid ACOF formulatio I rectagular rectagular polar Stochastic ACOF with Reewable Eergy Resources Computatioal Results Coclusios

Operatig the electric power grid Coceptually simple demad geeratio losses t T I practice, this is complicated by Dyamic costraits (ramp rates Trasmissio limitatios Security (reliability requiremets Market costraits

How do utilities implemet this i practice?

Day-Ahead Uit Commitmet (SCUC D-8h Day-Ahead Eergy Market (DAEM or DAM Clears demad bids ad supply offers at 1600h o the day prior to the operatig day roduces: Hourly schedules for the ext operatig day for market participats (i.e., geeratio ad demad Hourly iterchage schedules Hourly day-ahead Locatioal Margial rices (LMs Some markets optimize reserve requiremets simultaeously

Reliability Uit Commitmet RUC (SCUC D-h Reliability Assessmet (Reserve Adequacy Aalysis - RAA Miimize additioal start-up ad o load costs to provide sufficiet capacity to satisfy the forecasted load plus the operatig ad replacemet reserve requiremets Clears ISO forecasted load at 00h DAM commitmets are respected roduces: Additioal commitmets Updated geerator dispatch poits

Look-Ahead Ecoomic Dispatch (H-1h SCED with ability to brig olie fast start resources It s iteded to meet itra-hour power fluctuatios Update load ad variable geeratio forecasts It produces: Geerator setpoits Schedule of fast start uits

ower flow i a electric grid (Real ower ( is the itegrated product of voltage ( ad curret (I For AC circuits with oly resistive loads = RMS * I RMS Reactive loads (iductors & capacitors shift the I- phase; worst case: = 0 Apparet ower (S is RMS * I RMS Govers geeratio costs ad trasmissio costraits (S ½ is the Reactive ower (

Solvig Ecoomic Dispatch I practice, operators liearize the uderlyig system Basis of Stochastic UC preseted i previous talk DCOF Small agle assumptio for voltage agle Omit reactive power 1, = B ( θ θ 1 Yields suboptimal (ad frequetly ifeasible solutios Actual system voltage agles ca exceed 30 degrees

Formulatig the ACOF ACOF: Detere,,, I for every lie, geerator i grid Explicitly capture etwork topology Miimize total geeratio costs Multiplicity of formulatios Choice of basis (for etwork costraits sigificat: polar: rectagular: I rectagular: Rectagular basis,,, θ,, v r, v j,, i r, i j, v r, v j si(w+ θ = v r si(w + v j si(w+π/., I i the complex plae

ACOF- (polar cos si ( si cos (.. ( m m m d mk m mk m mk m d mk m mk m mk m B G B G s t c

ACOF- (rectagular * * arcta( arcta( ( ( (.. ( m r m j m r j m j j r r k mk mk d d v v v v s q p j Y j s t j c

ACOF-I (rectagular arcta( arcta( ( ( ( (.. ( m r m j m r j m j j r r k mk d j r r j d j j r r v v v v i i I I I I Y I s t j c

Stochastic ACOF Motivatio Icreased peetratio of reewable (odispatchable eergy impacts fossil-fuel plats, i particular uit cyclig. Sudde availability / loss of reewable geeratio must be absorbed by dispatchable geeratio Curtailig reewables frequetly disallowed by cotract / regulatio 0% ad 30% reewable peetratio for a easy week i July. Reproduced from NREL 010 Wester Wid ad Solar Itegratio Study

Stochastic ACOF Motivatio Uits that are desiged for base-load operatio are beig employed to help absorb reewable plats variability Results i icreased O&M costs ad higher forced outage rates How to dispatch geeratio to imize the cyclig of slow geerators while takig ito accout the stochastic ature of wid ad solar plats? 0% ad 30% reewable peetratio for a challegig week i April. Reproduced from NREL 010 Wester Wid ad Solar Itegratio Study

Addig Reewables to ACOF Reewable plats are a sigificat source of ucertaity They do ot use fuel to provide power No structural chages to objective fuctio (zero cost coefficiets Reewables are modeled as must-take resources egative real power loads gg where ww ll 0 " Î N,t Î T g w mm dd G are the geerators at ode W are the reewable geerators at ode L are the lies coected ito ode M are the lies out of ode l m d

Addig Reewables to ACOF Moder reewable power plats ca cotrol their reactive power output, thus is a variable Dispatchable å gîg w g + w + å l - å m - å d = 0 å wîw lîl durig periods whe geeratio is ormally off-lie. Additioally, reactive power limits are give by their power electroics iterface ratig (MA w + w MA ratig mîm reactive power capability. Some iverters have the capability to a power, if eeded, provided that curret ad teral voltage ratigs a Cosiderig that iverter cost is related to curret ratig, provisio of output meas that the iverter eeds to be larger for the same plat M a higher cost compared to existig idustry practice. Figure 6 shows iverter based o curret limits oly. Normally, this iverter would b power factor operatio (1 i Figure 6. Iverters would be able to p "w Î W,t Î T " Î N,t Î T power whe it operates at a power levels lower tha 1 (e.g.,. I provide dîd reactive power support at zero power, similar to a STATCOM fuctioality is ot stadard i the idustry. iverters are typicall grid at ight, i which case the iverter-based reactive power capabil practice could, of course, be modified, if site coditios dictate the us + - 1 Figure 6. Reactive power capability of a iverter (red curve

Stochastic ACOF Formulatio Objective s. t. tt C gt gg ( t C gt a b ( g t 0 g g, ( t c ( g, ( t k( gt Add explicit rampig costs (k First-stage: geerator set poits ( g at the first time step Secod-stage: rampig of the uits through time Base load uits have high k values eakers have low k value

Case Studies IEEE Test cases 118-bus system 118 buses 186 lies 54 geerators 300-bus system 300 buses 411 lies 69 geerators 4 period model Implemetatio (rectagular yomo / ys Ipopt 3.10.1 w/ MA7 4x Itel Xeo E5-4650 (3 cores @.70GHz 51 GB RAM

Extesive form ot tractable Target: <<300 secods Disclaimer: Sigificat tuig opportuities remai Warm-startig / iitializatio Model presolve (ot curretly available i yomo

rogressive Hedgig arallel solutio MI o sigle machie Weak scalig Extremely tight covergece costrait (1E-5 10x 80x improvemet i wall clock time

rogressive hedgig: tuig ρ Surprisigly high rho required (118 bus, 16 sceario, 1e-5 covergece

Impact of ocovexity

Impact of ocovexity (

Summary Stochastic ACOF has the potetial to sigificatly impact grid operatios Itegratio of reewable (odispatchable geeratio ricig rampig products? Extesive Form covergece a cocer rogressive Hedgig (with appropriate tuig effective at gettig high-quality solutios i tractable time Sigificat tuig, formulatio, ad algorithmic opportuities remai