DSP Based Online Power Quality Events Detection and Classification Using Hilbert Huang Transform and Random Forest Method
DSP Based Online Power Quality Events Detection and Classification Using Hilbert Huang Transform and Random Forest Method
DSP Based Online Power Quality Events Detection and Classification Using Hilbert Huang Transform and Random Forest Method
Abstract—In this paper, empirical mode decomposition (EMD), spike, and DC offset which makes sustenance of power quality
Hilbert transform (HT), and random forest technique are com- a precarious issue. Disturbances in the power system are
bined to detect and classify the power quality events (PQEs) classified basically into stationary and non-stationary where
in real-time. Empirical mode decomposition (EMD) is used to
decompose the non-stationary PQEs into the mono-component the non-stationary disturbances are enormously present in
mode of oscillations, known as intrinsic mode functions (IMF). the power system network. These disturbances lead to the
The most efficacious features are extracted by applying Hilbert disintegration of the execution of transmission and distribution
transform on individual IMF. Four features such as standard system components, so indeed diminutive unsettling influences
deviation of magnitude, Hilbert energy, Tsallis entropy, and in the original signal needs to be checked so as to take
standard deviation of phase are computed from the AM-FM
signal of Hilbert transform to train the random forest classifier. appropriate preventive activity to avoid power quality issues.
Random forest (RF) classifier, a congregated form of numerous Primarily preprocessing, optimal feature extraction and clas-
self-determining decision trees that has a unique ability to deal sification are followed by recognition and detection of power
well with uneven data sets with missing variables and faster system problems. The time and frequency domain information
learning speed compared with other machine learning methods. from PQ disturbances are extracted by utilizing the meth-
The main objective of this paper is to assess the worth of
the proposed HHT-RF method implemented in the hardware ods of Fourier transform (FT) and wavelet transform (WT)
prototype using a digital signal processor (DSP) in the power through the process of preprocessing and feature extraction
quality events detection and classification in real-time. [1-2]. Fourier analysis is deficient of time information during
Index Terms—Empirical mode decomposition (EMD), random transforming to frequency domain so it does not have sufficient
forest (RF), Hilbert transform (HT), intrinsic mode function time information about the event on the other hand also
(IMF), non-stationary power quality events, real-time analysis,
Digital signal processor (DSP). bounded by stationary signals only. Most PQEs issues are
non-stationary in nature so for that there is a requirement
I. I NTRODUCTION of method like a short-time Fourier transform (STFT) and
wavelet transforms that would give frequency components
Online Power quality (PQ) signal disturbance’s recognition included with the timing of event of such disturbances [3-
and classification is the foundation of PQ problem manage- 7]. As non-stationary signals are the integration of stationary
ment in this rapidly growing power system with the multiple signals hence STFT [8] decomposes it into time-frequency
uses of nonlinear loads, fast-acting power electronic devices, domains within the sliding window which provide the timing
computers, lighting controls, and the unavoidable unbalanced of the event of such disturbances. STFT cannot adequately
power system. Quality of power is remarkably weakened analyze transient phenomena whereas, the WT with multi-
due to the event of unsettling influences such as transient, resolution capability proves to be an efficient analyzer of
harmonics, flicker, sag, momentary interruption, swell, notch, nonlinear and transient phenomena by disintegrating the signal
978-1-7281-7274-3/20/$31.00
c 2020 IEEE into multiple time-frequency levels and, therefore, maintaining
Authorized licensed use limited to: NATIONAL INSTITUTE OF TECHNOLOGY HAMIRPUR. Downloaded on August 28,2023 at 12:08:00 UTC from IEEE Xplore. Restrictions apply.
IEEE International Conference on Computational Intelligence for Smart Power System and Sustainable Energy (CISPSSE-2020), July 29-31, 2020, Odisha, India
T T
II. H ILBERT H UANG TRANSFORM 2
2
SD(n) = | gn−1 (t) − gn (t) | / | gn−1 (t) | < 0.05
The advanced signal processing algorithm namely empirical t=0 t=0
mode decomposition is the primary algorithm of Hilbert
Huang transform (HHT) to extract the intrinsic mode functions (1)
(IMFs) from the original non-linear and non-stationary power
quality events (PQEs) signals. The Hilbert transform (HT) where n illustrates the nth difference between the original
is applied to each and every IMFs individually to extract signal r(t) and the envelope mean v(t). Otherwise apply the
competent features of the PQEs signals. step-2 to extract the IM F2 , IM F3 , IM F4 unless until an
IM F = 0. The PQE decomposed signal EMD is shown in
A. Empirical Mode Decomposition Fig. 1. The original PQE signal r(t) is regenerated using the
equation given below.
The empirical mode decomposition (EMD) is applied on m
the PQE signal to decomposed into a number of mono-
r(t) = gi (t) + km (t) (2)
component modes named as intrinsic mode functions (IMFs).
i=0
Authorized licensed use limited to: NATIONAL INSTITUTE OF TECHNOLOGY HAMIRPUR. Downloaded on August 28,2023 at 12:08:00 UTC from IEEE Xplore. Restrictions apply.
IEEE International Conference on Computational Intelligence for Smart Power System and Sustainable Energy (CISPSSE-2020), July 29-31, 2020, Odisha, India
where m represents the IMF number and they are orthogonal The performance of the proposed method is tested using
to each other. fourteen PQEs of synthetic, real, and synthetic with noise
signals presented in Table I. Fig. 2 shows the magnitude
response, phase response, frequency response, and energy
B. Hilbert transform
spectrum of PQE signal using HT for the detection of multiple
The Hilbert transform is applied on the imaginary part of synthetic PQEs.
the extracted real-valued IMF signal into analytic form for
the purpose of feature extraction [14]. The complex-valued
analytic signal representation of IMFs is defined as,
Authorized licensed use limited to: NATIONAL INSTITUTE OF TECHNOLOGY HAMIRPUR. Downloaded on August 28,2023 at 12:08:00 UTC from IEEE Xplore. Restrictions apply.
IEEE International Conference on Computational Intelligence for Smart Power System and Sustainable Energy (CISPSSE-2020), July 29-31, 2020, Odisha, India
Authorized licensed use limited to: NATIONAL INSTITUTE OF TECHNOLOGY HAMIRPUR. Downloaded on August 28,2023 at 12:08:00 UTC from IEEE Xplore. Restrictions apply.
IEEE International Conference on Computational Intelligence for Smart Power System and Sustainable Energy (CISPSSE-2020), July 29-31, 2020, Odisha, India
Confusion Matrix
98 0 0 0 0 0 0 0 0 0 0 0 0 0 100%
1
7.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 100 0 0 0 0 0 0 0 0 0 0 0 0 100%
2
0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 0 100 0 0 0 0 0 0 0 0 0 0 0 100%
3
0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 0 0 100 0 0 0 0 0 0 0 0 0 0 100%
4
0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 0 0 0 100 0 0 0 0 0 0 0 0 0 100%
5
0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 0 0 0 0 100 0 0 0 0 0 0 0 0 100%
6
0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
PQE Output Class
0 0 0 0 0 0 100 0 0 0 0 0 0 0 100%
7
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 0 0 0 0 0 0 98 0 0 0 0 0 0 100%
8
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
0 0 0 0 0 0 0 0 100 0 0 0 0 0 100%
9
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0%
2 0 0 0 0 0 0 2 0 100 0 0 0 0 96.2%
10
0.1% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.1% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0% 3.8%
0 0 0 0 0 0 0 0 0 0 100 0 0 0 100%
11
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0% 0.0%
0 0 0 0 0 0 0 0 0 0 0 100 0 0 100%
12
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0% 0.0%
0 0 0 0 0 0 0 0 0 0 0 0 100 0 100%
13
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0% 0.0%
0 0 0 0 0 0 0 0 0 0 0 0 0 100 100%
14
0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 7.1% 0.0%
98.0% 100% 100% 100% 100% 100% 100% 98.0% 100% 100% 100% 100% 100% 100% 99.7%
2.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 2.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.0% 0.3%
10
11
12
13
14
1
Target Class
method for the recognition of the single as well as combined ing the most suitable twelve features in both the noise-free
PQE patterns in both the noisy and noise-free environments. and noisy environments. Both the simulation and laboratory
The proposed HHT-RF method can be implemented for the experiment prove the remarkable recognition performance of
superior recognition of real-time PQE patterns to avoid power the proposed HHT-RF method for both the synthetic and
system pollution. practical PQEs. The outstanding classification capability, short
event detection time, simplicity and practicability of HHT-RF
VI. C ONCLUSION method as compared to the existing prevalent method prove
The less computational complex empirical mode decom- that it is the utmost choice to monitoring the single as well as
position is used to extract most suitable IMFs for the com- combined PQE patterns in real-time.
putation of PQIs. Four efficacious PQIs are extracted from
the Hilbert transformed array for better classification accuracy. R EFERENCES
The magnitude response, frequency response, phase response
and Hilbert energy spectrum show the quick, clear and accu- [1] Bracewell, Ronald Newbold, and Ronald N. Bracewell, “The Fourier
rate detection capability for the combined PQE patterns. The transform and its applications,” vol. 31, no. 99, New York: McGraw-
Hill, 1986.
high-speed T M S320C6713 digital signal processor starter
[2] Wang, Zhongde. “Fast algorithms for the discrete Wavelet transform and
kit is more suitable to embed the proposed HHT-RF method for the discrete Fourier transform,” IEEE Transactions on Acoustics,
in MATLAB/Simulink with code composer studio software Speech, and Signal Processing, vol. 32, no. 4, pp. 803–816, 1984.
environments for the quick recognition of the complex PQE [3] Chen, J., Rostami, J., Peter, W.T. and Wan, X.,“The design of a novel
mother Wavelet that is tailor-made for continuous Wavelet transform in
patterns in real-time. The ensemble learning based RF clas- extracting defect-related features from reflected guided wave signals,”
sifier produces outstanding classification accuracy by import- vol. 110, pp. 176–191, 2017.
Authorized licensed use limited to: NATIONAL INSTITUTE OF TECHNOLOGY HAMIRPUR. Downloaded on August 28,2023 at 12:08:00 UTC from IEEE Xplore. Restrictions apply.
IEEE International Conference on Computational Intelligence for Smart Power System and Sustainable Energy (CISPSSE-2020), July 29-31, 2020, Odisha, India
Authorized licensed use limited to: NATIONAL INSTITUTE OF TECHNOLOGY HAMIRPUR. Downloaded on August 28,2023 at 12:08:00 UTC from IEEE Xplore. Restrictions apply.