Acceso abierto

Performance Evaluation of ST-Based Methods for Simulating and Analyzing Power Quality Disturbances

Husham I. Hussein

Husham I. Hussein

  • Electrical Engineering Department, University of MalagaMalaga, Spain
  • Electrical Power &Machines Department, College of Engineering, University of DiyalaIraq
Buscar este autor en
Sciendo|Google Scholar
Hussein, Husham I.
, 
Ahmed Alazawi

Ahmed Alazawi

  • Electrical Engineering Department, University of MalagaMalaga, Spain
  • Electrical Power &Machines Department, College of Engineering, University of DiyalaIraq
Buscar este autor en
Sciendo|Google Scholar
Alazawi, Ahmed
, 
A. Rodríguez

A. Rodríguez

Electrical Engineering Department, University of MalagaMalaga, Spain
Buscar este autor en
Sciendo|Google Scholar
Rodríguez, A.
 y 
F. Muñoz

F. Muñoz

Electrical Engineering Department, University of MalagaMalaga, Spain
Buscar este autor en
Sciendo|Google Scholar
Muñoz, F.
  
18 oct 2023
Performance Evaluation of ST-Based Methods for Simulating and Analyzing Power Quality Disturbances's Cover Image
Acerca de este artículo
Artículo anterior
Artículo siguiente
Cite
Descargar portada
Categoría del artículo: Article
Publicado en línea: 18 oct 2023
Recibido: 11 ago 2023
DOI: https://doi.org/10.2478/ijssis-2023-0011
Palabras clave
© 2023 Husham I. Hussein et al., published by Sciendo
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Figure 1:

PQD signal models covering all types of disturbance signals (stationary and nonstationary) included normal sine waves.
PQD signal models covering all types of disturbance signals (stationary and nonstationary) included normal sine waves.

Figure 2:

Flowchart of the PQD analysis and classification system.
Flowchart of the PQD analysis and classification system.

Figure 3:

Normal sine wave: (a) target disturbance – nonstationary signal (clean and stable sine wave without any disturbance) for all methods, (b) DOST, (c) DCST, and (d) DCT.
Normal sine wave: (a) target disturbance – nonstationary signal (clean and stable sine wave without any disturbance) for all methods, (b) DOST, (c) DCST, and (d) DCT.

Figure 4:

Interruption disturbance signal: (a) target disturbance – nonstationary signal (from 0.075 s to 0.2 s) for all methods. Intended coefficients of the time and frequency fluctuations, including amplitude and angle, for all methods; (b) DOST – the signal time interval corresponds to the target signal and works to accurately capture localized changes in signal characteristics; (c) DCST – some effects on the windows; and (d) DCT – much more limited frequency resolution compared with other time–frequency analytic methods but may capture a few fine-frequency details in the PQDs.
Interruption disturbance signal: (a) target disturbance – nonstationary signal (from 0.075 s to 0.2 s) for all methods. Intended coefficients of the time and frequency fluctuations, including amplitude and angle, for all methods; (b) DOST – the signal time interval corresponds to the target signal and works to accurately capture localized changes in signal characteristics; (c) DCST – some effects on the windows; and (d) DCT – much more limited frequency resolution compared with other time–frequency analytic methods but may capture a few fine-frequency details in the PQDs.

Figure 5:

Oscillatory transient disturbance signal: (a) target disturbance – nonstationary signal and its corresponding Stockwell transform-based methods, in which disturbance in all plots occurs at 0.066–0.08 s; (b) DOST – may exhibit better noise robustness than other time–frequency analytic techniques, resulting in cleaner and clearer plotting of PQDs even in the presence of noise or interference; (c) DCST – frequency magnitude increased to 400 Hz corresponding to the time disturbance interval in (b); and (d) DCT – rising multiples in these disturbances, from 0 Hz to 105 Hz, with inter-harmonic and harmonic features.
Oscillatory transient disturbance signal: (a) target disturbance – nonstationary signal and its corresponding Stockwell transform-based methods, in which disturbance in all plots occurs at 0.066–0.08 s; (b) DOST – may exhibit better noise robustness than other time–frequency analytic techniques, resulting in cleaner and clearer plotting of PQDs even in the presence of noise or interference; (c) DCST – frequency magnitude increased to 400 Hz corresponding to the time disturbance interval in (b); and (d) DCT – rising multiples in these disturbances, from 0 Hz to 105 Hz, with inter-harmonic and harmonic features.

Figure 6:

Sag disturbance signal: (a) target disturbance – stationary signal and its corresponding plots allow each method to perform excellently, with the disturbance occurring from 0.16 s to 0.23 s; the intended coefficients of both time and frequency fluctuations, including amplitude and angle, for all methods. (b) DOST – sag effects with a corresponding rate of 100%. (c) DCST – moderate resolution may not capture all of the fine details for certain PQDs with rapidly changing characteristics, leading to less clear plots. (d) DCT – operates on discrete data points and assumes the signal to be periodic; this discrete nature may lead to artifacts in the time–frequency representation, affecting the clarity of the plots.
Sag disturbance signal: (a) target disturbance – stationary signal and its corresponding plots allow each method to perform excellently, with the disturbance occurring from 0.16 s to 0.23 s; the intended coefficients of both time and frequency fluctuations, including amplitude and angle, for all methods. (b) DOST – sag effects with a corresponding rate of 100%. (c) DCST – moderate resolution may not capture all of the fine details for certain PQDs with rapidly changing characteristics, leading to less clear plots. (d) DCT – operates on discrete data points and assumes the signal to be periodic; this discrete nature may lead to artifacts in the time–frequency representation, affecting the clarity of the plots.

Figure 7:

Swell disturbance signal: (a) target disturbance – a stationary signal and represented by all methods, with disturbance occurring from 0.2 s to 0.28 s; the intended coefficients of both time and frequency fluctuations, including amplitude and angle, for all methods. (b) DOST, (c) DCST, and (d) DCT – depend on the specific characteristics of the simulated PQDs. For disturbances with complex or rapidly changing time–frequency components, the DCT may not be the most suitable method, despite its simple but precise plots.
Swell disturbance signal: (a) target disturbance – a stationary signal and represented by all methods, with disturbance occurring from 0.2 s to 0.28 s; the intended coefficients of both time and frequency fluctuations, including amplitude and angle, for all methods. (b) DOST, (c) DCST, and (d) DCT – depend on the specific characteristics of the simulated PQDs. For disturbances with complex or rapidly changing time–frequency components, the DCT may not be the most suitable method, despite its simple but precise plots.

Figure 8:

Harmonic disturbance signal: (a) target signal – a stationary signal; the intended coefficients of both time and frequency fluctuations, including amplitude and angle, for all methods. (b) DOST time–frequency representation – ability to represent stationary and nonstationary signals, contributing to the clear plotting of disturbances. (c) DCST time–frequency representation. (d) DCT – known for its energy compaction properties, in which a large portion of the signal energy is concentrated in a few DCT coefficients; this method can result in some aspects of the disturbances being well represented, while others are not.
Harmonic disturbance signal: (a) target signal – a stationary signal; the intended coefficients of both time and frequency fluctuations, including amplitude and angle, for all methods. (b) DOST time–frequency representation – ability to represent stationary and nonstationary signals, contributing to the clear plotting of disturbances. (c) DCST time–frequency representation. (d) DCT – known for its energy compaction properties, in which a large portion of the signal energy is concentrated in a few DCT coefficients; this method can result in some aspects of the disturbances being well represented, while others are not.

Figure 9:

Sag with harmonics disturbance signal: (a) target signal – a stationary signal (from 0.05 s to 0.15 s); the intended coefficients of both time and frequency fluctuation, including amplitude and angle, for all methods. (b) DOST representation – uses different windowing functions in which a selected appropriate window can enhance the clarity of the plot by minimizing spectral leakage and side lobe effects. (c) DCST representation. (d) DCT representation – the rising multiple effective in this type of disturbance (inter-harmonics and harmonics: 25, 50, 75, and 104 Hz).
Sag with harmonics disturbance signal: (a) target signal – a stationary signal (from 0.05 s to 0.15 s); the intended coefficients of both time and frequency fluctuation, including amplitude and angle, for all methods. (b) DOST representation – uses different windowing functions in which a selected appropriate window can enhance the clarity of the plot by minimizing spectral leakage and side lobe effects. (c) DCST representation. (d) DCT representation – the rising multiple effective in this type of disturbance (inter-harmonics and harmonics: 25, 50, 75, and 104 Hz).

Figure 10:

Swell with a harmonic disturbance signal: (a) target signal – a stationary signal from 0.1 s to 0.18 s; the intended coefficients of both time and frequency fluctuation, including amplitude and angle, for all methods. (b) DOST time–frequency representation. (c) DCST representation. (d) DCT representation.
Swell with a harmonic disturbance signal: (a) target signal – a stationary signal from 0.1 s to 0.18 s; the intended coefficients of both time and frequency fluctuation, including amplitude and angle, for all methods. (b) DOST time–frequency representation. (c) DCST representation. (d) DCT representation.

Figure 11:

Flicker disturbance signal: (a) original signal – a stationary signal; the intended coefficients of both time and frequency fluctuation, including amplitude and angle, for all methods. (b) DOST, (c) DCST, and (d) DCT – the rising multiple effectiveness in these types of disturbances (25, 75, and 99 Hz for inter-harmonics and harmonics).
Flicker disturbance signal: (a) original signal – a stationary signal; the intended coefficients of both time and frequency fluctuation, including amplitude and angle, for all methods. (b) DOST, (c) DCST, and (d) DCT – the rising multiple effectiveness in these types of disturbances (25, 75, and 99 Hz for inter-harmonics and harmonics).

Figure 12:

Disturbance classification-based statistical analysis (mean, variation, standard deviation, entropy, skewness, and kurtosis) of the DOST, DCST, and DCT coefficients.
Disturbance classification-based statistical analysis (mean, variation, standard deviation, entropy, skewness, and kurtosis) of the DOST, DCST, and DCT coefficients.

Figure 13:

Accuracy comparison between the neural network method and SVM, KNN, and DT.
Accuracy comparison between the neural network method and SVM, KNN, and DT.

Figure 14:

Statistical analysis features properties’ comparison (mean, variation, standard deviation, entropy, skewness, and kurtosis) of DOST, DCST, and DCT coefficients.
Statistical analysis features properties’ comparison (mean, variation, standard deviation, entropy, skewness, and kurtosis) of DOST, DCST, and DCT coefficients.

Mathematical models for PQDs and relevant parameters_

PQ disturbance Class Equations Parameters
Normal voltage C0 vt=A1±αut2−ut1sinωt v\left( t \right) = A\left\{ {1 \pm \alpha \left[ {u\left( {{t_2}} \right) - u\left( {{t_1}} \right)} \right]} \right\}{\rm{\;sin\;}}\left( {\omega t} \right) α < 0.1
Sag C1 vt=A1−αut2−ut1sinωt v\left( t \right) = A\left\{ {1 - \alpha \left[ {u\left( {{t_2}} \right) - u\left( {{t_1}} \right)} \right]} \right\}{\rm{\;sin\;}}\left( {\omega t} \right) 0.16 ≤ α ≤ 0.23
Swell C2 vt=A1+αut2−ut1sinωt v\left( t \right) = A\left\{ {1 + \alpha \left[ {u\left( {{t_2}} \right) - u\left( {{t_1}} \right)} \right]} \right\}{\rm{\;sin\;}}\left( {\omega t} \right) 0.2 ≤ α ≤ 0.28
Interruption C3 vt=A1−αut2−ut1sinωt v\left( t \right) = A\left\{ {1 - \alpha \left[ {u\left( {{t_2}} \right) - u\left( {{t_1}} \right)} \right]} \right\}{\rm{\;sin\;}}\left( {\omega t} \right) 0.075 < α ≤ 0.2
Harmonics C4 vt=Asinωt+h3sin3ωtut3h2−ut3h1+ h5sin5ωtut5h2−ut5h1+ h7sin7ωtut7h2−ut7h1 \matrix{ {v\left( t \right)} \hfill & = \hfill & {A{\rm{\;sin\;}}\left( {\omega t} \right) + {h_3}{\rm{\;sin\;}}\left( {3\omega t} \right)\left[ {u\left( {{t_{3h2}}} \right) - u\left( {{t_{3h1}}} \right)} \right]} \hfill \cr {} \hfill & {} \hfill & { +\, {h_5}{\rm{\;sin\;}}\left( {5\omega t} \right)\left[ {u\left( {{t_{5h2}}} \right) - u\left( {{t_{5h1}}} \right)} \right]} \hfill \cr {} \hfill & {} \hfill & { + \,{h_7}{\rm{\;sin\;}}\left( {7\omega t} \right)\left[ {u\left( {{t_{7h2}}} \right) - u\left( {{t_{7h1}}} \right)} \right]} \hfill \cr } 0 ″ hi ″ 0.15
Flicker C5 ft=1+αfsinβfωtsinωt {\rm{f}}\left( t \right) = \left( {1 + {\alpha _f}sin\left( {{\beta _f}\omega t} \right)} \right)sin\left( {\omega t} \right) αf = 0.1 − 0.15, βf = 5 − 7.5
Sag with harm C6 vt=A1−αutsag2−utsag1sinωt+ h3sin3ωtut3h2−ut3h1+ h5sin5ωtut5h2−ut55h1+ h7sin7ωtut7h2−ut7h1 \matrix{ {v\left( t \right)} \hfill & = \hfill & {A\left\{ {1 - \alpha \left[ {u\left( {{t_{sag2}}} \right) - u\left( {{t_{sag1}}} \right)} \right]} \right\}{\rm{\;sin\;}}\left( {\omega t} \right)} \hfill \cr {} \hfill & {} \hfill & { +\, {h_3}\sin \left( {3\omega t} \right)\left[ {u\left( {{t_{3h2}}} \right) - u\left( {{t_{3h1}}} \right)} \right]} \hfill \cr {} \hfill & {} \hfill & { +\, {h_5}\sin \left( {5\omega t} \right)\left[ {u\left( {{t_{5h2}}} \right) - u\left( {t{5_{5h1}}} \right)} \right]} \hfill \cr {} \hfill & {} \hfill & { +\, {h_7}\sin \left( {7\omega t} \right)\left[ {u\left( {{t_{7h2}}} \right) - u\left( {{t_{7h1}}} \right)} \right]} \hfill \cr } 0.05 ≤ α ≤ 0.15
Swell with harm C7 vt=A1+αut2−ut1sinωt+ h3sin3ωtut3h2−ut3h1+ h5sin5ωtut5h2−ut5h1+ h7sinωtut7h2−ut7h1 \matrix{ {v\left( t \right)} \hfill & = \hfill & {A\left\{ {1 + \alpha \left[ {u\left( {{t_2}} \right) - u\left( {{t_1}} \right)} \right]} \right\}{\rm{\;sin\;}}\left( {\omega t} \right)} \hfill \cr {} \hfill & {} \hfill & { +\, {h_3}\sin \left( {3\omega t} \right)\left[ {u\left( {{t_{3h2}}} \right) - u\left( {{t_{3h1}}} \right)} \right]} \hfill \cr {} \hfill & {} \hfill & { +\, {h_5}\sin \left( {5\omega t} \right)\left[ {u\left( {{t_{5h2}}} \right) - u\left( {{t_{5h1}}} \right)} \right]} \hfill \cr {} \hfill & {} \hfill & { + \,{h_7}\sin \left( {\omega t} \right)\left[ {u\left( {{t_{7h2}}} \right) - u\left( {{t_{7h1}}} \right)} \right]} \hfill \cr } 01. ≤ α ≤ 0.18
Oscillatory transient C8 vt=Asinωt+βexp−t−t1/τ⋅sinωnt−t1ut2−ut1 \matrix{ {v\left( t \right)} \hfill & = \hfill & {A{\rm{\;sin\;}}\left( {\omega t} \right) + \beta {\rm{\;exp\;}}\left[ { - \left( {t - {t_1}} \right)/\tau } \right]} \hfill \cr {} \hfill & {} \hfill & { \cdot \sin \left[ {{\omega _n}\left( {t - {t_1}} \right)} \right]\left[ {u\left( {{t_2}} \right) - u\left( {{t_1}} \right)} \right]} \hfill \cr } 0.1 < β < 0.8,0.066 < τ 0.008ωn = 2π fn, 500 ≤ fn ≤ 1500
Idioma:
Inglés
Calendario de la edición:
1 veces al año
Temas de la revista:
Ingeniería, Introducciones y reseñas, Ingeniería, otros
RSS Feed de revista