Power flow dynamics in electricity grids can be described by equations resembling a Kuramoto model of non-linearly coupled oscillators with inertia. The coupling of the oscillators or nodes in a power grid generally exhibits pronounced heterogeneities due to varying features of transmission lines, generators, and loads. In studies aiming at uncovering mechanisms related to failures or malfunction of power systems, these grid heterogeneities are often neglected. However, over-simplification can lead to different results away from reality. We investigate the influence of heterogeneities in power grids on stable grid functioning and show their impact on estimating grid stability. Our conclusions are drawn by comparing the stability of an Institute of Electrical and Electronics Engineers test grid with a homogenized version of this grid.

1.
S.
Backhaus
and
M.
Chertkov
,
Phys. Today
66
,
42
(
2013
).
2.
M.
Timme
,
L.
Kocarev
, and
D.
Witthaut
,
New J. Phys.
17
,
110201
(
2015
).
3.
S.
Auer
,
K.
Kleis
,
P.
Schultz
,
J.
Kurths
, and
F.
Hellmann
,
Eur. Phys. J. Spec. Top.
225
,
609
(
2016
).
4.
M.
Mureddu
,
G.
Caldarelli
,
A.
Damiano
, and
H.
Meyer-Ortmanns
,
Sci. Rep.
6
,
34797
(
2016
).
5.
C.
Schiel
,
P. G.
Lind
, and
P.
Maass
,
Sci. Rep.
7
,
11562
(
2017
).
6.
H.
Hähne
,
J.
Schottler
,
M.
Wächter
,
J.
Peinke
, and
O.
Kamps
,
Europhys. Lett.
121
,
30001
(
2018
).
7.
S.
Tamrakar
,
M.
Conrath
, and
S.
Kettemann
,
Sci. Rep.
8
,
6459
(
2018
).
8.
B.
Schäfer
,
M.
Matthiae
,
X.
Zhang
,
M.
Rohden
,
M.
Timme
, and
D.
Witthaut
,
Phys. Rev. E
95
,
060203(R)
(
2017
).
9.
B.
Schäfer
,
C.
Beck
,
K.
Aihara
,
D.
Witthaut
, and
M.
Timme
,
Nat. Energy
3
,
119
(
2018
).
10.
M.
Colombino
,
D.
Groß
, and
F.
Dörfler
, in 2017 IEEE 56th Annual Conference on Decision and Control (CDC) (IEEE, 2017), pp. 5690–5695.
11.
T.
Nishikawa
and
A.
Motter
,
New J. Phys.
17
,
015012
(
2015
).
12.
We note that load nodes are more realistically described by a constant impedance to the ground,24 which guarantees that they consume a given power at the equilibrium voltage.
13.
K.
Schmietendorf
,
J.
Peinke
,
R.
Friedrich
, and
O.
Kamps
,
Eur. Phys. J. Spec. Top.
223
,
2577
(
2014
).
14.
K.
Schmietendorf
,
J.
Peinke
, and
O.
Kamps
,
Eur. Phys. J. B
90
,
222
(
2017
).
15.
G.
Filatrella
,
A. H.
Nielsen
, and
N. F.
Pedersen
,
Eur. Phys. J. B
61
,
485
(
2008
).
16.
S.
Lozano
,
L.
Buzna
, and
A.
Díaz-Guilera
,
Eur. Phys. J. B
85
,
231
(
2012
).
17.
P.
Menck
,
J.
Heitzig
,
J.
Kurths
, and
H.
Schellnhuber
,
Nat. Commun.
5
,
3969
(
2014
).
18.
R.
Albert
,
I.
Albert
, and
G. L.
Nakarado
,
Phys. Rev. E
69
,
025103
(
2004
).
19.
J.-W.
Wang
and
L.-L.
Rong
,
Saf. Sci.
47
,
1332
(
2009
).
20.
P.
Hines
,
E.
Cotilla-Sanchez
, and
S.
Blumsack
,
Chaos
20
,
033122
(
2010
).
21.
C.
Grigg
,
P.
Wong
,
P.
Albrecht
,
R.
Allan
,
M.
Bhavaraju
,
R.
Billinton
,
Q.
Chen
,
C.
Fong
,
S.
Haddad
,
S.
Kuruganty
et al.,
IEEE Trans. Power Syst.
14
,
1010
(
1999
).
22.
Even in the 2006 European blackout, the largest frequency deviations were of the order of 1 Hz.41 
23.
P. J.
Menck
,
J.
Heitzig
,
N.
Marwan
, and
J.
Kurths
,
Nat. Phys.
9
,
89
(
2013
).
24.
C.
Mitra
,
T.
Kittel
,
A.
Choudhary
,
J.
Kurths
, and
R. V.
Donner
,
New J. Phys.
19
,
103004
(
2017
).
25.
J.
Machowski
,
J.
Bialek
, and
J.
Bumby
,
Power System Dynamics
(
John Wiley & Sons
,
NJ
,
2008
).
26.
R.
Zimmerman
and
C.
Murillo-Sánchez
, Matpower 6.0b1 User Manual (Power Systems Engineering Research Center (PSerc), 2016) (Online), see www.pserc.cornell.edu/matpower/manual.pdf.
27.
The steady-state values of these admittances can be affected in transient time intervals by the armature reaction effect,25 which will not be considered here.
28.
F.
Dörfler
and
F.
Bullo
,
IEEE Trans. Circuits Syst. I: Regul. Pap.
60
,
150
(
2013
).
29.
The shift of π/2 in this definition of the phase angle γjk is introduced to yield a sine function in the subsequent equation (7), which then has the common form of a Kuramoto model of second order.42,43,44
30.
M.
Rohden
,
A.
Sorge
,
M.
Timme
, and
D.
Witthaut
,
Phys. Rev. Lett.
109
,
064101
(
2012
).
31.
“UCTE” operation handbook (2004), see https://www.entsoe.eu.
32.
D.
Mehta
,
N.
Daleo
,
F.
Dörfler
, and
J.
Hauenstein
,
Chaos
25
,
053103
(
2015
).
33.
R.
Delabays
,
T.
Coletta
, and
P.
Jacquod
,
J. Math. Phys.
57
,
032701
(
2016
).
34.
D.
Manik
,
M.
Timme
, and
D.
Witthaut
,
Chaos
27
,
083123
(
2017
).
35.
An alternative method for determining the fixed point is the “holomorphic embedding load flow method” suggested in Ref. 45.
36.
A.
Motter
,
S.
Myers
,
M.
Anghel
, and
T.
Nishikawa
,
Nat. Phys.
9
,
191
(
2013
).
37.
D.
Witthaut
and
M.
Timme
,
New J. Phys.
14
,
083036
(
2012
).
38.
D.
Jung
and
S.
Kettemann
,
Phys. Rev. E
94
,
012307
(
2016
).
39.
Let us note that a corresponding balance under consideration of losses, i.e., when keeping the resistances of the lines, can usually be ensured only if one generator node is allowed to exhibit a mechanical power different from the other generator nodes. Accordingly, in that case, one would not consider a fully homogenized grid.
40.

If one would be interested only in the grid return times (see below), an alternative way to obtain smoothly varying times is to consider the energy E=j12Hjωj2+jk14Kjkcos(θjθkγjk)[(θjθj)(θkθk)]2 of coupled harmonic oscillators, which characterizes the system in the neighborhood of the fixed point, see Eq. (7). If the phase space trajectory {[θj(t),ωj(t)]} enters the nearest neighborhood, this energy decreases monotonically with time t.

41.
C.
Li
,
Y.
Sun
, and
X.
Chen
, in 2007 International Power Engineering Conference—IPEC 2007 (IEEE, NJ, 2007), pp. 939–944.
42.
Y.
Kuramoto
, in Mathematical Problems in Theoretical Physics, Lecture Notes in Physics, edited by H. Araki (Springer Verlag, Berlin, 1975), Vol. 39, pp. 420–422.
43.
J. A.
Acebrón
,
L. L.
Bonilla
,
C. J.
Pérez Vicente
,
F.
Ritort
, and
R.
Spigler
,
Rev. Mod. Phys.
77
,
137
(
2005
).
44.
F. A.
Rodrigues
,
T. K. D.
Peron
,
P.
Ji
, and
J.
Kurths
,
Phys. Rep.
610
,
1
(
2016
).
45.
A.
Trias
, in 2012 IEEE Power and Energy Society General Meeting (IEEE, 2012), pp. 1–8.
You do not currently have access to this content.