Control of Power Converters in Low-Inertia Power Systems
Florian D¨
- rfler
Control of Power Converters in Low-Inertia Power Systems Florian D - - PowerPoint PPT Presentation
Control of Power Converters in Low-Inertia Power Systems Florian D orfler Automatic Control Laboratory, ETH Z urich Acknowledgements ! ! ! ! Marcello Colombino Jean Sebastien Brouillon Dominic Gro Irina Subotic Further: Gab-Su
!
! ! !
1
Hz *10 sec BEWAG UCTE
2
Hz *10 sec BEWAG UCTE
December 7, 1994
3
2Mω2 as safeguard against disturbances
4
2Mω2 as buffer
5
issues broadly recognized by TSOs, device manufacturers, academia, agencies, etc.
UPDATE REPORT ! BLACK SYSTEM EVENT IN SOUTH AUSTRALIA ON 28 SEPTEMBER 2016
AN UPDATE TO THE PRELIMINARY OPERATING INCIDENT REPORT FOR THE NATIONAL ELECTRICITY MARKET. DATA ANALYSIS AS AT 5.00 PM TUESDAY 11 OCTOBER 2016.
ERCOT is recommending the transition to the following five AS products plus one additional AS that would be used during some transition period:
ERCOT CONCEPT PAPER Future Ancillary Services in ERCOT
PUBLIC The relevance of inertia in power systems Pieter Tielens n, Dirk Van Hertem
ELECTA, Department of Electrical Engineering (ESAT), University of Leuven (KU Leuven), Leuven, Belgium and EnergyVille, Genk, Belgium Contents lists available at ScienceDirect journal homepage: www.elsevier.com/locate/rserRenewable and Sustainable Energy Reviews
Renewable and Sustainable Energy Reviews 55 (2016) 999–1009MIGRATE project: Massive InteGRATion of power Electronic devices
Frequency Stability Evaluation Criteria for the Synchronous Zone
– Requirements and impacting factors – RG-CE System Protection & Dynamics Sub Group
However, as these sources are fully controllable, a regulation can be added to the inverter to provide “synthetic inertia”. This can also be seen as a short term frequency support. On the other hand, these sources might be quite restricted with respect to the available capacity and possible activation time. The inverters have a very low
Impact of Low Rotational Inertia on Power System Stability and Operation
Andreas Ulbig, Theodor S. Borsche, Göran Andersson ETH Zurich, Power Systems Laboratory Physikstrasse 3, 8092 Zurich, Switzerland ulbig | borsche | andersson @ eeh.ee.ethz.ch
!
!"#$%% "&'()*%")+,-.)'%/),-)0% 1"2%/).3**)456(-34'%
!"#$%&!&$!&'"!()*!+$,,-&&""!
6
(Invited Paper)
Federico Milano
University College Dublin, Ireland email: federico.milano@ucd.ie
Florian D¨
ETH Z¨ urich, Switzerland emails: dorfler@ethz.ch, ghug@ethz.ch
David J. Hill∗ and Gregor Verbiˇ c
University of Sydney, Australia
∗ also University of Hong Kong
emails: dhill@eee.hku.hk, gregor.verbic@sydney.edu.au
include key features without burdening the model (whether for analytical or computational work) with uneven and excessive detail;
devices and time-scales associated with CIG, new loads and use of storage;
guidelines including data-based approaches;
mitigate the high rate of change of frequency in low inertia systems;
very fast control system, which are nearly antipodal characteristics to those of a synchronous machine. Thus, one should critically reflect the control of a converter as a virtual synchronous machine; and
(and cannot) be fixed by simply “adding inertia back” in the systems. The later sections contain many suggestions for further work, which can be summarized as follows:
[Milano, D¨
7
8
8
xb(t) xc(t)
xb(t + T) xc(t + T)
1 T
0 xi(t)dt = 0
π 2π −1 1
δ xabc
sin(δ(t) − 2π
3 )
sin(δ(t) + 2π
3 )
π 2π −1 1
δ xabc
sin(δ0 + ω0t − 2π
3 )
sin(δ0 + ω0t + 2π
3 )
π 2π −1 1
δ xabc
2 ) 90◦rotation
9
◮ quasi-steady state algebraic model ∼ diffusive (synchronizing) coupling
i1 . . . in
nodal injections
. . . ... . . . ... . . . −yk1I2 · · · n
j=1 ykjI2
· · · −yknI2 . . . ... . . . ... . . .
v1 . . . vn
nodal potentials ◮ salient feature: local measurement reveal global information
10
network idc
DC port modulation LC output filter AC port control to power grid (3-phase)
◮ passive DC port port (idc, vdc) for energy balance control
◮ modulation ≡ lossless signal transformer (averaged)
2vdcu with u ∈ [−1, 1] ◮ LC filter to smoothen harmonics with R, G modeling filter/switching losses
11
k io,k = p⋆ k ,
k R( π 2 ) 90◦rotation
k
k, q⋆ k}
12
v1 v2 θ⋆
12
v⋆ v1 v2
ω0 ω0
k, p⋆ k, q⋆ k) & stabilization of a limit cycle
13
!""" #$%&'%(#!)&' )& *)+"$ ','#"-'. /)01 23. &)1 2. -%, 2456 5676
!89:;8;<=><? />@=AB: !<;@=>B >< CD!EFGBH;I +><I *JK;@ E;<;@B=>J<
!"#$%&'"'() %* +$,(-.'() /'-#%(-' .( 0.1$%2$.3- 4-.(2 5.$)6,7 !('$).,
8.".-9 :%(.! "#$%&'# (&)*&+! ,---; :6$<,(,$,<,(, =%%77,! (&)*&+! ,---; ,(3 06>67 ?@ ?9,(3%$>,$! (&)*&+! ,---
Virtual synchronous generators: A survey and new perspectives
Hassan Bevrani a,b,⇑, Toshifumi Ise b, Yushi Miura b
a Dept. of Electrical and Computer Eng., University of Kurdistan, PO Box 416, Sanandaj, Iran b Dept. of Electrical, Electronic and Information Eng., Osaka University, Osaka, Japan!"#$%&'()*+,-+#'"(./#0*/1(2-33/*04($(5&*0-$1(( 6#+*0&$(7*/8&9+9(:"(!&;0*&:-0+9(<#+*="(20/*$=+( 0/(6;/1$0+9(7/>+*(2";0+%;((
?$-0@&+*(!+1&11+A(!"#$"%&'()))A(B*-#/()*$#C/&;A(*"+,-%'!"#$"%&'()))A($#9(?&11+;(D$1$*$#=+(
!""" #$%&'%(#!)&' )& *)+"$ ','#"-'. /)01 23. &)1 2. -%, 2456
!789:;< "=>?<:;@7 (@7:9@? ':9<:8AB C@9 /'(DE/F( #9<7G=;GG;@7 'BG:8=G
H;8I8; JK>. (<=LI8?? F1 M@@:K. N9<;7 *1 %O<=. %7O98P H1 $@GQ@8. <7O (K9;G N1 M9;AK:
!"#$%&#'$%()*+'",'"%-#,.%/#",012%3#*',#4% 5,)"16'%
7898:%+1*%;'<'*=''4>?%@%58;8A8%$'%A11*?@%!"#$%&'("()"&*'+,,,@%98%/1"'21B%1*$%38%/#<<4.'"C@% !"#$%&'("()"&'+,,,%
◮ PD control on ω(t) : M d dt ω(t) + D (ω(t) − ω0) = Pgeneration(t) − Pdemand(t) ◮ there are smarter implementations at the cost of algorithmic complexity
14
tracking control (cascaded PIs)
3 2 4 4 1
reference synthesis
15
[Chandorkar, Divan, Adapa, ’93] ◮ frequency control by mimicking p − ω
◮ voltage control via q − v droop heuristic: d dtv = −c1(v − v⋆) − c2(q − q⋆)
P2 P1 P ! !* !sync
ωsync ω p(t) − p∗ ω0
filtering
logic for sync droop tracking controllers tricks
16
Virtual synchronous generators: A survey and new perspectives
Hassan Bevrani a,b,⇑, Toshifumi Ise b, Yushi Miura b
a Dept. of Electrical and Computer Eng., University of Kurdistan, PO Box 416, Sanandaj, Iran b Dept. of Electrical, Electronic and Information Eng., Osaka University, Osaka, JapanContents lists available at ScienceDirect
Electrical Power and Energy Systems
journal homepage: www.elsevier.com/locate/ijepes Abstract- The method to investigate the interaction between a Virtual Synchronous Generator (VSG) and a power system is presented here. A VSG is a power-electronics based device that
To better study and witness the effects of virtual inertia, the hardware of a real VSG should be tested within a power
a power system is not easy as a power system cannot be
Real Time Simulation of a Power System with VSG Hardware in the Loop
Vasileios Karapanos, Sjoerd de Haan, Member, IEEE, Kasper Zwetsloot Faculty of Electrical Engineering, Mathematics and Computer Science Delft University of Technology Delft, the Netherlands E-mails: vkarapanos@gmail.com, v.karapanos@tudelft.nl, s.w.h.dehaan@tudelft.nl
Frequency Stability Evaluation Criteria for the Synchronous Zone
– Requirements and impacting factors – RG-CE System Protection & Dynamics Sub Group
However, as these sources are fully controllable, a regulation can be added to the inverter to provide “synthetic inertia”. This can also be seen as a short term frequency support. On the other hand, these sources might be quite restricted with respect to the available capacity and possible activation time. The inverters have a very low
[Jouini, Arghir, & D¨
17
droop control machine emulation PBC + matching
slow convergence
[Tayebi, D¨
18
0v + g(v) = io
’13], [Dhople, Johnson, D¨
19
P2 P1 P ! !* !sync ωsync ω p(t) − p∗ ω0
v!" idc m iI v LI τm θ, ω vf v if τe is Lθ C!" M rf rs rs G!" RI
R C L g(v) v +
k
L
˙
xαβ
y a constant: R C
L
n n
∗
[Colombino, Groß, Brouillon, & D¨
[Seo, Subotic, Johnson, Colombino, Groß, & D¨
20
d dt vk(t) = uk(vk, io,k)
io,k to network
dti + ri ≈
jk}
k, p⋆ k, q⋆ k)
θ⋆
12
v⋆ v1 v2
ω0 ω0 21
synchronization
v1 v2
ω0 ω0 eθ,1(v) eθ,2(v) e v,1(v1) e v,2(v2)
j=1 wjk
jk)vk
22
jk)vk)
jk
jk))vk
k = v⋆2 j rjk(1−cos(θ⋆
jk))−ω0ℓjk sin(θ⋆ jk)
r2
jk+ω2 0ℓ2 jk
k = −v⋆2 j ω0ℓjk(1−cos(θ⋆
jk))+rjk sin(θ⋆ jk)
r2
jk+ω2 0ℓ2 jk
global parameters
k
k
k
k
23
−ω0 ω0
j=1 wjk(vj − R(θ⋆ jk)vk)
−ω0 ω0
v⋆2
k
p⋆
k
−p⋆
k q⋆ k
24
◮ consistent v⋆, p⋆ k, and q⋆ k satisfy AC power flow equations ◮ magnitude control slower than synchronization control ◮ power transfer “small enough” compared to network connectivity
1 2 λ2(L) > max k
j=1
k
vk≈1 ω0 + c1 (p⋆ k − pk)
vk≈1 c1 (q⋆ k − qk) + c2 (v⋆ − vk)
25
2dist(v, S)2 + c2 v⋆2
Z{02N } 0-stable manifold
sync set S amplitude set A
T
target set T
02N
t→∞v(t)T ∪02N = 0
t→∞v(t)T = 0
26
1 2 3 v1 v2 v3 4 8 6 5 9 7
27
5 10 15 0.5 1 1.5 2 pk [p.u.] 5 10 15 0.99 1 1.01 time [s] ω [p.u.] 5 10 15 0.5 1 vk [p.u.] 5 10 15 0.5 1 1.5 2 time [s] io,k [p.u.]
28
◮ inverter control interferes
◮ controller needs to be
◮ recognized problem [Vorobev, Huang, Hosaini, & Turitsyn,’17]
◮ communication through
◮ very inductive lines delay
◮ the controller must be slow
29
d d t v = fv(v, i)
ia = h(v) −h(v)
d d t i = fi(v, i)
v i v y = i − h(v)
◮ slow system: V (v) for d d tv = fv(v, h(v)) ◮ fast system: W(y) for d d ty = fi(v, y + h(v))
d d tv = 0 & coordinate y = i − h(v)
◮ ν(x) = dW(i − h(v)) + (1 − d)V (v)
◮ d d tν(x) is decaying under stability condition
◮ T ′ ∪ {0n} globally attractive & T ′ stable ◮ Z{0n} has measure zero
30
5 10 15 20 10−5 10−4 10−3
3 · 1
−2
3 · 10
− 2
6 · 10−2 6 · 10−2 8 · 10−2 8 · 10−2 9.5 · 10−2 9.5 · 10−2
linear instability certified stability region constraints violated damping ratios
2 4 1 2 vk [p.u.]
31
d dt v(t) = u(v, io)
io
vdc i L R C v G io 1 2vdc u
◮ idea: invert LC filter so that v ≈ 1 2vdcu
◮ analysis: repeat proof via singular perturbation Lyapunov functions
[Subotic, ETH Z¨ urich Master thesis ’18] ◮ ...similar steps for control of vdc in a more detailed model
32
33
[Seo, Subotic, Johnson, Colombino, Groß, & D¨
black start of inverter #1 under 500 W load (making use of almost global stability) 250 W to 750 W load transient with two inverters active connecting inverter #2 while inverter #1 is regulating the grid under 500 W load change of setpoint: p⋆ of inverter #2 updated from 250 W to 500 W
34
dynamics on grid-forming dispatchable virtual oscillator control.
ization of coupled oscillators with application to the control of grid-forming power inverters.
35