ECEN 667 Power System Stability Lecture 20: Oscillations, Small Signal Stability Analysis

Transcription

1 ECEN 667 Power System Stability Lecture 20: Oscillations, Small Signal Stability Analysis Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, 1

2 Announcements Read Chapter 7 Homework 6 is due today Final is as per TAMU schedule. That is, Friday Dec 8 from 3 to 5pm 2

3 Oscillations An oscillation is just a repetitive motion that can be either undamped, positively damped (decaying with time) or negatively damped (growing with time) If the oscillation can be written as a sinusoid then t cos sin cos t e a t b t e C t 2 2 b where C A B and tan a And the damping ratio is defined as (see Kundur 12.46) 2 2 The percent damping is just the damping ratio multiplied by 100; goal is sufficiently positive damping 3

4 Power System Oscillations Power systems can experience a wide range of oscillations, ranging from highly damped and high frequency switching transients to sustained low frequency (< 2 Hz) inter-area oscillations affecting an entire interconnect Types of oscillations include Transients: Usually high frequency and highly damped Local plant: Usually from 1 to 5 Hz Inter-area oscillations: From 0.15 to 1 Hz Slower dynamics: Such as AGC, less than 0.15 Hz Subsynchronous resonance: 10 to 50 Hz (less than synchronous) 4

5 Example Oscillations The below graph shows an oscillation that was observed during a 1996 WECC Blackout 5

6 Example Oscillations The below graph shows oscillations on the Michigan/Ontario Interface on 8/14/03 6

7 Fictitious System Oscillation Movie shows an example of sustained oscillations in an equivalent system 7

8 Forced Oscillations in WECC (from [1]) Summer hour data: 0.37 Hz oscillations observed for several hours. Confirmed to be forced oscillations at a hydro plant from vortex effect data: Another 0.5 Hz oscillation also observed. Source points to hydro unit as well. And 0.7 Hz. And 1.12 Hz. And 2 Hz. Resonance possible when system mode poorly damped and close. Resonance observed in model simulations. 1. M. Venkatasubramanian, Oscillation Monitoring System, June

9 Forced Oscillations in WECC (from [1]) Summer hour data: 0.37 Hz oscillations observed for several hours. Confirmed to be forced oscillations at a hydro plant from vortex effect data: Another 0.5 Hz oscillation also observed. Source points to hydro unit as well. And 0.7 Hz. And 1.12 Hz. And 2 Hz. Resonance possible when system mode poorly damped and close. Resonance observed in model simulations. 1. M. Venkatasubramanian, Oscillation Monitoring System, June

10 Observing Modes and Damping With the advent of wide-scale PMU deployments, the modes and damping can be observed two ways Event (ringdown) analysis this requires an event Ambient noise analysis always available, but not as distinct Image Source: M. Venkatasubramanian, Oscillation Monitoring System, June

11 Resonance with Interarea Mode [1] Resonance effect high when: Forced oscillation frequency near system mode frequency System mode poorly damped Forced oscillation location near the two distant ends of mode Resonance effect medium when Some conditions hold Resonance effect small when None of the conditions holds 1. M. Venkatasubramanian, Oscillation Monitoring System, June

12 Medium Resonance on 11/29/ MW 0.26 Hz Forced Oscillation in Alberta Canada 200 MW Oscillations on California-Oregon Inter-tie System mode 0.27 Hz at 8% damping Two out of the three conditions were true. 1. M. Venkatasubramanian, Oscillation Monitoring System, June

13 An On-line Oscillation Detection Tool Image source: WECC Joint Synchronized Information Subcommittee Report, October

14 Small Signal Stability Analysis Small signal stability is the ability of the power system to maintain synchronism following a small disturbance System is continually subject to small disturbances, such as changes in the load The operating equilibrium point (EP) obviously must be stable Small system stability analysis (SSA) is studied to get a feel for how close the system is to losing stability and to get additional insight into the system response There must be positive damping 14

15 Model Based SSA Assume the power system is modeled in our standard form as x f x, y 0 = g(x,y) The system can be linearized about an equilibrium point Δx = AΔx BΔy Eliminating Dy gives 0 = CΔx + DΔy -1 Δx = A BD C Δx A Δx If there are just classical generator models then D is the power flow Jacobian;otherwise it also includes the stator algebraic equations sys 15

16 Model Based SSA The matrix A sys can be calculated doing a partial factorization, just like what was done with Kron reduction SSA is done by looking at the eigenvalues (and other properties) of A sys 16

17 SSA Two Generator Example Consider the two bus, two classical generator system from lectures 18 and 20 with X d1 '=0.3, H 1 =3.0, X d2 '=0.2, H 2 =6.0 GENCLS Bus 1 X=0.22 Bus 2 GENCLS slack Deg pu 0.00 Deg pu Essentially everything needed to calculate the A, B, C and D matrices was covered in lecture 19 17

18 SSA Two Generator Example The A matrix is calculated differentiating f(x,y) with respect to x (where x is d 1, D 1, d 2, D 2 ) dd1 D1. pus dt dd1, pu 1 dt 2H P P D D 1 2 M 1 E pu dd2 D2. pus dt dd2, pu 1 dt 2H P P D D P M 2 E pu 2 2 E E V G E E V G E V E V B Ei Di Di Di i Qi Qi Qi i Di Qi Qi Di i E je E cosd j sind Di Qi i i i 18

19 Giving SSA Two Generator Example A B, C and D are as calculated previously for the implicit integration, except the elements in B are not multiplied by Dt/ B

20 SSA Two Generator Example The C and D matrices are C, D Giving A sys A - BD C

21 SSA Two Generator Calculating the eigenvalues gives a complex pair and two zero eigenvalues The complex pair, with values of +/- j11.39 corresponds to the generators oscillating against each other at 1.81 Hz One of the zero eigenvalues corresponds to the lack of an angle reference Could be rectified by redefining angles to be with respect to a reference angle (see book 226) or we just live with the zero Other zero is associated with lack of speed dependence in the generator torques 21

22 SSA Two Generator Speeds The two generator system response is shown below for a small disturbance Notice the actual response closely matches the calculated frequency gfedcb Speed, Gen Bus 1 #1 gfedcb Speed, Gen Bus 2 #1 22

23 SSA Three Generator Example The two generator system is extended to three generators with the third generator having H 3 of 8 and X d3 '=0.3 GENCLS Bus 1 Bus 2 X=0.2 GENCLS slack 3.53 Deg pu X=0.2 X= Deg pu Bus 3 GENCLS Deg pu 200 MW 0 Mvar 23

24 SSA Three Generator Example Using SSA, two frequencies are identified: one at 2.02 Hz and one at 1.51 Hz The oscillation is started with a short, self-clearing fault Shortly we ll discuss modal analysis to determine the contribution of each mode to each signal PowerWorld case B2_CLS_3Gen_SSA 24

25 Visualizing the Oscillations with PowerWorld Visualization of results can be key to understanding and explaining power system dynamics The PowerWorld transient stability contour toolbar allows for the rapid creation of a time-sequenced oneline contours of transient stability results These displays can then be made into a movie by either Capturing the screen as the contours are creating using screen recording software such as Camtasia Or having Simulator automatically store the contour images as jpegs and then creating a movie using software such as Microsoft Movie Maker 25

26 Contour Toolbar The contour toolbar uses stored transient stability results, and hence it is used only after the transient stability solution has finished It can be used with either RAM or hard drive results It requires having a oneline with objects associated with the desired results (buses, generators, substations, etc) The contour toolbar is shown by either Selecting Add-ons, Stability Case Info, Show Transient Contour Toolbar On the Transient Stability Analysis Form select the Show Transient Contour Toolbar button at the bottom of the form 26

27 Contour Toolbar Buttons and Fields Several buttons and fields control the creation of the contour images The Options menu is used to specify three options Contour Options is used to display the Contour Options dialog; it must be used to specify the contouring options, including the maximum/nominal/minimum values When the Contour Options dialog is shown the Contour Type and Value fields are disabled since these values are specified on the toolbar Other options can be set to customize the contour 27

28 Contour Toolbar Contour Options These fields are disabled The Maximum, Nominal and Minimum values need to be set. These should be set taking into account how the values vary throughout the transient stability run. 28

29 Contour Toolbar Export Option Export is used to display the TS Contour Export Options dialog, which provides the option of saving the contours as jpeg format files On the dialog select Save Images as JPEGs to save the files Each jpeg file will have a base file name with the associated time in seconds appended 29

30 Contour Toolbar Value Meaning Option The Value Meaning option is used to indicate which values will actually be contoured Actual Value: Good for things like frequency Percent of Initial: Contours the percent of the initial value; useful for values with widely different initial values, like generator MW outputs Deviation from Initial: Contours the deviation from the initial; good for voltages Percent Deviation: Percentage deviation Contour limits should be set appropriately 30

31 Creating a Contour Sequence A typical way to use the contour tool bar is to Solved the desired transient stability contingency, making sure to save the associated contour fields Display the toolbar Set the options Select the play button which will automatically create the sequence of contours Contours can be saved with either screen recording software or as a series of jpegs Images can be combined as a movie using a tool such as Microsoft Movie Maker 31

32 Three Bus Example (With H 1 set to 6.0 and H 3 to 4.0) GENCLS Transient Stability Time (Sec): Bus 1 Bus 2 X=0.2 GENCLS slack 3.53 Deg pu X=0.2 X= Deg pu Bus 3 GENCLS Deg pu 200 MW 0 Mvar A low resolution copy of the movie is on the website 32

33 Comtrade Format (IEEE Std. C37.111) Comtrade is a standard for exchanging power system time-varying data Originally developed for power system transient results such as from digital fault recorders (DFRs), but it can be used for any data Comtrade is now being used for the exchange of PMU data and transient stability results Three variations on the standard (1991, 1999 and 2013 format) PowerWorld now allows transient stability results to be quickly saved in all three Comtrade Formats 33

34 Three Bus Example Comtrade Results The 1991 format is just ascii using four files; the 1999 format extends to allow data to be stored in binary format; the 2013 format extends to allow a single file format 34

35 Large System Studies The challenge with large systems, which could have more than 100,000 states, is the shear size Most eigenvalues are associated with the local plants Computing all the eigenvalues is computationally challenging, order n 3 Specialized approaches can be used to calculate particular eigenvalues of large matrices See Kundur, Section 12.8 and associated references 35

36 Single Machine Infinite Bus A quite useful analysis technique is to consider the small signal stability associated with a single generator connected to the rest of the system through an equivalent transmission line Driving point impedance looking into the system is used to calculate the equivalent line's impedance The Z ii value can be calculated quite quickly using sparse vector methods Rest of the system is assumed to be an infinite bus with its voltage set to match the generator's real and reactive power injection and voltage 36

37 Small SMIB Example As a small example, consider the 4 bus system shown below, in which bus 2 really is an infinite bus GENCLS Bus 4 X=0.1 Bus 1 Bus 2 X=0.2 Infinite Bus slack Deg pu X=0.1 Bus 3 X= Deg pu 4.46 Deg pu 0.00 Deg pu To get the SMIB for bus 4, first calculate Z Ybus j Z44 j Z 44 is Z th in parallel with jx' d,4 (which is j0.3) so Z th is j

38 Small SMIB Example The infinite bus voltage is then calculated so as to match the bus i terminal voltage and current V V Z I inf where i i i P i jq V In the example we have inf inf i i * * * P4 jq 4 1 j j V j V 1. 0 I i.. (. ). V j0 220 j j0 328 While this was demonstrated on an extremely small system for clarity, the approach works the same for any size system 38

39 Calculating the A Matrix The SMIB model A matrix can then be calculated either analytically or numerically The equivalent line's impedance can be embedded in the generator model so the infinite bus looks like the "terminal" This matrix is calculated in PowerWorld by selecting Transient Stability, SMIB Eigenvalues Select Run SMIB to perform an SMIB analysis for all the generators in a case Right click on a generator on the SMIB form and select Show SMIB to see the Generator SMIB Eigenvalue Dialog These two bus equivalent networks can also be saved, which can be quite useful for understanding the behavior of individual generators 39

40 Example: Bus 4 SMIB Dialog On the SMIB dialog, the General Information tab shows information about the two bus equivalent PowerWorld case B4_SMIB 40

41 Example: Bus 4 SMIB Dialog On the SMIB dialog, the A Matrix tab shows the A sys matrix for the SMIB generator In this example A 21 is showing D4, pu 1 PE4, cos d4 2H4 d

42 Example: Bus 4 SMIB Dialog On the SMIB dialog, the Eigenvalues tab shows the A sys matrix eigenvalues and participation factors (which we'll cover shortly) Saving the two bus SMIB equivalent, and putting a short, self-cleared fault at the terminal shows the 1.89 Hz, undamped response 42