Infrared Fault Diagnosis in Transformers
Infrared Fault Diagnosis in Transformers
DOI: 10.1049/hve2.12405
ORIGINAL RESEARCH
- -Revised: 8 August 2023 Accepted: 16 September 2023
- High Voltage
1
Department of Electrical Engineering, Faculty of Abstract
Engineering at Shoubra, Benha University, Cairo,
Egypt
Condition monitoring of the insulating system within power transformers has a massive
2
importance according to the electrical utilities. Dissolved gas analysis (DGA) is frequently
Department of Electrical Engineering and
Automation, School of Electrical Engineering, Aalto used for this purpose. However, DGA lacks the necessary level of accuracy to identify all
University, Espoo, Finland equipment faults, particularly in their initial stages of degradation. Also, it does not have
3
Department of Electrical Power Engineering, the capability for real‐time monitoring and relies on manual sampling and laboratory
Egypt‐Japan University of Science and Technology testing, causing potential delays in fault identification. Additionally, the interpretation of
(E‐JUST), New Borg El‐Arab City, Egypt
DGA data necessitates specialised expertise, which may pose difficulties for smaller
4
Department of Electrical Power and Machines entities that have limited access to resources. Therefore, the contribution of this research
Engineering, Faculty of Engineering, Tanta
University, Tanta, Egypt
is to use infrared spectroscopy measurements as a new effective technique substituting
the DGA method for fault diagnosis in insulating oil. The inception faults that were
Correspondence considered in this study were the electrical fault (discharges of high energy) and the
Mohamed M. F. Darwish and Matti Lehtonen. thermal fault (300°C < Temperature < 700°C). Regarding that, two test cells were crafted
Email: [email protected] and especially for serving the simulation processes inside the laboratory for both types of
[email protected]
inception faults. Subsequently, six samples of pure paraffinic mineral oil were taken to be
degraded in the laboratory. Following that, all of them besides another sample that were
Associate Editor: Qing Yang
not subjected to any kind of faults were taken to be examined by Fourier transform
Funding information infrared (FTIR) spectroscopy to obtain an overview of the oil's behaviour in each fault
Science, Technology & Innovation Funding case. After that, the FTIR analysis was initially verified utilising the DGA method. Then,
Authority, Egypt, Grant/Award Number: 31292 for further affirmation, the dielectric dissipation factor (DDF) for all samples was
measured. In the final analysis, the verification tests provide experimental evidence about
the outperformance of this new optical technique in detecting the transformer's inception
faults in addition to proving its potential for being a superior alternative to the well‐
known traditional diagnostic techniques.
-
This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is
properly cited.
© 2024 The Authors. High Voltage published by John Wiley & Sons Ltd on behalf of The Institution of Engineering and Technology and China Electric Power Research Institute.
(DGA). This technique is capable of predicting the type of since only one case of degradation severity was considered for
inception fault that occurred based on the concentrations each fault type. Also, the FTIR analysis was confirmed with
of certain fault gases dissolved in the oil. The concentrations of only one diagnostic technique.
such gases are measured via a device called gas chromatog- According to the plentiful survey mentioned above and the
raphy (GC), and then the analysis is conducted utilising one or authors' knowledge, no one in the previous works has deeply
more of the methods that lie under the umbrella of the DGA discussed the capability of monitoring the transformer incep-
technique. Generally, there are several traditional DGA tion faults based on the oil insulation's spectrum utilising any
methods such as the Doernenburg ratio, Duval triangles, key of the various optical spectroscopy techniques. In this regard,
gas, the international electrotechnical commission (IEC) ratio, the new proposed investigation has covered the lack in this
Duval pentagons, and Roger's ratio [4–6]. However, the ac- research area by utilising the FTIR spectroscopy technique to
curacy of the mentioned DGA methods is somewhat low comprehensively study the effects of various types of trans-
compared to the other diagnostic techniques. Fortunately, after former inception faults on the transformer oil's spectral
the implementation of the internet of things (IoT), artificial behaviour.
intelligence (AI), and machine learning (ML) in such traditional Therefore, the contribution of this paper is as follows:
DGA methods, diagnostic accuracy was much improved [7–
10]. Despite the high dependence on the DGA technique in � Degrading six samples of pure paraffinic mineral oil elec-
detecting transformer inception faults, it has many drawbacks. trically and thermally inside the laboratory utilising special
Actually, it is a destructive expensive technique that requires an test setups.
experienced operator to conduct the measurements; also, it is � Examining all of the degraded samples besides a reference
time‐consuming since the GC device takes a large amount of sample by FTIR spectroscopy according to the ASTM
measuring time [11, 12]. Further, it is not an eco‐friendly E1252‐98 standard in order to study the oil's spectral
technique by the way as it utilises toxic organic solvents be- behaviour in each fault type [23].
sides other resources. These disadvantages have forced � Confirming the FTIR analysis initially through the DGA
immediately the researchers to find other alternatives to be method utilising a GC according to the ASTM D3612‐02
utilised in monitoring the condition of power transformers. standard [24]. Thereafter, for further verification, the
Lately, optical spectroscopy techniques have attracted most of DDF for all of the oil samples was estimated with an LCR
the attention towards efficient monitoring of the insulating metre according to the IEC 60247 standard [25].
system quality within power transformers. Since such tech- � Checking the reliability of this new optical technique in
niques are non‐destructive, easy to be conducted, cost‐ discriminating between transformer inception faults besides
effective, and eco‐friendly besides having high accuracy in being an alternative to the well‐known traditional diagnostic
material characterisation [13, 14]. methods.
In this regard, the authors in Ref. [15] have tried to predict
the type of fault via the Duval triangle method based on the
intensities of the dissolved gases absorbance peaks utilising the 2 | SAMPLES PREPARATION AND
Fourier transform infrared (FTIR) spectroscopy. Furthermore, EXPERIMENTAL PROCEDURES
in Ref. [16], the transformer oil's quality was evaluated based
on its molecular structure with the aid of FTIR spectroscopy. For ensuring the applicability of the FTIR spectroscopy for
Also, the spectrum of such optical techniques can be correlated being an alternative to the well‐known conventional diagnostic
with the output of some traditional diagnostic techniques such techniques, six samples each of which was 400 mL of pure
as total acid number (TAN) and dielectric dissipation factor mineral oil (Shell Diala S2 ZU‐I Dried) were degraded in the
(DDF) as in Refs. [17, 18]. Moreover, ultraviolet‐visible (UV‐ laboratory utilising experimental setups that were crafted
Vis) spectroscopy was utilised in detecting the transformer oil especially for conducting the following types of faults. For
quality based on the by‐products formed within the oil due to more details, three of the samples named S2, S3, and S4 were
thermal ageing [19]. However, in Ref. [20], the transformer degraded electrically through different numbers of repeated
health index was calculated based on the number of ageing by‐ high‐voltage impulses (discharges with high energy). On the
products that were formed within each oil sample utilising UV‐ other hand, the rest of the samples named S5, S6, and S7 were
Vis spectroscopy. Concerning Raman spectroscopy, it was subjected to various periods of thermal overheating (300°C <
utilised also in determining the thermal ageing stage of oil‐ Temperature < 700°C) to be degraded thermally.
paper insulation [21]. Last but not least, the role of optical
spectroscopy techniques in detecting transformer inception
faults was introduced by our research group in Ref. [22]. The 2.1 | Electrical degradation
FTIR spectra of mineral oil samples that have been electrically
as well as thermally degraded were utilised in the discrimination According to the IEC and the institute of electrical and elec-
between those types of inception faults. However, this inves- tronics engineers (IEEE) standards, the electrical faults that
tigation has some drawbacks that have been fulfiled in this commonly happen within power transformers were classified
current work. These drawbacks are as follows: the oil's spectral into three classes depending on the energy level of the elec-
behaviour for each fault type was not studied comprehensively trical discharges [5, 26]. For more details, the classes under the
23977264, 2024, 2, Downloaded from https://ietresearch.onlinelibrary.wiley.com/doi/10.1049/hve2.12405 by Cochrane Peru, Wiley Online Library on [29/06/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
DARWISH ET AL.
- 321
umbrella of electrical fault type starting from the one with the positive polarity through a needle‐to‐plate electrode configu-
lowest energy level until reaching the one with the highest ration of 1 mm gap distance keeping a time interval of 30 s
energy level are partial discharge, discharges of low energy between impulses. Then, by reaching the end of the degrada-
(D1), and discharges of high energy (D2), respectively. There- tion process, the plunger of the syringe was pressed gently to
fore, to comprehensively study the impact of electrical fault on guarantee that most of the formed gases were dissolved into
mineral oil's behaviour, the worst electrical fault type which is the oil. Furthermore, the impulse voltage applied during this
D2 was conducted in the laboratory via the test platform accelerated impulsive degradation was a standard lightning
shown in Figure 1 according to IEC 60060‐1 [27]. impulse of 70 kV peak value as well as a front time and tail
time of 1.2 and 50 μs, respectively, according to the IEC
60060‐1 [27].
2.1.1 | Test cell specifications Anyway, the three samples S2, S3, and S4 that were utilised
in this kind of degradation were subjected to different numbers
It is clear from the scheme of the electrical fault experimental of repeated electrical impulses. The samples have been
setup depicted in Figure 1 that the test cell is composed of only degraded by 30, 45, and 90 impulses to study the behaviour of
two parts. The first part (upper part) is a PA6 Nylon cover mineral oil comprehensively when subjected to this kind of
consisting of a needle high‐voltage electrode made from steel fault. In addition, the waveforms of all impulses in each elec-
with a 0.01 mm tip radius and a 1.5 mm diameter. Besides, a trical fault condition were monitored carefully during the test
conical funnel with a 20° slope angle is already carved in the utilising the Tektronix DPO 7104C Digital Phosphor Oscil-
cover and connected to a syringe of 100 mL via an air‐bleeding loscope of 10X attenuation based on the measuring system that
pipe with a three‐way valve. The second part (lower part) is a will be described in the next sentences for ensuring that the
PA6 Nylon foundation where at its centre, there is a brass specifications of the high‐voltage impulse were nearly constant
circular plate electrode with a 20 mm diameter earthed through during the deterioration process for all oil samples. For more
an earthing system of 0.5 Ω. Also, upon the foundation the details, the measuring system that was utilised in this work was
body of the test cell in which the oil is placed exists. It is a calibrated to achieve the measuring requirements in the IEC
transparent Perspex vessel of 100 mm inner diameter, 5 mm 60,060‐2 and is composed of the following components.
thickness, 70 mm height, and a total capacity of 550 mL. Also, Firstly, an impulse voltage divider (high voltage capacitor) with
there are two O‐Rings for preventing the evolved gases from a 255 voltage dividing ratio is connected directly in parallel
escaping out of the vessel and the external gases from entering. with the test cell as shown in Figure 1 to reduce the peak of the
Note that the clearance and creepage distances were taken into impulse wave from 70 kV to about 275 V. Secondly, a coaxial
consideration during the design phase for minimising the risk cable (transmission system) which is connected to the divider's
of flashover [16, 28]. low voltage arm for transferring the output signal to an
attenuator of 10X attenuation further reduces the peak level
from 275 V to about 28 V before entering the oscilloscope to
2.1.2 | Experiment procedure prevent overloading the oscilloscope's input circuitry. Note
that the coaxial cable, the attenuator, and the oscilloscope's
For simulating the electrical fault of class D2, firstly, the test cell input each of them have a surge impedance of 50 Ω to provide
was cleaned very well; then, 400 mL of pure mineral oil was approximately a zero‐reflection coefficient (proper impedance
poured slowly into the Perspex vessel. Following that, the test matching throughout the system) and prevent the measured
cell was well‐sealed with bolts and nuts. After the test cell was wave from reflection that can lead to signal distortion, standing
sealed completely, the degradation was performed through the waves, and inefficient power transfer. Subsequently, for
impulse voltage test system shown in Figure 1. The oil within ensuring that the entire waveform is visible and well‐scaled on
the test cell was subjected to repeated high‐voltage impulses of the screen without any distortion or clipping, the oscilloscope
was properly configured by setting up the voltage and time per
division to 10 V and 2 µs, respectively [29].
generating thermal faults in the laboratory was employed to is sufficient to initialise vapour nucleation. Therefore, the heat
conduct the T2 thermal fault with the aid of the experimental energy is dissipated via vaporisation and strong oil convection.
setup depicted in Figure 2 for studying the thermal fault effect Since vaporisation normally absorbs large amounts of heat
on the mineral oil's behaviour. Actually, the basis for choosing energy, as the input power increases, the surface heat flux in-
the T2 thermal fault among the three types of transformer creases but with a minimal increase in the heating element's
thermal faults is as follows. Firstly, the T2 thermal fault is surface temperature. Thus, the heating element's temperature
widely adopted in the literature instead of the T1 and T3 in the nucleate boiling region remains almost constant ‘satu-
thermal faults. Regarding the T1 thermal fault, it was widely rated’ even with input power largely rising and has an upper
adopted in the literature as its effect on the oil's molecular limit temperature at the critical surface heat flux point known
structure is very low, and since the main goal of this work is to as the Nukiyama temperature. This temperature is 330ºC (al-
study the thermal effect on the mineral oil's molecular struc- ways about 40ºC higher than the liquid's boiling point). Hence,
ture behaviour, the T1 thermal fault will not be suitable with the temperature in the nucleate boiling region approximately
our research goal. Concerning the T3 thermal fault, it is diffi- remains between 320ºC and 330ºC. Following this stage, the
cult in the laboratory to go beyond the T2 thermal fault level by liquid should be in the transition boiling stage; however, in a
the immersed heating method. The underlying reason behind practical situation if the input power to the heater is individ-
that returns to the pool boiling phenomenon that limits the ually controlled, the heating process will follow the first and
temperature of the immersed heating element as reported in second stages and once the Nukiyama temperature is reached.
Refs. [28, 30, 31]. Where the input power approaches 1200 W or surface heat flux
In general, pool boiling is defined as boiling from a heated of 102 W/cm2, any small increment of the input power causes
surface submerged in a large volume of stationary liquid. a sudden jump to the heating element's surface temperature to
Subsequently, the immersed heating technique resembles a the film boiling region, and the transition boiling region is
pool boiling process. Theoretically, the pool boiling process missed.
passes through four stages which are natural convection, Consequently, the transition boiling's temperature range is
nucleate boiling, transition boiling, and film boiling, respec- missed in the immersed heating process. Within the film
tively. However, practically, the temperature range within the boiling stage, the heating element is covered by a stable vapour
transition boiling stage will be missed as illustrated with the aid film through which the heat energy is transferred by radiation
of Figure 3. For more details, the natural convection starts heat transfer. Accordingly, the heating element's surface tem-
when the input power is below approximately 290 W or the perature would significantly increase at a temperature normally
surface heat flux is below 25 W/cm2, where the heating ele- higher than 1000ºC as the input power (surface heat flux) in-
ment's temperature is below the mineral oil's boiling temper- creases. Note that based on many previous investigations, this
ature which is about 290ºC. In this stage, the heat energy is extremely high temperature may lead to inhomogeneous tem-
dissipated by the oil's light natural convection and as the input perature distribution along the heating element, and the hot-
power increases the heating element's temperature also in- spot temperature may reach the burnout temperature. In this
creases in a nearly linear way. Once the input power reaches way, 330ºC (Nukiyama temperature) is the highest stable
above 290 W, the nucleate boiling stage is initialised where the temperature that can be achieved by the immersed heating
heating element's temperature becomes greater than the oil's element in mineral oils which is located within the temperature
boiling temperature. Certainly, this difference in temperatures range of the so‐called T2 thermal fault (300–700ºC). While the
temperature range between the Nukiyama temperature (330ºC) under test. So, when the oil temperature approached 110°C,
and the burnout temperature (normally higher than 1000ºC) the power supply was disconnected until the oil was cooled
could not be reached, hence the T3 thermal fault cannot be down to room temperature (nearly 25°C), and then the heating
generated in the laboratory. In other words, heating beyond could be restarted.
1200 W is less controllable as compared to the situation below This was repeated until the total time of heating reaches
1200 W. In this regard, the experimental setup was utilised to the designated thermal fault duration. Note that after each
generate specifically the T2 thermal fault. heating period, the test cell was placed in cold water keeping its
upper part above the water level beside a fan directed towards
it to decrease the cooling period as much as possible. More-
2.2.1 | Test cell specifications over, after each cooling period and before the start of the next
heating period, the plunger of the syringe was compressed
It is obvious from the scheme of the thermal fault experi- gently to ensure that most of the formed gases were dissolved
mental setup depicted in Figure 2 that the test cell is composed into the oil. Regarding the above condition, the designated
of two parts. The first part which is the upper cap is a PA6 heating period for each oil sample was subdivided into small
Nylon lid consisting of the following: the power supply con- heating periods. For more details, in each small heating period,
nections and copper clamps that hold a heating element in a the heating element was powered for about 3 min followed by
70‐turn coil shape (the coil size is 5 mm in diameter, 70 mm in a 30‐min cooling period. In this investigation, the three oil
length, and made from Kanthal A‐1 resistive wire of 0.6 mm samples S5, S6, and S7 were subjected to localised overheating
diameter). Moreover, there was a terminal for the thermom- of 40, 60, and 90 min total heating durations, respectively, to
eter, a type K thermocouple with a temperature range from comprehensively study the behaviour of mineral oil when
−200°C to 1260°C, and a reverse funnel shape of 20° slope subjected to this fault type, and bearing in mind that the pre-
angle, which is already excavated in the cap itself and con- viously mentioned heating durations were achieved after 14,
nected to a syringe of 100 mL via an air bleeding pipe with a 20, and 30 cycles of heating and cooling, respectively.
three‐way valve. The second part (lower part) consists of the
following: a PA6 Nylon foundation upon a transparent Perspex
vessel of 100 mm inner diameter, 5 mm thickness, 70 mm 3 | CHARACTERISATION OF OIL
height, and a total capacity of 550 mL. Also, there are two O‐ SAMPLES
Rings for preventing the evolved gases from escaping out of
the vessel and the outside gases from entering. The O‐Rings After the degradation processes of the oil samples whose
are placed as follows: one of them is attached to the down- conditions are summarised in Table 1, test specimens were
side of the top cap to block the gap between the cap and the extracted from the test cells via a 100 mL glass gas‐tight sy-
vessel, whereas the other one is attached to the downside of ringe according to the ASTM D3613‐98 after the test cells
the vessel to block the gap between the vessel and the Nylon remained undisturbed for about 16 h. This was done to allow
foundation [28, 32]. for complete mixing and obtain a homogeneous sample be-
sides allowing the fault gases in the system to reach a quasi‐
equilibrium state [33]. For more details, the procedure
2.2.2 | Experiment procedure considered in taking the test specimens will be illustrated in the
next sentences. First of all, the upper cover of the test cell was
In the beginning, the vessel was well cleaned, then 400 mL of opened and then the syringe was positioned in the middle of
pure mineral oil was poured slowly into the test cell. Subse- the Perspex vessel to ensure that the test specimen obtained is
quently, the test cell was well sealed by bolts and nuts. After representative of the bulk from which it was taken. Following
that, the test was conducted to generate the T2 thermal fault as
seen in Figure 2. Firstly, the thermal fault was generated by
injecting controlled input power into the heating element TABLE 1 The degradation condition of each oil sample.
during the heating period ranging from 295 to 305 W. It has Number of
been found that this range of power lets the heating element Sample impulses/
reach the onset of the saturation region in the input power‐ name Fault type heating duration
heating element temperature relationship shown in Figure 3 S1 Normal case Not degraded
[30]. The previously mentioned range of power has been
attained by approximately adjusting the voltage of a 50‐Hz S2 Electrical 30 impulses
fault
single‐phase variable AC power supply (autotransformer) in S3 45 impulses
the range from 41.6 to 42.25 V which has let a passage of
S4 90 impulses
current through the heating element nearly from 7.11 to 7.23
A, respectively. Secondly, for keeping the vessel at a safe S5 Thermal fault 40 min heating
temperature away from its melting point since the oil was in S6 60 min heating
direct contact with the vessel's inner walls, a thermometer was
S7 90 min heating
utilised to monitor the surface temperature of the oil sample
23977264, 2024, 2, Downloaded from https://ietresearch.onlinelibrary.wiley.com/doi/10.1049/hve2.12405 by Cochrane Peru, Wiley Online Library on [29/06/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
324
- DARWISH ET AL.
that, the piston was pushed back slowly to pull out only 10 mL be divided into three main subregions which are as follows: the
of oil to avoid air bubbles being formed. near‐IR (NIR) region, the mid‐IR (MIR) region, and the far‐IR
Subsequently, the syringe was kept vertical and then the (FIR) region. In order to conduct the FTIR qualitative mea-
piston was carefully depressed to eject any air bubbles present surements for comprehensively studying the spectral behaviour
within the oil. Indeed, the previous steps were repeated until of the oil samples in the MIR range (from 2.5 μm (4000 cm−1)
the syringe was filled to its full capacity. Afterwards, the syringe to 50 μm (200 cm−1)) and comparing them together, a con-
was sealed with a Luer lock fitting immediately to minimise its stant volume of 20 μL was taken from each sample [36]. Then,
exposure to the atmosphere and protect the test specimen the characterisation process was carried out via the Agilent
from environmental contaminants. Actually, for each oil Cary 630 FTIR with its diamond attenuated total reflectance
degradation condition, there were two oil specimens preserved (ATR) according to the ASTM E1252‐98 standard [23]. The
in two separate syringes, one syringe was kept for the DGA reason behind choosing the MIR range specifically among the
test, while the other one was utilised to fill glass test tubes three main subregions was as follows. Most of the molecules at
fitted with screw caps to be utilised in the FTIR and DDF room temperature are in the vibrational ground state, and the
tests. Thereafter, all the test specimens were properly labelled transition from this vibrational state to the next vibrational
and stored in a cardboard box to be protected from sunlight to state is known as the fundamental transition (mode). The
avoid photodegradation. fundamental transition is of great importance for IR spec-
In this way, the specimens were ready to be examined as troscopy since it is the highest in intensity. In fact, such a
soon as possible via the following diagnostic techniques. transition within the molecules of mineral oils is always
Firstly, the chemical structure of the oil specimens was char- monitored in the MIR range [37].
acterised via FTIR spectroscopy according to ASTM E1252‐98
[23]. This characterisation was conducted to observe the
changes that occurred in the chemical structure of each oil 3.1.1 | Principle of operation
sample to find a methodology that helps in discriminating
between transformer inception faults via this optical method. As illustrated in Figure 4a, FTIR spectroscopy is composed of
Then, to initially confirm the FTIR spectroscopy analysis, the four main parts which are a source of polychromatic IR radia-
samples were analysed through the DGA method utilising a tion, a spectral analyser (Interferometer), a sample compart-
GC according to ASTM D3612‐02 [24]. Secondly, for further ment, and lastly, a radiation detector. Indeed, the MIR source
verification, the DDF test was performed according to IEC commonly used is the Globar whose typical operating tem-
60247 for each oil sample [25]. perature is about 300 K [37]. The IR beam incoming from the
source is immediately collimated and then directed towards the
interferometer. This Michelson's interferometer consists of
3.1 | FTIR spectroscopy technique three main parts, the first part is a beam splitter which is placed
at an angle of 45° and acts as a partially reflecting mirror, it splits
FTIR spectroscopy is a non‐destructive optical instrument the whole IR radiation into two equal beams, one half is
mainly utilised for particular purposes such as measuring the transmitted to the movable mirror and the second half is
basic properties of a molecule, quantitative analysis of known transmitted to the fixed mirror. After that, the two beams are
species within the material, compound identification, and reflected to the splitter again where they are recombined and
finally, structural elucidation [34]. This optical technique de-
pends on measuring the Infrared (IR) radiation transmitted or
absorbed by a certain material, which is mainly due to the
vibrational and rotational absorption bands of its molecules, as
a function of wavelength (λ) or its reciprocal which is termed
as wavenumber (ν). Furthermore, the plot resulting from such
a measurement is so called the FTIR spectrum, and it is
commonly presented with the Transmittance (T ) of the IR
radiation on the ordinate and the wavenumber on the abscissa.
T is the ratio between the radiant power transmitted by a
certain sample (I ) and the total radiant power incident on it (I0)
as clear from Equation (1) [35].
I
T¼ ð1Þ
I0
directed towards the sample compartment. Subsequently, por- analyser according to method C (Headspace Sampling) in the
tions from certain wavelength components in the whole IR ASTM D3612‐02 standard. The process of determining
radiation are absorbed by the sample's molecules, and the rest of the concentrations of dissolved gases was performed via the
the radiation is transmitted and directed towards the detector headspace GC (DANI GC1000 DPC Digital Gas Chromato-
[38]. The methodology through which the final spectrum is graph along with the DANI HSS 86.50 static headspace
obtained from the detector will be discussed more deeply. sampler). For more details, only the concentrations of the
Firstly, when the two beams are reflected from the mirrors hydrocarbon gases were considered for each oil sample, while
and recombined again at the splitter, a phase difference is the concentration of the Hydrogen gas (H2) which is an air gas
observed between them, since the paths they travelled are not was neglected for more accurate analysis as it may have a
of the same length. However, the main reason for this phase misleading impact on the whole study [24].
difference is the movable mirror movement. If the movable
mirror's distance from the splitter is equal to that of the fixed
mirror, the two beams will be in a phase resulting in 3.3 | Dielectric dissipation factor test
constructive interference, and the light power signal detected at
the detector will be of its maximum value. On the other hand, The tan‐delta (tan δ) or the dielectric dissipation factor (DDF)
if the movable mirror is displaced by λ/4 consequently the is an important electrical dimensionless ratio through which
phase difference becomes λ/2, and destructive interference the condition of the dielectric material can be determined.
between the two beams will result and the light power signal Basically, the dielectric loss of the insulation is the power
detected will be of its minimum value. But in fact, the movable dissipated via the insulating medium once it is subjected to an
mirror is always in continuous motion, forming at the detector alternating voltage. Good dielectrics usually have a very low
an oscillating light power signal (Interferogram) for each DDF. On the other hand, a high DDF may indicate that there
wavelength component within the whole IR radiation [14, 35]. are problems in the insulation's structure. Anyway, the DDF
For more information, the interferogram is a plot of the light test in this investigation was performed on all of the faulted oil
power intensity versus the optical path difference between the samples according to the IEC 60,247 besides the equipment's
two mirrors which is known as retardation. Finally, the detector user guides utilising the test setup which is composed of the
receives the overall interferogram which is due to the combi- Agilent E4980A precision LCR metre along with the Agilent
nation of all of the interferograms (cosine waves) produced 16452A liquid test fixture [25, 40].
from each wavelength component. Then, the overall interfer- While taking measurements, approximately 3.4 mL from
ogram is converted into an electrical signal, and by computing the oil sample was injected into the liquid test fixture. Then,
the Fourier transform of such signal, the FTIR spectrum is the LCR metre was set to the parallel circuit mode. This is to
simply generated [39]. measure the parallel equivalent capacitance (CP) in (F) besides
measuring the parallel equivalent resistance (RP) in (Ω) for the
oil dielectric. After selecting the circuit mode, the test was
3.1.2 | Experiment procedure performed for each oil sample in the frequency ( f ) range from
20 Hz to 2 MHz to study the oil dielectric's behaviour
The characterisation has been performed as follows: firstly, a comprehensively. Finally, for each oil sample, the DDF values
volume of 20 μL from the oil sample was placed in the ATR (represented by DDDF in the formula) corresponding to the
accessory as obvious in Figure 4b. Next to that, the IR beam whole range of frequency were simply calculated with the aid
was allowed to travel in the crystal which is a transparent of Equation (2) [25].
medium of a higher refractive index than that of the oil sample
until it struck the interface between the crystal and the oil 1
DDDF ¼ ð2Þ
medium at an angle greater than the critical angle. At that 2 � π � f � RP � C P
moment, attenuated total reflectance occurred, and the IR
beam was reflected internally into the crystal (denser medium).
Despite the internal reflection being total, an electromagnetic
wave called an evanescent wave penetrated the oil medium at 4 | RESULTS AND DISCUSSIONS
the point where the IR beam contacts the interface and was
partially absorbed by the oil sample. Finally, the reflected beam 4.1 | FTIR spectroscopy
carried the absorption information of the sample and was
collected at the detector for obtaining the FTIR spectrum of 4.1.1 | Mineral oil composition
the oil under test.
In general, the mineral oil's molecule is mainly composed of
three main hydrocarbon compounds with around 25 carbon
3.2 | Dissolved gas analysis method atoms per molecule as illustrated in Figure 5 [41, 42]. The first
one is the paraffinic compound, which is an open‐chain
For confirming the reliability of FTIR spectroscopy results, the aliphatic saturated hydrocarbon compound whose carbon
mineral oil samples each of 100 mL were examined via the GC atoms are bonded together through (C–C) sigma bonds (a type
23977264, 2024, 2, Downloaded from https://ietresearch.onlinelibrary.wiley.com/doi/10.1049/hve2.12405 by Cochrane Peru, Wiley Online Library on [29/06/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
326
- DARWISH ET AL.
a very weak intensity as the previous aromatic peak. The reason matching with the progression by which the oil samples were
for that is during the refining process, the aromatic and pol- degraded in both fault types. Similarly, the same situation
yaromatic contents are reduced to nearly 5% to enhance the occurred with aromatic peaks located at the wavenumber
stability and dielectric properties of the oil as most of these 3175 cm−1 in case of the thermal fault. Indeed, this is because
compounds are detrimental [43]. The fourth peak (at in some cases where a molecule contains multiple normal vi-
1457 cm−1) is mainly due to the bending mode of the C–H bration modes, some of these modes may occur at the same
sigma bonds in the CH2 groups within the aliphatic com- wavenumber which is known as degenerate modes. Accord-
pound. Similarly, the judgement on the oil type cannot be ingly, the transmittance intensities of the individual peaks
based on this peak only, since both paraffinic and naphthenic corresponding to the stretching and bending modes of the
compounds have CH2 groups in their structure [14, 36]. bonds within the different functional groups in a certain oil
Anyway, by entering the fingerprint region, the scene became specimen under test are varied randomly from another spec-
more clear. The fifth peak noticed (at 1377 cm−1) is mainly due to imen of the same volume taken from the same oil sample
the bending mode of the C–H sigma bonds in the CH3 groups (same degradation condition). For example, if two oil speci-
within the aliphatic compound. Based on the intensity of this mens of the same volume are taken from the same oil sample
peak, which is nearly half the intensity of the previously (same degradation condition), they may not have identical
mentioned peak, a decision was taken that the utilised mineral oil spectra, but certainly, the overall intensity of aliphatic (or ar-
is paraffinic since the CH3 groups are only found in the paraffinic omatic) peaks of both specimens will be the same [37].
compounds. Bearing in mind that, if the concentration of the Hence, the variation of the intensities is not a problem
naphthenic compounds is greater, the intensity of this peak will anymore, since this variation cannot alter the total concentra-
not be found or found with very weak intensity. The last peak tion of bonds already existing within the oil sample under
which is an important peak (at 724 cm−1) represents the bending characterisation [23]. Moreover, it is not accurate to compare
mode of the C–H sigma bonds in the CH2 groups within a long the changes that occurred in the molecular structure of the oil
open chain aliphatic compound [36, 38]. Certainly, the decision based on the alteration that happened to the intensity of certain
that was taken was confirmed by this peak. individual peaks. As such, a change does not represent the
actual change that occurred in the molecular structure of the
oil and is not a condition to be compatible with the gradual
4.1.3 | Spectral response due to electrical and increase in the degradation severity. Otherwise, for an equitable
thermal faults comparison, the overall alteration that happened in all of the
individual peaks (for each sample's spectrum) must be
The spectra of all electrical fault cases along with the no‐fault considered in the comparison instead. Based on that, the
case are depicted in Figure 6b. It is obvious from the spectra overall alteration that occurred in the aliphatic, as well as the
that there is a noticeable significant random reduction in the aromatic peaks, was considered in this investigation for
aliphatic peaks (not completely compatible with the degrada- studying the alteration that happened in the molecular struc-
tion sequence), which means a reduction in the aliphatic bonds ture of the oil due to degradation.
(aliphatic compounds decreased) within the oil molecules. Fortunately, it is observed, after utilising this methodology,
Concerning the two aromatic peaks that appeared at wave- that the overall alteration in both fault types is completely
numbers 3175 cm−1 and 1607 cm−1, respectively, there is a compatible with the gradual increase in degradation degree.
continuous decrement in the peaks (completely compatible For more details about this methodology, the total trans-
with the degradation sequence) until they completely dis- mittance intensity of all individual peaks representing the
appeared in S4 [14]. aliphatic C–H bonds within each oil sample was taken as an
Regarding the degradation via the thermal fault, the spectra indication of the concentration of the aliphatic compounds
of the samples together with the no‐fault case are presented in presented. In the same way, the total transmittance intensity of
Figure 6c. It is apparent from the spectra that there is a random all individual peaks representing the aromatic C–H bonds and
slight reduction in the aliphatic peaks (not completely the aromatic C═C bonds within each oil sample was taken as
compatible with the deterioration sequence). Subsequently, the an indication of the concentration of the aromatic compounds
aromatic peaks shown at the wavenumbers 3175 cm−1 and found.
1607 cm−1, there is a visible peak increment at the wave- Furthermore, additional analysis was performed based on
number 1607 cm−1 (inversely matches with the degradation the same concept in which the total area under all of the in-
sequence). For the peaks at the wavenumber 3175 cm−1, there dividual peaks representing the aliphatic C–H bonds within
is also an increment in the peaks (inversely quasi‐compatible each oil sample was taken as an indication of the concentration
with the deterioration sequence) [36]. of the aliphatic compounds. Also, the total area under all of the
individual peaks representing the aromatic C–H bonds and the
aromatic C═C bonds within each oil sample was taken as an
4.1.4 | Analysis interpretation indication of the concentration of the aromatic compounds
found. This analysis was conducted to validate the reliability of
First of all, it is obvious from the above spectra that the analysing the variation that occurred in the FTIR spectra by
reduction in most of the aliphatic peaks is not completely comparing the overall alteration that happened in all of the
23977264, 2024, 2, Downloaded from https://ietresearch.onlinelibrary.wiley.com/doi/10.1049/hve2.12405 by Cochrane Peru, Wiley Online Library on [29/06/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
328
- DARWISH ET AL.
peaks instead of comparing the alteration that occurred in were increased in the case of thermal degradation (total
certain individual peaks. Note that the total transmittance of transmittance decreased below the normal case).
aliphatic peaks as well as aromatic peaks besides the total area Also, Figures 7e,f show the rates of increment, as well as
under such peaks is presented in Table 3. the decrement in the total transmittance. The reason behind
the opposite alteration that occurred in the aromatic com-
Analysis via total intensity method pounds will be clarified in depth next. Concerning the degra-
By just looking at Figure 7, which illustrates the influence of dation via electrical impulses, the bond breakage that was
each fault type on the total transmittance intensities of the responsible for the reduction in the aromatic compounds
aliphatic as well as the aromatic compounds, it can be noticed occurred in two stages. First, the flow of the electrical arc
from the bar graph shown in Figure 7a that the total trans- (electrons of high energy) established an interruption to the
mittances of the aliphatics in the electrical fault cases have continuous pathways of the overlapping p‐orbitals within the
much higher values than that of thermal fault cases, which aromatic ring. This interruption made the 6 pi electrons unable
means that there is a higher reduction in the aliphatic com- to delocalise over the entire conjugate system and became
pounds within the oil samples that were subjected to electrical localised electrons. Consequently, the aromatic ring converted
deterioration. Secondly, the rate by which the total trans- from a stable resonance hybrid structure with 6 C═C bonds of
mittance increases is high concerning the cases of the electrical the same length between the carbon atoms to an unstable
fault type which is clear from Figures 7b,c. It is obvious from contributing resonance structure with only 3 C═C bonds (due
those figures that the aliphatic compounds are reduced at a to the localised pi electrons) and 3 C–C bonds. At this stage,
much higher rate in the case of electrical degradation, by the probability of the ring being thermally decomposed
approximately 2.4 from the rate of reduction in the case of the became higher than before. Secondly, the high heat energy
thermal fault. In fact, this returns to the high temperature of generated during the fault was sufficient to break most of the
the electrical degradation which is greater than 1000°C [5, 22]. ring's C–C bonds (90.2 kcal/moL), C–H bonds (110.7 kcal/
This temperature was responsible for generating much heat moL), and some of the ring's C═C bonds (174.5 kcal/moL).
energy that was capable of breaking a large amount of the C–C This results in a high decrement in the aromatic compounds
bonds (90.2 kcal/mol) and C–H bonds (101.1 kcal/moL) until they completely disappear [4, 36].
resulting in a significant reduction in the aliphatic compounds For the samples degraded by overheating, the low heat
with a very high rate. energy generated during the fault has a very low probability of
Conversely, the low temperature of the thermal fault type breaking the bonds of the stable aromatic ring which are the
(approximately 320°C) was responsible for generating lower C═C bonds (174.5 kcal/moL) and the C–H bonds (110.7 kcal/
heat energy which was only capable of breaking a small amount moL). Instead, each carbon atom within most of the cyclo-
of the C–C bonds (mostly within paraffinic compounds) and hexanes (naphthenic ring with 6 carbon atoms) breaks down
very few quantities of the C–H bonds that result in a moderate one of its two aliphatic C–H bonds. Consequently, additional
reduction in the aliphatic compounds with a very low rate [36, weak pi bonds were formed between the naphthenic carbon
44]. Furthermore, the bar graph in Figure 7d represents the atoms that were already joined together with strong sigma
alteration that happened in the aromatic compounds which bonds. This results in the formation of aromatic rings that are
provides valuable information that can be utilised in differ- more stable utilising only 49.8 kcal/moL [36, 45]. In this way,
entiating between the two kinds of faults. Concerning the an increment occurred in the aromatic compounds but in very
electrical degradation, there is a reduction that occurred in few quantities, since there were only small quantities of the
the aromatic compounds (total transmittance increased above naphthenic compounds within the utilised mineral oil. For
the normal case). On the contrary, the aromatic compounds more details, Figures 8a,b show the impacts of the electrical
TABLE 3 Results of all diagnostic techniques presented as follows: FTIR spectroscopy, gas chromatography (GC), and LCR metre.
Total transmittance of
peaks Total area under peaks
Fault type Aliphatics Aromatics Aliphatics Aromatics
Electrical fault 97.2729% 91.4572% 94.3096% 88.5425%
F I G U R E 1 3 Flowchart of the methodology proposed for utilising the In future work, the reliability of other spectroscopy tech-
FTIR spectra in the discrimination between transformer inception faults.
niques such as PL spectroscopy and UV‐Vis spectroscopy in
discriminating between transformer inception faults will be
transmittance of the aliphatics is highly increased, this may considered. Besides making some important modifications in the
point out that there is an electrical fault, otherwise the fault setups utilised in conducting the inception faults to be suitable
type may be a thermal fault. for performing online monitoring via FTIR spectroscopy.
In the previous stage, the fault type is initially determined,
and to confirm such a decision, the following step should have A C K N OWL E D G E M E N TS
proceeded. If the total transmittance of the aromatics is The authors thank the Science, Technology & Innovation
decreased, the decision will go to the thermal fault or else an Funding Authority (STIFA), Egypt, under a grant (31292) for
electrical fault will be predicted. In this way, the inception fault their support.
type within a power transformer can be determined precisely
in an easy procedure and low cost unlike the traditional C O N FL I C T O F I N TE R E S T S TA TE M E N T
techniques utilised. From the authors' point of view, the The authors declare no potential conflict of interest.
proposed methodology may be more effective if online
monitoring on the site is considered instead of offline moni- DA TA AVA I L A B I L I T Y S TA TE M E N T
toring [14]. The data that support the findings of this study are available
from the corresponding author upon reasonable request.
4. Meira, M., et al.: Dissolved gas analysis in mineral oil and natural ester 24. ASTM Standard (D3612‐02): Standard Test Method for Analysis of
liquids from thermal faults. IEEE Trans. Dielectr. Electr. Insul. 28(4), Gases Dissolved in Electrical Insulating Oil by Gas Chromatography,
1317–1325 (2021) pp. 1–22. American Society for Testing and Materials (2017)
5. IEEE Standard (C57.104‐2019): IEEE guide for the interpretation of 25. IEC Standard (IEC 60247): Insulating Liquids ‐ Measurement of Rela-
gases generated in mineral oil‐immersed transformers. In: Institute of tive Permittivity, Dielectric Dissipation Factor (Tan D) and d.c. Re-
Electrical and Electronics Engineers Inc, pp. 1–98. NY (2019) sistivity, 3rd ed. The International Electrotechnical Commission, Geneva
6. Ali, M.S., et al.: Conventional methods of dissolved gas analysis using oil‐ (2004)
immersed power transformer for fault diagnosis: a review. Elec. Power 26. IEC Standard (IEC 60599): Mineral Oil‐Filled Electrical Equipment in
Syst. Res. 216, 109064 (2023) Service – Guidance on the Interpretation of Dissolved and Free Gases
7. Taha, I.B.M., Ibrahim, S., Mansour, D.A.: Power transformer fault Analysis, 3rd ed. The International Electrotechnical Commission,
diagnosis based on DGA using a convolutional neural network with Geneva (2015)
noise in measurements. IEEE Access 9, 111162–111170 (2021) 27. IEC Standard (IEC 60060‐1): High‐voltage Test Techniques ‐ Part 1:
8. Taha, I.B.M., Mansour, D.A.: Novel power transformer fault diagnosis General Definitions and Test Requirements, 3rd ed. The International
using optimized machine learning methods. Intell. Automation Soft Electrotechnical Commission, Geneva (2010)
Comput. 28(3), 739–752 (2021) 28. Maher, A., et al.: Dissolved gas analysis and dissipation factor mea-
9. Ghoneim, S., et al.: Enhancing diagnostic accuracy of transformer faults surement of mineral oil‐based nanofluids under thermal and electrical
using teaching‐learning‐based optimization. IEEE Access 9, 30817–- faults. High Volt. 8(3), 1–11 (2023)
30832 (2021) 29. IEC Standard (IEC 60060‐2): High‐voltage Test Techniques ‐ Part 2:
10. Tran, M.‐Q., et al.: Reliable deep learning and IoT‐based monitoring Measuring Systems, 3rd ed. The International Electrotechnical Com-
system for secure computer numerical control machines against cyber‐ mission, Geneva (2010)
attacks with experimental verification. IEEE Access 10, 23186–23197 30. Wang, X.F., et al.: Dissolved gas analysis of thermal faults in transformer
(2022) liquids simulated using immersed heating method. IEEE Trans. Dielectr.
11. Badawi, M. et al.: Reliable estimation for health index of transformer oil Electr. Insul. 25(5), 1749–1757 (2018)
based on novel combined predictive maintenance techniques. IEEE 31. Przybylek, P., Gielniak, J.: Analysis of gas generated in mineral oil, syn-
Access 10, 25954–25972 (2022) thetic ester, and natural ester as a consequence of thermal faults. IEEE
12. Prasojo, R.A., et al.: Precise transformer fault diagnosis via random forest Access 7, 65040–65051 (2019)
model enhanced by synthetic minority over‐sampling technique. Elec. 32. Salah‐Eldin, A.M., et al.: Dissolved gas analysis of oil‐based TiO2
Power Syst. Res. 220, 109361 (2023) nanofluids under thermal faults. In: IEEE 21st International Conference
13. Qian, Y., Zhao, Y., Wang, Q.: Application of infrared spectroscopy in on Dielectric Liquids (ICDL), pp. 1–4. Sevilla (2022)
monitoring the aging process of transformer oil. In: IEEE 7th Asia 33. ASTM Standard (D3613‐98): Standard Practice for Sampling Insulating
Conference on Power and Electrical Engineering (ACPEE), pp. 844–848. Liquids for Gas Analysis and Determination of Water Content, pp. 1–4.
Hangzhou (2022) American Society for Testing and Materials (2017)
14. Agnello, S.: Spectroscopy for Materials Characterization, 1st ed. John 34. Mehmood, M.A., et al.: Analyzing the health condition and chemical
Wiley and Sons, Inc., NJ, USA (2021) degradation in field aged transformer insulation oil using spectro-
15. Dutta, A., Karmakar, S., Kalathiripi, H.: A novel method of fault pre- scopic techniques. In: IEEE 2018 International Conference on Di-
diction in transformer oil using infrared spectroscopy. In: IEEE 4th agnostics in Electrical Engineering (Diagnostika), pp. 1–4. Pilsen
International Conference on Condition Assessment Techniques in (2018)
Electrical Systems (CATCON), pp. 1–6. Chennai (2019) 35. Mayerhöfer, T.G., Pahlow, S., Popp, J.: The bouguer‐beer‐lambert law:
16. Karmakar, S., Dutta, A., Kalathiripi, H.: Investigation of the effect of shining light on the obscure. ChemPhysChem 21(18), 2029–2046 (2020)
high voltage impulse stress on transformer oil by infrared spectroscopy. 36. Bruice, P.Y.: Organic chemistry, 8th ed. Pearson, NJ, USA (2017)
In: IEEE International Conference on High Voltage Engineering and 37. Kaufmann, E.N.: Characterization of Materials, vol. 3, John Wiley and
Technology (ICHVET), pp. 1–5. Hyderabad (2019) Sons, Inc., NJ, USA (2012)
17. Smith, G.M., Holmes‐Smith, A.S., McMeekin, S.G.: FTIR‐spectroscopic 38. Pavia, D.L., et al.: Introduction to Spectroscopy, 5th ed Cengage learning,
and electrical analysis of thermally stressed liquid dielectrics. In: IEEE CT, USA (2014)
20th International Conference on Dielectric Liquids (ICDL), pp. 1–4. 39. Harris, D.C.: Quantitative chemical Analysis, 8th ed. Macmillan, NY
Roma (2019) (2010)
18. Zandbaaf, S., Khorrami, M.R.K., Afshar, M.G.: Prediction of dielectric 40. IEEE Standard (62‐1995): IEEE Guide for Diagnostic Field Testing of
dissipation factor by ATR‐FTIR spectroscopy based on multivariate Electric Power Apparatus ‐ Part 1: Oil Filled Power Transformers,
calibration methods for transformer oil samples in power industry. Regulators, and Reactors, pp. 1–64. Institute of Electrical and Electronics
Infrared Phys. Technol. 128, 104528 (2023) Engineers Inc., NY (2005)
19. Hadjadj, Y., et al.: Assessing insulating oil degradation by means of 41. Kaplan, I.R., Rasco, J., Lu, S.‐T.: Chemical characterization of trans-
turbidity and UV/VIS spectrophotometry measurements. IEEE Trans. former mineral‐insulating oils. Environ. Forensics 11(1‐2), 117–145
Dielectr. Electr. Insul. 22(5), 2653–2660 (2015) (2010)
20. Darwish, M.M.F., et al.: A new method for estimating transformer health 42. Pagger, E.P., et al.: Biological Insulating Liquids, 1st ed. Springer Nature,
index based on ultraviolet‐visible spectroscopy. In: IEEE 23rd Interna- Cham, Switzerland (2023)
tional Middle East Power Systems Conference (MEPCON), pp. 1–5. 43. Harlow, J.H.: Electric Power Transformer Engineering. CRC press, FL,
Cairo (2022) USA (2004)
21. Yang, D., et al.: Identification of the aging stage of transformer oil‐paper 44. Moldoveanu, S.C.: Pyrolysis of Organic Molecules: Applications to
insulation via Raman spectroscopic characteristics. IEEE Trans. Dielectr. Health and Environmental Issues. Elsevier, Inc., Oxford (2018)
Electr. Insul. 27(6), 1770–1777 (2020) 45. Zhang, J., et al.: Pyrolysis kinetics and determination of organic com-
22. Darwish, M.M.F., et al.: Application of infrared spectroscopy for ponents and N‐alkanes yields of Karamay transformer oil using TG,
discrimination between electrical and thermal faults in transformer oil. FTIR and Py‐GC/MS analyses. Fuel 306, 121691 (2021)
In: IEEE 9th International Conference on Condition Monitoring and 46. Meira, M., et al.: Dissolved gas analysis differences between natural esters
Diagnosis (CMD), pp. 255–258. Kitakyushu (2022) and mineral oils used in power transformers: a review. IET Generation,
23. ASTM Standard (E1252‐98): Standard Practice for General Techniques Transm. Distribution 13(24), 5441–5448 (2019)
for Obtaining Infrared Spectra for Qualitative Analysis, pp. 1–12. 47. Boudrâa, S., et al.: Study of spectral response of transformer oil under
American Society for Testing and Materials (2021) low electrical discharge and thermal stress. In: IEEE 18th International
23977264, 2024, 2, Downloaded from https://ietresearch.onlinelibrary.wiley.com/doi/10.1049/hve2.12405 by Cochrane Peru, Wiley Online Library on [29/06/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
DARWISH ET AL.
- 335