Filtering NumericalProtectionRelays
Filtering NumericalProtectionRelays
Filtering NumericalProtectionRelays
Ilhan Kocar
21/03/2022
Evolution of Technology
3
Basic Components of a Numerical Relay - II
4
Analog Input Acquisition
The analog input is the waveforms of voltages and currents
Input from power system is through Current Transformers (CT) and Voltage Transformers
(VT)
Protective relays require reasonably accurate reproduction of the normal, tolerable, and
intolerable conditions in the power system for correct sensing and operation1.
Typical Sampling
▪ Three voltage signals: phase voltages
▪ Three current signals: line currents
5
Fundamental Notions on Signal Processing - I
Periodic signal
▪ Same magnitude after a time interval of T called period. The voltage and current signals in North America have a
period of 16.66 ms. The period is often simply called cycle.
Frequency of a periodic signal
▪ 1/T, 60 Hz or cycles per second in North America
Sampling frequency
▪ Number of samples over a time interval of 1 second. In modern numerical relays: 240, 480, 720, 960 or 1920 Hz.
Period of Sampling
▪ Inverse of sampling frequency
Number of samples per cycle
▪ Sampling frequency over network frequency
Analysis window
▪ The time interval during which the signal is analyzed to extract a frequency component. In numerical relays, it is often the
fundamental frequency (60 Hz) but in certain cases harmonics (120, 180, 240 Hz etc.)
▪ Analysis window is independent form sampling frequency (if sampling is 960 Hz and window one cycle then 16 samples per window,
window half cycle than 8 samples per window)
Sampling theory: A signal with a frequency band limited to Fc can be numerically reconstructed with 2xFc sampling
6
Fundamental Notions on Signal Processing - II
For example, sampling at 960 Hz or 16 samples per second
▪ we assume the signal is limited to a frequency band of 480 Hz.
▪ All frequency components below 480 Hz can be reconstructed.
▪ All the components above 480 Hz need to be eliminated using an ideal low pass filter (anti-aliasing filter).
Sampling at 240 Hz
▪ The signal islimited to 120 Hz
▪ Above 120 Hz need to be eliminated. Not possible e.g., to reconstruct third harmonic at 180 Hz.
7
Filtering System with a Window of Two Samples - I
Consider a sinusoidal current or voltage waveform
x (t ) A sin( t )
It can also be expressed with
x (t ) A cos( )sin( t ) A sin( )cos( t ) or in matrix form
x (t1 ) sin( t1 ) cos( t1 ) A cos( )
x (t1 t) sin( t1 t ) cos( t1 t ) A sin( )
If let us say samples are taken at each 2pi/12 radians (12 samples per cycle, sampling frequency of 720 Hz at 60 Hz
8
Filtering System with a Window of Two Samples - II
y1 A cos
y2 A sin
Provided by non-recursive filters expressed as
y1 1.7321 x (0) 2.0 x ( t )
y2 1.0 x (0)
At each period of sampling, the amplitude of A can be simply obtained by the square root of the sum of squares
X y1 j y2
9
Filtering System using a Window with Many Samples
Assume we have more samples, say 12 per cycle with a step of dt
x (0) 0 1
x( t) sin( t ) cos( t) A cos( ) Expressed in
X = BY 0 1
A sin( ) matrix form
x (11 t ) sin(11 t ) cos(11 t)
0.5 0.866
0.866 0.5
More equations than unknows, rearranging to solve Y
1 0
X = BY 0.866 0.5
Results in
BT X = BT BY T
1
T
0.5 0.866
Y B B B X B
T
1
T T
1
T 0 1
B B B X= B B B BY
0.5 0.866
For a given network frequency and fixed rate of sampling, B is known a priori. 0.866 0.5
Consider 60 Hz with 720 Hz rate of sampling (12 samplings per cycle), B is given on the right.
Define also matrix C to simplify the final system of equation 1 0
T
1
T
0.866 0.5
C B B B Given B on the right, what is the dimension of C?
0.5 0.866
10
Sinusoidal Temporal and Frequency Response Analysis
Time domain response of filters depend on the sampling rate and number of samples per window (or window size),
When we deal with sinusoidal functions, the filters are Fourier filters with variable window size
Consider below 12 samples per cycle and varying window size
6 samples
window size
2 samples
3 samples
12
samples
11
Sinusoidal Temporal and Frequency Response Analysis
Time domain response of filters depend on the sampling rate and number of samples per window (or window size),
When we deal with sinusoidal functions, the filters are Fourier filters with variable window size
Consider below 12 samples per cycle and varying window size
window size
24 samples
120 samples
12
Fourier Filter 1/2
In case of a window size corresponding to one network cycle, we end up with full-cycle Fourier filter.
See the matrix C
k2
cos
T
1 2 N
C B B BT k 0,1, 2, N This is to extract the fundamental component (60 Hz)
N k2
sin
N
N 1
2 3n 2
x (n T ) cos
A cos( ) N n 0 N
Y N 1
A sin( ) 2 3n 2
x (n T ) sin
N n 0 N
This is to extract the phasor of the
fifth harmonic component (180 Hz)
N 1
2 5n 2
x (n T ) cos
A cos( ) N n 0 N
Y N 1
A sin( ) 2 5n 2
x (n T ) sin
N n 0 N
14
Fourier Filter Exercise 1/2
Full cycle filter with 12 samples per cycle, extraction of fundamental component (60 Hz Phasor)
% Full cycle Fourier Filter
Nsamples=12; %number of samples
t=0:1/60/Nsamples:(Nsamples-1)/60/Nsamples; %for signal sampling
x=1.0*sin(2*pi*60*t)+0.3*sin(3*2*pi*60*t); % sampled sinusoidal signal
%filter
y11=0; y21=0; %cosine and sine filter,
HC=1; %fundamental component
for k=1:Nsamples
kk=k-1;
y11=y11+x(k)*cos(HC*2*pi*kk/Nsamples);
y21=y21+x(k)*sin(HC*2*pi*kk/Nsamples);
end
y11=y11*2/Nsamples;
y21=y21*2/Nsamples;
y=sqrt(y11^2+y21^2); %amplitude is 1.0
15
Fourier Filter Exercise 2/2
Full cycle filter with 12 samples per cycle, extraction of the third harmonic
% Full cycle Fourier Filter
Nsamples=12; %number of samples
t=0:1/60/Nsamples:(Nsamples-1)/60/Nsamples; %for signal sampling
x=sin(2*pi*60*t)+0.3*sin(3*2*pi*60*t); % sampled sinusoidal signal
%filter
y11=0; y21=0; %cosine and sine filter,
HC=3; %third frequency component
for k=1:Nsamples
kk=k-1;
y11=y11+x(k)*cos(HC*2*pi*kk/Nsamples);
y21=y21+x(k)*sin(HC*2*pi*kk/Nsamples);
end
y11=y11*2/Nsamples;
y21=y21*2/Nsamples;
y=sqrt(y11^2+y21^2); %amplitude is 0.3
16
Fourier Filter – Typical Relay Application
Normalized frequency response of cosine filter Normalized frequency response of sine filter
Combined response
17
Fourier Filter – Typical Relay Implementation
18
DC Component Problem
Fault current has an exponentially decaying dc component
19
DC Component Problem
Fault current has an exponentially decaying dc component
Circuit representing
fault inception
Solutions: Mimic impedance, least squares filter, cosine filter
20
DC Component Problem, Cosine Filter
Fourier filter consists of two filters: cosine and sine
Cosine filter has a better performance in rejecting dc component
21
Numerical Protection Relays
Claudio Rimada
Senior Subject Matter Expert
April 2021
Program
1. Relay Types 6. Relay digital input/output
- Relay Trip Signals.
2. Protective relays - Input / Output.
-Time Synchronization.
3. Protection zones - Communications.
- Serial / Ethernet comms.
4. Relay analog inputs
- Current Measurements. 7. MiCOM Px40 Platform
- Voltage measurements. (Hardware)
5. Characteristics of numerical
relays
- Introduction.
- Relay Architecture.
A Alarm
ac or AC Alternating current
B Bus, battery, blower
BP Bypass
BT Bus tie
C Current, close, control, capacitor, compensator, case
CC Closing coil, coupling capacitor, carrier current
CS Control switch, contactor switch
CT Current transformer
CCVT Coupling capacitor voltage device
D D Down, direct, discharge
dc or DC Direct current
E Exciter, excitation
F Field, feeder, fan
G* Ground, generator
GND Ground
70 - Rheostat
72 - DC Circuit Breaker
73 - Load-Resistor Contactor
74 - Alarm Relay
●21
●87T
●50/51
● Input from the power system is usually through current transformer (CT)
and voltage transformer (VT).
● Bar Type Current Transformer. One that has a fixed, insulated, straight
conductor in the form of a bar, rod, or tube that is a single primary turn passing
through the magnetic circuit and is assembled to the secondary, core and
winding.
● Window Type Current Transformer. One that has a secondary winding insulated
from and permanently assembled on the core, but has no primary winding as an
integral part of the structure. Complete or partial insulation is provided for a
primary winding in the window through which one or more turns of the line
conductor can be threaded to provide the primary winding.
● Bushing Type Current Transformer. One that has an annular core and a
secondary winding insulated from and permanently assembled on the core, but
has no primary winding or insulation for a primary winding. This type of current
transformer is for use with a fully insulated conductor as the primary winding.
3000 spires
1A
245 kV
1A
400 kV
63 kV
36 kV
20 kV
10 kV
P1 P2
S1 S2
P1 P2
Ip
S1 S2
P1 ➔ P2
S1 ➔ S2 Externally
P1 P2
Ip
Test A
If Ammeter deflects + X ➔ S1
If Ammeter deflects - X ➔ S2
245 kV
100 V
P1
S1
P2 S2
Np.Us = Ns.Up
Protection
R S T
P1 S1
P2 S2
Protection
R S T
P1 S1
P2 S2
P2 S2
● Self-diagnosis.
● Relay software is able to conduct test routines to determine the correct operation of
hardware components. In case anomalies are found WD will be activated.
● Adaptive protection.
● Due to communication facilities relay can be adaptive. This feature enables the relay
setting to be changed depending on the operating condition of the network.
● Measurements.
● All the possible measurements taken formerly by panel meters (Voltages, currents, frequency,
power, energy, cos phi, etc)
● Control.
● The ability to operate primary equipment such as breakers and disconnectors by means of Input
and outputs.
● Fault record.
● Once a protection issue a trip, its very useful to have a simple indication of: What protection
element tripped, what was the fault value. This will give the operator a clear view about the fault
nature and a clue how to amend the problem and re-establish the substation operation.
● Disturbance recording.
● The ability to record the samples belonging to: current, voltages and I/O. This serves
two purposes: Fault analysis and Fault reproduction using digital injection devices.
● Logic.
● The capacity to program Boolean logic between input/output to solve a particular
problem such as: interlocking, blocking, intertripping, etc.
● Communications.
● One or more communication channels are available for: remote setting, remote
measurements, control, etc.
Com 2
IRIG-B
Signal Conditioner
Digital Inputs
Serial Comms.
CT Signal Conditioner
Signal Conditioner
RTC
Keyboard,
MMI Display,
Analog Inputs
LEDs
VT Signal Conditioner Multiplexor
Processor
Digital Outputs
Signal Conditioner
Sample and Hold A/D converter
Signal Conditioner
RAM Memory
WD system WD
Ethernet /
Flash memory Comms. processor
Fx
Parallel port
Frecuency
Frecuency Frecuency
-Fmax Fmax -fs -FL FL fs
Frecuency
Adaptavive Gain
CT Anti-aliasing filter
Amplifier
Amplitude
Multiplexor
data
Sample and Hold A/D converter
Time
EOC
SOC
S/H
Mux-Control
4
Amplitude
Time
Sample
Conversion time
Hold time
𝑵−𝟏
𝟐𝝅 𝒌 𝒏
−𝒋
𝑿𝒌 = 𝒙𝒏 𝒆 𝑵
𝒏=𝟎
𝑵−𝟏
𝟐𝝅 𝒌 𝒏
−𝒋
𝑿𝒌 = 𝒙𝒏 𝒆 𝑵 (1)
𝒏=𝟎
𝑵−𝟏
𝟐𝝅 𝒌 𝒏 𝟐𝝅 𝒌 𝒏 Bk
𝑿𝒌 = 𝒙𝒏 [𝒄𝒐𝒔 − + 𝒋 𝒔𝒊𝒏 − ] (2)
𝑵 𝑵
𝒏=𝟎
𝑩
𝒋 𝒕𝒂𝒏−𝟏 𝑨𝒌
𝑿𝒌 = 𝑨𝒌 + 𝒋 𝑩𝒌 = 𝑨𝒌 𝟐 + 𝑩𝒌 𝟐 𝒆 𝒌 (3)
Ak
Real
● If we take 8 samples beginning at t=0, the sampled values (xn) will be:
x0=0; x1=0.707; x2=1; x3=0.707; x4=0; x5=-0.707; x6=-1; x7=-0.707
● (X4; X5; X6; X7 ) should be discarded because they are over the Nyquist limit. (5)
Because of this the (Xk) will result in:
X0=0; X1=-8j; X2=0; X3=0; X4=0
● As we took 8 samples to make the calculation we average the value 8/8=1 of (6)
the coefficients.
● X0=0; X1=-1j; X2=0; X3=0; X4=0
4 8
Magnitude
(6)
1
Frequency [Hz]
X0 X1 X2 X3 X4
(0Hz) (1Hz) (2Hz) (3Hz) (4Hz)
● The 23rd is the first predominant harmonic that is not attenuated by the
Fourier filter and this is known as an ‘Alias’. However, the Alias is
attenuated by approximately 85% by an additional, analogue, ‘anti-
aliasing’ filter (low pass filter).
● Fourier components are stored in memory so they can be accessed by all the
protection elements algorithm.
● Unprocessed samples are used for: wave form recorder, RMS calculation,
metering etc.
AC protected circuit
52
CT
Trip
I
Protective
Relay
V
VT
Fault
Protection Relay
I Func. Delay DO Delay
ANSI
Function
V
Prot. trip
PS
CB Delay
CB
-P +P 52-A
CB Status
TC
-B
PS PS
CB
-P1 +P1 52-A +P2 -P2
52-A
TC1 TC2
-B1 -B2
Protection Relay
I
ANSI
Function
V
PS
CB
52-A
86
TC
-B
-P
Schneider Electric, Claudio Rimada, April 2021 59
DC trip, Open/Close circuits with 86
function
+VB (CB voltage) +VP (Prot. voltage)
Control-Close
Local-Close
Local-Open
Control-Open
52a 86
● Outputs
● Electromechanical relays (operation time 5-7ms)
● Solid State relays (SSD) consisting of:
1. Combination of IGBT with a varistor.
2. Electromechanical relay.
3. The operation time is 0.2ms. This response is suitable for sub-cycle
relays (Ex: Line protections).
IRIG-B IRIG-B
Bullet Antenna
IRIG-B
Mod
IRIG-B
Unmod
Switch (Fx or Eth)
SNTP SNTP
Client Client
Main Bck-Up
SEL Px40
ANT
Amplifier P594
Item. Description
1 Satellite Antena Power [W] = 26 W
2 Satellite Antena Input Power [dbm] = 44.14973 dbm
3 Satellite Antena Gain [db] 13 db
4 Satellite Antena Output Power [dbm] = 57.14973 dbm
● RS-232 / RS-485 ports are available for legacy protocols such as:
●MODBUS
●DNP3.0
●IEC60870
●Courier
● Ethernet Communications.
●100Mbits/s Fiber optic + 10/100Mbits/s Copper.
●A special processor is used to implement various protocols:
● IEC61850
● DNP3.0 over TCP/IP
PC running
MiCOM Studio
Courier
`
RS232 to RS485
converter
Settings, PSL,
DR, Events
Laptop running
MiCOM Studio
Settings, PSL,
MCL, DR, Events
RS232 to RS485
converter
Courier
Intranet
Substation Substation Serial Device
Router Switch Server
PC running
MiCOM Studio
PLC (uploading
measurements and status)
RS232 to RS485
converter
PACiS OI
PACiS GATEWAY
TX-RJ45 Copper
Ethernet
RS1600T (14TX) or RS8000T (6TX)
FX-Optical Fiber
Ethernet
●Gateway
●Ethernet standard:
●Eth Switch
●BCU
●IED ●IED
●IED ●IED
●IED
●Gateway
●Eth ●Eth
●Ethernet non-standard:
- Basée en RSTP ou e-RSTP
- Mécanisme d’auto-cicatrisation
Configurable
Interface
Dual Rated Programmable
CT Inputs LEDs
Flexible
Communications
Comprehensive
Recording
Programmable
Scheme Logic
MiCOM40-78
Menu Navigation
Alarm viewer Keys
(Read Key)
Programmable
Battery LEDs
back-up
25 Pin Download/
9 Pin Local
Monitor port
communications
Security seal