This section aims at evaluating the performance of the EMIF method to localize PD sources in power transformers. The section is divided into three subsections. In
Section 3.1, a simple 2D model is used to represent the transformer tank. The windings inside the transformer tank are modeled using metallic cylinders. Using a 2D model leads to a reduced computation time for the forward and backward phases needed in the proposed EMIF method. In
Section 3.2, the performance of the EMIF method is evaluated using a simple 3D model for the transformer. The assumed model is simulated using the commercial software CST-MWS. In both subsections, comparisons between the proposed method and classical EMTR are also provided. Finally, in
Section 3.3, a laboratory transformer model implemented in the CST-MWS software is used to assess the performance of the proposed method to localize a PD source inside a realistic power transformer tank.
3.1. 2D Transformer Model
Figure 2 shows the adopted 2D model to simulate a three-phase transformer. The transformer tank is modeled by a rectangular shape. Its length and width,
l and
w, are equal to 1000 mm and 500 mm, respectively. The material of the transformer walls is considered to be a perfect electric conductor. The three solid cylindrical windings inside the transformer are modeled as three metallic
a/2 = 100 mm radius circles. The circular barriers in
Figure 2 represent the outer radius of the three-phase transformer windings. The center of the middle winding coincides with the center of the tank. The edge-to-edge distance between two adjacent windings is
b = 50 mm. The considered model for the 2D transformer is symmetric about the
x and the
y axes. The transformer tank is considered to be empty. In other words, no insulating equipment (such as paper pressboard) or magnetic circuit (core) is considered in the transformer tank. Note that, even though the model used here is a simplified transformer model, it has been shown that the performance of the time-reversal method is improved in more complex and heterogeneous media due to the effect of multiple reflections [
25].
In
Figure 2, the locations of the PD source (point #1) and sensor (point #2) are denoted by a black circle and a triangle, respectively. The shown locations of the PD source and sensor in
Figure 2 are arbitrary and they are shown only for illustration purposes. In all the 2D simulations, one sensor and one PD source are considered. To implement the proposed EMIF method, a cosine-modulated Gaussian pulse is considered for the PD source signal, as shown in
Figure 3a. The magnitude of the frequency-domain transform of this cosine-modulated Gaussian pulse is shown in
Figure 3b,c in linear and logarithmic scale, respectively. The center frequency and bandwidth for the signal in
Figure 3 are 3 GHz and 3 GHz, respectively. Depending on the case study, different values for the bandwidth and frequency center are used in this paper.
For the 2D simulations, five case studies (CS#1 to CS#5) are considered with different locations of the source and the sensor (see
Table 1). In the first and second case studies, the windings (see
Figure 2) are not included in the model.
The 2D finite-difference time-domain (FDTD) technique was used to solve Maxwell’s equations in the transverse magnetic mode. In other words, one electric field component (Ez) and two magnetic field intensity components (Hx and Hy) will be used to solve Maxwell’s equations. In the FDTD technique, equally spaced mesh cells along the x and y axes are considered with a length of 2.5 mm. The time step used in the simulations, rounded to one decimal place, is 5.3 picoseconds. It should be noted that the sensor is modeled using a single-point field probe that records the electric field calculated by the FDTD method.
3.1.1. CS#1: Transformer Tank Model without Windings
In this case study, the source and sensor locations are respectively at positions (0.25 m, −0.2 m) and (−0.25 m, 0.2 m). The excitation signal for the source is shown in
Figure 3. The number of iterations in the FDTD simulation is 8192 time steps (corresponding to a time window of 43 ns).
Figure 4a,b show the time- and normalized frequency-domain waveforms of the electric field intensity recorded by the sensor. The resonance shown in
Figure 4a is due to the multiple reflections from the walls of the transformer tank. To apply the proposed EMIF method, the inverted version of the signal received by the sensor in the frequency domain is calculated using Equation (8) with
. The normalized inverted frequency-domain signal received by the sensor is shown in
Figure 4c. This inverted signal is transformed to the time domain using the inverse fast Fourier transform and then back injected into the transformer tank model.
After doing the backward simulation using the FDTD method with the same parameters as in the forward phase, the distribution of the normalized maximum electric field power over the whole time interval is shown in
Figure 5a. In all the figures, the colors represent the intensity of the represented quantity. In that figure, the red circle (o) and black cross (×) show the actual and estimated source locations, respectively. The location of the sensor is shown by the pink square (□) in
Figure 5a. An expanded view of the region close to the source location is shown in the figure inset. The localization error for the proposed EMIF method is zero. In
Figure 5, to estimate the location of the source, the maximum power criterion proposed in [
18] is used. The length of the square mask window around the sensor (used to exclude the sensor from the maximum power search domain) is 25 cells along the
x- and
y-axes.
For the sake of comparison, the distribution of the normalized maximum electric field power over the computational domain obtained using the classical EMTR method is shown in
Figure 5b. An expanded view for the source region is also shown in the figure inset. The red circle (o) and the black cross (×) show the actual and estimated source locations, respectively. The localization error with the EMTR method is 2.5 mm, namely, equal to a mesh cell.
It can be seen that both the EMIF and the EMTR methods can easily localize the source inside the transformer tank with only one sensor. The accuracy of the EMIF method appears to be better than that of the EMTR method for this specific case.
3.1.2. CS#2: Localization of a PD Source Close to a Metallic Wall
In the second case study, CS#2, the performance of the EMIF method is compared to the EMTR method when the PD is located in the immediate vicinity of the metallic transformer wall. In this case, the source is located at point (0.3 m, 0.24 m), 10 cm away from the metallic wall, while the sensor is located at (−0.35 m, −0.2 m). The other parameters are the same as in CS#1.
Figure 6a,b show the distribution of the normalized maximum electric field power over the computational domain obtained by using the EMIF and the EMTR methods. In this figure, the red circle (o) and the black cross (×) show the actual and estimated source locations, respectively. The location of the sensor is denoted by the pink square (□). The length of the square mask windows used in the EMIF and EMTR methods is 20 cells [
18]. An expanded view of the region around the source is shown in the insets to highlight the performance of both methods. The localization errors are, respectively, 10.3 mm and 17.5 mm for the proposed EMIF method and the classical EMTR method. This shows that EMIF performs better than EMTR in locating sources close to metallic walls.
3.1.3. CS#3: Localization of a PD Source between Two Metallic Windings
In the third case study, CS#3, three circular barriers representing the transformer windings are included in the transformer model. The windings are considered perfect electric conductor circles. The considered source is located between two windings. The coordinates of the PD source and the sensor are (0.125 m, 0.0 m) and (−0.25 m, −0.2 m), respectively. All the other parameters of the FDTD simulations are the same as in the first case study, CS#1.
The distributions of the normalized maximum electric field power over the computational domain obtained by the EMIF and EMTR methods are shown in
Figure 7a,b, respectively. In this figure, the location of the sensor is denoted by the pink square (□). The red circle (o) and the black cross (×) show the actual and estimated source locations. An expanded view of the region around the source is shown in the EMIF figure inset.
The localization errors are 2.5 mm and 75.0 mm for the proposed EMIF method and EMTR, respectively. The length of the square mask windows used in both the EMIF and EMTR methods is 20 cells [
18]. This case study shows that the proposed EMIF method is able to solve one of the main problems to localize sources located between the two adjacent metallic objects [
17]. In this case study, the EMIF method shows a better performance compared to the classical EMTR method again.
3.1.4. CS#4: Localization of PD Source behind the Metallic Windings
In this case study, in addition to the presence of circle barriers or windings in the transformer tank, the PD source is behind the farthest metallic winding from the sensor. Note that, in this case study, compared to CS#3, no direct line of sight exists between the sensor and the PD source. In other words, the source of PD is hidden from the view of the sensor. The PD location coordinates are (0.25 m, −0.125 m), and the sensor is located at the point (−0.25 m, 0.2 m). The PD source signal in this case study is assumed to be a cosine-modulated Gaussian pulse with a bandwidth and a frequency center of 1.5 GHz. Due to the reduction in bandwidth in this example and to the location of the PD source behind the transformer winding, the localization is more challenging compared to the previous cases.
The number of iterations for CS#4 is 16,384 time steps, corresponding to around 87 ns. All the other parameters are the same as CS#1.
Figure 8a,b show the simulation results obtained by the EMIF method and the EMTR method, respectively. In this figure, the location of the sensor is denoted by the pink square (□). The red circle (o) and the black cross (×) show the actual and estimated source locations. An expanded view of the region around the source for the EMIF method is shown in the figure inset.
The localization errors for the proposed EMIF method and EMTR are 12.7 mm and 25.7 mm, respectively. Again, the EMIF method shows a better performance compared to the classical EMTR method. The lengths of the square mask windows used in both the EMIF and EMTR methods are 30 and 20 cells along the x- and y-axes, respectively [
18].
3.1.5. CS#5: Localization of a PD Source between the Metallic Windings
In the fifth case study, the location of the PD source and the sensor are (0.125 m, 0.0 m) and (−0.25 m, 0.2 m), respectively. The considered center frequency and bandwidth of the PD source waveform are 2.0 GHz and 1.5 GHz, respectively. All the other parameters are the same as in CS#4. Compared to CS#3, the center frequency and bandwidth of the excitation signal are reduced in this case study.
Figure 9a,b show the simulation results obtained by the EMIF and EMTR methods. In this figure, the location of the sensor is denoted by the pink square (□). The red circle (o) and the black cross (×) show the actual and estimated source locations. An expanded view of the region around the source for the EMIF method is shown in the figure inset.
The localization errors are 25.1 mm and 57.5 mm for the proposed EMIF method and for EMTR, respectively. Again, the EMIF method outperforms the classical EMTR method. The lengths of the square mask windows used in the EMIF and the EMTR methods are 30 and 20 cells along the x- and y- axes, respectively [
18].
The degradation in the accuracy of the EMTR localization when the PD is located near a metallic plane or object (within
, where
is the minimum wavelength of the source) has been discussed in [
17]. For the same reason, the EMIF location accuracy is also reduced but not as much as in the EMTR method, as indicated by the results associated with case studies CS#2, CS#3, and CS#5. This is essentially because the upper 3 dB cutoff frequency used for the EMIF method is higher than for the EMTR method (more energy in the higher frequencies), making
In all the case studies, the obtained localization error by the EMIF method is significantly lower than that obtained using the EMTR method. The values of the localization error by both methods are summarized in
Table 2.
3.2. A Simple 3D Transformer Model
In this section, a simple 3D model for the transformer tank is considered as shown in
Figure 10. The size of the considered 3D tank is 1000 × 500 × 500 mm
3 in the direction of the x-axis, the y-axis, and the z-axis, respectively. The tank walls are assumed to be perfect electric conductors. Here, the 3D simulation commercial software CST-MWS is used to implement both the forward and backward simulations in the EMIF and EMTR methods. CST-MWS uses the finite integration technique (FIT) to solve Maxwell’s equations in the time domain.
To model the sensor, a dipole antenna with a center frequency of 1 GHz is considered. The length, diameter, and gap parameter for this dipole antenna are 135 mm, 5 mm, and 20 mm, respectively. As shown in
Figure 10, the direction of the dipole antenna is along the z-axis. A 72-Ohm resistance is used at the input of the dipole antenna. The center of the dipole antenna is placed at point (−0.25 m, −0.1 m, 0.15 m) inside the transformer tank.
To model the PD source, a discrete port is modeled inside the transformer tank. The center of the PD source is (0.3 m, 0.16 m, −0.1 m). A Gaussian pulse with a bandwidth of 3 GHz and 30 ns delay with respect to the one shown in
Figure 3a is used to excite the PD source. The minimum and maximum frequencies in the CST-MWS are 0 and 3 GHz, respectively. The length and direction of the PD source are 20 mm and along the z-axis, respectively.
In the backward phase for both EMIF and EMTR, the electric field intensities over three orthogonal planes containing the PD source are monitored. To localize the PD source, the maximum field power criterion [
18] is used over these three planes. In the backward phase, the PD source is removed from the simulation model. The number of mesh cells and the time step used in the CST-MWS software are 1,292,000 and 4.89 ps, respectively.
To implement the EMIF method, the simulation in CST-MWS is done for the forward phase. The recorded signal at the output of the dipole antenna is exported from CST-MWS and transformed to the frequency domain.
Figure 11a,b show the time-domain waveform recorded by the dipole antenna and its normalized frequency-domain magnitude spectrum. The reason for the time fluctuations in the time-domain signal is the multiple reflections from the metallic walls of the transformer tank.
The signal received by the dipole antenna is inverted according to Equation (8) and transformed back into the time domain. The normalized magnitude spectrum of the inverted signal is shown in
Figure 12c. Finally, in the backward phase, the time-domain version of the inverted signal is back injected into the CST-MWS model of the transformer without the discrete port used to model the PD source.
Figure 12a–c show the distribution of the normalized maximum electric field power over the computational domain obtained by the EMIF method in the three orthogonal cut-planes. In these figures, the estimated PD location is shown by the black cross (×).
If the center of the discrete port is considered as the actual location of the PD source, the 3D localization error in the EMIF method is 5 mm. The 3D localization error is defined as the 3D distance between the center of the discrete port used for the PD source and the estimated location obtained by the maximum field criterion. As in the 2D case studies presented earlier, here we use 0.9 as the value of the parameter in the EMIF method.
Figure 13a–c show the distribution of the normalized maximum electric field power over the computational domain in three orthogonal cut-planes, obtained by the classical EMTR method. In these figures, the estimated PD location is shown by the black cross (×).
The 3D localization error in the EMTR method is about 15.1 mm. As for the 2D simulations, EMIF appears to be more accurate than EMTR in locating PD sources in 3D.
It is obvious from
Figure 11b that the received signal has no significant frequency content beyond 3 GHz. On the other hand, as can be seen from
Figure 11c, the inverted signal in the frequency domain has significant frequency content up to 4.5 GHz. This increase in the frequency bandwidth is one of the main reasons explaining why EMIF outperforms EMTR.
A comment is in order on the choice of the parameter
in Equation (8). So far, a value
was adopted in the EMIF method simulations. This choice was based on reference [
24], in which by changing the value of
, the temporal focus is calculated and it is shown that, by setting
, the maximum temporal focusing can be achieved. To evaluate the effect of the parameter
on the performance of the EMIF method, the last 3D case studies are repeated for the EMIF method for different values of
. In these simulations, the parameter
is varied from 0.1 to 100.
Figure 14 shows the distribution of the normalized maximum electric field power over the whole time interval in the computational domain for different values of
from 0.1 to 100 in the EMIF method (see Equation (8)). It can be seen that the variation of
can affect the focal spot and the localization error in the EMIF.
By increasing the parameter , both the focal spot sharpness and the accuracy are degraded. Based on the authors’ experience and multiple simulations (not included in this paper), the range for that yields the best EMIF results is (0.1–5). It should be noted that more investigation is needed to determine this range for other applications.
3.3. A More Realistic 3D Transformer Model
In the last case study, a more realistic transformer model is adopted. The model consists of a transformer tank [
26,
27] with a size of 1200 × 800 × 730 mm, as shown in
Figure 15a. The model allows for the installation of up to five electromagnetic antennas to record the PD-emitted signals: one on the sidewall, two in the oil drain, and two on the top cover, as shown in
Figure 15b.
Since the PD in the laboratory tests was produced by the electrode system installed in the center of the tank, the PD source modeled in CST-MWS was also located in this position. A Gaussian signal with a bandwidth of 3 GHz was used to excite the PD in this location. The PD signal was injected through a discrete port in CST-MWS in the forward phase.
Due to the requirement of high memory resources for the backward phase to record the electric field inside the tank, the monitoring region was restricted to a 40-cm3 cubic region around the electrode system.
Figure 16 shows the distribution of the normalized maximum electric field power over the whole time interval for a square with a length of 40 cm around the electrode system obtained by the EMIF method. In this figure, the red crosses (×) show the estimated locations of the PD source obtained by the maximum field power criterion. It can be seen that the localization error is near zero (5 mm). The actual PD location is considered in the center of the gap used to excite the electrode system. To excite the electrode system shown in
Figure 17, a discrete edge port is used in the CST-MWS model. It should be noted that the same results have been obtained for the EMTR method. Therefore, the EMTR results have not been presented here.