Recent Advances in Systems Science

Research on Differential Protection Algorithms of Power Transformers

Using PSCAD/EMTDC Software
ADAM SMOLARCZYK, EMIL BARTOSIWICZ
Institute of Power Engineering
Warsaw University of Technology
ul. Koszykowa 75, 00-662 Warsaw
POLAND
[email protected], [email protected]
Abstract: - The paper describes proposed and PSCAD/EMTDC modeled test systems for testing the differential
relay algorithms used for protection of power transformers. The proposed systems can be used to test new
differential protection algorithms and real differential protection devices. The paper contains the description of
construction and operation principles of an example numerical differential relay for protection of power
transformer against the consequences of the internal faults. What is more, an example differential protection
device has been modeled in PSCAD/EMTDC software and the chosen results of its operating during simulated
faults are also presented in this paper.

Key-Words: - power system relays, current differential protection of power transformers, modeling of power
system relays, protection algorithms testing, testing of power system relays, PSCAD/EMTDC software
advantage of this type of modeling is, however, full
control of the modeled system and use of the only
necessary parameters for mapping the modeled
phenomenon.
In case of need for mapping the asymmetric
phenomena (phase-to-ground or phase-to-phase
short circuits, conductor openings or cut-offs, etc.)
in modeled parts of PS, more accurate models of PS
components are needed. There are many simulation
programs used for modeling of dynamic quicklyvarying (electromagnetic) and slowly-varying
(electromechanical) transients phenomena. One of
them is PSCAD/EMTDC software. [5]. This
program has well-equipped library of PS
components and elements of both the process of the
simulation and modeling elements of PS controls. It
is used to map both slowly- and quickly-varying
phenomena (for example to test the travelling-wave
fault location algorithms [6]).
Power system relays and operation algorithms
used in their protection functions must be tested
according to specific procedures. Extensive split of
the tests which the protection devices should be
subjected to is shown in [7]. Among the tests that
the relays should be subject to there are functional
performance tests and scheme performance tests.
These tests are related to validation of the relay
operation correctness in particular PS work
conditions and for the particular fault process in PS.
It can be done in such a way that the time-courses of

1 Introduction
Simulation software environments give the
opportunities to represent the entire power systems
(PS) or selected fragments thereof. They allow to
analyze: switching overvoltage, the correctness of
the design and operation of protection devices and
automatics, failures, electrical energy quality, drive
system
control,
FACTS
system
control.
MATLAB/Simulink environment is often used for
the study of control strategies or new protection
algorithms. In this program, particular elements of
the constructed part of PS (generators and control
systems, transmission lines, power transformers,
etc.) can be modeled independently or some readyto-use elements available in the SimPowerSystems
library (toolbox) can be used. For example, in [1] all
the elements of the PS for the test systems (singleand multi-machine) has been independently
modeled to check the fast generator resynchronization automatics, in [2] to verify the control
strategy of braking resistors for damping of power
swings, in [3] to check the control of FACTS
elements to dampen power swings and in [4] - to
check the control of power transmission through
direct current lines (HVDC). Unfortunately, singlehandedly modeling of power system components is
very time-consuming and usually allows to model
the elements used for mapping the phenomena of
symmetric character (power swing, symmetrical
faults, changes of working system load, etc.). The

ISBN: 978-960-474-314-8

137

Recent Advances in Systems Science

generated during modeling distinct types of system

faults.

electrical quantities just before and during the

disturbances can be obtained using a computer
simulation of the entire PS or corresponding
fragment thereof (e.g. using the PSCAD/EMTDC
software). Microprocessor tester forces the recorded
voltage and current time-courses by converting
digital signals to analog and enhancing it to the
proper level for analog inputs of tested relay. These
signals are injected into the relay and its operation is
observed by recording the pick-up and tripping
signals. In the case of systemic performance tests
the real time simulations are also used for checking
the power system relays. [8, 9].
In order to verify the operation correctness of the
differential relays for power transformer protection
and of the protection functions algorithms
implemented therein the microprocessor testers
software can be used. [10]. This software usually
gives an opportunity to validate the operating
characteristics of the differential function by using
the not-distorted sine waves of fundamental
frequency. For a more sophisticated tests of
differential relays and protection functions
algorithms implemented therein it is the appropriate
to use test-scheme models made in simulation
programs
and
voltage/current
time-courses

2 Structure of numerical differential

protection for power transformer
For the numerical differential protection, before a
direct comparison of differential currents Idiff and
biased/stabilized currents Ibias (for each phase), there
must be realized a mathematical processing of the
measured
current
signals.
Algorithm
for
measurement and adjustment of the measured
currents performed in the modern numerical
differential protection (for example for Siemens
devices) is shown in Fig. 1 [11]. Each module of the
differential relay is in charge for a particular
mathematical operation. It should be noted that all
the calculations are done relative to the reference
winding which for the Siemens differential relays
(e.g. 7UT512 [12]) is always the first winding
(WINDING 1 in Fig. 1). This means that the
elimination of the phase shift of currents is always
realized on the side of second winding.
(WINDING 2 in Fig. 1).

WINDING 1 (W1)
CT1
L1

WINDING 2 (W2)

Z1

Z2

I L1

I L1

I L2

I L2

I L3

* I L3

CT2
L1

I L3_CT1_NORM

I L1_CT1
I L2_CT1
I L3_CT1

Unit 5
I L1_CT2
I L2_CT2
I L3_CT2

L3

Unit 4
I L1_CT2
I L2_CT2
I L3_CT2

Unit 2

I L1_CT2_NORM
I L2_CT2_NORM
I L3_CT2_NORM

Currents normalizations

I L1_CT1_NORM
I L2_CT1_NORM

Unit 6

I0 elimination

I L3_CT1

Unit 3

I0 elimination

I L1_CT1
I L2_CT1

Currents normalizations

Unit 1

Vector group adaptation

L3

L2

Comparision
Idiff and Ibias

L2

I L1_CT2
I L2_CT2
I L3_CT2

Numerical differential relay

Fig. 1. The processing of current signals in differential relay; CT1 and CT2 current transformers, Unit 1 6
current-processing units; WINDING 1, WINDING 2 protected transformer windings; basis on [11]
Computational algorithm used in numerical
differential relay presented in Fig. 1 can be divided
into the following steps [11]:
- determination of amplitude matching factors on
the basis of rated currents of protected power
transformer (for lower- and higher-voltage
sides) and using rated currents of each current
transformer primary windings (Unit 1, Unit 2);
- the adjustment of phase currents amplitude,
having regard to the current ratio of current

ISBN: 978-960-474-314-8

138

transformers and amplitude matching factors

(Unit 1, Unit 2);
elimination
of
zero-sequence
current
symmetrical component I0 of normalized phase
currents. Zero-sequence current of the phase
currents is eliminated only when the protected
transformer winding is wye-connected and its
neutral point is solidly grounded (Unit 3, Unit
4);
vector group adaptation of the WINDING 2
phase currents to phase shift of the reference

Recent Advances in Systems Science

- a scheme of two-winding transformer

(110 kV/15 kV) with YNd11 vector group scheme 3,
- a scheme of auto-transformer (230 kV/120 kV) scheme 4,
- a scheme of two-winding transformer
(22 kV/250 kV) with Dyn11 vector group scheme 5,
- a scheme of two-winding transformer
(22 kV/250 kV) with Dyn5 vector group, build of
three single-phase units - scheme 6.

winding (WINDING 1) of the transformer, by

using additional phase shift (elimination of
phase shift basing on the transformer vector
group) via programmable coefficients matrix
(Unit 5);
computing of differential currents Idiff and
biased (stabilized) currents Ibias for each phase
separately according to
the
relations
Idiff =II+III and Ibias =II+III(formulas for
Siemens relays) wherein the L1 phase in
accordance with Fig. 1 current II = I*L1-CT1 and
**
current III = I L1-CT2 (Unit 6);
checking the position of operating point
described by Idiff and Ibias currents (for each
phase separately) with respect to the pick-up
characteristic Idiff = f(Ibias) and decision (after
taking into account the conditions associated
with the 2nd and 5th harmonic content in the
Idiff current) with tripping the pulse for opening
a circuit breakers located on both sides of the
protected transformer (Unit 6).

Classical approach of transformer modeling has

been used in schemes 1, 2, 3 and 4. This approach
allows the reproduction of windings placed on the
same column, which does not take into account the
interactions between the individual phases [5]. For
this kind of models, three-phase two-winding
transformer is equivalent to a combination of singlephase two-winding transformers. In scheme 5 there
has been exploited the model of 5-limb transformer,
that uses the unified magnetic equivalent circuit.
Transformer model of this type as distinct from
conventional transformer models, does take into
consideration the transformer core geometry, the
magnetic coupling between the phase windings and
internal coupling of each phase winding [5]. In
system 6, due to a lack of available proper model in
PSCAD/EMTDC main library, the three-phase twowinding transformer with Dyn5 vector group has
been built of three single-phase units. Depending on
the type of modeled transformer, and thus its
application, the proper feeding (one- or two-way)
has been envisaged.
Fig. 2 shows the appearance of a typical primary
circuit for the project of power transformer with
YNd11 vector group (scheme 3). In addition to the
elements shown in the figure, in the simulation
project there are blocks which contain the following
elements:
- elements of modeled differential relay,
- diagrams of analog signals (currents and voltages
on higher- and lower-voltage sides of the
transformer and signals made available in the
transformer model) and binary signals (states of
switches),
- elements related to the control (choice of location
and type of fault, breakers control, control of
saving the analog variables time-courses to
COMTRADE-format files).

3 Modeled schemes for testing the

transformer differential protection
algorithms
In [13] thesis there are presented the descriptions of
simulation systems projects containing transformers
and modeled differential relays for their protection.
These projects, made in the PSCAD/EMTDC
software, reproduce the power transformers with
different vector groups and rated powers, lines (as
PI-type two-terminal-pair networks), power source
and the load on lower- and higher-voltage sides. For
the modeled systems there are possibilities of
switches control in both a manual and automatic (by
simulating differential relay) way.
In addition, depending on the transformer vector
group and its rated power, there have been used
proper algorithms for appropriate operation of
modeled transformer differential protection in
simulation projects. The signals injected to the
differential relay may be brought from actual or
ideal current transformers (depending on the
requirements). In separate simulation projects,
following test schemes have been prepared:
- a scheme of two-winding transformer
(15 kV/6 kV) with Yy0 vector group - scheme 1,
- a scheme of two-winding transformer
(400 kV/110 kV) with YNy0 vector group scheme 2,

ISBN: 978-960-474-314-8

139

Recent Advances in Systems Science

Fig. 2. Primary circuit of simulation model using the 110 kV/15 kV transformer with YNd11 vector group scheme 3; derivate from [13]
thus the differential function correctly eliminates
the zero-sequence component of the current),
- switching on the unloaded transformer (whether
there is no unnecessary pick-up of differential
function during the inrush of magnetizing current
and appearance of a high content of 2nd harmonic
in differential current),
- transformer overfluxing (whether there is no
unnecessary pick-up of differential function due to
the large rise of supply voltage and the emergence
of a high content of 5th harmonic in the
differential current).

In all the modeled systems with transformers it is

possible to visualize all the courses of currents and
voltages on the higher- and lower-voltage sides and
in the fault point. The observed voltage and current
time-courses at relaying points, on the secondary
side of current and voltage transformers, can be
saved as the COMTRADE format.
As a result this time-courses can be retraced by
microprocessor testers, and thus be used to test the
real differential relays (functional performance tests
[7]).
Additionally, there is a possibility in modeled
systems to observe variation of differential current
(Idiff) as a function of stabilized current (Ibias) for the
A, B, C phases, occurring during disturbances, with
respect to the beginning characteristics Idiff = f(Ibias)
for modeled differential protections.
With the modeled test schemes can be checked
(within the functional performance tests) the
operation correctness of the differential protection
functions in case of:
- internal and external single and multi-phase faults,
with current transformers saturation absence and
for various pre-load of transformer (the impact of
the appearance of non-periodic current component
in fault current),
- metallic and resistive internal faults,
- various level of saturation of the CTs (on one or
both sides of the transformer) and various pre-load
of transformer, in the case of external faults
(whether there is no unnecessary pick-up of
differential function),
- external single-phase faults on the transformer side
of a grounded star-point (whether there is no
unnecessary pick-up of differential function, and

ISBN: 978-960-474-314-8

4 A scheme of modeled transformer

differential protection
Differential relays, designed to protect the
transformers from the schemes described in Sec. 3
of this paper, have been modeled in thesis [13]. For
modeling the protection device, shown in Fig. 1
block diagram of signal processing in numerical
differential protection has been used. In addition to
the stabilized differential function, the unstabilized
differential function has been modeled as operating
in case of exceeding the high value of differential
current (e.g. 10In).
The protection modeling process has been
focused on the protection operation logic, not on
digital signal processing in order to computing the
1st, 2nd and 5th harmonics of phase currents. For
the calculation of particular harmonics of phase
currents (phasors), the On-Line Frequency Scanner
module (from PSCAD/EMTDC main library, [5])
has been used. The description of direct calculation
method for measurement signals (current phasors),

140

Recent Advances in Systems Science

states: closed breakers CB_SEE1, CB_GN, CB_DN

and CB_OBC2 (marked with red squares) and
opened breakers CB_OBC1 and CB_SEE2 (marked
with green squares) - what is equivalent to one-way
supply of transformer on its higher-voltage side and
applying the load on its lower-voltage side.
Metallic B-C fault has been simulated (internal,
on the higher-voltage transformer side). Timecourses of differential (I_DIFF_H1) and stabilized
currents (I_BIAS_H1) for the A, B, C phases of
fundamental frequency, including time-courses of
the tripping signals (OW_A_H1, OW_B_H1 and
OW_C_H1) modeled relay, are shown in Fig. 4.
During the two-phase B-C fault, differential currents
of these phases (IB_DIFF_H1, IC_DIFF_H1) are
very similar (47.46 A) and nearly twice as high as
biased currents (IB_BIAS_H1, IC_BIAS_H1), which
amount to 24.16 A and 23.77 A. Differential current
value for phase A (IA_DIFF_H1) is close to zero.

without using the On-Line Frequency Scanner

module, may be found in [14].
Fig. 2 shows the location of differential relay
model in the simulation project. Block diagram of
the modeled relay includes signal processing shown
in Fig. 1, and due to the limited volume the article
does not show its construction in PSCAD/EMTDC.
Fig. 3 shows the operation (stabilized)
characteristics Idiff = f(Ibias) of modeled relay,
implemented by using the Dual Slope Current
Differential Relay module from PSCAD/EMTDC
differential function.
Idiff
Operating area

k1

k2
Restrain area

IS1

IS2
Ibias
Fig. 3. Stabilized characteristic of modeled relay,
implemented by using the Dual Slope Current
Differential Relay module; IS1 beginning value of
differential current; IS2 inflection-point value of
biased current ; k1 slope factor of first segment; k2
slope factor of second segment

Formulas for differential current Idiff and

stabilized current Ibias (separately for each phase),
used in Dual Slope Current Differential Relay
module and in modeled relay, have the following
form:
I diff = I I + I II
(1)

I bias = ( I I + I II ) / 2

Fig. 4. Time-courses of differential and stabilized

currents of fundamental frequency (for each phase)
and binary tripping signals of modeled relay for
internal metallic B-C fault on the higher-voltage
side of transformer

(2)

where: II, III transformed phase (A, B, C) currents

for the transformer sides of higher- and lowervoltage.

5 Example results of functional

performance tests for modeled
differential protection

Relay operating
characteristic

This section presents selected results of disturbance

simulations for modeled test schemes. The research,
described in greater detail in [13], have been aimed
to monitor and analyze the signals that occur during
chosen disturbances related to transformers and to
verify the operation of modeled transformer
differential protection.
The scheme simulation project is illustrated in
Fig. 2. In the figure are shown the default breakers

ISBN: 978-960-474-314-8

Fig. 5. Relay operating characteristic and calculated

currents of modeled relay for internal metallic B-C
fault on the higher-voltage side of transformer

141

Recent Advances in Systems Science

Fig. 5 illustrates the relay operating characteristic

of modeled differential relay and differential current
as a function of biased current (for each phase
separately). It shows exceeding of the relay
operating characteristic by the working points,
computed basing on differential and biased
(stabilized) currents of faulted phases (B, C), which
causes pick-up of modeled differential relay and
sending the tripping signals (OW_C_H1 and
OW_B_H1), shown in Fig. 4. The amount of
differential current in phase A does not reach the
boot value (Fig. 5), so the tripping signal for this
phase (OW_A_H1) is set to zero (Fig. 4). So
operation of the modeled differential protection is
consistent with the expectations.

References:
[1]
J.Machowski,
A.Smolarczyk,
J.Bialek,
Coordinated fast valving and excitation control as
an effective method to enhance transient stability
after severe faults, CIGRE Symposium Working
Plant and Systems Harder, London, 7-9 June
1999.
[2]
J.Machowski,
A.Smolarczyk,
J.W.Bialek,
Damping of power swings by control of braking
resistors, International Journal of Electrical Power
& Energy Systems 23 (2001), pp.539-548.
[3] S.Robak, Selection of the structure, type and
parameters of robust FACTS-based stabilizer.
Electrical Review, ISSN 0033-2097, Vol. 85, Issue:
8, 2009, pp. 76-79.
[4] .Nogal, Enhancing Stability of a Weakly
Connected AC Power System By Supplementary
Control of Embedded HVDC Lines, Internatio-nal
Conference on Environment and Electrical
Engineering, Wroclaw, Poland, 5 - 8 May 2013,
pp. 79-86, IEEE Conference Record #21360, IEEE
Catalog Number: CFP1351I-CDR.
[5] Manitoba HVDC Research Centre, PSCAD v.4.2.1
Online Help.
[6] K.Glik, D.Rasolomampionona, Laboratory Test
Results of Travelling Wave Fault Location
Scheme, PowerTech, Grenoble France, 16-20 June
2013.
[7] CIGRE: Analysis and guidelines for testing
numerical protection schemes, CIGRE Report No
159, Working Group 34.10, August 2000.
[8] RTDS Technologies Inc., Real Time Digital
Simulators, Available on-line: www.rtds.com.
[9] A.Smolarczyk, A modern laboratory test bench for
protection relays testing in real-time phenomena,
(in Polish), Electrical Review, Vol.: 88, Issue: 3A,
pp. 142-147, Publ.: 2012.
[10] OMICRON electronics: Test Universe Software,
Available on-line (2013): www.omicron.at.
[11] G.Ziegler, Numerical Differential Protection,
Principles and Applications, Publicis Corporate
Publishing, Erlangen 2005.
[12] Siemens, Numerical Differential Protection Relay
7UT51 v3.0, Instruction manual Order No
C53000-G1176-C99-6.
[13] K.Pyrzak, Application of PSCAD/ETMDC
software for modeling of transformers and
differential protective relays, The MSc thesis, (in
Polish), Warsaw 2011.

6 Summary
The PSCAD/EMTDC software can be successfully
used to modeling the disturbances in schemes with
transformers and differential relays protecting
transformer from effects of this disturbances
occurrence. There can be created the simulation
projects of schemes with power transformers, which
allow to analyze various types of disturbances that
may occur in this schemes, such as internal and
external faults, switching on the unladed
transformers, voltage spikes on the transformer
terminals, faults with current transformer saturation,
etc.
A lot of tests have been made for modeled test
schemes. Example results of the tests, described in
the article, have shown the usefulness of
PSCAD/EMTDC for verifying the operation
correctness of differential protections for power
transformers. It is important that the calculations are
in this program performed in the time domain.
Thanks to this, instantaneous values of phase
voltages and currents may be injected into the
modeled protection device. Modeled in the thesis
[13] test schemes can be used for testing the
protection algorithms of real protection devices
(thanks to the possibility of saving the time-courses
of analog and binary signal to COMTRADE files)
and for designing the new differential protection
solutions.
Modeled in the thesis [13] differential relays, for
protecting power transformers of different rated
powers and different vector groups, and example
tests results of its operating have proved that
PSCAD/EMTDC can be used for testing new
protection function algorithms and verifying
operating correctness of the existing algorithms (e.g.
for educational or training purposes).

ISBN: 978-960-474-314-8

[14] M.Sztykiel, C.L. Bak, Line Differential Protection

Scheme Modelling for Underground 420 kV Cable
Systems, EMTDC/PSCAD Relays Modeling,
OMICRON electronics, International Protection
Testing Symposium (IPTS), Brand 2011.

142