A novel method of fault line selection based on IOS is presented. Firstly, the IOS is established by using math model, which adopted TZSC signal to replace built-in signal of duffing chaotic oscillator by selecting appropriate parameters. Then, each line’s TZSC decomposed by db10 wavelet packet to get CFB with the maximum energy principle, and CFB was solved by IOS. Finally, maximum chaotic distance and average chaotic distance on the phase trajectory are used to judge fault line. Simulation results show that the proposed method can accurately judge fault line and healthy line in strong noisy background. Besides, the nondetection zones of proposed method are elaborated.
1. Introduction
When SPG fault occurs in power distribution network, fault line selection method is always study emphasis. For the arc suppression coil grounding system, the fault current is very weak, and judging the fault line accurately becomes more difficult. Therefore, power distribution network in china is also called small current to ground system. For a long time, many scholars did a lot of studies about fault line selection and came up with many novel methods of different theories, which obtained some achievements [1–7].
The signal processing tools, such as wavelet transform (WT), S-transform, Prony method, and Hough Transform, which used to obtain fault signals feature in many methods. Then, many papers establish fault line selection criterion with artificial neural network [8] or support vector machine [9]. In [1], the paper adopted WT to study SPG fault, it used good localization ability in time and frequency domain to extract the characteristic of transient signal on different scales, but WT is easily influenced by noise. Besides, different wavelet-based functions have different results for characteristic of transient signal. In [2], TZSC data of first T/4 is used, and S-transform is carried out to determine the main characteristic frequency. Therefore, fault line selection result is detected by comparing S-transform energies of each feeder line. To avoid TAs density effect of magnetic saturation, paper [3] puts forward a piecewise Prony method; it improved the precision of the traditional Prony method, and the proposed method was used to match TZSC signals to get fault information. In [4], Hough transform was used to achieve fault line selection; actually, Hough transform has an excellent performance in the case of not exceeding two-dimensional parameter space, but if parameter space increased, calculation amount will rise sharply, meanwhile, costly storage space and time-consumption also will be surged. Paper [8] presents a method based on rough membership neural network to classify 10 kinds of fault types about transmission line. In [9], support vector machine (SVM) has advantages of solving small sample size, nonlinear and high dimensional pattern recognition problems; however, its identifying ability is influenced by self-parameters easily.
Chaotic oscillator detection method has some advantages, such as it can effectively extract weak signals under the noise background, it is sensitive to small periodic signals, and it has certain immunity to noise. Therefore, the application areas of the method are very wide., for example, arcing faults protection [10], fault diagnosis of wind turbine gearbox [11], faults line detection [12], and islanding detection [13].
This paper proposed a novel method of fault line selection based on IOS. Firstly, in the aspect of extracting feature signal, TZSC is decomposed by db10 wavelet packet. Next, according to the maximum energy principle, CFB is selected and is inputted to IOS. Finally, fault line selection criterion based on maximum chaotic distance and average chaotic distance of phase trajectory is presented to judge the fault line.
The remaining of this paper is arranged as follows. In Section 2, we analyzed the small current to ground fault. In Section 3, we indicated duffing chaotic oscillator disadvantage that cannot detect the damped oscillation signal. In Section 4, we introduced the IOS work principle. In Section 5, we proposed the fault line selection method and introduced the wavelet packet theory and detection factor setting. In Section 6, we simulated the small current to ground fault and verified the effectiveness of proposed method. In Section 7, the nondetection zone of proposed method is presented. A brief summary concludes the paper in Section 8.
2. Analysis of Small Current to Ground Fault
When SPG fault occurs in small current to ground system, the current flow distribution is shown in Figure 1.
Current flow of SPG fault.
In Figure 1, the power’s neutral point attached an arc suppression coil, so current flow has a change by the coil, but capacitor voltage amplitude and distribution are both the same as nonsolidly grounded system. However, the difference is that inductive current I·L flow from line to ground point. Therefore, the total current is
(1)I·f=I·L+I·CΣ,
where I·CΣ and I·L were capacitance current to ground and inductance current, respectively. Using L to represent inductance, so I·L=-E·A/jωL.
Phase difference between I·CΣ and I·L is 180°, so I·f will decrease due to arc suppression coil compensation. When the coil runs with over compensation stage, the residual current to ground is inductive. Reactive power flow direction of capacitive is from bus line to fault line, it is same as healthy line. Therefore, the fault information is less.
Transient zero-sequence equivalent circuit of SPG fault is used to analyze TZSC. In Figure 2, C0, L0 are zero-sequence capacitance and zero-sequence inductance, respectively; Rg is transition resistance; Rp, Lp are equivalent resistance and inductance of arc suppression coil respectively; e(t) is zero-sequence voltage.
Transient zero-sequence equivalent circuit.
In Figure 2, TZSC i0.t is [14]
(2)i0.t=i0L.t+i0C.t=ILmcosφe-t/τL+ICm(ωfωsinφsinωtddddddddddddddddii-cosφcosωft(ωfωsinφsinωt)e-δt,
where i0L.t, i0C.t are inductive current and capacitive current of TZSC and their initial values are ILm and ICm, respectively (ICm=UphmωC, ILm=Uphm/ωL), Uphm is phase voltage amplitude, ω is angular frequency of power frequency, ωf, δ are oscillation angular frequency and attenuation coefficient of TZSC, τL is decay time constant of inductive current, and φ is initial phase angel.
When SPG fault occurs, ground point will flow transient capacitive current and transient inductive current; the former has fast attenuation speed and the later has slow speed, and i0C.t determined the amplitude and frequency of i0.t. With power lines length increasing, the self-oscillation frequency decreases, and the free oscillation component amplitude of i0C.t also decreases; meanwhile, free oscillation duration will be reduced to half power frequency cycle. Moreover, the distribution network is a small current to ground system, so fault signal is weak, and it influenced by noise easily. Therefore, TZSC difference between fault line and healthy line is small, and the fault information is less; it was difficult to achieve fault line selection for distribution networks.
In recent years, duffing chaotic oscillator detection method is used widely, because it is sensitive to small periodic signal and has immunity to noise. Therefore, it can be used to detect weak signals in strong noisy background for power distribution networks.
3. Duffing Chaotic Oscillator Disadvantage
The Holmes-type duffing chaotic oscillator has two characteristics.
System motion state is not influenced by large amplitude of irregular noise.
System motion state changed by the periodic signal, its frequency is the same as the system, even if the signal is very weak.
Holmes-type duffing equation is given as follows:
(3)x′′+kx′-x+x3=γf1(t),
where x, k, and f1(t) are output signal, damping ratio, and built-in signal, respectively, γ is amplitude factor of built-in signal, and f1(t) is equal to cos(t) generally.
If γ is equal to 0.8 and 0.82625, respectively, we can get chaotic state and periodic state; it is shown in Figure 3.
Different states of phase trajectories.
Chaotic state
Periodic state
In Figure 3, the system has different states by setting different γ values. Therefore, duffing chaotic oscillator detecting method usually used system’s different states to detect different signals.
Considering the external signal variability, (3) is changed by scale variation. Then, add measured signal, it is shown as follows:
(4)x..+ωkx.-ω2x3+ω2x5=γcos(ωt)+f2(t).
Here, ω is angular frequency of measured signal, f2(t) is external signal.
We choose damped oscillation signal to study, it contains a range of frequency components; it was shown in Figure 4; Figure 5 is obtained by calculating (4).
Damped oscillation signal.
Phase trajectories of different amplitudes.
Original amplitude
Amplitude reduced 10 times
Amplitude reduced 100 times
Amplitude reduced 1000 times
It showed that duffing oscillator detection system is not influenced by mutative amplitude of external signal in Figure 5. The reason is that external signal has different frequencies, but the system has only one frequency. Therefore, the damped oscillation signal cannot make duffing chaotic oscillator turn to other states, in other words, the duffing oscillator detection system has weak detection ability for damped oscillation signal.
4. IOS Work Principle
TZSC signal of small current to ground system is shown in Figure 6.
TZSC signal.
In Figure 6, TZSC signal has oscillation and attenuation states, and (2) shows that different TZSC have different oscillation and attenuation states. Therefore, TZSC math model is proposed. Consider the following:
(5)i(t)=A1cos(2πω0t+φ1)+A2eδ1tcos(2πω0k1t+φ2)+A3eδ2tcos(2πω0k2t+φ3)+A4eδ3t.
Equation (5) consists of four components, where Ai, k1, k2, δ and φj are component amplitude, natural number, decimal number, attenuation coefficient, and initial phase, respectively, i=1,2,3,4, and j=1,2,3. Each component, respectively, represents fundamental component, integer harmonic, noninteger harmonics, and attenuation component.
In order to make up for the disadvantage that cannot detect damped oscillation signal, the built-in signal f1(t) is replaced by (5); therefore, IOS is proposed and shown as follows:
(4)x′=y,y′=-ky+x3-x5iiiiiii+γ(+e-tcos(2π×ω0×1.1×t)+e-t)cos(2π×ω0t)iiiiiiiiiiiiii+e-tcos(2π×ω0×1×t)iiiiiiiiiiiiii+e-tcos(2π×ω0×1.1×t)+e-t)iiiiiii+β×f2(t),
where Ai=1(i=1,2,3,4), ω0=50, k1=1, k2=1.1, δ=-1, φj=0, γ is amplitude factor of built-in signal, and β is detection factor.
4.1. Amplitude Factor Impact
If f2(t)=0, k=0.5, and γ is equal to 10−4, 0.025, 0.03, and 0.035, respectively, the IOS phase trajectories are shown in Figure 7.
Phase trajectories of different γ.
γ=10-4
γ=0.025
γ=0.03
γ=0.035
Comparing Figure 3 with Figure 7, in this paper, the state of Figure 7(a) is named chaotic-like state (CLS). Similarly, Figures 7(b) and 7(c) are named critical state (CS), Figure 7(d) is named periodic-like state (PLS).
In Figure 7(a), when IOS is in CLS, points movement are small on phase trajectory, motion tendency is obscure, and it is integrally in chaotic state. When γ increased to a certain value, the phase trajectory turned to PLS; its movement scale is big and movement tendency is clear; the overall movement state can be thought as circular movement with reduced radius; it is shown in Figure 7(d).
CS is between CLS and PLS, and the overall movement state can be thought as a curve segment keeping away from the origin point; it shows in Figures 7(b) and 7(c).
4.2. Damping Ratio Influence
We used the PLS parameters and set damping ratio k as equal to 0.01, 0.05, 0.5, and 1, respectively. By solving IOS, the phase trajectories are shown in Figure 8.
Phase trajectories of different k.
k=0.001
k=0.05
k=0.5
k=1
According to Figure 8, with damping ratio kincreasing, movement scale of phase trajectory becomes smaller; overall movement state changed from circular movement with decreasing radius to a clear curve segment. However, when k decreased, movement scale became larger and overall movement state is from a clear curve segment turn to a closed curve. In other words, when k increased from small to large, PLS phase trajectory followed: sparse→dense→sparse.
All in all, rational k can make PLS phase trajectory become dense and make the difference between CS and CLS be more obvious; the phenomenon is more beneficial to identifying PLS. But, lower and bigger k would make PLS turn to other states. Therefore, choosing the rational damping ratio is very important.
4.3. IOS Example
To demonstrate IOS detection ability, transient zero-sequence signal from paper [15] is chosen and inputted to IOS (k=0.05, γ=0.002, β=0.001). By changing amplitudes and frequencies of transient zero-sequence signal, Figure 9 is obtained.
Phase trajectories of different signals.
Change frequency
Change amplitude 1
Change amplitude 2
No change
Figure 9(a) results are obtained by the changed frequency signal that comes from paper [15], on the contrary, Figures 9(b) and 9(c) results are obtained by changing the signal amplitude. Similarly, Figure 9(d) results are obtained by importing the signal directly without change.
According to Figure 9, the example signal in the paper can make IOS have different phase trajectories. Therefore, different oscillation damping signals can generate different phase trajectories; it shows that IOS has superior ability to detect oscillation damping signals.
5. Fault Line Selection Method5.1. Wavelet Packet Transform Theory
Wavelet packet transform (WPT) is an improved method of wavelet transform (WT), and it has better time-frequency resolution, therefore, it can effectively extract fault feature of TZSC signal. WPT is described as follows.
Defining Ujn is the closure space of un(t), so Uj2n is closure space of u2n(t), and un(t) meet the following two-scale equation:
(7)u2n(t)=2∑l∈Zh(l)un(2t-l),u2n+1(t)=2∑l∈Zg(l)un(2t-l),
where g(l)=(-1)lh(1-l), h(l) and g(l) have orthogonal relationship; when n=0, (7) can be expressed. Consider the following:
(8)u0(t)=2∑k∈Zh(l)u0(2t-l),u1(t)=2∑k∈Zg(l)u0(2t-l).
When n is equal to Z+, it can generalize Uj+1n=Ujn⊕Uj2n+1, j∈Z, n=Z+, orthogonal wavelet packet is determined by basis function: u0(t)=φ(t). Since φ(t) is determined by hl uniquely, so {un(t)} is called the orthogonal wavelet packet about sequence {hl}. About wavelet packet decomposition, it is the band-pass or low-pass filter actually. The filter bandwidth is [fs(l-1)/2j,fsl/2j], where j is level number of wavelet decompositions, l is the lth contacts of wavelet decomposition, and fs is input signal frequency; WPT decomposition structure is shown in Figure 10.
WPT decomposition structure.
A is low-frequency, D is high-frequency. The decomposition has the following equation. (9)S=AAA3+DAA3+ADA3+DDA3+AAD3+DAD3+ADD3+DDD3.
According to power signal feature, WPT can choose appropriate CFB to match signal spectrum adaptively; it improved time-frequency resolution and had a wide prospect. Because daubechies series wavelet has some merits, such as orthogonal characteristic, compact support, and N-1 vanishing moments; therefore, db10 wavelet packet is chosen to analyze TZSC signal.
5.2. Detection Factor Setting
Each line’s TZSC is decomposed by db10 wavelet packet; then, the CFB is chosen by maximum energy principle, with the change of grounding resistance, the different CFB are obtained; it is shown in Figure 11.
CFB of different resistances.
Small faults resistance
Big faults resistance
In Figure 11, CFB is also damped oscillation signal, and its maximum amplitude decreased with fault resistance increasing. In other words, when fault situation changed, CFB amplitude will be changed. In order to make IOS have a better applicability, the detection factor (β) should be chosen appropriately. Therefore, CFB was input to IOS, and the equation is proposed as follows:
(10)x′=y,y′=-ky+x3-x5y′=+γ(cos(2π×50t)+e-tcos(2π×50×1×t)y′=+γ+e-tcos(2π×50×1.1×t)+e-t)y′=+β×ξk(j)(n),
where ξk(j)(n) is CFB signal, k=0.05, γ=0.002, and β=0.000286. The value ofβis chosen by a large number of experiments [12].
5.3. Fault Line Selection Criterion Based on Chaos Distance
Compared with paper [12], it is difficult to accurately distinguish between CS and PLS only by the phase trajectory; the reason is that their phase trajectories are very similar. Therefore, chaotic distance and average chaotic distance are put forward to distinguish the phase trajectory state.
Chaotic distance (DC):DC is square distance from each point to origin point on phase trajectory; it can represent the size of each point movement scale:
(11)DC=(X-0)2+(Y-0)2.X is the vector which is acquired by solving (4), Y is the first-order differential vector of X, and the maximum of DC is named MDC.
Average chaotic distance (DA):DA is the average of DC; it can represent the overall movement scale:
(12)DA=∑DCm.m is the number of points on phase trajectory.
By using the chaos distance, a range of CFB signals is chosen and inputted into IOS, by solving IOS, we get Figure 12.
Chaotic distances of different states.
Chaotic-like state
Chaotic-like distance
Criticality state
Criticality distance
Period-like state
Period-like distance
In Figures 12(a) and 12(b), CLS scale is very small, and its points always move near the origin point, so the DC of CLS is also very small. The DC of CS increases with time increasing, which complies with curvilinear movement that is away from origin point; it is shown in Figures 12(c) and 12(d). But DC of PLS shows oscillation damping state, which complies circular movement with reducing radius; it is shown in Figures 12(e) and 12(f). Chaotic distance can clearly indicate movement scale and tendency of each state; it avoids error distinction. MDC and DA values of Figure 12 are shown in Table 1.
MDC and DA of Figure 12.
State
MDC
DA
PLC
1.4690
0.0356
CS
1.4670
0.0061
CLS
0.0019
0.0002
MDC value of CLS is very small, and its magnitude order stabilized near z1 (in this paper, z1=10-3), but MDC of CS and PLS are both to stabilize near z2 (in the paper, z2=1.4). For CS and PLS, their DA has a larger difference. In order to distinguish CS and PLS accurately, sj is proposed (sj=MDC-z2). While |sj|<a (in the paper, a=0.1), CS and PLS states cannot be distinguished by MDC solely. Therefore, it proposed hj:
(13)hj=max(DAj)DAj,
where max(DAj) is the maximum value of DA in all phase trajectories, j is the number of phase trajectories.
If hj>b (in the paper b=1.5), the phase trajectory that has max(DAj) is judged PLS, others are judged CS. If hj≤b, all phase trajectories are judged CS.
In addition, phase trajectory of magnitude order near the z1 is judged CLS.
From above analysis, CLS, CS, and PLS can be accurately distinguished by DA and MDC. When SPG fault occurs, CFB with detection factor is inputted to IOS, and it used DA and MDC to select the fault line. Therefore, the novel criterion of fault line selection method is proposed.
Step 1.
If only one line has |sj|<a (this line is called L), and the phase trajectory magnitude of other lines stabilized near z1. In other words, L makes IOS turn to PLS or CS and other lines make IOS turn to CLS. Therefore, L is judged as fault line and other lines are judged as healthy line.
Step 2.
If many lines have |sj|<a, calculated hj, If hj>b (b=1.5), the line with maximum DA makes IOS turn to PLS (this line is called L), other lines make IOS turn to CLS or CS. Therefore, L is judged as fault line and other lines are judged as healthy line.
Step 3.
If all lines do not comply with Steps 1 and 2, the fault is judged as bus fault.
6. Simulation and Verification6.1. Simulation Model
In this paper, SPG fault simulated by ATP-EMTP, and the simulation model is shown in Figure 13. There are 4 lines, such as: S1, S2, S3, and S4; their lengths are 13.5 km, 24 km, 17 km, and 10 km, respectively. Simulation parameters are
Transformer. 110/10.5 kV, leakage impedances of high voltage side: 0.40 + j12.20 Ω, leakage impedances of low voltage side: 0.006 + j0.183 Ω, excitation current: 0.672 A, magnetizing flux: 202.2 Wb, magnetic resistance: 400 kΩ.
Load. All are delta connected, ZL=400+j20Ω.
Arc Suppression Coil. LN=1281.9mH, RN=40.2517Ω.
6.2. SPG Fault Analysis
Set fault resistance and initial angle as equal to 0.0001 Ω, and 0°, respectively; fault line voltage of phases A, B, and C was shown in Figure 14, and each line zero-sequence current was shown in Figure 15.
Phase voltage of fault line.
Zero-sequence current of each line.
It shows that phase voltage of healthy phases B and C are increasing and their amplitude is equal, but the phase voltage of fault phase A is equal to 0 approximately in Figure 14. This phenomenon meets the boundary condition of A phase to ground fault. Furthermore, all lines have same voltage value, because each line is linked with the same bus line in parallel connection. Therefore, this paper chose zero-sequence current to study the fault line selection.
In Figure 15, when SPG fault occurs, zero-sequence current value of each line is not equal to 0, its current waveforms have damped oscillation state, and the attenuation speed is fast. The gap of zero-sequence current between fault line and healthy line is smaller after 0.03 s, but the gap of TZSC between fault line and healthy line is larger in 0.02-0.03 s. In addition, healthy line TZSC is similar. Therefore, TZSC is chosen to study fault line selection method in this paper.
6.3. Changing Initial Phase Angle and Resistance
When fault occurred in S1, fault location is 5 km away from bus line, and the initial angles are 0°, 30°, 60°, and 90° respectively. Therefore, each TZSC of T/2 decomposed by db10 wavelet packet when fault occurred. According to paper [16], IOS parameters were set as k=0.05, γ=0.002, and β=0.000286.
To show phase trajectory difference, this paper gives 2 examples, one is that initial angle is equal to 30° and fault resistance is equal to 1100 Ω, to express conveniently, this paper uses (30°, 1100 Ω) to represent the fault condition; the other example is (90°, 2000 Ω). The TZSC and CFB of the two examples were shown in Figures 16 and 17, respectively.
No 1 of TZSC and CFB.
Transient zero-sequence current
Characteristic frequency band
No 2 of TZSC and CFB.
Transient zero-sequence current
Characteristic frequency band
It shows that TZSC and CFB are both damped oscillation signal in Figures 16 and 17, different fault conditions have different oscillation damping degrees. Furthermore, CFB can describe time-frequency characteristics of transient signal more clearly in some frequency bands; it provides a basis for fault feature information processing [17]. Therefore, CFB is analyzed by IOS for fault line selection. According to new criterion, fault line selection result is shown in Table 2.
Fault line selection results of different resistances and angles.
Fault line
Fault situation
Distance of phase trajectory
S1
S2
S3
S4
Selection result
S1
(30°, 90 Ω)
MDC
1.4666
0.0005
0.0006
0.0017
S1
(0°, 400 Ω)
MDC
1.4691
1.4718
0.0005
0.0029
S1
DA
0.0598
0.0072
0.0001
0.0004
(0°, 2000 Ω)
MDC
1.4689
0.0013
0.0005
0.0007
S1
(30°, 400 Ω)
MDC
1.4679
0.0004
0.0019
0.0053
S1
(30°,1100Ω)
MDC
1.4690
1.4670
0.0014
0.0019
S1
DA
0.0356
0.0061
0.0002
0.0003
(30°, 2000 Ω)
MDC
1.4701
0.0035
0.0005
0.0012
S1
(60°, 90 Ω)
MDC
1.4667
0.0006
0.0009
0.0060
S1
(60°, 400 Ω)
MDC
1.4663
1.4668
0.3589
0.0043
S1
DA
0.0521
0.0088
0.0030
0.0007
(60°, 2000 Ω)
MDC
1.4677
0.0024
0.0005
0.0006
S1
(90°, 400 Ω)
MDC
1.4650
0.0014
0.0005
0.0019
S1
(90°, 2000 Ω)
MDC
1.4674
0.0004
0.0015
0.0006
S1
Under the fault condition: (30°, 1100 Ω), according to Table 2, MDC magnitude orders of S3 and S4 are both 10-3; it judges S3, S4 as healthy line. sj of S1 and sjof S2 are both less than 0.1; it indicated that S1, S2 are PLS or CS, combining DA, hj=0.0356/0.0061=5.84>1.5; therefore, S1 makes phase trajectory turn to PLS, and S2 makes phase trajectory turn to CS, which is same as Figure 18. Then, S1 is judged as fault line; other lines are judged as healthy line.
Phase trajectory of (30°, 1100 Ω).
S1
S2
S3
S4
Under the fault condition: (90°, 2000 Ω), the sj of S1 meets sj=1.4674-1.4=0.0674, because sj is smaller than 0.1 and the MDC magnitude orders of other lines are all 10−3, it judges that S1 makes IOS turn to PLS and other lines make IOS turn to CLS. Therefore, the S1 is judged as fault line, and other lines are judged as healthy line. The selection result is the same as the simulation setting. It is as shown in Figure 19.
Phase trajectory of (90°, 2000 Ω).
S1
S2
S3
S4
From above analysis, fault line selection results are same as simulation settings, which indicate that the method of this paper can adapt to different resistances and different phase angles.
6.4. Different Lines Fault
Refer to paper [18]; S3 is hybrid line and S4 is pure cable line, therefore, simulate the different lines occur fault, and the results are shown in Table 3. With cable lines taking part in simulation, TZSC decay process becomes shorter [16], but it does not affect fault line selection results. Table 3 shows that the method has better results when different lines fault occurred.
Different lines fault result.
Fault line
Fault situation
Distance of phase trajectory
S1
S2
S3
S4
Selection result
S4
(0°, 200 Ω)
MDC
0.0009
1.4187
1.4681
1.4678
S4
DA
0.0006
0.0048
0.0090
0.0863
(0°, 700 Ω)
MDC
0.0006
0.0845
0.0090
1.4684
S4
(0°, 2000 Ω)
MDC
0.0007
0.0004
0.0004
1.4658
S4
S3
(0°, 200 Ω)
MDC
0.0007
0.0777
1.4648
0.1994
S3
(0°, 700 Ω)
MDC
0.0016
0.0009
1.4692
0.0027
S3
(0°, 2000 Ω)
MDC
0.0003
1.4664
1.4691
0.1934
S3
DA
0.0001
0.0069
0.0126
0.0028
6.5. Fault with Strong Noise
Considering actual fault situation, fault signal usually has noise; therefore, add 0.5 db random Gaussian white noise to TZSC to validate the method’s noise immunity. The waveform was shown in Figure 20, it shows that the TZSC signal is covered by noise completely. Fault line selection results are shown in Table 4.
Strong noise fault results.
Fault line
Fault situation
Distance of phase trajectory
S1
S2
S3
S4
Selection result
S1
(30°, 90 Ω)
MDC
1.4684
0.0002
0.0005
0.0006
S1
(0°, 200 Ω)
MDC
1.4667
0.0004
0.0003
0.0018
S1
(0°, 2000 Ω)
MDC
1.4689
0.0013
0.0005
0.0007
S1
S3
(0°, 200 Ω)
MDC
0.0007
0.0777
1.4648
0.1994
S3
(0°, 2000 Ω)
MDC
0.0003
1.4664
1.4691
0.1934
S3
DA
0.0001
0.0069
0.0126
0.0028
S4
(0°, 200 Ω)
MDC
0.0009
1.4187
1.4681
1.4678
S4
DA
0.0006
0.0048
0.0090
0.0863
(0°, 2000 Ω)
MDC
0.0007
0.0004
0.0004
1.4658
S4
TZSC with 0.5 db noise.
Under various fault conditions, the proposed method can also achieve better result, and it indicates that the method has good noise immunity. Furthermore, the fault line selection result is better than paper [19] which adds 20 db noise.
6.6. Different Distances Fault
In fact, fault location has different distances from bus line when fault occurs. In Section 6.6, simulate S1 line with different distance fault, and the fault resistance is 300 Ω, fault distances are 5 km, 30 km, 35 km, 40 km, and 45 km, respectively. Fault line selection results are shown in Table 5.
Different distance fault results.
Fault line
Fault situation
Distance of phase trajectory
S1
S2
S3
S4
Selection result
S1
(0°, 300 Ω, 5 km)
MDC
0.0009
1.4187
1.4681
1.4678
S1
(0°, 300 Ω, 10 km)
MDC
1.4690
0.0004
1.4665
0.0022
S1
DA
0.0583
0.0002
0.0061
0.0004
(0°, 300 Ω, 30 km)
MDC
0.3808
0.0007
0.0288
0.0005
S1
(0°, 300 Ω, 35 km)
MDC
1.4692
0.0003
1.3629
0.0004
S1
(0°, 300 Ω, 40 km)
MDC
1.4687
0.0011
1.4687
0.0003
S1
DA
0.0274
0.0008
0.0154
0.0005
According to Table 5, the result is same as simulation setting, which indicates that the proposed method adapts to different distance faults.
7. Nondetection Zone of Proposed Method7.1. Small Resistance to Ground Fault
Set fault resistance between 0 Ω and 80 Ω. According to proposed method, partial results of fault line selection are shown in Table 6.
Small resistance to ground fault result.
Fault line
Fault situation
Distance of phase trajectory
S1
S2
S3
S4
Selection result
S1
(0°, 5 Ω)
MDC
1.4675
1.4643
1.4663
1.4677
Failure
DA
0.0127
0.0126
0.0379
0.0366
(0°, 10 Ω)
MDC
0.0571
0.0016
1.4677
1.4704
Failure
DA
0.0024
0.0010
0.0393
0.0359
(0°, 30 Ω)
MDC
0.0006
0.0005
1.4654
0.0004
Failure
(0°, 60 Ω)
MDC
0.0003
0.0007
1.4644
0.0030
Failure
(0°, 80 Ω)
MDC
0.0024
0.0003
0.0008
0.0021
Failure
The proposed method in this paper does not work for small resistance to ground fault in Table 6. In addition, when fault resistance value is bigger than 10 Ω, cable line and cable hybrids line usually make IOS turn to PLS, and overhead line makes IOS turn to CLS, but when fault resistance value is smaller than 10 Ω, overhead line can also make IOS turn to PLS; for example, the fault condition is (0°, 5 Ω). To express result obviously, phase trajectory of (0°, 10 Ω) is given in Figure 21.
Phase trajectory of (0°, 10 Ω).
S1
S2
S3
S4
It shows that S1 makes IOS turn to CS, S2 makes IOS turn to CLS, S3, and S4 make IOS turn to PLS in Figure 21; it contradicts the proposed criterion in this paper (Steps 1, 2, and 3). The reason is that fault signal amplitudes are large and damped oscillation states are obvious in this fault situation, it makes all total driving force of IOS too large. There are two conditions, one is that fault line makes IOS directly cross PLS and turn to CLS, the other is that healthy line makes IOS directly cross CLS and turn to PLS or remain CLS. When fault resistance continues to decrease, fault line and healthy line can both make IOS turn to PLS, such as, fault condition: (0°, 5 Ω). In a word, the most essential thing is that the detection factor is determined by experiment, and it is not suitable for small resistance to ground fault.
7.2. High Resistance to Ground Fault at the End of Line
The length of S1 is 13.5 km, simulate fault location at the end of S1, and fault resistance is set to 1.5 kΩ to 2 kΩ. Its result is shown in Table 7.
High resistance to ground fault result.
Fault line
Fault situation
Distance of phase trajectory
S1
S2
S3
S4
Selection result
S1
(0°, 1500 Ω)
MDC
0.0546
1.4653
1.4695
1.4687
Failure
DA
0.0027
0.0076
0.0362
0.0376
(0°, 1700 Ω)
MDC
0.1191
0.0052
1.4749
1.4676
Failure
DA
0.0025
0.0017
0.0342
0.0548
(0°, 1800 Ω)
MDC
0.0014
1.4679
1.4679
1.4672
Failure
DA
0.0001
0.0010
0.0175
0.0402
(0°, 1900 Ω)
MDC
0.0114
1.4658
1.4704
1.4675
Failure
DA
0.0017
0.0060
0.0241
0.0271
(0°, 2000 Ω)
MDC
0.0569
0.0007
0.0550
1.4669
Failure
According to Table 7, when each line CFB was solved by IOS, the healthy line usually makes IOS turn to PLS or CS, and fault line makes IOS turn to CLS, the phase trajectory of (0°, 1900 Ω) is shown in Figure 22.
Phase trajectory of (0°, 1900 Ω).
S1
S2
S3
S4
It is indicated that S1 makes IOS turn to CLS, S2 makes IOS turn to CS, and S3 and S4 make IOS turn to PLS in Figure 22. This result contradicts the proposed criterion in this paper (Steps 1, 2, and 3). This condition occurs due to the fact that fault location is in the end of line, and the fault resistance is larger. It makes voltage and current decrease and the fault current signal become weaker. Therefore, the proposed method is not suitable for this fault condition.
In conclusion, the proposed method is not suitable for small resistance or high resistance to ground fault in this paper; in other words, detection factor has a range, so how to choose better detection factor to adapt extensive fault conditions is to be further studied.
8. Conclusions
SPG fault line selection method based on IOS is proposed in this paper, and the conclusions are as follows.
(1) IOS has a better recognition for TZSC signal, to a certain extent; it makes up for deficiency of the duffing chaotic oscillator which did not detect oscillation damping signal. In addition, TZSC signal decomposed by db10 wavelet packet, it decreased CFB quantities, and improved IOS solving speed.
(2) PLS phase trajectory can be seen as a circular movement with reduced radius; CS phase trajectory can be seen as a curve segment far away from origin point, but CLS phase trajectory is generally seen as chaotic state of movement. In addition, the MDC of CS is similar to PLS, and their value difference usually is about 0.1. But the DA between CS and PLS has a significant gap which usually is equal to 1.5 for their ratio. Therefore, phase trajectory state is usually judged by DA, and MDC of CLS is very small. Besides, magnitude order of MDC is to stabilize at 10−3 for CLS.
(3) There are some nondetection zones of the proposed method: detection factor is determined by experiment, it can cause some errors, and the method has problem of insufficient sensitivity when small and high resistance to ground occur fault.
NotationsSPG:
Single phase to ground
TZSC:
Transient zero-sequence current
CFB:
Characteristic frequency band
SVM:
Support vector machine
IOS:
Improved oscillator system
WT:
Wavelet transform
WPT:
Wavelet packet transform
CLS:
Chaotic-like state
CS:
Critical state
PLS:
Periodic-like state
DC:
Chaotic distance
MDC:
Maximum of DC
DA:
Average of chaotic distance.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this article.
Acknowledgments
This work is supported by Science and Technology Research (12B470002, 14A470004) and Control Engineering Lab Project (KG2011-15) of Henan Province, China, and Youth Foundation (Q2012-28, Q2012-43A) of Henan Polytechnic University, China.
DongX. Z.ShiS. X.Identifying single-phase-to-ground fault feeder in neutral noneffectively grounded distribution system using wavelet transform20082341829183710.1109/TPWRD.2008.9179242-s2.0-54049149895ShuH.PengS.A fault line detection algorithm for distribution network of overhead line and underground cable mixed lines using S-transform energy from short window data200924101521592-s2.0-70749090977WangX. W.WuJ. W.LiR. Y.A novel method of fault selection based on voting mechanism of Prony relative entropy theory20134615965ShuH. C.GaoL.DuanR. M.CaoP.ZhangB.A novel hough transform approach of fault line selection in distribution networks using total zero-sequence current2013379172-s2.0-84878472146SunX.GaoM.LiuD.HouX.YuanR.Analysis and processing of electrical fault signals with modified adaptive prony method2010302880872-s2.0-77958592693NayakP. K.PradhanA. K.BajpaiP.A fault detection technique for the series-compensated line during power swing201328271472210.1109/TPWRD.2012.22318862-s2.0-84875665666ZhangZ. X.LiuX.ZailinP.Fault line detection in neutral point ineffectively grounding power system based on phase-locked loop201482273280LinS.HeZ.ZangT.QianQ.Novel approach of fault type classification in transmission lines based on rough membership neural networks2010302872792-s2.0-77958517040ShuZ.HeZ. Y.WangQ.LinS.Fault line selection of resonant grounding system based on the characteristics of charge-voltage in the transient zero sequence and support vector machine2013411271782-s2.0-84879867833YangJ. G.ZhangR. C.FangH. Y.Application of duffing oscillator signal detection to arcing faults protection2006302469722-s2.0-33846804551SunZ. Q.ChenC. Z.GuY. L.LiuH.Incipient fault diagnosis of large scale wind turbine gearbox based on chaos theory and sampling integral technology20133291131172-s2.0-84879178242ZhangS.ZhaiX.DongX.LiL.TangB.Application of EMD and Duffing oscillator to fault line detection in un-effectively grounded system201333101611672-s2.0-84877676520VahediH.GharehpetianG. B.KarrariM.Application of duffing oscillators for passive islanding detection of inverter-based distributed generation units20122741973198310.1109/TPWRD.2012.22122512-s2.0-84866669535ShuH. C.2008Beijing, ChinaChina Machine PressKangX.LiuX.SuonanJ.MaC.WangC.YangL.New method for fault line selection in non-solidly grounded system based on matrix pencil method20123612889310.3969/j.issn.1000-1026.2012.12.0162-s2.0-84865059083ZhangH.HeZ.ZhangJ.Frequency spectrum characteristic analysis of single-phase grounding fault in resonant grounded systems2012366798410.3969/j.issn.1000-1026.2012.06.0142-s2.0-84859488064ShuH. C.ZhuM. M.Faulty line selection based on time-frequency characteristics of transient zero-sequence current201333916KangZ.LiD.LiuX.Faulty line selection with non-power frequency transient components of distribution network2011314162-s2.0-79955974337GuoM. F.LiuS. D.YangG. J.A new approach to detect fault line in resonant earthed system based on Hilbert spectrum band-pass filter and transient waveform recognition20133236774