Faulty Line Selection Based on Modified CEEMDAN Optimal Denoising Smooth Model and Duffing Oscillator for Un-Effectively Grounded System

As the un-effectively grounded system fails, the zero-sequence current contains strong noise and nonstationary features. -is paper proposes a novel faulty line selection method based on modified complete ensemble empirical mode decomposition with adaptive noise (MCEEMDAN) and Duffing oscillator. Here, based on multiscale permutation entropy, fuzzy c-means clustering, and general regression neural network for abnormal signal detection, theMCEEMDAN is proposed.-e endpoint mirror method is used to suppress the endpoint effect problem in the decomposition stage. -e proposed algorithm is able to decompose the original signal into a series of intrinsic mode functions, which can complete the first filtering. -e research shows that it can efficiently suppress the mode confusing phenomenon of empirical mode decomposition (EMD) and is also more complete and orthogonal than ensemble empirical mode decomposition (EEMD) and complementary ensemble empirical mode decomposition (CEEMD). -e optimal denoising smooth model is established for choosing optimal intrinsic mode functions to complete the second filtering. It can ensure that the reconstructed filtered signal has better smoothness and similarity. -e optimal denoising smooth model of MCEEMDAN can not only keep useful details of the original signal but also reduce the noise and smooth signal. -e bifurcation characteristic of the chaotic oscillator is applied in weak signal detection. -e zero-sequence current’s denoising result is extracted as the input signal of the Duffing system.-e faulty line could be selected by observing the phase diagram of the system. -e research results verify the usability and effectiveness of the proposed method.


Introduction
e un-effectively grounded system is used in the distribution network. Data show that 80% of the fault occurring in the un-effectively grounded system is single-phase grounding fault. When a single-phase grounding fault occurs, the fault current is pretty small, and the zero-sequence current is nonstationary and nonlinear [1,2]. If the faulty line remains unsolved, it could cause a serious threat to the insulation of the distribution network. It is important for stable operation and safe of the distribution network to select the faulty line accurately [3,4]. But the adaptability and accuracy of traditional faulty line selection methods are not enough due to the weak fault signal and the complex working condition.
Most existing faulty line selection methods use the injection signal method, steady-state information method, or transient information method. e injection signal method needs additional signal injection equipment, and its price is high, and its engineering operation is complex. e characteristic of a fault signal is weak in the steady-state information method, and it is greatly affected by the arc suppression coil, which makes the faulty line selection unreliable [5]. e characteristic of the fault signal in the transient information method is larger than the steady characteristic, and it is not affected by the arc suppression coil [6] and it does not require additional equipment. Liu et al. [7] used wavelet analysis to extract fault information of the zero-sequence current for faulty line selection, but the wavelet analysis is easily affected by noise, and the fault feature band which is extracted by wavelet analysis is inefficient, and the method needs to select appropriate wavelet basis function. Shu et al. [8] used time-frequency characteristics to the analysis of transient zero-sequence current for the faulty line selection, but it is invalid while a high-resistance grounding fault. Fault characteristics as an initial abrupt wave, real power, and transient are employed in the proposed fusion method [9], but its calculation amount is large. e paper [10] shows that the 5 th harmonic fault signal at the exit of the faulty line is larger than any nonfaulty line 5 th harmonic signal and towards the opposite, but the 5 th harmonic fault signal is very weak, and the fault background signal interferes with the 5 th harmonic fault signal. Zhang et al. [11] used EMD to decompose the zero-sequence current, and the first intrinsic mode function is extracted as the characteristic parameter of the zero-sequence current. But when the signal contains pulse signal, intermittent signal, and noisy signal, the mode mixing of EMD would arise and cause wrong judgment [12,13]. Ensemble empirical mode decomposition (EEMD) can restrain modal aliasing of EMD at a certain level; however, the calculation amount grows, and the completeness loses due to residual noise [14,15]. Complementary ensemble empirical mode decomposition (CEEMD) can decrease the reconstruction error by adding white noise to the target signal and using the ensemble mean to extract intrinsic mode functions [16,17]. However, its calculation speed is slow, and the efficiency is low. Some effective 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) adds adaptive white noise to the original signal and calculates unique residual signal to obtain every intrinsic mode function (IMF) [18,19]. e CEEMDAN algorithm can improve computational efficiency and decrease the reconstruction error. But CEEMDAN still has noise existing in the IMF region, and some IMFs may appear in hysteretic IMFs.
In order to overcome the above deficiencies, modified complete ensemble empirical mode decomposition with adaptive noise (MCEEMDAN) is proposed. e modified method based on combing CEEMDAN, general regression neural network, fuzzy c-means clustering, and multiscale permutation entropy is established. e decomposition stage uses an endpoint mirror method to suppress end effects. MCEEMDAN, which restrains the mode mixing problem, has excellent adaptability. e approach can also sift out the pseudocomponents generated by the signal decomposition. It is complete and orthogonal than EEMD and CEEMD. e optimal denoising smooth mathematical model, which weights the contradiction between similarity and smoothness of the filtered signal, is established for completing signal reconstruction and effectively extracting useful information from the original signal. e zero-sequence current in each line, which is processed by the MCEEMDAN optimal denoising smooth mathematical model, is extracted as the periodic external dynamic of the Duffing system. e trisection symmetry phase estimation method is applied for searching the critical phase, and each input signal is moved according to the critical phase. e faulty line could be selected by observing the phase diagram of the Duffing system.

Multiscale Permutation Entropy
e permutation entropy (PE) can describe the random and dynamic behavior of the time series [20]. It is based on a comparison of neighborhood values that PE measures signal complexity, and it does not consider the specific value of the time series. PE includes many advantages, such as simple algorithms, high computational efficiency, good robustness, strong anti-interference ability, and so forth. e PE algorithm is described below. Given a time series {x(i) i � 1, 2, . . ., N} with the length N, the time series which is constructed by using the following equation is defined as [21] where τ represents the time delay and m represents the embedding dimension. us, X(i) is arranged in ascending sorting, and the permutation is defined as For Any vector X(i) has a unique symbol sequence as follows: where g � 1∼k and k ≤ m!. e ascending sequence is one of the m! permutations of m-tuple vector. e occurrence frequency of each permutation is defined as P i (i � 1, 2, . . ., k). e PE is defined by the Shannon entropy, which is calculated as follows: where 0 ≤ H p ≤ ln(m!), and the normalized PE can be defined as where 0 ≤ PE ≤ 1. PE is an effective method for reacting a slight change in the time series. In order to solve the singlescale problem of PE, multiscale permutation entropy (MPE) is proposed by the paper [22]. e most important step of calculating MPE is acquiring coarse-grained time series by using the following equation: 2 Mathematical Problems in Engineering where s is the scale factor and y s (j) is a coarse-grained time series. PE of each coarse-grained time series is calculated as MPE.
Before computing MPE, three important parameters, including the embedding dimension m, the time delay τ, and the scale factor s are needed to set. If the embedding dimension m is too small, it is difficult to distinguish the randomness of the signal. e method will be meaningless. If the m is too large, the algorithm could take a long time to calculate, and it is difficult to observe the subtle change of the signal. Zheng et al. [23] suggested that the embedding dimension m is set to 6 in this paper. e time delay τ is set to 1 in this paper because it has little effect on the permutation entropy. If the maximal scale factor s max is >10, the MPE can reflect the important information of the signal [24]. e scale factor s is set to 12 in this paper.

Fuzzy C-Means Clustering.
Many areas, such as data mining and pattern recognition, use fuzzy c-means clustering (FCM) to solve their problems. FCM, which is an unsupervised learning, uses a membership degree to determine the value of each sample objective, which belongs to a cluster. FCM continuously modifies the cluster centers and optimizes the cost function by calculating the Euclidean distance between the data points and the cluster center until the termination criterion is satisfied. In FCM, membership values range between 1 and 0. e label of each sample objective is based on the maximum coefficient of the membership. Similar sample objects are assigned to the same cluster.
X � {x 1 , x 2 , . . ., x n } is a data set that is composed of n samples of dimension. It is divided into C � {c 1 , c 2 , . . ., c c } fuzzy clusters. e cost function is defined as [25] where m ∈[1, ∞) is the fuzzy weighted exponent, u ij is the membership matrix of the jth sample point to the ith cluster, and d ij is the Euclidean distance between x j and c i . We use Lagrange multipliers λ � [λ 1 , λ 2 , . . ., λ n ] with the length n to form a new cost function for equation (7) to reach a minimum as follows: where J is a new cost function, and the necessary conditions for equation (7) to reach its minimum are We repeatedly update equation (9). If either J's value is below a certain threshold value or J's improvement value minus its previous iteration value is below a certain threshold, the constraint condition is satisfied, and the cyclic iteration is stopped. e cluster to which the individual belongs is determined by the membership degree of individual relative to each cluster center, and each cluster is labeled for representing the rule degree of the fuzzy cluster signal.

General Regression Neural Network.
e generalized regression neural network (GRNN) is a critical branch of radial basis function neural network. GRNN includes many advantages, such as the strong nonlinear mapping ability, good global convergence, high fault tolerance, fast training speed, and so forth. e schematic diagram of the GRNN is shown in Figure 1.
GRNN consisted of three layers, including the input layer, radial basis layer, and specific linear layer [26]. e input layer is composed of neurons and does not participate in the operation. It only transfers the input vectors to the radial basis layer, and the number of layers is equal to the dimension of the input vector. e number of radial basis layers is equal to the number of the learning samples, and the neuron transfer function is Gaussian function p i as follows: where n is the number of radial basic neurons, X is the input variable, X i is the sample observation corresponding to the ith neuron, and δ is the smoothness factor. e transfer function of the output layer is a linear function purelin. e final network output value is defined as To sum up, the structure of GRNN is simple, and the smoothness factor δ will affect the execution effect of the GRNN. If the smoothness factor δ is selected randomly by hand, it could have greater subjectivity and affect the output Mathematical Problems in Engineering 3 result of GRNN. In this paper, the fruit fly optimization algorithm (FOA) [27] is used to select the optimal smoothness factor δ adaptively.

FCM-GRNN Identifies Abnormal Signal Model.
We collect 1600 signals. ese signals include 492 regular signals and 1108 random signals. ere are white noise signals, Gaussian white noise signals, intermittent signals, and other signals in random signals.
e MPEs of the signals are clustered by the FCM algorithm. Because these signals have two states, the clustering center is set to 2. e weighted index is set to 2, and the iterative termination tolerance is set to 10e − 3 . e input is MPE of each signal. e clustering center of each cluster and fuzzy membership matrix are obtained. We calculate the within-class means of two categories and obtain the distance matrix from all samples in the corresponding category to within-class mean. e 100 samples with the smallest distance are, respectively, selected from the distance matrix of the two categories, and the corresponding network outputs are set to 1 and 2. e MPE of 12 × 200 groups of signals is taken as the input of GRNN, and the output is the corresponding signal category. e FOA algorithm optimizes the smoothness factor of GRNN. e iteration number is set to 100, and the population size is set to 50. e optimization process diagram of FOA-GRNN is shown in Figure 2. As shown in Figure 2, when the smoothness factor is 0.1212, the minimum RMSE is 0.0145. e smoothness factor 0.1212 is brought into the GRNN, and the above MPE of 12 × 200 sample signals is trained to obtain the prediction model. e prediction model is used to analyze the MPE. e MPE as the eigenvector cannot be visualized in a high dimension. In order to show the effect of the model, the principal component analysis (PCA) algorithm is used to reduce the dimensionality of the 12-dimensional MPE and the clustering center vector. Figure 3 shows the 2-dimensional spatial distribution and 2-dimensional countermap of the random signals and regular signals after PCA dimension reduction. e MPE Input layer Radial basis layer Specific linear output layer  eigenvectors of signals are distributed around the two clustering centers according to regularity and randomness. e distinguishing effect of different types of signals is ideal, and the characteristics of the same type of signals are closely related. is shows that the signal's MPE as the characteristic quantity combining with FCM unsupervised clustering and GRNN tutor learning can complete the random detection of the signal, and the detection rate is maintained at 100%. is method can be used to identify abnormal signals.

Modified CEEMDAN Optimal Denoising
4.1. Modified CEEMDAN Algorithm. EEMD and CEMMD make the distribution of the extreme points of the original signal more uniform by adding white noise, and the added noise covers the abnormal signals, such as high frequency, intermittent, or noise in the original signal. erefore, the problem of modal confusion can be reduced. However, the algorithms need to limit the number of iterations, and the decomposed components do not necessarily satisfy the definition of the IMF. From the perspective of the application, it is meaningless to sacrifice the accuracy of the component and the physical meaning of the instantaneous frequency of the component for the adaptive decomposition.
In fact, after the abnormal signal is decomposed, it is not necessary to use EMD to decompose the noises completely, and as long as the decomposition completeness of the added noise signal can be ensured. In order to solve these problems, MCEEMDAN is proposed in this paper. e MCEEMDAN can avoid unnecessary ensemble mean in EEMD and CEEMD. It can make the decomposition results satisfy the definition of inherent mode function. e proposed method can reduce the reconstruction error caused by adding white noise and ensure the completeness of decomposition. e key step of MCEEMDAN is the detection of abnormal signals. Given the advantages of the MPE-FCM-GRNN, MCEEMDAN can be designed. Here, based on MPE-FCM-GRNN for abnormal signal detection, the adaptability of MCEEMDAN also can be improved.
Supposing the operator E k (·) is the kth intrinsic mode function, which is decomposed by EMD. IMF q ′ is the qth intrinsic mode function, which is decomposed by CEEM-DAN. I(t) is the original signal. n i (t) stands for the ith white noise with the mean value of zero and the variance of one (i � 1, 2, . . ., NR, and NR is the times of performing EMD); β i− 1 stands for adaptive coefficient (i � 1, 2, . . ., NR). EMD adopts the end point mirror image method during the decomposition stage. e steps of MCEEMDAN are described as follows: (1) e first mode of CEEMDAN is obtained by EMD NR times to decompose the signal I(t) + β 0 n i (t) as follows: (2) e first residue is as follows: (3) e second mode of CEEMDAN is obtained by EMD NR times to decompose the signal r 1 (t) + E 1 (n i (t)) as follows: (4) When k � 2, 3, . . ., N, the kth residue and the (k + 1) th mode are as follows: Step (4) is repeated until either the mean value of the upper and lower extreme envelope of the residue is zero or the difference between the number of the zero points and the number of the extreme points is less than two. (6) e MPE of each intrinsic mode function of CEEMDAN decomposition is analyzed by using trained GRNN. When the output classification result is the same as the random signal, it is considered to be an abnormal signal and is removed from the original signal. Supposing the pth mode is an abnormal signal, then the residue signal R(t) is defined as follows: (7) e R(t) is decomposed by EMD to obtain the MCEEMDAN's intrinsic mode functions, which are arranged in the order from high frequency to low frequency.
e signal x is composed of sinusoidal signal x 1 � 8sin (2π40t + π/2), amplitude modulation signal x 2 � (t + 1) sin (2π10t + π/3), and intermittent random signal x 3 . e sampling frequency is 1000 Hz, and the total time is 2 s. en, the waveforms of simulation signal and its components are shown in Figure 4. Because the results of CEEMD and EEMD are the same, only the result of EEMD decomposition is given. e x is, respectively, decomposed by EMD, EEMD, CEEMDAN, and MCEEMDAN, and the results are shown in Figure 5. It can be seen from waveform in Figure 5 that there is an obvious problem of modal aliasing in the EMD algorithm. EEMD and CEEMD can restrain modal aliasing of EMD at a certain level, but there are still a lot of pseudocomponents. e decomposition effect of the MCEEMDAN algorithm is better. e noise signal with strong randomness is filtered out. Both the residual      e end effect problem can be suppressed. e decomposition results are regular and stable. Figure 6 shows the Hilbert-Huang spectrum of four decomposition methods. As shown in Figure 6, EMD, EEMD, and CEEMDAN have the problem of instantaneous frequency instability, but the MCEEMDAN decomposition method can effectively extract the frequency information of the signal, and the frequency is more stable. e quantitative analysis of the four decomposition methods is based on the correlation between the components of the original signal and the intrinsic mode functions of four decomposition methods as follows: where j is the numbering of intrinsic modal function and i � 1, 2, 3 denotes signal x 1 , x 2 , x 3 , and C is a covariance matrix. As the correlation degree |corr| is close to 1, the correlation is high. e IMF correlation analysis of four decomposition methods can be known from Tables 1-4. EMD cannot decompose the effective components of the original signal. Although x 1 , x 2 can be extracted by EEMD, both IMF 3 -IMF 4 and x 2 have a high correlation, indicating that there is modal aliasing in EEMD. e CEEMDAN has a high correlation in both IMF 4 -IMF 5 and x 1 , which indicates that CEEMDAN has the problem of model mixing. MCEEMDAN can effectively extract x 1 and x 2 from x, and the only IMF 1 is highly correlated with x 1 , and the only IMF 2 is highly correlated with x 2 . As described above, MCEEMDAN can restrain modal aliasing. e orthogonality index can evaluate the endpoint effect of the algorithm. e smaller the orthogonal index is, the better the orthogonality of each intrinsic modal function is, and the lower the influence of the endpoint effect on the algorithm is. Table 5 shows orthogonality indexes for different methods. As shown in Table 5, the endpoint effect of the MCEEMDAN decomposition method is effectively suppressed.
e simulation results show that the proposed algorithm has better decomposition results than the existing algorithms. MCEEMDAN can restrain model mixing and overcome the problem of the endpoint effect. Its frequency of the Hilbert-Huang spectrum is stable. e orthogonality of the intrinsic mode function is good. ese advantages provide a good basis for signal reconstruction.

Optimal Denoising Smooth
Model. Different combinations of intrinsic mode functions can construct different reconstructed signals. e optimal denoising smooth model is constructed by the similarity degree and smoothness index [28][29][30]. e standard deviation of the difference between the reconstructed signal HP(t) and the original signal x(t) is calculated as the similarity degree as follows: where X(t) is the difference between x(t) and HP(t), N is the number of sampling points, and Amse represents the standard deviation. e smaller the Amse is, the closer the reconstructed signal is to the original signal. When point x 1 is on the curve f(x), the left and right curvatures at the point x 1 are as follows: When f(x) is smooth and differentiable at point x 1 , the left and right curvature at point x 1 are equal as follows: e second derivative is expanded according to the discrete formula as follows: where h is the step length. en, the smoothness SN| x�x 1 of the curve f(x) at point x 1 is defined as Supposing smoothness index Asmse is the standard deviation of SN of all points excluding the left and right endpoints in the sample. e smaller the Asmse is, the smoother the signal is.
Considering the smoothness of the filtering algorithm and the error standard deviation between the reconstruction result and the original signal, a weight coefficient μ(0 ≤ μ ≤ 1) is set, and then the objective function is defined as When Aminf is the minimum value, the objective function obtains the optimal solution. e algorithm model corresponding to the optimal solution is the optimal denoising smooth model.

Duffing Oscillator Model.
Duffing equation, which is widely used in weak signal detection, is a typical mathematical model to describe the chaotic property of the system [31,32]. e chaotic system is sensitive to the periodic driving force. When the total periodic driving force of the system changes, the system will show different trajectories.
e Duffing equation is defined as where k is the damping ratio, r and ω are the amplitude and the angular frequency of the inner periodic driving force, respectively, and − x 3 + x 5 represents the nonlinear restoring force. Equation (26) is expressed as the differential equation: When the damping ratio k and the angular frequency ω of the system are fixed, the trajectory of the system changes with the change of driving force. When r exceeds a certain threshold r c , the trajectory of the system enters the chaotic state. If r continues increasing over a certain threshold r d , the trajectory of the system enters the large-scale periodic motion state [33,34]. When the signal is detected by the Duffing oscillator, the amplitude of the inner periodic driving force is adjusted to r d . e system is in a critical state of chaos, and then the cofrequency signal to be tested is added to the Duffing system as follows: x ″ + kx ′ − x 3 + x 5 � r d cos(ωt) + a cos(ωt + φ), (26) where acos(ωt + φ) is a signal to be measured. Duffing's total driving force F(t) is composed of the internal periodic driving force and external periodic driving force as follows: � r cos(ωt + θ). (27) en, the following equations come into existence: Multiply the signal to be measured by a detection factor, and the inequation a ≪ r d can be obtained. e value of θ will tend to 0, and its impact on the system will be negligible. At this time, the state of the system is mainly determined by the value of φ. When the value of φ is in the following range as the r < r d and the system does not transition to the largescale periodic motion state. If the value of φ is not in this range, the r >r d and the system can transition to the largescale periodic motion state. After the fault of the distribution power network, the phases diagram trajectories of the fifth harmonic zero-sequence current of the normal line and faulty line are opposite, so the faulty line can be selected by the phase diagram trajectory of the chaotic oscillator.
e above-discussed case is that the internal power angular frequency is ω � 1 rad/s. But the frequency of the signal to be detected is a high-frequency signal in the power grid, and it is necessary to transform the scale of the Duffing system as follows: e detection equation can be defined as Equation (31) is expressed as the differential equation: x ′ � ωy, e above change process is an equal amount of mathematical transformation, but only the scale of the chaotic system is different, and the characteristic of the system equation is not changed. After the change of scale, the state and change rule of the chaotic oscillator are not changed. As long as the value of ω is adjusted, the detection of the high-frequency signal can be completed. Because the field operation is complex and subject to electromagnetic interference, the phase diagram trajectory of the chaotic system will become rough, and excessive noise will affect the detection of the chaotic system. e signal needs to be processed by the optimal denoising smooth model of MCEEMDAN before it is input to the Duffing system.

Trisection Symmetry Phase Estimation.
e trisection symmetry phase estimation is used to determine the phase of the signal to be measured relative to the internal driving force [35]. e phase is divided into three equal parts within a 2π circle. e MCEEMAND optimal denoising smooth model of any line is selected as the signal to be measured. When the signal is input into the chaotic oscillator by phase shift at 0, 2π/3, and 4π/3, there must be a phase diagram state which is different from the other two phase diagram states. For example, the Duffing system is a periodic state after the phase shift is 0°and 4π/3, but it is a large-scale periodic motion state after the phase shift is 2π/3, and then the critical Mathematical Problems in Engineering 11 phase must be within the (0, 2π/3) or (2π/3, 4π/3). e dichotomy is used to reduce the search area of the critical phase within the (0, 2π/3) and (2π/3, 4π/3) until the critical phase is searched. e results of MCEEMDAN optimal denoising of other lines that are moved according to the critical phase are input into the chaotic system to complete the signal detection. When the output state diagram of one of the lines is different from that of other lines, it is considered to be the faulty line. When all output state diagrams are the same, the bus bar is deemed to be the fault.

Case Study
MATLAB-SIMULINK is used in this paper to simulate the un-effectively grounded system. e simulation model consists of four lines. e length of a single line is generally 2 km to 20 km in the un-effectively grounded system. Line l 1 is set as 16 km overhead line; line l 2 is set as 13 km cableoverhead line including 5 km overhead line and 8 km cable line; line l 3 is set as 20 km overhead line; line l 4 is set as 10 km cable line. Table 6 shows unit length line parameters. e overcompensation method is adopted in the arc suppression coil. When the compensation degree is set to 8% in this paper, the inductive reactance of the arc suppression coil is 0.4765H. e resistance value of the arc suppression coil takes 2% of the inductive reactance and is calculated to be 2.9943 Ω.
Case 1. (low-resistance grounding). We use l 3 , 0.02 s, and 150 Ω to indicate the fault occurrence in line number 3 when its faulty resistance is 150 Ω, and the fault occurs at 0.02 s. Fault of line l 3 occurs at the point 6 km from the bus. In order to perfectly simulate the actual situation in the field, Gaussian white noise is artificially added to the zero-sequence current in each line, and the noise intensity is − 20 dBW. Zero-sequence current with noise in each line while a low-resistance grounding fault is shown in Figure 7.
e zero-sequence current in each line is analyzed by Fast Fourier Transform (FFT), and the result is shown in Figure 8. Figure 8 shows that due to the compensation effect of the arc suppression coil, the fundamental component of the zerosequence current of the fault line is not the largest component, and the zero-sequence current fundamental amplitude method for faulty line selection fails. In addition, the fifth harmonic is little affected by the arc suppression coil, but its content is less and the amplitude is small, which makes it difficult to select the faulty line. In a word, the application of FFT is limited.
e zero-sequence current of four cycles after fault is decomposed by MCEEMDAN proposed in this paper. Set the number of adding noise as 100 times, the amplitude of adding noise as 0.3, and the maximum number of iterations as 100 times. e results of MCEEMDAN decomposition in each line are shown in Figure 9. It can be known from Figure 9 that the decomposed intrinsic mode functions are stable, and the noise component is suppressed effectively.
According to the optimal denoising smooth model, the optimal filtering algorithm is selected. e results of optimal denoising smooth model in each line are shown in Figure 10.
It can be known from Figure 10 that the optimal denoising signal is smooth and consistent with the change of the original signal. In short, not only can the optimal denoising smooth model of MCEEMDAN filter out noise but it also retains the useful information of the original signal.
When the external driving force is not added to the Duffing system, the chaotic oscillator is adjusted to the critical state. Set the damping ratio of the chaotic oscillator as 0.3; set the amplitude of internal driving force as 0.464602; set the angular frequency of internal driving force as 500π rad/s, and the calculated step size is 5 × 10 − 5 . When the amplitude of the external driving force signal is much smaller than that of the internal driving force of the system, the θ value of formula (29) can be ignored. e external driving force needs to be multiplied by a detection factor before it is input into the Duffing system. Set the detection factor as 0.001 by the experiment many times. e fourthorder Runge-Kutta method is used to solve the equation of state of a chaotic system. It can be seen from the critical state diagram in Figure 11 that the critical state is a chaotic state.
It is first input into the Duffing system that MCEEM-DAN optimal denoising result of line l 4 is multiplied by the detection factor 0.001. e critical phase is 0.208π, which is calculated by the trisection symmetry phase estimation. e results of MCEEMDAN optimal denoising of other lines are moved according to the critical phase, and then they are, respectively, input into the Duffing system. e phase diagram trajectory while a low-resistance grounding fault is shown in Figure 12.
It can be known from Figure 12 that the phase diagram trajectory of line l 3 changes from the chaotic state to the large-scale periodic motion state, moving around the focal point (±1, 0) and the saddle point (0, 0). e phase diagram trajectories of other lines change slightly relative to the critical state, but they are still chaotic state. erefore, it can be judged that the faulty line is line l 3 , which is consistent with the experimental setting.
Case 2. (high-resistance grounding fault). We use l 2 , 0.02 s, and 1000 Ω to indicate the fault occurrence in line number 2 when its faulty resistance is 1000 Ω, and the fault occurs at 0.02 s. Fault of line l 2 occurs at the point 2 km from the bus. In order to perfectly simulate the actual situation in the field, Gaussian white noise is artificially added to the zero-sequence current in each line, and the noise intensity is − 1 dBW. Zero-sequence current with noise in each line while a high-resistance grounding fault is shown in Figure 13. It can be known from Figure 13 that the zero-sequence current signal is submerged in noise, and it is difficult to distinguish the changing trend of the signal.
In order to reflect the performance of MCEEMDAN optimal denoising, it is compared with the wavelet threshold noise reduction method, which is widely used at present. e wavelet function which uses the soft threshold method adopts db2 wavelet, and the number of decomposition layers is three. e waveforms of different denoising methods are shown in Figure 14. It can be known from Figure 14 that there is still noise in wavelet threshold denoising, but the    Mathematical Problems in Engineering effect of the MCEEMDAN optimal denoising in this paper is better, and the noise content is less than the result of wavelet threshold denoising, and the waveform surface is smooth. It is calculated that the correlation coefficient between MCEEMDAN optimal denoising and zero-sequence current without noise is 0.9932. ere are still some problems in the wavelet threshold denoising method, such as the selection of appropriate wavelet basis function and the number of decomposition layers. In a word, the method proposed in this paper has some advantages over the wavelet threshold noise reduction method. It is first input into the Duffing system that MCEEMDAN optimal denoising of line l 4 is multiplied by the detection factor 0.001. e critical phase is π/3, which is calculated by the trisection symmetry phase estimation. e final phase diagram trajectory is shown in Figure 15. It can be known from Figure 15     is line l 2 , which is consistent with the experimental setting. In summary, the proposed method is also suitable for highresistance grounding fault.  Figure 16. It is first input into the Duffing system that the MCEEMDAN optimal denoising result of line l 3 is multiplied by the detection factor 0.001. e critical phase is π/6, which is calculated by the trisection symmetry phase estimation. e results of MCEEMDAN optimal denoising of other lines are moved according to the critical phase, and then they are, respectively, input into the Duffing system. e phase diagram trajectory while an intermittent arc grounding fault is shown in Figure 17. It can be known from Figure 17 that the phase diagram trajectory of line l 3 is the chaotic state, but the phase diagram trajectories of other lines are of the large-scale periodic motion state. erefore, it can be judged that the faulty line is line l 3 , which is consistent with the experimental setting. In summary, the proposed method is also suitable for intermittent arc grounding fault.  Table 7. As can be seen in Table 7, the method proposed in this paper can accurately select fault lines under different fault conditions.

Conclusion
In this paper, MCEEMDAN, which is based on combining GRNN, FCM, and MPE, has strong adaptability. e MCEEMDAN is more complete and orthogonal than the commonly used EMD, EEMD, and CEEMDAN in the adaptive decomposition stage. e modified method can also suppress the endpoint effect problem and mode confusing phenomenon at a certain level. Different intrinsic mode functions are obtained by using the MCEEMDAN scheme, and a series of filtering algorithms are obtained, respectively. e superior filtering algorithm is selected by the optimal denoising smooth model. Not only can the optimal denoising smooth model of MCEEMDAN filter out noise, but it also retains the useful information of the original signal. In addition, its denoising ability is better than that of the wavelet threshold denoising method. Duffing system trisection symmetry phase estimation is used to detect MCEEMDAN optimal denoising results in each line. e faulty line can be selected by observing the trajectory of the phase diagram. Compared with the traditional faulty line selection method, the result of this paper is more intuitive. A large number of experimental studies show that the proposed method can accurately select the faulty line under different fault situations. e research results verify the usability and effectiveness of the proposed method. It has a good application prospect and engineering value in an uneffectively grounded system single-phase grounding protection that this paper proposes the faulty line selection method.

Abbreviations
MCEEMDAN: Modified complete ensemble empirical mode decomposition with adaptive noise GRNN: Generalized regression neural network FCM: Fuzzy c-means Multiscale permutation entropy EMD: Empirical mode decomposition EEMD: Ensemble empirical mode decomposition CEEMDAN: Complete ensemble empirical mode decomposition with adaptive noise CEEMD: Complementary ensemble empirical mode decomposition PE: Permutation entropy PCA: Principal component analysis IMF: Intrinsic mode function FOA: Fruit fly optimization algorithm.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.