CÁTEDRA ENDESA DE LA UNIVERSIDAD DE SEVILLA

Transcription

1 Detection of System Disturbances Using Sparsely Placed Phasor Measurements Ali Abur Department of Electrical and Computer Engineering Northeastern University, Boston CÁTEDRA ENDESA DE LA UNIVERSIDAD DE SEVILLA Escuela Técnica Superior de Ingeniería de Sevilla 29 de mayo de 2015

2 Greetings from Northeastern University Location: Boston, Massachusetts, USA Web: Student population Undergraduate: ~16,000 Graduate: ~5000 International: 15% (125 countries) Colleges and schools College of Arts, Media, and Design College of Business Administration College of Computer and Information Science College of Engineering Bouvé College of Health Sciences College of Professional Studies College of Science College of Social Sciences and Humanities School of Law

3 Talk Outline Phasor Measurement Units (PMUs) External Network Model Error Identification Fault Detection and Location Closing remarks

4 Phasor measurement units (PMU) Synchronized Phasor Measurements and Their Applications by A.G. Phadke and J.S. Thorp Invented by Professors Phadke and Thorp in They calculate the real-time phasor measurements synchronized to an absolute time reference provided by the Global Positioning System. PMUs facilitate direct measurement of phase angle differences between remote bus voltages in a power grid.

5 Phasor Measurement Units (PMU) Phasor Data Concentrators (PDC) [*] PMU 2 PMU 1 Antenna CT GPS CLOCK PT PMU 3 PDC CORPORATE PDC REGIONAL PDC PMU 4 PMU K DATA STORAGE & APPLICATIONS DATA STORAGE & APPLICATIONS [*] IEEE PSRC Working Group C37 Report

6 Reference phasor +ϴ 1 --ϴ 2 --ϴ 3 +ϴ 1 --ϴ 2 Reference Reference --ϴ 3

7 Measurements provided by PMUs I phasors PMU V phasor ALL 3-PHASES ARE TYPICALLY MEASURED BUT ONLY POSITIVE SEQUENCE COMPONENTS ARE REPORTED V A V B VC = [T] V 0 V + V V + = 1 3 [V A+ V B + 2 V C ] = e jjπ 3

8 Solving for the Unknown Sparse Vector Underdetermined set of equations: C X = [b] P x N N x 1 = P x 1 X is known to be k sparse, i.e. k out of N entries are known to be significantly larger than the remaining (N-k) entries.

9 Two Examples Line Outage Detection Few equations, large number of candidate lines, one of which is disconnected. Fault detection and location Few equations, large number of candidate branches, one of which is faulted. Further complication: unknown fault location along the suspect branch.

10 Line Outage Detection Problem: Detect and identify topology changes (line switching) in the external system by monitoring the internal system. Assume limited access to external system PMUs. EXTERNAL SYSTEM [2] INTERNAL SYSTEM [1] REAL-TIME INTERNAL SYSTEM STATE REAL-TIME INTERNAL TOPOLOGY BASE CASE EXTERNAL STATE BASE CASE EXTERNAL TOPOLOGY

11 Problem Formulation BEFORE AFTER AFTER

12 Problem Formulation EXTERNAL SYSTEM [2] INTERNAL SYSTEM [1]

13 Problem Formulation EXTERNAL SYSTEM [2]?? INTERNAL SYSTEM [1]

14 Linear Approximation: Problem Formulation Partitioned into internal 1 and external 2 networks:

15 Line Switching Modeled by Bus Injections LINE IS OPEN k Pkm m P km - P km k m P km = 0

16 Line Switching Modeled by Bus Injections Bθ 0 = P 0 ( B + B) θ 1 = P 0 Bθ 1 = P 0 P Sparse Bus Injection Vector Pre-outage topology Post-outage state of the system k k b m -b [ 0 0 p 0 0 p 0 ] T P = 0 ΔB= m -b b

17 Linearized Problem Formulation

18 Linearized Problem Formulation J Internal Only Buses Boundary Buses

19 Linearized Problem Formulation

20 Conversion to Binary Variables

21 Optimization Problem min j u j + v j s. t. J = Λ. p. X + u v X j = 2 j e bb = u j ii e bb > 0 e bb = v j ii e bb < 0 v j = u j = 0 ii e bb = 0

22 Incorporating Losses ρ = min (u j + v j ) j s. t. J = Λ. p. A + Λ. q. B + u v A j = 2 j B j = 1 j j, A j B j j, u j + v j < p > 0 q > 0 p 10

23 Example: 118-Bus System / Two Zones

24 Example: 118-Bus System / Two Zones Internal System External System

25 Identified External Line Outages [DC Approximation, No Load Changes] Line outage P (Pre-outage flow on removed line) Solution found for X -1 at bus +1 at bus

26 Lossy Case Simulation Results [Different Line Outages in the External System] Line outage p Solution found for A q e +1@ bus -1@ bus

27 Lossy Case Simulation Results [Simultaneous topology and load/generation changes] Line outage p Solution found for A q e +1@ bus -1@ bus @

28 Fault location problem Fault occurs, typically along a transmission line, less frequently right at a bus / substation Determine the faulted line or line segment Determine the distance between one of the line terminals to the point of fault Branch l Distance x

29 Fault location methods Power frequency / impedance based methods Record fault transients at line terminals Filter HF signals, estimate impedance to fault Traveling wave based methods High frequency sampling ( > 20KHz ) Capture wave front arrival instants using synchronized sensors AI / pattern recognition / machine learning based methods

30 Motivation for this work +ΔV k V k V m +ΔV m I f V k +ΔV k I k I m V m +ΔV m

31 Virtual injections and fault location I 1 + I 2 = I 3 I 2 I 1 = Z 13 Z 23

32 Fault Detection and Location Can line faults be detected and located using a sparse set of phasor voltage measurements? Can all types of faults in three-phase systems be identified? How do we ensure that there are sufficient measurements to detect and identify faults on any line in the power grid?

33 Faulted Line 2-5 AB-G Fault 10-bus Example

34 Network Partitioning U Bus without Measurement M Measured Bus

35 Estimation Model TI u = RHS

36 Complex to Real Conversion [ T ][ I ] = [ RHS] = [ I Y V ] u m red m [ X][ β ] = [ y]

37 Least Angle Regression and Shrinkage (LARS) [*] [*] R. Tibshirani, Regression shrinkage and selection via the lasso, Journal of the Royal Statistical Society. Series B (Methodological), pp , B. Efron, T. Hastie, I. Johnstone, R. Tibshirani et al., Least angle regression, The Annals of statistics, vol. 32, no. 2, pp , 2004.

38 Uniqueness of Solution Define bus groups (column indices) for line terminals: g 1, g 2...g L Then the uniqueness condition will be: rank{ T } = full for i = 1, 2,, L g i That is: columns of each g i must be linearly independent.

39 Measurement Placement Flowchart

40 Positive Sequence Model 118 buses 41 voltage phasor measurements 131 fault cases required the use of LARS T I = [(100 + j200) (300 + j600)] pu.. Results Summary for Noise-Free Cases Ave. Running Time Ave. No. of Iterations Ave. Error s s 5.76 e-14

41 Positive Sequence Model

42 Three Phase Model 30-bus three-phase system All five types of faults (SLG, LL, DLG, 3P, 3PG) are simulated 16 voltage phasor measurements 23 line faults of all 5 types required the use of LARS

43 Three-Phase Model

44 Remarks and Conclusions PMUs can be used to track network changes (line outages / faults) without the need for full network observability. Measurement noise remains an issue which requires further work. Work is currently being extended to detection of faults in T-circuits.

45 Acknowledgements National Science Foundation/DoE NSF/ERC CURENT Thank You Collaborators: Any Questions? Guangyu Feng, Graduate student Roozbeh Emami, Ph.D. 2011