Next Article in Journal
Manipulation of Unknown Objects to Improve the Grasp Quality Using Tactile Information
Next Article in Special Issue
Integrated Method for Personal Thermal Comfort Assessment and Optimization through Users’ Feedback, IoT and Machine Learning: A Case Study
Previous Article in Journal
Demonstration of a Speckle Based Sensing with Pulse-Doppler Radar for Vibration Detection
Previous Article in Special Issue
Detection of Anomalous Noise Events on Low-Capacity Acoustic Nodes for Dynamic Road Traffic Noise Mapping within an Hybrid WASN
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radio-Frequency Localization of Multiple Partial Discharges Sources with Two Receivers

by
Guillermo Robles
*,
José Manuel Fresno
and
Juan Manuel Martínez-Tarifa
Department of Electrical Engineering, Universidad Carlos III de Madrid, Avda, Universidad, 30, Leganés, 28911 Madrid, Spain
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(5), 1410; https://doi.org/10.3390/s18051410
Submission received: 3 April 2018 / Revised: 27 April 2018 / Accepted: 28 April 2018 / Published: 3 May 2018

Abstract

:
Spatial localization of emitting sources is especially interesting in different fields of application. The focus of an earthquake, the determination of cracks in solid structures, or the position of bones inside a body are some examples of the use of multilateration techniques applied to acoustic and vibratory signals. Radar, GPS and wireless sensors networks location are based on radiofrequency emissions and the techniques are the same as in the case of acoustic emissions. This paper is focused on the determination of the position of sources of partial discharges in electrical insulation for maintenance based on the condition of the electrical equipment. The use of this phenomenon is a mere example of the capabilities of the proposed method but it is very representative because the emission can be electromagnetic in the VHF and UHF ranges or acoustic. This paper presents a method to locate more than one source in space with only two receivers, one of them in a fixed position and the other describing a circumference around the first one. The signals arriving from the different sources to the antennas are first separated using a classification technique based on their spectral components. Then, the individualized time differences of arrival (TDOA) from the sources collected at different angles describe a function, angle versus TDOA, that has all the geometric information needed to locate the source. The paper will show how to derive these functions for any source analytically with the position of the source as unknown parameters. Then, it will be demonstrated that it is possible to fit the curve with experimental measurements of the TDOA to obtain the parameters of the position of each source. Finally, the technique is extended to the localization of the emitter in three dimensions.

1. Introduction

The localization of emitting sources is an interesting research topic which covers several fields, from earthquake detection, [1], to RADAR and global positioning systems applied to the remote control of autonomous vehicles [2] or indoor positioning systems [3]. All these fields can use multilateration techniques applied to radio-frequency, acoustic and vibration signals which are acquired with several fixed sensors. Generally, when a pulse is radiated from a source, the wave will arrive to two receivers at different times. One of the advantages of measuring these time differences of arrival or TDOA is that a common clock is not required as in other localization techniques based on the time of arrival (TOA) of the pulse to the receiver. With two sensors, all the possible points in space that would give the same TDOA describe a hyperboloid [4]. Using two more independent receivers and calculating the intersection of the hyperboloids using all possible pairs of antennas would give the position of the source in three dimensions. Therefore, the spatial localization of emitters using multilateration techniques can be solved with at least four receivers.
In this work, the pulsed sources under study will come from a phenomenon known as partial discharges (PD) which are small-energy ionizations taking place within insulation materials revealing accelerated ageing [5]. Nevertheless, the results obtained here are not restricted to this type of pulsed source, since the same strategy could be applied to other types of signals whose source needed to be localized. It is well known that the detection of PD can be used to prevent premature failures of power cables and electric machines [6,7,8,9]. The interest in this technique is increasing since it can give indications on the insulation status with the equipment in service. In the particular case of open-air substations, the radio-frequency (RF) radiation emitted by PD can be used to detect and locate the imperfections, giving valuable information for condition monitoring [10,11,12].
There are other methods for localization based on the received signal strength (RSS) [13,14,15,16] or the angle of arrival (AOA) [17,18,19]. The use of the received signal strength indicator (RSSI) requires a previous mapping of the area to survey using a known source of discharges which sometimes would require putting the substation out of service. The sparker is placed in all points of a predefined grid and the power of the signal is recorded with a set of receivers. When a PD occurs in that mapped area, the signal acquired by the receivers is compared with the grid to set the approximate position of the source. The accuracy on the localization of the source depends on the density of data recorded in the grid and the previous elimination of all partial discharges that occurred during the characterization of the area. Concerning the calculation of the AOA, the usual approach is to use four receivers in the VHF and UHF bands or ultrasonic range forming and array of 2 × 2 elements. This array is virtually expanded to increase the resolution and then, the direction of the emission of the PD is computed with a multiple signal classification (MUSIC) algorithm [19].
Typically, the localization of PD sources using RF signals with multilateration techniques has been done through a system of four omnidirectional antennas [20]. The position of the receivers and TDOA of the pulses arriving to these sensors form a system of non-linear equations that can be solved to locate the emitter [21]. However, in the RF range the sampling frequencies ought to be high to have sufficient spatial resolution so the acquisition systems have a fairly high cost. It seems evident that a reduction in the number of receivers could help in the reduction of the cost of the instrumentation, mainly in sensors and high-speed digitizing channels. In this paper, the authors face this problem of moving at least one of the receivers with the aim of reducing the number of sensors using the known position of their relative movement to obtain the needed TDOA for the localization of the sources. Therefore, the authors propose a system for spatial localization of PD sources by means of two sensors mechanically linked by a rod. In order to get the necessary number of different TDOA for the localization, one of the sensors is rotated fixing the other one as a geometric reference. The mathematical function relating the relative distance between the two sensors, the rotation angle of the moveable one, and the TDOA will be derived in the paper showing good agreement with the experimental results. This approach was firstly applied in [22] obtaining good results for one source in the same plane as the antennas. In this paper, the method has been extended to two simultaneous PD sources and the localization of the sources in 3D. The casuistic when two or more sources are active is already complicated with four antennas, so reducing the number of antennas poses an additional challenge. This is tackled in the manuscript separating the origin of the signals based on their spectral signature using power ratio (PR) maps, a classification method based on the differences of power in certain bands of frequency [12,23].
The paper is structured in the following sections. First, the mathematical and geometrical principles of the method are explained in Section 2. Then, the problem of localizing more than one source is faced in Section 3 using the classification maps based on the spectral signature of the signals introduced in Section 4. Section 5 describes the measuring setup, the instrumentation and the experimental measurements. The algorithm and the method are extended to the localization of more than one source in three dimensions in Section 6. Finally the conclusions and insights of the work done in the paper are drawn in the last section.

2. Localization Basis

The proposed method is based on the information of the arrival of pulsed signals to two sensors, in this case, two simple omnidirectional antennas [24], as a function of their relative position with respect to the source. The TOA from the source to the antenna i, t i , is not known since the position of the source has not yet been determined. Therefore, the localization relies on the time differences of arrival, t 12 = t 1 - t 2 of the emission to the antennas as in all multilateration techniques. The novelty of this work is presented henceforth. The distance between antennas is fixed, attaching them to the ends of a shaft or rod. Then, both antennas can be moved jointly when the rig is rotated around the vertical axis passing through the middle point between the antennas or the antenna 1 (▲) can be fixed at the origin of a Cartesian coordinate system while the antenna 2 (▼) rotates around the antenna 1 as shown in Figure 1 where the source (⬤) is placed at ( - 1 , 2 ) m. Though both possibilities are valid, the latter configuration is preferable to moving both antennas simultaneously since one of the receivers is always placed in the same place. This will be the method followed in the rest of the paper. Moreover, this stationary antenna will give the needed reference for the spectral characterization of the pulses explained in Section 4. Additionally, there are two possible methods to determine the position of the source. One of them is based on the intersection of two lines perpendicular to the rod when t 12 = 0 and the other is based on the determination of the dependence of the rotation angle with the TDOA and fitting the curve to the experimental data. These two procedures are explained in the next sections.

2.1. Intersection of Lines

There are four characteristic positions of the antennas depending on certain values of the TDOA: two zero crossings, and maximum and minimum values. The maximum corresponds to the case when the rod is aligned with the direction of the source, therefore, antenna 1, antenna 2 and source are in the same line. Figure 1 shows this particular case when the antenna 2 is placed at 116 . 6 , then t 1 > t 2 and t 12 > 0 . There is also a minimum at the opposite angle, - 63 . 4 , because there are two positions where the three elements are aligned. In this case, the relative locations of antenna 1 and antenna 2 are swapped and the TDOA is negative, t 2 > t 1 and t 12 < 0 . The position of the antennas also crosses zero twice and this happens when the rod is perpendicular to the direction of the source because, in that case, both times of arrival to the antennas are the same, t 1 = t 2 and t 12 = 0 , black radii in Figure 1. These cases are more helpful for this technique because plotting the perpendicular lines at the angular positions defined by the two zero crossings and calculating the intersection of the two lines will determine and solve the position of the source. Figure 1 shows these locations at 34 . 9 and - 161 . 8 and the lines perpendicular to the rod that intersect at the source.
In summary, the procedure to determine the position of the source is straightforward. First, the TDOA at different angles in one turn are calculated analyzing the onset of the pulses of the incoming signals to both antennas using any of the algorithms presented in [25] which, in this work, is based on the Hinkley criterion. Then, the angles where the TDOA are zero, ζ 1 and ζ 2 , are calculated interpolating the closest angles above and below zero TDOA because the rotation is done in discrete angles. The next step traces the line perpendicular to the rod at those angles passing through its middle point. Finally, the position of the source is calculated with the intersection of those lines.

2.2. Analytical Expression of the Periodic Function

The analytical expression in two dimensions can be derived from the setup in Figure 2 where d i is the distance from the antenna i to the source. The antenna 1 is located at the origin of coordinates and the positions of the source and the antenna 2 are defined by their radius and azimuth in cylindrical coordinates ( d 1 , α ) and ( r , θ ) , respectively. The TDOA in meters can be calculated subtracting the distances, d = d 1 - d 2 , which is better done in Cartesian coordinates:
( x s , y s ) = ( d 1 cos α , d 1 sin α ) ( x 2 , y 2 ) = ( r cos θ , r sin θ ) d 2 = ( d 1 cos α - r cos θ , d 1 sin α - r sin θ )
where ( x s , y s ) are the coordinates of the source, ( x 2 , y 2 ) the position of the antenna 2 and r the constant distance between antennas. Therefore, d is calculated as:
d = d 1 - ( d 1 cos α - r cos θ ) 2 + ( d 1 sin α - r sin θ ) 2
or, rearranging the terms, as:
d = d 1 - r 2 + d 1 2 - 2 r d 1 cos ( α - θ )
Equation (3) depends on the position of the source, determined by d 1 and α ; the distance between antennas, r and the rotation angle, θ . The only unknown parameters are d 1 and α , so fitting the median of the experimental data for every angle to the analytical expression with a non-linear least squares algorithm will determine these parameters and, therefore, the position of the source.
Then, the process would consist of measuring the TDOA, d, at different angles completing a total turn to obtain the theoretical function plotted in Figure 3. The units of the TDOA are usually seconds but they can also be represented in meters multiplying time by the speed of light c, or samples, dividing time by the sampling time T s . Concerning the choice of samples for the plot in the figure, one should bear in mind that the sampling frequency and the length of the rod are related and affect the sensitivity, as will be explained in Section 2.3.
There is an alternative formulation of Equation (3) represented in Cartesian coordinates that would give directly the coordinates of the source in this system:
d = x s 2 - y s 2 - ( x s - x 2 ) 2 - ( y s - y 2 ) 2

2.3. Choosing the Separation Distance

The separation between antennas is a critical parameter that will determine the accuracy of the calculated position of the source. The maximum time difference of arrival d m a x , which is also the maximum and minimum of the analytic expression in Equation (3), can be derived considering that cos ( α - θ ) = 1 , or α = θ . In that case, the square root would be minimum and the maximum TDOA in meters would be:
d m a x = d 1 - r 2 + d 1 2 - 2 r d 1 = d 1 - ( d 1 - r ) 2 = r
This result is expected since when the antennas are aligned ( α = θ ) and in the same plane as the source, the relative distance between the antennas and the source would be maximum and equal to the distance between them, r. This parameter is directly related to the sampling frequency and the sensitivity of the rig. In effect, since the sampling time determines the minimum measurable distance, the minimum step of d would be one sample or c T s meters. There is no point in selecting a small step in the rotation angle when the sampling frequency is low because the resolution in the vertical axis of Figure 3 would also be low. We can overcome this issue setting the length of the rod considering that we need to have a number of discrete steps of TDOA in the vertical axis in Figure 3 such as we can reach a sufficiently high resolution in the fitting of the curve. It has been considered that having 10 steps or samples up to the maxima of the curve is adequate, therefore, d m a x > 10 c T s or r > 10 c T s . This figure of merit together with the discrete steps in the rotation angle would ensure that the zero crossings can be determined with good accuracy.

3. Separation and Localization Strategy for Two Sources

Both possibilities of rotating the antennas commented on in Section 2 would give similar results for the position of the source, however, the difficulty arises when there is more than one source in different positions. Since the number of antennas has been reduced to two, there is only one time difference of arrival of the pulse to the receivers. If the measured signals come from two or more different sources, there is no easy way to know which TDOA would correspond to which sources. This problem has already been studied in [10] where the authors study the probability that a set of TDOA come from a certain source, however, we propose another approach based on the analysis of the shape of the pulse which simplifies the localization of multiple sources. It is well known that the shape of the signal received by the antenna depends on the radio-frequency pulse itself, the path from the source to the antenna, the radiometric characteristics of the antenna and interferences. In industrial environments, there can be metallic structures that distort the emitted pulse and the received signal would have a voltage increase due to the front wave and oscillations due to reflections and the antenna response. This poses additional problems in the determination of the pulse onset that have already been commented and solved in [21,26,27]. Since the fingerprint of the signal arriving to the antennas is highly dependent on the relative position of the source and the sensors and the environmental characteristics, the differences in the power spectral density of the signals can be used to determine their source of origin. If one of the antennas is stationary it would always receive signals with fingerprints associated to the sources resulting in a pattern database that will help in the association of signal and source. On the other hand, the mobile antenna would have the information of the time difference of arrival because it receives signals from the sources with a delay that is related to the path followed by the emission. Then, every acquisition has pairs of signals with information about the number of sources reported by antenna 1 and the TDOA of the emission of certain source reported by antenna 2. Since it is possible to classify the signals arriving to antenna 1 depending on their origin, the paired signals received by antenna 2 are automatically assigned to the same origin. The result is as many groups of pairs of signals as sources exist. Finally, once the signals are separated, every group is analysed individually to determine the position of the corresponding source applying any of the two methods presented in Section 2.
The alternative method introduced in Section 2 to acquire the signals moving both antennas would pose an important difficulty in determining the origin of the pulse since the information about the number of sources reported by the stationary antenna would be missing. In this case, there would be differences in terms of power spectral density between the signals received by the two antennas because the path followed from the source would be different. Then, the pair of signals should be grouped based on a statistical analysis of the TDOA finding similar time differences of arrival and hence determining the number of sources and their positions. As partial discharges are a stochastic phenomenon, an antenna can receive events in random times which could lead to mistakes in the association to the appropriate source. Unfortunately, this procedure is less reliable and produces worse results than if one of the antennas is stationary.
In summary, the process would start leaving antenna 1 fixed and rotating antenna 2 acquiring signals with both antennas every angle step. Then, classifying the signals in as many classes as emitting sources are using the information in the power spectral density with the power ratio maps classification method described in [23,28] and summarized in the next section. Finally, signals are extracted from the same source to determine their origin using multilateration techniques.

4. Power Ratio Maps Enhanced with Particle Swarm Optimization (PSO)

The power ratios map is applied to separate signals through their spectral power finding those bands of frequency where their spectra is different [28]. Let [ f 1 L , f 2 L ] and [ f 1 H , f 2 H ] be the limits for the low and high frequency intervals, respectively. The spectral power is calculated in those bands, referred to as the total power of the signal, so every signal would be parameterized with a power ratio at low frequencies, or P R L , and a power ratio at high frequencies, P R H :
P R L = f = f 1 L f 2 L G ( f ) 2 f = 0 f T G ( f ) 2
P R H = f = f 1 H f 2 H G ( f ) 2 f = 0 f T G ( f ) 2
where G ( f ) is the Fourier transform of the signal g ( t ) and f T is the highest frequency of interest of g ( t ) .
Those signals corresponding to the same event would have similar spectra and so would be the parameters P R L and P R H . These signals would form a packed cluster when represented in a two dimensions scatter plot P R L versus P R H . Other events may present differences in these parameters, so the clusters would be plotted separately from the first one.
The selection of the frequency limits for the intervals can be done choosing the best set of frequencies that would group the same type of signals in k packed clusters and gives the largest separation between clusters. The clusters are selected using the k-means algorithm because it is manageable and easy to implement. Then, the distance between the centroids of clusters i and j and the spread of the clusters are simultaneously evaluated using the Mahalanobis distance m i j . The frequency intervals are selected maximizing the minimum distance between the centroids so the objective function M would be:
M = max min i j k m i j
In the case of two sources, k = 2 , and the Equation (8) can be simplified to M = max m 12 .
Particle swarm optimization (PSO) is then applied to obtain M through the selection of a valid set of intervals. PSO is a metaheuristic algorithm with hive intelligence based on the behaviour of animals such as birds or fish in search for food. The algorithm starts deploying a set of particles in the solutions space and moving them towards an optimum with information of the achievements of the swarm. The dimension of the particles is related to the five parameters used to define the PR map, [ f 1 L , f 2 L ] , [ f 1 H , f 2 H ] and f T . In every iteration l the algorithm stores the position of every particle for its personal best solution for M so far, P n , b , and the position of the particle that has achieved the overall best solution for M, P b . The movement of particle n towards the solution for M, v n , is calculated with the following standard PSO algorithm:
v n ( l + 1 ) = v n ( l ) + c 1 U 1 ( 0 , 1 ) [ P n , b ( l ) - P n ( l ) ] + c 2 U 2 ( 0 , 1 ) [ P b ( l ) - P n ( l ) ] , P n ( l + 1 ) = P n ( l ) + v n ( l + 1 ) ,
where U 1 ( 0 , 1 ) , U 2 ( 0 , 1 ) [ 0 , 1 ] 5 are five-dimensional random vectors, with each component independently drawn from a uniform distribution between 0 and 1. Both U 1 ( 0 , 1 ) , U 2 ( 0 , 1 ) randomize the movement of the particles towards their own best and the swarm’s best, respectively. The operator ⊗ multiplies the random numbers by the five coordinates component by component. The parameters c 1 and c 2 describe the balance between the personal influence of the particle and the social influence in the search of the solution. After a certain number of iterations the algorithm converges towards a set of frequencies that maximizes the distance between the two clusters. The outcome of the application of these algorithms is k sets of signals coming from k different sources that can be studied individually to determine the time differences of arrival to the two antennas.

5. Setup and Analysis of the Measurements

Partial discharges were activated winding a grounded enameled copper wire around an unshielded high-voltage cable, see Figure 4. When the voltage is applied there would be a high electric field divergence in that area that would create surface discharges. This procedure is repeated in another section of the cable so there would be two similar sources of partial discharges. The antennas are two simple omnidirectional monopoles 10 cm in length that have good response in the frequency range where partial discharge emits energy which for these test objects is below 400 MHz [24]. In the measurements in two dimensions, the stationary antenna is placed at the origin of a Cartesian coordinate system. The positions of the two emitters are S 1 = ( - 1 , 2 ) m and S 2 = ( 0 , 2 ) m in Cartesian coordinates and S 1 = 2 . 24 116 . 6 m and S 2 = 2 90 m in polar coordinates, Figure 5. The antennas are in the same plane attached to a rod r = 0 . 65 m long and connected to a high-speed digitizing oscilloscope through coaxial cables 5 m long. The sampling frequency was set to 5 GS/s so the sampling time is T s = 200 picoseconds. The distance between the antennas fulfills the relationship r > 10 c T s since 10 c T s = 0 . 6 m so there will be at least 10 TDOA in samples between the signals arriving to the antennas when the antennas and the source are in the same bearing. The antennas were rotated 360 degrees every 10 degrees. To have a statistically reliable TDOA, the oscilloscope acquired 500 signals for every angle step.
Figure 6 presents two PD pulses, one for each source, captured with the antenna 1. Once the signals have been acquired, the analysis process starts calculating their spectral power density arriving to the stationary antenna. Bear in mind that this antenna would receive signals in which the signature of the environmental characteristics of the installation is stable and distinctive for every source while the moveable antenna will be used only to find the TDOA. Figure 7 shows a PR map with two clusters corresponding to the signals arriving to this antenna from two different sources. The average power spectral density of the signals from each cluster is represented in Figure 8 together with the selected intervals in the low and the high frequency ranges to separate the sources. These bands, PRL and PRH, are selected automatically using PSO to maximize the separation of the clusters as explained in Section 4 and results f 1 L = 39 MHz, f 2 L = 50 MHz, f 1 H = 137 MHz, f 2 H = 254 MHz and f T = 280 MHz. Both intervals select more power from the dashed line than from the solid one. This places the dots in the PR map with higher PRH and PRL values obtaining the grey cluster. Additionally, the black cluster has only around 5% of the total power in the low frequency interval, PRL, which places it at the left side of the other cluster. Notice that the differences in the spectra are due to the fact that the paths that follow the signals are not the same since the physical effect of the discharges is quite similar.
The polygon in Figure 7 represents the boundaries of the first cluster selected with Matlab for the analysis of the position of one of the sources. As explained in Section 3 the signals received by antenna 1 are paired with the signals of antenna 2 thanks to the two channels in the acquisition system; so, selecting one of the clusters corresponding to one of the sources will give the signals from that source received by antenna 2.
The TDOA for the 500 pulses in each cluster are plotted versus the angle of rotation in a box plot in Figure 9 and Figure 10. The bottom and top of the boxes in blue represent the boundaries for 25% to 75% of all the values obtained for the TDOA while the red line inside the box represents the median value. The red crosses outside the box are estimations of the TDOA given by the algorithm for some of the pulses that do not have clear onsets. These cases give TDOA outside the boundaries of the box so they can be considered statistical outliers though they are also considered in the calculation of the median value and the boundaries of the box. Admitting that it may be thought that the number of outliers is high, it can be clearly seen that the boxes are very narrow which means that most of the results are concentrated around the median. Figure 11 shows the fitting of the periodic functions considering only the median of all TDOA in every angle of rotation. The numerical results for the source 1 after the fitting with Equation (4) are x ^ s 1 = - 0 . 04 m and y ^ s 1 = 2 . 11 m while the exact values, considering the layout in Figure 5 were x s 1 = 0 m and y s 1 = 2 m. The shift in Euclidean distance between the actual position of the source and the estimation is 12 cm. The fitting in Figure 11 is assessed with the mean square error considering the median of the TDOA yielding ϵ 1 = 0 . 1061 samples or ϵ 1 = 0 . 64 in centimeters which are also very low. The results for the source 2 are x ^ s 2 = - 0 . 99 m and y ^ s 2 = 2 . 07 m while the actual values are x s 1 = - 1 m and y s 1 = 2 m which are also very close to the estimated position. The error between the actual source position and the estimated one is 7 . 1 cm while the mean square error of the fitting considering the median of the TDOA is ϵ 2 = 0 . 3627 samples or ϵ 2 = 2 . 2 cm. As additional information, the estimated zero crossings obtained with the method presented in Section 2.1 for source 1 and 2, respectively, are ζ 11 = 34 . 17 , ζ 12 = - 161 . 18 and ζ 21 = 8 . 18 , ζ 22 = 167 . 5 which are very close to the theoretical ones in Figure 5. Though the error is well below 1 in both angles in the case of source 1, the intersection of the lines is shifted from its actual position: x ^ s 1 = - 1 . 09 m and y ^ s 1 = 2 . 18 m which corresponds to an error of 20.1 cm. In the case of the source 2, x ^ s 2 = 0 . 08 m and y ^ s 2 = 1 . 81 m placing the source at position shifted 20.6 cm from the actual one.

6. Localizing in Three Dimensions

The extension of the method to locate sources in three dimensions is straightforward, placing one of the antennas at a different height so there can be a shift of the TDOA received by that antenna with the information in the Z axis. The easiest setup would result in elevating the stationary antenna as shown in Figure 12. The method in Section 2.1 will be discarded henceforth since the results obtained in 2D for two sources have not been very accurate. Moreover, if one of the antennas is in a different plane, the zero crossings in Figure 3 would not correspond to a position in which the rod is perpendicular to the direction of the source.
The equations derived in Section 2.2 have to be modified to use three components. The starting point also considers that the TDOA in meters can be calculated subtracting the distances, d = d 1 - d 2 , where d 1 and d 2 are again the distances from antennas 1 and 2 to the source, respectively. From Figure 13, d s can be defined as the distance of the source to the vertical axis Z, d s = d 1 cos β , with β = arcsin [ ( z 1 - z s ) / d 1 ] . Then, the equations can be written as:
( x s , y s , z s ) = ( d s cos α , d s sin α , z s ) ( x 1 , y 1 , z 1 ) = ( 0 , 0 , z 1 ) d 1 = ( d s cos α , d s sin α , z s - z 1 ) ( x 2 , y 2 , z 2 ) = ( r cos θ , r sin θ , 0 ) d 2 = ( d s cos α - r cos θ , d s sin α - r sin θ , z s )
where ( x s , y s , z s ) are the coordinates of the source, ( x 1 , y 1 , z 1 ) and ( x 2 , y 2 , z 2 ) the positions of antennas 1 and 2 and r is the length of the rod that connects antenna 2 to the base of antenna 1. Therefore, d is calculated as:
d = d 1 - ( d s cos α - r cos θ ) 2 + ( d s sin α - r sin θ ) 2 + z s 2
or,
d = d 1 - d s 2 + r 2 - 2 r d s cos ( α - θ ) + z s 2
In these equations, there are only three unknown parameters: d 1 , α and z s since d s is related to d 1 and θ is the rotation angle.
The next steps are the same as in the case of the localization in 2D: The signals from two sources received at the stationary antenna are separated using the PR classification map. This sets the correspondence between the events captured by the stationary antenna and the sources. Since every signal in the PR map is paired to another signal acquired by the movable antenna, it is possible to determine the TDOA for each source. These time differences are again estimated with the Hinkley criterion at different positions of antenna 2 and plotted versus the angle of rotation in a box plot. The results are presented in Figure 14 and Figure 15. Then, Equation (12) is used to fit the plot of median values of the TDOA for a complete rotation of the rig to obtain the plots shown in Figure 16. The resulting fitting parameters yield the Cartesian coordinates of the source using Equation (10).
There is also an alternative formulation of Equation (12) using directly Cartesian coordinates:
d = x s 2 + y s 2 + ( z s - z 1 ) 2 - ( x s - x 2 ) 2 + ( y s - y 2 ) 2 + ( z s - z 2 ) 2
In this case, the position of the source ( x s , y s , z s ) is readily obtained from the fitting. The results for the median TDOA have been fitted again with Matlab with a non-linear least squares algorithm using the relations in Equations (10) and (13).
The localization is satisfactory as shows Table 1 and the mean square error of the fitting is very low ϵ 1 = 0 . 2277 samples or ϵ 1 = 1 . 37 cm for the source 1 and ϵ 2 = 0 . 1247 samples or ϵ 2 = 0 . 75 cm for the source 2.

7. Conclusions

This paper proposes a new method to determine the position of pulsed sources such as the emission of partial discharges. The localization is based on the fact that the signals that arrive to a stationary antenna from different sources have different signatures. The fingerprint has information about the direct path from the source to the antenna, reflections and attenuation. This information is used to separate the sources and study the time difference of arrival to the receivers to determine their position. The other antenna is rigged to the stationary one and describes a circumference to have different time differences of arrival. The function that relates the TDOA with the rotation angle can be derived mathematically based on the geometry of the setup and the position of the sources. Then, the experimental data is fitted to this function to obtain the position of the sources. The usability of the method has been demonstrated for two sources in two and three dimensions obtaining good results. The TDOA are determined with the Hinkley criterion based on the change of the accumulated energy in the signals and this step is critical to have accurate results. The value of the selected TDOA in every rotation angle is based on the median value because there is a certain number of outliers that would render the average value useless. These outliers correspond to signals where the Hinkley criterion fails to find a correct value because the first frontwave of the pulse is not clear or the noise triggers an early picking of the onset of the pulse. In any case, the number of TDOA close to the correct value is high and the experimental functions are fitted accurately. A previous denoising of the signals with the discrete wavelet transform or the PR maps (such as in [12]) may help in improving the accuracy in the determination of the position of the source since the pickings of the onsets with the Hinkley criterion would be better. Another source of random error is the accurate positioning of the movable antenna because the movement is not still automatized. We expect to have this resolved with a stepper motor capable of rotating accurately less than one degree per step.
The practical application of this setup with two antennas in an open-air substation is straightforward since the area that occupies the rotating antenna is not large and would not interfere with the assets. The rig can be placed anywhere in the installation as long as the distances to the farther equipment stay inside the range of the antennas so the signal-to-noise ratio is appropriate to detect the discharges. It might be necessary to measure in several places if the signals acquired are highly distorted due to metallic structures in the line-of-sight. In each place, a complete rotation of the mobile antenna is recommended to have a good fitting with the theoretical curve. If the separation between sources is very short, the proposed setup would also be able to detect them as long as it is larger than the minimum TDOA which in the case of this paper is 6 cm. The problem when two sources are very close would be having two different propagation paths to distinguish the signals with the PR method. Therefore, two very close PD in the same side of an insulator, for instance, might follow a similar propagation path and the separation method may fail; although, two close PD on different sides of the insulator would be clearly identified.

Author Contributions

The idea of paper was conceived by G.R. and J.M.F. The methodology and the experimental validation was proposed by G.R., J.M.M.-T. and J.M.F. during the different stages of the acquisition process and analysis. The analytical solution, the programming of the algorithms and the formal analysis of the outcomes was done by G.R. and J.M.F. G.R. prepared the original first version of the manuscript. All authors contributed to the revision and edition of the manuscript. Funding and resources were obtained by G.R. and J.M.M.-T.

Funding

The work done in this paper has been funded by the Spanish Government (MINECO) and the European Regional Development Fund (ERDF) under contract DPI2015-66478-C2-1-R (MINECO/FEDER, UE).

Acknowledgments

Tests were done in the High-Voltage Research and Test Laboratory (LINEALT) at Universidad Carlos III de Madrid.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, K.; Kwon, H.; You, K. Laser-Interferometric Broadband Seismometer for Epicenter Location Estimation. Sensors 2017, 17, 2423. [Google Scholar] [CrossRef] [PubMed]
  2. Moreau, J.; Ambellouis, S.; Ruichek, Y. Fisheye-Based Method for GPS Localization Improvement in Unknown Semi-Obstructed Areas. Sensors 2017, 17, 119. [Google Scholar] [CrossRef] [PubMed]
  3. Xiong, R.; van Waasen, S.; Rheinlnder, C.; Wehn, N. Development of a Novel Indoor Positioning System With mm-Range Precision Based on RF Sensors Network. IEEE Sens. Lett. 2017, 1, 1–4. [Google Scholar] [CrossRef]
  4. Zou, Y.; Liu, H.; Xie, W.; Wan, Q. Semidefinite Programming Methods for Alleviating Sensor Position Error in TDOA Localization. IEEE Access 2017, 5, 23111–23120. [Google Scholar] [CrossRef]
  5. Kreuger, F. Partial Discharge Detection in High-Voltage Equipment; Butterworths: Kent, UK, 1989. [Google Scholar]
  6. IEEE Guide for Field Testing and Evaluation of the Insulation of Shielded Power Cable Systems Rated 5 kV and Above; IEEE-400:2012, Standard; IEEE: New York, NY, USA, 2012.
  7. Rotating Electrical Machines—Part 27-2: On-Line Partial Discharge Measurements on the Stator Winding Insulation of Rotating Electrical Machines; IEC-TS-60034-27-2, Standard; International Electrotechnical Comission: Geneva, Switzerland, 2012.
  8. Stone, G.C.; Sedding, H.G.; Chan, C. Experience With Online Partial-Discharge Measurement in High-Voltage Inverter-Fed Motors. IEEE Trans. Ind. Appl. 2018, 54, 866–872. [Google Scholar] [CrossRef]
  9. AJ, C.; Salam, M.; Rahman, Q.; Wen, F.; Ang, S.; Voon, W. Causes of transformer failures and diagnostic methods—A review. Renew. Sustain. Energy Rev. 2018, 82, 1442–1456. [Google Scholar] [CrossRef]
  10. Li, P.; Zhou, W.; Yang, S.; Liu, Y.; Tian, Y.; Wang, Y. A Novel Method for Separating and Locating Multiple Partial Discharge Sources in a Substation. Sensors 2017, 17, 247. [Google Scholar] [CrossRef] [PubMed]
  11. Mishra, D.K.; Sarkar, B.; Koley, C.; Roy, N.K. An unsupervised Gaussian mixer model for detection and localization of partial discharge sources using RF sensors. IEEE Trans. Dielectr. Electr. Insul. 2017, 24, 2589–2598. [Google Scholar] [CrossRef]
  12. Robles, G.; Fresno, J.M.; Martínez-Tarifa, J.M. Separation of Radio-Frequency Sources and Localization of Partial Discharges in Noisy Environments. Sensors 2015, 15, 9882. [Google Scholar] [CrossRef] [PubMed]
  13. Xu, Y.; Zhou, J.; Zhang, P. RSS-Based Source Localization When Path-Loss Model Parameters are Unknown. IEEE Commun. Lett. 2014, 18, 1055–1058. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Upton, D.; Jaber, A.; Ahmed, H.; Khan, U.; Saeed, B.; Mather, P.; Lazaridis, P.; Atkinson, R.; Vieira, M.F.Q.; Glover, I.A. Multiple source localization for partial discharge monitoring in electrical substation. In Proceedings of the 2015 Loughborough Antennas Propagation Conference (LAPC), Loughborough, UK, 2–3 November 2015; pp. 1–4. [Google Scholar]
  15. Li, Z.; Luo, L.; Zhou, N.; Sheng, G.; Jiang, X. A Novel Partial Discharge Localization Method in Substation Based on a Wireless UHF Sensor Array. Sensors 2017, 17, 1909. [Google Scholar] [CrossRef] [PubMed]
  16. Li, Z.; Luo, L.; Sheng, G.; Liu, Y.; Jiang, X. UHF partial discharge localisation method in substation based on dimension-reduced RSSI fingerprint. IET Gener. Transm. Distrib. 2018, 12, 398–405. [Google Scholar] [CrossRef]
  17. Gao, S.; Zhang, Y.; Xie, Q.; Kan, Y.; Li, S.; Liu, D.; Lü, F. Research on Partial Discharge Source Localization Based on an Ultrasonic Array and a Step-by-Step Over-Complete Dictionary. Energies 2017, 10, 593. [Google Scholar] [CrossRef]
  18. Wielandt, S.; Strycker, L.D. Indoor Multipath Assisted Angle of Arrival Localization. Sensors 2017, 17, 2522. [Google Scholar] [CrossRef] [PubMed]
  19. Yongfen, L.; Xiaohu, X.; Fei, D.; Xiao, T.; Yanming, L. Comparison of DOA Algorithms Applied to Ultrasonic Arrays for PD Location in Oil. IEEE Sens. J. 2015, 15, 2316–2323. [Google Scholar] [CrossRef]
  20. Zhu, M.X.; Wang, Y.B.; Liu, Q.; Zhang, J.N.; Deng, J.B.; Zhang, G.J.; Shao, X.J.; He, W.L. Localization of multiple partial discharge sources in air-insulated substation using probability-based algorithm. IEEE Trans. Dielectr. Electr. Insul. 2017, 24, 157–166. [Google Scholar] [CrossRef]
  21. Robles, G.; Fresno, J.M.; Sánchez-Fernández, M.; Martínez-Tarifa, J.M. Antenna Deployment for the Localization of Partial Discharges in Open-Air Substations. Sensors 2016, 16, 541. [Google Scholar] [CrossRef] [PubMed]
  22. Fresno, J.M.; Robles, G.; Martínez-Tarifa, J.M. Planar Localization of Radio-Frequency or Acoustic Sources with Two Receivers. Proceedings 2018, 2, 103. [Google Scholar] [CrossRef]
  23. Fresno, J.M.; Ardila-Rey, J.A.; Martínez-Tarifa, J.M.; Robles, G. Partial discharges and noise separation using spectral power ratios and genetic algorithms. IEEE Trans. Dielectr. Electr. Insul. 2017, 24, 31–38. [Google Scholar] [CrossRef]
  24. Robles, G.; Sánchez-Fernández, M.; Albarracín Sánchez, R.; Rojas-Moreno, M.; Rajo-Iglesias, E.; Martínez-Tarifa, J. Antenna Parametrization for the Detection of Partial Discharges. IEEE Trans. Instrum. Meas. 2013, 62, 932–941. [Google Scholar] [CrossRef] [Green Version]
  25. Robles, G.; Fresno, J.M.; Giannetti, R. Ultrasonic bone localization algorithm based on time-series cumulative kurtosis. ISA Trans. 2017, 66, 469–475. [Google Scholar] [CrossRef] [PubMed]
  26. Boker, E.; Levy, A.; Goldreich, R. Mitigation of Multipath Distortions for TDOA-Based Geolocation. U.S. Patent 9,841,489, 12 December 2017. [Google Scholar]
  27. Bai, F.; Gagar, D.; Foote, P.; Zhao, Y. Comparison of alternatives to amplitude thresholding for onset detection of acoustic emission signals. Mech. Syst. Signal Process. 2017, 84, 717–730. [Google Scholar] [CrossRef]
  28. Robles, G.; Fresno, J.M.; Martínez-Tarifa, J.M.; Ardila-Rey, J.A.; Parrado-Hernández, E. Partial Discharge Spectral Characterization in HF, VHF and UHF Bands Using Particle Swarm Optimization. Sensors 2018, 18, 746. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Different positions representative of antenna 2 marked with triangles (▼) when turned around the fixed antenna 1 (▲).
Figure 1. Different positions representative of antenna 2 marked with triangles (▼) when turned around the fixed antenna 1 (▲).
Sensors 18 01410 g001
Figure 2. Source and antenna positions to determine the analytic expression of the periodic function θ vs. TDOA.
Figure 2. Source and antenna positions to determine the analytic expression of the periodic function θ vs. TDOA.
Sensors 18 01410 g002
Figure 3. Periodic function representing the TDOA vs. the rotating angle. The zero-crossings correspond to the positions of the rig perpendicular to bearing of the source so the distances from the antennas to the source are the same. The maximum and minimum correspond to positions where the rig is parallel to this bearing so one of the antennas is at the farthest relative position with respect to the source and the other at the closest one.
Figure 3. Periodic function representing the TDOA vs. the rotating angle. The zero-crossings correspond to the positions of the rig perpendicular to bearing of the source so the distances from the antennas to the source are the same. The maximum and minimum correspond to positions where the rig is parallel to this bearing so one of the antennas is at the farthest relative position with respect to the source and the other at the closest one.
Sensors 18 01410 g003
Figure 4. Detail of the section of the high voltage cable (in red) with the wound wire to ground.
Figure 4. Detail of the section of the high voltage cable (in red) with the wound wire to ground.
Sensors 18 01410 g004
Figure 5. Positions of source 1 (left) and source 2 (right) and the theoretical zero-crossing and bearing angles.
Figure 5. Positions of source 1 (left) and source 2 (right) and the theoretical zero-crossing and bearing angles.
Sensors 18 01410 g005
Figure 6. Two characteristic PD signals from the two sources acquired by the antenna 1.
Figure 6. Two characteristic PD signals from the two sources acquired by the antenna 1.
Sensors 18 01410 g006
Figure 7. Signals received by antenna 1 represented in a PR map where the differences in the power spectral density corresponding to two sources distributes the signals into clusters.
Figure 7. Signals received by antenna 1 represented in a PR map where the differences in the power spectral density corresponding to two sources distributes the signals into clusters.
Sensors 18 01410 g007
Figure 8. Average power spectral density of the signals captured by the antenna 1 and the selected intervals that maximize the distance between clusters in the PR map. The solid line corresponds to dark cluster and the dashed plot to the grey cluster.
Figure 8. Average power spectral density of the signals captured by the antenna 1 and the selected intervals that maximize the distance between clusters in the PR map. The solid line corresponds to dark cluster and the dashed plot to the grey cluster.
Sensors 18 01410 g008
Figure 9. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 1 (black cluster).
Figure 9. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 1 (black cluster).
Sensors 18 01410 g009
Figure 10. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 2 (grey cluster).
Figure 10. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 2 (grey cluster).
Sensors 18 01410 g010
Figure 11. Median of all TDOA versus rotation angle and fitting function. Blue plot corresponds to source 1 and red plot to source 2.
Figure 11. Median of all TDOA versus rotation angle and fitting function. Blue plot corresponds to source 1 and red plot to source 2.
Sensors 18 01410 g011
Figure 12. Antenna layout for three dimensional localization.
Figure 12. Antenna layout for three dimensional localization.
Sensors 18 01410 g012
Figure 13. Geometric representation of the setup with the main parameters of the equations.
Figure 13. Geometric representation of the setup with the main parameters of the equations.
Sensors 18 01410 g013
Figure 14. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 1 in 3D localization.
Figure 14. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 1 in 3D localization.
Sensors 18 01410 g014
Figure 15. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 2 in 3D localization.
Figure 15. Box plot of the time differences of arrival versus angle of rotation for the pulses generated by source 2 in 3D localization.
Sensors 18 01410 g015
Figure 16. Median of all TDOA versus rotation angle in 3D localization with the fitting function. Blue plot corresponds to the fitting function for source 1 and red plot to source 2.
Figure 16. Median of all TDOA versus rotation angle in 3D localization with the fitting function. Blue plot corresponds to the fitting function for source 1 and red plot to source 2.
Sensors 18 01410 g016
Table 1. Results of the estimation of parameters in the 3D localization. The theoretical position of the sources are: ( 0 , 2 , 0 ) m and ( - 1 , 2 , 0 . 3 ) m, for source 1 and 2, respectively.
Table 1. Results of the estimation of parameters in the 3D localization. The theoretical position of the sources are: ( 0 , 2 , 0 ) m and ( - 1 , 2 , 0 . 3 ) m, for source 1 and 2, respectively.
Source 1Source 2
x s = - 0 . 050 m x s = - 0 . 989 m
y s = 2 . 073 m y s = 2 . 027 m
z s = 0 m z s = 0 . 279 m
Euclidean error0.089 m0.036 m

Share and Cite

MDPI and ACS Style

Robles, G.; Fresno, J.M.; Martínez-Tarifa, J.M. Radio-Frequency Localization of Multiple Partial Discharges Sources with Two Receivers. Sensors 2018, 18, 1410. https://doi.org/10.3390/s18051410

AMA Style

Robles G, Fresno JM, Martínez-Tarifa JM. Radio-Frequency Localization of Multiple Partial Discharges Sources with Two Receivers. Sensors. 2018; 18(5):1410. https://doi.org/10.3390/s18051410

Chicago/Turabian Style

Robles, Guillermo, José Manuel Fresno, and Juan Manuel Martínez-Tarifa. 2018. "Radio-Frequency Localization of Multiple Partial Discharges Sources with Two Receivers" Sensors 18, no. 5: 1410. https://doi.org/10.3390/s18051410

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop