Research on Partial Discharge Source Localization Based on an Ultrasonic Array and a Step-by-Step Over-Complete Dictionary

Partial discharge (PD) in electrical equipment is one of the major causes of electrical insulation failures. Fast and accurate positioning of PD sources allows timely elimination of insulation faults. In order to improve the accuracy of PD detection, this paper mainly studies the direction of arrival (DOA) estimation of PD ultrasonic signals based on a step-by-step over-complete dictionary. The simulation results show that the step by step dictionary can improve the operation speed and save signal processing time. Firstly, a step-by-step over-complete dictionary covering all the angles of space is established according to the expression of the steering vector for a matching pursuit direction finding algorithm, which can save computation time. Then, the step-by-step complete dictionary is set up according to the direction vector, and the atomic precision is respectively set to 10◦, 1◦ and 0.1◦. The matching pursuit algorithm is used to carry out the sparse representation of the received data X and select the optimal atom from the step-by-step complete dictionary, and the angle information contained in atoms is DOA of the PD sources. According to the direction finding results, combined with the installation location of the ultrasonic array sensor, the spatial position of a partial discharge source can be obtained using the three platform array location method. Finally, a square ultrasonic array sensor is developed, and an experimental platform for the ultrasonic array detection of partial discharges is set up and used to carry out an experimental study. The results show that the DOA estimation method based on a step-by-step over-complete dictionary can improve the direction finding precision, thereby increasing the subsequent positioning accuracy, and the spatial position estimation error of the PD source obtained under laboratory conditions is about 5 cm, making this a feasible method.


Introduction
The deterioration of insulation in electrical equipment is mainly caused by partial discharge (PD) [1][2][3].Therefore, the timely identification of the partial discharge sources [4][5][6] in electrical equipment is very important for the stability of the whole power system.The principle of the ultrasonic array positioning method [7][8][9] is that a sensor array receives ultrasonic signals generated by PD sources and uses array signal processing techniques to determine the source positions.The ultrasonic array positioning method has advantages of strong anti-interference capability and high accuracy.
In the last century, Schmidt and Roy proposed the multiple signal classification method (MUSIC) [9][10][11][12][13] and the estimation of signal parameters via rotational invariance technique (ESPRIT) algorithm [14,15].The covariance matrix of the received signal is decomposed into two signal and noise subspaces, and then an orthogonal relation is used to estimate the direction of arrival (DOA) of the signal.Although these two methods achieve good results in high Signal to Noise Ratio SNR situations, the performance is more obvious when the number of sampling points is small or the signal to noise ratio is low.Later, a lot of new algorithms like the maximum likelihood (ML) algorithm [16,17] and weighted subspace fitting techniques [18], etc., have been commonly used.Compared to these improved algorithms and subspace-based algorithms, the estimation performance of the improved algorithms is greatly improved, but the computational burden is often too large to be practical.
In addition to the design and implementation of efficient and fast decomposition algorithms, the most important task is to reduce the complexity of the algorithm, while finding the appropriate atomic type and then construct a suitable over-complete dictionary is another key issue.Mallat and Zhang [19,20] pointed out that the over-complete dictionary of atoms is a set of Gauss functions that are scalable, time shifted, and frequency adjusted, known as the Gabor atom dictionary.Because this dictionary has very good time-frequency aggregation, therefore, this kind of dictionary sparse representation clearly shows the time-frequency characteristics based on the internal signal.
In this work, a square ultrasonic array sensor [21,22] is taken as an example.Firstly, the step over-complete dictionary [23] which can effectively reduce the operation time is constructed according to the direction finding need and the array flow pattern of the sensors.Then the matching pursuit (MP) algorithm [24][25][26] is used according to the step-by-step over-complete dictionary to search for the optimal DOA estimation.Moreover, using the spatial direction finding angle, the spatial position of the PD source is determined by using the three array platforms location method.Finally, the feasibility of the theory is verified by simulations and experiments.

Mathematical Model of Array Signal Processing
When the estimation of the DOA of the ultrasonic signals generated by a PD source is conducted based on signal sparse decomposition theory, the first step is to build a mathematical model of the array signals.For the PD ultrasonic array detection system, we suppose that the array sensor is composed of M ultrasonic sensors, there are k PD sources, and the ultrasonic signals on the array sensor are also k, the angle of the ultrasonic signal is ϕ i , i = 1, 2, ..., k.Moreover, the angle of the incident signal ϕ i is composed of two angles, the plane angle α and the pitch angle θ for the plane square array sensor.In the ideal case, considering that each element is in the same direction and independent from each other, the received signal of the m-th element at any time t is: where the g mi is the gain of the ith signal for the m-th array element, which can be normalized to 1. s i is the i-th PD source.n m (t) is noise for the m-th array element.mi is time delay of the i-th spatial signal incident on the m-th array element, which contains the spatial parameters of the signal source DOA.Because s i (t − τ) ≈ s i (t)e −jωτ , i = 1, 2, . . .k, the data received by the array sensor is arranged in a vector form, and the data can be obtained from the following expression: where ω = 2π f = 2π c λ , c is speed of sound, and λ is the wavelength.Equation ( 2) can be characterized as follows: A is a space array of M × k dimensional vector matrix and: where a i (ϕ i ) is oriented vector, ϕ i is the i-th space signal DOA (including the plane angle α and the pitch angle θ), and: is the oriented vector determined by the array mainfold.τ denotes the delay between the array elements.For a different array structure, τ will have a different expression.

Square Ultrasonic Array Structure
In a square array ultrasonic sensor composed of m 2 array elements, the array element is arranged in the form shown in Figure 1 (taking 9-yuan as example).Each edge has m elements, and according to the number of the icon of each array element number (the number ab means that the array element is located in the first a row and b column), the spacing between the array elements is d.where , c is speed of sound, and  is the wavelength.Equation ( 2) can be characterized as follows: A is a space array of M k  dimensional vector matrix and: where ( )   i i a  is oriented vector, i  is the i-th space signal DOA (including the plane angle α and the pitch angle θ), and: is the oriented vector determined by the array mainfold.τ denotes the delay between the array elements.For a different array structure, τ will have a different expression.

Square Ultrasonic Array Structure
In a square array ultrasonic sensor composed of m 2 array elements, the array element is arranged in the form shown in Figure 1 (taking 9-yuan as example).Each edge has m elements, and according to the number of the icon of each array element number (the number ab means that the array element is located in the first a row and b column), the spacing between the array elements is d.The ultrasonic propagation speed is set to c, the NO11 array element is taken as the reference element, where the position of the element is taken as the origin of the coordinates, and a threedimensional rectangular coordinate system is established.Suppose that the incident signal direction is (α, θ), azimuth angle α is the angle between the reference point to the source of the connection in the X-Y plane of projection and x axis, the pitch angle θ is the angle between the reference point to the source of the connection with the z axis.Then, the geometric relations can be known, and according to the reference array, the time delay of the ultrasonic signal arriving the ab NO array element is: where da is the distance projection between the ab NO array element and the reference array on the x axis.And db is the distance projection on the y axis.Equation ( 5) can be expressed as follows: The ultrasonic propagation speed is set to c, the NO 11 array element is taken as the reference element, where the position of the element is taken as the origin of the coordinates, and a three-dimensional rectangular coordinate system is established.Suppose that the incident signal direction is (α, θ), azimuth angle α is the angle between the reference point to the source of the connection in the X-Y plane of projection and x axis, the pitch angle θ is the angle between the reference point to the source of the connection with the z axis.Then, the geometric relations can be known, and according to the reference array, the time delay of the ultrasonic signal arriving the NO ab array element is: where d a is the distance projection between the NO ab array element and the reference array on the x axis.And d b is the distance projection on the y axis.Equation ( 5) can be expressed as follows: where f is the center frequency of the signal.

Ultrasonic Array Direction Finding and Location Based on
Step-by-Step Over-Complete Dictionary

DOA Estimation Based on the Matching Pursuit Algorithm
For the ultrasonic signal DOA estimation of k PD sources, the incidence angle of the signal can be searched in the N angle vector that has been set.Usually the number of sources is much smaller than the number of angles to be searched, that is k ≤ N.An angle search vector is constructed in the angle space, which has k components which are respectively equal to ϕ 1 , ϕ 2 , . . ., ϕ k .Constructing an angle search matrix is also an over-complete atomic dictionary: When the matching pursuit algorithm is used to decompose the received array signal, the signal in each direction has a maximum projection value on the corresponding atom.According to the mathematical projection theory, the signal and the maximum value of the inner product of the atom can be used to find the largest projection, which means to find the best matching atom vector.Suppose the sampling point of the signal sampling is T, then output signal X is a vector matrix of M × T dimensions.Atom α ν β(n) in the A s (β), n = 1, 2, ..., N, is the column vector for the M × 1 dimension.The collected array signal X forms the inner product with each atomic vector α ν β(n) in the matching pursuit algorithm.The most matched atom is selected by the size of the absolute value of the inner product, α ν (β γ0 ), γ0 ∈ {0, 1, ..., N}.The most matched atomic vectors should be consistent with: Then the received signal can be expressed in the projection of α ν (β γ0 ) and the residual signal: where P a v (β γ0 ) X is the projection part of the signal on the most matched atom.According to the definition of matrix projection, the projection part can be obtained from the following formula: Repeat the above steps for the optimal matched signal residual, the signal residual is small enough after K times and meets the error requirement, then the sparse decomposition result of the array signal X(t) is: After the completion of the reception signal decomposition, k optimal atoms are obtained.The angular information contained in the atom is composed of the DOA estimation matrix, which is

The Establishment of Step over Complete Atomic Dictionary
According to the need of realizing the signal DOA estimation, a complete vector atomic dictionary is constructed based on the vector matrix As(β).For the square array sensor, β is the vector constituted by the incident azimuth and pitch angle, which is β = [α, θ], Over a complete vector atomic dictionary is As(α, θ), the range of α and θ are all 0°-180°, so the direction finding accuracy acquires a uniform value, and the vector atom is determined by the values of α and θ.According to Equation ( 6), the following models are established for each vector atom in the atomic dictionary of the square array sensor: The process of the matching algorithm to find the optimal atom is the one by one inner product of atoms between the amount of decomposition and the over-complete dictionary, so the number of atoms in a complete dictionary directly determines the amount of computation needed.For the plane shaped ultrasonic array sensor, the DOA estimation can simultaneously measure the azimuth and pitch angle, which range from 0°-180°.If the accuracy in a complete dictionary, where the accuracy is the step size of two angles in the adjacent vector atoms is set to 0.1, then the number of atoms will reach 1800 × 1800 directly according to the expected accuracy needed to establish a complete dictionary which will make the subsequent search require a huge amount of computing, so it is necessary to set up a step over-complete dictionary.In the establishment of a step over-complete dictionary, the precision of the first established over-complete dictionary is 10°, the search ranges are

The Establishment of Step over Complete Atomic Dictionary
According to the need of realizing the signal DOA estimation, a complete vector atomic dictionary is constructed based on the vector matrix A s (β).For the square array sensor, β is the vector constituted by the incident azimuth and pitch angle, which is β = [α, θ], Over a complete vector atomic dictionary is A s (α, θ), the range of α and θ are all 0 • -180 • , so the direction finding accuracy acquires a uniform value, and the vector atom is determined by the values of α and θ.According to Equation ( 6), the following models are established for each vector atom in the atomic dictionary of the square array sensor: The process of the matching algorithm to find the optimal atom is the one by one inner product of atoms between the amount of decomposition and the over-complete dictionary, so the number of atoms in a complete dictionary directly determines the amount of computation needed.For the plane shaped ultrasonic array sensor, the DOA estimation can simultaneously measure the azimuth and pitch angle, which range from 0 • -180 • .If the accuracy in a complete dictionary, where the accuracy is the step size of two angles in the adjacent vector atoms is set to 0.1, then the number of atoms will reach 1800 × 1800 directly according to the expected accuracy needed to establish a complete dictionary which will make the subsequent search require a huge amount of computing, so it is necessary to set up a step over-complete dictionary.In the establishment of a step over-complete dictionary, the precision of the first established over-complete dictionary is 10 • , the search ranges are all 0 • ~180 • , the first step atom number is 18 × 18, and the first rough time estimation is obtained.Then the precision of the second established over-complete dictionary is 1 • , the search range is [η 1 − 10 • , η 1 + 10 • ], the second step atom number is 20 × 20, and the second more accurate estimate is obtained, where the search angle is η 2 .Similar to the above steps, the third over-complete dictionary is established, the degree is set to 0.1 • , the search range is [η 2 − 1 • , η 2 + 1 • ], the third step atom number is 20 × 20, and a high accuracy DOA estimation η 3 can be obtained with three precision searches, but the amount of computation required compared to a directly set precision degree like 0.1 • is much smaller.

Three Array Cross Positioning Principle
The location of the PD source can't be determined after obtaining the DOA using the MP method.In this paper, the direction finding results are used, and we combine the three array cross location method [10] to obtain the spatial position of the local source.The three array positioning principle is shown in Figure 3.
Then the precision of the second established over-complete dictionary is 1°, the search range is [η1 − 10°, η1 + 10°], the second step atom number is 20 × 20, and the second more accurate estimate is obtained, where the search angle is η2.Similar to the above steps, the third over-complete dictionary is established, the degree is set to 0.1°, the search range is [η2 − 1°, η2 + 1°], the third step atom number is 20 × 20, and a high accuracy DOA estimation η3 can be obtained with three precision searches, but the amount of computation required compared to a directly set precision degree like 0.1° is much smaller.

Three Array Cross Positioning Principle
The location of the PD source can't be determined after obtaining the DOA using the MP method.In this paper, the direction finding results are used, and we combine the three array cross location method [10] to obtain the spatial position of the local source.The three array positioning principle is shown in Figure 3.The space positions of the three groups of array sensors are set to A1 (x1, y1, z1), A2 (x2, y2, z2), A3 (x3, y3, z3), the direction finding and the position of the array sensor are obtained by using DOA estimation, then the direction finding line equation St can be obtained.The spatial coordinate of the signal source is set to S (x, y, z), and the three different direction finding lines should intersect in S (x, y, z) under ideal conditions, but there are many kinds of measurement errors, these lines are on different surfaces.Therefore, the sum function of the sum of the three straight lines can be constructed: The chaotic monkey [27] algorithm is used to search optimize the space, which makes the minimum sum of the distance.The point where the d value is the smallest, is regarded as the spatial location of the source.

Experimental Study
To start with, the PD signals are collected by using the PD detection experimental platform in the lab.Then the matching pursuit algorithm is prepared with the help of the MATLAB (2012 v7.0, Math Works, Natick, MA, USA) programming software to capture the array signal DOA estimation.Finally, after obtaining the space angle information of the PD source, combined with the sensor array location information, the spatial position of the PD source is determined by the three array cross positioning method.The space positions of the three groups of array sensors are set to A1 (x 1 , y 1 , z 1 ), A2 (x 2 , y 2 , z 2 ), A3 (x 3 , y 3 , z 3 ), the direction finding and the position of the array sensor are obtained by using DOA estimation, then the direction finding line equation St can be obtained.The spatial coordinate of the signal source is set to S (x, y, z), and the three different direction finding lines should intersect in S (x, y, z) under ideal conditions, but there are many kinds of measurement errors, these lines are on different surfaces.Therefore, the sum function of the sum of the three straight lines can be constructed: The chaotic monkey [27] algorithm is used to search optimize the space, which makes the minimum sum of the distance.The point where the d value is the smallest, is regarded as the spatial location of the source.

Experimental Study
To start with, the PD signals are collected by using the PD detection experimental platform in the lab.Then the matching pursuit algorithm is prepared with the help of the MATLAB (2012 v7.0, Math Works, Natick, MA, USA) programming software to capture the array signal DOA estimation.
Finally, after obtaining the space angle information of the PD source, combined with the sensor array location information, the spatial position of the PD source is determined by the three array cross positioning method.

Experimental System
The main components of the PD ultrasonic array detection hardware system include: simulation of discharge system, signal receiving and transmission system, synchronous multi-channel data acquisition and processing system, and the computer.This system architecture is shown in Figure 4.

Experimental System
The main components of the PD ultrasonic array detection hardware system include: simulation of discharge system, signal receiving and transmission system, synchronous multi-channel data acquisition and processing system, and the computer.This system architecture is shown in Figure 4.The model distribution transformer is a rectangular box body that is welded by steel plates.Three sides are composed of steel plates, and the other side is a glass plate, which makes it easy to observe the discharge source location.The body size is a length of 150 cm, width of 100 cm, and the height is 120 cm.A discharge tube (EPSON, Shenzhen, China) is used as a discharge device, and the critical discharge voltage is 230 V.The frequency of ultrasonic signals registered by the sensor array is in the range of 50 kHz-280 kHz, the center frequency is 150 kHz, and the equivalent velocity is 1500 m/s, while the wavelength is about 10 mm.
During the experimental process, the PD ultrasonic signals are acquired with a square array sensor [22,28], which is shown in Figures 5 and 6.In Figure 6 there are magnets on the top and bottom sides of the assembly.The ultrasonic array sensor assembly is tightly attached on the outer wall of the electrical equipment model.Moreover, the different order square array structures are changed by changing the positions in the assembly and the array elements, as shown in Figure 6.During the experiments, the environmental noise ratio was about 10 dB.The PD ultrasonic signal is difficult to collect, and a variety of interferences could affect the experimental results.Each of the PD source locations needs to collect more data.The model distribution transformer is a rectangular box body that is welded by steel plates.Three sides are composed of steel plates, and the other side is a glass plate, which makes it easy to observe the discharge source location.The body size is a length of 150 cm, width of 100 cm, and the height is 120 cm.A discharge tube (EPSON, Shenzhen, China) is used as a discharge device, and the critical discharge voltage is 230 V.The frequency of ultrasonic signals registered by the sensor array is in the range of 50 kHz-280 kHz, the center frequency is 150 kHz, and the equivalent velocity is 1500 m/s, while the wavelength is about 10 mm.
During the experimental process, the PD ultrasonic signals are acquired with a square array sensor [22,28], which is shown in Figures 5 and 6.In Figure 6 there are magnets on the top and bottom sides of the assembly.The ultrasonic array sensor assembly is tightly attached on the outer wall of the electrical equipment model.Moreover, the different order square array structures are changed by changing the positions in the assembly and the array elements, as shown in Figure 6.During the experiments, the environmental noise ratio was about 10 dB.The PD ultrasonic signal is difficult to collect, and a variety of interferences could affect the experimental results.Each of the PD source locations needs to collect more data.

Experimental System
The main components of the PD ultrasonic array detection hardware system include: simulation of discharge system, signal receiving and transmission system, synchronous multi-channel data acquisition and processing system, and the computer.This system architecture is shown in Figure 4.The model distribution transformer is a rectangular box body that is welded by steel plates.Three sides are composed of steel plates, and the other side is a glass plate, which makes it easy to observe the discharge source location.The body size is a length of 150 cm, width of 100 cm, and the height is 120 cm.A discharge tube (EPSON, Shenzhen, China) is used as a discharge device, and the critical discharge voltage is 230 V.The frequency of ultrasonic signals registered by the sensor array is in the range of 50 kHz-280 kHz, the center frequency is 150 kHz, and the equivalent velocity is 1500 m/s, while the wavelength is about 10 mm.
During the experimental process, the PD ultrasonic signals are acquired with a square array sensor [22,28], which is shown in Figures 5 and 6.In Figure 6 there are magnets on the top and bottom sides of the assembly.The ultrasonic array sensor assembly is tightly attached on the outer wall of the electrical equipment model.Moreover, the different order square array structures are changed by changing the positions in the assembly and the array elements, as shown in Figure 6.During the experiments, the environmental noise ratio was about 10 dB.The PD ultrasonic signal is difficult to collect, and a variety of interferences could affect the experimental results.Each of the PD source locations needs to collect more data.

Direction Finding and Location
During the experiment, the space position of the PD source is set up in advance.Then a plurality of preset fixed sensor locations are selected the outer wall of the tank.In order to ensure the experiments are performed under the same conditions and facilitate subsequent data processing, during the experiment, the PD source location is fixed, and only the position of the ultrasonic sensor array is changed, and thus we change the relative space position.In the simulation experiments, the PD source location coordinates are (65, 15, 60) cm.The PD ultrasonic array signal waveform in oil is collected by a square ultrasonic sensor array, as shown in Figure 7 The three array sensors are fixed at (45, 0, 20) cm, (80, 0, 20) cm, (0, 30, 40) cm, respectively.The three DOA theoretical values are: (36.9°, 32.0°), (135.0°,27.9°), (167.0°,73.3°).According to the principle of the matching pursuit (MP) algorithm for DOA estimation, the process of the array signal in the search for the optimal atom is to form the inner product with each atom, respectively.The absolute value of the inner product with the optimal atom is maximum.In order to represent this process visually, the absolute value of the inner product distribution from the perspective of space can be calculated.
In step-by-step over-complete dictionary, the atomic precisions are respectively 10°, 1° and 0.1°.Each finding process is divided into three steps.The ultrasonic sensor array position 1 is used as an example.The PD ultrasonic signal data is collected by an ultrasonic array sensor in position 1 and is used the MP algorithm for DOA estimation.The direction finding result is (37.5°,31.0°).The absolute values of the inner product in the space distribution are shown in Figures 8-10.

Direction Finding and Location
During the experiment, the space position of the PD source is set up in advance.Then a plurality of preset fixed sensor locations are selected the outer wall of the tank.In order to ensure the experiments are performed under the same conditions and facilitate subsequent data processing, during the experiment, the PD source location is fixed, and only the position of the ultrasonic sensor array is changed, and thus we change the relative space position.In the simulation experiments, the PD source location coordinates are (65, 15, 60) cm.The PD ultrasonic array signal waveform in oil is collected by a square ultrasonic sensor array, as shown in Figure 7.

Direction Finding and Location
During the experiment, the space position of the PD source is set up in advance.Then a plurality of preset fixed sensor locations are selected the outer wall of the tank.In order to ensure the experiments are performed under the same conditions and facilitate subsequent data processing, during the experiment, the PD source location is fixed, and only the position of the ultrasonic sensor array is changed, and thus we change the relative space position.In the simulation experiments, the PD source location coordinates are (65, 15, 60) cm.The PD ultrasonic array signal waveform in oil is collected by a square ultrasonic sensor array, as shown in Figure 7 The three array sensors are fixed at (45, 0, 20) cm, (80, 0, 20) cm, (0, 30, 40) cm, respectively.The three DOA theoretical values are: (36.9°, 32.0°), (135.0°,27.9°), (167.0°,73.3°).According to the principle of the matching pursuit (MP) algorithm for DOA estimation, the process of the array signal in the search for the optimal atom is to form the inner product with each atom, respectively.The absolute value of the inner product with the optimal atom is maximum.In order to represent this process visually, the absolute value of the inner product distribution from the perspective of space can be calculated.
In step-by-step over-complete dictionary, the atomic precisions are respectively 10°, 1° and 0.1°.Each finding process is divided into three steps.The ultrasonic sensor array position 1 is used as an example.The PD ultrasonic signal data is collected by an ultrasonic array sensor in position 1 and is used the MP algorithm for DOA estimation.The direction finding result is (37.5°,31.0°).The absolute values of the inner product in the space distribution are shown in Figures 8-10.The three array sensors are fixed at (45, 0, 20) cm, (80, 0, 20) cm, (0, 30, 40) cm, respectively.The three DOA theoretical values are: (36.9 • , 32.0 • ), (135.0 • , 27.9 • ), (167.0 • , 73.3 • ).According to the principle of the matching pursuit (MP) algorithm for DOA estimation, the process of the array signal in the search for the optimal atom is to form the inner product with each atom, respectively.The absolute value of the inner product with the optimal atom is maximum.In order to represent this process visually, the absolute value of the inner product distribution from the perspective of space can be calculated.
In step-by-step over-complete dictionary, the atomic precisions are respectively 10 • , 1 • and 0.1 • .Each finding process is divided into three steps.The ultrasonic sensor array position 1 is used as an example.The PD ultrasonic signal data is collected by an ultrasonic array sensor in position 1 and is used the MP algorithm for DOA estimation.The direction finding result is (37.5 • , 31.0 • ).The absolute values of the inner product in the space distribution are shown in Figures 8-10.In the same way, a step and step over-complete dictionary is used to measure the direction of the remaining two positions.The measurement results for array sensor position 2 and position 3 are as follows: (134.0°,26.7°), (167.6°,74.6°).Then, using the three DOA positions that have been measured, the spatial position of the PD source is obtained according to the three array cross positioning principle.The results are (67.7,17.6, 63.3) cm, and the error is 5.0 cm.The location diagram is shown in Figure 11.In the same way, a step and step over-complete dictionary is used to measure the direction of the remaining two positions.The measurement results for array sensor position 2 and position 3 are as follows: (134.0°,26.7°), (167.6°,74.6°).Then, using the three DOA positions that have been measured, the spatial position of the PD source is obtained according to the three array cross positioning principle.The results are (67.7,17.6, 63.3) cm, and the error is 5.0 cm.The location diagram is shown in Figure 11.In the same way, a step and step over-complete dictionary is used to measure the direction of the remaining two positions.The measurement results for array sensor position 2 and position 3 are as follows: (134.0°,26.7°), (167.6°,74.6°).Then, using the three DOA positions that have been measured, the spatial position of the PD source is obtained according to the three array cross positioning principle.The results are (67.7,17.6, 63.3) cm, and the error is 5.0 cm.The location diagram is shown in Figure 11.In the same way, a step and step over-complete dictionary is used to measure the direction of the remaining two positions.The measurement results for array sensor position 2 and position 3 are as

Conclusions
The DOA estimation method based on sparse representation of received data is studied, and the step-by-step over-complete dictionary is established.Firstly, the steering vector of the full array ultrasonic array sensor is derived, which is used as the model of the atoms in the dictionary.Then, the direction finding method based on the direct sparse representation of the received data is applied to the detection of the local discharge ultrasonic array.Finally, an experimental study shows that based on the step-by-step over-complete dictionary of transformer partial discharge detection, the PD source localization error is 5 cm, which meets the actual needs of the project.
realizes the signal DOA estimation.The ultrasonic array signal flow chart of DOA estimation based on matching pursuit algorithm is shown in Figure2.

Figure 2 .
Figure 2. The flow chart based on matching pursuit direction of arrival (DOA) estimation.

Figure 2 .
Figure 2. The flow chart based on matching pursuit direction of arrival (DOA) estimation.

Figure 3 .
Figure 3. Schematic diagram of the three array cross location.

Figure 3 .
Figure 3. Schematic diagram of the three array cross location.

Figure 4 .
Figure 4.The partial discharge (PD) detection experiment system structure diagram.

Figure 5 .
Figure 5. Model and physical aspect of a single array element.(a) The form composition of element of the PD ultrasonic array sensors; and (b) the picture of element of the PD ultrasonic array sensors.

Figure 4 .
Figure 4.The partial discharge (PD) detection experiment system structure diagram.

Figure 4 .
Figure 4.The partial discharge (PD) detection experiment system structure diagram.

Figure 5 .
Figure 5. Model and physical aspect of a single array element.(a) The form composition of element of the PD ultrasonic array sensors; and (b) the picture of element of the PD ultrasonic array sensors.Figure 5. Model and physical aspect of a single array element.(a) The form composition of element of the PD ultrasonic array sensors; and (b) the picture of element of the PD ultrasonic array sensors.

Figure 5 .
Figure 5. Model and physical aspect of a single array element.(a) The form composition of element of the PD ultrasonic array sensors; and (b) the picture of element of the PD ultrasonic array sensors.Figure 5. Model and physical aspect of a single array element.(a) The form composition of element of the PD ultrasonic array sensors; and (b) the picture of element of the PD ultrasonic array sensors.

Figure 8 .Figure 9 .Figure 10 .
Figure 8.The accuracy of array direction finding is 10° in position 1.(a) Three-dimensional map; and (b) contour map (the color is proportional to surface height, and dark red represents the highest point).

Figure 8 .Figure 8 .Figure 9 .Figure 10 .
Figure 8.The accuracy of array direction finding is 10 • in position 1.(a) Three-dimensional map; and (b) contour map (the color is proportional to surface height, and dark red represents the highest point).

Figure 9 .Figure 8 .Figure 9 .Figure 10 .
Figure 9.The accuracy of array direction finding is 1 • in position 1.(a) Three-dimensional ma; and (b) contour map (the color is proportional to surface height, and dark red represents the highest point).

Figure 10 .
Figure 10.The accuracy of array direction finding is 0.1 • in position 1.(a) Three-dimensional map; and (b) contour map (the color is proportional to surface height, and dark red represents the highest point).