Phase Angle Calculations:

Similar documents
Peak Reliability Delivering near real-time phase angle deltas using Inter-Control Center Communication Protocol (ICCP)

Real Time Phasor Data Processing using Complex Number Analysis Methods R.A. DE CALLAFON 1, C.H. WELLS 2 USA

Hunting for Anomalies in PMU Data

WECC Dynamic Probing Tests: Purpose and Results

Math 103 Day 13: The Mean Value Theorem and Tuesday How Derivatives OctoberShape 26, 2010 a Graph1 / 12

Real-Time Wide-Area Dynamic Monitoring Using Characteristic Ellipsoid Method

1 Modular Arithmetic Grade Level/Prerequisites: Time: Materials: Preparation: Objectives: Navajo Nation Math Circle Connection

Algebra Mat: Working Towards Year 6

How Corning Kennebunk Uses PI System Tools to Make Better Use of Data

ECE 241L Fundamentals of Electrical Engineering. Experiment 6 AC Circuits

SOLUTIONS to ECE 2026 Summer 2017 Problem Set #2

PMU-Based Power System Real-Time Stability Monitoring. Chen-Ching Liu Boeing Distinguished Professor Director, ESI Center

5.3 The Fundamental Theorem of Calculus. FTOC statement: Part 1. note: connection Increasing/decreasing functions and their derivatives

Physics Fall 2006 Laboratory 5: Rotational Dynamics

Enhancement of State Estimation Results using Phasor Measurements. Presented to: NASPI Work Group Meeting March 6-7, 2008, New Orleans, LA

Unit Trip April 20, 2011, 13:40:53 CDT. McDonald Angle Relative to U.T. Austin Second

Chapter 12 Variable Phase Interpolation

Systems of Linear Equations

PDQ Tracker High Level Requirements

5.5 Special Rights. A Solidify Understanding Task

Safety evaluations of level crossings

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING

Key Performance Task

Digital Power System Protection

On Computing Power System Steady-State Stability Using Synchrophasor Data

CpE358/CS381. Switching Theory and Logical Design. Class 16

Lecture 5: Using electronics to make measurements

A Supervisory Approach towards Cyber- Secure Generator Protection

LCR Series Circuits. AC Theory. Introduction to LCR Series Circuits. Module. What you'll learn in Module 9. Module 9 Introduction

Optimal PMU Placement

MATH 423/ Note that the algebraic operations on the right hand side are vector subtraction and scalar multiplication.

OSIsoft Data Compression Analysis

MATH 2331 Linear Algebra. Section 1.1 Systems of Linear Equations. Finding the solution to a set of two equations in two variables: Example 1: Solve:

Linear system of equations

Lecture 8: Signal Detection and Noise Assumption

Inner product spaces. Layers of structure:

Intro Vectors 2D implicit curves 2D parametric curves. Graphics 2012/2013, 4th quarter. Lecture 2: vectors, curves, and surfaces

Discovery Through Situational Awareness

Homework 2; due on Thursday, October 4 PY 502, Computational Physics, Fall 2018

New Approaches for the Design of Synchrophasor Estimation Systems

Using Algebra Fact Families to Solve Equations

Multi-Step Equations and Inequalities

Linear Algebra MATH20F Midterm 1

Numerical Methods I Solving Nonlinear Equations

BETWEEN PAPERS PRACTICE SET 3 of 4 (F&H)

Phasor mathematics. Resources and methods for learning about these subjects (list a few here, in preparation for your research):

Introduction EE 224: INTRODUCTION TO DIGITAL CIRCUITS & COMPUTER DESIGN. Lecture 6: Sequential Logic 3 Registers & Counters 5/9/2010

Power system modelling under the phasor approximation

PHY 221 Lab 3 Vectors and Motion in 1 and 2 Dimensions

Lecture 1: Systems of linear equations and their solutions

On-Line TSA and Control at EMS Control Centers for Large Power Grids

I. Signals & Sinusoids

Lecture 11. Andrei Antonenko. February 26, Last time we studied bases of vector spaces. Today we re going to give some examples of bases.

AP Calculus AB UNIT 1: PRECALCULUS REVIEW UNIT 2: BRIDGE TO CALCULUS LESSON 1: INTRO TO CALCULUS LESSON 2: FUNCTIONS

Chapter 5 Synchronous Sequential Logic

14:332:231 DIGITAL LOGIC DESIGN. Combinational Circuit Synthesis

IEEE Power and Energy Society General Meeting, July 2010, Minneapolis, MN USA

Lecture 5: Using electronics to make measurements

MATH240: Linear Algebra Review for exam #1 6/10/2015 Page 1

HPC and High-end Data Science

Music 270a: Complex Exponentials and Spectrum Representation

Name. Use Two-Color Counters to model each addition problem. Make pairs of red and yellow counters. Find the sum.

Math 360 Linear Algebra Fall Class Notes. a a a a a a. a a a

The General Linear Model. How we re approaching the GLM. What you ll get out of this 8/11/16

3 December Lesson 5.5

Linear Equations in Linear Algebra

Adding Integers with Different Signs. ESSENTIAL QUESTION How do you add integers with different signs? COMMON CORE. 7.NS.1, 7.NS.

Fluid Statics. Pressure. Pressure

SPEECH ANALYSIS AND SYNTHESIS

Chapter 7 PHASORS ALGEBRA

Finite Math - J-term Section Systems of Linear Equations in Two Variables Example 1. Solve the system

Designing Information Devices and Systems II Spring 2018 J. Roychowdhury and M. Maharbiz Discussion 6B

MA 242 LINEAR ALGEBRA C1, Solutions to First Midterm Exam

1 Invariant subspaces

Intro Vectors 2D implicit curves 2D parametric curves. Graphics 2011/2012, 4th quarter. Lecture 2: vectors, curves, and surfaces

Data Converter Fundamentals

Communication Signals (Haykin Sec. 2.4 and Ziemer Sec Sec. 2.4) KECE321 Communication Systems I

The ICESat 2 Mission Laser altimetry of ice, clouds and land elevation

The Gram-Schmidt Process

Putnam Pigeonhole Principle October 25, 2005

SVSM User Training and Example Sequences

CHAPTER 1 LINEAR EQUATIONS

Brief Review of Exam Topics

Unit #5 Applications of the Derivative Part II Homework Packet

2.0 COMPLEX NUMBER SYSTEM. Bakiss Hiyana bt Abu Bakar JKE, POLISAS BHAB 1

ω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the

Distortion Analysis T

Math 235: Linear Algebra

ROBERTO BATTITI, MAURO BRUNATO. The LION Way: Machine Learning plus Intelligent Optimization. LIONlab, University of Trento, Italy, Apr 2015

PATTERN RECOGNITION AND MACHINE LEARNING

Revision of Basic A.C. Theory

1. Introduction to commutative rings and fields

Contents lecture 6 2(17) Automatic Control III. Summary of lecture 5 (I/III) 3(17) Summary of lecture 5 (II/III) 4(17) H 2, H synthesis pros and cons:

Distributed Real-Time Electric Power Grid Event Detection and Dynamic Characterization

Computers also need devices capable of Storing data and information Performing mathematical operations on such data

Carleton College, winter 2013 Math 232, Solutions to review problems and practice midterm 2 Prof. Jones 15. T 17. F 38. T 21. F 26. T 22. T 27.

REAL TIME PMU-BASED STABILITY MONITORING ZIJIE LIN. A thesis submitted in partial fulfillment of the requirement for the degree of

On the Use of PMUs in Power System State Estimation

Detection theory. H 0 : x[n] = w[n]

a k cos kω 0 t + b k sin kω 0 t (1) k=1

Transcription:

Phase Angle Calculations: Considerations & Use Cases KEN MARTIN, EPG EATT PHASE ANGLE PAPER TEAM NASPI WORK GROUP MEETING OCTOBER 2016

Acknowledgment Jeff Anderson, BPA Dan Brancaccio, Bridge Jeff Dagle, PNNL Harold Kirkham, PNNL Jim Kleitsch, ATC Dmitry Kosterev, BPA Ken Martin, EPG Jim McNierney, NYISO Vincent Minden, BPA Ryan Nice, PJM Ryan Quint, NERC Alison Silverstein, NASPI Kyle Thomas, DVP Mani Venkatasubramanian, WSU Chuck Wells, OSIsoft

Phase Angle Formulation PMU reports measured values time-tagged to UTC Amplitude Phase Angle Frequency Rate of change of frequency (ROCOF) Signal can be modeled as x(t) and the Synchrophasor is calculated relative to the nominal system frequency synchronized to UTC: x(t) = X m (t) cos(2πf 0 t + φ(t)) The synchrophasor is: X (t) = (X m (t)/ 2)e j(φ(t)) The synchrophasor includes all the dynamic changes in amplitude and phase angle Off-nominal frequency appears as a change in phase angle

Frequency & Angle Example Off-nominal frequency f = f 0 + Δf: Signal formula is: x(t) = X m cos(2π(f 0 + Δf) t + φ) = X m cos(2πf 0 t + 2πΔf t + φ) Then the synchrophasor is: : X (t) = X m / 2 ej(2πδf t + φ) Constant f off-nominal results in angle increasing/decreasing linearly φ(t)

PMU Window Measurement A PMU examines the signal over short time (measurement window) Generally a few cycles It estimates the phase angle at each point and resolves the answer to a principal value θ where π < θ π. The PMU only estimates phasor angle, it does not track time error In a time plot, PMU may appear to wrap but that is not something a PMU does It is an artifact of the way we resolve angle to the principal value

Accumulated Phase As frequency varies off-nominal, phase angle will vary System frequency > nominal, the phase angle increases System frequency < nominal, the phase angle decreases Phase angle may decrease for period of time and then increase for some time Accumulated phase can be tracked based on some starting point

The Issues To calculate phase angle between two measurements, we have to be sure they do not span the 180 rollover As they rotate with off-nominal frequency, they do not roll-over at the same time and create a spike in the measurement We may want to record the accumulated system time error through the accumulated phase angle

Solution approaches Each time angle differences are calculated, resolve them into their principal value (subtract then unwrap) Monitor each angle from some starting point and add the full circle (360 ) to an accumulated phase angle value each time it rolls over (unwrap the angle). Calculate all angle differences from these unwrapped angles Starting points have to be unwrapped initially Both approaches assume there will be no real angle differences that are greater than 180. In fact, given any two angles, it is impossible to determine angle differences >180 without additional information Angles at intermediate points Flow calculations with equipment and line parameters (state estimation)

Angle Wrapping Definition

Angle Unwrapping Definition

Unwrap Then Subtract

Subtract Then Wrap

Vector Math This approach lets the computer math formulation resolve the angle. This works because vector (and complex number) math works in rectangular coordinates which do not have a roll over and resolves the answer to the principal value.

Phase Angle Difference Calculation Considerations

Accumulated Phase Angle Difference Phase angle reported by a PMU is defined as the principal value (PV) of the signal relative to the reference waveform synchronized to UTC. PMU is not designed to keep track of evolution of relative phase of two signals since some arbitrary time zero. A PMU could be built to track accumulated phase The user would have to establish a start time Accumulated phase can become arbitrarily large, so a reset mechanism is required Phasor value reporting would need modification to report very large values Accumulated phase can be calculated or obtained from the PMU signals by unwrapping the measurement over a period of time. This is effectively an alternative to calculating system time error used in SCADA systems

Phase Angle Distance Shows how phase angle can increase to large values over time

Rollover Counter (ROC) For keeping track of the number of roll-overs that have occurred

Concluding Remarks Applications must account for principal value of phase angle reported by PMU resolving to interval (+180,-180] in any calculation Applications should document and clearly define angle wrapping and unwrapping they use for industry reference Commonly used angle difference calculation methods explained in detail. Vector math is gold standard computationally complex but useful for offline applications Subtract then Wrap and Unwrap then Subtract both accurate and their pros and cons Simple arithmetic difference should be avoided Considerations include speed, accuracy, CPU usage, data storage, missing data, etc. Phase Angle Distance and Rollover Counter (ROC) introduced