Over Matlab
Over Matlab
RcaITimcSimuIation Hardvarc-In-thc-Loop,HIL)VaIidation
Muhammad Shoaib AImas, Rujiroj LccIaruji, and Luigi \anIrctti
Royal Institute of Technology (KTH)
Teknikringen 33, SE-lOO 44, Stockholm, Sweden
Email: [email protected]@kth.se [email protected]
Abstract-Digital microprocessor based relays are currently
being utilized for safe, reliable and efcient operation of power
systems. The overcurrent protection relay is the most exten
sively used component to safeguard power systems from the
detrimental efects of faults. Wrong settings in overcurrent relay
parameters can lead to false tripping or even bypassing fault
conditions which can lead to a catastrophe. Therefore it is
important to validate the settings of power protection equipment
and to confrm its performance when subject to difrent fault
conditions. This paper presents the modeling of an overcurrent
relay in SimPowerSystems (MATLAB/Simulink). The overcurrent
relay has the features of instantaneous, time defnite and inverse
defnite minimum time (IDMT) characteristics. A power system
is modeled in SimPowerSystems and this overcurrent relay
model is incorporated in the test case. The overall model
is then simulated in real-time using Opal-RT's eMEGAsim
real-time simulator to analyze the relay's performance when
subjected to faults and with different characteristic settings
in the relay model. Finally Hardware-in-the-Loop validation of
the model is done by using the overcurrent protection feature
in Schweitzer Engineering Laboratories Relay SEL-487E. The
event reports generated by the SEL relays during Hardware-in
the-Loop testing are compared with the results obtained from
the standalone testing and software model to validate the model.
A
CRONYMS
CT Current Transformer
HIT. Hardware-in-the-Loop
RT Real-Time
RTS Real-Time Simulation/Simulators
SIL Software-in-the-Loop
SPS SimPowerSystems
VT Voltage Transformer
I.
I
NTRODUCTION
A power system is vulnerable to faults, either due to
natural disasters or by misoperation of the system due to
operators' negligence. This can result in permanent damage
to power system components leading to considerable costs
for their replacement and in longer disconnections of power
supply to customers, which is highly undesirable. This sets
a requirement for a power system to sustain faults, while
protection systems should minimize the damage in important
components and the effect of faults as much as possible. This
is achieved by using power system protection techniques and
methodologies [1]. Microprocessor-based relays have their
own algorithms for monitoring the power system through
current and voltage inputs fom CTs and VTs respectively.
Protection relays are thoroughly tested using relay test sets
to confrm their reliability and safe operation before commis
sioning them in a substation.
With the advent of RTS and compact reconfgurable I/O
controllers , large power systems can be simulated and their
behavior can be analyzed in both steady state and faulted
conditions. Based on this analysis, protection system and
power system controllers can be developed for more effcient,
reliable and safe operation of power systems. RTS are cur
rently being used for HIT. validation, developing algorithms
for adaptive protection, designing system integrity protection
schemes (SIPS), remedial action schemes (RAS) and other
applications [2].
In this paper, a detailed model for overcurrent relays is
provided. To assess the performance of this model a test
system is designed in SPS and is simulated in RT using Opal
RT's eMEGAsim real-time simulator. Once the SIL model
is certifed, it is further validated by HIL simulation using
the overcurrent protection function of Schweitzer Engineering
Laboratories (SEL) relay SEL-487E. The results obtained
by the SIT. and the HIL are compared. SPS does not have
a dedicated library for protection functions. Therefore it is
necessary to validate the model with HIL so that the validated
software model can be used for accurate representation of
power system protection components. Similar approach for
protection relay model validation has been demonstrated in
[3] [4] using diferent modeling platform and RTS.
The remainder of the paper is arranged as follows. Section II
ofers a literature review on overcurrent relays. Section III
presents the details of the overcurrent relay modelled in
SPS. Section IV focuses on the design of test system in
SPS and incorporation of proposed relay model in the test
case. SIL simulation and results are discussed in Section V,
while standalone testing and HIL validation is explained in
Section VI. Finally in Section VII, conclusions are drawn and
future work is outlined.
II.
O
VERCURRENT
R
ELAYS
As the name states, an overcurrent relay provides protection
against over currents. This relay uses current inputs fom a
CT and compares the measured values with preset values.
Figure 1 shows the logical representation of an overcurrent
relay. If the input current value exceeds the preset value,
the relay detects an overcurrent and issues a trip signal
to the breaker which opens its contact to disconnect the
protected equipment. When the relay detects a fault, the
condition is called fault pickup. The relay can send a trip
978-1-4673-2421-2/12/$31.00 2012 111 4789
lHuI LuIIEHI
|L|uVuE
I 5gH
LOnIIOI
Fig. !. Logical representation of Over-Current Relay
signal instantaneously after picking up the fault (in the case
of instantaneous overcurrent relays) or it can wait for a
specifc time before issuing a trip signal (in the case of time
overcurrent relays). This time delay is also known as the
operation time of the relay, and is computed by the relay
on the basis of the protection algorithm incorporated in the
microprocessor [5].
. Classifcation of Over-Current Relays
Overcurrent relays are classifed on the basis of their
operation time, in the following three categories:
1) Instantaneous Overcurrent Relay: These relays instan
taneously send a trip command to the breaker as soon as the
fault is detected (input current greater than the preset value).
They do not have any intentional time delay. They are usually
implemented close to the source where the fault current level
is very high and a small delay in operation of relay can cause
heavy damage to the equipment. So an instantaneous relay is
used there to detect and respond to a fault in few cycles.
2) Defnite Time Overcurrent Relay: This type of over
current relay is used for backup protection (e.g. back up
protection for transmission line where primary protection is
distance relay). If the distance relay does not detect a line
fault and does not trip the breaker, then after a specifc time
delay, the overcurrent relay will send a trip command to the
breaker. In this case, the overcurrent relay is time delayed by
a specifc time which is just greater than the normal operating
time of the distance relay plus the breaker operation time.
3) Inverse Defnite Minimum Time (IDMT) Overcurrent
Relay: This relay has an inverse time characteristic. This
means that the relay operating time is inversely proportional
to the fault current. If the fault current is higher, the operating
time will be lesser [6]. It can be graded for a very large range
of operating times and fault currents [7]. The characteristics
of an IDMT overcurrent relay depend on the type of standard
selected for the relay operation. These standards can be
ANSI, IEEE, lAC or user defned. The relay calculates the
operation time by using the characteristic curves and their
corresponding parameters [8]. Any of the above mentioned
standards can be used to implement a characteristic curve
for an overcurrent relay. The overcurrent relay will then
calculate the operation time corresponding to that particular
characteristic curve.
In accordance with IEC 60255 or BS142, the characteristics
of IDMT relays are represented with the following equation:
where:
C
J
Z TMS
,` 1
s
(1)
artI art2 art3
Fig. 2. Functional Block Diagram of an Overcurrent Relay
J. Relay operation time.
C: Constant for relay characteristic.
Is: Current Set point.
I: Current Input to the relay.
0: Constant Representing Inverse Time Type (0 ` 0).
TMS: Time Multiplier setting controls the relay tripping time.
By using appropriate TMS settings, the grading of a protection
network system can be achieved [9]. The range of TMS is
normally 0.1 to 1.0.Different types of curves can be obtained
by varying 0 and C [10]. Table I below shows values for 0
and C corresponding to each curve:
TABLE I
OIIILRLN11YPLS CI !NvLRSLCHARAC1LRIS1ICSCURvLS
Relay Characteristic Tpe C C
Standard Inverse 0.02 0.14
Very Inverse 1 13.5
Extremely Inverse 2 80
Long Inverse 1 120
III. OVERCURRENT
R
ELAY
M
ODELING
F
OR
R
EAL-
T
IME
S
IMULATION
The detailed model of an overcurrent relay is implemented
in the SimPowerSystems (MATLAB/Simulink) Toolbox. The
reason for using the SPS Toolbox is because it operates in the
Simulink environment and is a dedicated tool for modeling
and simulating power systems [II]. Another motivation in
using SPS is because it is compatible with the OPAL-RT
eMEGAsim real-time simulator [12] which is used to simulate
the power system model in RT and for HIL validation of the
model.
The block diagram presented in Figure 2 shows the three
major functional components of an overcurrent relay. These
components and their functionalities are discussed next.
. Part 1
In order to present an accurate model of the microprocessor
based overcurrent relay, the current inputs which are analog
signals are converted to digital signals by Analog to Digital
Conversion. In Simulink this can be achieved by using a Zero
Order Hold Block. This ZOH block is available in the Discrete
Library of Simulink.
D. Part 2
A digital flter is used to extract the fundamental signal. A
down sampler is used in order to reduce the sampling rate.
Down sampling helps avoiding anti-aliasing [13] efects. The
digital flter used here is a low pass digital FIR flter [14].
4790
pOwQU
Fig. 5. Test Case Model developed in SimPowerSystems
C Part 3
Once the RMS value of the current is obtained, this
current is fed into the relay protection algorithm block. This
block compares the current value with the pickup value. If
the input current exceeds the pickup value, the relay will
compute the operation time for this scenario taking into
consideration the characteristic curve (Standard Inverse, Very
Inverse, Extremely Inverse or Long Inverse) and sends a trip
signal once the operation time is elapsed.
Figure 3 in p. 7 shows the overall Simulink model of the
overcurrent relay.
The protection algorithm implemented in the model is
shown in Fig. 4, p. 7. The overcurrent model is masked to
a single subsystem (block) and a graphical user interface is
provided where the user can select the type of overcurrent
relay (instantaneous, time defnite or IDMT) and can further
select the different characteristic curves for IDMT overcurrent
relay.
IV.
T
EST
C
ASE IN
S
IMPOWER
S
YSTEMS
Figure. 5 shows the test case model designed in SPS to
couple the overcurrent relay model and investigate its charac
teristics when subjected to faults. The major components of
the test case are:
Three phase voltage source, 50 Hz, 1 kV phase voltage.
Transmission Line (7-section), 100 km
Three phase fault block to introduce single phase to
ground, phase to phase and three phase to ground faults.
Circuit breaker to disconnect load with trip signals fom
the overcurrent relay model.
Three phase series RLC load of 1 MW.
Simulation time step " 50 j sec.
In order to guarantee a safe, reliable and fast operation
of overcurrent relays, settings such as pickup value, time
multiplier settings, characteristic curves, and others, have
to be chosen carefully. Otherwise, the relay can either not
detect the faulted conditions (and will not trip) or send a
false tripping command. Figure 6 shows some of the general
steps which are followed in selecting the pickup value of the
overcurrent relay. Following these steps, fault analysis of the
test system was conducted. The results are shown in Fig. 7.
Table II shows the settings of the overcurrent relay model
which are based on the results fom fault analysis (Fig. 7)
CaclulaIeu|| CalculaIeShorI Consider1urn3aIioIor
'oadCurrenI CircuiICurrenI CurrenI1ransIormer
SelecI|nsIanIaneous/ SelecIickupvalueoyusinIherelaIion
DeIiniIeor|DM1
l.2 xMox|mum|ooJCurrent:P|cku
st n un o
ro
~vu:,+=es
~>.rn u.u
Ratio (Input Current to Relay f Pickup Value)
Fig. 8. Different IDMT characteristic curves for TMS 0.5. The vertical
line at ratio 7.72 is the region of interest in the test case
X
ti:igoa| ttom:imu+ole|a,
Voae|is:ootiguteatoooeottne
0igita|Oututottne :imu|atot
ti:igoa|iste:ei.eaatooeottne0igita|ioutot
tne :imu|atotwni:nis :ootiguteaiotne:imuol
Voae|tosnowttiiogiotne:ooso|e
0iitalOutput
and 0iital
|nputare
looped oack
Fig. 9. Steps for Real-Time SIL Validation of the Relay Model
........
....
_
.... ....
Fig. 10. Results for Relay Model with Standard Inverse Curve
of the simulator is 20 mA. In order to make sure that the
current level never exceeds this limit, the three phase current
signals are attenuated by a factor of 1000. As the maximum
fault current for this test case model is 3088 A (see Table II)
and the CT ratio is set to 400: 1. So the maximum fault
i:
4
x+ts+
......
0 1 3 4 6
......
..
...... ......
Fig. 11. Results for Relay Model with Very Inverse Curve
........ .......
4
....
0.6
0.4
1
0.
:
L 0
....
.
....
Fig. 12. Results for Relay Model with Extremely Inverse Curve
TABLE III
CCMPARISCN CI 1HLRLSUL1SIRCMPRCPCSL RLLAYMCLLSIL
vALIA1ICN vSILALCALCULA1ICNS
Summar of Resuls From Proposed Relay Model
(ault Applied at t=2sec)
MATLAB}Simulink
Design models for real
time simulation
The Analog outputs of the
Simulator are fed into the CT
Inputs of the SEL-487E
Fig. 14. Model-to-Data Work-Flow for HIL simulation. Solid Lines
indicate hardwired connections, dashed lines indicate digital data streams
over ethernet
simulator are confgured to change the breaker status in the
SPS model in real-time. Figure l3 shows the test case model
with additional blocks which are highlighted. The overall
workfow from building the model to simulating it in real
time and executing HIL validation is shown in Fig. 14.
In order to enable the overcurrent feature of the SEL-487E
SEL AcSELerator Quickset [18] is used. The characteristi
curve chosen for the validation is extremely inverse curve. The
same settings were incorporated in SEL-487E (Figure 15) as
in the software model (see Table II). A fault is introduced at
t=2sec. The ACSELERATOR Analytic Assistant [] 9] is used
Fig. 16. Results from Standalone test using Freja-300. Showing the tripping
tIme of 0.698 sec wIth TMS=0.5, pIckup value = 1 and characteristic curve
of extremely inverse. Scale is 1.93 Adiv
@ o:apntoc:\uses\msa|mos\oasHo\ 1u.cv
.U;
v,s
s,o
<,
o,o
1
.,s
s,o
!.
O IQ O OOD, = tOO
Fig. 17. Results from SEL-487E event report generated by HIL simulation.
The operating time calculated here is 34.5 cycles which corresponds to 0.69
sec
to display event report data as oscillography. In addition Freja
300 (Relay Test System) [20] is used to inject the secondary
currents to SEL-487E to confrm the tripping times obtained
fom HIL simulation. For standalone testing, the sofware
interface for Freja 300 i.e. Freja Win [21] is used to give
secondary injections of 0.98 A for the frst 2 seconds and then
a step increase in current to 7.72 A. The results are presented
in Fig. 16 showing that before the fault is applied, the current
input is equal to 0.98 A which is lesser than the pickup value
of 1 A. As the fault is applied at t=2 sec, the current rises
up to 7.72 A. At this point the relay picks up the fault and
its IDMT characteristic is activated which fnally generates
a trip signal after 0.698 sec. As the fault is applied at t=2
sec., and operating time of relay is 0.698 sec., so the tripping
time is equal to 2.698 sec. At t=2.698 sec, the relay D/O
changes its contacts and thus it opens the breaker in the SPS
model and the load is protected against overcurrent. Figure. 17
shows the oscillography fom the event report of SEL-487E
resulted fom HIL simulation. As soon as the relay trips,
SEL-487E generates an event report which can be viewed by
using ACSELERATOR Analytic Assistant. Figure. 18 shows
the operating time for SIL and HIL validations for diferent
characteristic curves.
4793
L
0
W
0
E
1.
1.
1.
8
4
Comparison ofSll and Hil Simulation for Different Characterlctlc Curves
= VEIyhVEI$E
~10h@1hvEr$E
+ OE8 L8 80h
W blL
0 H1
1.2
o
C
O
1
0.8
.I I.I05 I.I1
.CW3@ Ouif j
ll5 UlCCK CCm3lC5 llC
CCl3llC llmC CCmulCO Wlll
llC llmC Cl3tCO UyllCllmCl
Cl
O
l
-
Cl
.C
3@
lll
Cl
CWmlC xll3CllglmCCI 3ull
lmCl
Ouif
HCl3llC3l
>CCC
LCl3lCl
Very |nverse
YES
3|erationTime
>=Timer
1kll
Fig. 4. Protection Algorithm Implemented in the Overcurrent Relay Model
4795
>CCC
ll5 UlCCK lmlCmCl5llCm3llCm3llC3l
Cu3llC lC CCmulC
IL-487E Dm| . 2Ul W
|
Fig. 15. Settings of SEL-487E relay using AcSELerator Quickset. (a) Inverse time overcurrent element enabling with extremely inverse curve, Pickup value
of 1 and TMS of 0.5, (b) Trip logic for the IDMT confguration, (c) Confguring D/O of relay to change its contacts for IDMT trip (d) Enabling CT of relay,
(e) Enabling Event Recording, (f Enabling IDMT feature for Terminal U input of SEL-487
4796