# Stepped Fault Line Selection Method Based on Spectral Kurtosis and Relative Energy Entropy of Small Current to Ground System

Authors: Xiaowei Wang; Xiangxiang Wei; Jie Gao; Yaxiao Hou; Yanfang Wei
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

This paper proposes a stepped selection method based on spectral kurtosis relative energy entropy. Firstly, the length and type of window function are set; then when fault occurs, enter step 1: the polarity of first half-wave extremes is analyzed; if the ratios of extremes between neighboring lines are positive, the bus bar is the fault line, else, the SK relative energy entropies are calculated, and then enter step 2: if the obtained entropy multiple is bigger than the threshold or equal to the threshold, the overhead line of max entropy corresponding is the fault line, if not, enter step 3: the line of max entropy corresponding is the fault line. At last, the applicability of the proposed algorithm is presented, and the comparison results are discussed.


---

## Body

## 1. Introduction

Most of 3~66 kV distribution networks are so-called small current to ground system in china, which includes unearthed neutral system, arc suppression coil compensated neutral system, and high resistance-grounded neutral system. Since the fault current is small and the arc is unstable instability, suitable fault line selection methods are lacking [1, 2].When fault occurs, the researchers are attracted by the rich transient fault features. At present, the main research methods are based on the transient fault signal, such as Wavelet transform [3, 4], S transform [5, 6], Prony algorithm [7], Hilbert-Huang transform (HHT) [8, 9], and correlation analysis [10, 11]. In [3, 4], the transient zero-sequence current (TZSC) is decomposed by wavelet, and the transient energy is calculated; according to the differences, the fault line is selected, when the high grounding resistance fault occurs, which can lead to misjudgment because of the compensation of arc suppression coil to the TZSC. Good time frequency of S transform can be used to select fault line; in [5, 6], the dominant frequencies of TZSC are obtained by S transform, which are used as the selection criterions; however, there is much information after S transform, and what should be further studied is how to take advantage of the phase angles. Prony algorithm can fit the low frequency transient signal well; in [7], Prony is used to fit the TZSC to select the dominant components, and then the relative entropy of dominant components is obtained to select the fault line, but what should be further studied is how to determine the model orders. In [8, 9], HHT is used to decompose the TZSC, and then the most high-frequency component of the intrinsic mode functions (IMF) can be obtained, and based on this, the selection criterion is built; however, the decomposition process may cause modal aliasing. In [10, 11], the poor waveform similarity of the TZSC between the fault line and the healthy line is noticed, and the coefficients are obtained to select the fault line; however, when the small-angle fault occurs, the similarity of the TZSC that flows the healthy cable line and the healthy overhead lines gets worse.Aiming at few features that are available in arc suppression coil to ground system, single fault line selection method is not very reliable, and this paper regards the ratios of the first half-wave extremes between neighboring lines and spectral kurtosis (SK) relative energy entropy as features to select the fault line, a stepped fault line selection method is proposed by utilizing two kinds of features. When fault occurs, enter step 1, and if the ratios of first half-wave extremes between neighboring lines are positive, the bus bar is selected as the fault line, and the line selection process is done. Otherwise, the SK relative energy entropies of branch lines are computed, then enter step 2, if the obtained entropy multiple is greater than or equal to the threshold, the overhead line is judged as the fault line, and the overhead line of max entropy corresponding is selected as the fault line, and the fault line selection process is done. If the entropy multiple is less than the threshold, enter step  3, and the line of max entropy corresponding is selected as the fault line.The remaining of this paper is organized as follows. In Section2, characteristics of TZSC are presented. In Section 3, basic theories about spectral kurtosis and wavelet denoising are presented. In Section 4, feasibility analysis of spectral kurtosis algorithm is presented. In Section 5, basic theories about fault line selection are presented. In Section 6, fault line selection criterions are proposed. In Section 7, simulation and verification are presented. In Section 8, applicability analysis is presented. In Section 9, two kinds of comparison results are discussed. In Section 10, the paper is completed with conclusions and future directions.


## 2. Characteristics of TZSC

Zero sequence equivalent circuit of single phase to ground is shown in Figure1, C
0 and L
0 are the zero-sequence equivalent capacitance and inductance, respectively, R
g is the ground transition resistance, R
p and L
p are the equivalent resistance and inductance of arc suppression coil, respectively, and G
(
t
) is the zero-sequence voltage.Figure 1
Zero sequence equivalent circuit.From Figure1, the TZSC i
0
(
t
) through the fault point is
(1)
i
0
(
t
)
=
i
0
L
(
t
)
+
i
0
C
(
t
)
=
I
L
m
cos
⁡
φ
e
-
t
/
τ
L
+
I
C
m
(
ω
f
ω
sin
φ
sin
ω
t
-
cos
⁡
φ
cos
⁡
ω
f
t
)
e
-
δ
t
,i
0
L
(
t
) and i
0
C
(
t
) are the inductor and capacitive current, respectively, I
L
m and I
C
m are the initial value of i
0
L
(
t
) and i
0
C
(
t
), respectively, ω is the angular frequency, ω
f and δ are the angular frequency of free oscillations and the attenuation coefficient, respectively, and δ = 1/τ
c, τ
c is the decay time constant of i
0
C
(
t
), τ
L is the decay time constant of i
0
L
(
t
), and φ is the initial phase of the fault line [12].From (1), i
0
L
(
t
) and i
0
C
(
t
) make up the fault TZSC, both of which include nonperiodic damped component and stable periodic component. At beginning, τ
L is greater than τ
c; therefore, i
0
L
(
t
) is deteriorated slower than i
0
C
(
t
).It can be known from the above analysis that the TZSC is mutant signal when fault occurs, and as time goes on, it would decrease gradually. Considering the good frequency-domain effect for unstable, nonperiodic, mutant, and damped signal of SK, SK algorithm is introduced to extract the feature information of TZSC in this paper.


## 3. Basic Theories

### 3.1. Definition of SK

Kurtosis is random 4-order cumulative amount and cannot reflect specific signal changes as a global index; therefore, kurtosis is improper to detect the signal with strong noise. In order to overcome the weakness of kurtosis, SK is proposed by Immovilli et al. [13–15].x
(
t
) is decomposed by Wold-Cramer under the nonstationary condition, and the result is as follows:
(2)
Y
(
t
)
=
∫
-
∞
+
∞
e
j
2
π
f
t
H
(
t
,
f
)
d
X
(
f
)
,H
(
t
,
f
) is the complex envelope at frequency f of Y
(
t
), and also H
(
t
,
f
) changes with time. In a practical system, H
(
t
,
f
) is random, X
(
f
) is the strict white noise spectrum, and ω is random with time varying of filter.Fourth-order spectrum cumulant ofY
(
t
) can be defined as follows:
(3)
C
4
Y
(
f
)
=
S
4
Y
(
f
)
-
S
2
2
Y
(
f
)
(
f
≠
0
)
,S
2
n
Y
(
f
) is the 2
n moments, which can be defined as follows:
(4)
S
2
n
Y
(
f
)
≜
E
{
|
H
(
t
,
f
)
d
X
(
f
)
|
2
n
}
d
f
.Therefore SK can be defined as follows:(5)
K
(
f
)
≜
C
4
Y
(
f
)
S
2
Y
2
(
f
)
=
S
4
Y
(
f
)
S
2
Y
2
(
f
)
-
2
,
(
f
≠
0
)
.The important properties of SK are as follows.(1)
The SK of pure stationary process is −1,f
≠
0.
(2)
The SK of a stationary Gaussian process is 0.The TZSCi
0
(
t
) can be expressed as follows:
(6)
i
0
(
t
)
=
Y
(
t
)
+
b
(
t
)
,Y
(
t
) is the actual TZSC and b
(
t
) is the noise. Y
(
t
) and b
(
t
) are independent; therefore, the SK of i
0
(
t
) is(7)
K
(
x
+
b
)
(
f
)
=
K
x
(
f
)
[
1
+
ρ
(
f
)
]
2
,ρ
(
f
)
=
(
S
b
(
f
)
/
S
i
0
(
f
)
) is the ratio of noise to signal, and S
b
(
f
) and S
i
0
(
f
) are power spectral densities of b
(
t
) and Y
(
t
), respectively.It can be known from (7) that K
(
x
+
b
)
(
f
)
≈
K
i
0
(
f
) when ρ
(
f
) is very small and K
(
x
+
b
)
(
f
)
≈
0 when ρ
(
f
) is very large. So the SK can search the whole frequency band, and the characteristic frequency bands can be detected easily.


### 3.2. SK Based on Short-Time Fourier Transforms (STFT)

i
0
(
n
) is the discrete form of i
0
(
t
),  n
=
1,2
,
…
,
N, and the calculation of SK based on STFT is as follows:
(8)
Y
w
(
k
P
,
f
)
=
∑
n
=
1
N
i
0
(
n
)
w
(
n
-
k
P
)
e
-
j
2
π
n
f
,w
(
n
) is the window function, whose length is N
w, P is the time step, and k is the number of time steps.The2
n-order spectrum distance of  Y
w
(
k
P
,
f
) is defined as follows:
(9)
S
2
n
Y
=
〈
|
Y
w
(
k
P
,
f
)
|
2
n
〉
k
,〈
·
〉
k is k-order time average.Considering the definition of SK in Section3.1, when n
=
1 and n
=
2, based on STFT, and SK can be described as follows [16]:
(10)
K
(
f
)
=
S
4
Y
(
f
)
S
2
Y
2
(
f
)
-
2
.It can be seen that thei
0
(
t
) can be characterized by SK based on STFT, and the frequency band of the max SK corresponding can be obtained finally.


### 3.3. Basic Principle of Wavelet Denoising

#### 3.3.1. Basic Principle of Wavelet Transform

The basis function of the wavelet transform can be defined as follows [17, 18]:
(11)
ψ
u
,
v
(
t
)
=
1
u
ψ
(
t
-
v
u
)
u
>
0
,
v
∈
R
,
where u is the scale parameters and v is the translation parameter.The continuous wavelet transform ofψ
(
t
) is a function defined as
(12)
W
T
f
(
u
,
v
)
=
〈
f
,
ψ
u
,
v
〉
=
1
u
∫
R
f
(
t
)
ψ
¯
(
t
-
v
u
)
d
t
.In order for an inverse wavelet transform to exist, the mother waveletψ
(
t
) should satisfy the admissibility condition: C
ψ
=
∫
R
(
|
ψ
(
ω
^
)
|
2
/
|
ω
|
)
d
ω
<
+
∞
(
ψ
(
ω
^
) is the complex conjugate of ψ
(
t
)). So ψ
(
t
) is an admissible wavelet. Choosing admissible wavelet, the original signals are reconstructed by
(13)
f
(
t
)
=
1
C
ψ
∫
∫
R
W
T
f
(
u
,
v
)
·
[
1
u
ψ
(
t
-
v
u
)
]
1
u
2
d
u
d
v
.In practice,u and v must be discrete. Let u
=
u
0
m and v
=
n
v
0
u
0
m, and (13) can be presented as
(14)
ψ
m
,
n
(
t
)
=
u
0
-
m
/
2
ψ
(
u
0
-
m
t
-
n
v
0
)
.Function family{
ψ
m
,
n
(
t
)
,
m
,
n
∈
Z
} is called discrete wavelet. When u
0
=
2 and v
0
=
1, function family {
ψ
m
,
n
(
t
)
} constructs a group orthonormal basis on which signals can be reconstructed.


#### 3.3.2. Process of Wavelet Denoising

The basic principle of wavelet denoising is based on the spectrum distribution characteristics of signals and noise. According to the principle, noise of each decomposition level is eliminated, and useful signal is preserved [19, 20].The substance of denoising is restraining the unwanted signal and boosting the useful signal. The one-dimensionalK
i
(
i
=
1,2
,
…
,
N
w
) can be expressed as follows:
(15)
K
i
=
f
i
+
ε
e
i
,K
i is the observable SK with noise, f
i is the real SK without noise, e
i is the gauss random signal, and ε is the standard deviation of noise.The denoising process of one-dimensional signals can be divided into three steps.Step I.
Mother wavelet function and the levels of decomposition are determined, and thenK
i can be decomposed by wavelet transform.Step II.
High frequency coefficients are processed by the threshold.
The threshold thr is(16)
thr
=
σ
2
log
⁡
e
⁡
(
N
w
)
,
whereσ is the standard deviation of K
i and N
w is the length of K
i.Step III.
Signals can be reconstructed by adding the low-frequency coefficients and high frequency coefficients, and the result isK
i
′.


## 3.1. Definition of SK

Kurtosis is random 4-order cumulative amount and cannot reflect specific signal changes as a global index; therefore, kurtosis is improper to detect the signal with strong noise. In order to overcome the weakness of kurtosis, SK is proposed by Immovilli et al. [13–15].x
(
t
) is decomposed by Wold-Cramer under the nonstationary condition, and the result is as follows:
(2)
Y
(
t
)
=
∫
-
∞
+
∞
e
j
2
π
f
t
H
(
t
,
f
)
d
X
(
f
)
,H
(
t
,
f
) is the complex envelope at frequency f of Y
(
t
), and also H
(
t
,
f
) changes with time. In a practical system, H
(
t
,
f
) is random, X
(
f
) is the strict white noise spectrum, and ω is random with time varying of filter.Fourth-order spectrum cumulant ofY
(
t
) can be defined as follows:
(3)
C
4
Y
(
f
)
=
S
4
Y
(
f
)
-
S
2
2
Y
(
f
)
(
f
≠
0
)
,S
2
n
Y
(
f
) is the 2
n moments, which can be defined as follows:
(4)
S
2
n
Y
(
f
)
≜
E
{
|
H
(
t
,
f
)
d
X
(
f
)
|
2
n
}
d
f
.Therefore SK can be defined as follows:(5)
K
(
f
)
≜
C
4
Y
(
f
)
S
2
Y
2
(
f
)
=
S
4
Y
(
f
)
S
2
Y
2
(
f
)
-
2
,
(
f
≠
0
)
.The important properties of SK are as follows.(1)
The SK of pure stationary process is −1,f
≠
0.
(2)
The SK of a stationary Gaussian process is 0.The TZSCi
0
(
t
) can be expressed as follows:
(6)
i
0
(
t
)
=
Y
(
t
)
+
b
(
t
)
,Y
(
t
) is the actual TZSC and b
(
t
) is the noise. Y
(
t
) and b
(
t
) are independent; therefore, the SK of i
0
(
t
) is(7)
K
(
x
+
b
)
(
f
)
=
K
x
(
f
)
[
1
+
ρ
(
f
)
]
2
,ρ
(
f
)
=
(
S
b
(
f
)
/
S
i
0
(
f
)
) is the ratio of noise to signal, and S
b
(
f
) and S
i
0
(
f
) are power spectral densities of b
(
t
) and Y
(
t
), respectively.It can be known from (7) that K
(
x
+
b
)
(
f
)
≈
K
i
0
(
f
) when ρ
(
f
) is very small and K
(
x
+
b
)
(
f
)
≈
0 when ρ
(
f
) is very large. So the SK can search the whole frequency band, and the characteristic frequency bands can be detected easily.


## 3.2. SK Based on Short-Time Fourier Transforms (STFT)

i
0
(
n
) is the discrete form of i
0
(
t
),  n
=
1,2
,
…
,
N, and the calculation of SK based on STFT is as follows:
(8)
Y
w
(
k
P
,
f
)
=
∑
n
=
1
N
i
0
(
n
)
w
(
n
-
k
P
)
e
-
j
2
π
n
f
,w
(
n
) is the window function, whose length is N
w, P is the time step, and k is the number of time steps.The2
n-order spectrum distance of  Y
w
(
k
P
,
f
) is defined as follows:
(9)
S
2
n
Y
=
〈
|
Y
w
(
k
P
,
f
)
|
2
n
〉
k
,〈
·
〉
k is k-order time average.Considering the definition of SK in Section3.1, when n
=
1 and n
=
2, based on STFT, and SK can be described as follows [16]:
(10)
K
(
f
)
=
S
4
Y
(
f
)
S
2
Y
2
(
f
)
-
2
.It can be seen that thei
0
(
t
) can be characterized by SK based on STFT, and the frequency band of the max SK corresponding can be obtained finally.


## 3.3. Basic Principle of Wavelet Denoising

### 3.3.1. Basic Principle of Wavelet Transform

The basis function of the wavelet transform can be defined as follows [17, 18]:
(11)
ψ
u
,
v
(
t
)
=
1
u
ψ
(
t
-
v
u
)
u
>
0
,
v
∈
R
,
where u is the scale parameters and v is the translation parameter.The continuous wavelet transform ofψ
(
t
) is a function defined as
(12)
W
T
f
(
u
,
v
)
=
〈
f
,
ψ
u
,
v
〉
=
1
u
∫
R
f
(
t
)
ψ
¯
(
t
-
v
u
)
d
t
.In order for an inverse wavelet transform to exist, the mother waveletψ
(
t
) should satisfy the admissibility condition: C
ψ
=
∫
R
(
|
ψ
(
ω
^
)
|
2
/
|
ω
|
)
d
ω
<
+
∞
(
ψ
(
ω
^
) is the complex conjugate of ψ
(
t
)). So ψ
(
t
) is an admissible wavelet. Choosing admissible wavelet, the original signals are reconstructed by
(13)
f
(
t
)
=
1
C
ψ
∫
∫
R
W
T
f
(
u
,
v
)
·
[
1
u
ψ
(
t
-
v
u
)
]
1
u
2
d
u
d
v
.In practice,u and v must be discrete. Let u
=
u
0
m and v
=
n
v
0
u
0
m, and (13) can be presented as
(14)
ψ
m
,
n
(
t
)
=
u
0
-
m
/
2
ψ
(
u
0
-
m
t
-
n
v
0
)
.Function family{
ψ
m
,
n
(
t
)
,
m
,
n
∈
Z
} is called discrete wavelet. When u
0
=
2 and v
0
=
1, function family {
ψ
m
,
n
(
t
)
} constructs a group orthonormal basis on which signals can be reconstructed.


### 3.3.2. Process of Wavelet Denoising

The basic principle of wavelet denoising is based on the spectrum distribution characteristics of signals and noise. According to the principle, noise of each decomposition level is eliminated, and useful signal is preserved [19, 20].The substance of denoising is restraining the unwanted signal and boosting the useful signal. The one-dimensionalK
i
(
i
=
1,2
,
…
,
N
w
) can be expressed as follows:
(15)
K
i
=
f
i
+
ε
e
i
,K
i is the observable SK with noise, f
i is the real SK without noise, e
i is the gauss random signal, and ε is the standard deviation of noise.The denoising process of one-dimensional signals can be divided into three steps.Step I.
Mother wavelet function and the levels of decomposition are determined, and thenK
i can be decomposed by wavelet transform.Step II.
High frequency coefficients are processed by the threshold.
The threshold thr is(16)
thr
=
σ
2
log
⁡
e
⁡
(
N
w
)
,
whereσ is the standard deviation of K
i and N
w is the length of K
i.Step III.
Signals can be reconstructed by adding the low-frequency coefficients and high frequency coefficients, and the result isK
i
′.


## 3.3.1. Basic Principle of Wavelet Transform

The basis function of the wavelet transform can be defined as follows [17, 18]:
(11)
ψ
u
,
v
(
t
)
=
1
u
ψ
(
t
-
v
u
)
u
>
0
,
v
∈
R
,
where u is the scale parameters and v is the translation parameter.The continuous wavelet transform ofψ
(
t
) is a function defined as
(12)
W
T
f
(
u
,
v
)
=
〈
f
,
ψ
u
,
v
〉
=
1
u
∫
R
f
(
t
)
ψ
¯
(
t
-
v
u
)
d
t
.In order for an inverse wavelet transform to exist, the mother waveletψ
(
t
) should satisfy the admissibility condition: C
ψ
=
∫
R
(
|
ψ
(
ω
^
)
|
2
/
|
ω
|
)
d
ω
<
+
∞
(
ψ
(
ω
^
) is the complex conjugate of ψ
(
t
)). So ψ
(
t
) is an admissible wavelet. Choosing admissible wavelet, the original signals are reconstructed by
(13)
f
(
t
)
=
1
C
ψ
∫
∫
R
W
T
f
(
u
,
v
)
·
[
1
u
ψ
(
t
-
v
u
)
]
1
u
2
d
u
d
v
.In practice,u and v must be discrete. Let u
=
u
0
m and v
=
n
v
0
u
0
m, and (13) can be presented as
(14)
ψ
m
,
n
(
t
)
=
u
0
-
m
/
2
ψ
(
u
0
-
m
t
-
n
v
0
)
.Function family{
ψ
m
,
n
(
t
)
,
m
,
n
∈
Z
} is called discrete wavelet. When u
0
=
2 and v
0
=
1, function family {
ψ
m
,
n
(
t
)
} constructs a group orthonormal basis on which signals can be reconstructed.


## 3.3.2. Process of Wavelet Denoising

The basic principle of wavelet denoising is based on the spectrum distribution characteristics of signals and noise. According to the principle, noise of each decomposition level is eliminated, and useful signal is preserved [19, 20].The substance of denoising is restraining the unwanted signal and boosting the useful signal. The one-dimensionalK
i
(
i
=
1,2
,
…
,
N
w
) can be expressed as follows:
(15)
K
i
=
f
i
+
ε
e
i
,K
i is the observable SK with noise, f
i is the real SK without noise, e
i is the gauss random signal, and ε is the standard deviation of noise.The denoising process of one-dimensional signals can be divided into three steps.Step I.
Mother wavelet function and the levels of decomposition are determined, and thenK
i can be decomposed by wavelet transform.Step II.
High frequency coefficients are processed by the threshold.
The threshold thr is(16)
thr
=
σ
2
log
⁡
e
⁡
(
N
w
)
,
whereσ is the standard deviation of K
i and N
w is the length of K
i.Step III.
Signals can be reconstructed by adding the low-frequency coefficients and high frequency coefficients, and the result isK
i
′.


## 4. Feasibility Analysis of Spectral Kurtosis Algorithm

The ideal signal is defined to show the SK algorithm characteristics, and the details are as follows:(17)
i
0
(
n
)
=
i
01
(
n
)
+
i
02
(
n
)
+
i
03
(
n
)
+
i
04
(
n
)
+
ε
(
n
)
,
i
01
(
n
)
=
5.6
cos
⁡
(
2
π
×
50
t
+
π
3
)
,
i
02
(
n
)
=
40
e
-
56
t
cos
⁡
(
2
π
×
250
t
+
π
6
)
,
i
03
(
n
)
=
72
e
-
102
t
cos
⁡
(
2
π
×
315
t
)
,
i
04
(
n
)
=
10
e
-
5.5
t
,
where i
01
(
n
) is the fundamental frequency component, i
02
(
n
) is the 5th harmonic component, i
03
(
n
) is noninteger harmonic component, i
04
(
n
) is the decaying direct current (DC) component, ε
(
n
) is gauss white noise signal, and SNR
=
-
13.5
dB.The ideal TZSCi
0
(
n
) is a nonperiodic signal, added by i
01
(
n
), i
02
(
n
), i
03
(
n
), i
04
(
n
), and ε
(
n
); what should be noted is that the ideal signal is generated according to the features of the TZSC [21, 22].


### 4.1. Set the Length of Window Function

The sampling frequencyf
c = 10000 HZ, the simulation time is 0.6 s, and SK algorithm is used to analyze the ideal TZSC. The adjacent frequency interval is f
c
/
2
N
w, as N
w are 64, 128, 256, and 512, and the frequency intervals are 78.13 Hz, 39.06 Hz, 19.53 Hz, and 9.77 Hz, respectively. The results of SK based on STFT are shown in Figure 2.Length of window function.
(a)
N
w
=
64
(b)
N
w
=
128
(c)
N
w
=
256
(d)
N
w
=
512As Figure2 shows, different lengths of window function make a difference for SK. When N
w is 64 or 128, to a certain extent, the ideal TZSC can be characterized; however, the frequency resolution is lower, fundamental frequency is not imprecise, and the 5th harmonic signal cannot be characterized. When N
w is 256, frequency resolution meets the requirements, and three types of signal frequency can be distinguished, because of the interval between adjacent SK. The frequency f causes errors, and the results are f
=
58.59
Hz, f
=
253.9 Hz, and f
=
312.5 Hz. When N
w is 512, frequency resolution is higher, and fundamental frequency is more accurate; however, SK contains much noise; for example, as Figure 2(d) shows, f
=
488.3 Hz belongs to spurious frequency components.


### 4.2. Set the Type of Window Function

After the length of window function is determined, the window functions are selected, such as Gaussian window (gausswin), Hamming window, Hanning window, and Kaiser window. Based on the SK algorithm, the results are shown in Figure3.Types of window function.
(a)
Gausswin
(b)
Hamming
(c)
Hanning
(d)
KaiserIt is known from Figure3 that different window functions can also characterize the TZSC’s frequency, and the SK of white noise is close to 0. Therefore, the types of window function have a little impact on the SK. Because of being less effective than other window functions, Kaiser should avoid being used.In conclusion, SK algorithm can distinguish the frequency components of the TZSC. In order to extract the features more perfectly, Hanning window is selected in this paper, whose length is 256.


## 4.1. Set the Length of Window Function

The sampling frequencyf
c = 10000 HZ, the simulation time is 0.6 s, and SK algorithm is used to analyze the ideal TZSC. The adjacent frequency interval is f
c
/
2
N
w, as N
w are 64, 128, 256, and 512, and the frequency intervals are 78.13 Hz, 39.06 Hz, 19.53 Hz, and 9.77 Hz, respectively. The results of SK based on STFT are shown in Figure 2.Length of window function.
(a)
N
w
=
64
(b)
N
w
=
128
(c)
N
w
=
256
(d)
N
w
=
512As Figure2 shows, different lengths of window function make a difference for SK. When N
w is 64 or 128, to a certain extent, the ideal TZSC can be characterized; however, the frequency resolution is lower, fundamental frequency is not imprecise, and the 5th harmonic signal cannot be characterized. When N
w is 256, frequency resolution meets the requirements, and three types of signal frequency can be distinguished, because of the interval between adjacent SK. The frequency f causes errors, and the results are f
=
58.59
Hz, f
=
253.9 Hz, and f
=
312.5 Hz. When N
w is 512, frequency resolution is higher, and fundamental frequency is more accurate; however, SK contains much noise; for example, as Figure 2(d) shows, f
=
488.3 Hz belongs to spurious frequency components.


## 4.2. Set the Type of Window Function

After the length of window function is determined, the window functions are selected, such as Gaussian window (gausswin), Hamming window, Hanning window, and Kaiser window. Based on the SK algorithm, the results are shown in Figure3.Types of window function.
(a)
Gausswin
(b)
Hamming
(c)
Hanning
(d)
KaiserIt is known from Figure3 that different window functions can also characterize the TZSC’s frequency, and the SK of white noise is close to 0. Therefore, the types of window function have a little impact on the SK. Because of being less effective than other window functions, Kaiser should avoid being used.In conclusion, SK algorithm can distinguish the frequency components of the TZSC. In order to extract the features more perfectly, Hanning window is selected in this paper, whose length is 256.


## 5. Basic Theories about Fault Line Selection

### 5.1. Steps of Denoising

The calculation process of SK can be interfered by various noise; firstly, wavelet analysis is introduced to eliminate the noises, whose frequency is approximate to the frequency of TZSC’s SKK
i
(
i
=
1,2
,
…
,
N
w
), since symlets wavelet are limited compactly supported wavelets and whose local capabilities in the time domain and frequency domain are good, considering SK algorithm has good local capabilities in the frequency domain, and the frequency characteristic of sym6 is good, whose frequency resolution can also meet the requirements. White noise belongs to a stationary random interference signal; according to the properties (
2
) in Section 3.1, the SK of white noise is 0, but the defects from the outside condition make the SK of white noise fluctuate around 0. Finally, the threshold method is proposed to correct SK, which can also eliminate the effects on the SK relative energy entropy caused by the white noise. Specific processes are as follows.(1)
K
i is decomposed into 4 layers by sym6.
(2)
Soft threshold process the high frequency coefficients of each decomposition scale in Section3.3.2 Step  II.
(3)
Reconstruct one-dimensional wavelet by adding the bottom of the low frequency coefficients and high frequency coefficients, and the result isK
i
′.
(4)
The waveform correction coefficientλ is introduced, and the equation is
(18)
ξ
=
λ
·
K
max
⁡
′
,
0
≤
λ
≤
1
,
where K
max
⁡
′ is the max K
i
′, and considering the value relation between K
i
′ and K
max
⁡
′, the coefficient λ is 20% in the paper.
(5)
K
i
′ is retained the same value as long as K
i
′
≥
ξ; otherwise, K
i
′
=
0, the result is K
^
i.After wavelet denoising and threshold correcting, the noise inK
^
i decreases significantly, and a premise of calculating the SK relative energy entropy is accurately provided.


### 5.2. The SK Relative Energy Entropy

In information theory, entropy represents uncertainty of average information and can also be used to estimate the complexity of random signals, it provides useful information about underlying the dynamic processes of signal, SK algorithm can magnify a local frequency characteristic of the TZSC, and small and short abnormalities are found in the frequency domain by SK energy entropy. The impact of fault conditions mainly reflected on time-frequency distribution of the TZSC [23–25], in order to reflect the changes of SK in the frequency domain, based on the relative entropy and energy entropy, and the SK relative energy entropy C
j is defined in the paper:
(19)
C
j
=
-
∑
i
=
1
N
w
u
j
i
lo
g
2
u
j
i
,
j
=
1,2
,
…
,
l
,
where u
j
i is the SK relative energy probability distribution, u
j
i
=
K
^
j
i
2
/
E
max
⁡, E
max
⁡
=
max
⁡
(
E
j
), E
j
=
∑
i
=
1
M
K
^
j
i
2, the number of overhead lines is l
1, and the SK relative energy entropies are C
a
(
a
=
1,2
,
…
,
l
1
). The number of hybrid lines and cable lines is l
2, and the SK relative energy entropies are C
b
(
a
=
1,2
,
…
,
l
2
).It can be known from the above definition thatC
j is the information entropy that the SK energy of line j is relative to the max energy of all lines, the SK relative energy probability distributions are based on the sum of the largest SK energy, the probability space of SK energy is united, u
j
i is a relative quantity, which has a unique sensitivity for the change of output frequency, and C
j can describe the relative uniformity for SK energy distribution in the frequency domain.


## 5.1. Steps of Denoising

The calculation process of SK can be interfered by various noise; firstly, wavelet analysis is introduced to eliminate the noises, whose frequency is approximate to the frequency of TZSC’s SKK
i
(
i
=
1,2
,
…
,
N
w
), since symlets wavelet are limited compactly supported wavelets and whose local capabilities in the time domain and frequency domain are good, considering SK algorithm has good local capabilities in the frequency domain, and the frequency characteristic of sym6 is good, whose frequency resolution can also meet the requirements. White noise belongs to a stationary random interference signal; according to the properties (
2
) in Section 3.1, the SK of white noise is 0, but the defects from the outside condition make the SK of white noise fluctuate around 0. Finally, the threshold method is proposed to correct SK, which can also eliminate the effects on the SK relative energy entropy caused by the white noise. Specific processes are as follows.(1)
K
i is decomposed into 4 layers by sym6.
(2)
Soft threshold process the high frequency coefficients of each decomposition scale in Section3.3.2 Step  II.
(3)
Reconstruct one-dimensional wavelet by adding the bottom of the low frequency coefficients and high frequency coefficients, and the result isK
i
′.
(4)
The waveform correction coefficientλ is introduced, and the equation is
(18)
ξ
=
λ
·
K
max
⁡
′
,
0
≤
λ
≤
1
,
where K
max
⁡
′ is the max K
i
′, and considering the value relation between K
i
′ and K
max
⁡
′, the coefficient λ is 20% in the paper.
(5)
K
i
′ is retained the same value as long as K
i
′
≥
ξ; otherwise, K
i
′
=
0, the result is K
^
i.After wavelet denoising and threshold correcting, the noise inK
^
i decreases significantly, and a premise of calculating the SK relative energy entropy is accurately provided.


## 5.2. The SK Relative Energy Entropy

In information theory, entropy represents uncertainty of average information and can also be used to estimate the complexity of random signals, it provides useful information about underlying the dynamic processes of signal, SK algorithm can magnify a local frequency characteristic of the TZSC, and small and short abnormalities are found in the frequency domain by SK energy entropy. The impact of fault conditions mainly reflected on time-frequency distribution of the TZSC [23–25], in order to reflect the changes of SK in the frequency domain, based on the relative entropy and energy entropy, and the SK relative energy entropy C
j is defined in the paper:
(19)
C
j
=
-
∑
i
=
1
N
w
u
j
i
lo
g
2
u
j
i
,
j
=
1,2
,
…
,
l
,
where u
j
i is the SK relative energy probability distribution, u
j
i
=
K
^
j
i
2
/
E
max
⁡, E
max
⁡
=
max
⁡
(
E
j
), E
j
=
∑
i
=
1
M
K
^
j
i
2, the number of overhead lines is l
1, and the SK relative energy entropies are C
a
(
a
=
1,2
,
…
,
l
1
). The number of hybrid lines and cable lines is l
2, and the SK relative energy entropies are C
b
(
a
=
1,2
,
…
,
l
2
).It can be known from the above definition thatC
j is the information entropy that the SK energy of line j is relative to the max energy of all lines, the SK relative energy probability distributions are based on the sum of the largest SK energy, the probability space of SK energy is united, u
j
i is a relative quantity, which has a unique sensitivity for the change of output frequency, and C
j can describe the relative uniformity for SK energy distribution in the frequency domain.


## 6. Fault Line Selection Criterions

In [26], it can be known from the simulation waveforms that the sign of first half-wave extremes between the fault branch line and the healthy branch lines is opposite. However, the signs of first half-wave extremes among the branch lines are the same when the bus bar fault occurs. According to this, the criterion of bus bar fault is obtained.


### 6.1. Step  1—The Bus Bar Fault Criterion

(1)
The first half-wave extremesi
0
z
j of branch lines are recorded when fault occurs, and then the ratios G
0
z
m of i
0
z
j are calculated:
(20)
G
0
z
m
=
i
0
z
j
i
0
z
(
j
+
1
)
,
m
=
1,2
,
…
,
l
-
1
.
(2)
IfG
0
z
m
>
0, the bus bar is selected as the fault line, and the selection process is done; otherwise, the following steps should be executed.


### 6.2. Step  2 and Step  3—The Branch Lines Fault Criterions

The distributed complexity and uncertainty of the TZSC are reflected by the energy entropy in the frequency domain [27]; when fault occurs, the SK relative energy entropy reflects the amount of information that can be used to select the fault line. Because of the interaction between the overhead lines and cable lines, the line with the max relative energy entropy may be not the fault line; in order to avoid the influence, the relative energy of the overhead line and cable line, hybrid line would be compared in this paper, respectively.The entropy multipleβ is introduced to judge the overhead line fault, and the calculation process is as follows:
(21)
β
=
C
a
max
⁡
C
a
s
max
⁡
,
β
≥
1
,
where C
a
max
⁡ is the max of C
a and C
a
s
max
⁡ is the second biggest of C
a. After a number of simulations, the threshold of entropy multiple β is set to 2 in this paper.When the branch line fault occurs, firstly, enter step 2, judge whether the overhead line is fault or not, according to (21), calculate β, if β
≥
2 is satisfied, the overhead line of C
a
max
⁡ corresponding is selected when the fault line is selected, and the selection process is done. Otherwise, enter step 3, the hybrid line or cable line of C
b
max
⁡ corresponding is selected as the fault line, and the selection process is done.The specific selection flowchart is shown in Figure4.Figure 4
Flowchart of fault line selection.


## 6.1. Step  1—The Bus Bar Fault Criterion

(1)
The first half-wave extremesi
0
z
j of branch lines are recorded when fault occurs, and then the ratios G
0
z
m of i
0
z
j are calculated:
(20)
G
0
z
m
=
i
0
z
j
i
0
z
(
j
+
1
)
,
m
=
1,2
,
…
,
l
-
1
.
(2)
IfG
0
z
m
>
0, the bus bar is selected as the fault line, and the selection process is done; otherwise, the following steps should be executed.


## 6.2. Step  2 and Step  3—The Branch Lines Fault Criterions

The distributed complexity and uncertainty of the TZSC are reflected by the energy entropy in the frequency domain [27]; when fault occurs, the SK relative energy entropy reflects the amount of information that can be used to select the fault line. Because of the interaction between the overhead lines and cable lines, the line with the max relative energy entropy may be not the fault line; in order to avoid the influence, the relative energy of the overhead line and cable line, hybrid line would be compared in this paper, respectively.The entropy multipleβ is introduced to judge the overhead line fault, and the calculation process is as follows:
(21)
β
=
C
a
max
⁡
C
a
s
max
⁡
,
β
≥
1
,
where C
a
max
⁡ is the max of C
a and C
a
s
max
⁡ is the second biggest of C
a. After a number of simulations, the threshold of entropy multiple β is set to 2 in this paper.When the branch line fault occurs, firstly, enter step 2, judge whether the overhead line is fault or not, according to (21), calculate β, if β
≥
2 is satisfied, the overhead line of C
a
max
⁡ corresponding is selected when the fault line is selected, and the selection process is done. Otherwise, enter step 3, the hybrid line or cable line of C
b
max
⁡ corresponding is selected as the fault line, and the selection process is done.The specific selection flowchart is shown in Figure4.Figure 4
Flowchart of fault line selection.


## 7. Simulation and Verification

In this paper, the ATP-EMTP is used to simulate a single phase-to-ground fault, and the simulation model is shown in Figure5, there are 4 lines (S1, S2, S3 and S4).Figure 5
Arc suppression coil grounding system.The parameters of simulation model are as follows.Overhead line positive-sequence parameters areR
1
=
0.17
Ω
/
km, L
1
=
1.2
mH
/
km, and C
1
=
9.697
nF
/
km, and zero-sequence parameters are R
0
=
0.23
Ω
/
km, L
0
=
5.48
mH
/
km, and C
0
=
6
nF
/
km.Cable line positive-sequence parameters areR
11
=
0.193
Ω
/
km, L
11
=
0.442
mH
/and km, C
11
=
143
nF
/
km, and zero-sequence parameters are R
00
=
1.93
Ω
/
km, L
00
=
5.48
mH
/and km, C
00
=
143
nF
/
km.Transformer is 110/10.5 kV, the leakage impedance of high voltage is(
0.40
+
j
12.20
)Ω, the leakage impedance of low voltage is (
0.006
+
j
0.183
)Ω, excitation current is 0.672 A, magnetizing flux is 202.2 Wb, and magnetic resistance is 400 kΩ. Load is all are delta, Z
L
=
400
+
j
20
Ω, and arc suppression coil is L
N
=
1281.9
mH, R
N
=
40.2517
Ω.When the overhead line S1 fault occurs, the initial fault angle is 0° and the grounding resistance is 50 Ω, and the TZSC waveforms are shown in Figure 6.Transient zero-sequence current.
(a)
Overhead line S1
(b)
Overhead line S2
(c)
Hybrid line S3
(d)
Cable line S4As Figure6 shows, the TZSC is mutant when fault occurs, then they are damped. The conclusions are consistent with the analysis about the characteristics of the TZSC in Section 2; therefore, SK algorithm is introduced to analyse the fault information contained in the TZSC, and the characteristic frequency bands can be selected.


### 7.1. Select Based on Step  1

The first half-wave extremesi
0
z
j of i
0
j
(
t
) (j
=
1, 2, 3, 4) are obtained when fault occurs, and the results are in Table 1.Table 1
First half-wave extremes.
Overhead line S1
Overhead line S2
Hybrid line S3
Cable line S4
178.7320
−15.0541
−72.0865
−149.6746Becausei
0
z
1
/
i
0
z
2
=
178.7320
/
(
-
15.0541
)
=
-
11.8726
<
0, according to step 1, it can be known that the bus bar is healthy in the power distribution system. In order to select the fault line, step 2 is needed.


### 7.2. Select Based on Steps  2 and  3

The definition of characteristic frequencyf
T is the frequency f that the max SK K
max
⁡ is corresponding in the specific condition.According to the above analysis of SK algorithm in Section4, Hanning window is selected, whose length is 256. And then K
i are calculated, the results are in Figure 7 (SNR
=
-
5.5623
dB), considering too high frequency without actual physical meaning, and only the frequency of SK below 10000 Hz is listed in this paper.Correspondingf
T of SK.
(a)
S1
(b)
S2
(c)
S3
(d)
S4The following conclusions are obtained from Figure7.(1)
Several different frequency components in the TZSC can be contained. For example, two kinds of different frequency within 2000 Hz are contained in S2 at least.
(2)
Because of the flaws of SK algorithm and the impact of auxiliary white noise, the values of SK are fluctuant around 0, which are not equal to the theoretical value 0.
(3)
Different frequency components of the transient zero sequence current signals are amplified by SK; among them, the most prominent frequency components aref
T
=
781.3 Hz for S1, because the amplitude of SK in 781.3 Hz is larger than any other amplitude of SK. Similarly, the most prominent frequency components are f
T
=
390.6 Hz for S2, the most prominent frequency components is f
T
=
585.9 Hz for S3, and the most prominent frequency components is f
T
=
585.9 Hz for S4.The conclusion(
3
) shows that the characteristic frequency f
T of fault line is larger than the f
T of healthy line in this fault condition. Therefore calculate SK of each branch line, and try to select the fault line by observing the characteristic frequency f
T directly [2]. According to the SK algorithm, f
T can be obtained with different initial angles θ and grounding resistances R, Hanning window is selected, whose length is 256, SNR
=
-
5.5623
dB, and the results are shown in Table 2.Table 2
f
T of each line with S1 fault.
Initial angle
θ
Grounding resistance
Characteristic frequencyf
T (Hz)
R/Ω
Overhead line S1
Overhead line S2
Hybrid line S3
Cable line S4
0°
20
585.9
390.6
585.9
585.9
200
3320
2344
3320
3125
150°
20
585.9
390.6
585.9
585.9
200
3320
2344
3320
3125The following conclusions can be obtained from Table2.(1)
S1 is the fault line, and f
T of S1 is larger than the f
T of S2, and f
T between S3 and S4 are either equal or close.
(2)
Since the inductance of the cable lines is less far than the inductance of overhead lines, the capacitance to ground of the cable lines is larger several times than the latter, the parameter of cable lines can make a difference in the spectral characteristics of TZSC [28], and it shows that the f
T of fault overhead line S1 is either equal or close to S3 and S4.
(3)
When S1 is fault, the initial phase angle θ is from 0° up to 150° and the grounding resistance R is fixed. The f
T of same line is unchanged; however, if the grounding resistance R is from 20 Ω up to 200 Ω and the initial phase angle θ fixed, the f
T of S1 is from 585.9 Hz up to 3320 Hz, and the f
T of S2 is from 390.6 Hz up to 2344 Hz. In conclusion, the initial phase angle θ has little effect on f
T, but the grounding resistance has great effect on the f
T, because of the effect of the grounding resistance making the f
T nonunique.It can be seen that thef
T would be nonunique with different resistances; therefore, the healthy line and the fault line cannot be distinguished only by f
T,. Considering that the more obvious of the fault features are, the larger difference of entropy between the healthy line and fault lines is, so the SK relative energy entropy is introduced to distinguish the fault line and the healthy lines in this paper.To eliminate the noises, whose frequency is approximate to the frequency ofK
i, wavelet analysis is introduced to denoising before calculating the relative energy entropy, and the results are in Figure 8.Denoising spectral kurtosis.
(a)
S1
(b)
S2
(c)
S3
(d)
S4From Figure8, the noise signals decrease obviously after denoising, the values of SK fluctuate around 0 in the high frequency (over 5000 Hz), the waveforms are more smoothly, and the result is K
i
′.To decrease the interference from the SK of white noise to the SK relative energy entropy of TZSC, the threshold method is introduced to correctK
i
′ in this paper, and the finally results K
^
i are shown in Figure 9.Revising spectral kurtosis.
(a)
S1
(b)
S2
(c)
S3
(d)
S4From Figure9, most of the frequencies are below 5000 Hz in K
^
j, and the SK of white noise is set to 0 artificially by the threshold method actually. The results show that the effect of white noise is decreased, and a premise that calculates C
j accurately is provided. The results are in Table 3.Table 3
SK relative energy entropy.
Branch line
S1
S2
S3
S4
C
j
g
g
g
2.2788
0.9179
1.6239
2.2834Based on the fault line selection method of step 2, according to (21), calculate β, β
=
C
a
max
⁡
/
C
a
s
max
⁡
=
2.2788
/
0.9179
=
2.48
>
2 is satisfied, although C
4 is larger than C
1, C
2, and C
3, in view of β
>
2, and S1 is selected as the fault line.For S1 belongs to the overhead line, it can be known that the bus bar is healthy after step 1, then enter step 2, calculate β, β
>
2 is satisfied, and therefore S1 is judged as the fault line without step 3.


### 7.3. Results of Different Fault Types

To verify the accuracy of the proposed stepped fault line selection method further, the results are given when the bus bar fault occurs, the overhead line is S1, the hybrid line is S3, and the cable line is S4 with different θ and R.


#### 7.3.1. Bus Bar Fault

When the bus bar fault occurs with different initial phase anglesθ and grounding resistances R, according to step 1, the selection results are in Table 4.Table 4
First half extremes of the TZSC with bus bar faults.
Fault line
θ
R/Ω
First half extremes of TZSCi
0
z
j
G
0
z
m
>
0?
Results
S1
S2
S3
S4
Bus bar
0°
20
−36.7519
−23.2945
−157.7988
−229.0609
Y
Bus bar
200
−9.7630
−8.8730
−35.6166
−68.9971
Y
Bus bar
150°
20
22.1604
20.3579
98.0850
150.7914
Y
Bus bar
200
6.0437
5.8571
23.0452
44.5919
Y
Bus barWhen the bus bar fault occurs, it can be known from Table4 that R, θ would make a difference for i
0
z
j, which can be described as the larger R is, the lower i
0
z
j is, i
0
z
j is from negative to positive when θ is from 0° up to 150°, and the absolute value of i
0
z
j decreases. When θ
=
0°, i
0
z
j
<
0 is satisfied; however, when θ
=
150°, i
0
z
j
>
0 is satisfied, the conclusion that whatever R and θ change, G
0
z
m
>
0 is right permanently, the criterion of bus bar fault are satisfied. Therefore, the fault bus bar can be selected by step 1, the theory of bus bar fault selection method is simple, the results are correct and credible, and the influence from R and θ is avoided.


#### 7.3.2. Overhead Line Fault in S2

When the overhead line S2 fault occurs with different θ and R, according to step 1 and step 2, the selection results are in Table 5, phase A to grounding fault occurs, and the fault location is 5 km away from the bus bar.Table 5
SK relative energy entropy with S2 fault.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S2
0°
20
0.5324
3.9787
1.8911
2.4950
7.47 > 2
S2
200
1.1521
4.4164
3.5851
3.2250
3.83 > 2
S2
150°
20
0.5761
2.9607
2.2864
2.2134
5.14 > 2
S2
200
2.9918
7.0234
7.8605
5.0687
2.35 > 2
S2From Table5, the SK relative energy entropy of fault line is larger than the healthy lines in most cases. According to step 2, calculate β, because of β
=
C
a
max
⁡
/
C
a
s
max
⁡
≥
2, the overhead line is judged as fault line preliminary, because C
2 belongs to the C
a
max
⁡, S2 is selected as the fault line finally, and the selection results are accurate without step 3. When the overhead line fault occurs, the results show that utilizing the stepped fault line selection method to select the fault line only to enter step 1 and step 2, and without step 3. The selection process is simplified.From the results of Tables4 and 5, the fault line selection method not only ensures the accuracy of fault line selection but also can save the time of fault line selection by steps 1 and 2. Besides, the fault line selection method cannot be influenced by R and θ.


#### 7.3.3. Hybrid Line S3 and Cable Line S4 Fault

When the hybrid line S3 or cable line S4 fault occurs with different θ and R, respectively, according to steps 1, 2, and 3, phase A to ground fault occurs, the fault location is 5 km away from the bus bar, and the stepped fault line selection results are in Table 6.Table 6
SK relative energy entropy with S3 and S4 faults.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Hybrid line S3
0°
20
1.4287
1.5408
4.9297
4.8318
1.08 < 2
S3
200
1.0894
1.1655
3.5634
2.8949
1.07 < 2
S3
150°
20
1.1355
1.1165
4.4099
3.8247
1.02 < 2
S3
200
0.7395
0.4809
3.5649
2.1245
1.54 < 2
S3
Cable line S4
0°
20
0.8222
0.7218
3.7140
5.1669
1.14 < 2
S4
200
1.7250
1.1495
2.5230
4.0801
1.50 < 2
S4
150°
20
0.4447
0.4522
3.3416
4.1139
1.02 < 2
S4
200
2.8497
4.0614
3.1244
3.6334
1.43 < 2
S4From Table6, according to steps 1 and 2, calculate β, because of β
=
C
a
max
⁡
/
C
a
s
max
⁡
<
2, the fault line belongs to S3 or S4 preliminary, then step 3 is needed, the line of C
b
max
⁡ corresponding is selected as the fault line, the fault line S3 or S4 can be selected after step 3 finally. From the above, the fault hybrid line or the fault cable line can be selected accurately after steps 1, 2, and 3.


## 7.1. Select Based on Step  1

The first half-wave extremesi
0
z
j of i
0
j
(
t
) (j
=
1, 2, 3, 4) are obtained when fault occurs, and the results are in Table 1.Table 1
First half-wave extremes.
Overhead line S1
Overhead line S2
Hybrid line S3
Cable line S4
178.7320
−15.0541
−72.0865
−149.6746Becausei
0
z
1
/
i
0
z
2
=
178.7320
/
(
-
15.0541
)
=
-
11.8726
<
0, according to step 1, it can be known that the bus bar is healthy in the power distribution system. In order to select the fault line, step 2 is needed.


## 7.2. Select Based on Steps  2 and  3

The definition of characteristic frequencyf
T is the frequency f that the max SK K
max
⁡ is corresponding in the specific condition.According to the above analysis of SK algorithm in Section4, Hanning window is selected, whose length is 256. And then K
i are calculated, the results are in Figure 7 (SNR
=
-
5.5623
dB), considering too high frequency without actual physical meaning, and only the frequency of SK below 10000 Hz is listed in this paper.Correspondingf
T of SK.
(a)
S1
(b)
S2
(c)
S3
(d)
S4The following conclusions are obtained from Figure7.(1)
Several different frequency components in the TZSC can be contained. For example, two kinds of different frequency within 2000 Hz are contained in S2 at least.
(2)
Because of the flaws of SK algorithm and the impact of auxiliary white noise, the values of SK are fluctuant around 0, which are not equal to the theoretical value 0.
(3)
Different frequency components of the transient zero sequence current signals are amplified by SK; among them, the most prominent frequency components aref
T
=
781.3 Hz for S1, because the amplitude of SK in 781.3 Hz is larger than any other amplitude of SK. Similarly, the most prominent frequency components are f
T
=
390.6 Hz for S2, the most prominent frequency components is f
T
=
585.9 Hz for S3, and the most prominent frequency components is f
T
=
585.9 Hz for S4.The conclusion(
3
) shows that the characteristic frequency f
T of fault line is larger than the f
T of healthy line in this fault condition. Therefore calculate SK of each branch line, and try to select the fault line by observing the characteristic frequency f
T directly [2]. According to the SK algorithm, f
T can be obtained with different initial angles θ and grounding resistances R, Hanning window is selected, whose length is 256, SNR
=
-
5.5623
dB, and the results are shown in Table 2.Table 2
f
T of each line with S1 fault.
Initial angle
θ
Grounding resistance
Characteristic frequencyf
T (Hz)
R/Ω
Overhead line S1
Overhead line S2
Hybrid line S3
Cable line S4
0°
20
585.9
390.6
585.9
585.9
200
3320
2344
3320
3125
150°
20
585.9
390.6
585.9
585.9
200
3320
2344
3320
3125The following conclusions can be obtained from Table2.(1)
S1 is the fault line, and f
T of S1 is larger than the f
T of S2, and f
T between S3 and S4 are either equal or close.
(2)
Since the inductance of the cable lines is less far than the inductance of overhead lines, the capacitance to ground of the cable lines is larger several times than the latter, the parameter of cable lines can make a difference in the spectral characteristics of TZSC [28], and it shows that the f
T of fault overhead line S1 is either equal or close to S3 and S4.
(3)
When S1 is fault, the initial phase angle θ is from 0° up to 150° and the grounding resistance R is fixed. The f
T of same line is unchanged; however, if the grounding resistance R is from 20 Ω up to 200 Ω and the initial phase angle θ fixed, the f
T of S1 is from 585.9 Hz up to 3320 Hz, and the f
T of S2 is from 390.6 Hz up to 2344 Hz. In conclusion, the initial phase angle θ has little effect on f
T, but the grounding resistance has great effect on the f
T, because of the effect of the grounding resistance making the f
T nonunique.It can be seen that thef
T would be nonunique with different resistances; therefore, the healthy line and the fault line cannot be distinguished only by f
T,. Considering that the more obvious of the fault features are, the larger difference of entropy between the healthy line and fault lines is, so the SK relative energy entropy is introduced to distinguish the fault line and the healthy lines in this paper.To eliminate the noises, whose frequency is approximate to the frequency ofK
i, wavelet analysis is introduced to denoising before calculating the relative energy entropy, and the results are in Figure 8.Denoising spectral kurtosis.
(a)
S1
(b)
S2
(c)
S3
(d)
S4From Figure8, the noise signals decrease obviously after denoising, the values of SK fluctuate around 0 in the high frequency (over 5000 Hz), the waveforms are more smoothly, and the result is K
i
′.To decrease the interference from the SK of white noise to the SK relative energy entropy of TZSC, the threshold method is introduced to correctK
i
′ in this paper, and the finally results K
^
i are shown in Figure 9.Revising spectral kurtosis.
(a)
S1
(b)
S2
(c)
S3
(d)
S4From Figure9, most of the frequencies are below 5000 Hz in K
^
j, and the SK of white noise is set to 0 artificially by the threshold method actually. The results show that the effect of white noise is decreased, and a premise that calculates C
j accurately is provided. The results are in Table 3.Table 3
SK relative energy entropy.
Branch line
S1
S2
S3
S4
C
j
g
g
g
2.2788
0.9179
1.6239
2.2834Based on the fault line selection method of step 2, according to (21), calculate β, β
=
C
a
max
⁡
/
C
a
s
max
⁡
=
2.2788
/
0.9179
=
2.48
>
2 is satisfied, although C
4 is larger than C
1, C
2, and C
3, in view of β
>
2, and S1 is selected as the fault line.For S1 belongs to the overhead line, it can be known that the bus bar is healthy after step 1, then enter step 2, calculate β, β
>
2 is satisfied, and therefore S1 is judged as the fault line without step 3.


## 7.3. Results of Different Fault Types

To verify the accuracy of the proposed stepped fault line selection method further, the results are given when the bus bar fault occurs, the overhead line is S1, the hybrid line is S3, and the cable line is S4 with different θ and R.


### 7.3.1. Bus Bar Fault

When the bus bar fault occurs with different initial phase anglesθ and grounding resistances R, according to step 1, the selection results are in Table 4.Table 4
First half extremes of the TZSC with bus bar faults.
Fault line
θ
R/Ω
First half extremes of TZSCi
0
z
j
G
0
z
m
>
0?
Results
S1
S2
S3
S4
Bus bar
0°
20
−36.7519
−23.2945
−157.7988
−229.0609
Y
Bus bar
200
−9.7630
−8.8730
−35.6166
−68.9971
Y
Bus bar
150°
20
22.1604
20.3579
98.0850
150.7914
Y
Bus bar
200
6.0437
5.8571
23.0452
44.5919
Y
Bus barWhen the bus bar fault occurs, it can be known from Table4 that R, θ would make a difference for i
0
z
j, which can be described as the larger R is, the lower i
0
z
j is, i
0
z
j is from negative to positive when θ is from 0° up to 150°, and the absolute value of i
0
z
j decreases. When θ
=
0°, i
0
z
j
<
0 is satisfied; however, when θ
=
150°, i
0
z
j
>
0 is satisfied, the conclusion that whatever R and θ change, G
0
z
m
>
0 is right permanently, the criterion of bus bar fault are satisfied. Therefore, the fault bus bar can be selected by step 1, the theory of bus bar fault selection method is simple, the results are correct and credible, and the influence from R and θ is avoided.


### 7.3.2. Overhead Line Fault in S2

When the overhead line S2 fault occurs with different θ and R, according to step 1 and step 2, the selection results are in Table 5, phase A to grounding fault occurs, and the fault location is 5 km away from the bus bar.Table 5
SK relative energy entropy with S2 fault.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S2
0°
20
0.5324
3.9787
1.8911
2.4950
7.47 > 2
S2
200
1.1521
4.4164
3.5851
3.2250
3.83 > 2
S2
150°
20
0.5761
2.9607
2.2864
2.2134
5.14 > 2
S2
200
2.9918
7.0234
7.8605
5.0687
2.35 > 2
S2From Table5, the SK relative energy entropy of fault line is larger than the healthy lines in most cases. According to step 2, calculate β, because of β
=
C
a
max
⁡
/
C
a
s
max
⁡
≥
2, the overhead line is judged as fault line preliminary, because C
2 belongs to the C
a
max
⁡, S2 is selected as the fault line finally, and the selection results are accurate without step 3. When the overhead line fault occurs, the results show that utilizing the stepped fault line selection method to select the fault line only to enter step 1 and step 2, and without step 3. The selection process is simplified.From the results of Tables4 and 5, the fault line selection method not only ensures the accuracy of fault line selection but also can save the time of fault line selection by steps 1 and 2. Besides, the fault line selection method cannot be influenced by R and θ.


### 7.3.3. Hybrid Line S3 and Cable Line S4 Fault

When the hybrid line S3 or cable line S4 fault occurs with different θ and R, respectively, according to steps 1, 2, and 3, phase A to ground fault occurs, the fault location is 5 km away from the bus bar, and the stepped fault line selection results are in Table 6.Table 6
SK relative energy entropy with S3 and S4 faults.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Hybrid line S3
0°
20
1.4287
1.5408
4.9297
4.8318
1.08 < 2
S3
200
1.0894
1.1655
3.5634
2.8949
1.07 < 2
S3
150°
20
1.1355
1.1165
4.4099
3.8247
1.02 < 2
S3
200
0.7395
0.4809
3.5649
2.1245
1.54 < 2
S3
Cable line S4
0°
20
0.8222
0.7218
3.7140
5.1669
1.14 < 2
S4
200
1.7250
1.1495
2.5230
4.0801
1.50 < 2
S4
150°
20
0.4447
0.4522
3.3416
4.1139
1.02 < 2
S4
200
2.8497
4.0614
3.1244
3.6334
1.43 < 2
S4From Table6, according to steps 1 and 2, calculate β, because of β
=
C
a
max
⁡
/
C
a
s
max
⁡
<
2, the fault line belongs to S3 or S4 preliminary, then step 3 is needed, the line of C
b
max
⁡ corresponding is selected as the fault line, the fault line S3 or S4 can be selected after step 3 finally. From the above, the fault hybrid line or the fault cable line can be selected accurately after steps 1, 2, and 3.


## 7.3.1. Bus Bar Fault

When the bus bar fault occurs with different initial phase anglesθ and grounding resistances R, according to step 1, the selection results are in Table 4.Table 4
First half extremes of the TZSC with bus bar faults.
Fault line
θ
R/Ω
First half extremes of TZSCi
0
z
j
G
0
z
m
>
0?
Results
S1
S2
S3
S4
Bus bar
0°
20
−36.7519
−23.2945
−157.7988
−229.0609
Y
Bus bar
200
−9.7630
−8.8730
−35.6166
−68.9971
Y
Bus bar
150°
20
22.1604
20.3579
98.0850
150.7914
Y
Bus bar
200
6.0437
5.8571
23.0452
44.5919
Y
Bus barWhen the bus bar fault occurs, it can be known from Table4 that R, θ would make a difference for i
0
z
j, which can be described as the larger R is, the lower i
0
z
j is, i
0
z
j is from negative to positive when θ is from 0° up to 150°, and the absolute value of i
0
z
j decreases. When θ
=
0°, i
0
z
j
<
0 is satisfied; however, when θ
=
150°, i
0
z
j
>
0 is satisfied, the conclusion that whatever R and θ change, G
0
z
m
>
0 is right permanently, the criterion of bus bar fault are satisfied. Therefore, the fault bus bar can be selected by step 1, the theory of bus bar fault selection method is simple, the results are correct and credible, and the influence from R and θ is avoided.


## 7.3.2. Overhead Line Fault in S2

When the overhead line S2 fault occurs with different θ and R, according to step 1 and step 2, the selection results are in Table 5, phase A to grounding fault occurs, and the fault location is 5 km away from the bus bar.Table 5
SK relative energy entropy with S2 fault.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S2
0°
20
0.5324
3.9787
1.8911
2.4950
7.47 > 2
S2
200
1.1521
4.4164
3.5851
3.2250
3.83 > 2
S2
150°
20
0.5761
2.9607
2.2864
2.2134
5.14 > 2
S2
200
2.9918
7.0234
7.8605
5.0687
2.35 > 2
S2From Table5, the SK relative energy entropy of fault line is larger than the healthy lines in most cases. According to step 2, calculate β, because of β
=
C
a
max
⁡
/
C
a
s
max
⁡
≥
2, the overhead line is judged as fault line preliminary, because C
2 belongs to the C
a
max
⁡, S2 is selected as the fault line finally, and the selection results are accurate without step 3. When the overhead line fault occurs, the results show that utilizing the stepped fault line selection method to select the fault line only to enter step 1 and step 2, and without step 3. The selection process is simplified.From the results of Tables4 and 5, the fault line selection method not only ensures the accuracy of fault line selection but also can save the time of fault line selection by steps 1 and 2. Besides, the fault line selection method cannot be influenced by R and θ.


## 7.3.3. Hybrid Line S3 and Cable Line S4 Fault

When the hybrid line S3 or cable line S4 fault occurs with different θ and R, respectively, according to steps 1, 2, and 3, phase A to ground fault occurs, the fault location is 5 km away from the bus bar, and the stepped fault line selection results are in Table 6.Table 6
SK relative energy entropy with S3 and S4 faults.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Hybrid line S3
0°
20
1.4287
1.5408
4.9297
4.8318
1.08 < 2
S3
200
1.0894
1.1655
3.5634
2.8949
1.07 < 2
S3
150°
20
1.1355
1.1165
4.4099
3.8247
1.02 < 2
S3
200
0.7395
0.4809
3.5649
2.1245
1.54 < 2
S3
Cable line S4
0°
20
0.8222
0.7218
3.7140
5.1669
1.14 < 2
S4
200
1.7250
1.1495
2.5230
4.0801
1.50 < 2
S4
150°
20
0.4447
0.4522
3.3416
4.1139
1.02 < 2
S4
200
2.8497
4.0614
3.1244
3.6334
1.43 < 2
S4From Table6, according to steps 1 and 2, calculate β, because of β
=
C
a
max
⁡
/
C
a
s
max
⁡
<
2, the fault line belongs to S3 or S4 preliminary, then step 3 is needed, the line of C
b
max
⁡ corresponding is selected as the fault line, the fault line S3 or S4 can be selected after step 3 finally. From the above, the fault hybrid line or the fault cable line can be selected accurately after steps 1, 2, and 3.


## 8. Applicability Analyses

### 8.1. Different Fault Distances

To verify the proposed method adaptability with different fault distances,θ and R are fixed, and the simulations are as follows. Taking the fault line S1 and S4 as examples, the fault locations are 10 km and 2 km away from the bus bar, respectively; the results are in Table 7.Table 7
SK relative energy entropy with different fault distances (θ
=
0
°, R
=
20
Ω).
Fault line
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S1
2.7267
1.2616
2.1665
2.8215
2.16 > 2
S1
Cable line S4
1.6140
1.1455
3.9845
4.3558
1.41 < 2
S4From Table7, when S1 fault occurs with different fault distances, based on step 2, because of β
1
=
2.7267
/
1.2616
=
2.16
>
2, step 3 is not needed, and S1 is selected as the fault line finally. Similarly, based on step 2 when S4 fault occurs, because of β
2
=
1.6140
/
1.1455
=
1.41
<
2, step 3 is needed, since 4.3558
>
3.9845, and S4 is selected as the fault line finally. The results show the stepped fault line selection method still applies to the changing fault distances.


### 8.2. Different Fault Phases

To verify the proposed method adaptability about different fault phases, taking the fault line S1 and S4 as examples, the results are as shown in Table 8.Table 8
SK relative energy entropy with different fault phases (θ
=
0
°, R
=
20
Ω).
Fault line
Fault phase
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S1
B
3.0161
1.3442
1.7646
2.9313
2.24 > 2
S1
C
2.1819
0.4247
1.0312
1.4207
5.14 > 2
S1
Cable line S4
B
0.6462
0.7357
3.5225
5.2236
1.14 < 2
S4
C
2.0681
3.6244
1.2483
1.4300
1.75 < 2
S4From Table8, when S1 fault occurs with different fault phases, based on step 2, because of β
1
=
2.7267
/
1.2616
=
2.16
>
2 and β
2
=
2.1819
/
0.4247
=
5.15
>
2, step 3 is not needed, and S1 is selected as the fault line finally. Similarly, based on step 2 when S4 fault with the changing fault phases, since β
3
=
0.7357
/
0.6462
=
1.14
<
2 and β
4
=
3.6244
/
2.0681
=
1.75
<
2, step 3 is needed, because of 5.2236
>
3.5225 and 1.4300
>
1.2483, and S4 is selected as the fault line finally. The results show the stepped fault line selection method still applies to the changing fault phases.


### 8.3. Simulate the Practical Data

Because the conditions of practical environment are poor when fault occurs, therefore, the practical fault data contains much noise. In this paper, in order to analyze the differences between the practical data and the simulation data intuitively, the different intensity noises are injected into the TZSC to simulate the practical data, when SNR is −1.5678 dB and −13.5678 dB, and the results are as shown in Figure10.Simulate the practical fault TZSC.
(a)
TZSC
(b)
SNR = −1.5678 dB
(c)
SNR = −13.5678 dBFrom Figure10, it can be known that the simulation data can simulate the practical data with the noise background. That is to say, it is possible to verify the proposed algorithm using the practical data. Taking the fault lines S1 and S4 as examples, the results are shown in Table 9.Table 9
SK relative energy entropy with different intensity noise (θ
=
0
°, R
=
20
Ω).
SNR/dB
Fault line
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
−1.5678
Overhead line S1
3.0665
1.2711
1.3520
2.7091
2.41 > 2
S1
Cable line S4
0.9898
0.8089
4.5347
5.6922
1.22 < 2
S4
−5.5678
Overhead line S1
2.9061
1.1725
1.8192
2.7349
2.48 > 2
S1
Cable line S4
0.9375
0.5388
3.7573
5.2826
1.74 < 2
S4
−10.5678
Overhead line S1
2.8449
1.0359
1.7698
2.4071
2.75 > 2
S1
Cable line S4
0.8267
0.4481
3.4586
3.9673
1.84 < 2
S4
−13.5678
Overhead line S1
2.6814
0.6951
2.0311
2.1538
3.86 > 2
S1
Cable line S4
1.1979
0.7646
4.5726
5.8559
1.57 < 2
S4From Table9, when overhead line S1 fault occurs, no matter what SNR is equal to, according to steps 1 and 2, β
>
2 establish permanently, step 3 is not needed, and S1 is selected as the fault line finally. Similarly, when S4 fault occurs, no matter what SNR is equal to, based on steps 1, 2, and 3, β
<
2 establish all the time, step 3 is needed, since C
4
>
C
3, and S4 is selected as fault line finally. It can be seen that the different fault lines can be selected correctly with different intensity noise, and the results show the stepped fault line selection method can apply to the simulative practical fault data.


### 8.4. Incipient Faults

Compared with the overhead line, the cable line is safer and more economic, so the cable lines are widely used in modern cities. However, the cables are more easily in the state of incipient faults before they fail into permanent faults. Such as the cable line S4 is in the state of incipient faults in this paper. Usually, incipient faults in power cables gradually result from the aging process, where the localized deterioration in insulation exists. Incipient faults are normally characterized as the faulty phenomena with the relatively small fault currents and the relatively short duration ranging from one-quarter cycle to multicycle. These short lasting current variations cannot be detected by the traditional distribution protection schemes because of their short duration and low increment in magnitude. However, such faults must be detected at the early stage to avoid the consequent catastrophe induced by the degradation themselves [29].The incipient faults of underground cables are similar to the intermittent arc faults; therefore, the arc model is essential to the analysis of the incipient faults. Paper [30] introduces the arc model; after theoretical analysis and practical verification, arc has nonlinear and time-varying characteristics, and the high-frequency components can be produced. Due to the fact that the arc model is easy, the nonlinear and time-varying characteristics can be expressed, and the arc model is introduced in this paper.The arc model is built by the energy balance of arc column, and the differential equations of conductivity are used to express(22)
d
p
d
t
=
1
τ
(
P
-
p
)
,
where τ is the time constant of arc, p is the instantaneous conductivity of arc, and P is stable conductivity of arc.Stable conductivity can be defined as(23)
P
=
|
i
arc
|
u
st
,
u
st
=
u
0
+
r
0
·
|
i
arc
|
,
where i
arc is the instantaneous current of arc, u
st is static voltage of arc, u
0 is the voltage characteristics of arc, r
0 is the arc resistance, and the characteristic parameters u
0 and r
0 depend on the length l
arc of arc, which can be calculated by
(24)
u
0
=
(
0.9
l
arc
m
+
0.4
)
kV
(25)
r
0
=
(
40
l
arc
m
+
8
)
m
Ω
.Equation (25) is generalized arc equation. In the resonance grounding system, the parameters of the small current arc change depend on the length of arc, which is the arc elongation e
l
o
n
g
s
p
d. e
l
o
n
g
s
p
d can be defined as
(26)
e
l
o
n
g
s
p
d
=
d
l
arc
d
t
=
7
l
0
(
0.2
v
th
/
v
max
⁡
)
+
0.2
,
where v
th is the instantaneous value of transient initial voltage and v
max
⁡ is the max normal voltage.The length of arc is defined as a time constant before simulating. The time constant of arc is defined as(27)
τ
=
τ
0
·
(
l
arc
l
0
)
α
,
where τ
0 is the initial time constant, l
0 is the initial length of arc, α is the negative coefficient, and the general value of α is −4.The length of arc is(28)
l
arc
=
(
e
l
o
n
g
s
p
d
·
(
t
-
t
f
)
+
1
)
·
l
0
,
where t is the simulation time and t
f is the moment that the arc occurs.In underground cables, the incipient fault is one type of transient, which is prone to an intermittent arc fault. The typical incipient faults are composed of two types: subcycle and multicycle incipient fault. The subcycle incipient fault always occurs near a voltage peak where the arc is ignited, lasts around one-quarter cycle, and self-clears when the current crosses zero. Figure11(a) shows the three-phase feeder currents when a subcycle incipient fault occurs between phase A to ground. The multicycle incipient fault also likely occurs near a voltage peak, lasts 1–4 cycles, and self-clears when the arc is quenched. The waveforms of the currents are shown in Figure 11(b).Feeder currents of incipient faults.
(a)
Subcycle
(b)
MulticycleTo verify the adaptability of the proposed method when the incipient faults occur, the TZSC of incipient faults is obtained, and then, according to the steps 1, 2, and 3, the selection results are as shown in Table10.Table 10
SK relative energy entropy with incipient faults.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Cable line S4
0°
1
1.1021
0.6623
4.7252
5.6418
1.66 < 2
S4
50
0.3872
0.3602
3.5047
4.1438
1.07 < 2From Table10, based on step 2 when S4 occur incipient faults, it can be known that β
1
=
1.1021
/
0.6623
=
1.66
<
2 and β
2
=
0.3872
/
0.3602
=
1.07
<
2, step 3 is needed, and since 5.6418
>
4.7252 and 4.1438
>
3.5047, S4 is selected as the fault line finally. The results show that the stepped fault line selection method can also adapt to the incipient faults.


## 8.1. Different Fault Distances

To verify the proposed method adaptability with different fault distances,θ and R are fixed, and the simulations are as follows. Taking the fault line S1 and S4 as examples, the fault locations are 10 km and 2 km away from the bus bar, respectively; the results are in Table 7.Table 7
SK relative energy entropy with different fault distances (θ
=
0
°, R
=
20
Ω).
Fault line
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S1
2.7267
1.2616
2.1665
2.8215
2.16 > 2
S1
Cable line S4
1.6140
1.1455
3.9845
4.3558
1.41 < 2
S4From Table7, when S1 fault occurs with different fault distances, based on step 2, because of β
1
=
2.7267
/
1.2616
=
2.16
>
2, step 3 is not needed, and S1 is selected as the fault line finally. Similarly, based on step 2 when S4 fault occurs, because of β
2
=
1.6140
/
1.1455
=
1.41
<
2, step 3 is needed, since 4.3558
>
3.9845, and S4 is selected as the fault line finally. The results show the stepped fault line selection method still applies to the changing fault distances.


## 8.2. Different Fault Phases

To verify the proposed method adaptability about different fault phases, taking the fault line S1 and S4 as examples, the results are as shown in Table 8.Table 8
SK relative energy entropy with different fault phases (θ
=
0
°, R
=
20
Ω).
Fault line
Fault phase
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Overhead line S1
B
3.0161
1.3442
1.7646
2.9313
2.24 > 2
S1
C
2.1819
0.4247
1.0312
1.4207
5.14 > 2
S1
Cable line S4
B
0.6462
0.7357
3.5225
5.2236
1.14 < 2
S4
C
2.0681
3.6244
1.2483
1.4300
1.75 < 2
S4From Table8, when S1 fault occurs with different fault phases, based on step 2, because of β
1
=
2.7267
/
1.2616
=
2.16
>
2 and β
2
=
2.1819
/
0.4247
=
5.15
>
2, step 3 is not needed, and S1 is selected as the fault line finally. Similarly, based on step 2 when S4 fault with the changing fault phases, since β
3
=
0.7357
/
0.6462
=
1.14
<
2 and β
4
=
3.6244
/
2.0681
=
1.75
<
2, step 3 is needed, because of 5.2236
>
3.5225 and 1.4300
>
1.2483, and S4 is selected as the fault line finally. The results show the stepped fault line selection method still applies to the changing fault phases.


## 8.3. Simulate the Practical Data

Because the conditions of practical environment are poor when fault occurs, therefore, the practical fault data contains much noise. In this paper, in order to analyze the differences between the practical data and the simulation data intuitively, the different intensity noises are injected into the TZSC to simulate the practical data, when SNR is −1.5678 dB and −13.5678 dB, and the results are as shown in Figure10.Simulate the practical fault TZSC.
(a)
TZSC
(b)
SNR = −1.5678 dB
(c)
SNR = −13.5678 dBFrom Figure10, it can be known that the simulation data can simulate the practical data with the noise background. That is to say, it is possible to verify the proposed algorithm using the practical data. Taking the fault lines S1 and S4 as examples, the results are shown in Table 9.Table 9
SK relative energy entropy with different intensity noise (θ
=
0
°, R
=
20
Ω).
SNR/dB
Fault line
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
−1.5678
Overhead line S1
3.0665
1.2711
1.3520
2.7091
2.41 > 2
S1
Cable line S4
0.9898
0.8089
4.5347
5.6922
1.22 < 2
S4
−5.5678
Overhead line S1
2.9061
1.1725
1.8192
2.7349
2.48 > 2
S1
Cable line S4
0.9375
0.5388
3.7573
5.2826
1.74 < 2
S4
−10.5678
Overhead line S1
2.8449
1.0359
1.7698
2.4071
2.75 > 2
S1
Cable line S4
0.8267
0.4481
3.4586
3.9673
1.84 < 2
S4
−13.5678
Overhead line S1
2.6814
0.6951
2.0311
2.1538
3.86 > 2
S1
Cable line S4
1.1979
0.7646
4.5726
5.8559
1.57 < 2
S4From Table9, when overhead line S1 fault occurs, no matter what SNR is equal to, according to steps 1 and 2, β
>
2 establish permanently, step 3 is not needed, and S1 is selected as the fault line finally. Similarly, when S4 fault occurs, no matter what SNR is equal to, based on steps 1, 2, and 3, β
<
2 establish all the time, step 3 is needed, since C
4
>
C
3, and S4 is selected as fault line finally. It can be seen that the different fault lines can be selected correctly with different intensity noise, and the results show the stepped fault line selection method can apply to the simulative practical fault data.


## 8.4. Incipient Faults

Compared with the overhead line, the cable line is safer and more economic, so the cable lines are widely used in modern cities. However, the cables are more easily in the state of incipient faults before they fail into permanent faults. Such as the cable line S4 is in the state of incipient faults in this paper. Usually, incipient faults in power cables gradually result from the aging process, where the localized deterioration in insulation exists. Incipient faults are normally characterized as the faulty phenomena with the relatively small fault currents and the relatively short duration ranging from one-quarter cycle to multicycle. These short lasting current variations cannot be detected by the traditional distribution protection schemes because of their short duration and low increment in magnitude. However, such faults must be detected at the early stage to avoid the consequent catastrophe induced by the degradation themselves [29].The incipient faults of underground cables are similar to the intermittent arc faults; therefore, the arc model is essential to the analysis of the incipient faults. Paper [30] introduces the arc model; after theoretical analysis and practical verification, arc has nonlinear and time-varying characteristics, and the high-frequency components can be produced. Due to the fact that the arc model is easy, the nonlinear and time-varying characteristics can be expressed, and the arc model is introduced in this paper.The arc model is built by the energy balance of arc column, and the differential equations of conductivity are used to express(22)
d
p
d
t
=
1
τ
(
P
-
p
)
,
where τ is the time constant of arc, p is the instantaneous conductivity of arc, and P is stable conductivity of arc.Stable conductivity can be defined as(23)
P
=
|
i
arc
|
u
st
,
u
st
=
u
0
+
r
0
·
|
i
arc
|
,
where i
arc is the instantaneous current of arc, u
st is static voltage of arc, u
0 is the voltage characteristics of arc, r
0 is the arc resistance, and the characteristic parameters u
0 and r
0 depend on the length l
arc of arc, which can be calculated by
(24)
u
0
=
(
0.9
l
arc
m
+
0.4
)
kV
(25)
r
0
=
(
40
l
arc
m
+
8
)
m
Ω
.Equation (25) is generalized arc equation. In the resonance grounding system, the parameters of the small current arc change depend on the length of arc, which is the arc elongation e
l
o
n
g
s
p
d. e
l
o
n
g
s
p
d can be defined as
(26)
e
l
o
n
g
s
p
d
=
d
l
arc
d
t
=
7
l
0
(
0.2
v
th
/
v
max
⁡
)
+
0.2
,
where v
th is the instantaneous value of transient initial voltage and v
max
⁡ is the max normal voltage.The length of arc is defined as a time constant before simulating. The time constant of arc is defined as(27)
τ
=
τ
0
·
(
l
arc
l
0
)
α
,
where τ
0 is the initial time constant, l
0 is the initial length of arc, α is the negative coefficient, and the general value of α is −4.The length of arc is(28)
l
arc
=
(
e
l
o
n
g
s
p
d
·
(
t
-
t
f
)
+
1
)
·
l
0
,
where t is the simulation time and t
f is the moment that the arc occurs.In underground cables, the incipient fault is one type of transient, which is prone to an intermittent arc fault. The typical incipient faults are composed of two types: subcycle and multicycle incipient fault. The subcycle incipient fault always occurs near a voltage peak where the arc is ignited, lasts around one-quarter cycle, and self-clears when the current crosses zero. Figure11(a) shows the three-phase feeder currents when a subcycle incipient fault occurs between phase A to ground. The multicycle incipient fault also likely occurs near a voltage peak, lasts 1–4 cycles, and self-clears when the arc is quenched. The waveforms of the currents are shown in Figure 11(b).Feeder currents of incipient faults.
(a)
Subcycle
(b)
MulticycleTo verify the adaptability of the proposed method when the incipient faults occur, the TZSC of incipient faults is obtained, and then, according to the steps 1, 2, and 3, the selection results are as shown in Table10.Table 10
SK relative energy entropy with incipient faults.
Fault line
θ
R/Ω
SK relative energy entropyC
j
Entropy multipleβ
Results
S1
S2
S3
S4
Cable line S4
0°
1
1.1021
0.6623
4.7252
5.6418
1.66 < 2
S4
50
0.3872
0.3602
3.5047
4.1438
1.07 < 2From Table10, based on step 2 when S4 occur incipient faults, it can be known that β
1
=
1.1021
/
0.6623
=
1.66
<
2 and β
2
=
0.3872
/
0.3602
=
1.07
<
2, step 3 is needed, and since 5.6418
>
4.7252 and 4.1438
>
3.5047, S4 is selected as the fault line finally. The results show that the stepped fault line selection method can also adapt to the incipient faults.


## 9. Compared with Other Fault Line Selection Methods

### 9.1. Wavelet Relative Energy Entropy

The wavelet mother function “db10” is used to decompose the TZSC [12], and the decomposed layers are 4. Select the frequency bands [4  0], [4  1], [4  2], and [4  3], and the results are shown in Figure 12.Wavelet decompose the TZSC.
(a)
(b)From Figure12, the lowest frequency band [40] is removed, the biggest energy band [4  1] is selected as the characteristic band, because the characteristic band contains the main features of TZSC. And then wavelet denoising and the threshold are used to correct the characteristic band coefficients, successively. Finally, the wavelet relative energy entropies C
j can be obtained, and the selection results are as shown in Table 11.Table 11
Wavelet relative energy entropy.
Fault line
θ
R/Ω
Wavelet relative energy entropyC
j
Entropy multipleβ
Results
Right?
S1
S2
S3
S4
Overhead line S2
0°
20
2.1015
2.4464
2.6801
4.1830
1.16 < 2
S4
N
200
1.2593
2.1661
3.6776
3.3619
1.72 < 2
S3
N
150°
20
2.1669
2.5476
2.5867
4.4014
1.18 < 2
S4
N
200
1.1975
1.9841
3.6561
3.6506
1.66 < 2
S3
NFrom Table11, the method based on wavelet relative energy entropy cannot select the fault line accurately. When the overhead line S2 fault occurs, according to (21), the entropy multiple β can be calculated, C
1
>
C
2, while β
<
2 is established, the selection results are incorrectly, because of the big differences on the distribution parameters between the cable line and the overhead line, which make the healthy cable line’s entropies bigger than the fault overhead line. Moreover, before the wavelet transform, the appropriate mother wavelet function should be selected, while the selection process is difficult. It can be seen that the accuracy of the method based on wavelet relative energy entropy is not guaranteed.


### 9.2. EMD Relative Energy Entropy

The local characteristics of the TZSC can be represented by empirical mode decomposition (EMD); after decomposition, IMF components can be obtained [8], and the waveforms of the mode function IMF1~IMF4 are as shown in Figure 13.EMD decompose the TZSC.
(a)
TZSC
(b)
IMF1
(c)
IMF2
(d)
IMF3
(e)
IMF4In view of IMF1 component belonging the highest frequency component of the TZSC, containing rich transient information, so IMF1 is selected as the characteristic component, and then wavelet denoising and the threshold are used to correct the characteristic component coefficients. Finally, the relative energy entropies of IMF1 can be obtained, and the selection results are as shown in Table12.Table 12
EMD relative energy entropy.
Fault line
θ
R/Ω
EMD relative energy entropyC
j
Entropy multipleβ
Results
Right?
S1
S2
S3
S4
Overhead line S2
0°
20
1.9112
4.5143
5.1857
4.9335
2.36 > 2
S2
Y
200
0.8032
3.8554
5.9764
4.8981
4.80 > 2
S2
Y
150°
20
2.0273
4.4264
5.7486
4.3812
2.18 > 2
S2
Y
200
0.9719
4.1269
6.8780
6.2598
4.25 > 2
S2
Y
Hybrid line S3
0°
20
0.7617
1.1209
4.8099
2.8481
1.47 < 2
S3
Y
200
0.6714
0.9501
4.7789
2.9695
1.12 < 2
S3
Y
150°
20
1.0806
1.5968
5.3013
4.0003
1.48 < 2
S3
Y
200
0.7927
1.3501
5.0642
2.7225
1.70 < 2
S3
Y
Cable line S4
0°
20
0.5502
0.9428
5.4257
4.2543
1.71 < 2
S4
Y
200
0.2029
0.7532
2.0443
3.7942
3.71 > 2
S2
N
150°
20
0.2706
0.3861
4.8885
4.0831
1.43 < 2
S3
N
200
0.3037
1.5814
4.2905
5.5391
5.21 > 2
S2
NFrom Table12, when the cable line S4 fault occurs, after decomposition by EMD, the EMD relative energy entropy of each branch line can be obtained, and then the entropy multiple β can be calculated, β is bigger than the threshold, according to steps 1 and 2, the cable line S4 belongs to the healthy line, and the selection results are incorrect. Moreover, the decomposition process may cause modal aliasing, which can make a big difference on IMF1 component to represent the TZSC. Therefore, it can be seen that the accuracy of the method based on EMD relative energy entropy is not guaranteed.From Tables5, 6, 11, and 12, when the branch line fault occurs, because of the big difference of distribution parameters between the cable line and overhead line, the relative energy entropy of each branch line can be obtained, according to the fault line selection steps 1, 2, and 3, and the accuracy of the method based on wavelet relative energy entropy and EMD relative energy entropy is not guaranteed. Moreover, wavelet transform needs to select the appropriate mother wavelet function, while the selection process is difficult and the modal aliasing may be caused by EMD decomposition process, while SK algorithm can decompose the TZSC in time-frequency domain adaptively, which contributes to the accuracy of the selection results. Therefore, the accuracy of selection method based on SK relative energy entropy is higher than the methods based on the wavelet relative energy entropy and EMD relative energy entropy.


## 9.1. Wavelet Relative Energy Entropy

The wavelet mother function “db10” is used to decompose the TZSC [12], and the decomposed layers are 4. Select the frequency bands [4  0], [4  1], [4  2], and [4  3], and the results are shown in Figure 12.Wavelet decompose the TZSC.
(a)
(b)From Figure12, the lowest frequency band [40] is removed, the biggest energy band [4  1] is selected as the characteristic band, because the characteristic band contains the main features of TZSC. And then wavelet denoising and the threshold are used to correct the characteristic band coefficients, successively. Finally, the wavelet relative energy entropies C
j can be obtained, and the selection results are as shown in Table 11.Table 11
Wavelet relative energy entropy.
Fault line
θ
R/Ω
Wavelet relative energy entropyC
j
Entropy multipleβ
Results
Right?
S1
S2
S3
S4
Overhead line S2
0°
20
2.1015
2.4464
2.6801
4.1830
1.16 < 2
S4
N
200
1.2593
2.1661
3.6776
3.3619
1.72 < 2
S3
N
150°
20
2.1669
2.5476
2.5867
4.4014
1.18 < 2
S4
N
200
1.1975
1.9841
3.6561
3.6506
1.66 < 2
S3
NFrom Table11, the method based on wavelet relative energy entropy cannot select the fault line accurately. When the overhead line S2 fault occurs, according to (21), the entropy multiple β can be calculated, C
1
>
C
2, while β
<
2 is established, the selection results are incorrectly, because of the big differences on the distribution parameters between the cable line and the overhead line, which make the healthy cable line’s entropies bigger than the fault overhead line. Moreover, before the wavelet transform, the appropriate mother wavelet function should be selected, while the selection process is difficult. It can be seen that the accuracy of the method based on wavelet relative energy entropy is not guaranteed.


## 9.2. EMD Relative Energy Entropy

The local characteristics of the TZSC can be represented by empirical mode decomposition (EMD); after decomposition, IMF components can be obtained [8], and the waveforms of the mode function IMF1~IMF4 are as shown in Figure 13.EMD decompose the TZSC.
(a)
TZSC
(b)
IMF1
(c)
IMF2
(d)
IMF3
(e)
IMF4In view of IMF1 component belonging the highest frequency component of the TZSC, containing rich transient information, so IMF1 is selected as the characteristic component, and then wavelet denoising and the threshold are used to correct the characteristic component coefficients. Finally, the relative energy entropies of IMF1 can be obtained, and the selection results are as shown in Table12.Table 12
EMD relative energy entropy.
Fault line
θ
R/Ω
EMD relative energy entropyC
j
Entropy multipleβ
Results
Right?
S1
S2
S3
S4
Overhead line S2
0°
20
1.9112
4.5143
5.1857
4.9335
2.36 > 2
S2
Y
200
0.8032
3.8554
5.9764
4.8981
4.80 > 2
S2
Y
150°
20
2.0273
4.4264
5.7486
4.3812
2.18 > 2
S2
Y
200
0.9719
4.1269
6.8780
6.2598
4.25 > 2
S2
Y
Hybrid line S3
0°
20
0.7617
1.1209
4.8099
2.8481
1.47 < 2
S3
Y
200
0.6714
0.9501
4.7789
2.9695
1.12 < 2
S3
Y
150°
20
1.0806
1.5968
5.3013
4.0003
1.48 < 2
S3
Y
200
0.7927
1.3501
5.0642
2.7225
1.70 < 2
S3
Y
Cable line S4
0°
20
0.5502
0.9428
5.4257
4.2543
1.71 < 2
S4
Y
200
0.2029
0.7532
2.0443
3.7942
3.71 > 2
S2
N
150°
20
0.2706
0.3861
4.8885
4.0831
1.43 < 2
S3
N
200
0.3037
1.5814
4.2905
5.5391
5.21 > 2
S2
NFrom Table12, when the cable line S4 fault occurs, after decomposition by EMD, the EMD relative energy entropy of each branch line can be obtained, and then the entropy multiple β can be calculated, β is bigger than the threshold, according to steps 1 and 2, the cable line S4 belongs to the healthy line, and the selection results are incorrect. Moreover, the decomposition process may cause modal aliasing, which can make a big difference on IMF1 component to represent the TZSC. Therefore, it can be seen that the accuracy of the method based on EMD relative energy entropy is not guaranteed.From Tables5, 6, 11, and 12, when the branch line fault occurs, because of the big difference of distribution parameters between the cable line and overhead line, the relative energy entropy of each branch line can be obtained, according to the fault line selection steps 1, 2, and 3, and the accuracy of the method based on wavelet relative energy entropy and EMD relative energy entropy is not guaranteed. Moreover, wavelet transform needs to select the appropriate mother wavelet function, while the selection process is difficult and the modal aliasing may be caused by EMD decomposition process, while SK algorithm can decompose the TZSC in time-frequency domain adaptively, which contributes to the accuracy of the selection results. Therefore, the accuracy of selection method based on SK relative energy entropy is higher than the methods based on the wavelet relative energy entropy and EMD relative energy entropy.


## 10. Conclusions

This paper proposes a stepped fault line selection method based on SK relative energy entropy, and the conclusions are as follows.(1)
SK belongs to higher order statistics with good ability to restrain the noise and can characterize the TZSC in frequency domain as well. The relative homogeneity of SK energy distribution is described quantitatively by the SK relative energy entropy. In addition, utilizing the max SK relative energy entropy as the criterion to select the fault line, the transient universal characteristic quantity is selected in fact. However, the SK is limited by the compromise problem of time-frequency domain, and the type of window function and the length of window should be determined at the beginning.
(2)
The ratios of the first half-wave extremes between neighboring lines and SK relative energy entropy are regarded as fault features to select the fault line. When the bus bar occur fault, the ratios of the first half-wave extreme between neighboring lines are positive. If the entropy multiple is greater than or equal to the threshold, judge the overhead line fault, if not, judge the hybrid line or the cable line fault.Future work should concentrate on how to determine the appropriate auxiliary noise intensity according to the characteristics of TZSC. Moreover, higher adaptability of the method based on SK relative energy entropy should be advanced with the high resistance.



---

*Source: 726205-2014-07-22.xml*