Optimal Denoising and Feature Extraction Methods Using Modified CEEMD Combined with Duffing System and Their Applications in Fault Line Selection of Non-Solid-Earthed Network

: As the non-solid-earthed network fails, the zero-sequence current of each line is highly non-stationary, and the noise component is serious. This paper proposes a fault line selection method based on modified complementary ensemble empirical mode decomposition (MCEEMD) and the Duffing system. Here, based on generalized composite multiscale permutation entropy (GCMPE) and support vector machine (SVM) for signal randomness detection, the complementary ensemble empirical mode decomposition is modified. The MCEEMD algorithm has good adaptability, and it can restrain the modal aliasing of empirical mode decomposition (EMD) at a certain level. The Duffing system is highly sensitive when the frequency of the external force signal is the same as that of the internal force signal. For automatically identifying chaotic characteristics, by using the texture features of the phase diagram, the method can quickly obtain the numerical criterion of the chaotic nature. Firstly, the zero-sequence current is decomposed into a series of intrinsic mode functions (IMF) to complete the first noise-reduction. Then an optimized smooth denoising model is established to select optimal IMF for signal reconstruction, which can complete the second noise-reduction. Finally, the reconstructed signal is put into the Duffing system. The trisection symmetry phase estimation is used to determine the relative phase of the detection signal. The faulty line in the non-solid-earthed network is selected with the diagram outputted by the Duffing system.


Introduction
Nowadays, the non-solid-earthed network is widely used in the distribution network.When a single phase-to-ground occurs, the line voltage of the system remains symmetrical, and it can run for 2 to 3 hours after the fault occurs.However, the phase voltage of the non-fault line could rise, which would be a severe threat to the insulation of the system.It is necessary to remove faults quickly for safety in operation.There are many difficulties in fault line selection, such as weak fault current signal and complex actual operating environment [1][2][3].Due to the compensation effect of the arc suppression coil, the fundamental component of the zero-sequence current cannot select the fault line.However, the arc suppression coil has little influence on the high-frequency component, but the highfrequency component is non-stationary and nonlinear, and the traditional faulty line selection method would fail.It is thus necessary to conduct further research in the area.
At present, faulty line selection methods are roughly divided two ways: fault steady-state characteristic method and fault transient characteristic method [4,5].The steady-state characteristic method mainly uses the magnitude and direction of the fundamental zero-sequence current to form the criterion.However, due to the compensation effect of the arc suppression coil, the fundamental zero-sequence current direction of the faulty line and the non-fault line may be the same, and the fundamental zero-sequence current of the faulty line may not be the maximum.The result of the steady-state characteristic method may be unreliable.The fault transient characteristic method is more reliable because the transient characteristics are abundant, and it is not affected by the compensation of the arc suppression coil.Previous studies [6,7] used wavelet transform to decompose signals and extract fault features for faulty line selection.However, the wavelet basis function needs to be accurately selected, and the frequency band selected by wavelet transform under intense noise may not be valid.Zhang, Liu, and Piao [8] used a phase-locked loop to extract the fifth harmonic zero-sequence current, and the faulty line was selected by the fifth harmonic zero-sequence current characteristics.However, if the component of other harmonic signals is significant, it will affect the accuracy of faulty line selection.Zhang and co-workers [9] used empirical mode decomposition (EMD) to analyze transient characteristics.However, when there is intensive noise in the signal, the EMD algorithm would have the problem of mode confusion and cause the failure of faulty line selection.Kang and colleagues [10] used the matrix pencil method to analyze the characteristic frequency band of transient zero-sequence current, and the faulty line was selected according to the frequency component.However, this method did not take account of the effect of noise.Li and co-workers [11] used variational mode decomposition (VMD) to deal with nonlinear signals.VMD does not need to use the recursive decomposition process.It alleviates some deficiencies in EMD and has been applied in the field of signal processing.However, it is limited by the selection of the number of decomposition layers.Ensemble empirical mode decomposition (EEMD) improves the mode confusion problems of the EMD algorithm by adding different white noises to the original signal repeatedly [12][13][14][15].However, the added white noises can cause reconstruction error.Complementary ensemble empirical mode decomposition (CEEMD) can decrease the reconstruction error of the EEMD by adding opposite white noises to the target signal [16][17][18].However, if the amplitude of the added white noise and the number of iterations are not appropriate, the result would be many pseudorandom components.Some useful characteristics will be lost when the high-frequency signal is filtered, which will affect the denoising effect.Complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) has been applied for denoising of the nonlinear signals [19][20][21].However, CEEMDAN still has noise existing in the result of intrinsic mode functions, and some intrinsic mode functions could exist in hysteretic intrinsic mode functions.
This paper extracts the transient characteristics of the zero-sequence current accurately based on modified complementary ensemble empirical mode decomposition (MCEEMD) and the Duffing system.The pseudorandom components of CEEMD's intrinsic mode functions (IMF) are eliminated according to generalized composite multiscale permutation entropy (GCMPE) and support vector machine (SVM).The EMD decomposes the remaining intrinsic mode functions.The results are arranged in the order from high frequency to low frequency, namely, the MCEEMD algorithm.MCEEMD is a process of smoothing the signal, which can adaptively decompose the signal into a series of intrinsic mode functions according to the local scale characteristics of the signal.MCEEMD algorithm can suppress the mode confusion of the EMD algorithm and has better completeness and orthogonality than the CEEMD algorithm.The optimal intrinsic mode functions are selected to ensure that the filtered signal has preferable similarity and smoothness by an optimal smooth denoising model.The filtered signal can reduce the noise content and contain the useful components of the original signal.Wang and colleagues [22] proposed that the Duffing system is highly sensitive to a weak signal.In this paper, a new method of identifying chaotic characteristics based on texture features of the phase diagram is proposed.It has such advantages as visibility, a simple algorithm, and a small computational amount.Every line's filtered signal is input into the Duffing system as a characteristic value, and the faulty line is selected by the change of the phase diagram of the system.

Permutation Entropy
Permutation entropy (PE) can detect the randomness and dynamic behavior of a time series.It has a strong anti-interference ability and is suitable for dealing with nonlinear data [23,24].Given a time series {x(i) i = 1, 2, …, N} with length N, the time series is reconstructed as follows: where m is the embedded dimension, and τ is a time delay.The m elements in X(i) are rearranged in ascending order as ( ( 1) ) ( ( 1) , the permutation is ordered as where 0 1 p H ≤ ≤ .When the value of Hp is larger, the time series is rather more complex.On the contrary, the time series is rather more regular.

Multiscale Permutation Entropy
Since PE is only a single-scale analysis of time series, it has some limitations.In order to solve this deficiency, multiscale permutation entropy (MPE) is proposed.The most crucial step of this algorithm is to make the time series coarse-grained [25].Given a time series {x(i) i = 1, 2, …, N} with length N, the coarse granulation's process is calculated as where s is the scale factor, and ( ) s j y is a new time series at a scale factor, namely, coarse-grained series.The PE of each coarse-grained time series can be calculated.MPE is made of these PE values: Although MPE can overcome the single time scale problem of PE, the method still has some deficiencies.The process of coarse granulation is unique, and some vital information would inevitably be lost [26][27][28].The results of MPE may show an unstable trend.

Generalized composite multiscale permutation entropy
Zheng [29] proposed the GCMPE algorithm solving the deficiencies of MPE.The algorithm is listed as follows: Firstly, the PE values of all coarse-grained sequences in the same scale factor are calculated, and all entropy values are averaged.Secondly, the first moment (average) can be extended to the second moment (variance), namely the GCMPE method.The generalized coarse-grained sequence can be expressed as where =  .For scale factor s, we calculate s generalized coarse-grained sequences' PE values.
Then, GCMPE( , , , ) X s m τ can be expressed as The coarse granulation's processes of MPE and GCMPE with s = 3 are shown in Figure 1.The MPE can only decompose one sequence, so some crucial information would inevitably be left out.GCMPE can decompose three sequences.It can synthesize the info of multiple coarse-grained sequences at the same scale and extend the first moment (average) to the second moment (variance).Theoretically, GCMPE is better than MPE [29].
The coarse granulation's processes of multiscale permutation entropy and generalized composite multiscale permutation entropy with s = 3.

Parameter Selection and Comparative Analysis
To study the influence of the embedding dimension on the calculation results, MPE and GCMPE were used to analyze the Gaussian white noise under the conditions of m = 4, 5, 6, 7. The result is shown in Figure 2. In Figure 2, it can be seen that when the embedded dimension m increased, the fluctuation and deviation of MPE increased.However, as the embedded dimension m increased, the changing trend of GCMPE is relatively smooth, and the swing is small.The comparison results show the superiority of GCMPE.As the embedded dimension m is small, the value of GCMPE does not change obviously, which cannot reflect the advantages of multi-scale analysis.On the other hand, as the value of m is large, the reconstructed time series will be homogenized, and GCMPE cannot reflect the subtle changes of the sequence [29].Therefore, the embedded dimension m is set to 6 in this paper.If the maximal scale factor smax is >10, the GCMPE can reflect the crucial information of the signal.The time delay τ has little effect on GCMPE.Consequently, the maximal scale factor smax is set to 16, and the time delay τ is set to 1.

Signal Randomness Detection
The randomness tests of white noise, Gaussian white noise, intermittent signal, sine signal, amplitude modulation signal, and frequency modulated signal are analyzed by GCMPE.The result is shown in Figure 3. From Figure 3, it is easy to know that the GCMPE values of white noise, Gaussian white noise, and intermittent signals are larger, indicating that these signals are more random.In contrast, the GCMPE values of other signals are smaller, meaning that these signals are more regular.The obtained results are consistent with the actual situation.Namely, the GCMPE value of the random signal is greater than that of the regular signal.The SVM algorithm, which is suitable for small sample classification, is introduced to realize the intelligent detection of abnormal signals.We collected 239 signals, including 70 regular signals and 169 random signals.After obtaining their GCMPE values, the GCMPE values of 120 signals were randomly selected as the training samples of SVM, and the rest of the GCMPE values of the signals were used as the test samples.The penalty parameter c = 0.25 and kernel function parameter g = 0.0625 of SVM were calculated by K-fold cross-validation [30].The output result was 1, which represents the regular signal; the output result was 2, which represents the abnormal signal.The recognition rates of signal types are 100%.
Through the simulation experiment, it can be seen that GCMPE-SVM can be used to detect abnormal signals.

The Steps of MCEEMD
After the abnormal signal is decomposed, it is not necessary to use EMD to decompose the added noises completely.Given the advantages of the GCMPE-SVM, MCEEMD can be designed.Firstly, the original signal is decomposed by CEEMD.Secondly, the IMF decomposed by CEEMD is detected by GCMPE-SVM, and the abnormal signal is eliminated from the original signal.Finally, the residual signal is decomposed by EMD, namely MCEEMD algorithm.
The detailed steps of the modified method are as follows: Step 1: Add K pairs of opposite white noises to the original signal where ni(t) (i=1, 2, …, K) is white noise, and its mean value is 0; ai is noise amplitude, and its value is 0.2 times the standard deviation of the original signal I(t).
Step 2: EMD decomposes each signal in the set.Each signal can get a series of intrinsic mode function components.The j th intrinsic mode function of the i th signal is IMFij.
Step 3: The corresponding intrinsic mode functions are averaged to eliminate white noise as much as possible: where IMFj is the j th intrinsic mode function, which is obtained by CEEMD.
Step 4: The GCMPE of each intrinsic mode function of CEEMD decomposition is analyzed by using SVM.When the output value is 2, the intrinsic mode function is considered to be an abnormal signal and is removed from the original signal.Supposing p modes are abnormal components.Then the residue signal R(t) is defined as Step 5: The R(t) is decomposed by EMD to obtain the MCEEMD's IMF.The results are arranged in the order from high frequency to low frequency.

Analysis of the Simulation Signal Using EMD, CEEMD, and MCEEMD
In order to demonstrate the advantages of the MCEEMD, the simulation experiments are carried out in this section.The artificial signals are given by Equation ( 13    In Figure 5a, it is evident that the EMD algorithm has the problem of mode confusion.It is also seen in Figure 5b that the CEEMD algorithm produces false components.As can be seen in Figure 5c, the MCEEMD has no spurious components, and it removes the noise components.The correlation coefficient [1,16] can be used to analyze the degree of linear association between variables.The correlation coefficient is close to 1, indicating a higher correlation between the two variables.In addition, Tables 1-3 show the correlation coefficients of the three decomposition methods.As can be seen in Table 1, the EMD algorithm cannot extract x2 effectively.In Table 2, it is evident that CEEMD's IMF3-IMF4 and x2 have a higher correlation.Here, the CEEMD decomposition encounters the mode confusion at a certain level.In Table 3, MCEEMD can effectively extract x1 and x2, and it restrains mode confusion.The results show that the method proposed in this paper can effectively extract the useful components of the original signal, and it suppress mode mixing at a certain level.The index of orthogonality (Ort) can be used as a valid index to evaluate the inhibition of the end effect [31].The smaller the Ort, the better the decomposition effect.Ort is shown in Table 4 where the Ort of MCEEMD is -5

10 ×
, and the value is almost close to 0. The results show that the Ort of MCEEMD is better than that of the other two algorithms.

10 ×
Hilbert-Huang transformation gains its popularity to analyze signals in time and frequency domains [32]. Figure 6 provides the Hilbert-Huang spectrums for three decomposition methods.As can be seen in Figure 6a, the result of EMD has a visible phenomenon of frequency mixing at 10 Hz, and the instantaneous frequency has an apparent fluctuation.Figure 6b shows the decomposition result of CEEMD has many pseudo components between 50 Hz to 90 Hz, and its instantaneous frequency is also unstable.However, the result in Figure 6c shows that MCEEMD can extract the effective components of the original signal, and there are no pseudo components in the Hilbert-Huang spectrum, and its instantaneous frequency is stable.To summarize, it is evident that the proposed algorithm has better decomposition results and overcomes mode mixing.

Optimal Smooth Denoising
The optimal smooth denoising model is composed of the MCEEMD and the objective function.The objective function takes into account the similarity and smoothness of different filtering algorithms.MCEEMD decomposes the original signal I(t), and the signal reconstruction SR(t) can be obtained as where SRk(t) represents the k th signal reconstruction, and i is the number of the IMF.
The standard deviation of the difference value between the signal reconstruction SR(t) and the original signal I(t) is taken as the similarity standard [33], like where N is the number of sample points, and Amse represents the similarity standard.The smaller Amse, the closer the reconstructed signal is to the original signal.
If the function f(x) on the interval (a, b) is assumed to have a continuous first-order derivative, its figure is a curve with tangents everywhere, and the tangents rotate continuously with the movement of the tangent point.The curve is defined as a smooth curve.In particular, it is obtained that a necessary and sufficient condition for function derivable at a certain point is that the left and right derivatives of the point coexistent and are equal.Suppose the point x1 lies on the curve f(x).The left curvature and the right curvature at point x1 are as follows: where K − and K + are the left curvature and the right curvature, respectively.If f(x) is smooth and differentiable at x = x1, the left curvature and the right curvature at x1 are equal as follows: )" f x − and 1 ( )" f x + are expanded to the difference forms as follows: where h represents the sampling step.It is evident, therefore, that if f(x) is smooth and differentiable The smoothness of the signal reconstruction at x = x1 can be defined as where SN is the smoothness index of the signal reconstruction at x = x1.The standard deviation of SN values of all points is taken as the smoothness index of the whole signal reconstruction as Asmse std( ), SN = (20) where Asmse represents the smoothness standard.The smaller Asmse, the smoother the signal reconstruction SR(t).Thus, the objective function Aminf can be defined as Aminf Amse (1 )Asmse, where μ is the weight factor, and 0 1 μ ≤ ≤ .The factor μ can be valued by balancing the requirements of the similarity and smoothness.It is set to 0.6 in this paper.As Aminf gets the minimum, the objective function can achieve an optimal solution, and then match the optimal solution with the filtering algorithm, namely, the optimal smooth denoising model.

Optimal Smooth Denoising Experiment of Simulation Signal
In order to validate the effectiveness of the optimal smooth denoising model, it is utilized to analyze simulation signals.The simulation signals are given as follows: where x1 represents the clear signal, and x is the noisy signal.The x1 is shown in Figure 7a; x is composed of x1 and Gaussian white noise as shown in Figure 7b.The noisy signal x is decomposed by the MCEEMD.The results are shown in Figure 8. From Figure 8 it is easy to know that the MCEEMD decomposition results can get eight components.
The similarity standard Amse, the smoothness standard Asmse, and the objective function Aminf of each filtering algorithm are listed in Table 5.As shown in Table 5, based on these results, the filtering algorithm SR8 is thought to be optimal when the similarity standard Amse and the smoothness standard Asmse are comprehensively considered.In order to prove the excellent performance of the optimal smooth denoising model, the filtering algorithm SR8 is compared with the wavelet threshold denoising.The results are shown in Figure 9.As shown in Figure 9a, it can be seen that the waveform still contains noises.The signal-to-noise ratio of the wavelet threshold denoising is 10.9047 dB.The wavelet threshold denoising also confronts the problems of how to choose the wavelet base and wavelet decomposition level.For the optimal smooth denoising model, its signal-to-noise ratio is 11.0150 dB, and the similarity coefficient between its result and the clear signal x1 is 0.9971.To summarize, the proposed denoising method not only can retain the accurate information of the original signal but also can reduce the noise and smooth the original signal.

Duffing System Model
The Duffing system is widely used in weak signal detection [34][35][36].When it is in the critical state, a minimal change of its parameters can cause the change of the system state.For the fifth harmonic, the amplitude of the fault line is large relatively, and the phase of the fault line is opposite to that of the non-fault lines.The fifth harmonic components of the fault line and the non-fault line can change the Duffing system into different state trajectories.The normal form of the Duffing equation can be defined as 3 cos( ), x kx x x r t ′′ ′ + − + = (24) where k represents the damping ratio, and rcos(t) denotes the internal driving force.
The damping ratio k is fixable, and it is set to 0.5 in this paper.As the internal driving force amplitude r varies from small to big, the Duffing system can have different motion states.If r = 0, the system will have three singular points; (0, 0) is the saddle point and (±1, 0) is the center point.When r exceeds a certain threshold rc, the trajectory of the system turns into a chaotic state.If r continues increasing over a certain threshold rd, the trajectory of the system changes from the chaotic state to the large-scale periodic motion state.The state at rd is used to detect the signal.rd can be previously determined by the computer simulation experiment [37].Based on the experiments, we chose rd = 0.8253.For detecting the high-frequency signal, we must do time scale transformation.Defining where cos( ) a ωτ ϕ + is the input signal.The total driving force can be expressed as Since a<<rd, θ is very small, and its impact on the system can be negligible.If arccos( / 2 ) arccos( / 2 ) ≤ , the trajectory of the Duffing system is still in a chaotic state.If θ is not in this range, it makes the transition from the chaotic state to the large-scale periodic motion state.Making use of the characteristic of the Duffing system can select the fault line in the non-solid-earthed network.In this paper, the angular frequency of the internal driving force is set as 500π rad/s, and the calculated step size is set as 5 × 10 -6 .For ignoring the effect of θ , the input signal shall be multiplied by a detection factor.We chose different detection factors for experiments, and 0.01 had the best effect.

Trisection Symmetry Phase Estimation
Trisection symmetry phase estimation (TSPE) can determine the phase of the input signal relative to the internal driving force [38].The detection process of TSPE is shown in Figure 10.In Figure 10,   1  d φ and 2 d φ are the critical phase sites.Dividing the 2π area into three equal parts, we shift the phase of the input signal by 0, 2π/3, and 4π/3, respectively.The phase-shifted signals are input into the Duffing system in turn.There must be a phase's motion state, which is different from the other two phases' motion states.From Figure 10, it can be seen that 2 φ is in the chaotic state, however, 1 φ and 3 φ are in the large-scale periodic motion states.Hence, it is concluded that critical sites are located in A line's zero-sequence current, which is denoised by the optimal smooth denoising model of MCEEMD, is the selected signal.TSPE determines the critical phase.Moving all line signals synchronously according to the critical phase, we input them into the Duffing system.The fault line and the non-fault line can make the Duffing system in different state trajectories.

A Method of Distinguishing Chaotic Characteristics
In this paper, a method of automatic recognition of chaotic characteristics based on texture features is proposed.According to this method, we can easily and quickly judge that the system is in a chaotic state or a stable large-scale periodic state.The gray level co-occurrence matrix (GLCM) can describe texture features [39].We construct the GLCM of four directions [40], namely 0 , 45 , 90 and 135 .
We calculated the mean value of the same feature in four directions.The four texture features were defined as a texture feature vector.The large-scale periodic motion state is shown in Figure 11.We calculated the texture parameters of this phase diagram.The results are shown in Figure 12.The mean result is defined as the reference vector.We calculated the Euclidean distance between this reference vector and the texture parameter vectors of other phase diagrams.After many experiments, it was concluded that if the Euclidean distance is less than 1.5, the phase diagram is in the large-scale periodic motion state.Otherwise, the phase diagram is in a chaotic state.The method can avoid the involvement of human factors and automatically identify chaotic nature by computers.

Introduction of the fault line selection steps
When a single phase-to-ground fault occurs in the non-solid-earthed network, the system can generate asymmetric zero-sequence currents.The zero-sequence currents of all lines are processed as follows: Step 1: The zero-sequence current is decomposed by MCEEMD.We can get a series of intrinsic mode functions.
Step 2: The superior filtering algorithm can be selected by the optimal smooth denoising model, and it is used as the input signal of the Duffing system.
Step 3: The Duffing system is adjusted to the critical state.Every line's denoising result is input into the Duffing system.
Step 4: We calculate the texture parameters of all phase diagrams.The state of the system is judged by Euclidean distance.The phase diagram state of the fault line is different from that of the non-fault line.If all phase diagrams are in the same state, the bus bar is the faulty line.

Simulation Study
The MATLAB software was used in this paper to simulate the non-solid-earthed network.The simulation model is shown in Figure 13.The three-phase voltage source is set as 110 kV.The transformer ratio is 110 kV/10.5 kV, and its nominal capacity is 20 MVA.Table 6 also shows the impedance parameters and length of the lines.The equivalent inductance of the arc-suppression coil is 1.2102 H.  Suppose there is a single-phase-to-ground fault in phase A of line L4 at 0.02 s, 5 km away from the bus bar.The ground resistance Rg is 10 Ω.The zero-sequence currents collected in the actual non-solid-earthed network often carry noise with them.After the simulation, we add the Gaussian white noise to the zero-sequence currents in lines 1-4.The results are shown in Figure 14 from 0 to 0.1 s.Every line's zero-sequence current is done with the help of fast Fourier transform.The result of fast Fourier transform analysis is shown in Figure 15.As seen from Figure 15, the zero-sequence current contains minor amounts of the fifth harmonic.In a study by Zhang and colleagues [8], the method of fault line selection failed by comparing the amplitude of the fifth harmonic.MCEEMD is used to decompose the post-fault zero-sequence current.The decomposition results are shown in Figure 16.As seen from Figure 16, the results of each line are arranged in descending sequence by frequency.Different intrinsic mode functions can be obtained by using MCEEMD, and a series of filtering algorithms can be obtained, respectively.The optimum filtering algorithm can be determined by the optimal smooth denoising model.Every line's optimal smooth denoising result is shown in Figure 17.As shown in Figure 17, the proposed algorithm not only achieves noise reduction effect, but also reflects well the available information of the original signal.Every line's denoising result can be used as the input signal of the Duffing system.The Duffing system is adjusted to the critical state.Every line's denoising result is input into the Duffing system.The equation of the Duffing system is solved by the fourth-order Runge-Kutta.The Duffing system's phase diagram when adding each line's denoising result is shown in Figure 18.As can be seen in Figure 18, the phase diagram of line L4 is in a chaotic state while the others are in a large-scale periodic motion state.It can be seen that the trajectory of the large-scale periodic motion state moves around (±1, 0) and (0, 0).The fifth harmonic components of the fault line and the nonfault line can make the Duffing system in different state trajectories.Therefore, line L4 has a singlephase-to-ground fault.For automatically identifying chaotic nature, we extract texture features of phase diagrams.As introduced in Section 5.3, the reference vector of the texture feature is (0.9086, 0.2638, 3.3849, 0.1029).The chaotic nature can be identified according to the Euclidean distance between the reference vector and texture feature of the phase diagram.The results are shown in Table 7.As shown in Table 7, only the Euclidean distance of line L4 is more than 1.5, and therefore its phase diagram is in a chaotic state.However, the other phase diagrams are the large-scale periodic state.The results mean the method can select the fault line accurately.8.As can be seen in Table 8, the Duffing systems of the fault line and the non-fault line show different state trajectories.Table 8 also shows that the state of the system can be accurately identified by the texture feature recognition of the phase diagram.If Euclidean distance is less than 1.5, the phase diagram is the large-scale periodic motion state.Otherwise, the phase diagram is a chaotic state.To summarize, the proposed method is able to select the fault line with different situations.The optimal smooth denoising model can be determined by balancing the similarity and smoothness of different filtering algorithms.The superior algorithm can retain the useful features of the original signal and reduce the noise and smooth the original signal.(3) A novel method of identifying chaotic nature based on texture features of the phase diagram is presented.This way can avoid the involvement of human factors and automatically identify chaotic nature by computers.(4) The fault line in the non-solid-earthed network can be selected with the diagram outputted by the Duffing system.A large number of experimental studies show that the proposed method can accurately select the fault line under different fault situations.The research provides a novel train of thought for the fault line in the non-solid-earthed network.

Figure 3 .
Figure 3. GCMPE values for various types of signals.
), where x is composed by x1, x2, and x3, and the sampling frequency is 1000 Hz.The original signals are shown in Figure 4.

Figure 8 .
Figure 8.The MCEEMD results of simulation signal x.

Figure 9 .
Figure 9.The denoising results for different methods.(a) The denoising result of the wavelet threshold denoising.(b) The denoising result of the optimal smooth denoising model.

3 .
linear restoring force.The equivalent form of the Duffing equation can be expressed as

Figure 13 .
Figure 13.Simulation model of the non-solid-earthed network.

Figure 14 .
Figure 14.Every line's zero-sequence current with noise.(a) The result of line L1; (b) the result of line L2; (c) the result of line L3; (d) the result of line L4.

Figure 15 .
Figure 15.The result of fast Fourier transform analysis.

Case 2 :
Since fault conditions are different in the actual non-solid-earthed network, we carried out different fault situations, such as fault line, fault resistance, and fault distance.The selection results of fault in different conditions are shown in Table

Table 1 .
The correlation coefficients of EMD decomposition.IMF: intrinsic mode functions.

Table 2 .
The correlation coefficients of CEEMD decomposition.

Table 3 .
The correlation coefficients of MCEEMD decomposition.

Table 6 .
The impedance parameters and length of the lines.

Table 7 .
Texture feature and Euclidean distance of phase diagram.

Table 8 .
Selection results of fault in different situations.