Insulator Partial Discharge Localization Based on Improved Wavelet Packet Threshold Denoising and Gxx−β Generalized Cross-Correlation Algorithm

TL;DR

This paper introduces a new method to locate partial discharges in insulators using improved signal processing techniques and a specialized algorithm for better accuracy in noisy environments.

Contribution

The novel contribution is the combination of improved wavelet packet denoising and a modified Gxx−β generalized cross-correlation algorithm for PD localization.

Findings

The improved wavelet packet denoising algorithm effectively separates PD signals from noise with low distortion.

The modified Gxx−β weighting function significantly improves the accuracy of time difference estimation between sensors.

Abstract

Partial discharge (PD) in insulators will not only lead to the gradual degradation of insulation performance but even cause power system failure in serious cases. Because there is strong noise interference in the field, it is difficult to accurately locate the position of the PD source. Therefore, this paper proposes a three-dimensional spatial localization method of the PD source with a four-element ultra-high-frequency (UHF) array based on improved wavelet packet dynamic threshold denoising and the Gxx−β generalized cross-correlation algorithm. Firstly, considering the field noise interference, the PD signal is decomposed into sub-signals with different frequency bands by the wavelet packet, and the corresponding wavelet packet coefficients are extracted. By using the improved threshold function to process the wavelet packet coefficients, the PD signal with low distortion rate and high signal-to-noise ratio (SNR) is reconstructed. Secondly, in order to solve the problem that the amplitude of the first wave of the PD signal is small and the SNR is low, an improved weighting function, Gxx−β, is proposed, which is based on the self-power spectral density of the signal and is adjusted by introducing an exponential factor to improve the accuracy of the first wave arrival time and time difference calculation. Finally, the influence of different sensor array shapes and PD source positions on the localization results is analyzed, and a reasonable arrangement scheme is found. In order to verify the performance of the proposed method, simulation and experimental analysis are carried out. The results show that the improved wavelet packet denoising algorithm can effectively realize the separation of PD signal and noise and improve the SNR of the localization signal with low distortion rate. The improved Gxx−β weighting function significantly improves the estimation accuracy of the time difference between UHF sensors. With the sensor array designed in this paper, the relative localization error is 3.46%, and the absolute error is within 6 cm, which meets the requirements of engineering applications.

Figures

Full text

1. Introduction

Insulators play a vital role in power systems, and their insulation performance has an important impact on ensuring the safe and reliable operation of power systems [1,2,3]. However, the insulator may acquire defects due to aging, pollution, or damage during long-term operation, which will cause its insulation performance to decline, and the PD phenomenon to occur [4,5]. As an early signal of insulator failure, PD has an extremely important early-warning role. If no timely treatment or repair is conducted, PD may further escalate and eventually lead to insulator flashover or even breakdown [6]. Therefore, through PD detection in insulators, the defect position in insulators can be found and determined in time, which is of great significance to prevent the occurrence of insulator failure and ensure the safety of the power system [7,8,9].

At present, the commonly used methods of insulator PD detection and localization include the pulse current method, ultrasonic method and UHF method [10,11,12]. The pulse current method shows high sensitivity only in power outage measurement and is relatively rarely used in charged or inline detection [13]. Because the ultrasonic wave is strongly affected by the medium in space, attenuation is fast, and the wave speed is unstable; so, the ultrasonic method is only suitable for accurate localization in a small range [14]. The UHF method can detect the electromagnetic wave signal generated by PD through the UHF sensor, which can not only realize online detection but also be widely used because of the electromagnetic wave’s strong penetration ability and long propagation distance [15,16,17]. The localization of the PD source is mainly based on the time difference of arrival (TDOA) of the signals from different sensors. So, the accuracy of the time difference directly determines the precision of PD source localization [18]. The commonly used time difference estimation algorithms include the threshold method, energy accumulation method, cross-correlation method, and so on. The threshold method takes the moment when the signal exceeds the set threshold for the first time as the first wave arrival moment of the PD signal. This method is simple in principle and convenient in operation, but the influence of noise and the selection of the threshold value will affect the judgment accuracy [19]. The energy accumulation method generates the energy accumulation curve and takes the inflection point of the accumulation curve as the starting moment of PD. Although this method can avoid the interference of small-amplitude noise, the selection of the inflection point often does not conform to a unified standard and is greatly influenced by human factors [20]. The cross-correlation method, which calculates the correlation function between two signals and considers the moment corresponding to the peak point as the time difference between the two signals, is currently the most widely used method [21,22]. However, in practical application, due to the large interference noise in the field, and even the fact that the PD signal is completely drowned by the noise, both the threshold method and the energy accumulation method cannot accurately obtain the characteristic moments of the PD signal. Although the cross-correlation method can suppress the influence of noise to a certain extent, there might be more fake peaks or the main peak might not be obvious when the noise is large, which leads to a wrong judgment of the peak point, thus increasing the difficulty of PD positioning.

Therefore, this paper proposes a three-dimensional spatial localization method with a four-element UHF array based on improved wavelet packet threshold denoising and the $[eqn]$ generalized cross-correlation algorithm. For the interference of field noise, wavelet packet technology is used to extract the wavelet packet coefficients of the corresponding sub-signals in each frequency band of the PD signal, and the wavelet packet coefficients are processed and reconstructed using an improved threshold function to obtain the denoised PD signal. On the premise of reducing the waveform distortion of useful signals, the method can effectively suppress the noise in the PD signal and improve the SNR of the PD signal. To solve the problem of inaccurate time difference estimation of the PD signal, an improved weighting function, $[eqn]$ , is proposed. The function adjusts the self-power spectral density of the signal by introducing an exponential factor, which solves the problem that it is difficult to accurately calculate the signal arrival moment and time difference due to the small amplitude and low SNR of the PD signal. In addition, this paper also analyzes the influence of different sensor array shapes and the position of the PD source on the localization results, finds a reasonable arrangement scheme, and finally achieves spatial position localization with a four-element three-dimensional rectangular array.

2. Improvement of the Insulator PD Localization Method

2.1. Improved Wavelet Packet Threshold Localization Signal Denoising Technology

The PD signal detected by a UHF sensor array is easily interfered by various noises in the complex noise environment in the field. These noises may come from the operation of surrounding equipment, electromagnetic interference, or other random noise signals. Because the amplitude of the PD signal is usually small, and the duration is short, the presence of noise further masks the PD characteristics, thus reducing the sensitivity and accuracy of detection. Therefore, it is necessary to design appropriate denoising technology to improve the reliability of the signal [23,24,25].

Wavelet transform is a mathematical tool commonly used in signal processing that can effectively analyze and process signals with non-stationary characteristics by breaking them down into components of different frequencies and time scales, which makes it widely used in the field of PD localization. The wavelet denoising process is shown in Figure 1.

The traditional wavelet transform usually provides only limited frequency resolution, while wavelet packet decomposition can explore the frequency content of a signal more deeply. The wavelet packet decomposition formula is:

$$[eqn]$$

where $[eqn]$ represents the high-pass filter coefficients corresponding to the orthogonal wavelet function $[eqn]$ , and $[eqn]$ the low-pass filter coefficients corresponding to the orthogonal scale function $[eqn]$ , satisfying $[eqn]$ , with L being the filter length. The recursively defined function $[eqn]$ , n = 0, 1, 2, … is called the wavelet packet and is determined by the orthogonal scaling function $[eqn]$ . The effect of wavelet packet decomposition is shown in Figure 2.

In the process of denoising using wavelet packet technology, the selection of the threshold function is extremely critical. There are two traditional threshold functions:

Hard threshold function

$$[eqn]$$

Soft threshold function

$$[eqn]$$

where $[eqn]$ is the wavelet coefficient obtained by decomposition, $[eqn]$ is the wavelet coefficient after threshold processing, $[eqn]$ is the threshold—usually a fixed threshold $[eqn]$ is used— $[eqn]$ is the variance of noise, and $[eqn]$ is the length of the signal.

It is not difficult to see that the hard threshold function is discontinuous at the piecewise point, which will cause the signal to be easily shaken when reconstructed. Although the soft threshold function is continuous at the piecewise point, the wavelet coefficients larger than the threshold undergo a constant linear contraction, which causes the signal to be deformed to a certain extent during processing, possibly introducing bias. In view of the problems existing in the hard threshold function and soft threshold function, this paper improves the threshold function and proposes a new adjustable threshold function, which is defined as follows:

$$[eqn]$$

where $[eqn]$ is the introduced regulatory factor. From the point of view of mathematics, $[eqn]$ , namely, the improved threshold function, is continuous in piecewise points, overcoming the shortcoming of the discontinuity of the hard threshold function at piecewise points. At the same time, the function satisfies the equations

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

From Equations (5)–(7), it can be seen that when $[eqn]$ , the output is subtracted from a fixed threshold, as with the soft threshold function, and when $[eqn]$ , the output remains unchanged and completely passes, as with the hard threshold function. When $[eqn]$ , the output amplitude is subtracted twice the threshold, which makes the processing of high-amplitude signals more aggressive. In practical application, the appropriate $[eqn]$ can be selected according to the characteristics of the signal to achieve the purpose of flexible application of the threshold function. Meanwhile, when $[eqn]$ , the deviation between the wavelet coefficient processed by the improved threshold function and the original wavelet coefficient is as follows:

$$[eqn]$$

and satisfies the equation:

$$[eqn]$$

When the parameter $[eqn]$ is constant, as $[eqn]$ increases, $[eqn]$ decreases, and when $[eqn]$ , $[eqn]$ . It can be seen that the improved threshold function overcomes the shortcoming of linear contraction of the wavelet coefficient when the soft threshold function is larger than the threshold value, so as to achieve a better signal denoising effect. A comparison of the threshold function curves is shown in Figure 3.

Accurate PD source localization is based on a clear signal foundation. By using the improved threshold function to denoise the PD signal, not only the interference of the surrounding noise signal is reduced and signal recognizability is improved, but also a good signal basis for time delay estimation is provided, thus improving the accuracy of time delay estimation, which directly determines the accuracy of the localization of the PD source. Thus, the improved threshold function has an important influence on the accuracy of the localization of the PD source.

2.2. Improved Time Delay Estimation Method for PD Localization Signals

Time delay estimation is one of the key factors to determine the localization accuracy of a PD source. Due to the small amplitude of the first wave and the low SNR of the UHF signal, the accuracy of the first wave arrival time and the time difference estimation is greatly affected. Therefore, it is necessary to research and apply an efficient time delay estimation algorithm to improve the localization accuracy of PD signals [26,27].

Suppose that the PD signals received by UHF sensor No. 1 and UHF sensor No. 2 are, respectively

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ and $[eqn]$ are signal attenuation coefficients, $[eqn]$ indicates the PD signals, $[eqn]$ and $[eqn]$ are non-interfering white noise signals, and $[eqn]$ and $[eqn]$ represent the time for the PD signals to spread to sensor 1 and 2, respectively, then $[eqn]$ is the time delay. The cross-correlation function of $[eqn]$ and $[eqn]$ is

$$[eqn]$$

By substituting $[eqn]$ , $[eqn]$ , Equation (12) can be turned into

$$[eqn]$$

If $[eqn]$ and $[eqn]$ , $[eqn]$ are unrelated, $[eqn]$ , $[eqn]$ , $[eqn]$ ; so, Equation (13) can be turned into

$$[eqn]$$

From the property of the correlation function $[eqn]$ , when $[eqn]$ , $[eqn]$ reaches its maximum, that is, $[eqn]$ also reaches its maximum. In this case, the $[eqn]$ value is the time delay. Therefore, the time delay for the two signals can be obtained by detecting the peak value of $[eqn]$ . Although the calculation amount of this method is small and the implementation is simple, the methos can easily be affected by noise and reverberation in practice, resulting in a difficult to detect peak value of $[eqn]$ , which makes detection difficult.

A generalized cross-correlation time delay estimation algorithm is proposed to overcome the shortcoming of basic cross-correlation in time delay estimation. By weighting the PD signal in the frequency domain, the interference of noise and reverberation can be effectively suppressed, and the actual signal can be highlighted, thus sharpening the peak value of the correlation function and improving the accuracy of time delay estimation. The time delay estimation process under the generalized cross-correlation algorithm is shown in Figure 4.

According to the Wiener–Hinchin theorem, the cross-correlation function and its cross-power spectral density are Fourier transform pairs of each other; then, the cross-correlation function of $[eqn]$ and $[eqn]$ can be expressed as

$$[eqn]$$

where $[eqn]$ is the cross-power spectral density function of $[eqn]$ and $[eqn]$ , whose mathematical expression is as follows:

$$[eqn]$$

$[eqn]$ represents the conjugate of $[eqn]$ . Therefore, the generalized cross-correlation function of the signals $[eqn]$ and $[eqn]$ can be expressed as

$$[eqn]$$

In Equation $[eqn]$ , $[eqn]$ is the generalized cross-correlation weighting function. According to prior knowledge of signal and noise, different weighting functions can be set to achieve the signal enhancement effect for different noise and reverberation situations. In practical applications, the selection of the weighting function is difficult and is also of key importance to achieve accurate time delay estimation. Conventional weighting functions include ROTH weighting, PHAT weighting, etc. However, when the signal energy is small, and the noise is relatively large, the denominator of the conventional weighting function will tend to zero, thus increasing the error. In order to suppress noise interference and sharpen the main peak of the inter-correlation function so to improve the accuracy of delay estimation, this paper proposes the improved weighting function

$$[eqn]$$

The influence of each parameter in Equation (18) on the localization of the PD source is as follows:

• (1) $[eqn]$ and $[eqn]$ represent the self-power spectral density of the signals $[eqn]$ and $[eqn]$ , respectively, which responds to the distribution of the energy of the signals at each frequency, and the frequency bands where the energy is concentrated usually correspond to the eigenfrequencies of the PD source. By analyzing the self-power spectral density, the frequency bands with higher SNR can be selected to reduce the influence of interference on PD localization.
• (2) $[eqn]$ and $[eqn]$ are the frequency-domain amplitudes of the two signals, respectively, and their sum is used as the denominator to adjust the overall weighted frequency-domain characteristics. This means that the difference in amplitude between the two signals will be balanced. The larger amplitude signal will not dominate the results, while the smaller amplitude signal will have some influence, thus improving the accuracy of the time delay estimation.
• (3)The exponential adjustment factor $[eqn]$ is used to control the degree of attenuation and sensitivity of the weighting function $[eqn]$ . A smaller value of $[eqn]$ will make the weighting function more sensitive to amplitude changes, while a larger value of $[eqn]$ will weaken the response to amplitude changes. By adjusting $[eqn]$ , the influence of certain frequency components in the frequency domain can be enhanced or suppressed, thus improving the sharpness of the peaks and the accuracy of the time delay estimation.

The $[eqn]$ weighting function balances the intensity differences of the signals by using the exponential adjustment factor $[eqn]$ to regulate the energy contribution of different frequency bands, emphasizes the information of representative frequency bands, and enhances the sensitivity and anti-jamming ability of the localization algorithm to improve the time delay estimation between the signals.

2.3. Mathematical Model of Time Difference-Based UHF PD Localization

The basic principle of UHF PD localization is to receive the signal of the PD source through the UHF sensors and use the time difference for the signal arriving at each sensor to locate the PD source. Taking the UHF sensors No. 1 and No. 2 as an example, the received signal and time delay diagram are shown in Figure 5.

Assume that the coordinates of the PD source are $[eqn]$ , the coordinates of the four UHF sensors are $[eqn]$ and $[eqn]$ , respectively. Considering the spatial geometric distance, we have

$$[eqn]$$

where $[eqn]$ is the speed of electromagnetic wave propagation in air, whose value is $[eqn]$ , $[eqn]$ is the time difference between the time when the $[eqn]$ sensor receives the signal and the time when the first sensor receives the signal, and $[eqn]$ represents the distance from the PD source to the $[eqn]$ sensor, namely:

$$[eqn]$$

By solving Equation (19), the position of the PD source can be obtained, thus realizing the localization of the PD source. The principle of PD source localization is shown in Figure 6.

3. Simulation Validation of the Improved Insulator PD Localization Method

3.1. Simulation Validation of the Improved Wavelet Packet Thresholding for the Denoising Algorithm of PD Signals

Research has shown that the waveform features of the PD signal mainly include exponential attenuation and exponential oscillation attenuation. In this paper, we chose to use double exponential oscillation attenuation to simulate PD signals. Its mathematical expression is as follows:

$$[eqn]$$

where $[eqn]$ is the amplitude of the PD signal, $[eqn]$ is the attenuation coefficient, $[eqn]$ is the central oscillation frequency, and $[eqn]$ is the signal duration, with bandwidth of 300 MHz–800 MHz and system sampling rate of 10 GSa/s.

Considering the complexity of the noise environment around the insulators, Gaussian white noise with SNR = −5 dB was added to the original signal to construct a noise interference environment close to reality. Figure 7 shows the time-domain waveform and frequency spectrum characteristics of the signal after adding noise.

Wavelet packet technology was used to decompose the generated PD signal. The number of the decomposition layers was three, and the wavelet packet coefficient obtained by decomposition is shown in Figure 8.

As can be seen from Figure 8, the wavelet packet coefficients of useful signals were mainly concentrated in nodes (3.1), while those of noisy signals were concentrated in other nodes. The obtained wavelet packet coefficients were processed using the hard threshold function, soft threshold function, and improved threshold function, and the signal was finally reconstructed. The reconstructed signal is shown in Figure 9.

In order to further demonstrate the effect of the improved threshold function, the SNR, root-mean-square error (RMSE), and normalized correlation coefficient (NCC) were used to evaluate the effect. They are defined as:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where $[eqn]$ represents the original signal, $[eqn]$ represents the signal after noise removal, and $[eqn]$ represents the length of the signal. The SNR, RMSE, and NCC for different denoising threshold functions are shown in Figure 10.

As can be seen from Figure 10, the SNR and NCC of the signal after the improved threshold processing were higher than those obtained using the hard and soft thresholds, and the RMSE was also reduced compared to that obtained with the hard and soft thresholds; so, the improvement was verified.

3.2. Simulation Validation of the Improved Time Delay Estimation Method for PD Signal Localization

The PD signal with noise with an SNR of −5 dB was subjected to a 10 ns time delay, and the signal delay with an SNR of 0 dB and 5 dB is also shown in Figure 11.

In order to verify the effectiveness of the $[eqn]$ weighting function designed in this paper, different weighting functions were used to estimate the time delay signals under different SNRs, and the results are shown in Figure 12.

The time delay estimation results of signals with different SNRs under different weighting functions are shown in Table 1.

From Figure 12 and Table 1, it can be seen that when the SNR was −5 dB and 0 dB, the time delay estimation under ROTH weighting and PHAT weighting was not as accurate as under $[eqn]$ weighting. When the SNR was 5 dB, although ROTH weighting could approximate the delay, and PHAT weighting could also accurately estimate the delay, the peak highlighting effect was not as obvious as with $[eqn]$ weighting.

The extent to which different values of $[eqn]$ affected the time delay estimation is shown in Figure 13.

As can be seen in Figure 13, in the range of $[eqn]$ ∈[0,1], as the exponential adjustment factor $[eqn]$ decreased, although the pseudo-peak amplitude around the peak point gradually decreased, the number of pseudo-peaks gradually increased, which made it difficult to accurately judge the peak point. When $[eqn]$ ∈[0.4,0.8], the effect was more satisfactory, and $[eqn]$ was set to 0.6 in this paper. Therefore, the $[eqn]$ weighting function proposed in this paper could improve the accuracy of delay estimation and ultimately the accuracy of PD source localization by adjusting the exponential adjustment factor $[eqn]$ , based on overcoming the defects of the traditional weighting functions.

3.3. Influence of Array Shape and PD Source Position on the Localization Result

In practical applications, the shape of the sensor array and the position of the PD source have important effects on the reception and processing of the signal. By designing sensor arrays with different shapes, this paper systematically discusses the influence of different array configurations on signal positioning accuracy and performance.

We assumed that the position of the PD source was $[eqn]$ , and the placement coordinates of the four UHF sensors were $[eqn]$ and $[eqn]$ . Taking the A1 sensor as a reference, the four sensors were on the same plane and formed a Y shape, as shown in Figure 14.

In this case, the time delay between the arrival of the PD signal to the reference sensor and to the A2, A3, and A4 sensors was $[eqn]$ , respectively. According to a theoretical numerical calculation, $[eqn]$ .

According to the time delay results obtained by theoretical calculation, the corresponding four groups of PD signals were simulated, and the SNR was −5 dB. The improved threshold function in this paper was used to denoise the signal, and the $[eqn]$ generalized cross-correlation algorithm was used to obtain the time difference between the reference array and the other array elements. The three groups of time differences obtained were substituted into Equation (19). The localization results obtained are shown in Table 2.

As can be seen from Table 2, in the five groups of simulation data, the average time delay for the three groups obtained with the Y-type array was $[eqn]$ and $[eqn]$ . The average position of the localization was (1.18 m, 1.51 m, 1.58 m), and the relative localization error was 4.85%. It can be seen that in the case of the Y-shaped array, the localization results on the x-axis and y-axis were more accurate, while the localization results on z-axis showed relatively large errors.

The position of the PD source, along with the coordinates of sensors A1, A2, and A3, remained unchanged, and the position’s coordinates of the A4 sensor became (0 m, 2.0 m, 0 m). In this case, the four sensors were still on the same plane and formed a rectangle, as shown in Figure 15.

In this case, the theoretical time delay between the PD source and the A2, A3, and A4 sensors was $[eqn]$ , $[eqn]$ , and $[eqn]$ , respectively. Similarly, the corresponding four groups of PD signals were simulated according to the theoretical time delay, and then the method described in this paper was used for denoising and time delay estimation. The final localization results are shown in Table 3.

As can be seen from Table 3, in the five groups of simulation data, the average time delay of the three groups obtained under the rectangular array was $[eqn]$ , $[eqn]$ , and $[eqn]$ . The mean value of the localization results was (1.17 m, 1.52 m, 1.72 m), and the relative localization error was 2.76%. It can be seen that in the case of a rectangular array, the localization results on the x-axis and y-axis were not much different from the localization results obtained with the Y-type array, but the localization on the z-axis as better than that obtained with the Y-type array.

In order to further improve the accuracy of z-axis localization, the A2 sensor in the above rectangular array was raised by 0.3 m, and the coordinate position of the A2 sensor was set to A2(0 m, 2.0 m, 0.3 m). The position of the PD source and other sensor coordinates remained unchanged. In this case, the four sensors were distributed in three-dimensional space, and the projection on the xoy plane was still rectangular, as shown in Figure 16.

In this case, the theoretical difference between the arrival times of the PD signal to the A1 sensor and the A2, A3, and A4 sensors was $[eqn]$ , $[eqn]$ , and $[eqn]$ , respectively. Similarly, four groups of corresponding PD signals were generated, and then the proposed method was applied for denoising and time delay estimation. The final localization results are shown in Table 4.

As can be seen from Table 4, in the five groups of simulation data, the average time delay of the three groups obtained under the three-dimensional rectangular array was $[eqn]$ , $[eqn]$ , and $[eqn]$ , the average value of the localization results was (1.16 m, 1.53 m, 1.80 m), and the relative localization error was 1.9%. It can be seen that in the case of the three-dimensional rectangular array, the localization results on the x- and y-axis were still relatively accurate, and not only the localization accuracy on the z-axis was further improved, but also the overall relative localization error was lower.

On the premise that the three-dimensional rectangular array remained unchanged, the coordinate positions of the PD source were changed to P1 (1.2 m, 1.5 m, 2.8 m), P3 (1.2 m, 1.5 m, 0.8 m), P4 (2.2 m, 1.5 m, 1.8 m), and P5 (3.2 m, 1.5 m, 1.8 m), and the effect of the localization of the PD source at different positions in the array was explored. The position of the PD source is shown in Figure 17.

In the same way as above, the theoretical time delay between each PD source position and each sensor was calculated, and according to the theoretical time delay, the corresponding four groups of PD signals were simulated; then, the method described in this paper was used to denoise and estimate the time delay. Combined with the position of the PD source at P2 (1.2 m, 1.5 m, 1.8 m), the localization results for different PD source positions are shown in Table 5.

As can be seen from Table 5, when the horizontal position of the PD source remained unchanged, the localization accuracy became higher and higher as the PD source gradually approached the sensor array. As the PD source moved away from the sensor array, the localization accuracy decreased.

As can be seen from Table 2, Table 3, Table 4 and Table 5, in a certain area, the localization with the rectangular array with plane projection was better than that of the Y-shaped array. Appropriately raising the height of a sensor in the array to form a three-dimensional spatial distribution could significantly improve the localization accuracy. The smaller the distance between the PD source and the sensor array, the higher the localization accuracy.

4. Insulator PD UHF Localization Experiment

In order to test the effectiveness of the three-dimensional spatial insulator PD localization based on wavelet packet decomposition and the $[eqn]$ generalized cross-correlation algorithm proposed in this paper, a PD localization platform was built in the laboratory. The localization platform mainly consisted of a test model, a high-voltage testing platform, and a UHF localization device. The test site is shown in Figure 18.

The test model mainly consisted of a UHF sensor, an MTESB-10 kVA/100 kV high-voltage experimental transformer, an FXBW4-110 composite insulator, a signal conditioner, a four-channel high-sampling-rate oscilloscope, and a computer.

The positions of the four UHF sensors were A1 (0.97 m, 0 m, 0 m), A2 (0.97 m, 0.49 m, 1.08 m), A3 (1.91 m, 0 m, 1.08 m), and A4 (1.91 m, 0 m, 0 m). A UHF sensor with integrated design, matching impedance of 50Ω, bandwidth of 300 MHz–800 MHz, a simple structure, a compact size, high sensitivity, and anti-interference ability was used. The oscilloscope analog bandwidth was 0–2 GHz, the sampling rate was 5 GSa/s, the number of sampling points was 2500, and the length of the feed line connected to the signal conditioner and the oscilloscope was the same.

The discharge level was detected without defects on the insulator, and when the voltage was increased to 50 kV, no PD signal was detected. After setting a discharge defect at P (1.41 m, 1.79 m, 1.12 m) on the high-voltage end of the insulator, the experimental voltage also slowly increased, and when the voltage increased to 28 kV, the UHF sensor detected a large PD signal. The PD signals of the four channels are shown in Figure 19.

The collected PD signal was denoised using the hard threshold function, soft threshold function, and the threshold function improved in this paper, and the denoised PD signal is shown in Figure 20, Figure 21 and Figure 22.

As can be seen from Figure 20, Figure 21 and Figure 22, compared with the hard and soft threshold functions, the PD signal after denoising with the improved threshold function had not only significantly improved recognizability, but also significantly improved SNR, which provides favorable conditions for the subsequent time delay estimation, and thus improves the accuracy of the final PD source localization.

The denoised reference signal and the remaining three groups of signals were used to calculate the time delay using the ROTH weighting function, the PHAT weighting function, and the $[eqn]$ weighting function, and the results are shown in Figure 23, Figure 24 and Figure 25.

The results of time delay estimation and localization with different weighting functions are shown in Table 6.

As can be seen from Table 6, the improved weighting function in this paper was better than the traditional weighting functions, both in terms of time delay estimation and final localization error. The final localization results with the improved weighting function described in this paper were (1.44 m, 1.73 m, 1.17 m), the absolute localization errors on the x-, y- and z-axis were 3 cm, 6 cm, and 5 cm respectively, and the relative localization error was 3.46%. The localization results were more accurate than with the traditional weighting functions, and the localization accuracy was improved, allowing for determining the position of the PD source.

To further verify the localization performance, multiple repetitions of the localization experiments were carried out using the improved $[eqn]$ weighting function described in this paper as well as the traditional weighting functions, and the error of the localization results is shown in Figure 26.

As can be seen from Figure 26, in the multiple repetitive localization experiments, the $[eqn]$ weighting function proposed in this paper still had higher localization accuracy compared with the ROTH weighting function and PHAT weighting function, and the range of fluctuation of the localization error was small, which is reliable to a certain extent. The method described in this paper was further verified to meet the requirements of engineering applications.

5. Conclusions

In this paper, an insulator PD localization method based on improved wavelet packet threshold denoising and the $[eqn]$ generalized correlation algorithm is proposed. The effectiveness of this method was verified by theoretical derivation, simulation analysis, and laboratory field tests, and the following conclusions were obtained.

• (1)In an environment with complex noise, the improved wavelet packet threshold function overcomes the shortcomings of the traditional soft and hard threshold functions in denoising and effectively suppresses the noise components in the PD signal under the premise of reducing the distortion of the useful signal waveform, which improves the SNR and recognizability of the PD signal.
• (2)Considering the problems of large time delay estimation errors caused by small first-wave amplitude and low SNR of the PD signals, the $[eqn]$ weighting function can better highlight the characteristics of the signal by adjusting the self-power spectrum of the signal, effectively suppress the influence of noise, and improve the accuracy of time delay estimation in complex environments.
• (3)When the shape of the sensor array is fixed, the sensor array with a rectangular plane projection is more accurate than the sensor array with a Y-shape projection. Appropriately increasing the height of one of the sensors can improve the accuracy of PD source localization. The closer the PD source is to the sensor array, the better the localization.
• (4)Field tests showed that the method described in this paper can realize the accurate localization of insulator PD sources with the relative localization error of 3.46% and the absolute error within 6 cm. It meets the requirements of engineering applications.