Removal Of Decaying Offset; Phaselet Computation - GE UR Series L90 Instruction Manual

9.1 OVERVIEW
other relay terminals. Also, the sum of the squares of the raw data samples is computed for each phase, and transmitted
with the phaselets.
At the receiving relay, the received phaselets are combined into phasors. Also, ground current is reconstructed from phase
information. An elliptical restraint region is computed by combining sources of measurement error. In addition to the
restraint region, a separate disturbance detector is used to enhance security.
The possibility of a fault is indicated by the detection of a disturbance as well as the sum of the current phasors falling out-
side of the elliptical restraint region. The statistical distance from the phasor to the restraint region is an indication of the
severity of the fault. To provide speed of response that is commensurate with fault severity, the distance is filtered. For mild
faults, filtering improves measurement precision at the expense of a slight delay, on the order of one cycle. Severe faults
are detected within a single phaselet. Whenever the sum of phasors falls within the elliptical restraint region, the system
assumes there is no fault, and uses whatever information is available for fine adjustment of the clocks.
The inductive behavior of power system transmission lines gives rise to decaying exponential offsets during transient con-
ditions, which could lead to errors and interfere with the determination of how well measured current fits a sinewave.
The current signals are pre-filtered using an improved digital MIMIC filter. The filter removes effectively the DC compo-
nent(s) guaranteeing transient overshoot below 2% regardless of the initial magnitude and time constant of the dc compo-
nent(s). The filter has significantly better filtering properties for higher frequencies as compared with a classical MIMIC filter.
This was possible without introducing any significant phase delay thanks to the high sampling rate used by the relay. The
output of the MIMIC calculation is the input for the phaselet computation. The MIMIC computation is applied to the data
samples for each phase at each terminal. The equation shown is for one phase at one terminal.
Phaselets are partial sums in the computation for fitting a sine function to measured samples. Each slave computes phase-
lets for each phase current and transmits phaselet information to the master for conversion into phasors. Phaselets enable
the efficient computation of phasors over sample windows that are not restricted to an integer multiple of a half cycle at the
power system frequency. Determining the fundamental power system frequency component of current data samples by
minimizing the sum of the squares of the errors gives rise to the first frequency component of the Discrete Fourier Trans-
form (DFT). In the case of a data window that is a multiple of a half cycle, the computation is simply sine and cosine
weighted sums of the data samples. In the case of a window that is not a multiple of a half-cycle, there is an additional cor-
rection that results from the sine and cosine functions not being orthogonal over such a window. However, the computation
can be expressed as a two by two matrix multiplication of the sine and cosine weighted sums.
Phaselets and sum of squares are computed for each phase at each terminal as follows. For the real part, we have:
For the imaginary part, we have:
where: k is the present phaselet index,
N is the number of samples per cycle, and
9
p is the present sample index
The computation of phaselets and sum of squares is basically a consolidation process. The phaselet sums are converted
into stationary phasors by multiplying by a precomputed matrix. Phaselets and partial sums of squares are computed and
time stamped at each relay and communicated to the remote relay terminals, where they are added and the matrix multipli-
cation is performed. Since the sampling clocks are synchronized, the time stamp is simply a sequence number.
9-2

N 2 1
–

4
i
--- -
I

= 1_f_A k p
–
1_Re_A k  
N
p 0
=

N 2 1
–

4
i
I
--- -

= – 1_f_A k p
–
1_Im_A k  
N
p 0
=
L90 Line Current Differential System
9 THEORY OF OPERATION

9.1.3 REMOVAL OF DECAYING OFFSET

9.1.4 PHASELET COMPUTATION

 1 2  
2 p
 
+
 ------------------------------- -
cos

 
N
 1 2  
2 p
 
+

sin ------------------------------- -

 
N
(EQ 9.1)
(EQ 9.2)
GE Multilin