A Novel Partial Discharge Ultra-High Frequency Signal De-Noising Method Based on a Single-Channel Blind Source Separation Algorithm

: To effectively de-noise the Gaussian white noise and periodic narrow-band interference in the background noise of partial discharge ultra-high frequency (PD UHF) signals in field tests, a novel de-noising method, based on a single-channel blind source separation algorithm, is proposed. Compared with traditional methods, the proposed method can effectively de-noise the noise interference, and the distortion of the de-noising PD signal is smaller. Firstly, the PD UHF signal is time-frequency analyzed by S-transform to obtain the number of source signals. Then, the single-channel detected PD signal is converted into multi-channel signals by singular value decomposition (SVD), and background noise is separated from multi-channel PD UHF signals by the joint approximate diagonalization of eigen-matrix method. At last, the source PD signal is estimated and recovered by the l 1 -norm minimization method. The proposed de-noising method was applied on the simulation test and field test detected signals, and the de-noising performance of the different methods was compared. The simulation and field test results demonstrate the effectiveness and correctness of the proposed method.


Introduction
Partial discharge (PD), which is caused by insulation defects, can be utilized to evaluate the insulation state of high voltage devices [1][2][3]. Since ultra-high frequency (UHF) signals (0.3-3 GHZ) will be also produced and propagated outside through the non-metal shielding dielectrics when PD occurs, the UHF antenna sensor can be used to detect the PD of high-voltage (HV) apparatus [4]. However, since the PD UHF signal amplitudes are very low in the early insulation defects, and vulnerable to interference due to the numerous external electromagnetic waves, the detected PD signal will have serious waveform distortion, even submerged in background noise interference. Thus, the de-noising method has attracted extensive attention from researchers. narrow-band interference signal and Gaussian white noise signal in the PD UHF signal is linear mixed, the mathematical model of linear mixed BSS can be expressed as [19,20].
x(t) = As(t) (1) where A is the mixed matrix with M in N order. Then, the unknown signal vector s(t) can be separated from the detected signal, x(t), by solving the separation matrix, W. Assume y(t) = [y 1 (t), y 2 (t), . . . , y N (t)] is the output signal vector after BSS. The mathematical model of BSS is y(t) = Wx(t) = WAs(t) (2) where W is the separation matrix with N in M order. When M ≥ N, it is non-underdetermined blind source separation. The unknown PD signal vector s(t) can be directly solved by Equation (2). When M < N, it is underdetermined blind source separation. Since the solution of Equation (2) is not unique, the optimal estimate solution of s(t) can be only calculated by the estimation method. However, for PD UHF de-noising, a single sensor is always used to detect the PD UHF signal (M = 1), and the number of source signals (N) (includes the PD signal, white Gaussian noise signal and other narrow-band interference signals) will be larger than three (N ≥ 3), which is called single-channel BSS. Hence, the BSS based on matrix computation is not feasible. The one-dimensional signal vector, x(t), detected by single channel should be decomposed into multiple time-domain signals. The single-channel PD signals can be converted into multiple-channel signals, through establishing the virtual channels, meeting the requirements of non-underdetermined blind source separation.

Number Estimation of Source Signals
To establish the virtual channels, the number of source signals (N) is needed. Firstly, the one-dimensional time domain signal is mapped to a two-dimensional time-frequency domain through the time-frequency analysis of detected source signals. The amplitude-time source signals can be converted into amplitude-time-frequency signal. Hence, the number of source signals (N) can be obtained.
Due to the advantages of flexibility of time-frequency resolution and convenience, S-transform is used to convert the source signals. S-transform function can be expressed as [21] s(τ, where f represents the frequency, t and τ represent the time, and g(t − τ, f ) represents the Gaussian window function.
Since the frequency distribution of periodic narrow-band interference is concentrated and the amplitude is large, a clustering analysis method of frequency section characteristic parameters, based on the time-frequency matrix of S-transform, is proposed to estimate the number of signal sources. The proposed clustering analysis method is as follows: Step 1: Solve the time-frequency matrix based on S-transform, the frequency sections are extracted at each 10 MHz interval, such as 0~10 MHz, 10 MHz~20 MHz, etc.
Step 2: Compare the amplitudes of all time-frequency sampling points for each frequency section. The maximum amplitude (AM max ) is characteristic parameter 1, and the minimum amplitude (AM min ) is characteristic parameter 2.
Step 3: Solve the average amplitude (AM av ) for all time-frequency sampling points for each frequency section as characteristic parameter 3.
Step 4: According to the characteristic parameters extracted in Step 2 and Step 3, the fast search and find of the density peaks clustering (FSFDPC) algorithm is proposed to cluster [22]. The cluster number, which is the amount of periodic narrow-band interference (N pn ), can be obtained. The FSFDPC algorithm has better convergence capacities.
Step 5: According to the source signal statistical characteristics of blind source separation, the number of PD UHF signals and the Gaussian white noise signal are equal to 1. Hence, the number of the source signals is N = N pn + 2.

Multi-Channel Detected Signal Recombination
For single-channel BSS, the traditional method of constructing virtual channels is the empirical mode decomposition (EMD), which can decompose the single-channel signal into intrinsic mode function signals. However, since the PD UHF signal and periodic narrow-band interference source signal are independent, the traditional method has the modal aliasing problem which will negatively affect the final signal separation result. Therefore, this paper proposed a novel multi-channel detected signal recombination method based on SVD. The trajectory matrix (X) of detection signal x(t) can be expressed as [23,24] where X is the trajectory matrix, U and V are the orthogonal matrices with m by n order, respectively, Λ a = diag (a 1 , a 2 , . . . , a n ) is a diagonal matrix, a 1 , a 2 , . . . , a n are the singular values of mixed matrix A.
Since matrix A with rank k can be expressed as the sum of k submatrices with rank 1, where u i and v i are the singular value vectors of the i-th column of matrices U and V, respectively. X i is a submatrix. Each submatrix in SVD is multiplied by the singular value weight and two eigenvectors are derived from matrices U and V, respectively. The submatrix and the singular value (a i ) have a positive correlation. Hence, since the time-frequency distribution of the PD signal and periodic narrow − band interference is relatively concentrated, the amount of information of the singular values is large. When the time-frequency distribution of Gaussian white noise signal is dispersed, the amount of information of the singular values is smaller. In addition, the submatrices for different singular values (a i ) are statistically independent, the corresponding submatrices with the first N singular values can be obtained through the SVD calculation of the single-channel detection signal, x(t). Then, the virtual channels are constructed by signal reconstruction calculation, and the multi-channel signals can be recombined for blind source separation de-noising.
Hence, the proposed multi-channel detected signal recombination method is as follows: Step 1: Construct a trajectory matrix (X) of signal x(t), as follows: Step 2: The SVD calculation for trajectory matrix X, and the singular value sequence a 1 , a 2 , . . . , a N is obtained.
Step 3: Extract the first N singular value sequence (a 1 , a 2 , . . . , a N ) and the corresponding submatrices. Reconstruct the virtual channels.
Step 4: Inversely calculate the trajectory matrix (C j ), according to Step 1. The virtual channel signal, C j (t), can be obtained.
Step 5: Recombine the virtual channel signal, C j (t), and the original single-channel detected signal, x(t), as the novel multi-channel detected signal, D(t) = [x(t), C 1 (t), C 2 (t), . . . , C N (t)]. It will become the non-underdetermined blind source separation.

BSS De-Noising Method Based on the JADE Algorithm
The separation matrix (W) is the most important for BSS. This paper proposed the JADE algorithm to solve the separation matrix, W. The JADE algorithm introduces a fourth-order cumulant of multivariate data to estimate separation matrix W by characteristics decomposition and joint approximation diagonalization [25]. The blind source separation de-noising method based on the JADE algorithm can be divided into two steps to solve the separation matrix (W) whitening process and orthogonal transform processes, as shown in Figure 1. Step 4: Inversely calculate the trajectory matrix (Cj), according to Step 1. The virtual channel signal, Cj(t), can be obtained.

BSS De-Noising Method Based on the JADE Algorithm
The separation matrix (W) is the most important for BSS. This paper proposed the JADE algorithm to solve the separation matrix, W. The JADE algorithm introduces a fourth-order cumulant of multivariate data to estimate separation matrix W by characteristics decomposition and joint approximation diagonalization [25]. The blind source separation de-noising method based on the JADE algorithm can be divided into two steps to solve the separation matrix (W) whitening process and orthogonal transform processes, as shown in Figure 1. Firstly, the multi-channel detected signal, D(t), should be whitening processed to eliminate the second-order correlation between the components. The whitening process can be expressed as [25] QAs Cs (8) where Q is the whitening matrix, D′(t) is the signal after the whitening process, in which each component is relatively independent and the variance is 1, and matrix C is the product of matrices Q and A. Assume Rx is the correlation matrix of the multi-channel detection signal D(t), and D(t) can be eigenvalue decomposed as  (9) where matrix

RE Λ E
is the orthogonal matrix, in which the diagonal elements are the eigenvalues of matrix Rx, and the corresponding standard orthonormal eigenvector is the column vector of characteristic matrix E. The whitening matrix (Q) can be expressed as Since the component variances of s(t) and D′(t) are both 1, each element of s(t) is independent, and each element of D(t) is orthogonal; matrix C is an orthogonal normalized matrix and can be expressed as TT N  CC C C I (11) where IN is the N-order unit matrix.
The matrix, D′(t), after the whitening process needs the orthogonal transformation to make aure that the components after blind source separation are not correlated and the variance is 1. The orthogonal transformation process can be expressed as

JADE Algorithm
Separation matrix W Firstly, the multi-channel detected signal, D(t), should be whitening processed to eliminate the second-order correlation between the components. The whitening process can be expressed as [25] D (t) = QD(t) = QAs(t) = Cs(t) (8) where Q is the whitening matrix, D (t) is the signal after the whitening process, in which each component is relatively independent and the variance is 1, and matrix C is the product of matrices Q and A.
Assume R x is the correlation matrix of the multi-channel detection signal D(t), and D(t) can be eigenvalue decomposed as where matrix Λ λ 2 = diag(λ 2 1 , λ 2 2 , . . . , λ 2 n ) is the orthogonal matrix, in which the diagonal elements are the eigenvalues of matrix R x , and the corresponding standard orthonormal eigenvector is the column vector of characteristic matrix E.
The whitening matrix (Q) can be expressed as Since the component variances of s(t) and D (t) are both 1, each element of s(t) is independent, and each element of D(t) is orthogonal; matrix C is an orthogonal normalized matrix and can be expressed as where I N is the N-order unit matrix.
The matrix, D (t), after the whitening process needs the orthogonal transformation to make aure that the components after blind source separation are not correlated and the variance is 1. The orthogonal transformation process can be expressed as where matrix P is the orthogonal matrix, s (t) is the approximate s(t) by orthogonal transformation estimation. The fourth-order cumulant of D (t) is where 1 < i, j, h, l < L, k p is the fourth-order cumulant of each element and c ip , c jp , c hp , c lp are the elements of (i, p), (j, p), (h, p), (l, p) of matrix C, respectively. For any matrix (M) with L by L order, the fourth-order cumulant matrix, Q D (M), can be expressed as where m hl is the element (h, l) of matrix M. According to Equation (9), Equation (14) is where Q S (M) is the fourth-order cumulant matrix of s(t).
Since the components of s(t) are statistically independent, matrix Q S (M) is the diagonal matrix. Select a matrix, M = [M 1 , M 2 , . . . , M r ], the orthogonal matrix (P) can be solved by joint approximate diagonalization of matrix Q D (M). The quadratic sum of diagonal elements is defined as the evaluation index of the joint approximate diagonalization calculation, The orthogonal matrix P can be solved by minimizing E(P). Then, the whitening matrix (Q) and the separation matrix (W = QP) can be obtained.

PD Signal Recovery After De-Noising
Since the signal, s (t), after blind source separation, based on the JADE algorithm, is an approximate signal vector of s(t), the amplitude of s (t) is different from the source signal. To improve the accuracy of blind source separation, this paper proposed an l 1 -norm minimization method to recover the original approximate PD UHF signal, s (t). For a matrix with M by N order, C M N solutions can be obtained in the l 1 -norm minimization [26]. Hence, the minimum norm solution can be obtained through comparison with the solutions. The l 1 -norm minimization method is shown as follows: Step 1: Calculate the C M N submatrices with M by N order of mixture matrix A, which is denoted Step 2: Solve the minimum l 1 -norm solution at t = t 0 . The solution can be expressed as Step 3: Calculate the l 1 -norm J g of solution s (g) (t 0 ), as follows Step 4: Calculate the minimum norm as the approximate signal s (t 0 ) at t = t 0 as follows Step 5: Calculate s (t) at all time points, reconstruct the PD source signal, and the signal recovery of PD signal is complete.
The proposed flow chart for the proposed PD UHF signal de-noising method based on single-channel blind source separation is shown in Figure 2. Firstly, the single-channel detected signal is analyzed through time-frequency analysis by S-transform, the frequency section characteristic parameters based on the time-frequency matrix of S-transform is used for clustering analysis, and the number of source signals (N) can be obtained. Then, construct the trajectory matrix of the detected signal, and the singular value sequence of the corresponding submatrices can be obtained by SVD calculation. The virtual channel can be reconstructed, and novel multi-channel signals can be obtained by combining the virtual channel and original single-channel signals. In addition, the PD signal can be separated from the mixed signals by BSS based on the JADE algorithm. Finally, the PD signal can be recovered by the l 1 -norm minimization method. Step 4: Calculate the minimum norm as the approximate signal s′(t0) at t = t0 as follows Step 5: Calculate s′(t) at all time points, reconstruct the PD source signal, and the signal recovery of PD signal is complete.
The proposed flow chart for the proposed PD UHF signal de-noising method based on singlechannel blind source separation is shown in Figure 2. Firstly, the single-channel detected signal is analyzed through time-frequency analysis by S-transform, the frequency section characteristic parameters based on the time-frequency matrix of S-transform is used for clustering analysis, and the number of source signals (N) can be obtained. Then, construct the trajectory matrix of the detected signal, and the singular value sequence of the corresponding submatrices can be obtained by SVD calculation. The virtual channel can be reconstructed, and novel multi-channel signals can be obtained by combining the virtual channel and original single-channel signals. In addition, the PD signal can be separated from the mixed signals by BSS based on the JADE algorithm. Finally, the PD signal can be recovered by the l1-norm minimization method.

Simulation Test Signals
Different types of typical insulation defects are simulated in the gas insulated switchgear (GIS) laboratory test platform, respectively. The corresponding PD UHF signals can be detected by the step voltage test. Firstly, the background noise interference was measured when the PD signals did not generate in GIS. Then, gradually increase the voltage until the remarkably stable PD UHF signals can be detected, the PD signals are recorded. Table 1 shows the insulation defects and test voltages. The simulation method for each insulation defect is shown in Figure 3. The simulation test platform and schematic diagram of the simulation test are shown in Figures 4 and 5, where T1 is the

Simulation Test Signals
Different types of typical insulation defects are simulated in the gas insulated switchgear (GIS) laboratory test platform, respectively. The corresponding PD UHF signals can be detected by the step voltage test. Firstly, the background noise interference was measured when the PD signals did not generate in GIS. Then, gradually increase the voltage until the remarkably stable PD UHF signals can be detected, the PD signals are recorded. Table 1 shows the insulation defects and test voltages. The simulation method for each insulation defect is shown in Figure 3. The simulation test platform and schematic diagram of the simulation test are shown in Figures 4 and 5, where T1 is the voltage regulator, T2 is the PD test transformer (YDTW-25/100, Wuhan Retop Power Equipment Company, Wuhan, China), R1 = 1 MΩ is the protection resistance, R2 = 50 Ω, C1/C2 is the power-frequency voltage divider (TAWF-1000/600), the PD UHF antenna is the Diagnostic Monitoring Systems (DMS) PD UHF antenna sensor (working bandwidth: 0.5~3 GHZ, highest antenna gain: 2.7 dBi, DMS, Glasgow, UK), the oscilloscope which is used to detect the PD UHF signals is WavePro 640Zi (maximum sampling frequency: 20 Gs/s, maximum bandwidth: 4 GHZ, Teledyne LeCroy, Chestnut ridge, New York, NY, USA. ), and N 2 /SF 6 mixed gas is filled in GIS (volume ratio of N 2 /SF 6 : 1:4, 0.5 MPa).        The PD UHF pulse signals produced in four kinds of insulation defects are combined, as shown in Figure 6a. The Gaussian white noise and periodic narrow-band interference should be added into the PD UHF pulse signals. Since the periodic narrow-band interference mainly consists of the carrier communication and mobile phone signals, which are the sine wave signals, with frequencies of 470, 900, and 1800 MHZ. The mathematical models of signals can be expressed as Since the signals of wireless communication are always from the amplitude modulation type, Z3(t) is the amplitude modulation sine wave model. In addition, the PD signal detection will be interfered by Gaussian white noise; thus, Gaussian white noise with distribution (0, 0.4) is added. Figure 6b shows the final noisy PD signal. It can be seen that the signal is seriously disturbed, and time-domain characteristics of signals cannot be identified.

De-Noising Results and Discussion
The proposed single-channel blind source separation algorithm was used to de-noise the noisy PD UHF signal in Figure 6b. Four traditional methods were also used for de-noising for comparison,  The PD UHF pulse signals produced in four kinds of insulation defects are combined, as shown in Figure 6a. The Gaussian white noise and periodic narrow-band interference should be added into the PD UHF pulse signals. Since the periodic narrow-band interference mainly consists of the carrier communication and mobile phone signals, which are the sine wave signals, with frequencies of 470, 900, and 1800 MHZ. The mathematical models of signals can be expressed as      Z 1 (t) = 0.005 sin(2π · 470 × 10 6 · t) Z 2 (t) = 0.005 sin(2π · 900 × 10 6 · t) Z 3 (t) = 0.005 sin(2π · 1800 × 10 6 · t) · cos(2π · 2 × 10 5 · t) The PD UHF pulse signals produced in four kinds of insulation defects are combined, as shown in Figure 6a. The Gaussian white noise and periodic narrow-band interference should be added into the PD UHF pulse signals. Since the periodic narrow-band interference mainly consists of the carrier communication and mobile phone signals, which are the sine wave signals, with frequencies of 470, 900, and 1800 MHZ. The mathematical models of signals can be expressed as Since the signals of wireless communication are always from the amplitude modulation type, Z3(t) is the amplitude modulation sine wave model. In addition, the PD signal detection will be interfered by Gaussian white noise; thus, Gaussian white noise with distribution (0, 0.4) is added. Figure 6b shows the final noisy PD signal. It can be seen that the signal is seriously disturbed, and time-domain characteristics of signals cannot be identified.

De-Noising Results and Discussion
The proposed single-channel blind source separation algorithm was used to de-noise the noisy PD UHF signal in Figure 6b. Four traditional methods were also used for de-noising for comparison,  Since the signals of wireless communication are always from the amplitude modulation type, Z 3 (t) is the amplitude modulation sine wave model. In addition, the PD signal detection will be interfered by Gaussian white noise; thus, Gaussian white noise with distribution (0, 0.4) is added. Figure 6b shows the final noisy PD signal. It can be seen that the signal is seriously disturbed, and time-domain characteristics of signals cannot be identified.

De-Noising Results and Discussion
The proposed single-channel blind source separation algorithm was used to de-noise the noisy PD UHF signal in Figure 6b. Four traditional methods were also used for de-noising for comparison, as shown in Table 2. For the AF algorithm, the least mean square (LMS) adaptive filter algorithm is used. The order of the filter is 255, the transition width of the filter is 100 MHZ, and the minimum width of a stop band is 50 MHZ [27]. For the AWT algorithm, the db10 wavelet is selected as the wavelet basis, and the decomposition layer is 10 [9]. For the RS algorithm, the constructed Hankel matrix is 15,000 by order 5001, and the larger number of singular values is 18 [11]. For the GSMT algorithm, the adjusted parameter of the generalized S-transform is 0.3, and the number of effective singular values is 4 [13]. Firstly, S-transform is used for the time-frequency analysis of the signal, the time waveform of the signal is shown in Figure 7a, and the time-frequency distribution is shown in Figure 7b. Then the clustering analysis method of frequency section characteristic parameters is used to analysize the time-frequency distribution and estimate the number of source signal; the clustering results are shown in Figure 7c, the number of source signals is N = 5.  Table 2. For the AF algorithm, the least mean square (LMS) adaptive filter algorithm is used. The order of the filter is 255, the transition width of the filter is 100 MHZ, and the minimum width of a stop band is 50 MHZ [27]. For the AWT algorithm, the db10 wavelet is selected as the wavelet basis, and the decomposition layer is 10 [9]. For the RS algorithm, the constructed Hankel matrix is 15,000 by order 5001, and the larger number of singular values is 18 [11]. For the GSMT algorithm, the adjusted parameter of the generalized S-transform is 0.3, and the number of effective singular values is 4 [13]. Firstly, S-transform is used for the time--frequency analysis of the signal, the time waveform of the signal is shown in Figure 7a, and the time-frequency distribution is shown in Figure 7b. Then the clustering analysis method of frequency section characteristic parameters is used to analysize the time-frequency distribution and estimate the number of source signal; the clustering results are shown in Figure 7c, the number of source signals is N = 5.  The reconstructed virtual channel signals by SVD calculation and trajectory matrix reconstitution are shown in Figure 8. Then, the de-noising PD signal was obtained by the blind source separation de-noising process based on the JADE algorithm and the l 1 -norm minimization method, as shown in Figure 9. In addition, four traditional methods in Table 2 were performed for de-noising the PD UHF signals. The de-noising results are shown in Figure 10. The reconstructed virtual channel signals by SVD calculation and trajectory matrix reconstitution are shown in Figure 8. Then, the de-noising PD signal was obtained by the blind source separation de-noising process based on the JADE algorithm and the l1-norm minimization method, as shown in Figure 9. In addition, four traditional methods in Table 2 were performed for de-noising the PD UHF signals. The de-noising results are shown in Figure 10.     Four evaluation indices are introduced to evaluate the PD UHF signal de-noising performance for each method: signal to noise ratio, root-mean-square error, waveform similarity coefficient, and variation trend parameter. Since the noise interference of the PD UHF signal in Figure 6 is very small, it can be approximately seen as the ideal non-noise original PD signal. A higher signal to noise ratio leads to a greater ability of the de-noising method to suppress interference. The root-mean-square error can show the difference between the de-noising PD and original PD signals. A lower value leads to a smaller distortion of the de-noising PD signal. When the waveform similarity coefficient is closer to 1, the de-noising PD and original PD signals are more similar. The variation trend parameter reflects the variation trend similarity between the de-noising PD and original PD signals. The variation trends for both signals are consistent when the value is closer to 1.
The evaluation results of various methods which were evaluated by the proposed indices are shown in Table 3. Four evaluation indices are introduced to evaluate the PD UHF signal de-noising performance for each method: signal to noise ratio, root-mean-square error, waveform similarity coefficient, and variation trend parameter. Since the noise interference of the PD UHF signal in Figure 6 is very small, it can be approximately seen as the ideal non-noise original PD signal. A higher signal to noise ratio leads to a greater ability of the de-noising method to suppress interference. The root-mean-square error can show the difference between the de-noising PD and original PD signals. A lower value leads to a smaller distortion of the de-noising PD signal. When the waveform similarity coefficient is closer to 1, the de-noising PD and original PD signals are more similar. The variation trend parameter reflects the variation trend similarity between the de-noising PD and original PD signals. The variation trends for both signals are consistent when the value is closer to 1.
The evaluation results of various methods which were evaluated by the proposed indices are shown in Table 3. It can be seen that the proposed de-noising method based on single-channel blind source separation can effectively suppress various noise interference without calculating the characteristic parameters of narrow-band interference, and the distortion of the de-noising PD signal is very small through the recovery method based on l 1 -norm minimization method, compared with traditional methods. For Method 2 (AF method), since the frequency of the PD UHF signal and periodic narrow-band interference signal overlap, part of the characteristic quantity of the PD signal will be removed in Method 2, leading to a large distortion of PD signals. This demonstrates that the proposed method has a good performance in suppressing the periodic narrow-band interference signal. For Method 3 (WT method), since the PD UHF signal is complicated and appropriate wavelet basis function is difficult to select, the distortion of the de-noising PD signals is apparent. In addition, the difference between the de-noising and original PD signals is larger in Method 4 (RS method), due to the unknown original PD signal. This demonstrates that, compared with the RS method, the amplitude of the de-noising PD UHF signal can be recovered by the proposed l 1 -norm minimization method. For Method 5 (GSMT method), since there is amplitude modulation sine wave periodic narrow-band interference, the submatrix modulus value method cannot be used to suppress the narrow-band interference, the suppression performance for periodic narrow-band interference is not good. Compared with the GSMT method, the proposed method can also effectively suppress the amplitude-modulation narrow-band interference. For Method 6 (UWT method), since the frequency distribution of the periodic narrow-band interference and PD UHF signals overlap, the suppressing performance is not good.

Field Test for De-Noising
The live PD detection was performed on a 500 kV open-type substation in operation, as shown in Figure 11. The devices are the same as those for the simulation test in Section 3. Figure 12 shows the detected PD UHF signal with noise interference. The PD signal is seriously interfered by the background noise and has obvious distortion. Six methods in Table 2 were performed to de-noise the detected PD UHF signal in the field test. The de-noising results are shown in Figure 13.  Since the original PD UHF signal with no noise interference cannot be obtained in the field test, the evaluation method in Section 3.2 cannot be used to evaluate the de-noising performance. Hence, the evaluation indices of the noise suppression ratio and the amplitude attenuation ratio are proposed for evaluating the de-noising performance in the field test. The noise suppression ratio denotes the suppression effect of the PD UHF signal after de-noising. The de-noising effect is positively proportional to the noise suppression ratio. The amplitude attenuation ratio denotes the attenuation degree of the PD UHF signal after de-noising. The distortion degree is positively proportional to the amplitude attenuation ratio. Table 4 shows the evaluation results of de-noising performance for various methods. The proposed de-nosing method can effectively suppress noise interference, and the amplitude attenuation of the de-noising PD signal is the smallest, compared with traditional methods. In addition, the proposed method can effectively reduce the calculation amount and computing time.   Since the original PD UHF signal with no noise interference cannot be obtained in the field test, the evaluation method in Section 3.2 cannot be used to evaluate the de-noising performance. Hence, the evaluation indices of the noise suppression ratio and the amplitude attenuation ratio are proposed for evaluating the de-noising performance in the field test. The noise suppression ratio denotes the suppression effect of the PD UHF signal after de-noising. The de-noising effect is positively proportional to the noise suppression ratio. The amplitude attenuation ratio denotes the attenuation degree of the PD UHF signal after de-noising. The distortion degree is positively proportional to the amplitude attenuation ratio. Table 4 shows the evaluation results of de-noising performance for various methods. The proposed de-nosing method can effectively suppress noise interference, and the amplitude attenuation of the de-noising PD signal is the smallest, compared with traditional methods. In addition, the proposed method can effectively reduce the calculation amount and computing time.  Since the original PD UHF signal with no noise interference cannot be obtained in the field test, the evaluation method in Section 3.2 cannot be used to evaluate the de-noising performance. Hence, the evaluation indices of the noise suppression ratio and the amplitude attenuation ratio are proposed for evaluating the de-noising performance in the field test. The noise suppression ratio denotes the suppression effect of the PD UHF signal after de-noising. The de-noising effect is positively proportional to the noise suppression ratio. The amplitude attenuation ratio denotes the attenuation degree of the PD UHF signal after de-noising. The distortion degree is positively proportional to the amplitude attenuation ratio. Table 4 shows the evaluation results of de-noising performance for various methods. The proposed de-nosing method can effectively suppress noise interference, and the amplitude attenuation of the de-noising PD signal is the smallest, compared with traditional methods. In addition, the proposed method can effectively reduce the calculation amount and computing time.

Conclusions
This paper proposed a novel PD UHF signal de-noising method, based on a single-channel blind source separation algorithm. This method utilizes the SVD calculation to convert the single-channel detected PD signals into multi-channel detected signals, and can effectively suppress

Conclusions
This paper proposed a novel PD UHF signal de-noising method, based on a single-channel blind source separation algorithm. This method utilizes the SVD calculation to convert the single-channel detected PD signals into multi-channel detected signals, and can effectively suppress the background noise in the PD UHF signals, combined with the JADE and l 1 -norm minimization methods. The conclusions are as follows: (1) The submatrix of the SVD decomposition of the original PD signal can convert the single-channel detected PD signal into multi-channel PD signals; the underdetermined problem of blind source separation can be effectively solved. (2) The l 1 -norm minimization method can effectively solve the large amplitude vibration problem after single-channel blind source separation, which is better for the subsequent signal processing and analysis. (3) Compared with traditional methods, the proposed method can effectively de-noise the Gaussian white noise and periodic narrow-band interference, and have small distortion.
The proposed method can be applied in the field partial discharge location and the insulation defect identification of HV apparatus. In addition, the field partial discharge location based on the proposed method will be studied in the future work.