CT Saturation
CT Saturation
CT Saturation
Abstract—Current transformer (CT) saturation, while a fairly understanding of CT saturation must extend beyond the
common occurrence in protection systems, is not often clearly knee point.
understood by protective relay engineers. This paper forgoes the
usual physics equations to describe how CTs saturate in a simple
and intuitive way. We explain the differences between II. CT SATURATION THEORY
symmetrical and asymmetrical saturation and how remanence To understand CT saturation, it is very important to
accumulates in the core of a CT. We then describe the CT understand the basic concept of how CTs work and what is
equivalent circuit and how it results in the familiar CT excitation actually happening when they saturate. This section describes
graph. ANSI ratings of CTs are explained, and we show how to
analyze the performance of CTs using simple equations and tools. what happens in the core of a CT during symmetrical saturation,
Finally, we explain how CT saturation can affect relay operation asymmetrical saturation, and remanence. It then explains how
and show how to detect CT saturation in protective relay event this core activity corresponds to the CT equivalent circuit,
reports. Real-world event reports are presented where correct ANSI voltage ratings, and the familiar CT excitation graph.
relay operation was compromised as a result of incorrect current
values from saturated CTs. A. How CTs Work
In its simplest form, a CT consists of two sets of wire
I. INTRODUCTION windings around an iron core, as shown in Fig. 1. The concept
Current transformer (CT) saturation is not a new topic, and is the same for a window or bushing CT, which consists of a
there have been many papers, books, application guides, and secondary winding around a core, with the primary winding
tutorials written on the subject. Sorting through this vast array being the primary conductor that passes through it.
of information to piece together a complete understanding of Transformers work based on the principle of electromagnetic
the topic is a time-consuming task and may not be realistic with induction. This principle states that an alternating magnetic flux
the schedules and demands placed on many practicing in the presence of a loop of wire induces a voltage across that
engineers. Because of this, engineers’ level of understanding is loop. Magnetic flux is simply the amount of magnetic field
often limited to the familiar CT excitation graph. The following passing through a material such as a transformer core.
is a list of common questions related to CT saturation: When alternating current IP flows in the primary winding of
Why does a CT saturate (Section II, Subsection B)? a transformer, it generates an alternating magnetic field H,
What is remanence, and do I have to worry about it
(Section II, Subsection C)? the transformer core. This alternating magnetic flux passes
What does it mean when a CT is a C800 (Section II, through the secondary winding. What happens next depends on
Subsection F)? the load connected to the secondary winding.
I have a C800 multiratio CT tapped at 400/5. Is it still If the secondary is connected to a burden, the alternating
a C800 (Section II, Subsection F)? magnetic flux in the core induces an alternating voltage VS
across the secondary winding. This causes a corresponding
How do I make sure my CT will not saturate for my
alternating current IS to flow in the secondary winding. The
worst-case fault current (Section III)?
alternating current in the secondary creates its own alternating
The knee point determines the saturation voltage of a
magnetic field and alternating magnetic flux that oppose those
CT, right (Section II, Subsection E)?
created by the primary winding. These primary and secondary
Will saturated CTs cause my relay to misoperate?
fluxes cancel, leaving a negligible amount of net flux in the
What if they just saturate a little bit (Section V)?
core. This occurs until the core becomes saturated.
After a misoperation, how do I know if CT saturation If the secondary is open-circuited, the alternating magnetic
was a cause (Section IV)? flux in the core induces a very high alternating voltage VS
Can modern relays prevent misoperations due to CT across the secondary winding. VS remains on the terminals with
saturation (Section VI)? no secondary current flowing, which is why it is very dangerous
The goal of this paper is to explain CT saturation to the to open-circuit an in-service CT. Because IS cannot flow, it
protective relay engineer and to answer these questions in a cannot create an opposing magnetic field, leaving a net flux in
clear and practical way. As this paper demonstrates, a proper the core equal to the flux created by the primary current.
2
The example in Fig. 3 shows the primary current decreasing saturates. This limit is defined by a symmetrical sine wave with
in magnitude every cycle. The point of this is to show that if a fixed voltage magnitude and fixed area under the curve in
primary current magnitude decreases, the CT is saturated for both the positive and negative directions. As long as the actual
less time. The lower magnitude in the second cycle of Fig. 3 CT waveform does not exceed this positive or negative volt-
generates a weaker magnetic field, requiring less flux density time area, the CT will not saturate. Consider the dc offset of the
to replicate the current correctly. Because fewer magnets are asymmetrical current in Fig. 4. This dc offset will result in an
used, IS reliably replicates IP for a longer time until all the accumulating positive volt-time area that eventually reaches the
magnets are aligned. In the third cycle, the magnitude of IP has maximum that the CT can handle at Point d, where saturation
been lowered to the point that the CT does not saturate and occurs.
replicates current correctly the entire time.
IS
IP
Time
a b c d
D. CT Equivalent Circuit when the ZE impedance is large (during normal conditions) and
We can represent the behavior of the C-class CT shown in large when the ZE impedance is small (during saturation). The
Fig. 8a with the simplified equivalent circuit shown in Fig. 8b. current measured by the load (IS) is what is left of the ratio
current (IST) after the error current (IE) is lost through the
(a) Secondary
magnetizing impedance (ZE). Therefore, the current the relay
sees is accurate during times of high ZE impedance (low error
current), and it is much lower than expected during times of low
ZE impedance (high error current). This behavior gives us the
Primary
classic saturation waveform shown in Fig. 3.
Fig. 9, adapted from [6], explains the saturation process in
the time domain. The graph in Quadrant III shows the primary
Fully Distributed Windings current IP being applied to the CT. The graph in Quadrant II
IST RS IS
shows the flux density over time. At time t0, when the CT is
(b)
first energized, the flux density starts to increase from whatever
IE
remanent value was left in the CT at the last de-energization.
Flux density continues to rise in the CT during the positive half
IP VE ZE VS ZB
cycle and decreases during the negative half cycle. Time t3 is
the point of maximum flux density, when the CT is fully
saturated. The graph in Quadrant IV shows the excitation
N1 N2
(error) current over time. We can see that the highest excitation
Fig. 8. Class C CT (a) and Simplified CT equivalent circuit (b) current occurs at the point of maximum flux density. The graph
IP (the primary current), VS (the secondary CT terminal in Quadrant I shows flux density versus excitation current and
voltage), and IS (the secondary load current) are is made up of the intersection of points between the graph in
self-explanatory. N1 and N2 represent the turns ratio of the CT, Quadrant II and the graph in Quadrant IV. The graph in
and the ratio current IST is the primary current divided by this Quadrant I is not time-dependent, but simply shows the amount
ratio. The circuit connected to the secondary of the CT is the of excitation current that is measured for a given flux density.
burden represented by ZB and includes the impedance of the This graph shows that when the flux density is low, the
relay and CT leads. The resistance RS represents the secondary excitation (error) current is low (because ZE is a high value). As
winding resistance of the CT. The magnetizing behavior of the flux density increases and the CT approaches saturation, the
core of the CT is represented by a varying reactance, ZE. This excitation current increases substantially (because ZE is a lower
impedance changes as the CT saturates, resulting in varying value), meaning less current will be available to the burden.
secondary excitation voltage VE and varying exciting branch Notice that the graph in Quadrant I corresponds directly to the
current IE. (Although ZE technically has a real and reactive B-H curve in Fig. 2, where the excitation current IE is
component, most practical CT analysis is close enough when proportional to the magnetic field strength H.
neglecting the angles of the impedances.) The flux density in Fig. 9 is directly proportional to both the
In the CT equivalent circuit in Fig. 8, the constant current flux and the excitation voltage VE. We can therefore interpret
source IP drives the total secondary current IST as determined by the graph in Quadrant I as the instantaneous relationship
the turns ratio. IST is the current that is measured by the load if between VE and IE, and use it to understand the changing
the CT is ideal. However, we must take into account the impedance of the CT excitation branch (ZE). The slope of the
magnetizing behavior of the CT core that is represented by the curve is the ratio of VE/IE, which Ohm’s law allows us to
varying impedance ZE. Because magnetization is nonlinear, characterize as ZE. We can easily visualize that the curve is
different values of impedance (ZE) must be used for different made up of two linear sections and that the knee point is where
states of operation (normal conditions vs. saturated conditions). the two linear sections transition. The steep part of the curve,
We can think of ZE as a varying impedance that changes with where the ratio of VE/IE is high (ZE is high), represents where
flux density. It has a very high value under normal conditions the core is not saturated. This is called the iron-core reactance
and a very low value (basically a short-circuit) when the CT is region of the curve. The flat part of the curve, where the ratio
saturated. The excitation current (IE) is the current that is lost of VE/IE is low (ZE is low), represents where the core is
through the magnetizing branch and corresponds to the CT saturated. This is called the air-core reactance region of the
error. We can also refer to this as “error current.” IE is small curve because the saturated iron core has no better permeability
than air.
6
2) Saturation
The area of saturation is along the low linear slope portion
of the excitation graph, where a small change in VE results in a
very large change in IE. We select a random point on the line in
this region as (IE E
the equivalent circuit, as shown in Fig. 12. We then solve for
ZE If we connect a
burden of 4 , we can solve for IS ST We
Fig. 10. Example excitation graph for a C400 CT [7] see that because the error current is so large, the ratio current
(IST) is not very close to the current the load actually sees (IS),
meaning the CT is behaving poorly.
7
increase in current delivered to the load. The definition also of the CT has a much lower terminal voltage at 10 A of error
states that assuming a 5 A nominal CT, 100 A is delivered to a current than the full 2000/5 tap. For multiratio CTs, the voltage
standard burden when VS matches the voltage rating of the CT. rating of a tap is directly proportional to the CT ratio
IEEE defines standard burdens for relaying CTs as 1, 2, 4, and corresponding to the tap divided by the full CT ratio (assuming
8 . Consequently, the standard voltage ratings for 5 A nominal the windings are fully distributed around the core). For
relaying CTs are C100, C200, C400, and C800 (20 times the example, if a C400, 1200/5 CT is operating on a 600/5 tap, the
5 A rated secondary current) [9]. For example, a C400 accuracy voltage rating at 600/5 is calculated as shown in (2).
class on a 5 A CT means that the ratio error will not exceed 600 / 5
10 percent for any current up to 100 A secondary (20 times the TAP
(2)
1200 / 5
standard burden. This means that the CT supports a voltage of 200 V at 20
From this definition, the voltage rating of a CT defines times rated secondary current of 100 A, which results in a
voltage VS, burden ZB, and current IS in the example shown in standard burden of 2 .
Fig. 13 for a C400, 5 A nominal CT. The definition also tells us A final note on CT ratings and the knee point: notice that the
that at this voltage, the ratio error (IE / IS) is a maximum of knee point of the graph in Fig. 10 is VE
10 percent. If we assume this worst case, then IE must equal the ANSI voltage rating of the CT (C400). In fact, 400 V is not
10 A. We can then back-calculate to find IST even on the high linear slope portion of the graph. Recall again
the primary current divided by the turns ratio. If we go above that the knee point only shows where the CT behavior starts to
this, we will start getting more than 10 percent error on the become nonlinear, while the ANSI rating defines the point of
output, and the current the connected load reads will not be 10 percent error. The ANSI rating, not the knee point, defines
reliable. Similarly, increasing the burden ZB to a value greater the practical threshold for saturation. Because the knee-point
than the standard burden will also result in an increased error. voltage is normally 46 percent of the saturation voltage, a
When applying CTs, a simple rule of thumb to avoid a popular rule of thumb to avoid saturation is to ensure that the
symmetrical saturation error over 10 percent is to use a burden ANSI rating is twice the terminal voltage developed by the
equal to or less than the standard burden and ensure that the maximum fault current [8]. This ensures operation near the
worst-case secondary fault current is less than 100 A. knee point for the maximum symmetrical fault current.
IST = 110 A RS IS = 100 A
III. ANALYZING CT PERFORMANCE
IE = 10 A
When a CT saturates, it provides distorted information to the
IP VS = 400 V ZB = 4
connected relay, which can cause the relay to behave
unexpectedly. When selecting a CT for a protective relay
application, it is important to check if the CT saturates for the
levels of fault current that will occur on the system. This section
Fig. 13. Example for a C400, 5 A nominal CT explains an equation that can be used to determine if a CT will
saturate for given fault conditions.
IEEE defines common CT voltage ratings as C100, C200,
C400, and C800. Because the ratings are defined at specific A. CT Steady-State Performance Analysis
voltage levels, not all CTs with the same voltage rating are Equation (3) is the criterion to avoid CT saturation for
created equal. For example, a CT with a terminal voltage of symmetrical fault currents. Real-world faults are rarely
810 V and a CT with a terminal voltage of 1,000 V would both symmetrical, so this equation should never be used practically
be classified as C800. The CT with the excitation graph shown to analyze CT performance. It is explained here and derived in
in Fig. 10 has an excitation voltage of 496 V at an error current Appendix B as a basis for the more practical equation in
of 10 A. To get the terminal voltage VS, we subtract the voltage Section III, Subsection B. Equation (3) can be used to
drop across the secondary winding resistance (RS determine the maximum allowable fault current for a given
the excitation voltage of 496 V. If the resulting terminal voltage burden or maximum allowable burden for a given fault current.
is over 400 V, we classify the CT as a C400 CT [2]. I FAULT ZB R S
It is important to note that these voltage ratings only apply (3)
IPRI ZB STD R S
to the full winding ratio of the CT, and tapping down a CT
reduces its accuracy. Because both windings are distributed The terms in (3) are defined as follows:
around the same core, they are each subject to the same flux and IFAULT is the maximum fault current in primary amperes
the voltage induced across every turn of any winding is the for a given fault.
same. That is, the maximum volts per turn the CT can support IPRI is the primary current rating of the CT (e.g., for a
without symmetrical saturation is fixed by the cross-sectional 2000/5 CT, IPRI is 2,000.)
area of the core. If we choose to tap a CT at anything other than ZB is the actual burden of the CT’s secondary circuit.
its full winding, there will be less voltage available at the RS is the internal resistance of the CT secondary winding.
terminals because not all of the available turns are connected to
ZB STD is the standard burden of the CT (e.g., for a C800
the load. In Fig. 10, for example, we can see that the 300/5 tap CT, the ZB STD .
9
ZB includes both the impedance of the connected relay and in less saturation. Therefore, a value over 20 in this equation
the impedance of the leads from the CT to the relay. does not mean that the CT will saturate for that fault current—
Microprocessor-based relays have a negligible burden, but just that it could, assuming the fault current had the worst-case
electromechanical relays do not. The impedance of the leads is dc offset.
determined by the wire gauge, the length of the leads, how the
C. What About Remanence?
CTs are connected, and the fault type. #10 AWG copper wire
has a resistance of 0.9989 per 1,000 feet. Reference [10] There is no way to predict the value of remanence that may
shows how to calculate the resistance of other wire gauges. exist in a CT at a given instant in time. As described in
When using wye-connected CTs and performing the calculation Section II, remanence can either impair or improve the CT
for a single-phase-to-ground fault, a two-way lead burden behavior for a given fault. The worst case is when a fault
should be used because the fault current has to make a full loop produces flux of the same sign as the remanence. In this case,
through the CT circuit (phase and neutral) and back to the relay. the flux change required to saturate the CT equals the difference
When using wye-connected CTs and performing the calculation between the core saturation flux and the remanence.
for a three-phase fault, a one-way lead burden should be used Regardless of how much remanence exists in a CT, recall
because the fault current for a single phase sums with the other from Section II that the effect remanence has on saturation only
two phases at the neutral point of the wye connection and no lasts about half a cycle. Modifications to (5) have been made to
current will flow through the neutral. This is derived in [5], account for the effects of remanence by derating the CT, but
which also shows how to calculate the burden when using delta- they are only valid for the first half cycle of the fault [2]. Using
connected CTs. these modifications to size CTs is unnecessarily conservative
RS, the secondary winding resistance of the CT, is specified and often yields impractically large CTs. Be aware that because
on the CT data sheet. The RS value is directly related to the of the random nature of remanence, even a CT that was properly
resistance of the wire and the number of turns and is therefore sized using (5) may still saturate momentarily. This saturation
dependent on the tap for multiratio CTs. will be short-lived and should have little effect on the
ZB STD is the standard burden of the CT and must be adjusted performance of standard relay algorithms.
if the CT is not tapped at its full ratio. In these instances, the Despite the best intentions, there are certain times when the
standard burden must first be multiplied by the ratio of the tap use of (5) to select CTs that will never saturate becomes an
rating to the full rating. For example, for a C800, 2000/5 CT exercise in futility. One such example is on applications near a
tapped at 1200/5, the standard burden is shown in (4). The generator bus. Here, the X/R ratios and fault current magnitudes
resulting C- are extremely high, and it is usually impossible to avoid
asymmetrical saturation even with the largest CTs.
800 1, 200 Reference [2] recommends selecting the highest practical rating
B STD (4)
for the CTs and always matching the CTs on the terminal and
B. CT Transient Performance Analysis neutral sides of the generator to each other (same ANSI voltage
rating, CT ratio, and connected burden). This way, even though
Section II shows that symmetrical fault currents are not the we know the CTs will saturate, they will saturate in the same
only risk for CT saturation. Fault currents with an exponentially way and the saturated waveforms will cancel each other out in
decaying dc offset, caused by the X/R ratio of the system, can the differential calculation. Because not all CTs with the same
produce significant CT saturation. Selecting CTs based only on voltage rating are created equal, it is important to make sure that
symmetrical fault current is not advised because it ignores the
both CTs have the same excitation curve, knee-point voltage,
risk of heavy CT saturation when the fault current includes dc
and terminal voltage at 10 A of excitation current. See
offset. Section III, Subsection D for an example illustrating this point.
To account for dc offset, we can improve (3) by including
Another example of not being able to avoid saturation is in
an extra X/R term, resulting in (5) [5]. This extra term is the
low- and medium-voltage switchgear applications, where fault
X/R ratio of the entire Thevenin equivalent of the system, from
currents are high and space requirements force the use of small
the source to the fault, through the particular CT. Equation (5)
and poorly-rated CTs. Refer to [11] and [12] for selection
only considers dc offset from sinusoidal fault current, not offset criteria for low-voltage applications. For more information on
from nonsinusoidal sources like transformer inrush.
selecting CTs for various protection applications, see [1], [2],
Appendix C gives an example of using (5) to analyze
[7], [10], [13], and [14].
asymmetrical saturation in a CT.
IFAULT ZB R S X
D. Transient Performance of Different CTs With the
(5) Same Rating
IPRI ZB STD R S R
To illustrate the point discussed in the previous section about
The initial magnitude of the dc offset that occurs is trying to make sure that the terminal and neutral-end CTs on a
determined by the fault incidence angle, and the time it takes generator differential have the same excitation curve, we can
for the dc offset to decay is determined by the X/R ratio of the look at two different C800 CTs manufactured to two different
system. Equation (5) is quite conservative and assumes a worst- designs. Fig. 14 shows two CTs that are both classified as C800
case, completely offset waveform, which is not often the case because
in the real world. Any fault other than the worst case will result 100 A secondary current. The construction of the two CTs,
10
18, 000 A 1.2 1 Fig. 16 shows the calculated operate and restraint values for
(6) a differential relay connected to these two CTs. We can see that
2,000 A 8 1
the false differential current reaches over 25 A secondary for
We get the following for the CT in Fig. 14b (larger RS): almost 1.5 cycles. This relay would have to have a fairly high
18, 000 A 1.2 2 slope ratio of 40 percent to ride through the transient differential
(7) current from these two CTs with the same rating.
2, 000 A 8 2
100
Even though the CT in Fig. 14b has more iron (VE is higher), 90
transiently, (5) predicts that it will saturate more severely. 80
70
IST = 110 A RS = 1 IS = 100 A Restraint
(a) 60
IE = 10 A 50
40
30
IP VE = 900 V VS = 800 V ZB = 8
20
10 Operate
0
0.000 0.017 0.033 0.050 0.067 0.083 0.010
Time (seconds)
(b) IST = 110 A RS = 2 IS = 100 A
Fig. 16. Operate and restraint values from a differential relay connected to
IE = 10 A the CTs in Fig. 14
IP VE = 1000 V VS = 800 V ZB = 8
E. CT Analysis Tools
Equation (5) shows how to verify if a CT will saturate for a
given fault current and burden. In addition to this equation,
other tools exist to analyze the behavior of a CT or sets of CTs.
Fig. 14. Two C800 CTs, one with more copper and less iron (a) and one 1) CT Saturation Theory and Calculator
with more iron and less copper (b)
Perhaps the most common is the “CT Saturation Theory and
Using a modeling tool like those described in Section III, Calculator,” an Excel® spreadsheet created by the IEEE Power
Subsection E, we can plot the response of both of these CTs and System Relaying Committee (PSRC) [15]. This spreadsheet is
compare the results. Fig. 15 shows the two CTs modeled in a shown in Fig. 17 and is well documented. Using the available
differential circuit, with the differential current being the tutorials, a level of proficiency can be gained in a short period
difference between the two CTs. We can see that the two CTs of time. The necessary inputs are a description of the CT
do indeed perform very similarly except in the second cycle, (e.g., ratio, accuracy class, and so on) and data on the available
where CT B saturates sooner than CT A as the evaluation of (5) fault current and X/R ratio. The calculator allows the settings
predicted. After the second cycle, the CTs have nearly the same and specification engineer to play “what if” scenarios with the
error and the errors mostly cancel (resulting in little differential CT. Being able to adequately predict the CT’s behavior during
current). fault conditions can help form well-reasoned decisions when
creating relay settings.
11
Fig. 17. Simulation in the “CT Saturation Theory and Calculator” spreadsheet for a C400 CT [15]
Fig. 18. Simulation in the “CT Saturation Theory and Calculator” spreadsheet for a C800 CT [15]
The simulations in Fig. 17 and Fig. 18 show the difference 2) Two CTs
between a C400 (Fig. 17) and C800 (Fig. 18) CT selection for The authors of [8] developed a program that simulates the
the same system and burden conditions. If a user runs this transient behavior of two CTs in a differential circuit. The
simulation for a C400 CT at a certain location, they can easily program was originally written in BASIC programming
see that the CT will saturate for the given fault current. At this language, but has since been converted to an executable
point, the user has some choices. Because the system fault program that can be run on a Windows® operating system
current and X/R ratio cannot be changed, they can try lowering (shown in Fig. 19). In addition to plotting the output of the CTs,
the burden or using a CT with a higher voltage rating or a higher the program produces ASCII and COMTRADE files that can
turns ratio. Fig. 18 shows that a C800 CT will still saturate but be used to play the signals into a relay for testing.
less quickly and severely.
12
100
Iunsat
60
0 1:ICW2
Isat 20
–20
–50 –60
0 2 4 6 8 –100
3 4 5 6 7 8 9 10
Cycles
Cycles
Fig. 21. CT saturation in a raw (a) and filtered (b) event report [19]
Fig. 23. C-phase CT saturates and recovers during a phase-to-phase fault
3) Which Type of Event Should I Download?
Realize that the currents the relay sees will vary depending
Both raw and filtered event types are useful and necessary
on what system condition is causing the saturation (fault,
in event analysis. The type of event report that should be used
inrush, and so on). In general, look for currents to not accurately
for analysis depends on whether the protection algorithm being
follow the waveform expected for the given system conditions.
analyzed uses raw or filtered data. The vast majority of
Take the system condition of inrush on a power transformer, for
protection algorithms use filtered data, and thus the filtered
example. A relay sees inrush currents when a downstream
event report is the most often analyzed. However, as shown in
transformer is first energized, and these inrush currents have an
Fig. 21, some power system conditions (CT saturation, inrush,
expected signature. Fig. 24 shows three examples of event
and so on) can only be detected in raw event reports because
reports during an inrush condition, each with varying levels of
the relay filters remove the characteristics that make them
CT saturation. Fig. 24a shows an inrush event with no CT
unique. Even the most seasoned event analysis expert will not
saturation—this is what inrush is supposed to look like [21].
have much luck looking for CT saturation in filtered event
Fig. 24b and Fig. 24c show a different inrush event. In Fig. 24b,
reports. Be sure to understand whether the data is raw or filtered
the B-phase CT is saturated slightly, and in Fig. 24c, there is
when analyzing a relay operation, and always download both
severe CT saturation on all three phases.
raw and filtered event reports after a fault.
In addition to whether the event report is raw or filtered, it is a) 20
15
important to be aware of the sampling rate of the data in the 10 2:IA69
event report. Fig. 22 shows CT saturation error between two 5 2:IB69
0 2:IC69
CTs that only lasts for a quarter cycle but manages to cause a –5
misoperation [20]. This could easily have been missed if the –10
–15
user was looking at a 4 sample-per-cycle event report. It is b) 600
important to always download the highest resolution data 400
1:IA_A
200
possible after a fault. 0 1:IB_A
–200 1:IC_A
1:IBW –400
8,000
1:IBX_NEG –600
4,000 –800
c) 600
0
400
200 1:IA
–4,000 1:IB
0
–200 1:IC
–8,000
25.030581 25.040581 25.050581 –400
–600
1,400 1:IB_SUM –800
1,000 Time
600
Fig. 24. Inrush currents with no CT saturation (a), slight CT saturation (b),
200 and severe CT saturation (c)
–200
25.030581 25.040581 25.050581
Time (s)
Fig. 22. Difference error from CTs lasting for a quarter cycle
14
2) Look for Frequency Changes compared with when a nonsaturated waveform gets passed
Depending on the level of CT saturation and how long it through. For example, the angle of the filtered current phasor
persists, it may not always be as obvious as what we have seen will start to go in the leading direction. It is important to note
so far. Light saturation, while capable of still causing relay that distortion in a filtered waveform may not be the result of
misoperations, does not always give us the typical sawtooth CT saturation—that is just one possible cause. Step changes in
waveforms shown previously. Fig. 25 shows a current magnitude or current reversals from sequential breaker opening
waveform from a saturated CT that might not be immediately (that result in a step change in angle), and evolving faults can
obvious. In cases like this, we can use the fact that saturated cause similar results. However, this distortion is a good hint to
CTs output currents rich in harmonics to our advantage. One download and analyze the raw event report and check for CT
way to detect higher-level harmonics is to look for frequencies saturation.
above 60 Hz.
3) Look for Harmonic Content
We can detect a change in frequency by using event analysis
Another way to detect if a CT has saturated is to look at the
software to measure and compare the periods of different parts
individual harmonics. Asymmetrical saturation produces a
of a signal in a raw event report. In Fig. 25, we placed two
large amount of even harmonics and symmetrical saturation
vertical lines to create a window that starts at one zero crossing
produces mainly odd harmonics [22]. Event analysis software
and ends at another zero crossing approximately 1 cycle later.
can be used to show the amount of harmonic content in the two
The software shows that the time period between these two
one-cycle windows of the previous event shown in Fig. 25 and
lines is 16.147 ms, which is very close to 16.666 ms (1 cycle,
Fig. 26. Fig. 27a shows the second-harmonic content during the
based on a 60 Hz nominal frequency). This means that the
first one-cycle window where the CT is not saturating (Fig. 25).
current signal between these two points in time is almost a true
Fig. 27b shows the second-harmonic content during the second
60 Hz signal and the CT is not experiencing much saturation.
one-cycle window where the CT is saturating (Fig. 26). Notice
20,000 that there is significantly more second-harmonic content in the
10,000 waveform when the CT has saturated.
0 a) b)
–10,000 100
1:IB
–20,000 80 37.58%
1:IB
–30,000 60 134.4199 (rms)
4.61% 120 Hz
–40,000 40 4.1226 (rms)
16.147 ms Second Harmonic
120 Hz
12:26:32.97 12:26:32.99 20 Second Harmonic
Time 0
0 120 240 360 0 120 240 360
Fig. 25. No frequency distortion exists when the period of a signal is Hz Hz
16.666 ms
Fig. 27. Second-harmonic content of current waveform when CT is not
Now, we move the time cursors further along the waveform saturated (a) and is saturated (b)
and compare the results. In Fig. 26, the period between the two
Just like with frequency changes, CT saturation is not the
zero crossings has decreased to 14.616 ms. This shows that the
only cause of harmonic content in current waveforms. The
signal has a frequency higher than 60 Hz, and it could be a sign
inrush condition shown in the previous section also generates
that the CT is saturating.
large amounts of harmonic content that the relay will measure.
20,000 Measuring harmonic content is just one method that can be used
10,000
when trying to determine if a CT has saturated.
0 4) Look for a Sudden Loss of DC Offset
–10,000 In cases where a fault current contains dc offset, the sudden
–20,000 dissipation of this offset can be a sign that CTs have saturated.
–30,000
Theory and testing that shows why dc disappears from the
–40,000
secondary current when saturation occurs is given in [4].
14.616 ms Fig. 28 gives an illustration of what a sudden loss in dc offset
12:26:32.99 12:26:33.00 looks like. Normally, we would expect a current with dc offset
Time to slowly decay along an exponential curve (shown in red in
Fig. 26. Frequency is greater than 60 Hz when the period of a signal is less Fig. 28). We can see that the waveform follows the exponential
than 16.666 ms curve for the first 2–3 cycles, at which point the dc offset
This same technique can be applied to filtered event reports completely goes away and the waveform is symmetrical around
to a certain extent. When a saturated waveform gets passed the 0 axis. If there is rapid dissipation of dc offset in a raw event
through the cosine filter, there is some distortion in the output report, suspect CT saturation.
15
current (calculated by the relay as the sum of all three phase Time (ms)
CTs or wired into a relay input in a residual connection) is a Fig. 30. External B-phase-to-C-phase fault on a two-terminal differential
sign of CT saturation—unexpected neutral current through a
Based on the principle of differential protection, it is
core-balanced CT is not. Whenever CT saturation is suspected
expected that for this external fault, IBW1 will be equal to and
from filtered event reports, the raw data should be downloaded
180 degrees out of phase with IBW2. The same is true for ICW1
to verify.
and ICW2. Fig. 31 shows the B-phase and C-phase currents
60 compared across the differential. The Winding 2 data points are
20
IAW1
IBW1
negated to make the difference between waveforms easier to
5ICW1
–20 see. If there had been no CT saturation, the curves would be on
10:01:17.615000 AM top of each other. However, we can see that the Winding 1
–60
14
currents experienced saturation around the 670 ms point in the
10 IRW1.Mag fault data, causing distortion and a difference from Winding 2.
6 The magnitudes of the differential signals resulting from the
2 saturation on Winding 1 can be found by adding the Winding 1
10:01:17.615000 AM
–2
and Winding 2 data points, as shown by the black traces
80
40 IAW2
in Fig. 31.
0 IBW2
5ICW2 Whenever phase currents across a differential are not equal
–40 and 180 degrees out of phase during an external fault condition,
10:01:17.615000 AM
–80 CT saturation is suspected. Remembering to account for CT
12 ratio differences and transformer and CT connections is
8 IRW2.Mag
important, as these things may also result in the phase currents
4
not being equal and 180 degrees out of phase during normal
0 10:01:17.615000 AM
conditions.
–2 –1 0 1 2 3 4 5 6
Cycles
Fig. 29. False residual current appears when CTs saturate during an external
fault [21]
16
In this example, because both sets of CTs see the same current
40 IBW1 with a 180-degree polarity difference, we would expect the
NEG_IBW2
20 5
PARTIAL_IBT
same inrush waveform to be mirrored on the neutral side of the
0 generator and the differential to be balanced for this external
–20 inrush event. The fact that this is not the case caused the
–40 differential element to measure false operate current and trip.
–60
TREF = 9:53:45 PM
60 ICW1
NEG_ICW2
40 5
PARTIAL_ICT
20 Generator
Switchgear
0
–20
–40 TREF = 9:53:45 PM
640 650 660 670 680 690 700 A1
Time (ms) Switchgear
zero for some period of time to allow the CT to partially bottom of Fig. 36 that this is not the case. We would also expect
recover. the nonfaulted phases to remain balanced, with 120-degree
Fig. 35 shows each phase current from the bus-side CTs separation—but, Fig. 36 shows that this is also not the case. We
compared with its corresponding phase on the neutral side. The can suspect CT saturation, but we will not know for sure until
bus-side CT signals have been flipped 180 degrees to account we look at the raw signals.
for the difference in CT polarity. We can see that the signals 10
start out completely equal to each other (which is expected) and 6 IAW1
then fall away as the neutral CTs saturate. C-phase is the worst 2 IBW1
5ICW1
–2
offender, and it is the phase that tripped the relay. –6
6:19:22.971000 AM
350 –10
250 IA 15
IA87_NEG
150 5 IAW2
50 IBW2
5ICW2
–5
–50 6:19:22.971000 AM
7:39:07 PM –15
–150
600 15
IB IAW1
400 5 IAW2
IB87_NEG
200
–5
0
7:39:07 PM 6:19:22.971000 AM
–200 –15
–3 –1 1 3 5
200 Cycles
0 IC
–200 IC87_NEG Fig. 36. Filtered event report for external A-phase-to-ground fault
–400
–600 7:39:07 PM
Fig. 37 shows the filtered differential event report. Here, we
170 190 210 230 250 270 see the restraint currents steadily decreasing throughout the
Time (ms) event, while operate currents increase and then decrease.
Decreasing operate current over time, while filtered phase
Fig. 35. Currents from bus-side CTs vs. neutral-side CTs
currents remain steady throughout a fault, is an indication of CT
Upon investigation, it was found that the CTs in fact were saturation [21].
sized the same on both sides of the generator—both 200/5, C20,
intentionally sized to the generator’s full-load ampere rating of 2
IOP1
150 A. With CTs rated similarly, it is unlikely that one set 1
IRT1
would perform so differently than the other. It was then found
that the length of the CT secondary leads between the CTs and 0
the relay were considerably different. The one-way lead length
of the bus-side CTs was 15 feet, while the neutral-side CTs 0.8 IOP2
IRT2
were 20 times longer—300 feet! This greatly increased the 0.4
burden of the neutral- neutral side vs.
0
bus side) and resulted in the difference in
saturation behavior. 1.2
IOP3
The ideal solution to this problem is to match the lead length 0.8 IRT3
burden on both sets of CTs to where they saturate similarly and 0.4
balance each other out. Until that could be accomplished, a 0
temporary measure was taken by forcing the relay into high –4 –2 0 2 4 6
Cycles
security mode for a short time during startup, when the
generator breaker first closed. High security mode was then Fig. 37. Filtered differential event shows decreasing operate current over
programmed to have a higher minimum operate current and a time
more secure slope setting than the standard settings. Fig. 38 shows the raw event report. Here we can see that the
B. Transformer Relay Operates During DC Offset A-phase currents on Winding 1 and Winding 2 have a
significant amount of dc offset at the beginning of the fault.
In March of 2017, a transformer differential relay
After about seven cycles, the dc offset goes away and the
misoperated on an external A-phase-to-ground fault on the low
waveforms balance out. Note that this dc offset cannot be seen
side of a 138/69 kV wye-wye autotransformer. The filtered
in the filtered event report because it is removed by the cosine
event report from this operation is shown in Fig. 36. The relay
filter. The dc offset itself does not cause the misoperation, but
and CTs are configured such that the taps on both sides are
it does cause the CTs to saturate. The bottom graph in Fig. 38
equal, and thus we would expect the A-phase currents on
shows the A-phase current on Winding 1 compared with the
Winding 1 and Winding 2 to be equal and 180 degrees out of
inverse of the A-phase current on Winding 2, and it proves that
phase with each other. We can tell from the waveforms on the
these signals are not equal and 180 degrees out of phase because
18
the CTs have saturated. Notice that the period of the event that IOP I1 I2
has the most difference between the current signals in Fig. 38 (8)
corresponds to the period in the event with the highest operate IRT I1 I2
current in Fig. 37. It is this difference in CT behavior that
causes false operate current and results in the relay tripping on Because the CTs around the zone are assumed to have
its differential element. A very similar event is shown in [21]— opposing polarities, the operate quantity (IOP) calculates as the
an external phase-to-ground fault on the wye side of a delta- difference between the current into and out of the zone. Under
wye transformer with significant dc offset. normal load or through-fault conditions, IOP should equal 0, and
15
under internal fault conditions, IOP calculates to a value.
5 IAW1 However, if one of the CTs enclosing the zone saturates during
IBW1
–5
5ICW1 an external fault, the currents from the CTs no longer cancel
–15
each other out exactly, and IOP calculates to some value. We do
–25
6:19:22.971000 AM not want to trip for this condition, so we must add a margin to
25
provide security for CT errors. This is normally done by
15 IAW2 comparing IOP with IRT times a slope (see (9))—instead of
IBW2
5
5ICW2 comparing it with a fixed value—which can be set to an
–5
appropriate margin to account for CT saturation errors. This
–15
6:19:22.971000 AM works well because when currents are high (high values of IRT)
15
and the chance for saturation is greater, IOP needs to overcome
IAW1 a larger threshold before the relay trips. Similar to the slope in
5 IAW2_NEG percentage-restrained differential relays, the angle of the
–5 characteristic of an alpha-plane differential relay also helps
–15 account for CT saturation differences between the two sets of
6:19:22.971000 AM CTs.
–25
–5 –3 –1 1 3 5
Cycles OP RT (9)
Fig. 38. Raw event report shows significant dc offset
B. External Fault Detector Enables High Security Mode
Several solutions could be implemented to prevent this type When an external fault occurs and causes the CTs to
of misoperation. In this case, the utility chose to temporarily saturate, there is a small window of time when the fault first
raise the minimum operate current threshold and slope setting occurs that the CTs have not yet saturated. During this time
in the differential relay. At a later date, they plan to increase the period, IRT increases because of the elevated fault current, but
taps of the CTs to allow for a higher voltage before saturation IOP does not. We can use this to our advantage to develop logic
(from 300/5 to 800/5 on the 138 kV side, and from 600/5 to that detects an external fault before the CTs have had a chance
1200/5 on the 69 kV side) and return the slope to its previous to saturate. When the relay detects a large change in IRT with no
setting. The external fault detection algorithm described in the (or very small) change in IOP, it declares an external fault. This
next section would also have prevented this misoperation is shown in Fig. 39. The success of this logic depends on CTs
because restraint current increased in all three phases providing valid output for some amount of time (2 ms to a half
significantly before the operate current increased. This would cycle) [14] [20].
have put the relay into a more secure mode with a higher slope
DIO PR
setting and kept it there until the dc offset decayed. +
DIO P
IOPR
DIO PP –
VIII. APPENDIX A Therefore, the actual secondary current coming out of the
To meet the IEEE C57.13-1993 standard, a CT must not CT must be greater than or equal to 90 percent of what we
exceed 10 percent ratio correction (PRC) [26]. The standard expect it to be. This means that only 10 percent of the expected
defines the following: signal can be lost as error.
The ratio correction factor (RCF) is the ratio of the IEEE C57.13-2016 10% ratio
true ratio to the marked ratio. The primary current is error to meet the standard [9]. Ratio error is defined similarly
equal to the secondary current multiplied by the to PRC in IEEE C57.13-1993, but with an added sine term to
marked ratio times the RCF. account for the fact that IE and IS are not exactly in phase for
The percent ratio correction (PRC) is the difference cases when the connected burden has a resistive component.
between the RCF and unity, expressed in percent. When the connected burden is not purely inductive, IE and IS
We can write the definition of RCF as shown in (10), where cannot simply be added algebraically, and this extra term is
CTR is the marked ratio of the CT. needed for accuracy.
IEEE C57.13-2016 also defines the composite error the
I P / IS same way IEEE C57.13-1993 defines PRC, and it says that this
RCF
CTR can be used in place of ratio error for ring-type cores of low
(10)
IST reactance. For more on CT measurement errors, see [27].
RCF
IS For simplicity in this paper, we define ratio error the same
as PRC, as shown in (16).
We can write the definition of PRC as shown in (11).
IE
(11) (16)
IS
Substituting (10) into (11), we get (12).
IST IX. APPENDIX B
IS Recall from Section II that the ANSI voltage rating of a CT
defines the minimum secondary voltage (VS) that the CT must
IST IS
deliver to a standard burden at 20 times rated secondary current
IS IS without exceeding a 10 percent ratio error [2]. If the fault
(12)
IST IS current through a CT is higher than 20 times the rated current,
or the connected burden is higher than the standard burden, we
IS
may risk going into saturation (over 10 percent error). We can
IE use the definition for the ANSI standard voltage rating to
IS develop an equation for the level of fault current (or connected
burden) that, if we go above, will result in more than 10 percent
IEEE C57.13-1993 states that (12) must not exceed
error. We can call this the saturation point for symmetrical
10 percent [26]. We want to know what this limit means in
current (symmetrical saturation).
terms of how the actual secondary current corresponds to the
We can write the definition of the ANSI standard voltage
expected, or ideal, secondary current. A more intuitive form of
rating as follows:
I I
calculating CT error would be S_actual , or S . We can write (17)
IS_expected IST STD S RATED B STD
Using the equivalent circuit in Fig. 41 and the definition of Because we are concerned with primary fault current values,
the ANSI voltage rating, we know that the ratio error will not we can incorporate the CT ratio and convert the first term in
exceed 10 percent as long as the secondary terminal voltage VS (27) to primary values. This results in (28).
is less than the secondary terminal voltage rating VSTD, as IFAULT ZB R S
shown in (18). (28)
I PRI ZB STD R S
VS VSTD (18)
where:
The voltages VS and VSTD are driven by the core excitation IFAULT is the maximum fault current in primary amperes
voltage, VE. Therefore, we can also write the relationship in for a given fault.
(18) as (19), which considers the voltage drop across the IPRI is the primary current rating of the CT (e.g., for a
secondary winding resistance RS. 2000/5 CT, IPRI is 2000.)
VE VE STD (19) ZB is the actual burden of the CT’s secondary circuit,
which includes both the impedance of the connected relay
We can write the definition of VE STD as: and the impedance of the leads from the CT to the relay.
E STD S RATED B STD S (20) RS is the internal resistance of the CT secondary winding.
ZB STD is the standard burden of the CT (e.g., for a C800
We can write the actual voltage VE as:
CT, the ZB STD
E S B S
(21) Equation (28) can be used to evaluate CT performance for
Therefore, based on the relationship in (19): symmetrical fault currents but should never be used in practice.
This derivation of (28) is shown simply as a basis for (5).
(22)
S B S S RATED B STD S
Equation (28) is not adequate for analyzing how a CT will
Rearranging terms, we obtain: behave in the presence of (typically asymmetrical) fault current,
IS ZB RS and (5) should always be used in practice. For even better
(23) accuracy, use a transient modeling program to model CT
IS RATED ZB STD RS
behavior in the time domain.
To make this equation useful, we must write it in terms of Note that the basis for (28) was first introduced in [1]. Later,
the secondary fault current, which is IST in Fig. 41. IST is the in [5], the original equation was adjusted to consider the CT
sum of IS and IE. For the worst case right at the saturation point, secondary resistance RS, resulting in the more accurate analysis
[9] defines the error as 10 percent (IE/IS Using this reflected here.
relationship, we can define:
(24) X. APPENDIX C
S ST
Equation (5), repeated in (29), is used to analyze CT
We can then write (23) as:
saturation resulting from fault current with dc offset.
ST B S
(25) A C800, 2000/5 CT with RS
IS RATED ZB STD RS burden (including relay and lead resistance). The system X/R
ratio is 12. What is the maximum primary three-phase fault
If we define IST as the secondary fault current If sec, (25)
current that can be applied to this CT without exceeding the
becomes:
maximum flux density (which is proportional to the volt-time
If sec ZB RS area)?
(26)
IS RATED ZB STD RS To determine the maximum fault current, set (5) to the
maximum of 20 and solve for IFAULT, as shown in (29) and
We can use (26) to solve for the maximum fault current (30).
(If sec) that the CT can handle before it goes into saturation. We
can also use it to solve for the maximum burden for a given IFAULT ZB R S X
(29)
fault current. IPRI ZB STD R S R
In development of (26), we assumed an inductive burden. I FAULT 1 0.5
Therefore, IS and IE are in phase and can be added algebraically (30)
2000 8 0.5
to get IST. This assumption also allows us to use the simplified
error equation for ratio error in Appendix A. Using modern I FAULT 17.4 kA
resistive burdens with these simplified equations produces
errors. To account for these errors, a more conservative form of XI. REFERENCES
(19) is used: [1] S. E. Zocholl, Analyzing and Applying Current Transformers.
Schweitzer Engineering Laboratories, Inc., Pullman, WA, 2004.
If sec ZB R S
(27) [2] IEEE Standard C37.110-2007, IEEE Guide for the Application of
IS RATED ZB STD R S Current Transformers Used for Protective Relaying Purposes.