Hybrid Time Domain Simulation: Application To Fault Induced Delayed Voltage Recovery (FIDVR)
Hybrid Time Domain Simulation: Application To Fault Induced Delayed Voltage Recovery (FIDVR)
Hybrid Time Domain Simulation: Application To Fault Induced Delayed Voltage Recovery (FIDVR)
Vijay Vittal
Arizona State University
PSERC Webinar
Jan. 20, 2015
Presentation background
2
Outline
1. Introduction
2. Part -1: OpenHybridSim: A new hybrid
simulation platform
3. Part -2: Application to FIDVR study
4. Conclusions
3
Introduction: Background
System
Zik
6
Time scale of transient phenomena
7
OpenHybridSim:
A new hybrid simulation tool
Overall design
• Loosely decoupled architecture
• Socket communication
• InterPSS, an open source power system simulator
• A generic interface to an EMT simulator, e.g.,
PSCAD, Maltab/SimPowerSystems
10
Three-phase Thévenin equivalents
of the external system
• Step 1: Calculate 3-sequence
Norton equivalents
zia +_ VTia
ziab
• Step 2: 3-sequence to 3-phase bus i ziac
transformation for each boundary bus zibc
zib +
_ VTib
120 −1
I Niabc = SI 120
Ni
yi
abc
= Syi S
zic + VTic
……
_
Detailed
……
where S is the transformation matrix, zikabc
is the bus primitive admittance system
zka
zkab
+
_ VTka
• Step 3: 1-phase Norton to zkac
Thévenin for each boundary bus bus k zkbc
zkb +_
VTkb
VTip = I Nip / yip zip = 1 / yip zkc +_
VTkc
where p stands for one of the three phases
Single-phase line Three-phase line
11
Detailed system represented by sequence
current sources in TS simulation
Need to extract three-sequence, fundamental frequency
currents from instantaneous, waveform values
Use well-established FFT algorithm
• PSCAD provides an FFT component
Directly used in the network solution (I=YV) in TS
simulation
Positive-sequence network
(1) solution and integration step (1)
I V TS ( t +∆T )
EMT ( t )
x (t ) = f ( x(t ), y ) (1)
(t )
x(t + ∆T ) = x(t ) + x (t ) ∆T
=0 g1 ( x(t + ∆T ), y((1)t + ∆T ) , I EMT
(1)
(t ) )
13
The interaction protocols
between two simulators
2
TS
(a) series 3
5
…
1
4
.. .. ..
…
. . . EMT
T T+ΔT T+2ΔT T+3ΔT
TS
…
(b) parallel
…
... ... ...
EMT
(c) combined T T+ΔT T+2ΔT T+3ΔT
TS
… …
… …
... ... ... ... ... ...
EMT
T T+ΔT T+2ΔT T+3ΔT T+4ΔT T+5ΔT T+6ΔT
I (1)
Positive-sequence network
solution and integration step (1)
VTS ( t +∆T )
InterPSS
x(t ) EMT ( t )
x (t ) = f ( x(t ), y((1)t ) )
x(t + ∆T ) = x(t ) + x (t ) ∆T x(t + ∆T )
V 120 Calculate 3-phase
TS (t + ∆T )
socket =0 g1 ( x(t + ∆T ), y((1)t + ∆T ) , I EMT
(1)
(t ) )
socket
server 120
Thévenin server
Negative- and zero- I EMT (t ) equivalent V abc
I (2)
VTS(2)(t +∆T ) T (t + ∆T ) VTS120(t +∆T )
VTS120(t ) EMT ( t ) sequence network solver
0 = g2 ( y (2)
( t + ∆T ) ,I (2)
EMT ( t ) )
(0)
I EMT VTS(0)(t +∆T ) Step (3)
120
(t )
0 = g0 ( y (0)
( t + ∆T ) ,I (0)
)
I EMT ( t )
EMT ( t )
V abc
Step (1) Step (2) Three-sequence TS simulation T (t + ∆T )
Step (4)
120 120
I EMT ( t ) I EMT ( t +∆T )
EMT EMT EMT Step (5) EMT EMT
socket
client step step step …………………………………….. step step
1 2 3 EMT N-1 N
t t + ∆T
15
Implementation with the parallel type protocol
t t + ∆T
16
Automatic protocol switching algorithm
Automatically switch interaction protocol based on the
system conditions, reflected by the maximum rate of
change of sequence current injections
( p) ( p)
I −I
120 EMT ( i , t ) EMT ( i , t − ∆T )
RI max( max ( )) / ∆T
EMT ( t ) (1)
i p ∈ (1, 2, 0) I
EMT ( i , t − ∆T )
120
1
I EMT ( t ) Maximum Y series
rate of
change
120
RI EMT >ε? N
I 120
EMT (t −∆T ) 120 Y 0
RI EMT N last step Delay
used series? parallel
time step
∆T
Logic of the protocol switching algorithm
17
Automatic creation and initialization
of the external network
Boundary
bus i
Base case Ii120 bus m
……
System system
to identify bus n
Boundary detailed Vk120
Configuration Ik120
system bus k
Boundary
bus i
bus m
1) Set buses and
Ii120
The remainder
branches within the
of the external detailed system out of Power flow
……
system
bus k service
bus n
2) initialize the external
Ik120 Vk120
system
18
Testing of the developed platform
Test system
Bus 8
230 kV
Bus2 Bus 7
18 kV 230 kV
0.0085+j0.072
100+j35MVA Bus 9 Bus 3
230 kV 13.8 kV (case 1)
0.0119+j0.1008
G2 G3
0.032+j0.161
0.039+j0.17
192MVA T2 T3 128MVA
XT=0.0625 XT=0.0586
Bus 5 Bus 6
230 kV 230 kV
0.017+j0.092
0.01+j0.085
125+j50MVA
90+j30MVA
Bus 4
230 kV
T1
XT=0.0576
Bus 3
16.5 kV
G1
247.5MVA
The IEEE 9 Bus system Modeling of Bus 5 and the external network
Thévenin equivalents in PSCAD
19
Case 1: The interaction protocols
Case 1 setting: 12
PSCAD
- EMT time step = 50 us 10 EMT-TS(series)
EMT-TS(parallel)
-TS time step = 5 ms
EMT-TS(combined)
8
current (kA)
Protocol control setting:
6
-Threshold ε = 0.004
- Delay = 2 cycles 4
230 kV/69kV
G2
0.0085+j0.072 0.0119+j0.1008
G3
ground fault 150 MVA
10%
• starting at 2.0 s
0.032+j0.161
69kV
0.039+j0.17
192MVA T2 T3 128MVA
XT=0.0625 XT=0.0586
• 4 cycles 69/12.47 kV
Bus 5 Bus 6 50 MVA 69/12.47 kV
230 kV 50 MVA
230 kV 8%
0.017+j0.092
0.01+j0.085
8%
125+j50MVA
90+j30MVA 12.47 kV 12.47 kV 12.47 kV
Bus 4 Equivalent
230 kV feeder
T1
XT=0.0576
Bus 3
Composite Composite Composite
16.5 kV
Load Load Load
G1
Model Model Model
247.5MVA
Detailed modeling of Bus 5
[1] Y. Liu, V. Vittal, J. Undrill, and J. H. Eto, "Transient Model of Air-Conditioner Compressor Single Phase Induction Motor," IEEE
Transactions on Power Systems, vol. 28, pp. 4528-4536, 2013.
21
Case 2: Performance of the developed
platform under unsymmetrical faults
PSCAD EMT-TS(parallel) EMT-TS(combined) protocol signal
200 2
PSCAD
100
Voltage (kV)
EMT-TS(parallel)
1.5
0 EMT-TS(combined)
protocol signal
-100
1
-200
Phase A
-300
1.95 2 2.05 2.1 2.15 2.2 0.5
200
Current (kA)
0
100
Voltage (kV)
0
-0.5
-100 0.7
-200 -1
Phase B 0.6
-300
1.95 2 2.05 2.1 2.15 2.2
-1.5 0.5
200
0.4
100 -2
Voltage (kV)
22
Takeaways of part-1
23
Fault induced delayed voltage recovery (FIDVR)
Fault
• distribution, sub-
transmission and
transmission
systems
Delayed voltage
recovery
• several to tens of
seconds
Root cause of
FIDVR
Voltage profile during a typical FIDVR event [2]
• air conditioner
(A/C) compressor
motor stalling and [2] D. N. Kosterev, A. Meklin, J. Undrill, B. Lesieutre, W. Price, D. Chassin, et al., "Load modeling
in power system studies: WECC progress update," in 2008 IEEE Power and Energy Society General
Meeting - Conversion and Delivery of Electrical Energy in the 21st Century, 2008, pp. 1-8.
prolonged tripping 24
Fault induced delayed voltage recovery (FIDVR)
Accurate FIDVR studies require
• Detailed A/C modeling
• Detailed network model down to distribution feeder level
Limitations of the conventional positive-sequence TS
simulator
• Distribution system configurations
• Response of single-phase devices when subjected to
unsymmetrical faults
FIDVR events are generally localized
• Detailed modeling can be limited to a small portion of a large
power system
FIDVR A large
Area power system
25
Determination of the boundary
of detailed system
Use a simple yet generic test case to quantify the voltage dip
threshold at a transmission bus causing A/C motor stalling
Criterion: A bus is included in the internal system if a single-
phase or three-phase fault at that bus causes a phase-to-neutral
voltage at buses with a large percentage of A/C motor loads to
drop below 0.75 pu.
x = 0.025 x = 0.05 x = 0.1 x = 0.15 x = 0.2
T T T T T
0.5 0.5
Total number
238
24086
15021
of buses
sub 15090
sub
500 kV 7
24092
24236
Number of 230 kV 37
24801 EAST
24138 buses of 161 kV 3
sub
28040 different 115 kV 68
15093
sub
24151
sub
voltage levels 92 kV 18
115 kV
500 kV <= 66kV 105
M M M
Load Load sub
sub Sub-system below 500 kV Total Load 11.9 GW
One-line diagram of the study region
27
Initialization of the detailed system with a
large percentage of load as induction motors
The built-in initialization function of PSCAD fails to
initialize the detailed system
A two-step initialization approach
Bus HV
(1) The switch K is turned to the position 0, LV
s
1 0
and the distribution systems and the CLMs
K Static
are energized by the fixed voltage sources equivalent
load
and initialized independently.
(2) After the CLMs are successfully Distribution network
initialized, the switch K is turned to the
position 1 such that distribution systems are
CLM CLM CLM
connected to the transmission system.
Voltage magnitude(pu)
1
is applied on the phase A of bus EMT-TS
InterPSS
24151 at t = 0.2 s 0.95
0.9
Bus 24151 positive sequence voltage
1.1
0.85
1.05 0 0.1 0.2 0.3 0.4 0.5
Voltage magnitude(pu)
Time (s)
0.8 0.19
0.75
0.17
0.7
0 0.1 0.2 0.3 0.4 0.5 0.15
Time (s) 0 0.1 0.2 0.3 0.4 0.5
Time(s)
29
FIDVR event triggered by a SLG fault
Detailed modeling of the region served by Bus 24151
Bus 24151
500 kV
415 +j*75 MVA 415 +j*75 MVA 403 +j*96 MVA 403 +j*96 MVA
22 + j*1.1
46.8 MVAr 0.12 pu
MVA
250 MVA
Bus 25423 115 /12.47 kV
115 kV 0.12 pu
250 MVA
115 /12.47 kV
M M M
f-2 f-3 f-4 f-5 f-6 f-7 f-8
Equivalent Feeder f-1
Equivalent Equivalent Equivalent Equivalent Equivalent Equivalent Equivalent Equivalent
feeder feeder feeder feeder feeder feeder feeder feeder
Model Model Model Model Model Model Model Model
202+j*(-8) MVA 202+j*(-8) MVA 202+j*(-8) MVA 202+j*(-8) MVA 202+j*(-8) MVA 202+j*(-8) MVA 202+j*(-8) MVA 202+j*(-8) MVA
Feeder 1
Z
A 3
Z
4 feeder 4 feeder
modeled in B
C
two sections M M
A/C
30
FIDVR event triggered by a SLG fault (cont’)
Terminal voltage of the A/C motors at the 1/4 length point
impedance 1
phase A
0.5 phase B
• Phase A and C were 0 phase C
directly affected by the 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3
Speeds of the A/C motors at the end of the feeder
1.4 1.5
propagated to non- 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3
Sequence voltages of the 115 kV bus# 24160
1.4 1.5
31
Effects of load composition on FIDVR
phase A Casephase
1 B phase C
Case 1.5
Load composition
#
Speed (pu)
1
1 75% 1-Ф A/C motor, Case 1
25% constant impedance 0.5
Speed (pu)
1
3 70% 1-Ф A/C motor, Case 2
15% 3-Ф NEMA type B induction 0.5
motor, 15% constant impedance 0
4 60% 1-Ф A/C motor, 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
1.5
15% 3-Ф NEMA type B induction
Speed (pu)
motor, 25% constant impedance 1
Case 3
5 50% 1-Ф A/C motor, 0.5
25% 3-Ф NEMA type B induction 0
motor, 25% constant impedance 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
Speed (pu) 1.5
• Propagation of A/C motor stalling to 1
unfaulted phase is consistent across Case 4
0.5
a substantial range of load 0
0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
compositions 1.5
• The impacts of SLG faults close to
Speed (pu)
1
certain regions of the system with Case 5
0.5
high A/C penetration could be more 0
severe than perceived 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
Time (s)
32
Effects of point-on-wave (POW) on FIDVR
• POW effects on the
Phase A Phase B Phase C
occurrence and
1 1
Speed(pu)
Speed(pu)
evolution of FIDVR (a1) (a2)
0.5 0.5
are apparent,
0 0
based on the 0.5 1 1.5 2 0.5 1 1.5 2
1 1
Voltage(pu)
Voltage(pu)
differences in the (b1)
(b2)
0.5 0.5
response of the
A/C motors of bus 0
0.5 1 1.5 2
0
0.5 1 1.5 2
Positive Seq.Voltage(pu)
Positive Seq.Voltage(pu)
24160 bus 24151
bus 24160 1
bus 24151
bus 24160
1
• The POW when the 0.8
(c1)
0.8
(c2)
35
Thank you!
36