Abstract

Since weak fault signals often lead to the misjudgment and other problems for faulty line selection in small current to ground system, this paper proposes a novel faulty line selection method based on variable scale bistable system (VSBS). Firstly, VSBS is adopted to analyze the transient zero-sequence current (TZSC) with different frequency variety scale ratio and noise intensity, and the results show that VSBS can effectively extract the variation trends of initial stage of TZSC. Secondly, TZSC is input to VSBS for calculation with Runge-Kutta equations, and the output signal is chosen as the characteristic currents. Lastly, correlation coefficients of every line characteristic current are used as the index to a novel faulty line selection criterion. A large number of simulation experiments prove that the proposed method can accurately select the faulty line and extract weak fault signals in the environment with strong noise.

1. Introduction

As an important part of the power system, distribution network is closely associated with its users and also has direct impact on the users. Data show that 80% of fault occurring in distribution network is single phase-to-ground fault. When single phase-to-ground fault occurs, the line voltage value is still symmetrical, the fault current is weak, and it could run 1 to 2 hours after fault occurs, which significantly improves the reliability of power supply. However, during the single phase-to-ground fault period, nonfault phase voltage could rise, which will threaten the system insulation and result in interphase shortage, protection tripping, power supply outage, and other problems. Because of the weak fault signal and the harsh working condition, faulty line selection becomes difficult. Therefore, it is necessary to carry out further research in this area [1, 2].

At present, scholars have put forward various faulty line selection methods. Based on different characteristic components, faulty line selection methods for single phase-to-ground could be divided into 3 categories, that is, signal injection method [3], steady-state component method [4], and transient component method [5, 6]. The signal injection method needs additional signal device and its engineering realization is complex. In steady-state component method the characteristic signal is weak, which makes the result unreliable, while the transient characteristics method is more reliable and applicable because the transient characteristics component is larger than steady component and it will not be influenced by the arc suppression coil and it will need no additional devices [711]. Papers [7, 8] use wavelet transform to extract characteristic information for faulty line selection, but the wavelet transform is easily influenced by noise, and the chosen characteristic frequency band may be nonvalid transient faulty component. In addition, different wavelet basis function would lead to different extraction results and thus lead to error judgment. Paper [9] adopts Prony algorithm to fit TZSC signal when fault occurs. This method not only effectively avoids the effect of current transformer saturation flux density on collected signals but also improves the overall Prony fitting precision to a certain extent; but its calculation amount is large, its fitting order is difficult to determine, and the antinoise ability is not strong. In [10], the support vector machine shows its advantages in solving the problem of small sample, nonlinear, and high dimensional pattern recognition, but the recognition ability is easily influenced by its own parameters. Paper [11] uses Empirical Mode Decomposition (EMD) of TZSC to extract the five harmonic components in characteristic components and input them into Duffing oscillator to achieve faulty line selection according to the change of phase diagram. But when TZSC is greatly interfered with, modal aliasing phenomenon of EMD would arise and cause error judgment. Paper [12] employs the ratios of the first half-wave extreme and Spectral Kurtosis relative energy entropy from TZSC to build the stepped faulty line selection method. Paper [13] uses the S-transform to obtain the modulus and phase information of electric components at each frequency point, and this information is employed to detect the faulty line. In [14], Hilbert-Huang transform 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. Paper [15] adopts evidence uncertainty reasoning and compared abnormal events to reduce computation amount and to improve the accuracy of faulty line selection. Paper [16] employs cross-correlation theory to calculate the integrated correlation coefficient of pure fault component of zero-modulus current for each line and takes the line with the smallest one as the faulty line. Paper [17] carries out the wavelet transform to decompose the transient zero-sequence current for each line, calculates the high and low frequency wavelet energies according to the wavelet coefficients, and selects the faulty line according to the maximum value of high or low frequency energy; however, in the strong noise background, the waveform and energy of weak TZSC will be affected.

In recent years, the research on stochastic resonance has made great progress. Stochastic resonance is a new practical technology which uses stochastic resonance principle to detect weak signal, and its research and application have spread into physical fields [18, 19], signal processing [20, 21], mechanical fault diagnosis [22], biology [23], neural network [24], and other academic fields; however, the research on this technology in power system is still needed. Therefore, with detailed study of the effect of TZSC on bistable system, this paper proposes a novel faulty line selection method for small current to ground system based on stochastic resonance theory. For signal feature extraction, the method employs VSBS to deal with TZSC and, then, choose the initial stage of output signal as characteristic current; for faulty line selection criterion, a novel faulty line selection criterion, which is based on cross-correlation coefficient sign, is proposed through calculating correlation coefficient of characteristic signal.

2. Characteristic Analysis of Single Phase-to-Ground Fault

The structure of small current to ground system is shown in Figure 1(a); when it experiences single phase-to-ground, the TZSC analysis circuit of faulty line is shown in Figure 1(b). In Figure 1, and are zero-sequence capacitance and inductance, respectively, is transition resistance of ground point, and are, respectively, equivalent resistance and inductance of arc suppression coil, and is zero-sequence voltage.

When distribution network fault occurs, from Figure 1(b), the TZSC flowing through the fault point is shown as [25]where and are inductive current and capacitive current of TZSC, and its initial values are and , respectively (, ), is phase voltage amplitude, is angular frequency of power frequency, and are oscillation angular frequency and attenuation coefficient of TZSC, is decay time constant of inductive current, and is initial phase angel.

From (1), when single phase-to-ground fault occurs in distribution network, the transient capacitance current has the characteristic of periodic attenuation oscillation. And [1] indicates that the free oscillation frequency of overhead line is within 300 Hz to 1500 Hz and the free oscillation frequency of cable lines is 1500 Hz~3000 Hz.

In addition, studies show that when single phase-to-ground fault occurs, the traveling wave pole is consistent with the overall changing trend of initial stage TZSC in transient process, so the mutation direction characteristic of initial stage TZSC can be used to replace the traveling wave polarity characteristic of TZSC, which can greatly reduce the hardware requirements and improve the reliability of faulty line selection [26]. Besides, whether it is big initial fault angle or small initial fault angle, the whole changing trend of faulty line is opposite to that of nonfaulty line in TZSC initial stage. But the introduction of arc suppression coil will greatly reduce ground fault current of distribution network, and when the fault occurs in voltage zero position, the changing trend of the initial stage TZSC is not easy to distinguish, which will make faulty line selection more difficult.

Hence, in some faulty conditions, it can be seen from the above analysis that TZSC of distribution network belongs to weak signal. As stochastic resonance (SR) theory has the unique advantage of amplification and detection of weak signals [27], it is helpful to employ the stochastic resonance theory to detect TZSC which is used to select the faulty line.

3. Signal-Detecting Ability of Variable Scale Bistable System

The bistable system for studying stochastic resonance is shown in [28]where is time, is input signal, is noise whose intensity is , is output signal, and is the speed of Brownian particle.

According to Fokker-Planck equation, the probability distribution function of is shown in (3) when is and and are equal to and 1, where and are the amplitude and frequency of periodic signal:

Since (3) has nonautonomous , it has no steady-state solution; that is, it can not have exact expression. However, in the adiabatic approximation condition with , , and , the output signal-to-noise ratio of bistable system can be obtained and shown in

Supposing is 1 and is 0.2, the curve of is shown in Figure 2 when changed. The feature of Figure 2 is that, with the rise of , presents a trend of increasing and begins to decrease when reached 0.13, which is the feature of stochastic resonance. So the bistable system can use noise to increase of signal; that is, the weak signal is amplified and detected.

Under the small parameters of adiabatic approximation condition, the theoretical analysis of stochastic resonance of bistable system coincides with the numerical simulation of the bistable system [28]. However, it is improper to apply the method of small parameters stochastic resonance directly to the processing of signals with large parameters. Reference [27] introduces variable scale stochastic resonance to the process and gets better results; however, when the signal is TZSC, what change will happen to stochastic resonance feature of bistable system, and what rules can we get? This section will focus on the influence of each parameter on bistable system and try to figure out the VSBS characteristics under the effect of TZSC.

3.1. Variable Metric Algorithm and Its Evaluation Index

The principle of variable metric algorithm is to transform high frequency into low frequency in order to make the large parameter signal close to or meeting small parameters conditions of stochastic resonance, which means that the frequency is compressed and then detected by bistable system.

The Calculation Process of Variable Scale. According to the frequency and sampling frequency of signals, a frequency compression-scale ratio (CR) is determined, based on which the compression sampling frequency is defined (). Then numerical calculation step () is obtained from , and, finally, the response output of bistable system is numerically calculated.

Since TZSC are generally broadband signals and their frequency range is not confined to one or a small number of frequencies but distributed in a wide frequency band, to which the traditional signal-to-noise ratio measurement cannot be effectively applied, it is necessary to develop other measurement indexes.

Reference [27] shows that although nonlinear Langevin equation cannot accurately predict the motion of Brownian particles, it can well predict the statistical properties of the particle orbits. So, this paper uses cross-correlation coefficient as a measure to describe the response of VSBS whose input signal is weak aperiodic signal. The covariance and cross-correlation coefficients of two signals are shown in

Additionally, in the initial faulty stage, the overall changing trend of TZSC of faulty line is opposite to that of nonfaulty line, so this paper focuses on the changing trend of input signal and output signal.

3.2. Simulation of Variable Scale Bistable System
3.2.1. Simulation of Nonintroducing Variable Scale

Supposing there is a set of measured signals, the sampling points are 500, the corresponding parameters of (2) are ,  A,  Hz, and  db, respectively, and the value of sampling frequency is 1000 Hz. Fourth-order Runge-Kutta algorithm is adopted to calculate (2). And the value of cross-correlation coefficients of and is −0.0078, whose results are shown in Figure 3. Figure 3(a) shows the result of with noise intensity as 100 db, Figure 3(b) shows without noise, and Figure 3(c) shows without being solved by variable scale.

It can be known from Figure 3 and that when both weak signal frequencies and are large parameters (larger than 1), the output and input of the system differ dramatically, and the information contained in the output signal will not be able to represent the original signal. That is why the stochastic resonance method with small parameters can not be directly applied to large parameters signal, so the detection is ineffective.

3.2.2. Simulation of Introducing Variable Scale

Bring in variable scale thought, choose CR as 100, ,  A,  Hz,  db, and the value of sampling frequency is 1000 Hz; calculate (2) with fourth-order Runge-Kutta algorithm and cross-correlation coefficients of and , will be obtained, and its value is 0.8088. The result is shown in Figure 3(d).

Figure 3(d) shows that after the treatment of VSBS, the waveform of output signal becomes orderly. Compared to Figure 3(c), the cross-correlation coefficients between and have obviously improved as well as the amplitude value of ; besides, and belong to strong correlation. Therefore, through frequency conversion, the disorganized large parameter signal is made clear and orderly; besides, is equal to 0.8088, which indicates that could better represent changing trend of submerged in noise, achieving the large parameter stochastic resonance or, exactly speaking, a kind of stochastic resonance.

3.3. Transient Zero-Sequence Current Detection

In order to test whether VSBS can detect the TZSC, the ideal TZSC [29] is defined as below:It can be seen that , which consists of 5 signals, has the characteristic of multifrequency and attenuation; therefore it is a nonperiodic signal. Input it into (2), and its corresponding parameters are and  db and sampling frequency is 100000 Hz. CR is equal to 1000, and then the results of its numerical simulation are shown in Figure 4.

Definition of Characteristic Current . Characteristic current is the output signal obtained by solving VSBS with TZSC by fourth-order Runge-Kutta algorithm. Choose nonnoises and to calculate cross-correlation coefficient , and the value is 0.7628.

When only the first 0.01 s of and is chosen, as shown in Figure 4(a), the noise causes strong disturbance in the initial stage of , which makes the changing trend not so clear as the original signal. It is known from Figure 4(c) that, after VSBS treatment, the changing trend of is similar to that of ; then, their is calculated, and the value has improved to 0.8909. Therefore, VSBS can effectively extract TZSC changing trend of the initial stage.

This method can be used to better extract the change trend of TZSC in the initial stage. This paper defines 0~0.01 s as the initial stage of TZSC, 0.01 s~∞ as noninitial stage, and signal length as the whole stage. To put it vividly, TZSC from (6) is chosen as the label, and the results are shown in Figure 4(d).

3.4. The Detection Adaptability of TZSC

In order to test detection adaptability of VSBS for TZSC, the paper will analyze frequency compression-scale ratio, noise intensity, the initial value, and signal amplitude, respectively.

3.4.1. Frequency Compression-Scale Ratio (CR)

Set and  db and sampling frequency equals 100000 Hz; set CR as 10, 100, 1000, and 5000, respectively, and the change of is shown in Table 1.

It can be seen from Table 1 that, with the increase of CR, between and first increases and then decreases; the reason is that the increase of CR can gradually compress the frequency band range of into VSBS’s detection range, and there may be a most suitable CR making the input and output most relevant, but when CR continues to increase, excessive frequency compression will also lead to decrease of gap between different frequencies of , showing the reduction of frequency species, which will further weaken the detection ability of VSBS. In addition, the calculation of cross-correlation coefficients of and in initial stage shows that has greatly improved, and the changing trend of is the same as that of , which verifies that VSBS can effectively extract the changing trend of TZSC in initial stage.

3.4.2. Noise Intensity (D)

Set , sampling frequency  Hz, and CR = 1000. Set noise intensity as 0 db, 50 db, 100 db, 1000 db, and 5000 db, respectively, and the change of cross-correlation coefficient is shown in Table 1. The change of characteristic current is shown in Figure 5(a) when the noise intensity of is 0 db, and the change of characteristic current is shown in Figure 5(b) when the noise intensity of is 100 db.

It can be seen from waveform in Figure 5 that the amplitude of increases with the increase of D, which means that part of the noise energy is transferred to [27]. In addition, the increases of made the latter part waveform disorderly; however, the changing trend of the waveform in initial stage is clear which shows no difference with the changing trend of nonnoise; therefore, this once again shows that the VSBS can extract initial stage TZSC. From in Table 1 we know that, within a certain range of noise intensity, the increase of shows little effect on of whole stage and of initial stage, which indicates that VSBS can well extract changing trend of TZSC in initial stage with the disturbance of strong noise. However, excessive noise intensity will produce a wide range of interference frequency components, which will affect the existence of the original signal, resulting in the decrease of cross-correlation coefficient and fuzziness of the changing trend in initial stage. It is worth noting that when the noise intensity is 0, the VSBS can also predict the changing trend of TZSC.

3.4.3. Signal Initial Value (IV)

Set ,  Hz, CR = 1000, and  db; set initial value of as 0 and 34.8, respectively, where 34.8 is initial value (IV) of . Then carry out numerical simulation and the change of characteristic current and cross-correlation coefficient is shown in Table 1 and Figure 5. Figure 5(c) is characteristic current when initial value of is 0, and Figure 5(d) is characteristic current when initial value of is 34.8.

It can be known from Figures 5(c) and 5(d) and Table 1 that when initial value of is 0, the change trend of is closest to that of , especially in the initial stage. Either increase or decrease of the initial value will decrease , because when initial value of VSBS is 0, any tiny disturbance is likely to cause it to move in the double potential well with large amplitude, so it can better reflect the moving trend of signals. However, when the initial value is too large, tiny disturbance may not be enough to cause a large amplitude motion in the double well potential or only small range of motion in a single potential well, therefore, it will weaken the detection ability of VSBS, and this is consistent with the decrease of in Table 1. Another reason for adopting TZSC to select faulty line in this paper is that when the initial value of TZSC is 0, VSBS and faulty line selection can be better combined [30].

3.4.4. Signal Amplitude

Set ,  Hz, CR = 1000, and  db; set the initial value as 0, increase the amplitude of , and set the amplitude factor as 1/100, 1/10, 1, 10, and 100, respectively.

It is known from Table 1 that, with the increase of amplitude , cross-correlation coefficient of whole stage and initial stage first increases and then decreases, but of noninitial stage increases. The reason is that belongs to damped oscillation signal, and, compared with that of noninitial stage, the amplitude values of initial stage are always larger, and the detection ability of VSBS on in noninitial stage is weaker. Therefore, appropriate increase of signal value would improve detection ability of VSBS.

Based on the above analysis, the features of BSVS detecting TZSC are summarized as follows:(1)Appropriate frequency compression ratio can improve the signal detection performance.(2)For small amplitude signal, appropriate increase of amplitude can improve signal detection performance.(3)For the signal with zero initial value, VSBS has a better detection performance of changing trend in its initial stage.(4)The cross-correlation coefficient of whole stage is always smaller than that of initial stage.

4. Faulty Line Selection Method

4.1. Parameter Setting

Based on the above analysis, this paper will select faulty line according to the following characteristics of VSBS detecting TZSC:The overall changing trend of TZSC in initial stage between faulty line and nonfaulty line is opposite.VSBS has excellent detection ability for the changing trend of TZSC in initial stage.When single phase-to-ground fault occurs in small current to ground system, the free oscillation frequency of overhead lines is generally 300 Hz~1500 Hz, while the free oscillation frequency of cable lines is 1500 Hz~3000 Hz. In addition, different fault conditions may cause the TZSC spectrum to transfer into low frequency band [25].Appropriate increase of signal amplitude helps to improve detection performance of VSBS.

According to ① and ②, this paper will focus on cross-correlation coefficient of different lines. Since the TZSC before failure is 0, when calculating cross-correlation coefficients, we will choose T/4 cycle after fault as the initial stage (0.02 s~0.025 s) in the paper; with ③ frequency varieties are compressed as much as possible to make frequency varieties into the frequency range which VSBS can detect, in order to enhance the adaptability of the method; therefore, the frequency compression-scale ratio (CR) is set as 1500. Based on ④ and simulation experiment, when the maximum amplitude of signal is less than 5, we first expand the amplitude by 10 times and then input it into VSBS. In addition, we find that TZSC amplitude before fault is not 0 but a very small value, which needs to be set to 0 before fault.

4.2. Pretreatment of Faulty Line Selection

Take line as an example:Choose TZSC of line number 1 from one cycle before fault to one cycle after fault as TZSC , and set the signal one cycle before fault as 0.Judge whether the maximum amplitude of is smaller than 5, if it is, carry out ③, and if it is not, carry out ④.Expand amplitude of by 10 times and input it into VSBS to calculate characteristic current .Directly input into VSBS to calculate characteristic current .

Calculate of all lines according to the above steps.

4.3. Steps of Faulty Line Selection

Take as an example to explain the steps of faulty line selection.

Step 1. Calculate cross-correlation coefficients between line number () and other lines () in initial stage.

Step 2. Count positive and negative signs of calculated :(1)If all the signs of are the same, the output is “” and line is judged as faulty line.(2)If all the signs of are not the same, the output is “1” and line is judged as nonfaulty line.

The remaining lines are also judged with the same steps.

5. Case Study

5.1. Simulation Model

ATP-EMTP is used in this paper to simulate the single phase-to-ground fault. The simulation model is shown in Figure 6 and the parameters of simulation model are the same as [31].

5.2. Simulation Results with Changing Phase and Resistance

Build the simulation model according to the parameters, make fault of line occur at the point 5 km from the bus, and change the initial fault angle (0°, 30°, 60°, and 90°) and ground resistance for simulation. It is known from [25] that when single phase-to-ground fault of the small current to ground system occurs, the fault resistance value is generally 0 kΩ to 2 kΩ; therefore, the maximum fault resistance is set as 2 kΩ in the paper. Then, select the faulty line according to the proposed method, and the parameters of VSBS are as follows: and CR = 1500. In addition, we use (, 0°, and 300 Ω) to indicate the fault occurrence in line number 1 when its initial angle is 0° and faulty resistance is 300 Ω. The results of their specific cross-correlation coefficients are shown in Table 2 in which represents cross-correlation coefficient of characteristic currents and .

The paper takes (, 90°, and 2000 Ω) as an example to explicitly show the consequence of VSBS. Under this fault condition, the TZSC of and of are shown in Figures 7(a) and 7(b); the cross-correlation coefficients are shown in Table 3.

The comparison of Figures 7(a) and 7(b) shows that, after the treatment of VSBS, the changing speed of in initial stage becomes slow, and oscillation part is reduced, which makes the changing trend of in initial stage easier to identify than that of TZSC in initial stage. It can be seen from Figure 7(b) that waveform of faulty line is steadier than that of nonfaulty lines, because frequency compression makes the part in low frequency band and with strong intensity easy to detect by VSBS, and the part in high frequency band and with weak intensity may be ignored, besides, the intensity of low frequency band fault transient component in the faulty line is much larger than that of nonfaulty line, and, therefore, the characteristic current waveform of faulty line is steadier than that of nonfaulty line. Besides, through simulation, we discover that because the transient free oscillation components and zero-sequence steady-state components are offset, the increase of fault resistance also makes waveform of each line steadier.

It can be seen from Table 3 that of and other lines are equal to −0.9732, −0.8092, and −0.7535, respectively; since all of them are the same sign, the output is “”; of and other lines are equal to −0.9732, 0.7408, and 0.6597, respectively; since all of them are not the same sign, the output is “1”; therefore, we judged as faulty line, and the judging result is consistent with simulation.

In summary, we know from Table 2 that TZSC of various fault conditions is, respectively, input into VSBS, and the output judging results are consistent with actual fault situations. Therefore, the proposed method can accurately achieve faulty line selection under different fault resistance and initial fault angle.

5.3. Simulation Results of Fault with Random Gauss White Noise Added

Since signals collected in actual system with fault often carry noise with them, to verify the antinoise performance of the method, we added 0.5 db or −0.5 db noise intensity to TZSC when fault in different line occurred and set signal before fault as 0. The selection results and specific cross-correlation coefficients are shown in Tables 4 and 5, respectively.

Signal-to-noise ratio equaling −0.5 db and fault situations as (, 60°, 1500 Ω) are taken as an illustration. The TZSC with noise of and of are shown in Figures 7(c) and 7(d), and the cross-correlation coefficients are shown in Table 6.

Firstly, from the faulty line selection method and Table 6, cross-correlation coefficients of and other lines are all negative, so the output is “” and of other lines are different, so the output is “1”; therefore, we judge line as faulty line, which is consistent with actual fault situation. Then, comparison of Figures 7(c) and 7(d) shows that, with the disturbance of strong noise, even if the TZSC of each line is submerged in strong noise, the proposed method is still able to effectively extract the changing trend of TZSC in initial stage and can accurately judge the faulty line. Finally, Tables 4 and 5 indicate that, with the disturbance of different noise intensity, the changing trends of characteristic currents in initial stage between faulty line and nonfaulty line still have a better discrimination after the treatment of VSBS, so we can say the method shows a good antinoise performance. The definition of signal-to-noise ratio [28] shows that the smaller the ratio is, the larger the noise intensity will be; therefore, antinoise performance of the method in this paper is much better than the one proposed in [8], which added 15 db and 0.5 db noise.

5.4. Adaptation Analysis of Faulty Line Selection Method
5.4.1. Different Faulty Lines

When fault occurs in , , and , respectively, we carry out faulty line selection method proposed in the paper to verify its adaptability, and the results and specific cross-correlation coefficients are shown in Table 7.

We know from [25] that, with the introduction of cable lines, although the attenuation process of fault transient current becomes shorter, the frequency spectrum principal component of transient component will move to low frequency band, which helps to detect VSBS. Therefore, different line fault conditions will not affect selection results of the method, and excellent results can also be obtained with different fault resistance.

5.4.2. Different Fault Distance

Since the distance of fault point is different in actual fault situations, we carry out simulation of line , with different distance from the bus line, and the fault distance is set as 4.5 km, 7.5 km, 10.5 km, and 13.5 km, respectively. Select the faulty line with the method and the results are shown in Table 8, with specific cross-correlation coefficients shown in Table 8.

It can be seen that the selection results are consistent with actual fault situation, which indicates that the method can also achieve faulty line selection of fault with different distance situation, especially with high ground resistance in the end of line.

5.4.3. Influence of Different Initial Stage Length on Selection Accuracy

From the moment of fault occurrence, choose time length of different initial stage as 0~0.005 s, 0~0.01 s, and 0~0.015 s, initial fault angle as 0°, 30°, 60°, and 90°, respectively, fault resistance as 10 Ω, 50 Ω, 100 Ω, 500 Ω, 1000 Ω, 1500 Ω, and 2000 Ω, and faulty line as , , , and , respectively, that is, a total of different fault conditions. Then, use the method proposed in the paper to carry out faulty line selection, and the selection results as shown in Table 9.

Table 9 shows that time length of the initial stage can affect the selection accuracy. The longer the length of initial stage is, the lower the selection accuracy will be. The main reasons are as follows:(1)VSBS can well detect the changing trend of TZSC in initial stage, besides, TZSC is an oscillation attenuation signal whose initial value is 0, therefore, with smaller time length, and the cross-correlation coefficient has better representation.(2)When single phase-to-ground fault occurs in distribution network, TZSC of each line will increase suddenly, and the TZSC mutation direction between faulty line and nonfaulty line is opposite. However, in the following T/4 time period, this situation will not happen, so the increase of signal length will affect the overall changing trend of the signal and and, then, cause wrong judgment.(3)It is known from Section 3 that, with the increase of time length, would also decrease, which means that the characteristic current can not well extract the changing trend of TZSC, and it would lead to wrong judgment.

5.5. Adaptation Analysis of Faulty Line Selection Method

In order to compare with other faulty line selection methods, choose TZSC with (, 90°, and 10 Ω) fault situation as an example, and demonstrate it from the following two cases, respectively: with noise and without noise. With the disturbance of noise, signal-to-noise ratio of the added noise is −0.5 db, and antinoise performances of existing methods are emphatically analyzed. At the end, in different faulty conditions, the selection results, which are from different faulty line selection methods, are given.

5.5.1. Without Disturbance of Noise

VSBS Method. According to the method in the paper, input TZSC of each line to VSBS, use fourth-order Runge-Kutta method for numerical simulation, and calculate cross-correlation coefficients of every line , the results of which are shown in Table 10. Choose characteristic current of and display its waveform in 0.019 s~0.021 s, as is shown in Figure 8(a).

Wavelet Packet Method. We use db10 wavelet packet to decompose TZSC of each line by four layers. Choose characteristic frequency band according to the maximum energy selection principle [32], restructure it with single branch, and calculate cross-correlation coefficient of every frequency band, whose results are shown in Table 10. Choose characteristic frequency band of and display its waveform in 0~50, as is shown in Figure 8(b).

Wavelet Method. We use db10 wavelet to decompose TZSC of each line by four layers. Choose the approximation coefficients of the four-layer wavelet of each line as characteristic signal, restructure it with single branch, and calculate cross-correlation coefficient after the restructuring of every characteristic signal, the results of which are shown in Table 10. Choose approximation coefficient waveform of line and display its waveform in 0.019 s~0.021 s which is shown in Figure 8(c).

EMD Method. We use EMD algorithm [33] to decompose TZSC of each line. Choose the first intrinsic mode components (IMF1 component) after treatment as characteristic mode component, calculate cross-correlation coefficient of each IMF1 component, and the results are shown in Table 10. Choose IMF1 component of line and display its waveform in 0.019 s~0.021 s, which is shown in Figure 8(d).

Firstly, for waveform, the changing trend of initial stage waveform in Figure 8(a) is clearer than that in Figure 8(b), indicating that VSBS can better describe the changing trend of TZSC compared to wavelet packet transform, because when the initial value is 0, the brown particles are in potential peak position of bistable system, and any small disturbance will make the brown particles of bistable system move drastically, so the bistable system can well track the signal changing trend. Oscillating components in Figure 8(b) are more abundant than that in Figure 8(a), which means that the characteristic signal processed by wavelet packet could contain more frequency components; the reason is that wavelet packet has such good capability of time-frequency analysis that it can elaborately divide the high frequency and low frequency of signals, while the frequency compression and transformation of VSBS will make some frequency components lost.

The waveform of Figure 8(a) changed after fault occurred, and the changing amplitude is larger than that before fault, while the waveform of Figure 8(c) changed before fault occurred, and the changing amplitude after fault occurred is smaller than that before fault, showing that VSBS can better reflect the changing time and trend of TZSC compared to wavelet transform. In addition, the oscillation degree of Figure 8(c) is smaller than that of Figure 8(b), indicating that although wavelet transform has good time-frequency localization, its high frequency resolution is poor.

The oscillation degree of Figure 8(d) is the strongest, because IMF component obtained by EMD contains frequency component which changes with the signal itself and is more suitable for nonstationary signals like TZSC. However, similar to Figure 8(c), Figure 8(d) also changed before fault occurred, indicating that EMD algorithm has a weaker ability to describe changing time and trend of TZSC compared to VSBS.

Then, from the cross-correlation coefficient and faulty line selection results we can see that, with the method proposed in this paper, after processing with VSBS and EMD, only the cross-correlation coefficients of TZSC between line and other lines are all the same, so line is judged as faulty line, which is consistent with actual situations. However, processed by wavelet packet, the cross-correlation coefficients between characteristic signal of line and other lines are equal to −0.4887, −0.0218, and −0.2531, respectively, all of which are the same negative sign, and, in the same way, the cross-correlation coefficients between and other lines are equal to −0.2531, −0.5990, and −0.9548, respectively, which are also the same sign, so and are judged as faulty line, but this result is not consistent with actual fault situation. And, then, processed by wavelet algorithm, none of the cross-correlation coefficients of characteristic signal between one line and other lines are the same sign, so all the lines are judged as healthy line, which, obviously, is not consistent with actual situation. This shows that wavelet transform and wavelet packet transform are not suitable for faulty line selection in this paper.

5.5.2. With Disturbance of Noise

With the same method and steps of Section 5.5.1, taking the waveform of line in 0.019 s~0.021 s as an example, we add a strong noise with signal-to-noise ratio as −0.5 db for simulation, the results of which are shown in Figure 9 and Table 11. Figure 9(a) is obtained by the process of VSBS, Figure 9(b) is obtained by the process of wavelet packet algorithm, Figure 9(c) is obtained by the process of wavelet transform algorithm, and Figure 9(d) is obtained by the process of EMD algorithm.

As to waveform, there are no obvious differences between other figures and Figure 8 except that Figure 9(d) is submerged in noise. As to cross-correlation coefficients, we will choose cross-correlation coefficient between line and line for analysis: without noise, processed in turn by VSBS, wavelet packet, and wavelet algorithm, the value is 0.9065, −0.4887, and 0.4483, respectively, while, with noise, processed in turn by VSBS, wavelet packet, and wavelet algorithm, the value is 0.9193, −0.4851, and 0.4328, respectively. Thus it can be seen that, with noise, the cross-correlation coefficient values by VSBS, wavelet packet, and wavelet algorithm are of little difference, so all of them have better antinoise ability. However, the cross-correlation coefficient processed by EMD algorithm without noise is 0.1968, and, with noise, the value is −0.4851, which changes from positive correlation to negative correlation, and the change is large, so combined with Figure 8(d) we can say that the antinoise ability of EMD algorithm is weak.

In summary, VSBS can extract the changing trend in initial stage of weak TZSC with the disturbance of strong noise, and its performance is better compared to wavelet transform, wavelet packet transform, and EMD algorithm; therefore, we choose VSBS to extract characteristic frequency band of TZSC in this paper.

5.5.3. Faulty Line Selection Results from Different Method

In strong noise background whose signal-to-noise ratio is −0.5 db, when different fault occurs including different lines, faulty resistance, and initial phase, the VSBS, wavelet packet, wavelet, and EMD are employed to select faulty line, respectively, and their faulty line selection results are shown in Tables 1215, respectively.

Table 12 shows that VSBS has no misjudgment in strong noise background and different faulty conditions; that is, the VSBS method can select faulty line correctly. However, there are many misjudgments in Tables 1315. These data indicate further that the antinoise performance of VSBS is better compared to wavelet transform, wavelet packet transform, and EMD algorithm.

6. Conclusions

This paper proposes a novel faulty line selection method for distribution network based on VSBS theory, and our research gets the following conclusions:(1)VSBS has better recognition for TZSC, which can effectively extract the changing trend of TZSC in initial stage under different fault situations, and the method can accurately judge the faulty line. In addition, VSBS has better antinoise ability, which helps extract the changing trend of weak TZSC with the disturbance of strong noise, and its antinoise performance is better than that of EMD algorithm and harmonic selection criterion.(2)The changing trend of TZSC in initial stage (0~0.005 s) is used to judge faulty line, which can reduce calculation time and the requirements for hardware. Besides, for the characterization capability of changing time and trend of TZSC in initial stage, the method in this paper is better than wavelet algorithm and wavelet packet algorithm.(3)The inadequacies of this paper are as follows: the frequency compression ratio is obtained through experiment, which might cause deviation. In addition, high resistance to ground fault with −10 db strong noise needs further study owing to the insufficient sensitivity of the present research.

Appendix

Build the simulation model according to the parameters, make fault of line occur at the point 5 km from the bus, and change the initial fault angle (0°, 30°, 60°, and 90°) as well as ground resistance for simulation. Then, with the proposed selection method, the cross-correlation coefficients of each line and faulty line selection results are shown in Table 2.

Add 0.5 db or −0.5 db noise intensity to TZSC when fault in different lines occurs. And set signal before fault to 0. The selection results and specific cross-correlation coefficients are shown in Tables 4 and 5.

In Figure 5, is cable-overhead line, and is pure cable line; we carry out faulty line selection with the method proposed in the paper to verify its adaptability, the results of which are shown in Table 7, and specific cross-correlation coefficients are shown in Table 7.

Since the distance of fault point is different in actual fault situations, we carry out simulation of line , with different distance from the bus line, and the fault distance is 4.5 km, 7.5 km, 10.5 km, and 13.5 km, respectively. Select the faulty line with the method and the results are shown in Table 8.

Notations

VSBS:Variable scale bistable system
TZSC:Transient zero-sequence current.

Competing Interests

The authors declare no conflict of interests.

Authors’ Contributions

Xiaowei Wang and Jie Gao conceived and designed the experiments; Jie Gao performed the experiments; Qiming Cheng analyzed the data; Guobing Song, Xiangxiang Wei, and Yanfang Wei contributed reagents/materials/analysis tools; Jie Gao wrote the paper.

Acknowledgments

This work was supported by National Natural Science Fund (61403127) of China, Science and Technology Research (12B470003, 14A470004, and 14A470001) of Henan Province, and Control Engineering Lab Project (KG2011-15, KG2014-04) of Henan Province, China, and Doctoral Fund (B2014-023) of Henan Polytechnic University, China.