**1. Introduction**

In the recent years compact ultra-wideband (UWB) antennas have received a great attention [1-3] due to the allocation of various frequency bands to UWB systems by Federal Communi‐ cations Committee (FCC) [4]. One of the most interesting aspects of these antennas is their aptitude of detecting electromagnetic (EM) transients with frequency content up to the very high frequency (VHF) range. EM transient phenomena with the above mentioned frequency content are those encountered, for example, in partial discharge (PD) detection. This diagnostic method is now widely used to identify defects taking place in the insulation systems. Wireless systems offer the possibility to achieve the whole shape of PD signals radiated from the source of the PD defects with less distortion due to the electric transmission path.

Modern diagnostic procedures require clear PD patterns for the identification of defects generating PD because different defects can affect differently the insulation reliability, [5-7]. Different pulsating sources can be simultaneously active during a PD measurement session and mixed PD patterns can be recorded. Thus, an effective separation of mixed PD patterns into sub-patterns each one pertinent to a specific noise or PD source typology, is a fundamental task to avoid wrong defect identification, [8-9]. One of the possible approaches in signal separation is based on the assumption that the same defect generates similar waveforms and features derived from signals grouped by similarity can be adopted to identify the defect or noise source. Thus, the dynamic characteristics of the antenna probe must be designed/ evaluated in order to fit the above mentioned requirements.

Both the time and the frequency domain can be used for separation purposes. Recently, the use of the Auto-Correlation Function (ACF) has been proposed for its ability to synthesize both the time and frequency domain features as well as to be less affected by superimposed highfrequency noise, random truncation of the pulse-tail (frequency leakage) and different trigger activation, [10].

© 2014 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The ACF has been adopted here to evaluate the performances of different EM probes in order to design the more suitable geometry to build a portable low-cost device for PD detection. Monopole, triangular and spherical antennas were simulated by means of the surface Method of Moments (MoM) in the frequency domain, [11-14]. The transmitting system is modeled by a power electronic device with a fault current between two metal plates. The ACF shapes of the simulated, transmitted and received signals, have been compared using the Pearson correlation index to verify which sensor is providing the best fidelity among the three.

the conduction band. Dielectric materials are divided into three categories named solid, liquid and gaseous. The main electrical properties that characterize the behaviour of insulating

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

179

**•** *Electric permittivity εr=ε/ε0*; where *ε* is the permittivity of a generic material and *ε0* is the

**•** *Loss factor*, defined as tgδ with δ=90°-φ, the loss angle of the material, with φ the phase shift between voltage and current, which, in the ideal case is equal to 90 ° while, in a real dielectric,

**•** *Dielectric Rigidity* (RD):value of the voltage gradient required to cause discharge through the dielectric under examination or the limit value of the electric field beyond which produces a conduction of electricity (electrical discharge) through the dielectric material. It is commonly measured in kV/mm. The discharge gradient is strongly influenced bya very high number of different parameters and by the test procedure, in particular by the shape

The Test procedures for these characteristic parameters depending on the dielectric materials state are often very different. To better understand the partial discharge phenomenon, the discharge base theory in gas is synthetically presented in the next section, while extended information are reported in [15, 16]. Gaseous insulation, especially ambient air, is the most

The occurrence of partial discharges in electrical equipment has been identified as one of the main causes of breaking devices and reduction in useful life of insulation. For many years, studies on the physical and chemical aspects of PD have been thoroughly addressed and news discharge models were introduced. The principle of discharge starting mechanism, for small distances, was introduced by Towsend [15, 16] and it is based on the coefficient of ionization of the gas subjected to an electric field. For distances more than a few millimetres, Towsend theory is wrong, since discharge develops in a times lower than necessary ones for electrons to cross the distance between electrodes and form an avalanche. To explain the formation of the discharge's mechanism over long distances, several assumptions were made based on phenomenon such as the formation of individual electronic avalanches or photo-ionization phenomena. Meek and Loeb, and independently Rather, proposed the Streamer theory. The equations that describe phenomena involving electrons and ions, avalanches and streamers,

It is possible to macroscopically represent the discharge phenomenon with a circuit model. This one is representing the system cavity-electrodes in the dielectric. The Whitehead model is exhibiting three capacitors to represent a cavity embedded in a dielectric. Fig. 1 details this

are too complex to be dealt with in a few lines and can be found in [15,16].

of electrodes, by the voltage applied and by the way the limit voltage is reached.

**•** *Intrinsic/extrinsic conductivityσ* and *γ* are about 10-5÷10-18 [S/m].

is less than 90 ° due to their leaks.

**2.1. Discharge theory**

model [17].

commonly used insulation in high voltage engineering.

vacuum dielectric permittivity with value 8∙85418781762 ×10-12 [F/m].

materials are:

This chapter focuses on the evaluation of the performances of different antenna sensors suitable for Partial Discharge (PD) measurements. Monopole, triangular and spherical antennas were simulated by means of the surface method of moments. The transmitting system is modeled by a power electronic device with a fault current between two metal plates. The shape of the simulated, transmitted and received signals, has been compared to verify the sensor that provides the best fidelity among the three. The auto-correlation function and the Pearson correlation index are adopted here for the comparison. A discussion on the dynamic characteristic of the different antenna probes and their use in different application is proposed.

#### **2. Partial discharges**

The increasing necessity of high voltage levels in high voltage (HV) power systems has driven to a stronger demand of insulating materials with high electrical performances at affordable costs. Since synthetic polymers partially meet these requirements [15, 16], they have been studied quite a bit in the most recent literature. However, they are affected by some aging problems also the PD phenomenon is a major cause of failure in the insulation systems of electrical apparatus. These discharges take place in consequence of unavoidable local defects produced by the industrial manufacturing process, and promote local erosion of the material that may cause electric breakdown of the component over time. The difficulty of performing reliable life predictions and suitable evaluations of their reliability in presence of such degradation phenomena, hindered a wide diffusion of HV components made of epoxy materials. The diagnostic techniques of recognition, identification and classification of discharge phenomena are therefore fundamental for the evaluation of the time reliability of the electrical equipment, especially for high voltages. In order to study and understand the behaviour of different materials, the first step is to perfectly know the dielectric characteris‐ tics.Materials are divided into three basic categories that are metals, ceramics and polymers. Furthermore, the combination of two or more of these elements gives rise to materials employed in common applications normally named composite materials. Successively, another division is based on the type of atomic bond and on the basis of the crystal structure that composes them. Besides, in solid state physics, the electronic band structure of a solid describes the range of energies that an electron of a certain material is "allowed" or "forbidden" to possess. The electronic structure of a material greatly is influencing the materials behaviour and the energy gap between the valence band. The conduction is used to classify the materials on the electronic characteristics, in conductors, semiconductors and insulators. The insulators materials or dielectric materials present a wide forbidden zone between the valence band and the conduction band. Dielectric materials are divided into three categories named solid, liquid and gaseous. The main electrical properties that characterize the behaviour of insulating materials are:


The Test procedures for these characteristic parameters depending on the dielectric materials state are often very different. To better understand the partial discharge phenomenon, the discharge base theory in gas is synthetically presented in the next section, while extended information are reported in [15, 16]. Gaseous insulation, especially ambient air, is the most commonly used insulation in high voltage engineering.

#### **2.1. Discharge theory**

The ACF has been adopted here to evaluate the performances of different EM probes in order to design the more suitable geometry to build a portable low-cost device for PD detection. Monopole, triangular and spherical antennas were simulated by means of the surface Method of Moments (MoM) in the frequency domain, [11-14]. The transmitting system is modeled by a power electronic device with a fault current between two metal plates. The ACF shapes of the simulated, transmitted and received signals, have been compared using the Pearson correlation index to verify which sensor is providing the best fidelity among the three.

This chapter focuses on the evaluation of the performances of different antenna sensors suitable for Partial Discharge (PD) measurements. Monopole, triangular and spherical antennas were simulated by means of the surface method of moments. The transmitting system is modeled by a power electronic device with a fault current between two metal plates. The shape of the simulated, transmitted and received signals, has been compared to verify the sensor that provides the best fidelity among the three. The auto-correlation function and the Pearson correlation index are adopted here for the comparison. A discussion on the dynamic characteristic of the different antenna probes and their use in different application is proposed.

The increasing necessity of high voltage levels in high voltage (HV) power systems has driven to a stronger demand of insulating materials with high electrical performances at affordable costs. Since synthetic polymers partially meet these requirements [15, 16], they have been studied quite a bit in the most recent literature. However, they are affected by some aging problems also the PD phenomenon is a major cause of failure in the insulation systems of electrical apparatus. These discharges take place in consequence of unavoidable local defects produced by the industrial manufacturing process, and promote local erosion of the material that may cause electric breakdown of the component over time. The difficulty of performing reliable life predictions and suitable evaluations of their reliability in presence of such degradation phenomena, hindered a wide diffusion of HV components made of epoxy materials. The diagnostic techniques of recognition, identification and classification of discharge phenomena are therefore fundamental for the evaluation of the time reliability of the electrical equipment, especially for high voltages. In order to study and understand the behaviour of different materials, the first step is to perfectly know the dielectric characteris‐ tics.Materials are divided into three basic categories that are metals, ceramics and polymers. Furthermore, the combination of two or more of these elements gives rise to materials employed in common applications normally named composite materials. Successively, another division is based on the type of atomic bond and on the basis of the crystal structure that composes them. Besides, in solid state physics, the electronic band structure of a solid describes the range of energies that an electron of a certain material is "allowed" or "forbidden" to possess. The electronic structure of a material greatly is influencing the materials behaviour and the energy gap between the valence band. The conduction is used to classify the materials on the electronic characteristics, in conductors, semiconductors and insulators. The insulators materials or dielectric materials present a wide forbidden zone between the valence band and

**2. Partial discharges**

178 Progress in Compact Antennas

The occurrence of partial discharges in electrical equipment has been identified as one of the main causes of breaking devices and reduction in useful life of insulation. For many years, studies on the physical and chemical aspects of PD have been thoroughly addressed and news discharge models were introduced. The principle of discharge starting mechanism, for small distances, was introduced by Towsend [15, 16] and it is based on the coefficient of ionization of the gas subjected to an electric field. For distances more than a few millimetres, Towsend theory is wrong, since discharge develops in a times lower than necessary ones for electrons to cross the distance between electrodes and form an avalanche. To explain the formation of the discharge's mechanism over long distances, several assumptions were made based on phenomenon such as the formation of individual electronic avalanches or photo-ionization phenomena. Meek and Loeb, and independently Rather, proposed the Streamer theory. The equations that describe phenomena involving electrons and ions, avalanches and streamers, are too complex to be dealt with in a few lines and can be found in [15,16].

It is possible to macroscopically represent the discharge phenomenon with a circuit model. This one is representing the system cavity-electrodes in the dielectric. The Whitehead model is exhibiting three capacitors to represent a cavity embedded in a dielectric. Fig. 1 details this model [17].

**2.2. Partial discharge characterization**

expressed in picocoulombs (pC).

*and time ti*

*Phase angle ϕ<sup>i</sup>*

interval and the duration of this time interval.

discharge pulse and T is the period of the test *voltage.*

A partial discharge is a stochastic phenomenon that involves many chemical and physical aspects of materials proprieties. The PD detection over time has become an indispensable method for the evaluation of insulation systems in the realization phase of new devices but also for evaluate reliability in work. The discharge manifests itself through charge displace‐ ment, acoustic waves, lights and electromagnetic field. The wired detection of charge dis‐ placement was the most frequently measurement method chosen before the introduction of antenna probes. The basic discharge detection circuit, in accordance with the International Electrotechnical Commission (IEC) Standard, is shown in fig. 3 [19]. The coupling capacitor provides a low impedance path for high-frequency discharge currents which cause voltage pulses over the measuring impedance. The pulses are then amplified and stored in a computer for further processing and displayed in different patterns as 3D patterns (phase-amplitudenumber of discharge) and 2D patterns (amplitude-number, phase-number, phase-amplitude).

In the same IEC Standard the main definitions related to PD are defined as follows:

between conductors and which can or cannot occur adjacent to a conductor.

which have been introduced into the test circuit for the purpose of the test.

*of occurrence of a PD pulse:* is *ϕ*<sup>i</sup>

fundamental components because they affect the shape of the impulses.

*Partial discharge (PD)*: Localized electrical discharge that only partially bridges the insulation

*Partial discharge pulse (PD pulse):* Current or voltage pulse that results from a partial discharge occurring within the object under test. The pulse is measured using suitable detector circuits,

*Apparent charge q:* Of a PD pulse is that charge which, if injected within a very short time between the terminals of the test object in a specified test circuit, would give the same reading on the measuring instrument ad the PD current pulse itself. The apparent charge is usually

*Pulse repetition rate n:* Ratio between the total number of PD Pulses recorded in a selected time

between the preceding positive going transition of the test voltage through zero and the partial

The pulses generated by PDs are detected as second order pulses; in Fig.4 the pulse was recorded by means of an oscilloscope and shows a signal proportional to the discharge current. It was detected by means of a 50 Ω resistance in series to the circuit where the defect was present. This kind of pulses are very fast (~50ns) and have a small amplitude (~5mV). As we will see in the following, not only the amplitude of the pulse but also its duration is a funda‐ mental parameter. The most important recent studies that have allowed PD measurements to become a reliable diagnostic system are based on the shape of the pulses, and in particular on the ratio between duration and frequency content. The acquisition techniques (bandwidth of the recording system) and the kind of sensor adopted (HFCT, capacitive couplers) become

=360(ti

/T) where ti is the time measured

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

181

Figure 1. Three capacity of the model of Whitehead. **Figure 1.** Three capacity of the model of Whitehead.

voltage.

Fig. 2 presents the *Vc* voltage applying on the c capacitor. It should be noticed that *Vc* is a fraction of the total voltage V. When *Vc* exceeds the so called inception voltage *Vci*, the discharge starts. After the discharge, *Vc* reaches the discharge extinction voltage value, *Vce*. The global behaviour of the system is affected and the PD activity can be monitored by employing current and voltage measures. Fig. 2 presents the *Vc* voltage applying on the c capacitor. It should be noticed that *Vc* is a fraction of the total voltage V. When *Vc* exceeds the so called inception voltage *Vci*, the discharge starts. After the discharge, *Vc* reaches the discharge extinction voltage value, *Vce*. The global behaviour of the system is affected and the PD activity can be monitored by employing current and voltage measures.

A field approach model that best represents the phenomenon is conceived by Pedersen [18]. Since the charge involved by the discharge inside the cavity, named true charge, cannot be measured, an auxiliary current, named apparent charge, is **Figure 2.** *V* voltage applied to the dielectric, *Vc* voltage across the cavity, *Vci* inception voltage of the discharge, *Vce* discharge extinction voltage.

measured outside of the specimen. The phenomenon is then traced to two successive stages: - the first defines the external current measured with the induced charge to the electrodes. - the second one correlates the induced charge to the electrodes with the charge inside the cavity. A field approach model that best represents the phenomenon is conceived by Pedersen [18]. Since the charge involved by the discharge inside the cavity, named true charge, cannot be measured, an auxiliary current, named apparent charge, is measured outside of the specimen. The phenomenon is then traced to two successive stages:


the realization phase of new devices but also for evaluate reliability in work. The discharge manifests itself through charge displacement, acoustic waves, lights and electromagnetic field. The wired detection of charge displacement was the most frequently measurement method chosen before the introduction of antenna probes. The basic discharge detection circuit, in accordance with the International Electrotechnical Commission (IEC) Standard, is shown in fig. 3 [19]. The coupling capacitor provides a low impedance path for high-frequency discharge currents which cause voltage pulses over the measuring impedance. The pulses are then amplified and stored in a computer for further processing

#### **2.2. Partial discharge characterization**

Figure 1. Three capacity of the model of Whitehead.

**Figure 1.** Three capacity of the model of Whitehead.

employing current and voltage measures.

**2.2. Partial discharge characterization** 

The phenomenon is then traced to two successive stages:

voltage.

discharge extinction voltage.

cavity.

voltage measures.

180 Progress in Compact Antennas

Fig. 2 presents the *Vc* voltage applying on the c capacitor. It should be noticed that *Vc* is a fraction of the total voltage V. When *Vc* exceeds the so called inception voltage *Vci*, the discharge starts. After the discharge, *Vc* reaches the discharge extinction voltage value, *Vce*. The global behaviour of the system is affected and the PD activity can be monitored by

Fig. 2 presents the *Vc* voltage applying on the c capacitor. It should be noticed that *Vc* is a fraction of the total voltage V. When *Vc* exceeds the so called inception voltage *Vci*, the discharge starts. After the discharge, *Vc* reaches the discharge extinction voltage value, *Vce*. The global behaviour of the system is affected and the PD activity can be monitored by employing current and

*a* 

*b s <sup>D</sup> <sup>d</sup> s D*;

*c* 0 *s d* ;

*S s <sup>D</sup>* 

*S D*;

Figure 2. *V* voltage applied to the dielectric, *Vc* voltage across the cavity, *Vci* inception voltage of the discharge, *Vce* discharge extinction

A field approach model that best represents the phenomenon is conceived by Pedersen [18]. Since the charge involved by the discharge inside the cavity, named true charge, cannot be measured, an auxiliary current, named apparent charge, is

A partial discharge is a stochastic phenomenon that involves many chemical and physical aspects of materials proprieties. The PD detection over time has become an indispensable method for the evaluation of insulation systems in the realization phase of new devices but also for evaluate reliability in work. The discharge manifests itself through charge displacement, acoustic waves, lights and electromagnetic field. The wired detection of charge displacement was the most frequently measurement method chosen before the introduction of antenna probes. The basic discharge detection circuit, in accordance with the International Electrotechnical Commission (IEC) Standard, is shown in fig. 3 [19]. The coupling capacitor provides a low impedance path for high-frequency discharge currents which cause voltage pulses over the measuring impedance. The pulses are then amplified and stored in a computer for further processing

measured outside of the specimen. The phenomenon is then traced to two successive stages:

**Figure 2.** *V* voltage applied to the dielectric, *Vc* voltage across the cavity, *Vci* inception voltage of the discharge, *Vce*

A field approach model that best represents the phenomenon is conceived by Pedersen [18]. Since the charge involved by the discharge inside the cavity, named true charge, cannot be measured, an auxiliary current, named apparent charge, is measured outside of the specimen.

**•** the first defines the external current measured with the induced charge to the electrodes.

**•** the second one correlates the induced charge to the electrodes with the charge inside the



A partial discharge is a stochastic phenomenon that involves many chemical and physical aspects of materials proprieties. The PD detection over time has become an indispensable method for the evaluation of insulation systems in the realization phase of new devices but also for evaluate reliability in work. The discharge manifests itself through charge displace‐ ment, acoustic waves, lights and electromagnetic field. The wired detection of charge dis‐ placement was the most frequently measurement method chosen before the introduction of antenna probes. The basic discharge detection circuit, in accordance with the International Electrotechnical Commission (IEC) Standard, is shown in fig. 3 [19]. The coupling capacitor provides a low impedance path for high-frequency discharge currents which cause voltage pulses over the measuring impedance. The pulses are then amplified and stored in a computer for further processing and displayed in different patterns as 3D patterns (phase-amplitudenumber of discharge) and 2D patterns (amplitude-number, phase-number, phase-amplitude).

In the same IEC Standard the main definitions related to PD are defined as follows:

*Partial discharge (PD)*: Localized electrical discharge that only partially bridges the insulation between conductors and which can or cannot occur adjacent to a conductor.

*Partial discharge pulse (PD pulse):* Current or voltage pulse that results from a partial discharge occurring within the object under test. The pulse is measured using suitable detector circuits, which have been introduced into the test circuit for the purpose of the test.

*Apparent charge q:* Of a PD pulse is that charge which, if injected within a very short time between the terminals of the test object in a specified test circuit, would give the same reading on the measuring instrument ad the PD current pulse itself. The apparent charge is usually expressed in picocoulombs (pC).

*Pulse repetition rate n:* Ratio between the total number of PD Pulses recorded in a selected time interval and the duration of this time interval.

*Phase angle ϕ<sup>i</sup> and time ti of occurrence of a PD pulse:* is *ϕ*<sup>i</sup> =360(ti /T) where ti is the time measured between the preceding positive going transition of the test voltage through zero and the partial discharge pulse and T is the period of the test *voltage.*

The pulses generated by PDs are detected as second order pulses; in Fig.4 the pulse was recorded by means of an oscilloscope and shows a signal proportional to the discharge current. It was detected by means of a 50 Ω resistance in series to the circuit where the defect was present. This kind of pulses are very fast (~50ns) and have a small amplitude (~5mV). As we will see in the following, not only the amplitude of the pulse but also its duration is a funda‐ mental parameter. The most important recent studies that have allowed PD measurements to become a reliable diagnostic system are based on the shape of the pulses, and in particular on the ratio between duration and frequency content. The acquisition techniques (bandwidth of the recording system) and the kind of sensor adopted (HFCT, capacitive couplers) become fundamental components because they affect the shape of the impulses.

The pulse shape of singular discharge phenomena present differences in term of time duration and frequency content. These differences are exhibited in Fig. 5 where an example of pulse

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

183

shape and frequency spectrum for each discharge family is showed.

**Figure 5.** Time duration and frequency spectrum of surface, corona and internal discharges.

patterns for each discharge family as shown in Fig.7, 8 and 9.

scale represents the intensity of the discharges.

Moreover, the acquisition of the pulse repetition on a time interval ΔT is displayed in a 3D histogram generally known as PD Pattern (phase-amplitude-number of discharge) and 2D PD patterns (amplitude-number, phase-number, phase-amplitude,) as displayed in Fig.6. Fur‐ thermore, the pattern allows the understanding of the kind of the defect.We obtain different

**Figure 6.** PD Pattern with 2D amplitude-number (right) and phase-number (below). In the 3D PD Pattern the colour

**Figure 3.** Discharge detection circuit.

**Figure 4.** Generic second order PD Pulse.

In the PD detection theory, three families of PD phenomena have been taken into account: namely internal, surface and corona. Internal discharges occur inside a cavity embedded in a dielectric. These are cavities within resin components, extruded plastic cables, joint, stator insulation and external particles such as dust or textile fibres embedded in insulation. Surface discharge may occur along dielectric surfaces with a high tangential field stress. Typical examples of high voltage components that present surface discharges are bushing, twisted pairs, ends of cables, the overhang of generator windings. Corona discharges occur in a needle plane configuration at sharp points at high voltage or ground potential [20].

The pulse shape of singular discharge phenomena present differences in term of time duration and frequency content. These differences are exhibited in Fig. 5 where an example of pulse shape and frequency spectrum for each discharge family is showed.

**Figure 5.** Time duration and frequency spectrum of surface, corona and internal discharges.

**Figure 4.** Generic second order PD Pulse.

**Figure 3.** Discharge detection circuit.

182 Progress in Compact Antennas

In the PD detection theory, three families of PD phenomena have been taken into account: namely internal, surface and corona. Internal discharges occur inside a cavity embedded in a dielectric. These are cavities within resin components, extruded plastic cables, joint, stator insulation and external particles such as dust or textile fibres embedded in insulation. Surface discharge may occur along dielectric surfaces with a high tangential field stress. Typical examples of high voltage components that present surface discharges are bushing, twisted pairs, ends of cables, the overhang of generator windings. Corona discharges occur in a needle

plane configuration at sharp points at high voltage or ground potential [20].

Moreover, the acquisition of the pulse repetition on a time interval ΔT is displayed in a 3D histogram generally known as PD Pattern (phase-amplitude-number of discharge) and 2D PD patterns (amplitude-number, phase-number, phase-amplitude,) as displayed in Fig.6. Fur‐ thermore, the pattern allows the understanding of the kind of the defect.We obtain different patterns for each discharge family as shown in Fig.7, 8 and 9.

**Figure 6.** PD Pattern with 2D amplitude-number (right) and phase-number (below). In the 3D PD Pattern the colour scale represents the intensity of the discharges.

Fig.7 shows a typical internal void rabbit-like PD pattern. Internal discharge are characterized by a quite symmetry for positive and negative discharges due to electron emission from the same materials (dielectric). In particular the showed pattern is referred to a spherical air void embedded in a dielectric subjected to a uniform electric field.

Fig.9 shows a typical corona PD pattern. Two characteristic differentiates it from the previous discharges. The first is the presence of discharge only in half period of AC source and the second is that the discharge occurs in correspondence of the voltage pick of AC source, in the other cases, differently, the discharges occur in presence of the maximum gradient of the

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

185

In many real cases PD patterns present very different shapes and sometimes with more than only one defect at the same time. For this reason pattern recognition techniques for the discrimination and classification of discharges represent a fundamental instrument for the reliability's evaluation of high voltage electrical apparatus. Many different methods have been proposed in literature to automate the identification of defects through PD pattern recognition. They differ either in presentation of PD data (PD pulse shape data, statistical parameters of these data and pulse density distribution) or in the classification methods (statistical algo‐ rithms or neural network (NN)) [20-22]. Recently the use of wavelet analysis and advanced models has provided interesting results in PD denoising processes, pattern recognition and

One limit of the classic PD measurement s is due to the wired circuit, signals are dumped or affected by the particular circuit path encountered. In order to overcome these limits wireless

The Method of Moments (MoM) in two dimensions has been successfully applied to describe the behaviour of three-dimensional objects that have a highly shielding external structure, for which the hypothesis of perfect electric conductors can be used [11]. Other applications of MoM involving three-dimensional objects composed of voxels or thin wire elements can be

voltage.

**Figure 9.** PD Pattern of corona discharges.

classification of multi-source PD patterns [23,26].

measurement have been proposed.

**3. Method of moments**

**Figure 7.** PD Pattern of internal discharges.

Fig.8 shows a typical surface PD pattern, where a substantial dissymmetry between positive and negative discharges can be noted. The inception phase angle, Δ*ϕ*, defined as the phase interval between the first and the last discharge (positive or negative) is smaller than the one observed for internal discharges.

**Figure 8.** PD Pattern of surface discharges.

Fig.9 shows a typical corona PD pattern. Two characteristic differentiates it from the previous discharges. The first is the presence of discharge only in half period of AC source and the second is that the discharge occurs in correspondence of the voltage pick of AC source, in the other cases, differently, the discharges occur in presence of the maximum gradient of the voltage.

**Figure 9.** PD Pattern of corona discharges.

Fig.7 shows a typical internal void rabbit-like PD pattern. Internal discharge are characterized by a quite symmetry for positive and negative discharges due to electron emission from the same materials (dielectric). In particular the showed pattern is referred to a spherical air void

Fig.8 shows a typical surface PD pattern, where a substantial dissymmetry between positive and negative discharges can be noted. The inception phase angle, Δ*ϕ*, defined as the phase interval between the first and the last discharge (positive or negative) is smaller than the one

embedded in a dielectric subjected to a uniform electric field.

**Figure 7.** PD Pattern of internal discharges.

184 Progress in Compact Antennas

observed for internal discharges.

**Figure 8.** PD Pattern of surface discharges.

In many real cases PD patterns present very different shapes and sometimes with more than only one defect at the same time. For this reason pattern recognition techniques for the discrimination and classification of discharges represent a fundamental instrument for the reliability's evaluation of high voltage electrical apparatus. Many different methods have been proposed in literature to automate the identification of defects through PD pattern recognition. They differ either in presentation of PD data (PD pulse shape data, statistical parameters of these data and pulse density distribution) or in the classification methods (statistical algo‐ rithms or neural network (NN)) [20-22]. Recently the use of wavelet analysis and advanced models has provided interesting results in PD denoising processes, pattern recognition and classification of multi-source PD patterns [23,26].

One limit of the classic PD measurement s is due to the wired circuit, signals are dumped or affected by the particular circuit path encountered. In order to overcome these limits wireless measurement have been proposed.

#### **3. Method of moments**

The Method of Moments (MoM) in two dimensions has been successfully applied to describe the behaviour of three-dimensional objects that have a highly shielding external structure, for which the hypothesis of perfect electric conductors can be used [11]. Other applications of MoM involving three-dimensional objects composed of voxels or thin wire elements can be found in [12]. We note *S* the surface of a metal object considered a perfect conductor (PEC) and *n* the normal to surface *S*. By considering an Electromagnetic (EM) stress, *Ei* denotes the incident electric field, which can be found in the absence of the object having *S* surface. The presence of an EM field generates surface current *JS* on *S*. The aim of the EM analysis is to determine the current in *S*. Since no analytical solution can be found for complex geometries only a numerical one can be obtained approximating the current *JS* on *S* into *Ji* currents on *<sup>S</sup> <sup>N</sup>* <sup>=</sup>⋃<sup>1</sup> *<sup>N</sup> Si* , (see Fig. 10). Triangular elementary surfaces are used to describe the whole *SN*, by introducing edge limiting each triangular faces. Vertexes are the points shared by different triangles. Each triangle, illuminated by external EM field, has its own current, which can conceptually be studied by examining the average current that passes through the three edges. The aim of the MoM is to evaluate each current *Ji* passing through an edge.

The current on the whole *S* can be approximated by:

$$
\vec{J}\_S = \sum\_{i=1}^N f\_i \vec{f}\_i \text{ (}\vec{r}\text{)}\tag{1}
$$


*A* → (*r* <sup>→</sup> )= *<sup>μ</sup>* <sup>4</sup>*<sup>π</sup> ∫ S J* → *S* (*r* → ) *e* - *jkR*

*φ*(*r* <sup>→</sup> )= <sup>1</sup> <sup>4</sup>*πε ∫ <sup>S</sup> σ*(*r* → ) *e* - *jkR <sup>R</sup> dS* '

*σ*, which can be related to *JS* by taking into account the following:

the one used in equation 1 is required, considering only the unknowns *Ji*

By considering an edge *i* shared by two triangles called *Ti*

distance from the free vertex of *Ti*

O

Ji

(M)

+

+ for *Ti* -

+ i r M

 and C- .

> T+ i

a

A

where *A*

with *k* =*ω εμ*.

Points (M) belonging to *Ti*

+

crosses the *i* side, leaving *Ti*

**O** or through ρ<sup>i</sup>

to the centroids C+

*tan* = - ( *jωA*

Such potentials are depending on the local surface current *JS* and on the local charge density

The set of equations described here takes the name of Electric Field Integral Equation (EFIE).

Equations 4 and 5 have no analytical solution for a complex *S*, so an approximation similar to

can be detected by means of the distance *r*

+

the current in each point M of the triangles, by using the basis functions(see Fig. 11a).

choice of the+and-polarities is due to the implementation of a surface current density *Ji*

Ti

B

**Figure 11.** Triangles sharing an internal edge. Focus on the current through the edge (a), focus on centroids.

A

O

C+ C-

M

ri c-

c+i

ri c+

<sup>→</sup> and *φ* are the vector and the scalar potential respectively:

<sup>→</sup> - <sup>∇</sup>*φ*)*tan*, (3)

*<sup>R</sup> dS*' (4)

∇ ∙ *JS* = - *jωσ* (6)

<sup>+</sup> and *Ti* -

. In Fig.11 the behaviour of the current in M is associated

B

ci

2

. Similar remarks can be applied for *Ti*

( c+i + ci ) <sup>1</sup>

b

, (5)

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

187

set in equation 1.

, it is possible to describe

<sup>→</sup> from an absolute center


, which

where N is the number of internal sides, and fi (r <sup>→</sup> ) are a suitable functions, called basis functions.

**Figure 10.** Loop antenna described with triangular surfaces

The basis functions to be used must meet two basic requirements. Indeed, they have to:


Integro-differential equation can be derived by imposing equality between the incident field and scattered field on the boundary condition on *S*,

$$\stackrel{\wedge}{m} \times \left(E\_i + E\_s\right) = 0\tag{2}$$

getting:

Building Partial Discharge Signal Wireless Probes http://dx.doi.org/10.5772/58840 187

$$-E\_{inc}^{\text{tan}} = -\left(j\omega\vec{A} - \nabla\,\varphi\right)^{\text{tan}},\tag{3}$$

where *A* <sup>→</sup> and *φ* are the vector and the scalar potential respectively:

$$
\vec{A}(\vec{r}) = \frac{\mu}{4\pi} \mathbb{f}\_S \vec{f}\_S(\vec{r}) \frac{e^{-j k \vec{\ell}}}{R} dS' \tag{4}
$$

$$\varphi(\vec{r}) = \frac{1}{4\pi\epsilon} \mathbb{J}\_S \sigma(\vec{r}) \frac{e^{-\vec{\beta}\cdot\vec{R}}}{R} dS \,\hat{\phantom{r}}\,\tag{5}$$

with *k* =*ω εμ*.

found in [12]. We note *S* the surface of a metal object considered a perfect conductor (PEC) and

incident electric field, which can be found in the absence of the object having *S* surface. The presence of an EM field generates surface current *JS* on *S*. The aim of the EM analysis is to determine the current in *S*. Since no analytical solution can be found for complex geometries

introducing edge limiting each triangular faces. Vertexes are the points shared by different triangles. Each triangle, illuminated by external EM field, has its own current, which can conceptually be studied by examining the average current that passes through the three edges.

(r

The basis functions to be used must meet two basic requirements. Indeed, they have to:

**•** Allow modeling with flat triangles. Further useful developments can be found in [27].

*n*

Integro-differential equation can be derived by imposing equality between the incident field

, (see Fig. 10). Triangular elementary surfaces are used to describe the whole *SN*, by

passing through an edge.

edge face

^ ×(*Ei* <sup>+</sup> *Es*) =0 (2)

<sup>→</sup> ), (1)

<sup>→</sup> ) are a suitable functions, called basis functions.

vertex

denotes the

currents on

*n* the normal to surface *S*. By considering an Electromagnetic (EM) stress, *Ei*

only a numerical one can be obtained approximating the current *JS* on *S* into *Ji*

*J* → *<sup>S</sup>* = ∑ *i*=1 *N Ji f* → *i* (*r*

*<sup>S</sup> <sup>N</sup>* <sup>=</sup>⋃<sup>1</sup>

getting:

*<sup>N</sup> Si*

186 Progress in Compact Antennas

The aim of the MoM is to evaluate each current *Ji*

where N is the number of internal sides, and fi

aperture

**Figure 10.** Loop antenna described with triangular surfaces

boundary edge

**•** be suitable to be used with an electromagnetic equation.

and scattered field on the boundary condition on *S*,

The current on the whole *S* can be approximated by:

Such potentials are depending on the local surface current *JS* and on the local charge density *σ*, which can be related to *JS* by taking into account the following:

$$
\nabla \bullet \mathbf{J}\_S = -\mathbf{j}\omega\sigma \tag{6}
$$

The set of equations described here takes the name of Electric Field Integral Equation (EFIE).

Equations 4 and 5 have no analytical solution for a complex *S*, so an approximation similar to the one used in equation 1 is required, considering only the unknowns *Ji* set in equation 1.

By considering an edge *i* shared by two triangles called *Ti* <sup>+</sup> and *Ti* - , it is possible to describe the current in each point M of the triangles, by using the basis functions(see Fig. 11a).

Points (M) belonging to *Ti* + can be detected by means of the distance *r* <sup>→</sup> from an absolute center **O** or through ρ<sup>i</sup> + distance from the free vertex of *Ti* + . Similar remarks can be applied for *Ti* - . The choice of the+and-polarities is due to the implementation of a surface current density *Ji* , which crosses the *i* side, leaving *Ti* + for *Ti* - . In Fig.11 the behaviour of the current in M is associated to the centroids C+ and C- .

**Figure 11.** Triangles sharing an internal edge. Focus on the current through the edge (a), focus on centroids.

The purpose of the basis function is to simplify the use of equations 4 and 5, projected in the numerical approximation. The integrals on *S* are solved for each triangle by taking into account three currents leaving the three edges. The basis function vector is defined by:

$$
\vec{\hat{f}}\_i(\vec{r}) = \begin{vmatrix}
\frac{l\_i}{2A\_i^+} \vec{\hat{\rho}\_i} & \vec{\text{r}} \, \text{in} \, T\_i^+ \\
\frac{l\_i}{2A\_i^-} \vec{\hat{\rho}\_i} & \vec{\text{r}} \, \text{in} \, T\_i^- \\
0 & \text{elsewhere}
\end{vmatrix} \tag{7}
$$

with *f*

*Einc <sup>p</sup>*,*tan*, *f* →

is possible to write equation 11 as follows:

*<sup>p</sup>* <sup>=</sup> *<sup>j</sup><sup>ω</sup> <sup>μ</sup>* <sup>4</sup>*<sup>π</sup> ∫ S J* → *S* (*r* → ) *e* - *jkR <sup>R</sup> dS*', *f* → *p tan* - <sup>1</sup> <sup>4</sup>*πε ∫ <sup>S</sup> σ*(*r* → ) *e* - *jkR*

on edge of length *lm* of *p* and edge of length *li*

a

and the one inside surface integral *∫*

+ and *Tm* -

*<sup>m</sup> tan* <sup>+</sup> *<sup>j</sup>ω<sup>A</sup>*

→ (*r* → *m c*- ), *f* → *<sup>m</sup> tan* <sup>=</sup>*<sup>∫</sup> T <sup>m</sup>* <sup>+</sup> *jωA* → (*r* → *m c*+ ) *f* → *<sup>m</sup>dS* ' +*∫ T <sup>m</sup>* - *jωA* → (*r* → *m c*- ) *f* → *<sup>m</sup>dS* ' =

1 <sup>2</sup> *jωA* → (*r* → *m c*+ )*l <sup>m</sup>ρ<sup>m</sup> <sup>c</sup>*<sup>+</sup> + 1 <sup>2</sup> *jωA* → (*r* → *m c*- )*l <sup>m</sup>ρ<sup>m</sup> c*-

*φ*(*r* → *m c*+ ), ∇ ∙ *f* → *m*

*A* → (*r* → *m c*± ) <sup>=</sup> *<sup>μ</sup>* <sup>4</sup>*<sup>π</sup> ∫ Ti* <sup>±</sup>*Ji f* → *i e* - *jk Rm* ± *Rm*

edge *m,* shared by *Tm*

for all triangle, we can write:

*Einc <sup>m</sup>* , *f* → *<sup>m</sup>* =*∫ T <sup>m</sup>* + +*T <sup>m</sup>* -*Einc m f* → *<sup>m</sup>dS* ' =*E inc m*, *c*+ *∫ T <sup>m</sup>* + *f* → *<sup>m</sup>dS* '

*Einc m*, *c*+ ( *l m* 2*Am* <sup>+</sup> *ρ<sup>m</sup> c*+ )*Am* <sup>+</sup> <sup>+</sup> *Einc m*, *c*- ( *l m* 2*Am* - *ρ<sup>m</sup> c*- )*Am* - <sup>=</sup> <sup>1</sup> <sup>2</sup> *Einc m*, *c*+ *l <sup>m</sup>ρ<sup>m</sup> <sup>c</sup>*<sup>+</sup> + 1 <sup>2</sup> *Einc m*, *cl <sup>m</sup>ρ<sup>m</sup> c*- ,

*∫ T <sup>m</sup>* <sup>+</sup>*φ*(*r* → *m c*+ ), ∇ ∙ *f* → *<sup>m</sup>dS* ' +*∫ T <sup>m</sup>* -*φ*(*r* → *m c*- ), ∇ ∙ *f* → *<sup>m</sup>dS* '

*jωA* → (*r* → *m c*+ ), *f* →

on the edge *i*:

r r'

O

**Figure 12.** Interactions between currents of different triangles

<sup>q</sup> <sup>p</sup>

Jp

<sup>→</sup> equals to zero on the boundary of both *S* and ∂*S*. Finally, by considering a triangle *p* it

To study the effects of the triangle *p* on the current due to the current of a triangle *q*(seeFig. 12.a), all the edge currents have to be considered. In Fig.12.b, the interaction between currents

> T+ m

> > T-

In equation 13 the inner product has to be considered twice: the formal inner product with *f*

*tan* + *φ*(*r* → *m c*- ), ∇ ∙ *f* → *m tan* =

The second inner product is inside the vector and scalar potentials, that are due to the current

<sup>±</sup> *dS*'= *<sup>μ</sup>*

<sup>8</sup>*<sup>π</sup> Ji l i ρi c*± *e* - *jk Rm* ± *Rm*

Jq

of *q*,is shown.

<sup>r</sup> r' <sup>l</sup><sup>i</sup> lm

<sup>m</sup> b

*<sup>S</sup> dS'*. By considering initially the first inner product, for

, and also considering value on the centroid as a constant value

+ *Einc m*, *c*- *∫ T <sup>m</sup>* - *f* → *<sup>m</sup>dS* ' =

=*φ*(*r* → *m c*+ )*l <sup>m</sup>* - *φ*(*r* → *m c*- )*l*

O

*<sup>R</sup> dS*', ∇ ∙ *f*

T+ i → *p tan*

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

Ii

Ti . (13)

189

→ *p*

(14)

(15)

*<sup>m</sup>*. (16)

<sup>±</sup> , (17)

Where *li* is the length of the edge taken into consideration, *Ai* + and *Ai* are the areas of the triangles. Also the following two features are required:

$$\nabla \bullet \vec{\hat{f}}\_i(\vec{r}) = \begin{cases} \frac{l\_i}{A\_i^+} & \min T\_i^+ \\ \frac{l\_i}{A\_i^-} & \min T\_i \\ 0 & \text{elsewhere} \end{cases} \tag{8}$$

*f* → *i* (*r* <sup>→</sup> )=0 <sup>→</sup> on the boundary of *S*.

It is not possible to solve equation 3 on the whole S, but it can be solved on *<sup>S</sup> <sup>N</sup>* <sup>=</sup>⋃<sup>1</sup> *<sup>N</sup> Si* , by employing the numerical projection on a N-dimension system. Similarly to the Fourier's sampling a particular function has to be used to project in the N-dimension space, instead of the delta of dirac function here are used the basis functions:

$$
\{\vec{X}, \,\vec{f}\} = \mathfrak{f}\_{\mathbb{S}} \vec{X} \bullet \vec{f} \, dS,\tag{9}
$$

which in the *N-space* becomes:

$$\{\vec{\mathbf{X}},\,\,\vec{f}\} \cong \sum\_{i=1}^{N} \vec{X}\_{i\,} \bullet \,\vec{f}\_{\,i} A\_{i\prime} \tag{10}$$

with *X* <sup>→</sup> a vector function. In such a way equation 3 becomes:

$$
\{\vec{E}\_{inc}^{\tan}, \vec{f}\} = \{j\omega \vec{A}, \,\vec{f}\}^{\tan} + \{\nabla \,\varphi, \,\vec{f}\}^{\tan}.\tag{11}
$$

In equation 11 the gradient is applied on the scalar potential, but considering the following identity (similar to the well known Stokes theorem), it is possible to apply it on the basis function:

$$\{\nabla\,\boldsymbol{\varphi},\,\vec{f}\} \mathbf{=} \{\nabla\,\boldsymbol{\varphi}\,\boldsymbol{\sigma}\,\boldsymbol{\tau}\} \mathbf{dS} = \mathfrak{f}\_{\sf S} \nabla\ \bullet \{\boldsymbol{\varphi}\} \mathbf{dS} \ \mathbf{-} \mathbf{\{\!}\boldsymbol{\varphi}\boldsymbol{\nabla}\bullet\ \vec{f}\} \mathbf{dS} = \mathfrak{f}\_{\sf S} \boldsymbol{\varphi}\boldsymbol{f} \mathbf{d} \mathbf{S} \ \mathbf{-} \mathbf{\{\!}\boldsymbol{\varphi}\boldsymbol{\nabla}\bullet\ \vec{f}\} \mathbf{dS} = \mathbf{-}\{\!\boldsymbol{\varphi}\,,\,\mathbf{\nabla}\bullet\ \vec{f}\} \mathbf{\!} \tag{12}$$

with *f* <sup>→</sup> equals to zero on the boundary of both *S* and ∂*S*. Finally, by considering a triangle *p* it is possible to write equation 11 as follows:

$$\left\{\mathbb{E}\_{\text{inc}}^{p,\text{tan}}, \begin{array}{c} \vec{f} \end{array} \right\} = \left\{j\omega\frac{\mu}{4\pi}\mathbb{f}\_{\text{S}}\vec{f}\right\}\_{\text{S}}(\vec{r})\xleftarrow{\varepsilon\cdot\mathbb{R}^{3}} dS', \left\{\vec{f} \right\}\_{\text{p}}^{\text{tan}}\cdot\left\{\frac{1}{4\pi\varepsilon}\mathbb{f}\_{\text{S}}\sigma(\vec{r})\frac{\varepsilon\cdot\mathbb{R}^{3}}{\mathbb{R}}dS', \left.\nabla\bullet\right\cdot\vec{f}\right\}\_{\text{p}}^{\text{tan}}.\tag{13}$$

To study the effects of the triangle *p* on the current due to the current of a triangle *q*(seeFig. 12.a), all the edge currents have to be considered. In Fig.12.b, the interaction between currents on edge of length *lm* of *p* and edge of length *li* of *q*,is shown.

**Figure 12.** Interactions between currents of different triangles

The purpose of the basis function is to simplify the use of equations 4 and 5, projected in the numerical approximation. The integrals on *S* are solved for each triangle by taking into account

three currents leaving the three edges. The basis function vector is defined by:

*l i* 2*Ai* - *ρi* <sup>→</sup> *rinTi* -

is the length of the edge taken into consideration, *Ai*

0 *elsewhere*

<sup>+</sup> *rinTi* +

0 *elsewhere*

employing the numerical projection on a N-dimension system. Similarly to the Fourier's sampling a particular function has to be used to project in the N-dimension space, instead of


It is not possible to solve equation 3 on the whole S, but it can be solved on *<sup>S</sup> <sup>N</sup>* <sup>=</sup>⋃<sup>1</sup>

, (7)

, (8)

<sup>→</sup> *dS*, (9)

, (10)

. (11)

<sup>→</sup> *dS* <sup>=</sup> - *<sup>φ</sup>*, <sup>∇</sup> <sup>∙</sup> *<sup>f</sup>*

<sup>→</sup> , (12)

are the areas of the

*<sup>N</sup> Si* , by

+ and *Ai* -

*f* → *i* (*r* → )={ *l i* 2*Ai* + *ρi* <sup>→</sup> *rinTi* +

∇ ∙ *f* → *i* (*r* → )={ *l i Ai*

the delta of dirac function here are used the basis functions:

*X* <sup>→</sup>, *<sup>f</sup>* <sup>→</sup> <sup>=</sup>*<sup>∫</sup> <sup>S</sup> X* <sup>→</sup> <sup>∙</sup> *<sup>f</sup>*

*X* <sup>→</sup>, *<sup>f</sup>* <sup>→</sup> <sup>≅</sup> <sup>∑</sup> *i*=1 *N X* → *<sup>i</sup>* ∙ *f* → *iAi*

<sup>→</sup> a vector function. In such a way equation 3 becomes:

<sup>→</sup> )*dS* - *<sup>∫</sup>*

<sup>→</sup> <sup>=</sup> *<sup>j</sup>ω<sup>A</sup>*

<sup>→</sup> , *<sup>f</sup>* <sup>→</sup> *tan*

*<sup>S</sup>φ*∇ ∙ *f*

+ ∇*φ*, *f*

In equation 11 the gradient is applied on the scalar potential, but considering the following identity (similar to the well known Stokes theorem), it is possible to apply it on the basis

<sup>→</sup> *dS* <sup>=</sup>*<sup>∫</sup>*

<sup>∂</sup>*Sφ f* <sup>→</sup> *dS* - *<sup>∫</sup>*

*<sup>S</sup>φ*∇ ∙ *f*

<sup>→</sup> *tan*

*Einc tan*, *f*

*<sup>S</sup>*∇ ∙(*φ f*

triangles. Also the following two features are required:

<sup>→</sup> on the boundary of *S*.

which in the *N-space* becomes:

*<sup>S</sup>*∇*φ* ∙ *f*

<sup>→</sup> *dS* <sup>=</sup>*<sup>∫</sup>*

Where *li*

188 Progress in Compact Antennas

*f* → *i* (*r* <sup>→</sup> )=0

with *X*

function:

∇*φ*, *f* <sup>→</sup> <sup>=</sup>*<sup>∫</sup>* In equation 13 the inner product has to be considered twice: the formal inner product with *f* → *p* and the one inside surface integral *∫ <sup>S</sup> dS'*. By considering initially the first inner product, for edge *m,* shared by *Tm* + and *Tm* - , and also considering value on the centroid as a constant value for all triangle, we can write:

$$\begin{split} \{\boldsymbol{E}\_{\rm inc}^{m}, \stackrel{\rightarrow}{\boldsymbol{f}}\_{m}\} &= \mathbb{J}\_{\Gamma\_{m}^{+}+\Gamma\_{m}^{-}} \boldsymbol{E}\_{\rm inc}^{m} \stackrel{\rightarrow}{\boldsymbol{f}}\_{m} d\boldsymbol{S} \stackrel{\rightarrow}{\boldsymbol{\ast}} = \boldsymbol{E}\_{\rm inc}^{m,c^{\circledast}} \boldsymbol{J}\_{\Gamma\_{m}^{-}} \stackrel{\rightarrow}{\boldsymbol{f}}\_{m} d\boldsymbol{S} \stackrel{\rightarrow}{\boldsymbol{\ast}} + \boldsymbol{E}\_{\rm inc}^{m,c^{\circledast}} \boldsymbol{J}\_{\Gamma\_{m}^{-}} \stackrel{\rightarrow}{\boldsymbol{f}}\_{m} d\boldsymbol{S} \stackrel{\rightarrow}{\boldsymbol{\ast}} \\ &= \boldsymbol{E}\_{\rm inc}^{m,c^{\circledast}} \left( \frac{l\_{m}}{2A\_{n}^{+}} \boldsymbol{\rho}\_{m}^{c^{\circledast}} \right) \boldsymbol{A}\_{m}^{+} + \boldsymbol{E}\_{\rm inc}^{m,c^{\circledast}} \left( \frac{l\_{m}}{2A\_{n}^{-}} \boldsymbol{\rho}\_{m}^{c^{\circledast}} \right) \boldsymbol{A}\_{m}^{-} = \frac{1}{2} \boldsymbol{E}\_{\rm inc}^{m,c^{\circledast}} l\_{m} \boldsymbol{\rho}\_{m}^{c\_{\circledast}} + \frac{1}{2} \boldsymbol{E}\_{\rm inc}^{m,c^{\circledast}} l\_{m} \boldsymbol{\rho}\_{m}^{c\_{\circledast}}. \end{split} \tag{14}$$

$$\begin{aligned} \left\{ \dot{j}\omega \overrightarrow{A} \Big| \left( \dot{r}\_{m}^{\circ \circ} \right)\_{\prime} \overrightarrow{f}\_{m} \right\}^{\text{tran}} + \left\{ \dot{j}\omega \overrightarrow{A} \Big| \left( \dot{r}\_{m}^{\circ \circ} \right)\_{\prime} \overrightarrow{f}\_{m} \right\}^{\text{tran}} &= \mathbf{f}\_{\Gamma\_{m}^{+}} \dot{j}\omega \overrightarrow{A} \Big( \dot{r}\_{m}^{\circ \circ} \right) \overrightarrow{f}\_{m} dS \overset{\circ}{\to} \mathbf{j}\_{\Gamma\_{m}^{-}} \dot{j}\omega \overrightarrow{A} \Big( \dot{r}\_{m}^{\circ \circ} \right) \overrightarrow{f}\_{m} dS \overset{\circ}{=} \\ &\quad \frac{1}{2} \left\{ \dot{j}\omega \overrightarrow{A} \Big| \left( \dot{r}\_{m}^{\circ \circ} \right) l\_{m} \rho\_{m}^{\circ \circ \circ} + \frac{1}{2} \left\{ \dot{j}\omega \overrightarrow{A} \Big| \left( \dot{r}\_{m}^{\circ \circ} \right) l\_{m} \rho\_{m}^{\circ \circ} \right. \end{aligned} \tag{15}$$

$$\begin{aligned} \{\rho(\stackrel{\cdot}{r}\_{m}^{c\*}), \,\,\nabla\,\bullet\,\vec{\boldsymbol{f}}\_{m}\}\_{m}^{\text{tan}} &+ \{\rho(\stackrel{\cdot}{r}\_{m}^{c\*}), \,\,\nabla\,\bullet\,\vec{\boldsymbol{f}}\_{m}\}\_{m}^{\text{tan}} =\\ \{\!\_{\boldsymbol{\Gamma}\_{m}}\!\!\rho(\stackrel{\cdot}{r}\_{m}^{c\*}), \,\,\nabla\,\bullet\,\vec{\boldsymbol{f}}\_{m}\}\_{m}^{\text{cis}} &\stackrel{\cdot}{\boldsymbol{\Gamma}}\_{\text{T}\_{m}^{c}}\!\!\rho(\stackrel{\cdot}{r}\_{m}^{c\*}), \,\,\nabla\,\bullet\,\vec{\boldsymbol{f}}\_{m}\stackrel{\cdot}{\boldsymbol{\Gamma}}\_{m}^{\text{cis}} &= \rho(\stackrel{\cdot}{r}\_{m}^{c\*})l\_{m} - \rho(\stackrel{\cdot}{r}\_{m}^{c\*})l\_{m}. \end{aligned} \tag{16}$$

The second inner product is inside the vector and scalar potentials, that are due to the current on the edge *i*:

$$\overrightarrow{A}\left(\overrightarrow{r'}\_{m}^{\rm ex}\right) = \frac{\mu}{4\pi} \mathbf{f}\_{T\_i} \mathbf{J}\_i \overrightarrow{f}\_i \overleftarrow{f}\_i \overleftarrow{f}\_i \overleftarrow{\frac{e^{\cdot \cdot k' R\_m^{\rm ex}}}{R\_m^{\rm ex}}} dS' = \frac{\mu}{8\pi} \mathbf{J}\_i \mathbf{l}\_i \mathbf{l}\_i \rho\_i^{\rm ex} \frac{e^{\cdot \cdot k' R\_m^{\rm ex}}}{R\_m^{\rm ex}} \tag{17}$$

where *Rm* ± is the distance of the centroids of triangles *Tm* + and *Tm* from the ones of *Ti* + and *Ti* - . In similar way:

$$\log \left( \vec{r}\_{m}^{c.c} \right) = -\frac{1}{\vec{\gamma}\omega 4\pi} \mathbf{f}\_{T\_{i}^{z}} \mathbf{f}\_{i} \nabla \cdot \bullet \vec{f}\_{i} \, \frac{e^{-\vec{\gamma} \cdot \mathbf{R}\_{m}^{z}}}{R\_{m}^{z}} dS' = \mp \frac{1}{\vec{\gamma}\omega 4\pi} \mathbf{f}\_{i} \mathbf{f}\_{i} \frac{e^{-\vec{\gamma} \cdot \mathbf{R}\_{m}^{z}}}{R\_{m}^{z}}.\tag{18}$$

The double inner product gives a symmetry to the system, the actions of the current *i* on *m* and of *m* on *i* are similar. The total effect on *m* is due to an extended summation of index *i*. In such a way the following equation system holds:

$$\left[\text{ZI}\right] = \text{V}\_{\text{A}} \tag{19}$$

*H* → *i* (*r* <sup>→</sup> )= *jk* <sup>4</sup>*<sup>π</sup>* (*m* → *<sup>i</sup>* ×*r* → ) 1 *<sup>r</sup>* <sup>2</sup> (1 <sup>+</sup>

the dipole model for the current crossing edge *i* (b).

**4. Antenna modeling with method of moments**

Figure 14. Discrete surface MoM model of the monopole antenna.

<sup>0</sup> 0.02 0.04


**Figure 14.** Discrete surface MoM model of the monopole antenna.





Figure 15. Discrete surface MoM model of the triangular and spherical antennas.

0.05

0.1

discharge are metallically linked by the discharge path.

0

Figure 16. Discrete surface MoM model of the emitting source.

the linking edge is studied.

tools also create the Delaunay triangulation.

unknown is set in equation (6).

*E* → *i* (*r* <sup>→</sup> )= *<sup>η</sup>* <sup>4</sup>*<sup>π</sup>* (( (*<sup>r</sup>* <sup>→</sup> ∙ *m* → *i* )*r* → *<sup>r</sup>* <sup>2</sup> - *m* → *i* )( *jk <sup>r</sup>* + 1 *<sup>r</sup>* <sup>2</sup> +

a

1 *jkr* )*e*- *jkr*

> (*r* <sup>→</sup> ∙ *m* → *i* )*r* → *<sup>r</sup>* <sup>4</sup> (1 <sup>+</sup>

1 *jk <sup>r</sup>* <sup>3</sup> ) + 2

**Figure 13.** A triangular antenna divided into triangular subparts, in which internal current are shown (a). Particular of

The aim of this study is to create a model to describe the behaviour of PD signal probe. Antennas can be easily designed with commercial software like Matlab® in which suitable

The simplest antenna to be designed is the monopole antenna, shown in Fig.14. In order to create the 20 cm long monopole two rectangles have been created, one thin in which only a line of triangle is needed, another bulky.By considering the possible interactions between the two triangles, problems arise. There is an edge shared between three triangles (red, yellow, blue). The current of this edge will be divided into 2 sub-currents and another row and

, (24)

Ti -

Building Partial Discharge Signal Wireless Probes

o

ri - ri +

r

. (25)

b

P

http://dx.doi.org/10.5772/58840

191

1 *jkr* ))*e*- *jkr*

li

Ti +

Similarly all non-planar antennas are treated. In Fig. 15 the triangular antenna and the spherical one are shown. The

<sup>l</sup><sup>i</sup> <sup>J</sup>i1

Ji2

a b

In Fig.16 an electronic power device is sketched: two metal plates, normally at different voltage, during a partial



0

0

0.05

The numerical model allows to feed all the antennas in two different ways. The first one uses a voltage generator, set between the antenna (upper object) and ground plate (bottom metallic plate). The other way takes into account the illuminating external field on the surface of the whole object. In order to simulate the transmission of fault EM signal and its receiving, the power device is fed by generator and the other antennas are used as probes, so the provided voltage on

triangular antennas simulates a planar 20 cm equilateral triangle. The sphere has 6 cm diameter.

0.05 0 0.02 0.04 0.06

with *Z* = *Zm*,*<sup>i</sup>* ,

$$\mathcal{L}\_{m,l} = \frac{1}{2} \left\| \omega \frac{\mu}{8\pi} \frac{e^{\cdot \cdot \cdot \mathbb{R}\_{m}^{+}}}{R\_{m}^{+}} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{m\boldsymbol{\ell}} \rho\_{\boldsymbol{m}}^{c\*} + \frac{1}{2} \left\| \omega \frac{\mu}{8\pi} \frac{e^{\cdot \cdot \cdot \mathbb{R}\_{m}^{-}}}{R\_{m}^{-}} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{m\boldsymbol{\ell}} \rho\_{\boldsymbol{m}}^{c\*} + \frac{1}{j\omega 4\pi} \frac{e^{\cdot \cdot \cdot \mathbb{R}\_{m}^{+}}}{R\_{m}^{+}} l\_{m\boldsymbol{\ell}} \mathbb{I}\_{\boldsymbol{i}} - \frac{1}{j\omega 4\pi} \frac{e^{\cdot \cdot \cdot \mathbb{R}\_{m}^{-}}}{R\_{m}^{-}} l\_{m\boldsymbol{\ell}} \mathbb{I}\_{\boldsymbol{i}} \right\| \rho\_{\boldsymbol{i}}^{c\*} l\_{m\boldsymbol{\ell}} \rho\_{\boldsymbol{i}}^{c\*} l\_{m\boldsymbol{\ell}} \rho\_{\boldsymbol{m}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{m}} \rho\_{\boldsymbol{m}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}} \rho\_{\boldsymbol{i}}^{c\*} l\_{\boldsymbol{i}$$

$$J = \begin{bmatrix} J\_1 \\ \vdots \\ J\_N \end{bmatrix} \tag{21}$$

$$V = \begin{bmatrix} \frac{1}{2} E\_{inc}^{1,c^{\ast\ast}} l\_1 \rho\_1^{c\ast} + \frac{1}{2} E\_{inc}^{1,c^{\ast\ast}} l\_1 \rho\_1^{c\ast} \\ \vdots \\ \frac{1}{2} E\_{inc}^{N\_{\cdot},c^{\ast\ast}} l\_N \rho\_N^{c\ast} + \frac{1}{2} E\_{inc}^{N\_{\cdot},c^{\ast\ast}} l\_N \rho\_N^{c\ast} \end{bmatrix} \tag{22}$$

The vector *V* contains in each row the effect of the incident electric field due to an external source; if the source is an internal voltage generator, all rows are zero except the one in which is enforced the electric field due to the generator, approximated by the imposed voltage divided by the distance between centroids.

Once defined equation 19, the unknown currents *J* can be evaluated on the whole surface of the object (see Fig.13.a). Each current starts and ends on a centroid. A dipole moment is definable for each current *Ji* crossing an internal edge of length *li* (see Fig.11.a):

$$
\vec{m}\_i = \mathbf{f}\_{T\_i^+ + T\_i^-} \mathbf{J}\_i \vec{f}\_i \, dS \, \vec{S} = \frac{1}{2} l\_i \mathbf{J}\_i \left( \vec{\rho}\_i^+ + \vec{\rho}\_i^- \right) = l\_i \mathbf{J}\_i \left( \vec{r}\_i^- - \vec{r}\_i^+ \right) \tag{23}
$$

The radiated magnetic and electric field, due the dipole of current *Ji* , in a generic point P (Fig. 12.b) located at a distance *r* <sup>→</sup> , are given by:

**Figure 13.** A triangular antenna divided into triangular subparts, in which internal current are shown (a). Particular of the dipole model for the current crossing edge *i* (b).

$$
\overrightarrow{H}\_{i}(\overrightarrow{r}) = \frac{jk}{4\pi} (\overrightarrow{m}\_{i} \times \overrightarrow{r}) \frac{1}{r^{2}} \left(1 + \frac{1}{jkr}\right) e^{-jkr} \,\tag{24}
$$

$$\overrightarrow{E}\_{i}(\overrightarrow{r}) = \frac{\eta}{4\pi} \left\{ \left( \frac{(\overrightarrow{r}\bullet\overrightarrow{m}\_{i})\overrightarrow{r}}{r^{2}} - \overrightarrow{m}\_{i} \right) \left( \frac{jk}{r} + \frac{1}{r^{2}} + \frac{1}{jkr^{3}} \right) + 2\frac{(\overrightarrow{r}\bullet\overrightarrow{m}\_{i})\overrightarrow{r}}{r^{4}} \left( 1 + \frac{1}{jkr} \right) \right\} e^{-jkr}.\tag{25}$$

#### **4. Antenna modeling with method of moments**

where *Rm* ±

In similar way:

190 Progress in Compact Antennas

with *Z* = *Zm*,*<sup>i</sup>* ,

<sup>2</sup> *<sup>j</sup><sup>ω</sup> <sup>μ</sup>* 8*π e* - *jk Rm* + *Rm* <sup>+</sup> *l i ρi c*+ *l <sup>m</sup>ρ<sup>m</sup> <sup>c</sup>*<sup>+</sup> + 1 <sup>2</sup> *<sup>j</sup><sup>ω</sup> <sup>μ</sup>* 8*π e* - *jk Rm* - *Rm* - *l i ρi cl <sup>m</sup>ρ<sup>m</sup> <sup>c</sup>*- + 1 *jω*4*π e* - *jk Rm* + *Rm* <sup>+</sup> *l ml <sup>i</sup>* - <sup>1</sup> *jω*4*π e* - *jk Rm* - *Rm* - *l ml i*

*Zm*,*<sup>i</sup>* <sup>=</sup> <sup>1</sup>

*φ*(*r* → *m c*± ) <sup>=</sup> - <sup>1</sup> *<sup>j</sup>ω*4*<sup>π</sup> ∫ Ti* ±*Ji* ∇ ∙ *f* → *i e* - *jk Rm* ± *Rm*

a way the following equation system holds:

is the distance of the centroids of triangles *Tm*

+ and *Tm* -

1 *<sup>j</sup>ω*4*<sup>π</sup> Ji l i e* - *jk Rm* ± *Rm*

*ZJ* =*V* , (19)

, (21)

(see Fig.11.a):

) (23)

, in a generic point P (Fig.

<sup>±</sup> *dS*'= ∓

The double inner product gives a symmetry to the system, the actions of the current *i* on *m* and of *m* on *i* are similar. The total effect on *m* is due to an extended summation of index *i*. In such

*J* =

*V* =

divided by the distance between centroids.

*m* → *<sup>i</sup>* =*∫ Ti* + +*T <sup>i</sup>* -*Ji f* → *i dS* ' = 1 <sup>2</sup> *liJi* (*ρ* → *i* <sup>+</sup> + *ρ* → *i* - ) =*liJi* (*r* → *i* - - *r* → *i* +

definable for each current *Ji*

12.b) located at a distance *r*

1 <sup>2</sup> *Einc* 1, *c*+ *l* <sup>1</sup>*ρ*<sup>1</sup> *<sup>c</sup>*<sup>+</sup> + 1 <sup>2</sup> *Einc* 1, *cl* <sup>1</sup>*ρ*<sup>1</sup> *c*-

The radiated magnetic and electric field, due the dipole of current *Ji*

<sup>→</sup> , are given by:

1 <sup>2</sup> *Einc N* , *c*+ *lN ρ<sup>N</sup> <sup>c</sup>*<sup>+</sup> + 1 <sup>2</sup> *Einc N* , *clN ρ<sup>N</sup> c*-

*J*1 ⋮ *J N*

⋮

The vector *V* contains in each row the effect of the incident electric field due to an external source; if the source is an internal voltage generator, all rows are zero except the one in which is enforced the electric field due to the generator, approximated by the imposed voltage

Once defined equation 19, the unknown currents *J* can be evaluated on the whole surface of the object (see Fig.13.a). Each current starts and ends on a centroid. A dipole moment is

crossing an internal edge of length *li*

from the ones of *Ti*

+ and *Ti* - .

, (20)

(22)

<sup>±</sup> . (18)

The aim of this study is to create a model to describe the behaviour of PD signal probe. Antennas can be easily designed with commercial software like Matlab® in which suitable tools also create the Delaunay triangulation.

The simplest antenna to be designed is the monopole antenna, shown in Fig.14. In order to create the 20 cm long monopole two rectangles have been created, one thin in which only a line of triangle is needed, another bulky.By considering the possible interactions between the two triangles, problems arise. There is an edge shared between three triangles (red, yellow, blue). The current of this edge will be divided into 2 sub-currents and another row and unknown is set in equation (6).

triangular antennas simulates a planar 20 cm equilateral triangle. The sphere has 6 cm diameter.

0.05 0 0.02 0.04 0.06

Figure 15. Discrete surface MoM model of the triangular and spherical antennas.

0.05

0.1

discharge are metallically linked by the discharge path.

0

Figure 16. Discrete surface MoM model of the emitting source.

the linking edge is studied.

Similarly all non-planar antennas are treated. In Fig. 15 the triangular antenna and the spherical one are shown. The

In Fig.16 an electronic power device is sketched: two metal plates, normally at different voltage, during a partial



0

0

0.05

The numerical model allows to feed all the antennas in two different ways. The first one uses a voltage generator, set between the antenna (upper object) and ground plate (bottom metallic plate). The other way takes into account the illuminating external field on the surface of the whole object. In order to simulate the transmission of fault EM signal and its receiving, the power device is fed by generator and the other antennas are used as probes, so the provided voltage on

Figure 14. Discrete surface MoM model of the monopole antenna. **Figure 14.** Discrete surface MoM model of the monopole antenna.




0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2


Similarly all non-planar antennas are treated. In Fig. 15 the triangular antenna and the spherical one are shown. The triangular antennas simulates a planar 20 cm equilateral triangle. The sphere has 6 cm diameter. Figure 14. Discrete surface MoM model of the monopole antenna. -0.04 -0.02 -0.04

triangular antennas simulates a planar 20 cm equilateral triangle. The sphere has 6 cm diameter.

Similarly all non-planar antennas are treated. In Fig. 15 the triangular antenna and the spherical one are shown. The

<sup>l</sup><sup>i</sup> <sup>J</sup>i1

As previously explained the source of the EM field illuminating the probes is given by the discharge current flowing through the discharge path, which is due to the ionized air gap inside of insulation, and embraces previously damaged materials and metallic surfaces. A full characterization of the exact discharge channel is beyond the scope of this work and also for sake of simplicity the discharge path is simulated with a metallic wire in which the discharge current flows unhindered. Then neglecting particular behaviour to the geometry of the discharge channel, scrupulous attention is paid to the waveform of the current, the harmonic

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

193


time [ms]

<sup>0</sup> <sup>20</sup> <sup>40</sup> <sup>60</sup> <sup>80</sup> <sup>100</sup> <sup>120</sup> <sup>0</sup>

frequency [MHz]

In Fig.17, a PD signal is represented, which is characterized by an underdamped oscillation. This latter propagates until the energy balance is once again stable, and the charges back to be imprisoned in the dielectric. In order to study the transient PD signal Discrete Fast Fourier Transform (DFFT) has to be used. The PD signals may have a faster or slower evolving, but roughly the spectrum of interest is contained in a [0-200MHz] window. This aspect allows characterizing the type of receiving probes. Unfortunately, in the environment in which PD signals are generated, other signals, that should affect in a limited way the diagnostic tool, are present. PD signal of Fig.17intentionally contains a part of noise that will interact with the

Once defined the behaviour of the PD current, it is possible to characterize the behaviour of the source system, focusing at 20 cm from the discharge, chosen as a compromise between the possibility of bringing a wireless sensor to the faulty part and the need to maintain a safe distance or in other cases distance imposed by the presence of screens, shields and spacers.The PD discharge path has a normal direction with respect to the ground plate and an unitary current source is applied within the range of [0-1000MHz]. The simulated z-component of the

**Figure 17.** Discharge current and Fourier spectrum. Noise is on 90-110 MHz range.

content, which will characterize the fidelity of the receiving probes, is studied.

0

2

Furier's com

probes receivers.

4

 ponents

<sup>6</sup> x 10-4

0.01 0.02

current [A ]

Ji2

a b

Figure 15. Discrete surface MoM model of the triangular and spherical antennas. **Figure 15.** Discrete surface MoM model of the triangular and spherical antennas.

<sup>0</sup> 0.02 0.04

In Fig.16 an electronic power device is sketched: two metal plates, normally at different voltage, during a partial discharge are metallically linked by the discharge path. In Fig.16 an electronic power device is sketched: two metal plates, normally at different voltage, during a partial discharge are metallically linked by the discharge path.

**Figure 16.** Discrete surface MoM model of the emitting source.

the linking edge is studied.

The numerical model allows to feed all the antennas in two different ways. The first one uses a voltage generator, set between the antenna (upper object) and ground plate (bottom metallic plate). The other way takes into account the illuminating external field on the surface of the whole object. In order to simulate the transmission of fault EM signal and its receiving, the power device is fed by generator and the other antennas are used as probes, so the provided voltage on the linking edge is studied.

As previously explained the source of the EM field illuminating the probes is given by the discharge current flowing through the discharge path, which is due to the ionized air gap inside of insulation, and embraces previously damaged materials and metallic surfaces. A full characterization of the exact discharge channel is beyond the scope of this work and also for sake of simplicity the discharge path is simulated with a metallic wire in which the discharge current flows unhindered. Then neglecting particular behaviour to the geometry of the discharge channel, scrupulous attention is paid to the waveform of the current, the harmonic content, which will characterize the fidelity of the receiving probes, is studied.

**Figure 17.** Discharge current and Fourier spectrum. Noise is on 90-110 MHz range.

Similarly all non-planar antennas are treated. In Fig. 15 the triangular antenna and the spherical one are shown. The triangular antennas simulates a planar 20 cm equilateral triangle. The

Figure 15. Discrete surface MoM model of the triangular and spherical antennas.

0.05

0.1

In Fig.16 an electronic power device is sketched: two metal plates, normally at different voltage,

The numerical model allows to feed all the antennas in two different ways. The first one uses a voltage generator, set between the antenna (upper object) and ground plate (bottom metallic plate). The other way takes into account the illuminating external field on the surface of the whole object. In order to simulate the transmission of fault EM signal and its receiving, the power device is fed by generator and the other antennas are used as probes, so the provided

discharge are metallically linked by the discharge path.

during a partial discharge are metallically linked by the discharge path.

**Figure 15.** Discrete surface MoM model of the triangular and spherical antennas.

0

Figure 16. Discrete surface MoM model of the emitting source.

the linking edge is studied.

voltage on the linking edge is studied.

**Figure 16.** Discrete surface MoM model of the emitting source.

triangular antennas simulates a planar 20 cm equilateral triangle. The sphere has 6 cm diameter.

0.05 0 0.02 0.04 0.06

Similarly all non-planar antennas are treated. In Fig. 15 the triangular antenna and the spherical one are shown. The

<sup>l</sup><sup>i</sup> <sup>J</sup>i1

Ji2

a b

In Fig.16 an electronic power device is sketched: two metal plates, normally at different voltage, during a partial



0

0

0.05

The numerical model allows to feed all the antennas in two different ways. The first one uses a voltage generator, set between the antenna (upper object) and ground plate (bottom metallic plate). The other way takes into account the illuminating external field on the surface of the whole object. In order to simulate the transmission of fault EM signal and its receiving, the power device is fed by generator and the other antennas are used as probes, so the provided voltage on

sphere has 6 cm diameter. Figure 14. Discrete surface MoM model of the monopole antenna.

<sup>0</sup> 0.02 0.04



192 Progress in Compact Antennas




In Fig.17, a PD signal is represented, which is characterized by an underdamped oscillation. This latter propagates until the energy balance is once again stable, and the charges back to be imprisoned in the dielectric. In order to study the transient PD signal Discrete Fast Fourier Transform (DFFT) has to be used. The PD signals may have a faster or slower evolving, but roughly the spectrum of interest is contained in a [0-200MHz] window. This aspect allows characterizing the type of receiving probes. Unfortunately, in the environment in which PD signals are generated, other signals, that should affect in a limited way the diagnostic tool, are present. PD signal of Fig.17intentionally contains a part of noise that will interact with the probes receivers.

Once defined the behaviour of the PD current, it is possible to characterize the behaviour of the source system, focusing at 20 cm from the discharge, chosen as a compromise between the possibility of bringing a wireless sensor to the faulty part and the need to maintain a safe distance or in other cases distance imposed by the presence of screens, shields and spacers.The PD discharge path has a normal direction with respect to the ground plate and an unitary current source is applied within the range of [0-1000MHz]. The simulated z-component of the electric field at a distance of 20 cm, is shown in Fig.18. The electric field profile is growing with the frequency, similar to the one of a dipole with its dimension equal to the discharge path.

In order to compare the different antennas behaviour a particular transfer function is consid‐ ered:

$$H = \frac{V\_{reiral}}{E\_{incident}} \left[ m \right] \tag{26}$$

<sup>0</sup> <sup>50</sup> <sup>100</sup> <sup>150</sup> <sup>0</sup>

<sup>0</sup> <sup>100</sup> <sup>200</sup> <sup>300</sup> <sup>400</sup> <sup>500</sup> <sup>600</sup> <sup>700</sup> <sup>800</sup> <sup>900</sup> <sup>1000</sup> <sup>0</sup>

Frequency [MHz]

monopole triangle sphere

monopole triangle sphere

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

195

Frequency [MHz]

In Fig.20, the different transfer functions are plotted in the frequency range going from 0 to 200 MHz. The spherical and triangular antennas show an almost flat behaviour while the monopole is characterized by a significant decrease in amplitude. It is easily to argue that these different profiles have a different impact on the shape of the recorded PD signals. The selection of a specific profile depends on the application, e.g. PD measurements in dc, ac and rectangular (adjustable speed drive applications) voltage supply. A procedure able to evaluate the signal distortion of the different antennas has been developed to select the antenna able to fit the

separation and identification requirements.

**Figure 20.** Transfer function of the three antennas for the range of 0-150MHz.

transfer function [m

 ]

**Figure 19.** Transfer function of the three antennas.

transfer function [m

 ]

where *Vreceived(ω)* is the voltage response of antennas, illuminated by the incident electromag‐ netic field.

The incident EM field due to the PD current, which is function of the signal frequency content, has to be computed. The received voltages can be computed by simulating the response of the antenna illuminated by the incident EM field. The voltage response is due to both source and antenna impedance behaviour. The three transfer functions corresponding to the three antenna typologies are plotted in Fig. 19considering the above mentioned frequency range.

As it can be seen, the monopole antenna provides the highest gain but its profile is character‐ ized by a resonant peak centered at about 450 MHz. The triangular antenna shows a resonant peak centered in about 400 MHz but with lower amplitude with respect to the former. On the contrary, the spherical antenna shows a flat bandwidth with the lowest amplitude among the three.

**Figure 19.** Transfer function of the three antennas.

electric field at a distance of 20 cm, is shown in Fig.18. The electric field profile is growing with the frequency, similar to the one of a dipole with its dimension equal to the discharge path.

<sup>0</sup> <sup>50</sup> <sup>100</sup> <sup>150</sup> <sup>200</sup> <sup>1</sup>

Frequency [M Hz]

**Figure 18.** *Ez* field generated by the discharge at a horizontal distance of 20 cm. Only the 0-200 MHz interval is shown.

In order to compare the different antennas behaviour a particular transfer function is consid‐

where *Vreceived(ω)* is the voltage response of antennas, illuminated by the incident electromag‐

The incident EM field due to the PD current, which is function of the signal frequency content, has to be computed. The received voltages can be computed by simulating the response of the antenna illuminated by the incident EM field. The voltage response is due to both source and antenna impedance behaviour. The three transfer functions corresponding to the three antenna

As it can be seen, the monopole antenna provides the highest gain but its profile is character‐ ized by a resonant peak centered at about 450 MHz. The triangular antenna shows a resonant peak centered in about 400 MHz but with lower amplitude with respect to the former. On the contrary, the spherical antenna shows a flat bandwidth with the lowest amplitude among the

*Eincident <sup>m</sup>* (26)

*<sup>H</sup>* <sup>=</sup> *Vreceived*

typologies are plotted in Fig. 19considering the above mentioned frequency range.

1.5 2 2.5 3

3.5 4 4.5 5

Ez [V /m]

194 Progress in Compact Antennas

ered:

netic field.

three.

5.5 x 10-4

**Figure 20.** Transfer function of the three antennas for the range of 0-150MHz.

In Fig.20, the different transfer functions are plotted in the frequency range going from 0 to 200 MHz. The spherical and triangular antennas show an almost flat behaviour while the monopole is characterized by a significant decrease in amplitude. It is easily to argue that these different profiles have a different impact on the shape of the recorded PD signals. The selection of a specific profile depends on the application, e.g. PD measurements in dc, ac and rectangular (adjustable speed drive applications) voltage supply. A procedure able to evaluate the signal distortion of the different antennas has been developed to select the antenna able to fit the separation and identification requirements.

#### **5. Transmission of Signals**

In order to test the different antennas, bursts of damped pulses were simulated with super‐ imposed high frequency noise in the frequency range of 0-200MHz. The performances of the three antenna configurations were evaluated considering the simulated received signals. An example of the normalized voltage response of the three antennas is plotted in Figure 21. According to Fig. 19, it was found an amplitude reduction of-10 dB and-14 dB for the triangular and spherical antennas with respect to the monopole probe.

where *s(t)* represents the time behaviour of the simulated PD signal in discrete time domain. In order to compare the shape of signals with different energy, the normalized ACF in discrete

> *s*(*i*)*s*(*i* + *k* ) *N*

Where *s(i)* is the signal digital sample, *N* is the number of digital samples, *ES* is the energy of

1

*ES* , (28)

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

*<sup>N</sup>* - <sup>1</sup> , (30)

, (31)

*d<sup>ρ</sup>* = 2(1 - *ρ*) (32)

(29)

197

*RS k* =

∑ *i*=1 *N*

*ES* = ∑ *i N s*(*i*)2 *N*

Fig.s22 and 23 show the NACF of the signals reported in Fig. 17 and 21, respectively.

The covariance between two NACFs can be evaluated by the following expression:

∑ *k* =1 *N RS* <sup>1</sup> (*k* ) - *R* -

*<sup>ρ</sup>* <sup>=</sup> *cov*(*RS* 1, *RS* <sup>2</sup>) *var*(*RS* 1) ∙ *var*(*RS* 2)

Previous investigations established a dρ<0.4 as a good threshold level in the PD shape com‐ parison. The differences between the transmitted and the received signals using the three different configurations, evaluated considering the distance of equation 32, are reported in

*cov*(*RS* 1, *RS* 2) =

where *var*(*RS* 1) and *var*(*RS* 2) are the variance of the NACFs.

The distance in the NACF metric can be assumed as:

*S* 2 are the mean value of NACFs.

The Pearson correlation index, ρ, is used to compare the different NACF with the emitted one.

*<sup>S</sup>* <sup>1</sup> ∙ *RS* <sup>2</sup>

(*k* ) - *R* - *S* 2

time domain has been here adopted:

the signal:

where *R* -

Table 1.

*S* 1 and *R* -

The Pearson correlation index is:

**Figure 21.** Time profiles of the received voltage.

Moreover, the high frequency noise has the highest influence on the pulse signal received by the triangular antenna since its transfer function reduces the amplitude of the pulse signals and amplifies the superimposed noise. The opposite considering the monopole configuration, while the flat bandwidth of the spherical antenna reduces both the amplitude and the super‐ imposed noise of the received pulse signal. A dedicated procedure has been developed in order to validate this empirical evaluation.

The proposed validation procedure is based on the classification of the autocorrelation function (ACF) of the PD signals, via the K-Means classification algorithm [10]. The ACF has been proposed for its ability to synthesize both the time and frequency domain features as well as to be less affected by superimposed high-frequency noise, random truncation of the pulsetail (frequency leakage) and different trigger activation. The ACF is evaluated by taking into account the following expression:

$$R\_S = \bigcup\_{\omega} \mathbf{s}(t)\mathbf{s}(t+\tau)dt\tag{27}$$

where *s(t)* represents the time behaviour of the simulated PD signal in discrete time domain. In order to compare the shape of signals with different energy, the normalized ACF in discrete time domain has been here adopted:

$$\left[R\_S \prod\_{i=1}^{N} k\right] = \frac{\sum\_{i=1}^{N} s(i)s(i+k)}{N} \frac{1}{E\_S \ \prime} \tag{28}$$

Where *s(i)* is the signal digital sample, *N* is the number of digital samples, *ES* is the energy of the signal:

$$E\_S = \frac{\sum\_{s(i)}^{N} s(i)^2}{N} \tag{29}$$

Fig.s22 and 23 show the NACF of the signals reported in Fig. 17 and 21, respectively.

The Pearson correlation index, ρ, is used to compare the different NACF with the emitted one. The covariance between two NACFs can be evaluated by the following expression:

$$cov\left(R\_{S1\nu}, R\_{S2}\right) = \frac{\stackrel{\circ}{\sum}\left[R\_{S1}(k) \cdot R\_{S1}\right] \bullet \left[R\_{S2}(k) \cdot R\_{S2}\right]}{N \cdot 1} \tag{30}$$

where *R* - *S* 1 and *R* - *S* 2 are the mean value of NACFs.

The Pearson correlation index is:

**5. Transmission of Signals**

196 Progress in Compact Antennas

0 1

0 1

0 1

**Figure 21.** Time profiles of the received voltage.

to validate this empirical evaluation.

account the following expression:

and spherical antennas with respect to the monopole probe.

In order to test the different antennas, bursts of damped pulses were simulated with super‐ imposed high frequency noise in the frequency range of 0-200MHz. The performances of the three antenna configurations were evaluated considering the simulated received signals. An example of the normalized voltage response of the three antennas is plotted in Figure 21. According to Fig. 19, it was found an amplitude reduction of-10 dB and-14 dB for the triangular


monopole

t, [ms]

triangle


t, [ms]

sphere


t, [ms]

Moreover, the high frequency noise has the highest influence on the pulse signal received by the triangular antenna since its transfer function reduces the amplitude of the pulse signals and amplifies the superimposed noise. The opposite considering the monopole configuration, while the flat bandwidth of the spherical antenna reduces both the amplitude and the super‐ imposed noise of the received pulse signal. A dedicated procedure has been developed in order

The proposed validation procedure is based on the classification of the autocorrelation function (ACF) of the PD signals, via the K-Means classification algorithm [10]. The ACF has been proposed for its ability to synthesize both the time and frequency domain features as well as to be less affected by superimposed high-frequency noise, random truncation of the pulsetail (frequency leakage) and different trigger activation. The ACF is evaluated by taking into

*s*(*t*)*s*(*t* + *τ*)*dt* (27)

*RS* = *∫* -∞ ∞

$$\rho = \frac{cov\left(R\_{\text{S}\,\text{V}} \parallel R\_{\text{S}\,\text{2}}\right)}{\sqrt{var\left(R\_{\text{S}\,\text{1}}\right)} \bullet \sqrt{var\left(R\_{\text{S}\,\text{2}}\right)}},\tag{31}$$

where *var*(*RS* 1) and *var*(*RS* 2) are the variance of the NACFs.

The distance in the NACF metric can be assumed as:

$$d\_{\rho} = \sqrt{2(1 - \rho)}\tag{32}$$

Previous investigations established a dρ<0.4 as a good threshold level in the PD shape com‐ parison. The differences between the transmitted and the received signals using the three different configurations, evaluated considering the distance of equation 32, are reported in Table 1.

the PD signal. Due to the different attenuation in the different frequency ranges (Fig. 19), the monopole antenna tends to modify differently the shape of signals generated by the same source. This can induce an over classification in separation procedure that is, the generation of more classes than PD sources. The spherical antenna, due to its flat bandwidth, seems to be

Building Partial Discharge Signal Wireless Probes

http://dx.doi.org/10.5772/58840

199

The table I shows that the most suitable probes for the identification of partial discharge signals are the monopole probe and the spherical probe. This evaluation was obtained using one of the simplest classification algorithms, certainly other more advanced algorithms may change this rating, and also the use of digital filters in the same implementation may overturn the comparison. The use of filters or the use of more elaborate techniques of signal processing goes beyond the intentions of this chapter. Therefore, the comments will be expressed only in

The monopole probe has a high gain, superior to the others. This is one of the key aspects in the evaluation of partial discharge phenomena that makes it attractive: the discharges take place within structures often shielded, for which the outgoing signals are weak. The use of too attenuating probe could make the measurement insufficient and the monopole can solve this issue. However, the use of this probe in the industrial environment is not wise since it could involve the absence of higher harmonics. Indeed, their contribution can modify the shape of the signal and even changing the initial instant of the discharge. Such behaviour, in the diagnosis of PD phenomenon, is not acceptable, since the determination of the type of discharge (internal, surface, corona, etc..) is based on the presence of a greater or lesser occurrence of discharges at certain angles of the sine wave function. One way to change the behaviour of the monopole is to change its length, but this length cannot be sufficiently modified to remove the effect of resonance. A simple way to try to mitigate the effect of the resonance is to insert appropriate RLC circuits. UWB antennas are actually filters, which introduce dependences between the transmitter and the receiver frequencies. Therefore such kind of antennas can be treaded. Generally, antennas have linear and passive behaviour, and input impedance can be represented by means of the canonical forms of Foster [28-30].The structure shown in Fig.24 is suitable to represent electric antenna as the electric dipole and monopole.In order to modify the behaviour of the antenna, the electric model of the antenna is required. The numerical implementation of the antenna gives the complex impedance used in the equation 19. The equivalent circuit has to show a similar behaviour for all the harmonic of interest. By following the steps in [28] the value of C0 is set by taking the value of the capacity presented by the antenna at low frequencies; the value of L0 is initially set to cancel the imaginary part of impedance at the first resonance frequency. It is required that the first RLC parallel branch has a resonant frequency corresponding to the first maximum of the real part

the most suitable wireless probe for ac applications.

relation to the analogical performance of probes.

**6. Discussion**

**6.1. Monopole antenna**

**Figure 22.** NACF of the PD current of Figure17.

**Figure 23.** NACFs of the transmitted signals of Figure21.


As it can be seen, the monopole and the spherical antenna preserve the pulse shape in almost similar way while the triangle's received signal is affected by noise and it is groped far from the PD signal. Due to the different attenuation in the different frequency ranges (Fig. 19), the monopole antenna tends to modify differently the shape of signals generated by the same source. This can induce an over classification in separation procedure that is, the generation of more classes than PD sources. The spherical antenna, due to its flat bandwidth, seems to be the most suitable wireless probe for ac applications.
