THE HOLOMORPHI C EMBEDDI NG LOAD FLOW METHOD ( HELM) Antonio Trias - - PDF document

the holomorphi c embeddi ng load flow method helm antonio
SMART_READER_LITE
LIVE PREVIEW

THE HOLOMORPHI C EMBEDDI NG LOAD FLOW METHOD ( HELM) Antonio Trias - - PDF document

THE HOLOMORPHI C EMBEDDI NG LOAD FLOW METHOD ( HELM) Antonio Trias EEI Transmission, Distribution, & Metering Conference Spring 2012 Newport, Rhode Island April 3, 2012 1 Contents 1 . Sum m ary 2 . The problem w ith iterative m


slide-1
SLIDE 1

–1 THE HOLOMORPHI C EMBEDDI NG LOAD FLOW METHOD ( HELM) Antonio Trias

1

EEI Transmission, Distribution, & Metering Conference Spring 2012 Newport, Rhode Island April 3, 2012

Contents

1 . Sum m ary 2 . The problem w ith iterative m ethods 3 . The HELM m ethod: 1 . Holom orphic em bedding

2

2 . The m ethod procedure 4 . Exam ples

slide-2
SLIDE 2

–2 Sum m ary

  • W HAT:

A l l th d f l i th l d fl A novel, general-purpose m ethod for solving the load flow equations of pow er system s of any size.

  • W HY:
  • NR shortcomings

3

  • HELM is deterministic and complete: it ensures that if there is a

solution the method will find it and, conversely, if there is no solution (voltage collapse) the method will unequivocally signal such condition as well.

Sum m ary

  • HOW :

The m ethod is based on a holom orphic em bedding procedure that extends the voltage variables into analytic procedure that extends the voltage variables into analytic functions in the com plex plane. This provides a fram ew ork to obtain and study the solutions using the full pow er of com plex analysis techniques:

  • The holomorphic embedding method provides a non-iterative

procedure for constructing the complex power series of voltages at a well-defined reference point, where it is trivial to

4

identify the correct branch of the multivalued problem

  • Then uses analytical continuation by means of algebraic

approximants to reach the solution. It can be proven that the continuation propagates the chosen branch to the maximal (in logarithmic capacity) possible domain on the complex plane. This proves that the method is deterministic and non- equivocal.

slide-3
SLIDE 3

–3 Sum m ary

  • Relevance:

1 . As an enabler of new real-tim e applications:

  • HELM allows to reliably implement real-time applications that

HELM allows to reliably implement real time applications that perform non-supervised exploratory load flows under all possible conditions

  • Examples: Contingency Analysis, Limit-Violations solvers,

Restoration plan builders

2 . As an enabler of new pow erful analysis tools:

  • new insights into the load flow problem

5

  • new magnitudes with practical application (e.g. new measures
  • f distance to voltage collapse)
  • comprehensive framework for computing the multiple

solutions to the problem

  • exact treatment of any general (algebraic) constraints under

the same methodology

Sum m ary

  • Relevance:

3 . Final com m ents

  • The mathematical foundations (complex analysis, geometric

The mathematical foundations (complex analysis, geometric function theory) are advanced, but the numerical implementation is straightforward

  • Performance characteristics that make it competitive and even

superior to fast-decoupled load flow algorithms, thus making it a good general-purpose load flow method.

  • The method has been implemented in industrial-strength EMS

products now operating at several large utilities in Europe,

6

Mexico, and the US, and has been granted two US Patents (7519506 and 7979239).

slide-4
SLIDE 4

–4 Contents

1 . Sum m ary 2 . The problem w ith iterative m ethods 3 . The HELM m ethod: 1 . Holom orphic em bedding

7

2 . The m ethod procedure 4 . Exam ples

The problem w ith iterative m ethods

HELM New ton-Raphson

Physical Solution P Q V Unstable Solution

8

W hile HELM alw ays finds the physical solution, N-R can converge to either the correct physical solution, an unstable solution, a w rong solution, or not converge at all. Physical Unstable Erroneous Non-Convergence Voltage Collapse Point P,Q

Note: Every point w ithin the circum ference represents a possible starting point for New ton-Raphson calculation.

slide-5
SLIDE 5

–5 The problem w ith iterative m ethods

For m ore inform ation about the “fractality” problem :

  • Bibliography:
  • R. Klump and T.Overbye, “A new method for finding low-

voltage power flow solutions,” in Power Engineering Society Summer Meeting, 2000. IEEE, vol. 1, 2000, pp. 593–597.

  • J. Thorp and S. Naqavi, “Load-flow fractals draw clues to

erratic behaviour,” IEEE Comput. Appl. Power, vol. 10, no. 1 pp 59–62 jan 1997

9

1, pp. 59 62, jan 1997.

  • Gridquant, “Convergence Issues with Newton-Raphson

Method”, 2012. http: / / www.gridquant.com/ technology/

Contents

1 . Sum m ary 2 . The problem w ith iterative m ethods 3 . The HELM m ethod: 1 . Holom orphic em bedding 2 . The m ethod procedure

1 0

4 . Exam ples

slide-6
SLIDE 6

–6 The HELM m ethod

1 1

The HELM m ethod Two bus system

V  V0  Z S * V * ; Z  R  jX ; S  P  jQ

The load flow equation for the tw o-bus system : I ntroducing dim ensionless variables:

  ZS* |V |2 ;  R  XQ  RP |V |2 ;  I  XP  RQ |V |2

U  V V0

1 2

|V0 | |V0 | |V0 |

U  1  U *

The load flow equation in its m ost essential form :

slide-7
SLIDE 7

–7 The HELM m ethod Two bus system

Exact solution( s) :

U 1  1 

2

U

subject to the condition:

U R  2  4   R   I

2

; U I   I

D  1 4   R   I

2  0

1 3

Rem em ber that the load flow problem is m ulti-valued in general. How ever, only one solution actually m akes sense in the real system . All the other so-called “low - voltage” solutions w ould be unstable in any w ell- designed pow er system .

4

Holom orphic Em bedding

Holom orphic em bedding:

  s

w ith

U(s)  U*(s*) U(s)  1 s U(s) U(s)  1 s *

Em bedded load flow equations:

1 4

U(s)  1 U(s)

This is now a pair of functional equations involving tw o holom orphic functions U, Û of the com plex variable s. At the lim it s= 1 , w e recover the original load flow equation At the lim it s= 0 , the system is trivially solvable

slide-8
SLIDE 8

–8 The m ethod procedure

Once the equations are under the proposed em bedding, this is how the m ethod w orks:

  • Since U(s) is holom orphic, consider the pow er series

expansion about s= 0. The embedded equations allow us to find all the coefficients of the power series as the solution to a succession of linear system s (1x1 in the two-bus case), order after order. The solution of each system yields the information to find the right-hand side for the next one.

1 5

  • Using analytic continuation techniques (Padé

Approximants), the solution at s= 1 can be constructed. Stahl’s theorem guarantees that the result is maximal: if there is a solution, the Padé Approximants will find it; if there is no solution, the Padé Approximants will signal this condition as well.

The m ethod procedure

Once the equations are under the proposed em bedding, this is how the m ethod w orks:

  • Furthermore the solution is analytically continued from the

well-defined reference solution at s= 0: it is guaranteed to be the correct operative solution.

1 6

slide-9
SLIDE 9

–9 References on HELM

For further reference:

  • Bibliography
  • “The Holomorphic Embedding Load Flow Method”, A.
  • Trias. Accepted for publication in IEEE Trans. Power

Syst., PES General Meeting, July 2012.

  • “Two-bus load flow: exact developments”, A. Trias.

Notes available at URL:

1 7

Notes available at URL: http: / / www.gridquant.com/ technology/

Contents

1 . Sum m ary 2 . The problem w ith iterative m ethods 3 . The HELM m ethod: 1 . Holom orphic em bedding 2 . The m ethod procedure

1 8

4 . Exam ples

slide-10
SLIDE 10

–10 10 Exam ples ( ALPS off-line Sim ulation)

1 9

Exam ples ( ALPS off-line Sim ulation)

2 0

slide-11
SLIDE 11

–11 11 HELM-based real-tim e applications

AGORA: A grid management system for Real‐Time operations

RT Restoration

  • m onitoring

HELM Load Flo Contingency State Estimator PV/QV Curves OPF Simulator Solver

  • Lim. Viol.

Solver

2 1

  • analysis
  • ptim al decision m aking
  • full integration w ith SCADA

Load Flow Analysis OPF

SCADA

AGORA: Advanced Grid Observation Reliable Algorithm s

A grid m anagem ent system for Real-Tim e operation “in Real-Tim e operation, w e can’t afford to have a convergence rate of less than 1 0 0 % ”

  • AGORA applications are based on new, ground-breaking

algorithmic advances

  • Based on HELM, a new, direct method (US Patent No. 7

519 506, April 2009) for solving the load-flow equations has been developed

  • Several industry-first applications:

2 2

  • Several industry first applications:
  • Advanced Topological State Estimation
  • PV-QV Curves on any node in Real time
  • Limits Violations Solver
  • Restoration Solver
slide-12
SLIDE 12

–12 12 AGORA: Advanced Grid Observation Reliable Algorithm s

A grid m anagem ent system for Real-Tim e operation “in Real-Tim e operation, w e can’t afford to have a convergence rate of less than 1 0 0 % ”

  • The solvers give answers in terms of real SCADA actions
  • And they are fully checked via load-flow computation on

the current estimated state

2 3

AGORA Restoration Solver

  • First system in the market

that computes restoration plans in real time

  • “Like a GPS navigator” for

the Operator Search for the Optim al Path

  • Define a state space where the initial

state is the disturbance condition and the final state is the pre-disturbance condition T ansitions a e allo able actions

the Operator

  • Monitors real-time network

state, and re-computes when necessary

  • Allows the user to simulate

and evaluate the plan

  • Transitions are allowable actions
  • Define a transition cost function
  • Consider admissible suboptimal

heuristic (guarantees convergence of the A* algorithm)

Final

2 4

I nitial State Final State

Allow able Maneuvers ( Actions)