Detection of Faults in Power System Using Wavelet Transform and
Independent Component Analysis

1
1
2

Dept. of Electrical and Electronics Engineering, IIIT, Bhubaneswar, India.
Dept. of Electrical Engineering, SOA University, Bhubaneswar, India.

3
3
4

Prakash K. Ray, 2B. K. Panigrahi, 2P. K. Rout

Asit Mohanty, 4Harishchandra Dubey
Dept. of Electrical Engineering, CET, Bhubaneswar, India.
Dept. of ECE, The University of Texas at Dallas, USA.

ABSTRACT: Uninterruptible power supply is the main motive of power utility companies that motivate them
for identifying and locating the different types of faults as quickly as possible to protect the power system
prevent complete power black outs using intelligent techniques. Thus, the present research work presents a
novel method for detection of fault disturbances based on Wavelet Transform (WT) and Independent
Component Analysis (ICA). The voltage signal is taken offline under fault conditions and is being processed
through wavelet and ICA for detection. The time-frequency resolution from WT transform detects the fault
initiation instant in the signal. Again, a performance index is calculated from independent component analysis
under fault condition which is used to detect the fault disturbance in the voltage signal. The proposed
approach is tested to be robust enough under various operating scenarios like without noise, with 20-dB noise
and variation in frequency. Further, the detection study is carried out using a performance index, energy
content, by applying the existing Fourier transform (FT), short time Fourier transform (STFT) and the
proposed wavelet transform. Fault disturbances are detected if the energy calculated in each scenario is
greater than the corresponding threshold value. The fault detection study is simulated in MATLAB/Simulink
for a typical power system.
1 INTRODUCTION
Modern power utilities require a efficient
protection scheme to protect the system itself as
well as the connected equipments for improving
better performance under normal as well as
abnormal/faulty operating scenarios. Now-a-days
the electromechanical relays are replaced by digital
relays because of their characteristics like faster
operation, accuracy and reliability. The fault
detection through digital relays and fault detector
(FD) is very vital for implementing any real-time
solutions (Phadke 1988, Dash 2000).
In the literature, so may work are for
monitoring and identification of faults in power
system. Fault situation can be identified by
comparing the difference between the value
between two consecutive cycles when higher than a
threshold value based on Phasor (Sidhu et al. 2002,
Sachdev et al., 1991). But, its demerit being the
modeling of fault resistance. Kalman filter is used
for fault detection based on estimation method
(Chowdhury et al. 1991, Zadeh et al. 2010, &
Girgis 1982). Wavelet transform (WT) is being

used for the detection of fault disturbances in power
system (Ukil & R. Živanović 2007) which detects
the changes in the signal. Then, adaptive filters and
wavelet transform in combination was used for fault
identification in power system (Ray et al. 2010).
But, these techniques affected by variation in
frequency, noise etc.
This paper proposes an algorithm for fault
detection using wavelet transform and independent
component analysis. The sudden changes are to be
detected by the detection techniques so that suitable
solutions can be adapted to protect the power
system from any type of disturbances. The proposed
method performances are being analyzed under
different operating scenarios like in presence of
noise, harmonics and frequency variations. Wavelet
transform (WT) and independent component
analysis (ICA) are suitable candidates for detection
of fault disturbances because of efficient timefrequency resolutions and their reliability to extract
the suitable features for identification purpose.
Again, the techniques are considered under different
operating conditions in order to assess the
robustness and accuracies (Pradhan et al. 2005).

This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by the authors or by the respective
copyright holders. The original citation of this paper is:
P. K. Ray, B. K. Panigrahi, P. K. Rout, A. Mohanty, H. Dubey, "Detection of Faults in Power System Using Wavelet Transform and
Independent Component Analysis", First International Conference on Advancement of Computer Communication & Electrical Technology, October 2016,
Murshidabad, India, DOI: 10.13140/RG.2.2.20394.82882

The techniques used for detection analysis is
presented in Section 2 followed by the algorithm of
implementation in Section 3. Then, Section 4
explains the result analysis both qualitatively and
quantitatively, followed by the conclusions in
Section 5.
2 TECHNIQUES FOR FAULT DETECTION
This section describes the techniques used for
identification of different types of faults in wind
system connected to grid power system. The WT
and ICA as detection techniques are presented
briefly with their mathematical modeling. Different
operating scenarios are taken as case studies to test
these methods for fault identifications which help in
assessing both normal as well as faulty operating
conditions. The details of these methods are as
follows.
2.1 Wavelet transforms (WT)

Figure 1: Wavelet decomposition tree

Considering a voltage signal of the power system
as ( ), the continuous wavelet transform (CWT) is
expressed as:


tN
V (t )* 
a 
 M

1



DWT (V , M , N ) 

1
M 0m


k

 n  kM 0m
V (k )* 
 Mm
0






(2)
Where M and N are replaced by M and kM ,
and k, mare integers. Scaling function in one stage
is expressed as sum of that of next stage and can be
given by:


 h ( n)

 (t ) 

2  (2t  n)

n 

(3)
Using the above equation, the original voltage
signal can be written as [8]:


V (t ) 



a

jo

(k )2 jo / 2  (2 jo t  k ) 



 d

j

(k )2 j / 2  (2 j t  k )

(4)
Where 0 is the coarse adjustment parameter in
the scaling function. The detail and approximate coefficients are written:
k 

k  j  jo



The wavelet transform is a signal processing
algorithm which is useful in detection of abnormal
operating conditions based on decomposition of the
power signals into different ranges of frequencies
by the help of a series of low-pass and high-pass
filters. This usually provides us a time-frequency
multi-resolution analysis that greatly useful for
identifying any short of abrupt variations in the
electrical parameters such as voltage, phase,
current, frequency etc. Here, Daubechies4 (dB4) is
being used as the mother wavelet basis function for
the fault detection analysis (Ukil & R. Živanović
2007). Usually, the signal is divided into a set of
approximate (a) and detail (d) co-efficient
representing the low-frequency and high frequency
bands respectively. The decomposition is as
presented below in Figure 1.

CWT (V , M , N ) 

Where
and
are called as dilation and
translation parameters and
is known as the
wavelet basis function. Now WT in discrete form
as:


 dt


(1)

a j (k ) 

a

j 1

(m)h(m  2k )
(5)

m 


d j (k ) 

c

j 1

(m)h1 (m  2k )

m 

(6)

2.2 Independent Component Analysis (ICA)
Independent Component Analysis (ICA) is an
advanced or modified version of principal
component analysis (PCA). It uses high order
statistical analysis based de-correlation for the
source input signal and provides important
information regarding the irregularities presence. In
ICA, the input data is modeled as linear coefficients
that are mutually independent to each other. Based
on the blind source separation algorithm, it can
efficiently and accurately transform the input
voltage signal into mutually and statistically
independent components, thereby helps the
detection process (Hyvärinen et al. 2001, Pöyhönen
et al. 2003, Lopez et al. 2011, & Dubey et al. 2011).
The mutual information nothing but the
individual independency measurement of the input
signal and its entropy is Gaussian in nature and can
be written as:
J ( y )  h(vgaussian )  h(v)
(7)
Differential entropy H of an input signal y with
density ( ) is given by:
h(v)    ph   log pv   d
(8)
Here, the negentropy is estimated based on the
estimation of probability functions of input signals
Then, rather the negentropy can be expressed as:

    e g  w x   e g  v

J  vi   J e wiT x

T
i

gaussian



2

(9)

Here, e is the statistical expectation and G is a nonquadratic factor. For n linear mixtures taken as
x , x , … . x for the n independent components, we
can express :
x j  a1s1  a2 s2  ....  an sn For all
(10)
Here, x is a random mixture and s is a random
independent variable. Assuming x as the random
vector with elements x , x , … . x and random
vector with elements s , s , … . s , if A is the matrix
having elements a , we can take column vectors;
x ,as the transpose of x. Then, we can write in
matrix form:
x  As
(11)
Where, A is a column matrix of elements a j and
then model is expressed as:
n

x   ai si

(12)

i 1

The ICA mathematical model is an iterative one
where the independent components are called latent
variables. Here the input mix matrix is considered
to be unknown and the components s are
statistically independent. Finally, once the matrix A
is estimated, we can determine the inverse, W, then
can evaluate the independent component by the
following expression:

s  Wx

(iv) Then, the principal components (PC) of the
input data are determined.
(v) In the current study, the independent
components are calculated based on fixed point
iteration [12]-[15].
(vi) Then, de-correlating
matrix, W and the
independent component, S of voltage signal is
determined from which a matrix, x is generated.
(vii) And finally, the pre-fault signal is considered
for constructing the matrix, x which can be used to
calculate the performance index and fault is
detected when the index is more than a set threshold
value. The performance index is given by
PI (k)  (normal(absolute(Wf (k )* xn (k )  s f (k )))2 ) (14)

4 SIMULATED RESULTS
This section explains the detection results and
corresponding discussions in a wind system
connected to grid. The power system is simulated in
MATLAB/Simulink environment and the grid is of
230 kV, 50 Hz rating consisting of thermal based
power plant. The considered grid-connected power
system is shown in Figure 2 and here, the sampling
frequency is taken as 2 kHz.

(13)
Figure 2: Grid-connected wind power System.

3 FAULT DETECTION METHODOLOGY
The fault detection methodology is explained in
this section for the grid-connected wind power
system as follows:
(i) A grid-connected wind power system is
simulated in MATLAB where different types of
fault are created. The voltage signal is extracted at
PCC under faulty conditions. The signal is then
passed through wavelet transform to get the timefrequency analysis for identification of faulty
situations.
(ii) The voltage signal is processed through to find
the mean and de-correlation. This is the 1st level of
processing. .
(iii) Then, X data size is reduced and filtered for
better data redundancy.

4.1 Fault detection using wavelet transform
The voltage signal at point of common coupling
(PCC) is taken offline under AG fault. The signal is
processed through WT which detects the fault
instant based on filtering through a set of low-pass
and high-pass filters. The voltage signal with its
detection result is shown in Figure 3. Similarly, the
detection result using AB and ABCG fault is shown
in Figure 4 and Figure 5 respectively. These results
clearly show that WT nicely detect the fault
initiation instant nicely. But, as noticed when the
fault is cleared, as shown in Figs., the recovery to

normal waveform is un-noticed

by wavelet

transform.

Figure 3: Detection of phase-to-ground fault in wind system
connected to grid (red- voltage signal; green-WT output)

4.3 Phase-to-phase fault detection using
Independent Component Analysis
Next, a phase-to-phase fault is near to the point of
common coupling and the corresponding voltage
signal at the PCC. Then, the index as explained in
ICA section is calculated under different operating
conditions like without noise, with 20-dB noise and
under variation in frequency condition. The
simulated of the performance index obtained from
ICA is shown in Figure 6 (b). It is observed that, the
fault initiated at near about 0.065 s detected through
sudden increase in the index value. Before the fault,
the index value is observed to be almost zero. Once,
the fault occurs, the value increases suddenly and if
we set a threshold value, then the can be detected.
But in this case, we have not set the threshold value.
Rather, we are assessing the faulty condition based
on sudden increase.

Figure 4: Detection of phase-to-phase fault in wind system
connected to grid (red- voltage signal; green-WT output)

(a)

Figure 5: Detection of three-phase-to-ground fault in wind
system connected to grid (red- voltage signal; green-WT
output)

4.2 Phase-to-ground fault detection using
Independent Component Analysis
A phase-to-ground fault is near to the point of
common coupling and the corresponding voltage
signal at the PCC is taken offline. Then, the index
as explained in ICA section is calculated under
different operating conditions like without noise,
with 20-dB noise and under variation in frequency
condition. The simulated of the performance index
obtained from ICA is shown in Figure 6 (a). It is
observed that, the fault initiated at near about 0.065
s detected through sudden increase in the index
value. Before the fault, the index value is observed
to be almost zero. Once, the fault occurs, the value
increases suddenly and if we set a threshold value,
then the can be detected. But in this case, we have
not set the threshold value. Rather, we are assessing
the faulty condition based on sudden increase.

(b)
Figure 6: Detection of (a) phase-to-ground fault (b) phase-tophase fault in wind system connected to grid

4.4 Detection of fault disturbances using the
performance index (PI)
This sub-section describes the detection
technique using a performance index called energy
content of the processed faulty signal. The voltage
signal at the point of common coupling is extracted
and processed through different signal processing
transforms like Fourier transform (FT), Short time
Fourier transform (STFT) and wavelet transform.
Then the output waveform after processing is used
to calculate energy which is compared with a
threshold value to know the faulty or normal
operating conditions in the power system. It is

observed that proposed WT provides larger energy
value so that a suitable threshold can be chosen to
accurately detect the faulty conditions. The index in
different cases is shown in Table 1.
Table 1 Energy content using FT, STFT,WT

Zadeh, R.A., A. Ghosh, & G. Ledwich (2010). Combination
of Kalman Filter and Least-Error Square Techniques
IEEE Trans. on Power System 25 (4) ,2868-2880.
Girgis A. (1982). A new Kalman ﬁltering based digital
distance relaying IEEE Trans. on Power Apparatus and
Systems 101,3471–3480.

Scenario
AG Fault
BG Fault
CG Fault
AB Fault
BC Fault
ABC Fault

FT
1.3672
1.4525
1.3324
2.2341
2.6342
3.1302

Energy content
STFT
1.7863
2.2414
1.8367
2.3532
3.1230
3.8225

WT
2.1663
2.7352
2.6538
3.2514
3.8724
4.2431

Ukil, A. & R. Živanović (2007). Abrupt change detection in
power system fault analysis using wavelet transform.
International Conference on Power Systems Transients,
Montreal, Canada.June. 19-23.
Ukil, A. and R. Živanović, (2007). Application of Abrupt
Change Detection in Power Systems Disturbance
Analysis and Relay Performance Monitoring.

5 CONCLUSIONS
This paper presented detection of fault
disturbances using Wavelet Transform and
Independent Component Analysis. Voltage signal
extracted at PCC is being processed through the
above techniques under different fault and operating
scenarios. It is observed that WT transform detected
the fault instant very accurately. But, in some cases,
it could not discriminate the disturbance condition
may be because of noise or frequency variation.
Therefore, the faults are accurately detected using
ICA under all operating scenarios

IEEE

Trans. on Power Delivery 22 (1),59-66.
Ray, P. K., H. C. Dubey, S. R.Mohanty, N. Kishor, & K.
Ganesh (2010). Power quality disturbance detection in
grid-connected wind energy system using wavelet and Stransform, IEEE ICPCES,1-4.
Pradhan, A. K. , A.Routray, & S. R.Mohanty (2005). A
Moving Sum Approach for Fault Detection of Power
Systems Electric Power Components and Systems 34 (4),
385 – 399.
Hyvärinen, A., J. Karhunen, & E. Oja (2001). Independent
Component Analysis. A Wiley Interscience Publication,

REFERENCES
Phadke, G. & J. S. Thorp (1988). Computer Relaying for Power
Systems, New York, John Wiley.
Dash, P. K., A. K. Pradhan, & G. Panda (2000). A novel fuzzy
neural network based distance relaying scheme. IEEE
Trans. on Power Delivery. 15, 902–907.

John Wiley & Sons, Inc.
Sanna, Pöyhönen.,

P. Jover, & H. Hyötyniemi (2003).

Independent component analysis of vibrations for fault
diagnosis of an induction motor Proceedings of IASTED
International

Conference

Circuits,

Signals,

and

Systems,19-21.
Lopez, M. G., H. M.Lozano, L. P. Sanchez, &L. N. O.

Sidhu, T. S., D. S. Ghotra, & M. S. Sachdev (2002). An

Moreno (2011). Blind Source Separation of Audio

adaptive distance relay and its performance comparison

Signals Using Independent Component Analysis and

with a ﬁxed data window distance relay. IEEE Trans. on

Wavelets

Power Delivery. 17,691–697.

Communications and Computers (CONIELECOMP)

Sachdev, M. S. & M. A.Nagpal (1991). Recursive least square

21st International Conference on Electrical

Mexico, 152-157.

algorithm for power system relaying and measurement

Dubey, H.C., S. R. Mohanty, & N. Kishore (2011). Abrupt

applications. IEEE Trans. on Power Delivery. 6,1008–

change detection of fault in power system using

1015.

independent

Chowdhury, F. N., J. P. Christensen, & J. L.Aravena (1991).

component analysis

Conference on Signal Processing, Communication,

Power system fault detection and state estimation using

Computing

Kalman ﬁlter with hypothesis testing. IEEE Trans. on

(ICSCCN),659 – 664.

Power Delivery. 6, 1025–1030.

International

and

Networking

Technologies

