Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2177 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … Evaluation of Moving Average Window Technique as Low-pass Filter in Microprocessor-Based Protecting Relays Naser Khodabakhshi-Javinani Electrical Engineering Department Amirkabir University of Technology Tehran, Iran naser@aut.ac.ir Hossein Askarian-Abyaneh Electrical Engineering Department Amirkabir University of Technology Tehran, Iran askarian@aut.ac.ir Abstract—Over the last decades, with the increase in the use of harmonic source devices, the filtering process has received more attention than ever before. Digital relays operate according to accurate thresholds and precise setting values. In signal flow graphs of relays, the low-pass filter plays a crucial role in pre- filtering and purifying waveforms performance estimating techniques to estimate the expected impedances, currents, voltage etc. The main process is conducted in the CPU through methods such as Man and Morrison, Fourier, Walsh-based techniques, least-square methods etc. To purify waveforms polluted with low- order harmonics, it is necessary to design and embed cutting frequency in a narrow band which would be costly. In this article, a technique is presented which is able to eliminate specified harmonics, noise and DC offset, attenuate whole harmonic order and hand low-pass filtered signals to CPU. The proposed method is evaluated by eight case studies and compared with first and second order low-pass filter. Keywords-Low-pass Filter; Power System Protection; Harmonics and noisy waveforms; Moving Average Window (MAW); Central Processing Unit (CPU) I. INTRODUCTION Disturbances are a main power quality issue in power systems. Voltage and current waveform distortions such as harmonics and noises are the main disturbances. Distorted waveform might be originated by sources, appliances, and loads [1-4]. Most of the harmonic sources such as power electronic converters, inject different levels of harmonics to the power grid. The harmonic sources can be divided into the following categories in a power grid [5-7]:  Electrical machines (in point of physically distortion creation) and non-linear loads are defined as non- semiconductor based harmonic sources.  Motor control devices, HVDC power transmission system, FACTS systems, electrical vehicle and cyclo-converters in low-speed and high torque machines are called semiconductor based harmonic sources. As introduced in [8], Adjustable Speed Driver (ASD), capacitors, circuit breakers, fuses, conductors, lighting, meters, transformers, rotating machines, telephone interference and protective relaying are affected by harmonics distortion. International standards for power quality such as IEEE519, IEC 61000, and EN 50160 established admissible harmonic levels and their emission limits [9-11]. A prior research on power system harmonics has been done in 1983 by an IEEE working group [12]. Along with 50 years of operational reports, the effect of harmonics has been categorized in some clauses. The 8th clause of this report points to relay malfunction by harmonics [5]. This topic is notable when the decision-making process depends on the rate of waveform distortion. This distortion affects the operating point of relays and consequently the decision-making system [13]. In protective relays, most of the measurement units determine the peak value and the zero crossing of the voltage and current waveforms. Therefore, the harmonics can cause mis-operation of the microprocessor- based relays [14]. This issue has notably been investigated in case of power distribution system protection conducted on static under-frequency and over-current relays [15]. Authors in [7, 16] presented analytical approaches based on relay characteristics and performed theoretical and experimental analysis to derive phase angles of low-order harmonics on the tripping time of overcurrent relays. Harmonic order and amplitudes are determinant. Authors in [17] indicated the effects of noisy electrical signals on digital relays and they proposed a filter capable of filtering selective harmonic components. A low-pass filter is a well known yet primitive approach to deal with power system harmonics. To eliminate low-order harmonics, low-pass filter with sharp cut-off frequency is required and that frequency should be set among the fundamental and the next order harmonics [1]. Setting the cut-off frequency in a narrow frequency band makes the design complex. Although the process of low-pass filter design is theoretically effortless, it is practically not economical in order to filter low order harmonics. Any distortion may lead to unexpected relay operation. This fact is apparent especially in overcurrent relays [17, 19]. The fast operation of the relays is influenced by harmonic current which is less or more than nominal [20, 21]. This deviation is limited to 15%. In over current relays and with the presence of low-order harmonics, they can operate up to 60% faster than Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2178 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … normal and standard condition [22]. Harmonics also have undeniable effects on frequency protection relays. In under frequency relays, the 10% THD, increases operating time twice. Therefore, it is recommended to limit the harmonic pollution to 5% [23]. Also, harmonics have a severe impact on differential relays. Figure 1shows the basic components of the typical digital relay which LP filter situates after surge protective circuits Although the prevalent low-pass filtering might be suitable for high-frequency harmonics, the performance and cost are considerable especially in filtering low-order harmonics [24]. There are many detection and calculation algorithms and techniques operated on filtered- discrete waveforms to estimate operation territory such: Mann & Morrison method [25], Rockefeller & Udren technique [26], Fourier and Walsh-based technique [27, 28], least-square methods [29], differential equation technique and travelling- wave-based methods but they are applied in Central Processing Units (CPU) [30]. The Moving Average Window (MAW) technique is a numeric-mathematical based transform function to simulate the waveform with selective resolution. Also, the averaging technique used to model discrete spaces to continuous one such as switches in converters [18]. In low pass filter design, setting sharp cutting frequency makes the design process complex and costly. The proposed method is an efficient technique to estimate signals with variable resolution. In this article, the MAW technique in two schemes is compared with the common first and second order low-pass filters which are applied in protective relays. Figure 1 describes MAW embedded in the protective relay as a low-pass filter. The MAW has the sampling frequency as the only tuning parameter and its response is comparable with the second order filter design with two tuning parameters. The design simplicity and the low computational burden is the two main advantages of the proposed method.  Mann & Morrison method,  Rockefeller & Udren technique,  Fourier and Walsh-based technique,  Least-square methods,  Differential equation technique  Travelling-wave-based methods  First order  Second order  Moving Averaging Window (The Proposed Method) Fig. 1. Basic components of a digital relay and the situation of Moving Averaging Window in filtering process. Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2179 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … II. PROPOSED METHOD In order to design a conventional low-pass filter, the optimal tuning of a sharp and accurate cutoff frequency is required, something complicated and cost-effective in implementation. The Moving Averaged Window (MAW) technique can estimate the fundamental only through selecting a proper sampling frequency. A. Moving Average Window Many studies have been published on MAW applications in signal processing and modeling scope. MAW applications are summarized in the estimation of waveforms with various degree of resolution and rate of processing power [18]. The MAW technique is capable of fault detection, direction estimation, and faulty phase identification [31-33]. This technique is less sensitive to noise for sampling frequency around the noise frequency. The proposed technique defines a transfer function in continuous time-domain mode. Equation (1) shows MAW transfer function wherein x(t) is under process waveform. ( )x t is averaged waveform along window length of TS. 1 ( ) ( ) S t t T S x t x d T      (1) The frequency response of (1) is shown in (2) which is the Fourier transform of impulse response of (1) and is called low- pass filter frequently [31]. 2( ) sinc( ) 2 STj STX j e      (2) where sinc( )=sin( )/   , Ts, and fs are sampling time and frequency respectively and fs defines the sampling/window length. The moving window is demonstrated in Figure 2. Figure 3 shows the MAW‘s frequency response. In Figure 3(a), ω frequency (nominal frequency) is rejected by the filter in π radian phase and other frequencies are attenuated except dc. Selecting TS, the filter output can be zero at fS. As shown in Figure 3(b), the filter’s phase response lags. The shown characteristic in amplitude and phase could lead up to article’s main purpose of the usage of MAW as a low-pass filter. B. Filter Response to Harmonic Polluted Signal By applying polluted waveform in (1) as the input variable, (3) is obtained, which shows mathematical-time domain of harmonically distorted signals with odd-order frequencies. 0 2 1 2 1 2 1 1 ( ) sin(2 ) n k k k k x t x x f t        (3) where, x1 is the fundamental frequency, x2k+1 and f2k+1 are odd- harmonic order amplitude and frequencies, respectively and φ2k+1 is corresponded phase angle. Applying (1) to (3), the averaged equation is derived as (4): 0 2 1 2 1 2 1 2 1 1 ( ) sin(2 ) n k k k k S k x t x a f t f T           (4) 2 1 2 1 2 1sinc( )k k k Sa x f T   (5) where πf2k+1Ts is phase lagging in filtered signal. In Figure 4, the amplitudes of 13 harmonic orders are illustrated by different averaging window frequencies. In Figure 4 we see that the amplitude of harmonics varies by sampling frequency. The main idea lies there. Since high sampling frequency leads to the amplifying of harmonics order, in some frequencies, harmonic orders tend to be zero. The zero crossing frequencies are an integer multiple of each harmonic frequency. Table I shows the amplitude of each zero-crossing frequency which plays an important role in MAW philosophy. Moving Window 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 Ts Time, sec A m pl itu de , p .u Fig. 2. Moving Window applied to distorted input waveform containing harmonics and Gaussian noise with variance 0.01 M ag ni tu de Ph as e, D eg Frequency, Hz (a) (b) Frequency, Hz Fig. 3. The proposed filter characteristics based on MAW. (a) Magnitude response. (b) Phase response 0 200 400 600 800 1000 1200 1400 1600 1800 2000 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 Sampling Frequency, Fs, Hz Amplitude P.u. Fig. 4. The amplitude of response of MAW by varying sampling frequency up to 13th harmonic order Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2180 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … In fact, averaging theory waived harmonic order which is swept with the corresponding frequency. In Table I, by selecting 150 Hz as MAW sampling frequency, weakening the related harmonic order (3rd and 9th) is achieved. This fact is carried out on other harmonic orders. Table II shows the phase lagging of each selected sample frequency. As it can be derived from Tables I and II, by applying fs=150 Hz, 300-degree phase lagging in fundamental and 17.32% decaying in magnitude are achieved. In section IV, two filtering strategies are discussed and the ability of MAW would be evaluated in comparison with first and second order low-pass filter. TABLE I. THE AMPLITUDE OF HARMONIC ORDER IN ZERO-CROSSING FREQUENCY Harmonic Order fs (Hz) 1st 3rd 5th 7th 9th 11th 13th 150 0.827 0 0.1654 0.1181 0 0.0752 0.0636 250 0.9355 0.5045 0 0.2162 0.1039 0.085 0.1164 350 0.9667 0.7241 0.3484 0 0.1935 0.1974 0.0744 450 0.9798 0.827 0.5642 0.263 0 0.1674 0.217 550 0.9864 0.882 0.6931 0.455 0.2103 0 0.1456 650 0.9903 0.9146 0.7738 0.5868 0.3784 0.1748 0 TABLE II. PHASE LAGGING DEGREE IN APPLYING MAW WITH VARIABLE SAMPLING FREQUENCY Harmonic Order fs (Hz) 1st 3rd 5th 7th 9th 11th 13th 150 300 180 60 300 180 60 300 250 324 252 180 108 36 324 252 350 334.28 282.85 231.43 180 128.57 77.14 25.71 450 340 300 260 220 180 140 100 550 343.64 310.91 278.18 245.45 212.73 180 147.27 650 346.15 318.46 290.77 263.07 235.38 207.69 180 III. SIMULATION AND RESULT A. DC Offset Elimination According to Figure 3 and to the frequency response of MAW, the DC component is amplified to 100%. The MAW technique is capable of eliminating DC offset of input waveforms easily through surviving MAW in fundamental frequency. Figure 5 shows the process of DC offset elimination which is proposed by MAW with a fundamental frequency sweep, e.g. fs=50 Hz, DC offset totally removed. MAW Fs=50 Hz + _ Inputs DC offset DC eliminated inputs Fig. 5. Proposed method for eliminating of DC offset on input waveforms before the main filtering B. Sampling Frequency, fs Sampling frequency has a crucial role in estimating fundamental waveform and harmonic’s content. The sampling frequency is obtained by inherent characteristics of local harmonic properties. In order to filter the harmonic polluted waveforms, two harmonic elimination strategies are defined. Figure 6 shows the flowchart in explanation of these strategies. As it shows, the fs has a key role in this filtering process. 1) Selective Harmonic Elimination Due to selective harmonic order elimination, it is enough to select corresponding harmonic frequency as the sampling frequency. Table I illustrates selective harmonic scheme. 2) General Harmonic Attenuation As it can be derived from Table I, by selecting fs=150 Hz, the harmonic reduction is started with the aim of maximum elimination of 3rd harmonic order. This strategy is defined as selective harmonic weakening by concentration on 3rd harmonic order elimination. In General Harmonic Attenuation strategy, the objective function sets on achieving minimum THD. In this strategy, fs tunes to minimize output THD by applying MAW. 3) MAW Technique Evaluation In order to prove MAW sufficiency, the simulation is performed by comparing two conventional low-pass filters and MAW. Table III shows 8 case studies of polluted waveforms. The first two cases of Table III present the most common harmonic existence with notable amplitude in the power system. The THD would be considered as a criterion pollution rate [32, 33]. The simulation was performed on first order low-pass filter (LPF-1), second order low-pass filter (LPF-2), MAW-General, and MAW-Selective. Table IV presents the result of filtered case studies which have been introduced in Table III. The first two columns are the result of the first and second order traditional low-pass filter. These filter characteristics are given in Appendix. The last two columns result applied general and selective MAW to these case studies. Methodologies have been discussed in section IV. According to first order filtering results in case 1, the THD decreases from 60.21% to 23.99% and fundamental weakens from 100% to 71%. The second order low-pass filter is set on high-quality factors and the full amount of fundamental amplitude. Two parameters (ωc, ξ) should be set and therefore, two poles add to system transfer function order level. The result shows that by applying second order LPF, fundamental is kept to 100% and THD factor decreased to 5.98%. Filtering with MAW-General strategy has been done by setting fs=144 Hz to achieve minimum THD. In the case of presence of various harmonic orders, the general MAW suggested for general decreasing harmonic content. In addition, Table IV shows the ability of MAW technique in both general and selective harmonic decaying. The result is that MAW has the best performance in comparison to traditional low-pass filtering methods. As for the various order of harmonics, the MAW-General scheme results better than MAW -Selective. Similarly, in the polluted waveform with specific harmonic order, the MAW-Selective has the best performance. Figure 7 illustrates the effect of LPF1, LPF2, MAW-General and MAW-Selective on case 1. FFT analysis is calculated with transient-free waveforms through THD, TIHD and TSHD. Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2181 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … Strategies and fs Determination MAW fs Filtered Waveform Amplitude DC offset elimination + _ MAW 50 Hz Is strategy selective? fs fs=144Hz Yes No Phase-shifting Calculation Start Input Waveform Output Decision Symbol Output Symbol Process Symbol Manual input Symbol Merge Symbol Merging Amplitude and Phase Fig. 6. Proposed technique flowchart TABLE III. 8 CASE STUDIES EVALUATING BY MAY TECHNIQUE Fundamental 3rd 5th 7th 9th 11th 13th THD C as es 1 100% 50% 30% 10% 8% 6% 5% 60.21% 2 100% 50% 50% 10% 0 0 0 71.41% 3 100% 50% 0 0 0 0 0 50% 4 100% 0 30% 0 0 0 0 30% 5 100% 0 0 30% 0 0 0 30% 6 100% 0 0 0 30% 0 0 30% 7 100% 0 0 0 0 30% 0 30% 8 100% 0 0 0 0 0 30% 30% IV. CONCLUSION The signal processing of digital relays has mainly focused on algorithms to find out and estimate the operation criteria. Many techniques and methods have been proposed. It is clear that pre-filtering produces better results when using a low-pass filter. The design of sharp cutting frequency among fundamental and first harmonic order, causes complexity and cost. Therefore the Moving Average Window was proposed to eliminate DC offset, noise and attenuate harmonic content. The proposed technique could satisfyingly substitute the traditional first and second order low-pass filter. The evaluation test was done in 8 case studies and the results are tabulated and illustrated. The MAW has the ability to eliminate specified harmonic order or attenuate whole harmonic order. Thus, two filtering strategies are defined and verified with the 8 case studies. Low computational burden and only one tuning parameter (fs) are two main advantages of MAW. TABLE IV. RESULT OF 8 CASE STUDIES FILTERING BY 4 STRATEGIES: LPF-1, LPF-2, MAW-GENERAL AND MAW SELECTIVE LPF-1 (%) LPF-2 (%) MAW-General (%) MAW-Selective (%) Fund. 3rd 5th THD Fund. 3rd 5th THD Fund. 3rd 5th THD Fund. 3rd 5th THD C as es 1 71 16 6 23.99 100 5.85 1.22 5.98 fs=144 Hz fs=150 Hz 82.7 0 4.97 6.2 81.315 1.994 4.05 5.81 fs=250 Hz 93.55 25.23 0 27.09 2 71 16 10 26.39 100 5.85 2.04 6.20 fs=136 Hz fs=150 Hz 82.7 0 8.27 10.1 79.20 4.58 4.21 8 fs=250 Hz 93.54 25.22 0 27.06 3 71 16 0 22.36 100 5.85 0 5.85 82.7 0 0 0 82.7 0 0 0 4 71 0 6 8.32 100 0 1.22 1.22 93.55 0 0 0 93.55 0 0 0 5 71 0 0 6 100 0 0 0.62 96.67 0 0 0 96.67 0 0 0 6 71 0 0 4.69 100 0 0 0.37 97.98 0 0 0 97.98 0 0 0 7 71 0 0 3.84 100 0 0 0.25 98.64 0 0 0 98.64 0 0 0 8 71 0 0 3.25 100 0 0 0.18 99.03 0 0 0 99.03 0 0 0 Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2182 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … a 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 5 10 15 20 25 30 0 0.2 0.4 0.6 0.8 1 THD% = 60.2084 TIHD% = 0.046798 TSHD% = 0.011681 -1.5 -1 -0.5 0 0.5 1 1.5 b -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 0 5 10 15 20 25 30 0 0.2 0.4 0.6 0.8 THD% = 23.9932 TIHD% = 0.023194 TSHD% = 0.023993 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 c -2 -1 0 1 2 -2 -1 0 1 2 0 5 10 15 20 25 30 0 0.5 1 THD% = 5.9823 TIHD% = 0.42009 TSHD% = 0.77405 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 d -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 0 5 10 15 20 25 300 0.5 1 THD% = 5.814 TIHD% = 1.5935e-13 TSHD% = 3.1673e-13 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 f -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 0 5 10 15 20 25 30 0 0.5 1 THD% = 27.0933 TIHD% = 2.272e-13 TSHD% = 2.7745e-13 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 Fig. 7. Simulation result of case 1: (a) Original waveform (b) LPF1 waveform both original and without transient-free waveform (c) LPF2 waveform both original and transient-free waveform (d) MAW-General with Fs=144 Hz both original and transient-free waveform (e) MAW-Selective with Fs=150 Hz both original and transient-free waveform (f) MAW-Selective with Fs=250 Hz both original and transient-free wavefor APPENDIX FIRST ORDER LOW-PASS FILTER 1 ( ) 1 c H s s    Where, ωc=2πfc and fc is the cutting frequency which is set on 50 Hz for achieving maximum THD SECOND ORDER LOW-PASS FILTER 2 2 1 ( ) 2 1 c c H s s s       Where ωc=2πfc and fc is the cutting frequency which is set on 50 Hz. ξ is damping ration which is set on 0.5 for achieving better quality factor. REFERENCES [1] W. Rebizant, J. Szafran, A. Wiszniewski, Digital signal processing in power system protection and control, Springer Science & Business Media, 2011 [2] M. H. J. Bollen, I. Y. H. Gu, Signal processing of power quality disturbances, John Wiley & Sons, 2006 [3] S. Upadhyaya, S. Mohanty, C. N. Bhende, “Hybrid Methods for Fast Detection and Characterization of Power Quality Disturbances”, Journal of Control, Automation and Electrical Systems, Vol. 26, No. 5, pp. 556- 566, 2015 [4] S. Khokhar, A. A. B. Mohd Zin, A. S. B. Mokhtar, M. Pesaran, “A comprehensive overview on signal processing and artificial intelligence techniques applications in classification of power quality disturbances”, Renewable and Sustainable Energy Reviews, Vol. 51, pp. 1650-1663, 2015 [5] R. Yacamini, “Power system harmonics. Part 1. Harmonic sources”, Power Engineering Journal, Vol. 8, No. 4, pp. 193-198, 1994 [6] J. Arrillaga, N. R. Watson, “Harmonic Sources”, in Power System Harmonics, John Wiley & Sons, pp. 61-142, 2004 [7] M. A. S. Masoum, A. S. Masoum, A. S. Masoum, “Influence of low order harmonic phase angles on tripping time of overcurrent relays”, Engineering, Technology & Applied Science Research Vol. 7, No. 6, 2017, 2177-2183 2183 www.etasr.com Khodabakhshi-Javinani and Askarian-Abyaneh: Evaluation of Moving Average Window Technique … 40th North American Power Symposium, pp. 1-5, Calgary, Canada, 2008 [8] V. E. Wagner, J. C. Balda, D. C. Griffith, A. McEachern, T. M. Barnes, D. P. Hartmann, D. J. Phileggi, A. E. Emannuel, W. F. Horton, W. E. Reid, R. J. Ferraro, W.T. Jewell, “Effects of harmonics on equipment”, IEEE Transactions on Power Delivery, Vol. 8, No. 2, pp. 672-680, 1993 [9] IEEE Recommended Practice and Requirements for Harmonic Control in Electric Power Systems, IEEE Std 519-2014 (Revision of IEEE Std 519-1992), pp. 1-29, 2014 [10] I. E. C., Electromagnetic compatibility (EMC), 2014 [11] H. Markiewicz, A. Klajn, Voltage characteristics of electricity supplied by public distribution networks, Copper Development Assosiation, 2004 [12] Power Harmonics, “Power System Harmonics: An Overview”, IEEE Transactions on Power Apparatus and Systems, Vol. PAS-102, No. 8, pp. 2455-2460, 1983 [13] W. F. Horton S. Goldberg, “The Effect of Harmonics on the Operating Points of Electromechanical Relays”, IEEE Power Engineering Review, Vol. PER-5, No. 5, pp. 45-46, 1985 [14] M. Shao, P. Poonpun, W. T. Jewell, “An advanced methodology for under-frequency load shedding relay testing”, Transmission and Distribution Conference and Exposition, pp. 1-4, Chicago, USA, 2008 [15] J. F. Fuller, E. F. Fuchs, D. J. Roesler, “Influence of harmonics on power distribution system protection”, IEEE Transactions on Power Delivery, Vol. 3, No. 2, pp. 549-557, 1988 [16] M. A. S. Masoum, S. M. Islam, K. Tan, N. X. Tung, “Impact of harmonics on tripping time of overcurrent relays”, 2007 Australasian Universities Power Engineering Conference, pp. 1-5, 2007 [17] J. Lazaro, J. F. Minambres, M. A. Zorrozua, “Selective estimation of harmonic components in noisy electrical signals for protective relaying purposes”, International Journal of Electrical Power & Energy Systems, Vol. 56, pp. 140-146, 2014 [18] R. W. Erickson, D. Maksimovic, Fundamentals of Power Electronics: Springer, 2001 [19] V. Gurevich, Cyber and Electromagnetic Threats in Modern Relay Protection, Crc Press, 2014 [20] J. M. Ho, C. C. Liu, “The effects of harmonics on differential relay for a transformer”, 16th International Conference and Exhibition on Electricity Distribution, Part 1: Contributions. CIRED, Vol. 2, pp. 5, 2001 [21] M. Kezunovic, B. Kasztenny, “Design optimization and performance evaluation of the relaying algorithms, relays and protective systems using advanced testing tools”, 21st IEEE International Conference on Power Industry Computer Applications, pp. 309-314, 1999 [22] P. M. Donohue, S. Islam, “The Effect of Nonsinusoidal Current Waveforms on Electromechanical and Solid-State Overcurrent Relay Operation”, IEEE Transactions on Industry Applications, Vol. 46, No. 6, pp. 2127-2133, 2010 [23] E. F. Fuchs, D. J. Roesler, M. A. S. Masoum, “Are harmonic recommendations according to IEEE and IEC too restrictive?”, IEEE Transactions on Power Delivery, Vol. 19, No. 4, pp. 1775-1786, 2004 [24] A. T. Johns, S. K. Salman, Digital protection for power systems: IET, 1997 [25] B. J. Mann, I. Morrison, “Digital calculation of impedance for transmission line protection”, IEEE Transactions on Power Apparatus and Systems, Vol. PAS-90, No. 1, pp. 270-279, 1971 [26] G. Rockefeller, “Fault protection with a digital computer”, IEEE Transactions on Power Apparatus and Systems, Vol. PAS-88, No. 4, pp. 438-464, 1969 [27] A. Johns, M. Martin, “Fundamental digital approach to the distance protection of EHV transmission lines”, Proceedings of the Institution of Electrical Engineers, Vol. 125, No. 5, pp. 377-384, 1978 [28] J. W. Horton, “Use of Walsh functions for high-speed digital relaying”, IEEE PES Summer Meeting, San Francisco, Paper No.A75-582-7, pp.1- 9. 1975 [29] A. Ranjbar, B. Cory, “An improved method for the digital protection of high Voltage transmission lines”, IEEE Transactions on Power Apparatus and Systems, , Vol. 94, No. 2, pp. 544-550, 1975 [30] Q.-H. Wu, Z. Lu, T. Ji, Protective relaying of power systems using mathematical morphology, Springer Science & Business Media, 2009 [31] S. Dambhare, S. A. Soman, M. C. Chandorkar, “Current Differential Protection of Transmission Line Using the Moving Window Averaging Technique”, IEEE Transactions on Power Delivery, Vol. 25, No. 2, pp. 610-620, 2010 [32] S. M. Hashemi, M. Tarafdar Hagh, H. Seyedi, “Transmission-Line Protection: A Directional Comparison Scheme Using the Average of Superimposed Components”, IEEE Transactions on Power Delivery, Vol. 28, No. 2, pp. 955-964, 2013 [33] A. K. Pradhan, A. Routray, S. R. Mohanty, “A Moving Sum Approach for Fault Detection of Power Systems”, Electric Power Components and Systems, Vol. 34, No. 4, pp. 385-399, 2006