arxiv: v1 [physics.soc-ph] 11 Dec 2017

Transcription

1 Cascadng Faures as Contnuous Phase-Space Transtons Yang Yang and Adson E Motter, Department of Physcs and Astronomy, Northwestern Unversty, Evanston, IL 68, USA Northwestern Insttute on Compex Systems, Northwestern Unversty, Evanston, IL 68, USA arxv:7453v [physcssoc-ph] Dec 7 In network systems, a oca perturbaton can ampfy as t propagates, potentay eadng to a argescae cascadng faure Here we derve a contnuous mode to advance our understandng of cascadng faures n power-grd networks The mode accounts for both the faure of transmsson nes and the desynchronzaton of power generators, and ncorporates the transent dynamcs between successve steps of the cascade In ths framework, we show that a cascade event s a phase-space transton from an equbrum state wth hgh energy to an equbrum state wth ower energy, whch can be sutaby descrbed n cosed form usng a goba Hamtonan-ke functon From ths functon we show that a perturbed system cannot aways reach the equbrum state predcted by quas-steadystate cascade modes, whch woud correspond to a reduced number of faures, and may nstead undergo a arger cascade We aso show that n the presence of two or more perturbatons, the outcome depends strongy on the order and tmng of the ndvdua perturbatons These resuts offer new nsghts nto the current understandng of cascadng dynamcs, wth potenta mpcatons for contro nterventons Cascadng processes undere a myrad of network phenomena [], ncudng backouts n power systems [, 3], secondary extnctons n ecosystems [4, 5], and compex contagon n fnanca networks [6, 7] In a such cases, an otherwse sma perturbaton may propagate and eventuay cause a szabe porton of the system to fa Varous system-ndependent cascade modes have been proposed [8 3] and used to draw genera concusons, such as on the mpact of nterdependences [4] and countermeasures [5] There are outstandng questons, however, for whch t s necessary to mode the cascade dynamcs startng from the actua dynamca state of the system In power-grd networks, the state of the system s determned by the power fow over transmsson nes and the frequency of the power generators, whch must be respectvey beow capacty and synchronzed under norma steady-state condtons Athough a oca perturbaton has a mted mpact on the connectvty of the network, t may trgger a cascade of faures and protectve responses that swtch off grd components and may aso ead generators to ose synchrony Much of our current understandng about ths process has been derved from quas-steady-state cascade modes [6 ], whch use teratve procedures to mode the successve nactvaton of network components caused by power fow redstrbutons, whe omttng the transent dynamcs between steady states as we as the dynamcs of the generators Further understandng has resuted from stabty studes focused on the synchronzaton dynamcs of power generators n the absence of fow redstrbutons [ 6] Yet, to date no theoretca approach has been deveoped to ncorporate at the same tme these two fundamenta aspects of power-grd dynamcs frequency change and fow redstrbuton n the modeng of cascadng faures [7] The goa of our study s to f ths gap and consder the nteracton between these two factors Our framework s nspred by energy functon anayss approaches consdered n the study of power system stabty [8, 9] and of bstabty of crcut eements [3] Specfcay, n ths Letter we ntroduce a tmecontnuous cascade mode that ncudes the dynamcs of the state varabes governed by the swng equatons of the generators, frequency dependence of oads, and power fow equatons as we as the dynamcs of the status varabes descrbng the on/off (e, operatona/dsabed) condton of the transmsson nes Wthn ths mode, the steady operatng states of the system correspond to stabe equbra, and a cascade event s a phase-space transton from one stabe equbrum to another We study these states and show that the stabe equbra are the oca mnma of an energyke functon Ths eads to numerous mportant mpcatons that have not been systematcay studed before In partcuar, t foows from the propertes of ths functon that a perturbed system cannot aways reach the equbrum state predcted by quas-steady-state modes, and may nstead approach an equbrum correspondng to a arger cascade; ths hghghts the mportance of the dynamcs between successve steps of a cascade, as consdered n our contnuous mode, whch s a factor that has remaned unexpored wth few exceptons [, 3 33] It aso foows that the equbrum energy does not depend monotoncay on the number of faures, and that cascades trggered by mutpe perturbatons depend strongy on the perturbaton order These resuts suggest the possbty of cascade mtgaton usng judcousy desgned perturbatons to steer the system to a preferred equbrum that woud not be reached spontaneousy We frst consder the protectve operaton, common to most power networks, that removes a transmsson ne when the fow on t exceeds ts capacty We assocate each ne wth a contnuous varabe η representng ts on/off status (as we as the contnuous process of swtchng between the two condtons) and a parameter

2 (a) (b) 4 6 f(`) - -4 (f) ` 5 ` ` (c) ` (n) ` 6``! F (`) 4 ` = ` = (f) ` 5 ` ` (c) ` = (n) ` Fgure Lne-status swtch mode (a) Functon f(η ) for a =, whose roots are the equbrum ponts of Eq () when λ = (for other vaues of λ, see Fg S n Suppementa Matera [34]) (b) Potenta functon φ(η ) = λ η F (η ), where the oca mnma for λ = correspond to the stabe equbra n (a) When λ s ncreased past, the oca mnmum η (n) merges wth η (c) and then dsappears λ ndcatng the fracton of the ne capacty used by the fow As shown beow, ths aows us to ncorporate the ne status nto the dynamca equatons by scang the power fow terms by η, wth η representng the norma status for λ < and the faed status for λ, where η s thus constraned to the unt nterva To mode the automatc remova of the overoaded nes, we can then defne the dynamcs of η as (c) ` (n) ` η = f(η ) λ, () where the rhs s defned to satsfy three physca condtons: (I) for λ <, there are three equbra η (f) < η (c) < η (n), where η (n) s a stabe equbrum representng the norma operaton status, η (f) s a stabe equbrum representng the faed status, and η (c) s an unstabe equbrum markng the crtca vaue beow whch η evoves to the faed status; (II) for λ, ony the equbrum η (f) remans stabe, whch s satsfed f the oca maxmum of f n (η (c), η (n) ) s ; (III) η (c) s aways cose to, snce a ne shoud be fuy operatona under norma condtons The dynamcs does not depend senstvey on the detas of functon f provded these condtons are satsfed Throughout, we use overdot to ndcate tme dervatve Here we defne f(η ) = a [ η ( η ) ] +aη 4 b, where a and b are postve parameters The terms η and ( η ) constran η above and beow, respectvey, as they ensure that f(η ) for η + and f(η ) for η The term η 4 aows f to have three roots correspondng to η (f), η (c), and η (n) for λ =, as shown n Fg (a) The parameters a and b are adjustabe to set η (f) cose to, to set η (c) and η (n) suffcenty cose to, and to set the oca maxmum of f to For ths choce of functon f, Eq () satsfes condtons (I)-(III) Moreover, the equaton can be rewrtten as a gradent system η = dφ(η )/dη, where φ(η ) = λ η F (η ), and df (η )/dη = f(η ) As shown n Fg (b), the stabe equbra of ths system correspond to the oca mnma of φ(η ) Foowng a perturbaton, the power fowng on transmsson nes can change dynamcay When the fow on ne reaches ts capacty (λ ), the system w experence a sadde-node bfurcaton and the status varabe η w evove to the stabe equbrum η (f), representng a ne swtch-off operaton Ths s a one-way acton, snce the equbrum η (f) s stabe for any vaue of λ Havng defned the dynamcs of the status varabes, we now ncorporate the system s protectve response nto the dynamca equatons governng the state of the network In a network of n nongenerator nodes, each such node s an eectrc pont where power s extracted by a oad, receved from generators, and/or redstrbuted among transmsson nes We denote by n g the number of generators, and by n the number of transmsson nes To proceed, we consder the extended representaton of the network [37] n whch each generator s now an addtona node connected to the network through a vrtua ne (not ncuded n n and not subject to faure), eadng to a network of n+n g nodes For notatona convenence, we rendex the generators as the frst n g nodes Assumng that the votage satsfes V (n per unt) for a nodes and that no rea power s ost on transmsson nes, we can defne the state of a power system as x = (ω, δ, η) Here, ω = (ω ) are the frequences of the generators reatve to the system s nomna frequency, δ = (δ ) are the phase anges of a other nodes reatve to a reference node (taken to be =, so that δ ), and η = (η ) are the status varabes of the (nonvrtua) transmsson nes L, where L The state of the system s sutaby determned by the foowng equatons: ω = D ω [ P + n g+n M M j=n g+ j sn δ j ], =,,, n g, δ = ω ω, =,, n g, δ = [P + n g T B j= j sn δ j + n g+n B ] j=n g+ j η -j sn δ j ω, = n g +,, n g + n, [ η -j = f(η -j ) B j ( cos δ j ) ], -j L W -j ()

3 3 x = J Ψ(x), (3) where Ψ(x) s an energy functon defned as Ψ(x) = Png +n e Png h M ω + j=n Bj ( cos δj ) = g + Png +n Png +n e + =ng + j=+ Bj ( cos δj )η`-j Png +n P + = P δ `-j L W`-j F (η`-j ), (4) and J s a matrx of the form J J J3 T J J = J T J33 3 J44 (5) In ths matrx, the off-dagona bocks are M M J = M M3 M Mng M, J3 = M, (6) (b) - - x (d) `` ` Fgure Smuated cascade event n Iceand s power grd (a) Dagram of the network, whch conssts of 35 generators ( ), 89 nongenerator nodes ( ), and 3 transmsson nes ( ) [4] The remova of the marked ne ( ) trggers a sequence of 6 subsequent ne faures (magenta) that separate the network nto 5 custers (coor coded) (b-c) Correspondng generator frequences ω (b) and ne-status varabes η` (c) as functons of tme [coor coded as n (a)] (d) Correspondng fracton λ` of the ne capacty used, shoud the ne overoaded at s (arrow) not be dsabed and the dagona bocks are Mn, M,, D g ), J = dag( M D D ng J33 = dag( Tn+, Tn+,, Tn +n ), g J44 = g dag( W, W, g, Wn 5 T T (c) -5 6` `` /! /! `!! (a) ` e s a symmetrc matrx wth Here, δj = δ δj and B e nonzero eements Bj = /x`-j, where x`-j s the transent reactance of a generator or s the reactance of a transmsson ne, dependng on whether the ne connectng and j s vrtua or not The frst two equatons are the swng equatons descrbng the dynamcs of the generators, where M s the generator rotor nerta, D s the rotor dampng rato, and P s the negatve of the (m) mechanca power nput P of the generator [38] The thrd equaton descrbes oads (and nongenerator nodes n genera, under the assumpton that they ncude some frequency-dependent power exchange) as frst-order rotors, where T s the oad frequency rato and P s the (d) power P demanded at the node We further assume Png +n that = P =, so that there exsts an equbrum pont at ω = and δ = cte Note that the term representng the power fow on ne `-j s mutped by the status varabe η`-j, whch automatcay turns off the (f) ne n the event of an overoad (when η` η` ) The ast equaton descrbes the dynamcs of the status varej ( cos δj ), abes, where λ`-j n Eq () s repaced by B the reactance energy stored n the transmsson ne `-j, dvded by W`-j, the maxmum reactance energy that ne `-j can hod The prefactor n ths equaton assures that the tme scae for ne faures s much shorter than that of the other dynamca changes n the network For more detas on the dervaton of Eq (), see Suppementa Matera [34] Importanty, we can show that Eq () can be derved from a Hamtonan-ke system of the form /! //! / (7) ) For detas on the dervaton of Eq (3), see Suppementa Matera [34] Crucay, the matrx J s the sum of a skew-symmetrc matrx and a dagona matrx wth nonpostve eements, from whch we can show that dψ(x(t))/dt = Ψ(x)T x = Ψ(x)T J Ψ(x) Moreover, because J s aso fu rank (whch foows from ts reduced row echeon form), we have that dψ(x(t))/dt = f and ony f Ψ(x) =, and hence f and ony f x = J Ψ(x) = Thus, when the network s perturbed, the energy functon Ψ(x) monotoncay decreases as the system evoves, and becomes constant agan ony when the system reaches an equbrum pont of Eq (3) [and hence of Eq ()] Such equbra represent stabe steady states, where the generators are synchronzed [ω (t) = ω (t) = = ωng (t)], the ange dfferences are fxed for a pars of nodes, and the fow s beow capacty for a operatng transmsson nes We frst ustrate our formasm on Iceand s powergrd network, shown n Fg (a) (for parameter settng, see Suppementa Matera [34]) The system s desgned to have a stabe steady state wth no addtona faures when any snge transmsson ne s mssng (provded

4 4 the network remans connected), whch s verfed n our smuatons We test whether such a cascade-free steady state s actuay reached foowng the remova of a ne when the transent dynamcs between steady states represented n our mode s taken nto consderaton Startng from the stabe steady state determned by Eq (), we smuate a 68 snge-ne remova perturbatons that keep the network topoogcay connected (performed by changng η to η (f) ) Of these, do not converge to the best avaabe stabe steady state and nstead undergo subsequent faures (Fg S3 n Suppementa Matera [34]) Insghts nto the underyng mechansm are provded by the exampe shown n Fgs (a)-(c), where a sequence of ne overoads separates the network nto 5 custers As shown n Fg (d), the system woud eventuay have approached the desgned steady state wth no addtona faures, but a ne overoad whose automatc swtch-off trggers subsequent overoads occurs before the system can reach that state In ths case, no feasbe trajectory exsts n the phase space connectng the nta state to the steady state predcted by quassteady-state modes Ths scenaro s common n genera, as shown for fve other systems n the 3rd coumn of Tabe S (Suppementa Matera [34]) When the network s subject to mutpe perturbatons, our framework shows that the cascade outcome w generay depend on the order and tmng of the perturbatons A natura measure to quantfy ths dfference s the sze C (e, number of nodes) of the argest connected custer n the postcascade stabe state As an ustraton, we consder the foowng three scenaros for two-ne remova perturbatons: () remove ne -j and then, after the stabe state s reached, remove ne -j ; () same as n () but for -j swapped wth -j ; () remove -j and -j concurrenty Consderng a, 7 pars of nes ( -j, -j ) that keep Iceand s network connected after ther remova (but not necessary after the resutng cascadng faures), our smuatons ndcate that 3% of these perturbatons ead to cascades n at east one of the scenaros above For ths subset of ne pars, we obtan that: (a) order matters n 79% of the cases, n that C dffers for at east one of the scenaros; (b) choosng between the orders n () and () eads to the argest C n 8% of the cases; (c) () and () ead to equay best C n 43% of the cases; (d) the concurrent remova scenaro () trumps () and () n the remanng 8% of the cases (for specfc exampes, see Fgs S4 and S5 n Suppementa Matera [34]) Smar trends are observed for a fve other systems consdered, as shown n Tabe S (Suppementa Matera [34]) Ths order dependence has potenta mpcatons for contro, as t can be expoted n proactve ne removas to prevent subsequent faures (Fg S6 n Suppementa Matera [34]) Ths reveas a sharp contrast between processes for whch order s mmatera, such as percoaton, and the cascades consdered here Energy, *(x) (a) 99 unperturbed state unperturbed state Energy (b) generators non-generators Fgure 3 Energy eves Ψ(x) of the stabe states n the 4- bus test system Each pane shows a combnatons of (eft coumn) to 7 (rght) successve ne removas that eave the network connected (a) A stabe states wthout addtona faures determned usng the MATLAB functon fsove (b) Subset of stabe states n (a) that the system actuay evoves to for the same ne removas as n (a) Aso marked are the fractons of perturbatons for whch a stabe state s dentfed (a) and the fractons of those stabe states actuay reached (b) The dagram on the rhs shows the topoogy of the network Takng the anayss one step further, our formasm offers unque nsght nto the reaton between ne remova perturbatons and energy eves Fgure 3(a) shows a energy eves for stabe steady states of the IEEE 4-bus test system (chosen n pace of Iceand s network to avod a cuttered pcture) for a combnatons of to 7 ne removas that keep the network connected Fgure 3(b) shows the states that the system actuay approaches foowng these successve ne removas the mssng states [compared to Fg 3(a)] are the ones not reached because the system undergoes a cascade Two major resuts foow from ths Frst, t confrms that upon perturbaton the system often does not reach the avaabe stabe steady state wth smaest number of faures (eg, for 7 ne removas, ths s so for 98% of a cases) Second, the range of energy eves wth k + ne removas overaps wth the range for k ne removas There are, for exampe, stabe steady states wth ony one ne faure at ower energy than many stabe states wth, 3,, 6 ne faures Ths shows that, foowng a perturbaton that coud eventuay ead to a stabe state wth mutpe faures, the system can n prncpe be steered to a ower-energy state whch has, nevertheess, a reduced number of faures Crucay, ths s possbe wthout an externa nput of energy, as the system tends to go spontaneousy to ower-energy states foowng a perturbaton In summary, the mode presented here accounts n a snge phase space for the nteracton between the fu dynamcs of a cascade (ncudng transents) and the changes to the underyng network structure mposed by the resutng faures The resuts expan the combnatora mpact of perturbatons, dentfy condtons under whch a cascade may deveop despte the presence of a stabe state that woud wthstand the perturbaton, and

5 5 suggest new opportuntes for cascade contro The authors thank Takash Nshkawa for feedback on the manuscrpt Ths work was supported by an ISEN Booster Award, NSF Grant DMS-578, Smons Foundaton Award 3496, and ARPA-E Award de-ar7 The vews and opnons of authors expressed heren do not necessary state or refect those of the Unted States Government or any agency thereof [] A E Motter and Y Yang, Physcs Today 7(), 3 (7) [] R Lscousk, et a, Fna Report on the August 4, 3 Backout n the Unted States and Canada: Causes and Recommendatons (4) [3] F Vandenberghe, et a, Fna Report of the Investgaton Commttee on the 8 September 3 Backout n Itay (4) [4] J A Dunne and R J Wams, Phos T R Soc B 364, 7 (9) [5] S Sahasrabudhe and A E Motter, Nature Commun, 7 () [6] P Ga and S Kapada, Proc Roy Soc A 466, 4 () [7] M Eott, B Goub, and M O Jackson, Am Econ Rev 4, 35 (4) [8] D J Watts, Proc Nat Acad Sc USA 99, 5766 () [9] A E Motter and Y-C La, Phys Rev E 66, 65(R) () [] K-I Goh, D-S Lee, B Kahng, and D Km, Phys Rev Lett 9, 487 (3) [] P Cructt, V Latora, and M Marchor, Phys Rev E 69, 454 (4) [] S V Budyrev, R Parshan, G Pau, H E Staney, and S Havn, Nature 464, 5 () [3] C D Brummtt, G Barnett, and R M D Souza, J R Soc Interface, 57 (5) [4] C D Brummtt, R M D Souza, and E A Lecht, Proc Nat Acad Sc USA 9, E68 () [5] A E Motter, Phys Rev Lett 93, 987 (4) [6] I Dobson, B A Carreras, V E Lynch, and D E Newman, Chaos 7, 63 (7) [7] M Anghe, K A Werey, and A E Motter, Proc 4th Ann HICSS (IEEE, New York, 7), Vo, p 3 [8] D Watts and H Ren, Proc 8 IEEE Int Conf Sustan Energy Tech (IEEE, New York, 8), p [9] R D Zmmerman, C E Muro-Snchez, and R J Thomas, IEEE T Power Syst 6, () [] D Wtthaut, M Rohden, X Zhang, S Haerberg, and M Tmme, Phys Rev Lett 6, 387 (6) [] A Moussaw, N Derzsy, X Ln, B K Szymansk, and G Kornss, Sc Rep 7, 79 (7) [] Y Susuk, I Mezć, and T Hkhara, J Nonn Sc, 43 () [3] M Rohden, A Sorge, M Tmme, and D Wtthaut, Phys Rev Lett 9, 64 () [4] A E Motter, S A Myers, M Anghe, and T Nshkawa, Nature Phys 9, 9 (3) [5] F Dörfer, M Chertkov, and F Buo, Proc Nat Acad Sc USA, 5 (3) [6] P J Menck, J Hetzg, J Kurths, and H J Schenhuber, Nature Commun 5, 3969 (4) [7] Theoretca understandng s especay needed gven that avaabe reports of emprca data on cascade backouts generay do not ncude detaed nformaton about generator behavor aong wth power fow nformaton [8] A R Bergen and D J H, IEEE T Power Ap Syst PAS-, 5 (98) [9] A Pa, Energy Functon Anayss for Power System Stabty (Sprnger, ) [3] C L DeMarco, IEEE Contr Syst Mag, 4 () [3] I Smonsen, L Buzna, K Peters, S Bornhodt, and D Hebng, Phys Rev Lett, 87 (8) [3] SP Corneus, WL Kath, and AE Motter, Nature Commun 4, 94 (3) [33] B Schäfer, D Wtthaut, M Tmme, and V Latora, arxv:7788 [nnao] (7) [34] See Suppementa Matera [ur], whch ncudes Refs [35, 36], for detas of the anayss and addtona exampes [35] H Zheng and C L DeMarco, Proc North Amercan Power Symposum, (IEEE, New York, ), p [36] F Mano, An open source power system anayss toobox, IEEE T Power Syst, 99 (5) [37] T Nshkawa and A E Motter, New J Phys 7, 5 (5) [38] These equatons foow from Newton s second aw apped T (e) to the generator rotor [39]: I ω = D ( ω + T (m) where I s the moment of nerta, D s the dampng coeffcent, T (m) s the mechanca torque, and T (e) s the torque due to eectrca oad n the network Thus, M = = [(ω + ω o)/p o]t (m) (ω o/p o)i, D = (ω o/p o) D, P (m) and n g+n B j=n g+ j sn δ j = [(ω + ω o)/p o]t (e) ),,, where the reatve frequency ω s assumed to be sma compared to the nomna frequency ω o, and the base power P o s used to transform the power terms nto per unt quanttes [39] J J Granger and W D Stevenson, Power System Anayss (McGraw-H, 994) [4] Iceand s transmsson network, wwwmathsedacuk/opt energy/networkdata/howtousehtm (accessed: )

6 Suppementa Matera Cascadng Faures as Contnuous Phase-Space Transtons Yang Yang and Adson E Motter DERIVATION OF EQUATION () IN THE MAIN TEXT In a power-grd network, we defne a nongenerator node as a bus and transmsson nes as eectrca connectons, ncudng transformers, between pars of such nodes We denote the number of nongenerator nodes by n, the number of generators by n g, and the number of transmsson nes by n The operatng condton of the network can be characterzed by the compex eectrca power S (e) = P (e) + j Q (e) at each node In steady state, a generators n an aternatng current network run at the same frequency and the dstrbuton of power fows through the network s determned by the compex votage V = V e jδ at each node, where δ s the votage ange n the reference frame Ths s determned through the power fow equatons: P (e) = Q (e) n V V j (G j cos δ j + B j sn δ j ), (S) j= = n V V j (G j sn δ j B j cos δ j ), (S) j= where δ j = δ δ j and Y j = G j + j B j defne a Lapacan-ke matrx In ths matrx, an off-dagona eement Y j s the negatve of the admttance of the ne -j that connects nodes and j Assumng that V pu for a nodes and that there s no rea power ost on the transmsson nes (e, G j = ), we can rewrte the rea power as P (e) = n j= B j sn(δ δ j ), where B s a symmetrc matrx wth off-dagona eements B j = /x and x s the reactance of ne These assumptons are vad throughout ths paper In genera, the state of both generators and oads can change n tme Durng a dsturbance, the generator rotors deceerate or acceerate wth respect to the nomna frequency (6 Hz n the US and 5 Hz n Europe, ncudng dω Iceand) The dynamcs of the generator rotor s governed by the swng equaton: M dt +D ω = P (m) P (e), where ω s the frequency (reatve to the nomna frequency), M s the rotor nerta, D s the rotor dampng rato, and P (m) s the net shaft power nput nto the generator Consderng Eq (S), we can combne the power fow equatons nto the swng equaton as dω M dt + D ω = P (m) n B j sn(δ δ j ) j= (S3) Here we choose the ange of the frst generator as the reference ange, and we defne δ = α α, where dα = ω dt In prncpe, the power consumed by oad nodes coud depend nonneary on the frequency at that node Assumng that the frequency at each node does not devate strongy from the nomna frequency, we can use a nearzed power-frequency reaton to descrbe the dynamcs of the node connected wth a oad: dδ T dt = P (d) n B j sn(δ δ j ) ω, j= (S4) where T s a postve constant and P (d) s the power requested by the oad To account for the nterna reactances of the generators, we adopted an extended representaton of the power grd As expaned n the man text, n ths extended representaton we add n g nodes connected to the network through vrtua ossess nes to represent the n g generators n the system (see Fg S) The reactance of a vrtua ne represents the transent reactance of the correspondng generator Accordngy, we modfy the B matrx n Eqs (S3) and (S4) as [ ] X B = X T, (S5) B

7 where X j s the recproca of the transent reactance of the th generator that connects the jth node n the network Combnng together the dynamcs of generators and oads, the equatons of moton take the foowng form: dω M dt + D ω dδ dt T dδ dt = P (m) n g+n j= B j sn(δ δ j ), =,,, n g, = ω ω, =,, n g, = P (d) n g+n j= B j sn(δ δ j ) ω, = n g +,, n g + n, (S6) where, for smpcty, the ast equaton s assumed to appy to a nongenerator nodes under the assumpton that they a ncude some frequency-dependent power exchange Addng up a the equatons n Eq (S6), at the fxed pont of the dynamcs (where dω /dt =, dδ /dt = ) we obtan ( n g = D + n)ω = n g = P (m) n g+n n P (d) g+ We assume that the rea power s baanced, e, n g = P (m) = n g+n n P (d) g+, and hence a generators operate wth the nomna frequency at the fxed pont To compete the dervaton of Eq (), for each transmsson ne -j we need to defne the fracton λ -j of the ne capacty used by the fow Ths quantty s determned by the average amount of reactve power stored n the transmsson ne, cacuated as (V V j )I j, where I j = j B j (V V j ) s the current on the ne Notng that V V j (both approxmatey equa to pu), we can determne ths reactve energy as B j ( cos δ j ), and repace λ -j n Eq () by ths reactve energy dvded by W -j (the maxmum reactve energy that ne -j can hod) By ncorporatng the dynamcs of the status varabes, as presented n the man text, nto the dynamcs of the power system n Eq (S6), we fnay obtan Eq () Ths equaton descrbes, at the same tme, the state varabes of the power system and the status of the transmsson nes PROOF OF EQUATION (3) IN THE MAIN TEXT We now expcty show the equvaence between Eqs (3) and () n the man text We reca that the state of a power system s defned as x = (ω, δ, η) As n the man text, the (reatve) frequences of the generators are represented by a vector ω of sze n g, the phase anges of a nodes reatve to the reference node are represented by a vector δ of sze (n g + n ), and the status varabes of the transmsson nes are represented by a vector η of sze n To factate our verfcaton of Eq (3), we further separate δ nto δ = (δ δ,δ δ ), where the vector δ = (δ,, δ ng ) represents the phase anges of the generator nodes and vector δ = (δ ng+,, δ ng+n) represents the phase anges of nongenerator nodes Accordng to Eq (4), the gradent of Ψ(x) can be decomposed as Ψ(x) = [ ω Ψ(x), δ Ψ(x), δ Ψ(x), ηψ(x) ], where ω Ψ(x) = [ M ω, M ω,, M ng ω ng ] T, (S7) δ δ δ Ψ(x) = δ j= j= j= B j sn δ j + P B 3j sn δ 3j + P 3 B ngj sn δ ngj + P ng, (S8) ng B j= (ng+)j sn δ (ng+)j + n g+n B j=n g+ (ng+)j sn δ (ng+)jη (ng+)-j + P ng+ ng B j= (ng+)j sn δ (ng+)j + n g+n B j=n g+ (ng+)j sn δ (ng+)jη (ng+)-j + P ng+ Ψ(x) =, (S9) ng B j= (ng+n)j sn δ (ng+n)j + n g+n B j=n g+ (ng+n)j sn δ (ng+n)jη (ng+n)-j + P ng+n

8 η Ψ(x) = B j ( cos δ j ) W -j f(η -j ) 3 (S) In the ast equaton, each component of η Ψ(x) corresponds to the dervatve of Ψ(x) wth respect to the status varabe of a transmsson ne -j For notatona convenence, n the expresson of δ Ψ(x) we ntroduce the constants η (ng+)-j assocated wth pars of nodes ( n g +, j ) that are not connected by a transmsson ne (these constants shoud not be confused wth the components of η, whch are varabes assocated wth pars of nodes that are connected) Makng use of the defnton of matrx J n Eqs (5)-(7) of the man text, we obtan D M ω ng M = ( n g+n j= B j sn δ j + P ) D M J ω Ψ(x) + J δ Ψ(x) + J ω M ( n g+n j= B j sn δ j + P ) 3 δ Ψ(x) = + ω ng M ng ( n g+n j= B j sn δ ngj + P ng ) + ng+n M D ng M ng =n ( n g B g+ j= j sn δ j + n g+n B j=n g+ j sn δ j η -j + P ) We note that n g+n B =n g+ j=n g+ j sn δ j η -j = gven that B j = B j and sn δ j = sn δ j Accordngy, the frst row of the rhs of Eq (S) can be rewrtten as (S) D M ω + = D M ω + = D M ω + = D M ω + = D M ω + = D M ω ( n g = j= B j sn δ j + n g+n ng =n M B g+ j= j sn δ j ) + = P M ( n g = j= B j sn δ j + n g+n ng =n M B g+ j= j sn δ j ) P M ( n g = j= B j sn δ j + n g+n ng j=n M B g+ = j sn δ j ) P M ( n g = j= B j sn δ j n g+n ng j=n M B g+ = j sn δ j ) P M ( n g = j= M M j= B j sn δ j n g+n ng B j= = j sn δ j ) P M B j sn δ j M P (S) In the above dervaton, we assume that the tota amount of power njecton nto the network equas the tota amount of power extracted from the network (e, n g+n = P = ), and we use the fact that B j = for, j =,,, n g Smar cacuaton appes to the other rows of Eq (S) Then, notng that the rhs of Eq (S) s the same as the rhs of the frst equaton n (), we concude that J ω Ψ(x) + J δ Ψ(x) + J dω 3 δ Ψ(x) = dt Next, we estabsh the equvaence between the other components of Eq (3) and the other equatons n () We obtan the second equaton n () usng that ω ω J T ω 3 ω ω Ψ(x) = = dδ δ dt, (S3) ω ng ω

9 4 the thrd equaton n () usng that J3 T ω Ψ(x) + J 33 δ Ψ(x) = ω ω + ω T ng+ ( ng j= B (ng+)j sn δ (ng+)j + n g+n j=n g+ B (ng+)j sn δ (ng+)jη (ng+)-j + P ng+ T ng+ ( n g j= B (ng+)j sn δ (ng+)j + n g+n j=n g+ B (ng+)j sn δ (ng+)jη (ng+)-j + P ng+) T ng+n ( n g j= B (ng+n)j sn δ (ng+n)j + n g+n = dδ δ dt, ) B j=n g+ (ng+n)j sn δ (ng+n)jη (ng+n)-j + P ng+n) (S4) and the fourth equaton n () usng that J 44 η Ψ(x) = B j ( cos δ j ) + f(η -j ) W -j = dη dt (S5) Combnng Eqs (S)-(S5), we have proved Eq (3) n the man text We note that a smar resut on a 3-bus network s reported n Ref [35] However, that study does not offer a framework to address a network wth an arbtrary number of buses, whch s derved here For competeness, we note that our energy-functon formuaton s aso naturay suted for stabty anayss Rewrtng Eq (3) as ẋ = g(x), f x s an equbrum state, by defnton we have g(x ) = The stabty of ths state s then determned by the egenvaues of the Jacoban matrx dg/dx x=x = JH(x ), where the components of matrx H are gven by H kk = Ψ x k x k GENERALIZED HAMILTONIAN-LIKE ENERGY WHEN THE NETWORK SPLITS Our anayss of the energy functon Ψ(x) assumes a baance of power n the network, whch s guaranteed when the network remans connected but can be voated when t spts durng a cascade Ths s easy remedated, however, by extendng the formasm to ntroduce a reference generator n each custer to mmc the system s operaton of rebaancng rea power As shown beow, the reference generator n each custer then serves as a sack bus that prevents mbaances between power nput and output, whch woud cause acceeraton or deceeraton of generators In our appcatons n the paper, the network remans connected or ese we consder each custer separatey We frst note that the choce of the reference generator does not mpact our descrpton of power-grd dynamcs, as ong as the power nput and output at each node s fxed Then we can show that n a power network satsfyng ng+n = P =, the vaue of the energyke functon Ψ(x) remans the same for dfferent choces of the reference node Specfcay, when the ange of generator s chosen to be the reference, the state s orgnay defned as x = (ω, δ, η), where the frequency s defned by ω = dα dt, the reference ange s defned by δ = (δ, δ 3,, δ ng ) (α α, α 3 α,, α ng α ), and α s the phase ange of node Here, f we change the reference node to be node r, the state vector needs to be redefned as x (r) = (ω, δ (r), η), where δ (r) = (δ (r), δ(r),, δ(r) r, δ(r) r+,, δ(r) n g+n) (α α r, α α r,, α r α r, α r+ α r,, α ng+n α r ) (S6) We note that, n the redefned phase space, the ange dfference between two nodes remans the same as n the orgna frame, e, δ (r) j δ (r) δ (r) j = (α α r ) (α j α r ) = α α j = δ δ j δ j Therefore, the energyke functon n the new reference frame, whch we denote Ψ(x (r) ), can be wrtten usng the expresson n Eq (4) but wth the term

10 5 = P δ repaced by n g+n = n g+n = r n g+n P δ (r) = = r P α r [ ng+n = r P δ (r) By makng use of n g+n = P =, we can further derve P ]α r = n g+n = r P α +P r α r = n g+n = P α +P α = n g+n = P α [ ng+n = P ]α = (S7) We can concude, therefore, that the energy functon remans unchanged when we choose a dfferent reference generator node In the same way we proved Eq (3), we can aso prove that where J (r) s a fu-rank matrx of the form J (r) = dx (r) dt n g+n = P δ = J (r) Ψ(x (r) ), (S8) J J (r) J (r) 3 [J (r) ]T [J (r) 3 ]T J 33 J 44 In ths matrx, the off-dagona bocks are M M J (r) = M r M r M r M r M r M r, J (r) M r M r+ M ng 3 = (S9) M r M r, (S) where J (r) 3 s an n g n matrx wth nonzero eements n row r, and the dagona bocks are the same as defned n Eq (7) of the man text To proceed, we consder the stuaton n whch a network G s spt nto k dsconnected custers G,, G k due to ne faures Wthout oss of generaty, we assume that G s the argest custer n the network and that G contans the reference node r (whch can be otherwse reassgned snce the choce of the reference node does not mpact the dynamcs or the vaue of Ψ) After rebaancng the power nput and output n G by settng G P =, we can consder the dynamcs n the correspondng subspace of the phase space of the system usng our formasm: dx (r) dt = J (r) Ψ(x (r) ) (S) Here, the substate x (r) contans ony the varabes for the nodes (generators and nongenerators) and transmsson nes n the custer G The matrx J (r) and functon Ψ(x (r) ) are defned on the custer G n the same way as J (r) and Ψ(x (r) ) were defned for the entre network G The same procedure can be used n each custer of the network, thereby eadng to a sef-consstent approach that can be apped to the genera case n whch the network spts nto dsconnected custers durng a cascade There s an exceptona case under whch our formasm w fa Recang that P r s the negatve of the mechanca power nput from the generator node r, the sum G r P = P r must be a nonnegatve number smaer than the recproca of the transent reactance of the generator The correspondng condton must hod true for each custer If for any custer we are not abe to seect a reference generator such that ths condton s satsfed, we decare t an unsovabe state Ths state corresponds to the stuaton n whch we are not abe to rebaance the power nput and output n a custer by adjustng the nput from any snge generator n the custer, and hence the cascade w necessary contnue to propagate n that custer Other operatons, such as sheddng power and adjustng the nput from mutpe generators, woud be generay needed n ths case to rebaance the system

11 6 POWER NETWORK DATA To the best of our knowedge, ths s the frst study on power grds to account for both ne overoad and generator dynamcs n the same cascade event Our anayss requres statc power fow data and dynamc data on the parameters of the generators The power fow data ncude the graph topoogy and eectrca parameters of the transmsson nes as we as the power-demand and generator-output data The power-grd network of Iceand consdered n ths work s the argest pubcy avaabe system known to us wth consstent statc and dynamc data, n whch the statc generator data can be matched to the generator data used n stabty tests Ths network has the addtona advantage of beng soated, and therefore nvoves no assumptons on possbe externa connectons The other networks used to compement our anayss are standard test systems often empoyed n power fow studes They consst of four IEEE test transmsson systems and the PEGASE 89-bus system, whch represents a transmsson network n Europe avaabe through Matpower 5 Iceand s network s shown n Fg S3 The statc and dynamc parameters of ths system are provded n Ref [4] The avaabe data ncude the rea power P (m) supped by each generator, the rea power P (d) demanded at each node, the reactance x of each transmsson ne, and the dynamc data (ncudng rotor nerta M and transent reactance) of each generator We choose the oad frequency rato to be T = for every nongenerator node and the rotor dampng rato to be D = 5 for every generator n the network The assumpton that these parameters are the same for the dfferent nodes s not essenta; other choces are possbe, and are consstent wth the mathematca assumptons underyng Eq (S6) and our formasm, provded that the parameters n the denomnators of Eqs (6) and (7) are nonzero If there are two nes connectng a par of nodes, we repace them by a snge ne wth the combned mpedance For one of the nes, whose actua reactance s not avaabe n the data, the reactance s assgned to be pu to assure that ths ne has a arge capacty To defne the capacty of each transmsson ne, we cacuate the steady-state power-fow soutons [determned by Eq (S)] for a possbe one-ne faure scenaros that keep the network connected We then assgn the capacty W of each ne to be % of the maxmum reactance energy that the ne stores n these soutons In ths way, the system has at east one stabe steady state when a transmsson ne s dsconnected (as ong as the network remans connected, as assumed) The basc propertes of the test systems are sted n Tabe S [36], where we aso ncude Iceand s network for competeness In partcuar, the IEEE 4-bus test system n Fg 3 conssts of 5 generator nodes, 4 nongenerator nodes (e, buses), and transmsson nes The nformaton avaabe on these systems ncudes data on the rea power P (m) supped by each generator, the rea power P (d) requested at each node, and the reactance x for each transmsson ne [36] In our cacuatons, we keep a these parameters unchanged except for the P (m) of the generators wth zero output of rea power To ncude those generators n our formasm, we assume that they have a sma rea-power output of MW We aso assume that a generator rotors have dentca dynamc parameters, wth rotor nerta M = 5 (n seconds) and rotor dampng rato D = 5 The transent reactance of each generator s chosen to be pu, whch guarantees that the ange dfference between the two ends of the ne connectng the generator node (and hence nstabty) remans sma durng cascades Fnay, the oad frequency rato T and ne capactes W are assgned n the same way as n Iceand s power grd

12 7 Tabe S Descrpton of the power systems used n ths work The coumns represent the number of buses (n), generators (n g), and power nes (n ), respectvey Power systems n n g n IEEE4 4 5 IEEE IEEE PEGASE IEEE Iceand Tabe S Cascades trggered by snge- and doube-ne remova perturbatons on the test power systems consdered Here, N s the number of perturbatons that keep the network ntay connected and N T s the number of such perturbatons that trgger cascades Aso shown s the number of doube-ne removas N OM for whch the resutng C depends on the perturbaton schedue ()-() (defned n the man text), as we as the breakdown nto the number for whch the argest C resuts from ether () or () (N E ), from both () and () (N B ), and from () (N C ) Smar resuts for Iceand s power grd are presented n the man text System Snge Doube N N T N N T N OM N E N B N C IEEE IEEE IEEE PEGASE IEEE f(`)! 6` ` = 6` = :5 6` = 6` = :5 5 ` Fgure S Equvaent of Fg (a) [e, the rhs of Eq ()] for dfferent vaues of λ

13 8 (a) (b) Fgure S Network representatons for a hypothetca 5-bus power system (a) Orgna representaton, where each node s assocated wth a oad (bue arrow) and/or a redstrbuton pont, whch may be connected to one or more generators (red arrows) The wde gray arrows ndcate the drecton of the fow on the nes (b) Extended representaton, where each generator s modeed as an addtona node (square) that njects power nto the orgna network through a vrtua ne (wde red arrow) Fgure S3 Network representaton of Iceand s power-grd system Generators are ndcated by squares and nongenerator nodes by crces Under the condtons consdered here, there are snge-ne remova perturbatons (red) that ead to cascades n ths system

14 9 ` (a) (b) `9- `36-56 ` ! C (c) ` `9- ` Fgure S4 Impact of perturbaton order n Iceand s power grd for the remova of nes `36-56 and `9- accordng to the three scenaros consdered n the man text: (a) remova of `36-56 foowed by the remova of `9- ; (b) remova of `9- foowed by the remova of `36-56 ; (c) concurrent remova of `36-56 and `9- The dynamcs of the power grd s represented by the status of the transmsson nes η` (top panes), the frequency of the generators ω (mdde panes), and the sze of the argest custer C (bottom pane) ` (a) (b) `56-6 `9- (c) `9- ` `56-6 ` ! C Fgure S5 Same as Fg S4 for the remova of nes `56-6 and `9- In ths case the concurrent remova of the two nes has the hghest mpact, whe n the exampe of Fg S4 a cascade can ony be trggered by the separate remova of the nes

15 !!! (a) ` /! / *(x) (b) 5 5 # 5 ` !! `36-6 `36-5 ` *(x) /! # 5!!!!s !!!!s 3 - Fgure S6 Rescue perturbaton n Iceand s power-grd system (a) The remova of ne 36-6 eventuay eads to the separaton of generator due to the overoad of ne 4-45, whch causes the network to ose 89 MW n power generaton (38% of ts tota) (b) The concurrent remova of ne 36-5 aong wth 36-6 prevents subsequent faures and power osses, keepng the network connected The top panes show the status of the transmsson nes η and the bottom panes show the frequency of the generators ω