1. Introduction
Partial discharge (PD) detection and location systems are indispensable for performing quality assurance and fault identification in high-voltage apparatus. PD can be detected and located through PD signals (electromagnetic waves), especially the ultra-high-frequency (UHF) component of PD signals [
1,
2,
3]. In recent years, PD detection and location methods based on UHF signals have been popular, with the advantages of strong anti-interference, high sensitivity and a stable transmitting speed. Thus, UHF detection methods, especially those using an omnidirectional UHF sensor array to detect and locate the PD signals in substations, have obtained fruitful results [
4,
5,
6,
7,
8,
9].
In some situations, the omnidirectional UHF sensor array is not a good choice because it need a small and light sensor array. For example, the intelligent inspection robot used for the substation inspection system can only carry a small sensor array [
10,
11]. However, PD detection methods using the time difference of arrival (TDOA) algorithm based on an omnidirectional UHF sensor array require a relatively large size sensor array to obtain an accurate TDOA. Therefore, a directional antenna array is a new choice in this situation. As a tentative research direction, the directional multiple signal classification (Dir-MUSIC) algorithm for direction of arrival (DOA) estimation based on signal strength, as represented by antenna gain array manifold, is proposed [
12]. A miniaturized spiral directional antenna array has been designed and developed to receive the PD signal. Using the signal strength information, the Dir-MUSIC algorithm can successfully and precisely find the direction of the PD source. However, in the Dir-MUSIC algorithm, the antenna gain is assumed to be approximatively equal at different frequencies of interest. In fact, when the antenna gain at different frequencies is too different, the direction error will be unacceptable. Therefore, it is necessary to improve a wideband Dir-MUSIC algorithm.
It is difficult to deal with the wideband signals directly using a subspace DOA estimation algorithm. Most wideband subspace methods decompose the sensor outputs into several narrowband signals using a filter bank or the discrete Fourier transform (DFT). Normally, they nclude the incoherent signal subspace method (ISSM) [
13] and the coherent signal subspace method (CSSM) [
14]. ISSM is one of the simplest wideband methods. It process the decomposed narrowband signals independently, applying narrowband techniques, and the results of all the frequency bins are averaged to obtain the final DOA estimation. CSSM calculates a focusing matrix between different frequency bins and estimates a coherently averaged sample covariance matrix, which can be obtained by applying narrowband techniques. CSSM requires initial values to find the focusing matrices, and the estimation performance of CSSM is sensitive to these initial values.
Wideband subspace techniques have many applications. An efficient implementation of CSSM on the field-programmable gate array (FPGA) has been proposed to achieve the desired performance [
15]. A generic subspace model has been proposed for characterizing a variety of radio frequency interference (RFI) types, which reveals a low-dimensional structure of the RFI subspace [
16]. For PD location in substations in oil, a PD location method based on an ultrasonic phased array and wideband array signal processing is proposed [
17].
In this paper, based on the Dir-MUSIC algorithm, we propose a novel wideband Dir-MUSIC algorithm using the strength proportion (DirSP), which can be regarded as another form of CSSM. In the proposed algorithm, PD signals will be decomposed into several narrowband signals using the DFT. We will choose some frequency bins to estimate an averaged sample covariance matrix at a certain frequency, and we will first estimate a focusing PD signal at this frequency. Then, we could apply the Dir-MUSIC algorithm on the averaged sample covariance matrix and the focusing PD signal.
2. Wide-Band Dir-MUSIC Algorithms
A PD signal is always a microsecond pulse signal, and the time interval between two pulses occurs randomly at the interval 50 us, 1 ms. Therefore, the probability of multi-source signal pulse superposition is very small, and this paper can only consider a single PD source in one time. We consider a uniform circular array of N elements (directional antenna), and the opening direction of each element is toward outside. The bandwidths of PD signals need not be identical, but there should be some frequency band where PD signals’ bandwidths overlap.
Most wideband subspace methods decompose the sensor outputs into several narrowband signals using a filter bank or the DFT. Generally, PD signals are always measured at some frequency band of interest, such as the designed resonance points of antennas. If
contains all the frequency bands of interest, then the output of the filter bank or DFT module can be written in vector form at
K frequencies, as follows:
where
where
is the output of the
jth antenna at frequency
for
,
and
.
is the component of the PD source and
is the antenna gain matrix at frequency
(
is a vector in this paper):
where
is the
jth antenna gain when the direction angle of the PD source is
at frequency
, and
holds.
2.1. The Normal Incoherent and Coherent Method
Next, a normal incoherent method is introduced. First, find the direction using the Dir-MUSIC algorithm at frequency
, as follows:
Then take an arithmetic mean of
as the estimated direction:
In the coherent method, an estimation of the coherently averaged sample covariance matrix is indispensable. Focusing matrices are effective in this operation, and can transform the constituent narrow-band components of the received signal into appropriate narrow-band representations at a certain given frequency. Next, we will introduce the rotational signal subspace (RSS) focusing matrix
, which satisfies the following constrained minimization problem:
subject to
where
is the Frobenius matrix norm, and
is the centre frequency. One solution to (
6) is
where the columns of
and
are the left and right singular vectors of
.
Then, the coherently averaged sample covariance matrix at
can be estimated by
It must be noted that directly applying the RSS focussing matrix is not a good choice.
and
in (
6) always have different modules, so
T can not have a theoretical solution that satisfies
. In the traditional CSSM, the steering vectors have the same modules, and this characteristic ensures the effectiveness of the traditional CSSM. Therefore, when the antenna gain vectors have quite different modules, CSSM in this paper will lose accuracy.
As an alternative, a wideband Dir-MUSIC algorithm using the strength proportion is proposed in the next subsection.
2.2. The Wide-Band Dir-MUSIC Algorithm Using the Strength Proportion
In this paper,
is only a real number when
is given. Then, from (
1) and (
3), we can obtain
Let
be the reference frequency; (
9) could be written as
Let
be the sum of the left side in (
10):
Then, the coherently averaged sample covariance matrix at
can be estimated by
where
.
Next, we will explain that this method is equivalent to CSSM when it is applied in the normal array signal processing.
We consider a single source and a linear array, and the incident angle is
. Then (
9) and (
10) convert to
The coherently averaged sample covariance matrix of this method and the CSSM are identical. Therefore, the proposed method is an application of CSSM in this paper.
Now, we can estimate the direction of the PD source using the Dir-MUSIC algorithm, and the steps are as follows.
We calculate the eigenvalues and eigenvectors of
via eigenvalue decomposition. It is clear that
is a real symmetric matrix, and the eigenvalues are all real numbers.
also consists of signal subspace
and noise subspace
:
where
is a simple expression of
. Therefore,
is a positive definite. Since
are the eigenvalues, and
are the related corresponding eigenvectors. Additionally, the eigenvectors are orthogonal to each other; that is,
This paper considers only one signal at one time; therefore, there is only one eigenvalue related to the signal. Certainly,
is the biggest eigenvalue related to the signal, and suppose
is one of the smaller eigenvalues. Then, we have
Because
is positive definite, we have
From (
3), we can see that
is also a vector function of
; therefore, we can search the direction of the PD source by traversing the possible value of
. Let
eigenvectors be related to the smaller
eigenvalue to construct a noise matrix:
Let the space spectrum be:
where
represents
.
When we take the maximum value of , it indicates that is closest to 0 in all ; then, we decide this is the estimated direction angle.
3. Simulations
In this section, simulations will be carried out to verify the advantages of the wideband Dir-MUSIC algorithm. A relatively simple but effective wideband PD signal is designed for simulations.
A double exponential oscillation attenuation function is used to simulate the PD signal [
18], and the specific expression is
where the parameters can be
. This is a narrow-band signal in the approximate frequency band
. Consider other sets of parameters
and
; the related simulated PD signals are written as
and
, and their approximate frequency bands are
and
. Meanwhile, the simulated PD signal related to the first set of parameters is written as
. Then, a wide-band PD signal is simulated as
The sampled pulse signal of
without noise at the sample rate
is shown in
Figure 1, and the frequency spectrum obtained from DFT is shown in
Figure 2.
From
Figure 1 and
Figure 2, this simulated PD signal is a wideband signal that contains three bands.
,
,
are the centre frequencies of the three bands. When the antenna gain is different among
,
and
, the Dir-MUSIC algorithm will lose its effectiveness.
Based on the actual antenna pattern, a linear combination of three Gaussian functions is enough to describe the two-dimensional antenna pattern function, as shown in
Suppose the set of parameters at
is
, and the other two set sof parameters are different at
and
. The normalized antenna gains of the three bands are shown in
Figure 3.
We provide simulation steps in
Figure 4.
Let the direction of the PD source
change from
to
, and SNR is set at 5, 0, −5, and the direction results of DirSP compared to Dir-MUSIC are shown in
Figure 5,
Figure 6 and
Figure 7.
From
Figure 5,
Figure 6 and
Figure 7, we can see when the antenna pattern functions among different frequency bins are quite different; the Dir-MUSIC algorithm only using one function will lose accuracy. Meanwhile, the DirSP algorithm using three functions performs very well.
Let the direction of the PD source
change from
to
; SNR is set at 5, 0,
, and the direction results of the wide-band methods are shown in
Figure 8,
Figure 9 and
Figure 10. Additionally, each result in the following figures is a mean of 50 random estimations.
From
Figure 8,
Figure 9 and
Figure 10, we can see that the errors of ISSM and DirSP are similar, but much less than CSSM. Specifically, the mean errors of ISSM, CSSM and DirSP are
,
and
at SNR
,
,
and
at SNR
,
,
and
at SNR
. This is ineffective for CSSM when directly applied in the wide-band Dir-MUSIC algorithm. On the other hand, errors in ISSM and DirSP are acceptable even when SNR is
.
It seems that ISSM is a little better than the proposed DirSP. The reason for this is that measurements are generated from the ideal model perfectly, and every Dir-MUSIC estimation is a perfect estimation of ISSM. Additionally, in DirSP, the strength proportion in (
10) is always not equal to 1, and the noise distribution changes in
in (
12). Therefore, ISSM performs best in the ideal model. However, the relationship between the actual measurements and the measured antenna pattern are not as ideal. Moreover, the proportion of
and
is always not equal in (
23); having the wrong proportion will cause mistakes in direction estimation. Next, we simulate a new PD signal as
From
Figure 11,
Figure 12 and
Figure 13, we can see that ISSM errors become higher when the proportion of frequency bins is wrong. Specifically, the mean errors of ISSM, CSSM and DirSP are
,
and
at SNR
,
,
and
at SNR
,
,
and
at SNR
. Another method could possibly be used to estimate the proportion; however, this will increase computational complexity and the accuracy will rely on the proportion estimation. Fortunately, errors in the new method are not influenced with changes in proportion and the errors of DirSP are quite acceptable.
At last, we should analyze why CSSM performs badly in the simulations.
and
in (
6) always have different modules, so
T could not have a theoretical solution. In the traditional CSSM, the steering vectors have the same modules, and this characteristic ensures the effectiveness of the traditional CSSM. Therefore, when the antenna gain vectors have quite different modules, CSSM in this paper will lose accuracy.
4. Experimental Data Processing
Experiments were carried out to verify the effectiveness of DirSP. The experimental platform consisted of a digital oscilloscope with storage function, four equal-length radio frequency coaxial cables, the developed uniform fan-shaped directional Vivaldi antenna array, and a lighter as a signal generator.
Figure 14 and
Figure 15 are two photos of the experiment environment.
To achieve a good direction, high gain and high angular resolution requirements, we used the developed uniform fan-shaped directional Vivaldi antenna array to receive signals. These four antennas were denoted from Antenna 1 to Antenna 4 in counter-clockwise order. The angle between the symmetry axes of Antenna 1 and Antenna 4 was , and the angle between the symmetry axes of adjacent antennas was . The location of the experimental site was calibrated by infrared range finders, and PD signals were generated by the lighter. We provide an example of the PD signal received from the lighter.
In
Figure 16, the signal contains 50 pulses. We used a high sample rate of 10 GHz to verify our algorithm. We calculated the frequency spectrum of a single pulse using DFT, as shown in
Figure 17. Normally, interference is stronger at low frequencies than at high frequencies. In this paper, we chose the frequency band 1–2 GHz for performance verification. Whether the frequency band 0.5–1 GHz can be used for the algorithms could be studied in subsequent research.
From
Figure 18, we chose three frequency bands, 1.3–1.5 GHz, 1.5–1.7 GHz and 1.7–1.9 GHz, to run the wideband Dir-MUSIC algorithm. We measured the antenna patterns for all four antennas every
in the band 1–2 GHz. Because it is a fan-shaped array, we only measured the antenna pattern every
, from
to
, and the
direction was the direction of the symmetry axis. The measured normalized patterns at 1.4 GHz, 1.6 GHz and 1.8 GHz of all four antennas are shown in
Figure 18 and
Figure 19.
Next, we will process the experimental data with the wideband Dir-MUSIC algorithm, and the mean error and standard deviation of the direction results will be shown in
Table 1. Additionally, it should be noted that the location area is the fan-shaped area, which means that the true azimuth is in
.
In
Table 1, the first column is the PD coordinate (distance and azimuth). At each position, angles in the first row are mean errors (between the average direction result and the true azimuth) and angles in the second row are standard deviation. The mean errors of DirSP are smaller than other methods in most situations. Only at (7.8 m, 30
) are the mean errors of ISSM and CSSM smaller than DirSP. The standard deviation of DirSP is also smaller than other methods in most situations. More importantly, the mean error and standard deviation are both smaller than
. Bad results can be obtained using other methods. In general, the experimental data-processing results are consistent with simulations and DirSP performs better than other methods.
5. Conclusions
In this paper, wideband Dir-MUSIC algorithms to find PD directions are discussed, and a novel wideband Dir-MUSIC algorithm using the strength proportion is proposed. Conclusions are drawn as follows.
When the antenna gain pattern functions among different frequency bins are quite different, the Dir-MUSIC algorithm only using one function will lose accuracy. Meanwhile, DirSP using three functions performs very well.
It is non-effective for CSSM to be directly applied to the Dir-MUSIC algorithm because the steer vectors have different modules, and DirSP is another application of CSSM in this paper.
ISSM may perform better when the signals are generated in an ideal model. However, DirSP performs better when the strength proportion among the frequency bins is unknown. DirSP is widely applicable in most situations.
Wideband Dir-MUSIC algorithms can avoid the influence of interference signals.
Finally, experimental data-processing has verified the effectiveness of DirSP. Compared to Dir-MUSIC, ISSM and CSSM, DirSP performs better in terms of both mean errors and standard deviation. More importantly, DirSP is more stable than other methods, while some bad results are obtained using other methods. Therefore, DirSP is effective for PD location.