Next Article in Journal
The Impact of Selected Biofuels on the Skoda Roomster 1.4TDI Engine’s Operational Parameters
Next Article in Special Issue
Short-Term Load Dispatching Method for a Diversion Hydropower Plant with Multiple Turbines in One Tunnel Using a Two-Stage Model
Previous Article in Journal
A Decentralized Architecture Based on Cooperative Dynamic Agents for Online Voltage Regulation in Smart Grids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Phasor Estimation Algorithm Based on a Least Squares Method

Department of Electrical Engineering, Myongji University, Yongin 17058, Korea
*
Author to whom correspondence should be addressed.
Energies 2019, 12(7), 1387; https://doi.org/10.3390/en12071387
Submission received: 11 March 2019 / Revised: 2 April 2019 / Accepted: 7 April 2019 / Published: 10 April 2019
(This article belongs to the Special Issue Electric Power Systems Research 2019)

Abstract

:
This paper proposes an adaptive phasor estimation algorithm based on a least square method that can suppress the adverse effect of an exponentially decreasing DC offset component in a phasor estimation process. The proposed algorithm is composed of three stages: a basic least squares model, a time constant calculation, and an adaptive least squares model. First, we use the basic least squares model to estimate the parameter of the DC offset component in the fault current signal. This model is designed to incorporate fundamental frequency, and harmonic and constant components. Second, we use the estimated parameter to calculate the time constant of the DC offset component. Third, we redesign a least squares model that incorporates fundamental frequency, harmonic components, and exponential function of the DC offset component. Since this model incorporates the exponential function of the DC offset component contained in the fault current signal, it estimates the phasor of the correct fundamental frequency component without influence of the DC offset component. We evaluated the performance of the proposed algorithm using computer generated signals and EMTP simulation signals. The evaluation results show that the proposed algorithm can effectively suppress the adverse influence of the exponentially decaying DC offset component.

1. Introduction

Accurate phasor estimation of current and voltage signals is very important for the protection and control of a power system. When a power system fault occurs, the current signal may include noise components such as an exponentially decaying DC offset and high-frequency due to traveling waves in addition to the fundamental frequency component. This noise degrades phasor estimation performance. The high-frequency components generated by a traveling wave may be both harmonic and non-harmonic. The most commonly used discrete Fourier transform (DFT) and least squares method (LSM) for phasor estimation can eliminate harmonic components. However, non-harmonic components are not completely removed using these phasor estimation methods and must be removed because they can cause errors in the phasor estimation result. The frequency band of the high-frequency components generated by a traveling wave is determined by the length of a power system’s transmission line [1]. Therefore, the cut-off frequency of the low-pass filter used for a power system’s protection relay should be chosen to be lower than the lowest frequency of high-frequency components generated by the traveling wave, thus suppressing the phasor estimation error caused by non-harmonic components generated by the traveling wave.
The exponentially decaying DC offset component contained in the fault current signal is not removed during the phasor estimation process and the filtering process by the low-pass filter, since the DC offset component is a non-periodic signal and its frequency spectrum includes all frequencies. Therefore, the DC offset component should be carefully taken into consideration in calculating the phasor of the fundamental frequency component.
Many techniques have been proposed for estimating the phasor of the correct fundamental frequency component without the adverse effect of the exponentially decaying DC offset component of the fault current signal. These can be classified into three main groups.
In the first group, the exponentially decaying DC offset component contained in the fault current signal is first filtered out and the phasor of the fundamental frequency component is estimated from the filtered signal using a DFT. In [2], a digital mimic filter is proposed to suppress the DC offset component in a fault current signal. The digital mimic filter can completely remove the DC offset component only when the actual time constant of the fault current signal’s DC offset component is the same as the mimic filter’s time constant. In [3], a Kalman filter is proposed to eliminate the DC offset component, only when the time constant of the DC offset component is the same as that modeled in a state transition matrix. In [4], an algorithm that uses the average value of the fault current signal for one cycle is proposed. In [5,6], algorithms that use the integral value of the fault current signal during one cycle are proposed. These algorithms [4,5,6] cannot completely eliminate the DC offset component if its time constant is small.
The second group uses a DFT to estimate the phasor of the fundamental frequency component and then calculate and remove the error due to the DC offset component in the DFT’s estimated phasor value. In [7,8,9,10,11], the error due to the DC offset component in the DFT’s estimated phasor value is calculated and eliminated by using three successive estimated phasor values. However, these methods take one cycle plus two additional samples to estimate the correct phasor and are susceptible to high frequency noise. In [12], the error due to the DC offset component in the DFT’s estimated phasor value is calculated and eliminated using the output of another DFT filter tuned at the n-th harmonic component. However, the performance of this algorithm is highly affected by random noise. In [13], a partial summation algorithm is proposed to eliminate the error due to the DC offset component in the DFT’s estimated phasor value. However, if the DC offset component has a small time constant, this algorithm cannot completely eliminate the error due to the DC offset component in the DFT’s estimated phasor value. In [14,15,16], the error due to the DC offset component in the DFT’s estimated phasor value is calculated and eliminated using the outputs of the even- and odd-sample-set DFTs. In [17], the samples used for the DFT computation are split into four groups. The DFT then obtains an estimated phasor value from each of the four groups. Then, the parameters of the DC offset component are obtained from the four phasor results. In [18], the error due to the DC offset component in the DFT’s estimated phasor value is calculated and eliminated by combining the DFT outputs from even and odd sample-sets obtained for a one cycle data window decimated by two and by four. In [19,20], algorithms using an average value over one cycle are proposed to estimate the error due to the DC offset component in the DFT’s estimated phasor value. However, these methods [14,15,16,17,18,19,20] are disadvantageous because the high frequency noise components are amplified by the subtraction process that calculates the error due to the DC offset component.
The third group uses the least squares methods. These methods find a coefficient that minimizes the error between the input signal and the predefined signal. Therefore, the predefined signal must be selected appropriately. As in the techniques of the first two groups, it is necessary to consider both the harmonic components and the exponentially decaying DC offset component to achieve an accurate phasor estimate. In [21,22,23,24,25,26], the DC offset component is approximated to the first two or three terms of its Taylor series expansion. However, as in the first group of techniques, these methods cannot accurately approximate the DC offset component when the time constant of the DC offset component is small, so the estimated phasor of the fundamental frequency component will include the error due to the DC offset component.
In this paper, we propose an adaptive phasor estimation method based on a least squares method. This method estimates the time constant of the exponentially decaying DC offset component contained in a fault current signal and then redesigns the least squares model using the estimated time constant. This approach estimates the phasor of the fundamental frequency component without the influence of the DC offset component. We tested the performance of the proposed algorithm by running it against computer-generated signals and EMTP simulation signals. To evaluate the algorithm’s performance, we compared the results produced by the proposed algorithm to the results obtained from the conventional DFT, the even/odd DFT [14], the modified DFT [10], and the conventional LSM [21].

2. Proposed Phasor Estimation Algorithm

2.1. Basic Least Squares Model

When a fault occurs in a power system, the current signal can be represented as the combination of a fundamental frequency component, a harmonic component, and an exponentially decaying DC offset component, as in:
y k = B 0 e k Δ t / τ + n = 1 M B n cos n 2 π N k + θ n   ,
where, B 0 and τ are the magnitude and the time constant of an exponentially component, B n and θ n are the amplitude and the phase angle of the n-th harmonic component, Δ t is the sampling interval, N is the number of samples per cycle, M is the maximum order of the harmonic components included in the current signal, and k is the sample index.
The harmonic components in Equation (1) can be expressed as follows using the trigonometric formula:
n = 1 M B n cos n 2 π N k + θ n = B 1 cos θ 1 · c o s 2 π N k B 1 sin θ 1 · s i n 2 π N k +   B M cos θ M · c o s M 2 π N k B M sin θ M · s i n M 2 π N k
The exponentially decaying DC offset component in Equation (1) can be expressed as the Taylor series expansion:
B 0 e k Δ t / τ = B 0 B 0 k Δ t τ + B 0 2 k Δ t τ 2  
The conventional least squares methods are designed using the first two or three terms of the Taylor series expansion of Equation (3) in order to consider the exponentially decaying DC offset component in the fault current signal. If the time constant of the DC offset component is small, this component cannot be approximated accurately.
In this paper’s proposed technique, we estimate the time constant of the DC offset component and then redesign the least squares model matrix to estimate an accurate phasor of the fundamental frequency component. In order to estimate this time constant, the basic least squares model is designed to incorporate fundamental frequency, harmonic and constant components as in:
y k N + 1 y k 2 y k 1 y k = cos 2 π N sin 2 π N cos M 2 π N sin M 2 π N 1 cos 2 π N N 2 sin 2 π N N 2 cos M 2 π N N 2 sin M 2 π N N 2 1 cos 2 π N N 1 sin 2 π N N 1 cos M 2 π N N 1 sin M 2 π N N 1 1 cos 2 π N N sin 2 π N N cos M 2 π N N sin M 2 π N N 1 B 1 cos θ 1 B 1 sin θ 1 B M cos θ M B M sin θ M B 0 k
Y A N X N N × 1 N × 2 M + 1 2 M + 1 × 1
where, Y is an N-by-1 vector of the input signals, A is an N-by-(2M + 1) model matrix for the basic least squares method, and X is a (2M + 1)-by-1 vector of coefficients.
Equation (4) can be briefly expressed as:
Y = A · X
The unknown vector of coefficients in Equation (5) can be obtained as follows:
X = A T · A 1 · A T · Y   ,
where, A T is the transpose of the model matrix A .
If the model matrix for the basic least squares method is designed as shown in Equation (4) and the input signal is equal to Equation (1), the estimated vector of coefficients will include the influence of the exponentially decaying DC offset component. This is because the DC offset component is not considered in the model matrix.

2.2. Time Constant Calculation

The coefficient of the constant component ( B 0 k ) among the coefficients estimated by the basic least squares model is expressed as:
B 0 k = B 0 N · e Δ t τ · k N + 1 e Δ t τ · k 1 e Δ t τ = B 0 N · 1 e Δ t τ · N 1 e Δ t τ e Δ t τ · k N + 1
As shown in Equation (7), the coefficient of the constant component is estimated as the mean value of the exponentially decaying DC offset component included in the data window of the basic least squares method. Appendix A contains proofs related to Equation (7).
The time constant ( τ ) of the exponentially decaying DC offset component is calculated as in Equation (8) using two successive coefficients of the constant component:
B 0 k B 0 k 1 = B 0 N · 1 e Δ t τ · N / 1 e Δ t τ B 0 N · 1 e Δ t τ · N / 1 e Δ t τ · e Δ t τ · k N + 1 e Δ t τ · k N + 2 = e Δ t / τ
Using Equation (8), the time constant of the exponentially decaying DC offset component can be obtained as:
τ = Δ t l n B 0 k B 0 k 1 = Δ t l n e Δ t / τ

2.3. Adaptive Least Squares Model

In order to estimate the exact coefficients without affecting the exponentially decaying DC offset component included in the current signal, the exponential function should be incorporated in the model matrix for the adaptive least squares method. Since the time constant of the DC offset component contained in the current signal is estimated in Equation (9), the adaptive least squares model can be designed as:
y k N + 1 y k 2 y k 1 y k = cos 2 π N sin 2 π N cos M 2 π N sin M 2 π N e Δ t τ cos 2 π N N 2 sin 2 π N N 2 cos M 2 π N N 2 sin M 2 π N N 2 e Δ t τ N 2 cos 2 π N N 1 sin 2 π N N 1 cos M 2 π N N 1 sin M 2 π N N 1 e Δ t τ N 1 cos 2 π N N sin 2 π N N cos M 2 π N N sin M 2 π N N e Δ t τ N B 1 cos θ 1 B 1 sin θ 1 B M cos θ M B M sin θ M B 0 k
Y A N X N N × 1 N × 2 M + 1 2 M + 1 × 1
where, Y is an N-by-1 vector of the input signals, A N is an N-by-(2M + 1) model matrix for new least squares method, and X N is a (2M + 1)-by-1 vector of new coefficients.
Equation (10) can be briefly expressed as:
Y = A N · X N
The unknown vector of coefficients in Equation (11) can be obtained as:
X N = A N T · A N 1 · A N T · Y
The newly designed least squares model takes the exponential terms in the model matrix into account so that the exponentially decaying DC offset component included in the input signal can be approximated accurately. Thus, the error due to the DC offset component is excluded from the phasor of the fundamental frequency component.
Figure 1 shows a flowchart of the proposed algorithm. A low-pass Butterworth filter prevents aliasing and removes non-harmonic components generated by traveling waves. The basic least squares model obtains the parameters of the exponentially decaying DC offset component contained in the fault current signal. The DC offset component’s time constant is calculated from these parameters. An adaptive least squares model is designed taking into account the time constant of the exponentially decaying DC offset. Finally, the phasor of the fundamental frequency component is accurately estimated using the adaptive least squares model.

3. Case Studies

We evaluated the proposed algorithm’s performance using computer-generated signals and EMTP-simulation signals. We compared these test results to those of the conventional DFT, the even/odd DFT [14], the modified DFT [10] and the conventional LSM [21].
To evaluate the accuracy of the proposed algorithm, we calculated the phasor estimation error to be:
Error = reference   RMS   value estimated   RMS   value reference   RMS   value × 100 % ,
where, the reference RMS value is the RMS value at 10 cycles after the fault.
In all cases, the sampling frequency was set to 7680 Hz (i.e., 128 sample per cycle in a 60 Hz System). The input signal was pre-processed by a second-order Butterworth low-pass filter with a gain of 0.1 at a stop-band cut-off frequency of 750 Hz in order to remove high-frequency noises and prevent aliasing errors. Although the sampling frequency of the algorithm is 7680 Hz, the high-frequency noise that is higher than 750 Hz must be rejected by a preprocessing filter because the frequency band generated by a traveling wave in the assumed 100 km transmission line is 750~1500 Hz [1].
In the conventional least squares method, the model matrix is designed to include harmonic and exponentially decaying DC offset components. The harmonic components are included up to the 12th order and the exponentially decaying DC offset is derived from the first two terms of the Taylor series expansion in Equation (3). The data window length of the conventional least squares method is one cycle. In the proposed algorithm, the model matrix of the basic least square is designed to include the harmonic and constant components. The harmonic components are included up to the 12th order. The data window length of the proposed algorithm is one cycle.

3.1. Computer-Generated Signals

We evaluated the proposed algorithm’s performance with computer-generated signals consisting of sinusoidal and exponentially decaying DC offset components. We investigated the algorithm’s sensitivity to DC offset variation by running the tests with two different time constants (0.5/ f 0 and 5/ f 0 ).
Case 1. Test for basic signals—the test signal was:
y k = 100 e k Δ t / τ 100 cos 2 π N k   ,
Figure 2 shows the test results for a time constant of 0.5/ f 0 . As shown in Figure 2a, the result of the conventional DFT shows the largest transient overshoot and the largest oscillatory response. In Figure 2b, the proposed algorithm is seen to more accurately remove the adverse influence of the exponentially decaying DC offset than the other methods.
Figure 3 shows the test results for a time constant of 5/ f 0 . These results show the same performance advantage demonstrated in Figure 2 for a time constant of 0.5/ f 0 .
Case 2. Test for Random Noise: In this case, the focus is on the noise immunity of the proposed algorithm. Gaussian white noise is added to the signal of Case 1. The signal-to-noise ratio (SNR) is set to 40 dB. The test signal used in this case is:
y k = 100 e k Δ t / τ 100 cos 2 π N k + Gaussian   white   noise
Figure 4 and Figure 5 show the robustness of the proposed algorithm in the presence of Gaussian white noise. Figure 4 shows the test results for a time constant of 0.5/ f 0 .
As seen in Figure 4a, the result of the conventional DFT and the even/odd DFT show a very large transient overshoot and a very large oscillatory response. In Figure 4b, the proposed algorithm be seen to more accurately remove the adverse influence of the exponentially decaying DC offset than the other methods. Since the even/odd DFT and the modified DFT intrinsically use difference equations, an error occurs if a signal contains Gaussian white noise.
Figure 5 shows the test results for a time constant of 5/ f 0 . As in the previous results for time constant 0.5/ f 0 , the results for this time constant demonstrate that the proposed algorithm can more accurately remove the adverse influence of the exponentially decaying DC offset than the other methods.
The test results for the proposed algorithm show stable and superior results over the other methods, even when the signal contains Gaussian white noise.

3.2. EMTP Simulation Signals

The model system for the simulations is evaluated for a-g faults on a 154 kV, 100 km overhead transmission line with sources at both ends, as shown in Figure 6. The transmission line parameters used in the simulations are given in Table 1.
The equivalent source parameters used in the simulation are given Table 2. The fault current signals are generated by the EMTP.
The magnitude and time constant of the exponentially decaying DC offset component varies depending on fault conditions such as fault location, fault resistance, and fault inception angle.
Figure 7, Figure 8 and Figure 9 show the estimated RMS values and percent error of the fault current. The magnitude and time constant of the DC offset components depend on the fault distance from the relaying point to the fault point, so the performance of the proposed algorithm is evaluated by changing fault distance. A fault inception angle of 0° was used to maximize the magnitude of the DC offset component. Figure 7 shows the test results for a solid single-phase to ground fault at 10% of the total line length from the relaying point. As seen in Figure 7, the proposed algorithm accurately removes the adverse influence of the exponentially decaying DC offset more accurately than the conventional DFT, the even/odd DFT, and the conventional LSM, and also converges faster than the modified DFT.
Figure 8 shows the test results for solid single-phase to ground faults at 50% of the total line length from the relaying point. These results show the same performance advantage demonstrated in Figure 7 for the 10% distance.
Figure 9 shows the test results for solid single-phase to ground faults at 90% of the total line length from the relaying point. These results show the same performance advantages demonstrated in Figure 7 and Figure 8 for their respective distances.
The magnitude and frequency band of the high-frequency components varies depending on fault conditions such as fault location, fault resistance, and fault inception angle.
Figure 10, Figure 11 and Figure 12 show the estimated RMS values and percent error of the fault current. The magnitude and frequency band of the high-frequency components depend on the fault distance from the relaying point to the fault point; thus, as in the tests of Figure 7, Figure 8 and Figure 9, the performance of the proposed algorithm is evaluated by changing fault distance. A fault inception angle of 90° was used to maximize the magnitude of the high-frequency component.
Figure 10 shows the test results for a solid single-phase to ground fault at 10% of the total line length from the relaying point. As shown in Figure 10b, the proposed algorithm can more accurately estimate the phasor of the fundamental frequency component than the conventional DFT, the modified DFT, the even/odd DFT, and the conventional LSM.
Figure 11 shows the test results for a solid single-phase to ground fault at 50% of the total line length from the relaying point. These results show the same performance advantages demonstrated in Figure 10 for the 10% distance.
Figure 12 shows the test results for a solid single-phase to ground fault at 90% of the total line length from the relaying point. These results show the same performance advantages demonstrated in Figure 10 and Figure 11 for their respective distances.
The EMTP test results show that, the proposed algorithm can both remove the exponentially decaying DC offset and accurately estimate the phasor of the fundamental frequency component when the signal contains high-frequency noise components.

4. Conclusions

In this paper, we proposed an adaptive phasor estimation algorithm based on a least squares method. The proposed algorithm estimates the time constant of the exponentially decaying DC offset component in the fault current signal first and redesigns the model matrix for the least squares model using the estimated time constant, and finally estimates the phasor of the fundamental frequency component using the redesigned least squares model. The result is that the phasor of the fundamental frequency component can be accurately estimated without the influence of the exponentially decaying DC offset component. The performance of the proposed algorithm was evaluated using computer-generated signals and EMTP simulation signals. Several computer-simulated signals, which take into account the fundamental frequency, exponentially decaying DC offset, and Gaussian white noise components, are used. The evaluation results show that the proposed algorithm can effectively suppress the adverse influence of the exponentially decaying DC offset and Gaussian white noise components and that the phasor of the fundamental frequency component can be accurately and quickly estimated. Future research will be focused on research to reduce the computational complexity of the least squares method.

Author Contributions

Conceptualization, W.-J.K. and S.-H.K.; Methodology, W.-J.K. and S.-H.K.; Software, W.-J.K.; Validation, W.-J.K., S.-R.N. and S.-H.K.; Formal Analysis, W.-J.K. and S.-H.K.; Investigation, W.-J.K. and S.-R.N.; Writing-Original Draft Preparation, W.-J.K.; Writing-Review & Editing, W.-J.K., S.-R.N. and S.-H.K.

Funding

This research was supported by Korea Electric Power Corporation, Grant Number R17XA05-2 and R18XA01.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In the least squares method, the straight line y = b is fitted through the given points x 1 , y 1 ,   ,   x n , y n so that the sum of the squares of the distances of those points from the straight line is at a minimum, where the distance is measured in the vertical direction (the y-direction) [27].
The point on the line with abscissa x j has the b . Hence, its distance from x j , y j is y j b in Figure A1 and that sum of squares is:
ε = j = 1 n y j b 2
ε depends on b . A necessary condition for ε to be minimum is
ε b = b j = 1 n y j b 2 = 2 b j = 1 n 1 2 j = 1 n y j = 0
The above equation can be expressed as:
b = j = 1 n y j j = 1 n 1
If the input value y j is an exponentially decaying DC offset component, the coefficient value b is as follows:
b = j = 1 n B 0 · e Δ t τ j j = 1 n 1 = B 0 n 1 e Δ t τ n 1 e Δ t τ
Figure A1. Vertical distance of a point x j , y j from a straight line y = b .
Figure A1. Vertical distance of a point x j , y j from a straight line y = b .
Energies 12 01387 g0a1

References

  1. Swift, G.W. The spctra of fault-induced transients. IEEE Trans. Power Appar. Syst. 1979, 98, 940–947. [Google Scholar] [CrossRef]
  2. Girgis, A.A. A new kalman filtering based digital distance relay. IEEE Trans. Power Appar. Syst. 1982, 101, 3471–3480. [Google Scholar] [CrossRef]
  3. Benmouyal, G. Removal of DC-offset in current waveforms using digital mimic filtering. IEEE Trans. Power Deliv. 1995, 10, 621–630. [Google Scholar] [CrossRef]
  4. Eisa, A.A.A.; Ramar, K. Removal of decaying DC offset in current signals for power system phasor estimation. In Proceedings of the 43rd International Universities Power Engineering Conference, Padova, Italy, 1–4 September 2008. [Google Scholar] [CrossRef]
  5. Cho, Y.-S.; Lee, C.-K.; Jang, G.; Lee, H.J. An innovative decaying DC component estimation algorithm for digital relaying. IEEE Tran. Power Deliv. 2009, 24, 73–78. [Google Scholar] [CrossRef]
  6. Abdoos, A.A.; Gholamian, S.A.; Farzinfar, M. Accurate and fast DC offset removal method for digital relaying schemes. IET Gener. Transm. Distrib. 2016, 10, 1769–1777. [Google Scholar] [CrossRef]
  7. Gu, J.-C.; Yu, S.-L. Removal of DC offset in current and voltage signals using a novel fourier filter algorithm. IEEE Trans. Power Deliv. 2000, 15, 73–79. [Google Scholar] [CrossRef]
  8. Yu, S.-L.; Gu, J.-C. Removal of decaying dc in current and voltage signals using a modified fourier filter algorithm. IEEE Trans. Power Deliv. 2001, 16, 372–379. [Google Scholar] [CrossRef]
  9. Nam, S.-R.; Kang, S.-H.; Park, J.-K. An analytic method for measuring accurate fundamental frequency components. IEEE Trans. Power Deliv. 2002, 17, 405–411. [Google Scholar] [CrossRef]
  10. Hwang, J.K.; Song, C.K.; Jeong, M.G. DFT-based phasor estimation for removal of the effect of multiple DC components. IEEE Trans. Power Deliv. 2018, 33, 2901–2909. [Google Scholar] [CrossRef]
  11. Kim, W.-J.; Kagn, S.-H. A DC offset removal algorithm based on AR model. J. Int. Counc. Electr. Eng. 2018, 8, 156–162. [Google Scholar] [CrossRef]
  12. Sidhu, T.S.; Zhang, X.; Albasri, F.; Sachdev, M.S. Discrete-fourier-transform-based technique for removal of decaying DC offset from phasor estimates. IEE Proc. Gener. Transm. Distrib. 2003, 150, 745–752. [Google Scholar] [CrossRef]
  13. Guo, Y.; Kezunovic, M.; Chen, D. Simplified algorithms for removal of the effect of exponentially decaying DC-offset on the fourier algorithm. IEEE Trans. Power Deliv. 2003, 18, 711–717. [Google Scholar] [CrossRef]
  14. Kang, S.-H.; Lee, D.-G.; Nam, S.-R.; Crossley, P.A.; Kang, Y.-C. Fourier transform-based modified phasor estimation method immune to the effect of the DC offsets. IEEE Trans. Power Deliv. 2009, 24, 1104–1111. [Google Scholar] [CrossRef]
  15. Rahmati, A.; Adhami, R. An accurate filtering technique to mitigate transient decaying DC offset. IEEE Trans. Power Deliv. 2014, 29, 966–968. [Google Scholar] [CrossRef]
  16. Silva, K.M.; Nascimento, F.A.O. Modified DFT-based phasor estimation algorithms for numerical relaying applications. IEEE Trans. Power Deliv. 2018, 33, 1165–1173. [Google Scholar] [CrossRef]
  17. Yu, C.-S.; Huang, Y.-S.; Jiang, J.-A. A Full- and Half-Cycle DFT-based Technique for Fault Current Filtering. In Proceedings of the 2010 IEEE International Conference on Industrial Technology, Vina del Mar, Chile, 14–17 March 2010. [Google Scholar] [CrossRef]
  18. Silva, K.M.; Kusel, B.F. DFT based phasor estimation algorithm for numerical digital relaying. Electron. Lett. 2013, 49, 412–414. [Google Scholar] [CrossRef]
  19. Dadash, M.R.; Zhang, Z. A new DFT-based current phasor estimation for numerical protective relaying. IEEE Trans. Power Deliv. 2013, 28, 2172–2179. [Google Scholar] [CrossRef]
  20. Jafarpisheh, B.; Madani, S.M.; Shahrtash, S.M. A new DFT-based estimation algorithm using high-frequency modulation. IEEE Trans. Power Deliv. 2017, 32, 2416–2423. [Google Scholar] [CrossRef]
  21. Sachdev, M.S.; Baribeau, M.A. A new algorithm for digital impedance relays. IEEE Trans. Power Apparatus Syst. 1979, 98, 2232–2240. [Google Scholar] [CrossRef]
  22. Rahman, M.A.; Dash, P.K.; Dowton, E.R. Digital protection of power transformer based on weighted least square algorithm. IEEE Trans. Power Appar. Syst. 1982, 101, 4204–4210. [Google Scholar] [CrossRef]
  23. Isaksson, A. Digital protective relaying through recursive least-squares identification. IEE Proc. C Gener. Transm. Distrib. 1988, 135, 441–449. [Google Scholar] [CrossRef]
  24. Sachdev, M.S.; Nagpal, M. A recursive least error squares algorithm for power system relaying and measurement applications. IEEE Trans. Power Deliv. 1991, 6, 1008–1015. [Google Scholar] [CrossRef]
  25. Belega, D.; Petri, D. Dynamic synchrophasor estimation in the presence of decaying DC offsets. In Proceedings of the 2014 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Montevideo, Uruguay, 12–15 May 2014. [Google Scholar] [CrossRef]
  26. Jafarian, P.; Sanaye-Pasand, M. Weighted least error squares based variable window phasor estimator for distance relaying application. IET Gener. Transm. Distrib. 2011, 5, 298–306. [Google Scholar] [CrossRef]
  27. Kreyszig, E. Advanced Engineering Mathematics, 10th ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2008; pp. 872–873. ISBN 978-0-470-45836-5. [Google Scholar]
Figure 1. Flowchart of the Proposed Algorithm.
Figure 1. Flowchart of the Proposed Algorithm.
Energies 12 01387 g001
Figure 2. Test results for Case 1 ( τ :0.5/ f 0 ): (a) Estimated RMS values; (b) errors.
Figure 2. Test results for Case 1 ( τ :0.5/ f 0 ): (a) Estimated RMS values; (b) errors.
Energies 12 01387 g002
Figure 3. Test results for Case 1 ( τ :5/ f 0 ): (a) Estimated RMS values; (b) errors.
Figure 3. Test results for Case 1 ( τ :5/ f 0 ): (a) Estimated RMS values; (b) errors.
Energies 12 01387 g003
Figure 4. Test results for Case 2 ( τ :0.5/ f 0 ): (a) Estimated RMS values; (b) errors.
Figure 4. Test results for Case 2 ( τ :0.5/ f 0 ): (a) Estimated RMS values; (b) errors.
Energies 12 01387 g004
Figure 5. Test results for Case 2 ( τ :5/ f 0 ): (a) Estimated RMS values; (b) errors.
Figure 5. Test results for Case 2 ( τ :5/ f 0 ): (a) Estimated RMS values; (b) errors.
Energies 12 01387 g005
Figure 6. Single-line diagram of the model system.
Figure 6. Single-line diagram of the model system.
Energies 12 01387 g006
Figure 7. Test results for fault distance: 10%, fault inception angle: 0°: (a) Estimated RMS values; (b) errors.
Figure 7. Test results for fault distance: 10%, fault inception angle: 0°: (a) Estimated RMS values; (b) errors.
Energies 12 01387 g007
Figure 8. Test results for fault distance: 50%, fault inception angle: 0°: (a) Estimated RMS values; (b) errors.
Figure 8. Test results for fault distance: 50%, fault inception angle: 0°: (a) Estimated RMS values; (b) errors.
Energies 12 01387 g008
Figure 9. Test results for fault distance: 90%, fault inception angle: 0°. (a) Estimated RMS values; (b) errors.
Figure 9. Test results for fault distance: 90%, fault inception angle: 0°. (a) Estimated RMS values; (b) errors.
Energies 12 01387 g009
Figure 10. Test results for fault distance: 10%, fault inception angle: 90°: (a) Estimated RMS values; (b) errors.
Figure 10. Test results for fault distance: 10%, fault inception angle: 90°: (a) Estimated RMS values; (b) errors.
Energies 12 01387 g010
Figure 11. Test results for fault distance: 50%, fault inception angle: 90°. (a) Estimated RMS values; (b) errors.
Figure 11. Test results for fault distance: 50%, fault inception angle: 90°. (a) Estimated RMS values; (b) errors.
Energies 12 01387 g011
Figure 12. Test results for fault distance: 90%, fault inception angle: 90°. (a) Estimated RMS values; (b) errors.
Figure 12. Test results for fault distance: 90%, fault inception angle: 90°. (a) Estimated RMS values; (b) errors.
Energies 12 01387 g012
Table 1. Transmission line parameters.
Table 1. Transmission line parameters.
SequenceParameterValueUnits
Positive & NegativeR1, R20.0419Ω/km
L1, L20.8799mH/km
C1, C20.0128μF/km
ZeroR00.2294Ω/km
L02.6666mH/km
C00.0043μF/km
Table 2. Equivalent source parameters.
Table 2. Equivalent source parameters.
SequenceSource 1Source 2
S [GVA]X/RR [GVA]X/R
Positive & Negative5.010.72.17.4
Zero11.55.16.34.9

Share and Cite

MDPI and ACS Style

Kim, W.-J.; Nam, S.-R.; Kang, S.-H. Adaptive Phasor Estimation Algorithm Based on a Least Squares Method. Energies 2019, 12, 1387. https://doi.org/10.3390/en12071387

AMA Style

Kim W-J, Nam S-R, Kang S-H. Adaptive Phasor Estimation Algorithm Based on a Least Squares Method. Energies. 2019; 12(7):1387. https://doi.org/10.3390/en12071387

Chicago/Turabian Style

Kim, Woo-Joong, Soon-Ryul Nam, and Sang-Hee Kang. 2019. "Adaptive Phasor Estimation Algorithm Based on a Least Squares Method" Energies 12, no. 7: 1387. https://doi.org/10.3390/en12071387

APA Style

Kim, W. -J., Nam, S. -R., & Kang, S. -H. (2019). Adaptive Phasor Estimation Algorithm Based on a Least Squares Method. Energies, 12(7), 1387. https://doi.org/10.3390/en12071387

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop