Newton s Power Flow algorithm

Similar documents
G S Power Flow Solution

Power Flow S + Buses with either or both Generator Load S G1 S G2 S G3 S D3 S D1 S D4 S D5. S Dk. Injection S G1

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

18.413: Error Correcting Codes Lab March 2, Lecture 8

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

Lecture Notes 2. The ability to manipulate matrices is critical in economics.

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

MOLECULAR VIBRATIONS

The number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Power Flow Analysis for IEEE 30 Bus Distribution System

Functions of Random Variables

Simple Linear Regression

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

General Method for Calculating Chemical Equilibrium Composition

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

PTAS for Bin-Packing

Lecture 9: Tolerant Testing

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3

Multiple Choice Test. Chapter Adequacy of Models for Regression

THE COMPLETE ENUMERATION OF FINITE GROUPS OF THE FORM R 2 i ={R i R j ) k -i=i

AN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Chapter 9 Jordan Block Matrices

TESTS BASED ON MAXIMUM LIKELIHOOD

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Feature Selection: Part 2. 1 Greedy Algorithms (continued from the last lecture)

CHAPTER 4 RADICAL EXPRESSIONS

Special Instructions / Useful Data

Simulation Output Analysis

Analysis of Lagrange Interpolation Formula

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

ESS Line Fitting

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

MA/CSSE 473 Day 27. Dynamic programming

X X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then

4 Inner Product Spaces

ENGI 4421 Propagation of Error Page 8-01

The Necessarily Efficient Point Method for Interval Molp Problems

α1 α2 Simplex and Rectangle Elements Multi-index Notation of polynomials of degree Definition: The set P k will be the set of all functions:

13. Parametric and Non-Parametric Uncertainties, Radial Basis Functions and Neural Network Approximations

Ahmed Elgamal. MDOF Systems & Modal Analysis

Ideal multigrades with trigonometric coefficients

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

Lecture 07: Poles and Zeros

Discrete Mathematics and Probability Theory Fall 2016 Seshia and Walrand DIS 10b

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

ECON 5360 Class Notes GMM

Third handout: On the Gini Index

Multiple Linear Regression Analysis

Chapter 5 Properties of a Random Sample

III-16 G. Brief Review of Grand Orthogonality Theorem and impact on Representations (Γ i ) l i = h n = number of irreducible representations.

1 Onto functions and bijections Applications to Counting

The Mathematical Appendix

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Bayes (Naïve or not) Classifiers: Generative Approach

Kernel-based Methods and Support Vector Machines

Statistics: Unlocking the Power of Data Lock 5

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

3D Geometry for Computer Graphics. Lesson 2: PCA & SVD

Chapter 11 Systematic Sampling

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

The Selection Problem - Variable Size Decrease/Conquer (Practice with algorithm analysis)

LINEARLY CONSTRAINED MINIMIZATION BY USING NEWTON S METHOD

End of Finite Volume Methods Cartesian grids. Solution of the Navier-Stokes Equations. REVIEW Lecture 17: Higher order (interpolation) schemes

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

Lecture Note to Rice Chapter 8

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Wu-Hausman Test: But if X and ε are independent, βˆ. ECON 324 Page 1

Point Estimation: definition of estimators

Objectives of Multiple Regression

Introduction to local (nonparametric) density estimation. methods

Bayes Interval Estimation for binomial proportion and difference of two binomial proportions with Simulation Study

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

1 Solution to Problem 6.40

ENGI 3423 Simple Linear Regression Page 12-01

UNIT 7 RANK CORRELATION

2SLS Estimates ECON In this case, begin with the assumption that E[ i

STA 108 Applied Linear Models: Regression Analysis Spring Solution for Homework #1

Summary of the lecture in Biostatistics

THE ROYAL STATISTICAL SOCIETY HIGHER CERTIFICATE

Loop-independent dependence: dependence exists within an iteration; i.e., if the loop is removed, the dependence still exists.

Chapter 3 Sampling For Proportions and Percentages

L5 Polynomial / Spline Curves

Chapter 14 Logistic Regression Models

4. Standard Regression Model and Spatial Dependence Tests

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Chapter 3. Linear Equations and Matrices

Multi-Step Methods Applied to Nonlinear Equations of Power Networks

A Study of the Reproducibility of Measurements with HUR Leg Extension/Curl Research Line

Lecture 3 Probability review (cont d)

Unsupervised Learning and Other Neural Networks

Chapter 8. Inferences about More Than Two Population Central Values

Beam Warming Second-Order Upwind Method

Analyzing Control Structures

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

Transcription:

Power Egeerg - Egll Beedt Hresso ewto s Power Flow algorthm

Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 2 Calculatos. For the referece bus #, we set : V = p.u. ad δ = 0 For all other buses, =2,3,, guess a soluto, whch could be : V = p.u ad δ = 0 degrees ad set the terato couter m = tlt tlt Defe P ad Q whch s a specfed ow et power jecto at each bus of the power system except the referece bus. Substtute the guessed values of voltages to the power flow equatos ad calculate the followg quattes step 2:

Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 3 Calculatos (2) re re 2. Calculate P ad Q for stace accordg to the equatos, Assumg a polar coordate decomposto all varables: re = cos = P V V y [ δ γ ] re = s = Q V V y [ δ γ ] for all buses except the referece bus,.e.: = 2, 3,...,. where j V = V δ y = y γ = y e γ

Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 4 Calculatos (3) We get the followg dfferece vector: Δy tlt re P2 P2 tlt re P3 P3 tlt re ΔP P P = = tlt re Q2 Q ΔQ 2 tlt re Q3 Q3 tlt re Q Q

Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 5 Calculatos (4) 3. ext calculate the Jacob-matrx : P P J J δ V 2 J = = J3 J4 Q Q δ V where the elemets of the matrx are show o the ext 2 sldes : J ad J 2 are the real power dervatve equatos J 3 ad J 4 are the reactve power dervatve equatos

Power Egeerg - Egll Beedt Hresso 6 The Jacob-Matrx P Wth the J sub-matrx s as follows: δ =, ( ) = V V y s δ δ γ ( ) P = V V y cos δ δ γ = P δ s ( ) = V V y δ δ γ The J 2 sub-matrx s as follows : P V = 2V y cos + V y cos ( γ ) ( δ δ γ ) j j j j j=, j P V = V y cos( δ δ γ )

Power Egeerg - Egll Beedt Hresso 7 The Jacob-Matrx (2) Wth ( ) Q = V V y s δ δ γ j j j j j= The J 3 sub-matrx s as follows : Q δ j=, j ( ) = V V y cos δ δ γ j j j j Q δ = V V y cos( δ δ γ ) The J 4 sub-matrx s as follows : Q V ( γ ) ( δ δ γ ) = 2V y s + V y s j j j j j=, j Q = V y s ( δ δ γ ) V

Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 8 Calculatos (4) 4. Solve the lear set of equatos: ( m+ ) ( m) tlt ( m) re P 2 2 2 P2 δ δ ( m+ ) ( m) tlt ( m) re P 3 3 3 P δ δ 3 P P ( m ) ( m tlt ( m) re + ) P P δ V δ δ tlt ( m) re = Q2 Q2 Q Q ( m+ ) ( m) V2 V tlt ( m) re 2 Q3 Q δ V ( ) ( ) 3 m+ m V3 V 3 tlt ( m) re Q Q ( m+ ) ( m) V V óþetur vetor með leðréttgum Δy P P δ V = Δx Q Q δ V

Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow Calculatos (5) 9 The last matrx ca also be wrtte: P P ΔP δ V Δδ = ΔQ Q Q Δ V δ V ΔP J J2 Δδ = ΔQ J 3 J4 Δ V 5. Calculate a correcto to the varables, or: 6. ad ( m+ ) ( m) ( m) δ δ δ = +Δ. If the chages Calculate the load flow system braches accordg to: * S = P + jq = V I * * * 2 * s p S = V V V Y + V Y ( m) Δ V betwee teratos are sgfcat (>ε) go bac to step # 2, otherwse go to step # 6 Calculate the power of the referece bus ( = ) accordg to the equatos: I ( m+ ) ( m) ( m) = +Δ V V V ad Y p Δδ ( m) Y p [ δ γ ] P = V V y cos re = re = s = S Y s S Q V V y [ δ γ ] I

Power Egeerg - Egll Beedt Hresso 0 Gauss/ewto Comparso of Methods Both are soluto methods to solve a set of olear equatos wth the use of terato. The ewto s method s more dffcult programmg, but gves a faster soluto. The Gauss-method s smple to program but gves a more accurate soluto ad requres a loger soluto tme.

Power Egeerg - Egll Beedt Hresso umber of Iteratos for Load Flow Solutos (Gauss ad ewto) From: Stagg, El-Abad: Computer Methods of Power System Aalyss

Computato Tme for Each Power Egeerg Iterato - Egll Beedt Hresso for 2 Load Flow Methods From: Stagg, El-Abad: Computer Methods of Power System Aalyss

Power Egeerg - Egll Beedt Hresso 3 DC Load Flow (A lear approxmato for stuatos where may estmates of load flow are eeded a short tme such as dyamc models)

Power Egeerg - Egll Beedt Hresso A Lear Approxmato to Load Flow (DC Load Flow) 4 The electrcal power system s ductve ad resstaces are comparatvely small,. e.: g 0 Phase agles dffereces are geerally small betwee adjacet buses a electrcal power system,.e.: δ δ s cosδ Voltages the system uder ormal operatg crcumstaces are close to uty, or.: V

Power Egeerg - Egll Beedt Hresso A Lear Approxmato to Load Flow (DC Load Flow) (2) By relaxg the last assumptos (where voltages are assumed costat ad phase agles small), we get the reactve power quadratc equatos as follows: 5 [ cos( δ δ ) s( δ δ )] P = V V g + b = [ s( δ δ ) cos( δ δ )] Q = V V g b = Q = V V b = However, the reactve power quadratc equatos are ot as terestg as the equatos for real power (see ext slde):

Power Egeerg - Egll Beedt Hresso A Lear Approxmato to Load Flow (DC Load Flow) (3) 6 By assumg that all of the above assumptos are vald, we get the followg DC-load flow equatos, where the relatos of real power ad phase agle s a lear set of equatos: [ cos( δ δ ) s( δ δ )] P = V V g + b = [ s( δ δ ) cos( δ δ )] Q = V V g b = P = b δ = b = = [ δ δ ]

Power Egeerg - Egll Beedt Hresso A Lear Approxmato to Load Flow (DC Load Flow) (4) 7 These equatos are aalogous to equatos for a drect curret crcut, f : P s aalogous to the curret the curret source (Kow for load buses) b s aalogous to the resstace or coductace (Kow for load buses) δ s aalogous to the voltage at each ode. (These are the uow varables) DC power flow s approxmately 0-00 tmes faster tha covetoal" load flow (Also called full AC - load flow), but may have assocated errors.

Power Egeerg - Egll Beedt Hresso 8 Fast Decoupled Load Flow (FDLF) Ths s aother approxmato to the full load flow equatos to crease speed)

Power Egeerg - Egll Beedt Hresso The dervatves wth rectagular 9 represetato of admttaces The load flow equatos wth rectagular represetato of admttaces are repeated here: y = g + jb 2 = + + = ( δ δ ) s ( δ δ ) ( cos ) P V g V V g b 2 = + = ( δ δ ) cos( δ δ ) ( s ) Q V b V V g b [ cos( δ δ ) s( δ δ )] P = V V g + b = [ s( δ δ ) cos( δ δ )] Q = V V g b = or wrtte alteratvely :

Power Egeerg - Egll Beedt Hresso 20 Full AC Dervatves Real Power dervatve equatos (J ad J 2 ) are (wth dex replacg former dex ad dex m replacg former dex ) P δ m ( δ δ ) cos( δ δ ) = V V g s b m m m m m P V m ( δ δ ) cos( δ δ ) = V g s b m m m m P ( δ δ ) cos( δ δ ) = V Vm gm s m bm m δ + m= m P = 2V g + V g cos( δ δ ) + b s( δ δ ) V m m m m m m= m Reactve Power dervatve equatos (J 3 ad J 4 ) are: Q δ m ( δ δ ) s ( δ δ ) = V V g cos b m m m m m Q V m ( δ δ ) cos( δ δ ) = V g s b m m m m Q ( δ δ ) s ( δ δ ) = V Vm gm cos m bm m δ + m= m Q = 2V b + V g s( δ δ ) b cos( δ δ ) V m m m m m m= m

Power Egeerg - Egll Beedt Hresso 2 Decoupled Power Flow Equatos Mae the followg assumptos: δ δ m V cos( δ δ m ) s( δ δ m ) 0 The dervates smplfy to: 0 r m << x m 0 g m P δ = m= m b m P δ m = b m P V = 0 P V m = 0 Q Q = 2 b + ( bm ) = b Q m = 0 V V δ m= m m How do we terpret the above dervatve terms? δ Q m = 0

Power Egeerg - Egll Beedt Hresso To do that we remember the Y bus Matrx 22 22 Composto Rule: The dagoal elemets of the Y bus matrx cosst of the sum of all the (seres ad shut) admttaces coected to the bus questo The off-dagoal elemets of the Y bus matrx cosst of the egatve value of the seres admttaces coectg the 2 buses questo. Therefore shut terms at bus # of : shut, m m= m are the sum of each row or colum of Y. The shut terms are show o the ext slde: Y bus bus y y y = +

V V 3 Power Egeerg - Egll Beedt Hresso A example of shut terms or shut elemets the Y bus matrx The shut terms at bus # are Y s,,3 y = Y + Y shut, p,,2 p,,3 23 I Y p,,3 Y p,3, I 3 Y p,,2 V 4 Y s,,2 Y s,2,4 Y p,4,2 Y p,2, V 2 The seres term betwee bus # ad bus #2 I 2 Y p,2,4 The shut terms at bus #2 are y = Y + Y shut,2 p,2, p,2,4

Power Egeerg - Egll Beedt Hresso 24 We chec out the Shut terms B bus g g2 g b b2 b g2 g22 g 2 b2 b22 b 2 Y j bus = + = Gbus + jb g g g b b b 2 2 Dagoal elemets the magary part of Y (that s B ) are composed of seres ad shut terms, or: b = b + b bus Shut part of b : b = b + b shut, m m= m shut. seres, bus bus Seres terms of b : b = ( b ) seres, m m= m

Power Egeerg - Egll Beedt Hresso 25 Defto of the B matrx: Therefore the shut terms at bus # of : b = b + b = b ( b ) shut, m m m= m= m m are the sum of all elemets each row or colum of B. Compare ths P to the above dervatves: = bm ad δ m m= m B P P = B = J. Sce = bm = ( b bshut, ) = bseres, we fd that: δ δ bus P δ = m= m b bus m ad defe the matrx B : B s the s the egatve of the magary part of the Y bus matrx (B bus ) oly wth seres terms or wth all the shut terms removed from B bus.

Power Egeerg - Egll Beedt Hresso 26 Defto of the B matrx: ow loo at the dervatves: Q V Q = 2 b + ( b ) ad = b ad defe the matrx B : m m m= Vm m Q = B = J. Sce b = b + b ad b = ( b ) ad sce V Q V 3 shut, m seres, m m= m= m m = 2 b + ( b We fd that: m= m ) = ( b + 2 b ) m seres, shut, B s the s the egatve of the magary part of the Y bus matrx (B bus ) wth all the shut terms couted twce (double couted) B bus.

Power Egeerg - Egll Beedt Hresso Comparso wth textboo verso (Saadat) 27 Wth o shut terms (all ls as serse mpedaces) B ad B are smply the egatve of B bus (page 242, Saadat). Ths s the case the textboo. Wth some voltages ot set equal to pu we get equatos (6.75) ad (6.76) Saadat. Wth Q = 0 (6.7) o shut elemets are preset B ad B ther textboo verso. It s a cruder approxmato tha the above where they are ether removed or double couted.

Power Egeerg - Egll Beedt Hresso 28 Soluto wth B ad B Matrces ow terate the decoupled equatos Δδ = ΔV = [ ] ' B ΔP [ B ] '' ΔQ

Power Egeerg - Egll Beedt Hresso 29 Fast Load Flow ΔP J J2 Δδ = ΔQ J2 J4 Δ V J2 0 J 0 3 2 depedet sets of lear equatos. The dmesos for each of them s half the orgal set of equatos ΔP B 0 Δδ ΔQ 0 B Δ V ' = '' ΔP = J Δδ ΔQ = J Δ V 4

Power Egeerg - Egll Beedt Hresso Fast Decoupled Load Flow (FDLF) 30 FDLF s based o a approxmato ad s faster but more accurate tha covetoal power flow The approxmato taes to accout: That voltage a power system s tghtly coected to reactve power That phase agles a power system are tghtly coected to actve power

- securty rule Power Egeerg - Egll Beedt Hresso 3 Fally we meto a very well ow practcal securty rule of thumb: The system ca loose oe ut wthout losg customers or a blacout occurrg the system Ths ut ca be a trasmsso le, a trasformer or a geerator, for stace We specfy that durg a dsturbace customers wll ot be terrupted Ths rule s called the - rule

Hstorc refereces Power Egeerg - Egll Beedt Hresso 32 B. Stott ad O. Alsac, Fast decoupled load flow, IEEE Trasactos o Power Apparatus ad Systems, Vol. PAS-93, Jue 974, pp. 859-869. W. F. Tey ad C. E. Hart, Power Flow Soluto by ewto s Method, IEEE Trasactos o Power Apparatus ad Systems, Vol. PAS-86, o., ov. 967, pp. 449-460. B. Stott, Revew of Load-Flow Calculato Methods, Proceedgs of the IEEE, Vol. 62, o. 7, July 974, pp. 96-929.