Complete Elimination of DC Offset in Current Signals for Relaying Applications
Jun-Zhe Yang
Chih-Wen Liu
Department of Electrical Engineering, National Taiwan University, Taipei, Taiwan
Abstract
-
This paper proposes a new method for exact estimation of DC offset in current signals such that DC offset can be completely eliminated for relaying applications. This new method, we called Smart Discrete Fourier Transform (SDFT), based on Discrete Fourier Transform (DFT) keeps all advantages of DFT, and smartly avoids the defects resulting from that frequency scales can’t match signal characteristics. Moreover, we use smoothing windows to filter noise in practice. The difference between SDFT and other methods, which use smoothing windows, is that SDFT can alleviate the difficulties of phase shift and amplitude decay caused by smoothing windows. Recursive and half-cycle computing can also be used in SDFT, soSDFT
is very suitable for use in real-time. To validate the claimed benefits of SDFT, the performances of SDFT algorithm are demonstrated with the data generated by EM” simulator. The simulation results are very encouraging.Keywords: Discrete Fourier Transforms (DFT), Frequency
estimation, phasor measurement.
I. Introduction
In power systems, high voltage transmission and distribution lines are vital links fiom the generating plants to the end users. From views of security, economics, and quality of power feeding, it is necessary for protecting lines using relay [ 11. When fault occurs, relays prevent involving area from growing. Hence correct action of relay is important in power systems. However, it is a fact that a line relay tends to overreach in the presence of DC offset component
in
fault current waveforms. Therefore, DC offset component has to be removed from these waveforms SDFT can estimate DC offset component exactly. In fact, SDFT obtains exact DC offset from fault current waveforms such that it can be eliminated. This is useful for fault detection and location. Three advantages of using SDFT are presented here. First, relays c q work correctly. Without influence of DC offset, relays ‘obtain accurate frequency and phasor components from SDFT. Second,L/R
time constant of DC offset can be used in line parameter estimation after fault occurs. Uncertainty about the line parameters, which usually results in a significant error in the calculated fault location, can also be resolved by SDFT. Third, SDFT obtains exact fault angle of DC offset. Therefore, one can use this extra information in protection and control.0-7803-5935-6/00/$10.00 (c) 2000 IEEE
Noise is unavoidable in power systems. Hence every algorithm applied to power systems has to take noise into consideration. Currently more prevalent methods are using smoothing windows. Generally, using smoothing windows is known to provide flexible solutions to filtering noise, but loses the ability of computing phasors. SDFT can allay the phase shift and amplitude decay programs caused by smoothing windows.
This paper is organized into four sections, the first of which is introduction. In section 11, we introduce the proposed algorithm, SDFT, used in estimating DC offset, and prove that SDFT can obtain phasors when using smoothing windows. We test SDFT using EMTP simulation data in Section 111. Finally, section
IV
is devoted to conclusions.11. The Proposed Digital Algorithm
This section proposes the algorithm of SDFT that estimates the DC offset, frequency, and phasor from a fault current signal. Since there are several components in a fault current signal [4-61, the algorithm first takes DC offset into consideration and uses smoothing windows to eliminate other components in a fault current signal. The main components of a fault current signal can be expressed as follows [7]:
x ( f ) = Xsin(of+&)+ Xsin(q$)exp(-ar) (1) where
X
: the amplitude of the signal,4,
: the phase angle of the signal, : the fault angle of the signal,1
-
= r : the time constant of the signal aSuppose that x ( t ) is sampled with a sampling rate (60*N) Hz to produce the sample set { x(k)}.
k k
x ( k ) = Xsin(a,-++l)+ Xsin(q$)exp(-cr-) 60 N 60 N
The signal x ( t ) is conventionally represented by a phasor (a complex number)
X
.
Z=Xe”s =Xcos~,+jXsin$, (3)
Then ~ ( t ) can be expressed as
(4) Fe,- -ye-”’
+
k X sin( +>) exp( -a -) 2 60 N x ( t ) = 1933where
*
denotes complex conjugate.Moreover, the fundamental fiquency (60Hz) component of DFT of { ~ ( k )
1
is obtained by(5) Combing Eq. (4) and Eq. (5) and taking frequency deviation (o = 2x(60
+
Af)) into consideration, we obtainWe rearrange Eq.
(6)
as the following:We use the following identity to simplify Eq. (7).
N9 I-1 sin
-
e sin - 2 I..Then Eq. (7) can be expressed as
r(2 + g,
Where 8 , = K , a n d o,l 6 0 .
60N N
Rearranging Eq. (9) fiuther, we obtain
Wedefine A,, B, and C, asfollows:
Then Eq. (1 0) can be expressed as
So far the development of the algorithm of SDFT is the same as the traditional DFT method. So the SDFT can keep all advantages of DFT such as recursive and half-cycle computing manner. But in the DFT, it doesn't take DC offset into consideration and it assumes that the frequency deviation is small enough to be ignored. It always considers ir = A r , so traditional DFT based methods incur error in estimating frequency and phasor when frequency deviates fiom nominal frequency (60 Hz) or DC
offset is present. If we want to obtain exact solution, we must take B, and C, into consideration. Then we define
and
And fiom Eq.
(lo),
we will find the following relations.Then
subtracting from Eq. (20) give La-?, =
(d
-l)(u-b) A A y, = X I + , -i,b=A,(u-b)+B,(u-'-b)1
N sin (0 ) X = & ( A , - sin (N 0,) (33)Similarly, we may obtain
(35)
U%,+'
-
u(i,
+
+
i,+,
c,
=( a
-
b)(d-
1) As same treatment with Eq. (14) and Eq. (20), we mayrearrange Eq. (22)' Eq. (23), and Eq. (24) to obtain
Furthermore, we take noise into consideration and use smoothing windows to filter noise. Consider a sampled set
{ X ( A ) } to be a filtered set { z(b) } by a smoothing window
{ sw(m)l s,.s,,.
.
.S. }, m is window size.Dividing Eq. (25) by Eq. (26), we get
Then substituting Eq. (23) and Eq. (24) into Eq. (27) and
eliminating variable 'a', we obtain as follows: Moreover, the DFT of { z(~) } is given by
Solve Eq. (28) to obtain 'b'. From the definition of 'b' in Eq. (16)' we can obtain the exact solution of the time constant.
1 60Nlogb
r=-
From the definition of Eq. (1 4), we can obtain: Then rearrange
Eq.
(27) as follows(39)
Solve Eq. (30) to obtain 'a'. From the defmition of 'a' in
Eq. (16)' we can get the exact solution of the fiequency. The relations of Eq. (17), Eq. (18), and Eq. (19) are still
kept in Eq. (39). Therefore, the same steps fiom Eq. (20) to Eq. (33) can be used in Eq. (39). Hence we can estimate time constant and frequency without modifying equations, but we have to do some change in Eq. (32) and Eq. (35)
when we estimate phasor and fault angle.
(31) 60N
f = 60
+
Af = cos-'(Re(a))- 27FFrom Eq. (29) and Eq. (3 I), it is observed that SDFT can
provide exact time constant and frequency using
5 ,
2,, ,x ~ + ~ , x,,
,
and i,, in the absence of noise. Moreover, we can estimate phasor and fault angle after getting exact time constant and f?equency by the following equations:A A
A I
where j, =
i,,,
-2,b and y.+, = z,+* -i,+lb(41)
aai,+, -a(;, +i,+,
(a
-
b)(ab-
l)(s,+
s2b + * . a + s,b’-’)c,
= 2.1159 Q 2.1159 2.1159 2.1159-.
9The phasor obtained from
Eq.
(40) and fault angle obtained fiom Eq. (41) will allay the phase shift and amplitude decay caused by smoothing windows. In the next section, we will test the SDFT method with the data generated by EMTP simulator.111. Simulation Results
The first example of
EMTP
simulation is a single-phase system and the model is shown Fig.1. It is a 100 kmtransmission line connected to generators at sending and receiving end, respectively. The transmission line model used in the simulation is the distributed parameters model.
Fig. 1
EMTP
simulation model-
(e60Hz)60.001
-76.17 I J
0 0.078
time(sec) Fig. 2(e)
A fault is simulated to occur at the peak value of voltage V I and we measured current 11 at BUS1. Computation results by the SDFT are showed in Fig. 2 and theoretical results are as follows:
(0.0003
+
O.OOO5 70) 11 (0.0003+
0.0005*
30) (0.0001*
70) 11 (0.0001*
30)+
1 a-= = 0.01065135 = 10.65135* lo-’ zl = 0.0001*
70+
j w*
(0.0003 + 0.0005*
70) 22 = 0.0001*
30+ j w*
(0.0003+
0.0005*
30) zl*V,+22*V, -v,
=-
677141-82.9962’ zl+ 22+ zl*
22The current of measurement point is I , =-=
-
’’
21 159L-
76.16 I 1’ (A) zl -3 x 10 0.653I
10.852 10.651 10.65 10.649 0.078 time(sec) Fig. 2(b) x l o 4r-
1
0.078 time(sec) Fig. 2(d) 80 60 40 20 0Fig. 2 Computation results by SDFTds
time(sec)
To distinguish different kind of SDFT [8], we add ‘d’ and ‘s’ in the suffix of SDFT means this kind of SDFT has taken DC offset and Smoothing windows into consideration. In this example, we used Blackman window and window size is 28 points and SDFTd moving window
size is 36 points. It means SDFTds gets frst computing results after fault occurred 2 cycles. From Fig.2, it is clear that SDFTds succeed in eliminating
DC
offset. Since the amplitude of DC offset is decaying with time, the computing result of time constant will become senseless.Next example, we add capacitance into model, change system frequency and change location and time of a fault. Computation results are showed in Fig. 4 and theoretical results are as follows:
Transmission Line
Measurement Point R=O.OoOl(ohm/lun,~
I I
281.7LVKV
Fig.3. EMTP simulation model (P60.5Hz)
60.52 1 1 80.47 I 1 0 0.078 &-“sec) Fig. 4(a) 1.827 I 0.078 1.825
‘
tlmo(o0a) Fig. 4(c) -72.5 L I 0 0.078 tlme(seo) Fig. 4 ( e ) (0.0003+
0.0005*
80) 11 (0.0003+
0.0005*
20) (0,0001*
80)11 (0.0001
*
20)+
1 r = = 0.00819025 = 8.19025 *IO-’I
281.7*10’1
0.039-0.063% 0.01 51+ 0.004%-0.015 I-0.004% -0.0593+0.2364 281.7* IO‘L- IO”
I , = 18259L-72.3607’ (A) Fig. 4@) 15000 Fig. 4@) 15000
P
g
5000 8 0 .._-__-.-__l.-L. 0 0.078 t l m e ( m e c ) Fig. 4(d) 0 Ums(.eo) Fig. 4(f) 0-7803-5935-6/00/$10.00 (c) 2000 IEEE IFig. 4. Computation results by SDFTds
1937
Since we added capacitance into the model, the transient current signal added high fiequency component that caused by resonance between line capacitance and line inductance [4-61. In this example, the high fiequency component of transient current signal doesn’t make serious affection to computing results of SDFTds. Actually, added window size of smoothing windows is helpful to resist the high fiequency component and noise of transient current signal, but how large of window size can be added depends on how long the line breaker open. Generally speaking, the fault will be cleared after 3 cycles, and it should be enough time for using SDFTds in fault location or other applications.
IV. Conclusion
In this paper we introduce the SDFT method and demonstrate its performance with the data generated by EMTP simulator. From the ExampIes presented above, SDFT succeed in eliminating DC offset. SDFT both keeps the advantages of DFT and also deals with the cause of DC offset errors, while taking smoothing windows into consideration. These aspects make SDFT an accurate and practical method in line relaying applications.
Reference
J. A. Jrang, Y. H. Lin, C. W Liu, J. Z Yang, T. M. Too, “An
Adaphve Fault Locator System For Transmission Lines,’’
JEEIYPESSummer Meeting, July 1999, pp. 930-936
J. C. Gu and S. L. Yu, “Removal of DC Offset in Current and Voltage Signals Using a Novel Fourier Filter Algorithm,” IEEE
Trans. on Power Delhrery, PE-OSZPRD(06-99).
Gabriel Benmouyal, “Removal of DC-OBet in Current Waveforms Using Digital Mimic Filtering,” IEEE Trans on Power Defwev,
Vol. 10, No. 2, April 1995, pp. 621-630.
Adly A. Girgis, A. A. Sallam, A. Karim El-Din, “An Adaptive Protection Scheme for Advanced Series Compensated (ASC) transmission Lines,” IEEE Trans. on Pmver L k f w e v , Vol. 13, No. 2. Apnl 1998, pp. 414-420.
Adly A. Girgrs and R Grover Brown, ”Modeling of Fault-Induced Noise Signals for Computer relaying Applications.” IEEE Trans. on Power Apparatus and Systems, Vol. PAS-102, No. 9, September Adly A Girgis and R. Grover Brown, “Application of Kalman Filtering in Computer Relaying,” IEEE Trans. on Power Apparatus
and Systems, Vol PAS-100, July 1981, pp. 3387-3397.
Norman T. Stringer, ‘The Effect of DC OBet on Current- Operated relays,” IEEE Trans. on Industry Applications, Vol. 34, No. 1. JanuarylFebruaty 1998, pp. 30-34.
J. Z Yang and C. W. Liu, “A Smart Method Makes DFT More Precise for Power System Frequency Estimation:’ IEEUPES
Winter Meetzng, Feb. 1999, pp. 909-913.
1983, pp 2834-2841.
Biographies
hih-Wen Liu was born in Taiwan in 1964. He received the B.S. degree in Electrical Engineering from National Taiwan University in 1987, Ph.D. degree in electrical
engineering from Cornell University in 1994. Since 1994, he has been with National Taiwan University, where he is associate professor of electrical engineering. His research interests include application of numerical methods to power system, motor control and