Next Article in Journal
A Unified Model of a Virtual Synchronous Generator for Transient Stability Analysis
Next Article in Special Issue
Simulation Environment for the Testing of Electrical Arc Fault Detection Algorithms
Previous Article in Journal
Half-Wave Phase Shift Modulation for Hybrid Modular Multilevel Converter with Wide-Range Operation
Previous Article in Special Issue
A Method for Power Flow Calculation in AC/DC Hybrid Distribution Networks Considering the Electric Energy Routers Based on an Alternating Iterative Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Harmonic State Estimation in Power Systems Using the Jaya Algorithm

by
Walace do Nascimento Sepulchro
* and
Lucas Frizera Encarnação
*
Department of Electrical Engineering, Federal University of Espírito Santo (UFES), Av. Fernando Ferrari, 514, Vitória 29075-910, Brazil
*
Authors to whom correspondence should be addressed.
Electronics 2024, 13(17), 3559; https://doi.org/10.3390/electronics13173559
Submission received: 16 August 2024 / Revised: 2 September 2024 / Accepted: 5 September 2024 / Published: 7 September 2024
(This article belongs to the Special Issue Compatibility, Power Electronics and Power Engineering)

Abstract

:
The increasing use of nonlinear loads in power systems introduces voltage and current components at non-fundamental frequencies, leading to harmonic distortion, which negatively impacts electrical and electronic devices. A common mitigation strategy involves identifying harmonic sources and installing filters nearby. However, due to the high cost of power quality (PQ) meters, comprehensive harmonic level monitoring across the entire power system is impractical. To address this, various methodologies for Harmonic State Estimation (HSE) have been developed, which estimate distortion levels on unmonitored system buses using data from a minimal set of monitored ones. Many HSE techniques rely on optimization algorithms with numerous tuning parameters, complicating their application. This paper proposes a novel methodology for fundamental frequency power flow and harmonic state estimation using the Jaya algorithm, which is characterized by fewer tuning parameters for easier adjustment. It also introduces a strategy to determine the minimal number of buses that need monitoring to achieve system observability. The methodology is validated on the IEEE-14 and IEEE-30 bus systems, demonstrating its effectiveness. The results of the proposed methodology are compared with those obtained using Evolutionary Strategies (ESs), highlighting its enhanced accuracy and computational efficiency.

1. Introduction

Over the decades, the increase in nonlinear loads, such as LED lighting systems, frequency inverters, and rectifiers, has introduced harmonic components into power systems. These components have significant negative effects on the operation of equipment supplied by waveforms distorted by harmonics [1].
The continuous evolution of power systems, along with the growing use of nonlinear loads, has highlighted the importance of harmonic distortion and its potential to affect system performance and reliability. This emphasizes the need to develop HSE methodologies to identify harmonic sources, allowing for their treatment and mitigation.
In 1970, initial efforts in the field of power system state estimation by [2,3,4] provided a fundamental framework to address these challenges, setting a benchmark for future research in HSE. The main advantage of this methodology lies in its ability to provide a comprehensive framework for system analysis, although its direct applicability to HSE was limited due to the lack of focus on harmonic distortions.
Subsequent advancements were made by [5] in 1989, who introduced a new state estimation technique specifically aimed at identifying harmonic sources, a crucial development in the field of HSE. However, the effectiveness of this approach was contingent upon the quality of the available measurement data.
In 1991, Beides and Heydt [6] expanded the scope of HSE by incorporating the Kalman filter for a dynamic estimation of harmonic states. This innovation allowed for adaptation to the dynamic variations of systems, although challenges related to computational complexity and the need for accurate modeling limited its applicability.
The exploration of the potential of neural networks by Hartana and Richards [7] in 1990 introduced an adaptive methodology for monitoring and identifying harmonic sources, noted for its flexibility. The main limitation, however, was the dependence on the quality of training data and the specific architecture of the neural network used.
Advancing HSE using the least squares method, a hybrid nonlinear method based on the least squares method and Kirchhoff’s current laws was introduced by [8]. One of the main disadvantages noted was that the accuracy of the estimate decreased with the harmonic order to be estimated.
In 1992, artificial neural networks were also used by [9] for the estimation of harmonic voltages. The results were compared with those obtained with known methods up to that time, yielding good results in terms of accuracy.
A methodology using the least squares method in conjunction with Singular Value Decomposition (SVD) was presented by Lobos et al. [10]. Despite satisfactory results, this technique requires greater computational effort and longer processing time than other techniques.
In 2004, the SVD technique was applied to the three-phase HSE of partially observable electrical systems [11]. Although the technique proved reliable and computationally stable, the SVD technique requires more time and processing cost than traditional techniques.
In 2005, a technique using neural networks for HSE from a reduced number of meters was proposed by [12]. In this technique, neural networks are initially used to provide pseudo-measurement values based on available measurements. These values are improved using state estimation with redundant measurements. The performance was satisfactory for the simulations carried out; however, obtaining good results is tied to the location of the meters.
Arruda et al. [13] advanced the use of ESs for harmonic distortion estimation, providing robustness and adaptability to the HSE process. However, this technique presents several parameters that need to be empirically adjusted. Thus, the selection and configuration of evolutionary parameters emerged as critical challenges for the effectiveness of this technique. A limitation of the presented methodology is that it only estimates the harmonic components of the network, assuming that the PQ meter data are synchronized and that the power flow for the entire electric power system (EPS) at the fundamental frequency is known.
Also in 2010, Arruda et al. [14] presented a study of HSE based on Evolutionary Strategies, however, performing a three-phase analysis of the power system. The results were compared with those obtained using the traditional Monte Carlo technique. The robustness of the proposed method was hindered by the adjustment of numerous HSE parameters for algorithm convergence, as well as the need for meter synchronization and prior knowledge of power flow in all EPS bars and phases.
Sepulchro, Encarnação, and Brunoro introduced an algorithm in 2014 [15] for estimating both harmonic distortion and power flow at the fundamental frequency, leveraging Evolutionary Strategies. This research extends the application of HSE by incorporating the estimation of fundamental power flow. However, it also encounters difficulties related to the fine-tuning of ES parameters.
In 2021, a hybrid algorithm for HSE in a transmission EPS from a reduced number of meters was proposed by [16]. A strategy called “HA” was used, which is a combination of ES algorithms with SADE. SADE is a type of modified Differential Evolution algorithm. The goal is to combine the high convergence speed of ES with the excellent accuracy of SADE. The technique was applied to the IEEE 14-bus and IEEE 57-bus systems, and the results obtained are satisfactory; however, the technique encountered a number of parameters involved in the combination of the two algorithms.
The introduction of the Jaya algorithm for solving multi-objective optimal power flow by Warid et al. in 2018 [17] brought a new perspective on the use of optimization algorithms with few parameters in various problems. This work highlighted the efficiency and simplicity of the Jaya algorithm as valuable features to tackle optimization problems in power systems, such as HSE. The simplicity of the code, underscored by its few tuning parameters, makes the Jaya algorithm a promising tool for HSE, particularly in larger electrical systems where complexity renders established methods like Evolutionary Strategies impractical.
In 2024, a technique for optimal power flow analysis to minimize power losses was proposed by [18]. This technique is based on two optimization methods, the Jaya algorithm, and the Teacher Learning Based Optimization (TLBO) algorithm. The absence of parameters in both algorithms makes the proposed technique less complex. The technique was applied to the IEEE 39-bus system. Also in 2024, Sepulchro and Encarnação extended the application of the Jaya algorithm to HSE [19]. In this paper, preliminary results using the Jaya algorithm were presented for harmonic state and power flow estimation for the IEEE 14-bus system. The main objective of this work was just to validate the Jaya algorithm on a small-scale system.
In this sense, as a continuation of the previous work, the authors expand the application of Jaya algorithm for harmonic state and power flow estimation for larger systems. The methodology specifically targets IEEE 14-bus and IEEE 30-bus systems, aiming to demonstrate the effectiveness and efficiency of the Jaya algorithm in resolving harmonic estimation challenges. Results are compared with those obtained using ESs, chosen as a benchmark due to their proven effectiveness in HSE, supported by empirically fine-tuned numerous parameters, contrasting with the fewer parameters of the Jaya methodology. This comparison will assess the algorithms’ efficacy in terms of accuracy and processing time. Furthermore, the authors also present a strategic approach for the allocation of PQ meters in power systems, based on the number of branches at each bus and their proximity, aiming to cover the fewest possible buses while maintaining system observability for power flow estimation and HSE.

2. Evolutionary Strategies

Evolutionary Strategies (ESs) are an optimization technique inspired by evolutionary principles, originally developed by Ingo Rechenberg and Hans-Paul Schwefel in the 1960s and 1970s [20,21,22,23], with subsequent refinements and developments proposed in later works [24]. These strategies were initially conceived as a way to mimic natural evolution processes in the context of solving complex engineering problems. In ESs, a population of potential solutions, represented by individuals, undergoes iterative improvement through genetic mutations and recombination. Mutations introduce variations in individual solutions by altering their parameters, while recombination allows the exchange of information between individuals to generate potentially better solutions. The fitness of each individual is evaluated based on its proximity to the optimal solution, driving the selection of individuals for the next generation. This iterative cycle continues until a predefined stopping criterion is met, such as achieving an acceptable level of solution quality or completing a maximum number of generations.
At the end of each generation, a selection process takes place where the fittest individuals—those closest to the optimal solution—are chosen to form the next generation. This cycle of mutation, recombination, and selection continues until a stopping criterion is met, which could be a predefined number of generations or the attainment of an acceptable solution.
Briefly, an evolutionary algorithm can be described as follows [25]:
  1:
t 0 ;
  2:
initialize P ( t ) ;
  3:
evaluate P ( t ) ;
  4:
while (stopping criterion) do
  5:      
P ( t ) variation P ( t ) ;
  6:      
evaluate P ( t ) ;
  7:      
Q ( t ) = f ( P ( t ) ) ;
  8:      
P ( t + 1 ) selection [ P ( t ) P ( t ) ] ;
  9:      
t t + 1 ;
10:
end while
In this context, P ( t ) represents a population with μ individuals. P ( t ) denotes a set of λ individuals generated through recombination and mutation from P ( t ) . Q ( t ) is a function of P ( t ) and can be zero or equal to P ( t ) . During the evaluation phase, each individual is scored based on its distance from the optimal solution. The best-performing individuals are then selected to form the initial population for the next generation. This iterative process continues until the stopping criterion is satisfied. The basic flowchart of the ES algorithm is shown in Figure 1.

3. Jaya Algorithm

3.1. Introduction

The Jaya algorithm, developed by [26], is a groundbreaking methodology in the field of mathematical optimization, distinguished by its simplicity and efficiency. Unlike many other optimization algorithms, Jaya requires only a few specific parameters for its operation, significantly reducing the complexity of its implementation and tuning.
The core principle of the Jaya algorithm is to steer solutions towards the best solution obtained so far, while simultaneously moving them away from the worst solution observed within the solution population. This “directional search” approach enables rapid convergence to an optimal or near-optimal solution, without the need for crossover or mutation operators typical of other optimization algorithms, such as genetic algorithms and Evolutionary Strategies.
The effectiveness and versatility of the Jaya algorithm have been demonstrated in various studies and applications, addressing complex optimization problems across different engineering fields [27]. The Jaya algorithm has shown great promise in solving power flow problems in electric power systems, particularly noted for its optimization capability without the need for tuning specific parameters. In a study conducted by [28], the Jaya algorithm proved effective in optimizing power flow, producing optimal solutions with rapid convergence, and demonstrating superior performance compared to other stochastic algorithms in terms of the optimality and feasibility of the proposed solutions. This research underscores the applicability of the Jaya algorithm in obtaining efficient and viable solutions for complex power flow problems.
A significant innovation was introduced by [17] with the development of the modified quasi-oppositional Jaya (QOMJaya) algorithm to solve multi-objective optimal power flow (MOOPF) problems. This method enhances convergence properties, exploration capabilities, and solution optimality by incorporating a learning strategy based on quasi-opposition.
Additionally, [29] proposed a hybrid Jaya–Powell’s Pattern Search approach to solve the multi-objective optimal power flow problem, integrating the Jaya algorithm for exploration and Powell’s Pattern Search method for exploitation.
However, considering the absence in the current literature of any applications of the Jaya algorithm in the HSE of electric power systems, the canonical version of Jaya as presented by Rao in 2016 will be applied.

3.2. Basic Equation

The basic update equation in the Jaya algorithm for optimizing a problem is given by [27]
X i , j ( n e w ) = X i , j + r 1 · ( X b e s t , j X i , j ) r 2 · ( X w o r s t , j X i , j )
where X i , j ( n e w ) is the new value of the j th variable for the i th candidate solution in the next iteration, X i , j is the current value of the j th variable for the i th candidate solution, X b e s t , j is the value of the j th variable for the best candidate solution in the current population, X w o r s t , j is the value of the j th variable for the worst candidate solution in the current population, and r 1 and r 2 are random numbers generated within the range [0, 1] for each iteration and each variable, providing stochastic elements to the update mechanism.
This equation is used to update each variable of each candidate solution in the population. The idea is to move the current solution towards the best solution and away from the worst solution, thereby iteratively improving the quality of solutions in the population.

3.3. Process Overview

As shown in Figure 2, the Jaya algorithm begins with initialization and initial settings, where a population of solutions is randomly generated within the search space. Following this is the evaluation step, where each solution is assessed based on the problem’s objective function, and the best and worst solutions in the population are identified, crucial for the subsequent steps. The update step uses Equation (1) to adjust each solution, propelling them towards the best solution and distancing them from the worst. This process is repeated in a loop, where after each update, the solutions are reevaluated, and the best and worst are identified again. The loop continues until a stopping criterion is satisfied, such as reaching the maximum number of iterations or a specified convergence condition [26]. The flowchart ends when the algorithm presents the best solution found as the result.

4. Problem Formulation

4.1. General Scheme

Firstly, the harmonic orders of interest are defined, and the admittance matrices for the fundamental frequency and harmonic frequencies are calculated from the power system data. Considering that the methodology proposed in this work aims to estimate power flow and HSE by assessing harmonic voltage distortion levels at buses without direct measurements, it is assumed that measurements are available only at select buses in the power system. These measured values serve as a database for estimating fundamental and harmonic frequency values at buses without direct measurements. The analysis and estimation of data for the fundamental frequency and each harmonic order will be conducted separately, employing either the Jaya algorithm or Evolutionary Strategies. Finally, by comparing the measured and estimated voltage values for the fundamental frequency and harmonic orders, the Total Harmonic Distortion (THD) can be determined. The proposed methodology is illustrated by the block diagram shown in Figure 3.

4.2. Reference Values

PQ meters installed at selected buses of the system deliver synchronized reference data for current and voltage, including magnitude and phase values, across the harmonic orders of interest as well as the fundamental frequency.

4.3. Harmonic Orders of Interest

The determination of which harmonic orders to examine is influenced by the availability of PQ meters and the specific objectives of the research. In this particular study, the emphasis is placed on identifying and analyzing the most frequently encountered harmonics in power systems. This includes the odd harmonic orders 3, 5, 7, 9, 11, and 13, which are often most significant in terms of their impact on system performance and equipment [30], enabling the calculation, from the values at the fundamental frequency, of the THD given by
THD = h = 2 n V h 2 V 1
where THD represents the Total Harmonic Distortion, V h signifies the magnitude of the h th harmonic component, V 1 denotes the magnitude of the fundamental component, and n represents the last harmonic component considered in the summation.

4.4. Individual Representation

In the HSE study, each individual represents a potential solution, embodying a possible state of the system. To address the study’s objectives effectively, it is essential to identify a distinct individual for each specific harmonic order h. This individual is responsible for representing the optimal solution for the system at that particular harmonic order. In this context, for a system consisting of t buses and m meters, and using the Jaya algorithm for optimization, each individual is represented by a matrix of ( t m ) rows and two columns.
Taking the IEEE 14-bus system as an example, the representation of an individual for Jaya is defined according to the expression in Equation (3). In this specific case, monitoring is limited to buses 2, 3, 8, 9, 12, and 14. The matrix columns, first and second, represent the magnitudes and phases of the voltages at each bus for the harmonic order h under analysis, respectively.
Individual h = V h 2 θ h 2 V h 3 θ h 3 V h 5 θ h 5 V h 7 θ h 7 V h 9 θ h 9 V h 11 θ h 11 V h 12 θ h 12 V h 13 θ h 13
where V h i denotes the absolute voltage at bus i in harmonic order h; θ h i indicates the phase of voltage at bus i in harmonic order h.
However, within the context of the ES methodology, an additional column is appended for each variable, considering that it is imperative to account for mutation steps aligned with the magnitudes and phases of voltages, thus necessitating the addition of two columns in this instance. As an illustrative example, the representation of an individual for the IEEE 14-bus network is defined as shown in (4), where measurements are taken exclusively at buses 1, 4, 6, 8, 10, and 14. The first and third columns represent the magnitudes and phases of voltage at each bus and harmonic order h, while the second and fourth columns denote the mutation steps.
Individual h = V h 2 σ V , h 2 θ h 2 σ θ , h 2 V h 3 σ V , h 3 θ h 3 σ θ , h 3 V h 5 σ V , h 5 θ h 5 σ θ , h 5 V h 7 σ V , h 7 θ h 7 σ θ , h 7 V h 9 σ V , h 9 θ h 9 σ θ , h 9 V h 11 σ V , h 11 θ h 11 σ θ , h 11 V h 12 σ V , h 12 θ h 12 σ θ , h 12 V h 13 σ V , h 13 θ h 13 σ θ , h 13 ,
where V h i denotes the absolute voltage at bus i in harmonic order h; θ h i indicates the phase of voltage at bus i in harmonic order h; σ V , h i represents the mutation step for the voltage parameter at bus i in harmonic order h (not applicable to the Jaya algorithm); and σ θ , h i signifies the mutation step for the phase parameter at bus i in harmonic order h (not applicable to the Jaya algorithm).

4.5. State Estimation for Each Harmonic of Interest

In the HSE process, following the selection of harmonic orders considered relevant for this study and the definition of individual representation, each is subjected to the state estimation algorithm proposed in this work, which is based on the Jaya algorithm or ES algorithm. The input data for the algorithm not only include the selected harmonic orders but also actual measurements of voltages and currents, coupled with the system’s admittance matrix. The execution of the algorithm for each harmonic of interest is meticulously carried out, and the results are compiled into a final report.

4.6. Evaluation

The evaluation process is designed to determine the fitness of each individual by quantifying how closely each solution approximates the established reference values. In this context, each individual is considered a potential solution to the HSE problem. Thus, for each harmonic order h selected for analysis, the voltage vector corresponding to each individual, combined with the system’s admittance matrix, results in the formulation of an estimated current vector I h i , as illustrated by
I h i = j = 1 n b Y h i j V h j
where V h j is the voltage at bus j, harmonic order h, Y h i j is the element ( i j ) of the admittance matrix at the frequency for harmonic order h, and n b is the number of monitored buses.
The accuracy of these estimates is verified by calculating the difference between the estimated current vector and the measured current vector I h M i , a process that generates the error vector e e i h , as shown by
e e h i = | I h i I h M i |
where e e i h is the estimation error at bus i, harmonic order h, I h i is the calculated current at bus i from an individual who represents the voltages at the buses for selected harmonic orders, and I h M i is the measured current, harmonic order h at bus i. This error vector quantitatively represents the discrepancy between the estimated condition and the actual system condition, providing a crucial metric for performance evaluation. After calculating the errors for all monitored buses, the fitness of each individual analyzed for the specified harmonic order is calculated as the inverse of the sum of the squares of the error vector elements, as defined by
F i t n e s s h = 1 i = 1 n b ( e e h i ) 2
where Fitness h is the individual fitness at harmonic order h.

4.7. Selection

In the case of the ES methodology, the selection is deterministic, ensuring that only the fittest individuals, i.e., those with the highest fitness, are chosen while maintaining a constant population size from generation to generation. Selection includes both parents and offspring to allow for elitism, ensuring that the best parent individuals are preserved throughout generations.
On the other hand, the Jaya algorithm methodology employs a unique selection methodology within its iterative process, as detailed in Figure 2. Originating from an initial population generated randomly within the search space, each individual in the population automatically generates a successor using the Equation (1). This successor then undergoes a fitness evaluation, where its fitness is compared to that of its progenitor. If the successor’s fitness surpasses that of its originator, it replaces the latter in the population; otherwise, it is discarded, thus optimizing the process by focusing resources on more promising candidates.
Furthermore, the Jaya algorithm identifies the best and worst individuals within the population throughout all iterations. The best individuals are those who consistently demonstrate the highest levels of fitness, effectively guiding the search toward more optimal solutions. On the other hand, the worst individuals represent the lowest levels of fitness, serving as benchmarks for the necessary improvements and adjustments in the algorithm’s search dynamics. Once the stopping criterion is met, the best solution for each harmonic order is considered the final solution for that specific harmonic order.

5. Application

5.1. Initial Considerations

In this work, the IEEE 14-bus and IEEE 30-bus systems were selected for the application of the proposed HSE methodology and are illustrated in Figure 4 and Figure 5, whose parameters and other information are obtained from [31,32], respectively. More informations can be found in Appendix A. The choice of these specific systems stems from their established usage in power system studies, offering a balance between complexity and manageability. Additionally, the availability of test data and widespread use in academic research make them suitable candidates for validating the proposed methodology.
Initially, both systems were simulated using PSCAD X4 software (version 4.5.0.0) at the fundamental frequency to validate the data. Since the operating conditions at the fundamental frequency considered in this study are the same as those of the official test models, the power flow results could be validated through [33,34] for the IEEE 14-bus and IEEE 30-bus systems, respectively.
Subsequently, both systems were simulated using the PSCAD software in a scenario featuring multiple harmonic sources of orders 3, 5, 7, 9, 11, and 13, which were randomly distributed.
Figure 4. IEEE 14-bus system (adapted from [31]).
Figure 4. IEEE 14-bus system (adapted from [31]).
Electronics 13 03559 g004
The power flow at the fundamental frequency and the harmonic state of the studied systems at all buses, including the values of current and voltage harmonic orders of interest, were obtained from PSCAD software and are provided in Appendix B.
Figure 5. IEEE 30-bus system (adapted from [35]).
Figure 5. IEEE 30-bus system (adapted from [35]).
Electronics 13 03559 g005

5.2. PQ Meters Allocation

The initial and fundamental step in the HSE methodology is to determine the minimum number and optimal locations of buses that need to be monitored to ensure system observability.
One of the selection criteria for measurement points involves considering the number of line segments connected to each bus. Buses with a higher number of connections are prioritized for measurement inclusion. Additionally, it is crucial to avoid selecting adjacent buses. This approach is based on the premise that voltages at neighboring buses can be estimated if one of them is already measured, thus minimizing the total number of PQ meters required.
To effectively select these buses, a comprehensive strategy was developed. This strategy, illustrated in the flowchart in Figure 6, involves several steps. First, all buses are mapped and the number of line segments branching from each bus is counted. Then, the buses are listed in descending order according to the number of branches. This hierarchical listing allows for prioritizing buses with the highest connectivity for measurement inclusion.
The strategy also incorporates a method to systematically check for adjacency, ensuring that directly connected neighboring buses are not selected consecutively. This step is crucial for maintaining the efficiency of the measurement system and reducing redundancy.
Moreover, the flowchart incorporates an iterative process to re-evaluate and adjust the selection as required. This ensures that the final group of selected buses offers comprehensive observability of the system with minimal PQ meters. Implementing this approach enhances the efficiency and effectiveness of the HSE process, resulting in improved accuracy and reliability in HSE.
The mapping results for the IEEE 14-bus system are presented in Table 1. According to the adopted strategy, monitoring is required for 6 specific buses: buses 1, 4, 6, 8, 10, and 14.
For the IEEE 30-bus system, the mapping yielded Table 2, and the adopted strategy resulted in the selection of 12 buses for comprehensive monitoring. Specifically, the chosen buses are buses 3, 5, 6, 11, 12, 17, 18, 20, 21, 24, 26, and 27. This strategic approach ensures extensive and effective coverage, enabling detailed and accurate observation of harmonic states across different sections of the system.

5.3. ES Methodology Tuning

Table 3 presents the various parameters considered characteristic of an algorithm based on ES. These values were empirically obtained, aiming for convergence in the shortest possible time and accuracy of estimated values.

5.4. Jaya Algorithm Tuning

Given that the Jaya algorithm in its original form includes only two parameters—the number of individuals in the initial population and the number of generations—these parameters were selected empirically. For the IEEE 14-bus system, the number of individuals in the initial population was set at 50, and the maximum number of generations was set at 2000. However, for the IEEE 30-bus system, these parameters were set at 150 and 5000, respectively, as listed in Table 4.
The stopping criterion was defined as either reaching the maximum number of generations or the stabilization of the fitness of the best individual, whichever occurs first. In this study, stabilization is determined when the fitness of the best individual remains unchanged for 10% of the maximum number of generations.

5.5. Implementation

The methodologies were implemented using the Python programming language, incorporating specific routines for computing the admittance matrix across various harmonic orders, reading reference values, and performing HSE using ES methodology and the Jaya algorithm.

6. Results and Discussion

Considering the premises previously presented, the proposed methodology was applied to the IEEE 14-bus and IEEE 30-bus systems using the proposed Jaya algorithm as well as the ES algorithm for a comparative study between the two methodologies.
Thus, the analysis involved the introduction of harmonic currents of order 3, 5, 7, 9, 11, and 13 in some buses, causing harmonic distortion in the system. This scenario presents a challenge for the system operator to identify these sources, considering a hypothetical situation where the operator is unaware of these sources and does not have PQ meters throughout the system.
As discussed, the ES algorithm required empirical tuning of various input parameters, which became a daunting task as the parameters are sensitive to each other. This necessitated running the codes multiple times to find suitable parameter tuning, a situation that is exacerbated when dealing with different problems, requiring specific parameters tuned for each analysis situation. This tuning is performed to allow both the diversity of the initial population of possible solutions and the rapid convergence to a global optimum. Conversely, the Jaya algorithm does not allow parameter adjustment freedom, containing only two parameters, namely the initial population size and the maximum number of generations, making it very simple and easy to tune.
With the initial parameters defined, based on voltage measurements at some selected buses, the algorithm receives system data, only the measured values at these buses and the initial parameters, and estimates the voltage amplitudes and phases at unmonitored buses for the fundamental frequency and the harmonic orders of interest.
To demonstrate that the voltage estimation errors are consistently kept within a minimal value range for each harmonic order, the proposed HSE methodology was executed 30 times, first applying the ES algorithm and then the Jaya algorithm, obtaining the results for each fundamental and harmonic order of interest at unmonitored buses. This allows us to obtain the THD and compare the time and consequently the processing cost for both algorithms. The results are presented and discussed below.
Considering that the amplitudes of the harmonic orders can be reasonably small and that any slight difference would result in considerable relative errors, the approach was to address the absolute errors between the estimated and real values, which allows for better analysis of the results. The average absolute errors of the obtained voltage magnitudes are presented in Table 5 and Table 6 using the ES and Jaya methodologies, respectively, for the IEEE 14-bus system, and in Table 7 and Table 8 using the ES and Jaya methodologies, respectively, for the IEEE 30-bus system. The same results are presented graphically in Figure 7 and Figure 8 for the IEEE 14-bus system and in Figure 9 and Figure 10 for the IEEE 30-bus system. Considering the estimates for the fundamental frequency, it is observed that the highest errors obtained using ES were 0.051 p.u. and 0.043 p.u. for the IEEE 14-bus and IEEE 30-bus systems, respectively, higher than the errors obtained with the Jaya algorithm, which were 0.003 p.u. and 0.011 p.u. for the two systems, respectively, demonstrating the efficiency of the Jaya algorithm for power flow estimation. For the harmonic components, considering the IEEE 14-bus system, the highest errors obtained are 0.010 p.u. using ES and 0.002 p.u. using Jaya, both in the third harmonic. In the case of the IEEE 30-bus system, the highest errors are 0.006 p.u. using ES and 0.007 p.u. using Jaya, both in the eleventh harmonic, which may be an insignificant difference with little impact on the THD estimation, as it is a higher-order harmonic whose magnitudes, in this case, are smaller than those of lower-order harmonics. Thus, the errors obtained with the Jaya algorithm, which are comparable or lower to those obtained with the ES algorithm for the harmonic components, demonstrate the effectiveness of the HSE using the proposed methodology based on the Jaya algorithm.
Considering the intrinsic characteristic of an evolutionary optimization algorithm like ES, which is based on selecting the fittest individuals to compose the initial population of the next generation, and an algorithm with few parameters like Jaya, which is directed by the best individual of each generation, the success in HSE is verified by the convergence of individuals or possible solutions to the optimal solution, that is, to the solution with the highest fitness and consequently the lowest error. In the present study, the lower error implies estimating the fundamental and harmonic components closer to the real values. To illustrate the error evolution curves’ behavior throughout the iterative process, Figure 11 and Figure 12 show the error evolution curves for the IEEE 14-bus network using, respectively, the ES and Jaya methodologies, and Figure 13 and Figure 14 show the error evolution curves for the IEEE 30-bus network using, respectively, the ES and Jaya methodologies. It is noted that the curves are raised for each harmonic order, considering that the algorithm is applied to each harmonic order at a time. The analysis of these curves allows us to conclude that they all move towards the lowest error plateau, stabilizing at certain values, proving the convergence of both methodologies. For harmonic frequencies, errors ranging from 0.000 to 0.061 p.u. for the IEEE 14-bus system and from 0.000 to 0.081 p.u. for the IEEE 30-bus system are observed. Some of the most significant error values are related to higher- and lower-order harmonics, whose impact on the total THD calculation is practically negligible. These errors for harmonic frequencies demonstrate that the algorithm is also effective in HSE.
From the estimated voltage values for the fundamental frequency and harmonic orders at the unmonitored buses, the THD values at these buses are determined using Equation (2). These values are listed in Table 9 and Table 10 using the ES and Jaya methodologies, respectively, for the IEEE 14-bus system, and in Table 11 and Table 12 using the ES and Jaya methodologies, respectively, for the IEEE 30-bus system. The obtained THD error values for both methodologies are also illustrated in Figure 15 and Figure 16. It is observed that the THD errors for the IEEE 14-bus system are slightly higher for the ES algorithm, with a maximum error of 0.990% compared to 0.190% for the Jaya algorithm. However, for the IEEE 30-bus system, the Jaya algorithm resulted in a maximum error of 0.609%, which is slightly higher than the maximum error of 0.446% obtained with the ES algorithm, although the average error values are within the same range. Finally, the comparative analysis of these THD results allows for the conclusion that the Jaya algorithm is as effective as the ES algorithm by presenting average and maximum THD errors that are comparable to or better than those resulting from the ES algorithm. These errors are irrelevant for HSE, considering that the obtained results allow a comprehensive analysis of the systems for identifying the buses with the highest and lowest THD and consequently allowing the location of harmonic sources.
Finally, considering the main characteristic of the Jaya algorithm compared to other methodologies existing in the literature, such as the ES algorithm, which has a simpler code structure and a minimal number of parameters, it is worth addressing the average processing times recorded for both methodologies when applied to the IEEE 14-bus system, shown in Figure 17, as well as when applied to the IEEE 30-bus system, presented in Figure 18. It is observed that in both systems and for all harmonic orders, the Jaya algorithm exhibited considerably lower processing times, finding results of practically the same quality as the ES algorithm in up to 50% less time. This proves that the simplicity of the code and the minimal number of parameters enable a quick solution using the Jaya algorithm.
Figure 19 shows the total average time to execute the proposed methodology in this work during the 30 simulations performed. This methodology, which estimates the power flow at the fundamental frequency and all harmonic orders of interest, is applied to both the IEEE 14-bus and IEEE 30-bus systems using the ES and Jaya algorithms, as illustrated in Figure 3. It is observed that with the Jaya algorithm and considering the IEEE 14-bus system, the proposed methodology takes about half the time to obtain the results compared to the ES algorithm. For the IEEE 30-bus system, the Jaya algorithm achieves the desired results in approximately 64% of the total time required by the ES algorithm. Even though the convergence of Jaya may be slower than the ES algorithm, meaning it may take a larger number of generations to converge, Jaya has a lower processing time, which can be attributed to its simplicity and fewer parameters. Thus, corroborating the objective proposed in this work, the Jaya algorithm is significantly faster in achieving equally satisfactory results as the ES algorithm, resulting in lower processing costs as a significant advantage.

7. Conclusions

The application of the proposed power flow and HSE methodology based on the Jaya algorithm to the IEEE 14-bus and IEEE 30-bus systems, especially at buses not equipped with PQ meters, using the measured values from other buses, resulted in estimated voltage values with minimal error or little impact on the THD calculation. This led to the determination of more realistic THD values, demonstrating the effectiveness of the proposed methodology based on the Jaya algorithm compared to the results obtained with the ES algorithm.
The adopted strategy for determining the minimum number of meters and selecting the buses to receive these meters proved to be accurate, as the results achieved with both the Jaya and ES algorithms demonstrated the effectiveness of the proposed methodology in estimating the fundamental and harmonic voltages at the other buses from the indicated buses, with minimal absolute errors.
In general, the implementation of the power flow and HSE methodology, based on the Jaya algorithm, proved to be effective in achieving satisfactory results in less time and with lower processing costs compared to those obtained with the ES algorithm. The dynamics of the Jaya algorithm, which aims to move candidate solutions away from the worst solutions and closer to promising solutions, demonstrated robustness in achieving convergence to solutions near to the reference values, allowing for a more precise calculation of THD at the unmonitored buses.
Even though it presented slower convergence compared to the ES algorithm, the absence of parameters and the simplicity of the Jaya algorithm provided a quicker response in the iterative optimization process, without the need for costly empirical parameter adjustments or a suitable set of parameters for each scenario or power system under study, making the proposed methodology applicable to various scenarios with minimal changes in the number of individuals in the initial population and the maximum number of generations.
Based on the obtained results, the presented methodology based on the Jaya algorithm emerges as a promising strategy of HSE, using a minimum number of meters allocated at specific buses through the outlined strategy in this study. This opens avenues for its application in larger, real-world systems, likely yielding equally satisfactory outcomes, through an algorithm that requires few tuning parameters. Such implementation could result in cost savings in studying and identifying harmonic sources in power systems, by reducing the required quantity of PQ meters.

Author Contributions

Conceptualization, W.d.N.S. and L.F.E.; methodology, W.d.N.S.; software, W.d.N.S.; validation, W.d.N.S.; formal analysis, W.d.N.S. and L.F.E.; investigation, W.d.N.S.; resources, L.F.E.; writing—original draft preparation, W.d.N.S.; writing—review and editing, W.d.N.S. and L.F.E.; visualization, W.d.N.S.; supervision, L.F.E.; project administration, L.F.E.; funding acquisition, L.F.E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Council for Scientific and Technological Development—CNPq (grant numbers 404857/2023-0 and 311848/2021-4) and Espírito Santo Research and Innovation Support Foundation—FAPES (grant number 1024/2022).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
EPSelectric power system
HSEharmonic state estimation
IEEEInstitute of Electrical and Electronics Engineers
PQpower quality
THDTotal Harmonic Distortion

Appendix A

Appendix A.1

Table A1 and Table A2 reproduce the original data of the IEEE 14-bus system, which were used in modeling the system in the electromagnetic transient software PSCAD X4.
Table A1. Line data—IEEE 14-bus system.
Table A1. Line data—IEEE 14-bus system.
Line NumberFrom BusTo BusResistance
(p.u.)
Reactance
(p.u.)
Half Line Charging
Susceptance (p.u.)
MVA
Rating
1120.019380.059170.02640120
2150.054030.223040.0219065
3230.046990.197970.0187036
4240.058110.176320.0246065
5250.056950.173880.0170050
6340.067010.171030.0173065
7450.013350.042110.0064045
84700.20912055
94900.55618032
105600.25202045
116110.094980.1989018
126120.122910.25581032
136130.066150.13027032
147800.17615032
157900.11001032
169100.031810.0845032
179140.127110.27038032
1810110.082050.19207012
1912130.220920.19988012
2013140.170930.34802012
Table A2. Bus data—IEEE 14 bus system.
Table A2. Bus data—IEEE 14 bus system.
Bus NumberBus VoltageGenerationLoadReactive Power Limits
Magnitude (p.u.)Phase Angle (Degree)Real Power (MW)Reactive Power (MVAR)Real Power (MW)Reactive Power (MVAR)Qmin (MVAR)Qmax (MVAR))
11.0600114.17−16.900010
21.045040.00021.712.7−42.050.0
31.01000094.219.123.440.0
4100047.8−3.9
510007.61.6
6100011.27.5
7100000
8100000
9100029.516.6
1010009.05.8
1110003.51.8
1210006.11.6
13100013.85.8
14100014.95.0

Appendix A.2

Table A3 and Table A4 reproduce the original data of the IEEE 30-bus system, which were used in modeling the system in the electromagnetic transient software PSCAD X4.
Table A3. Line data—IEEE 30-bus system.
Table A3. Line data—IEEE 30-bus system.
Line NumberFrom BusTo BusResistance
(p.u.)
Reactance
(p.u.)
Half Line Charging
Susceptance (p.u.)
MVA
Rating
1120.020.060.03130
2130.050.200.02130
3240.060.180.0265
4250.050.020130
5260.060.180.0265
6340.010.040130
7460.010.04090
841200.23065
9570.050.120.0170
10670.030.080130
11680.010.09032
126900.21065
1361000.56032
146280.070.060.0132
158280.060.200.0232
1691100.21065
1791000.11065
1810200.090.21032
1910170.030.09032
2010210.030.08032
2110220.070.15032
22121300.14065
2312140.120.26032
2412150.070.13032
2512160.010.12032
2614150.220.12016
2715180.110.22016
2815230.100.21016
2916170.080.19016
3018190.060.13016
3119200.030.07032
3221220.010.22032
3322240.110.18016
3524250.190.33016
3625260.250.38016
3725270.110.21016
3827290.220.4016
3927300.320.60016
40282700.4065
4129300.240.45016
Table A4. Bus data—IEEE 30 bus system.
Table A4. Bus data—IEEE 30 bus system.
Bus NumberBus VoltageGenerationLoadReactive Power Limits
Magnitude (p.u.)Phase Angle (Degree)Real Power (MW)Reactive Power (MVAR)Real Power (MW)Reactive Power (MVAR)Qmin (MVAR)Qmax (MVAR))
1100024.963−4.638−20150
21021.712.760.9727.677−2060
3102.41.20000
4107.61.60000
510000000
610000000
71022.810.90000
81030300000
910000000
10105.91920000
1110000000
121011.27.50000
1310003713.949−1544.7
14106.21.60000
15108.22.50000
16103.51.80000
171095.80000
18103.20.90000
19109.53.40000
20102.20.70000
211019.66911.200000
22100031.5940.34−1562.5
23103.21.622.28.13−1040
2410156.700000
25101.000.000000
26103.502.300000
27100028.9110.97−1548.7
2810000000
29103.6590.900000
301012.001.900000

Appendix B

Appendix B.1

Table A5, Table A6, Table A7, Table A8, Table A9 and Table A10 present the fundamental and harmonic component values of current and voltage derived from the power flow and harmonic state analysis for the IEEE 14-bus system using the electromagnetic transient software PSCAD X4. These values are based on the current injection harmonics in randomly selected buses.
Table A5. Reference values of magnitude in p.u. and phase in degrees for the fundamental current and the 3rd- and 5th-order current harmonics for the IEEE 14-bus system.
Table A5. Reference values of magnitude in p.u. and phase in degrees for the fundamental current and the 3rd- and 5th-order current harmonics for the IEEE 14-bus system.
BusI1hI3hI5h
MagnitudePhaseMagnitudePhaseMagnitudePhase
11.35680.84690.0675114.00840.01041.7873
21.061025.72890.076952.35770.01351.7903
30.3546−77.94380.0303124.66840.003814.0743
43.0349114.73110.1398131.07760.0539−5.4048
52.3628−113.96610.2465−71.81840.0822178.0589
61.0942−154.48070.0707−69.58130.0203179.2846
71.1543−161.33940.0645−96.90260.0129161.8508
80.4530−2.66870.024077.42070.0040−15.1816
91.254012.05370.061874.99230.0115−16.0899
101.3448−165.29710.0381−99.66080.0077167.0406
111.553410.31820.0320110.88680.0068−0.0682
120.138423.28140.0239102.98300.0126−1.8550
130.278530.95530.0490113.07340.0102−1.9336
140.2758−176.68540.0150−59.11080.0015165.4553
Table A6. Reference values of magnitude in p.u. and phase in degrees for the fundamental voltage and the 3rd- and 5th-order voltage harmonics for the IEEE 14-bus system.
Table A6. Reference values of magnitude in p.u. and phase in degrees for the fundamental voltage and the 3rd- and 5th-order voltage harmonics for the IEEE 14-bus system.
BusV1hV3hV5h
MagnitudePhaseMagnitudePhaseMagnitudePhase
11.03890.00000.019649.73810.000754.5349
21.0249−2.23240.024151.65910.000954.6651
31.0474−5.67070.018133.09860.002284.7141
40.9918−9.36220.028623.36990.00112.5361
51.0733−8.43360.047527.22490.0116−79.7017
61.0475−8.43120.037227.22480.0046−79.0903
71.0357−9.41910.035624.60040.0029−84.2933
81.0139−5.12640.024837.49280.000951.8321
91.0079−6.55820.027333.54270.000735.1230
101.0297−9.05000.032224.23130.0020−84.0741
111.03410.48520.025226.69730.000221.3587
121.0334−6.73450.020435.31040.007187.7918
131.0277−7.03530.021132.27010.003486.1968
141.0142−9.38670.031533.04500.000668.3699
Table A7. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order current harmonics for the IEEE 14-bus system.
Table A7. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order current harmonics for the IEEE 14-bus system.
BusI7hI9h
MagnitudePhaseMagnitudePhase
10.0042−73.30470.0014−159.6064
20.0046−97.27430.0038−40.3494
30.0069−22.28870.0040−2.4660
40.0164−128.36700.0235177.1767
50.031889.59320.017826.5938
60.0119170.70230.0157171.3318
70.0115−73.61660.0187165.3096
80.002897.43520.0067−17.4348
90.0116108.45480.0241−18.7536
100.0060−79.78520.0135163.6955
110.0025−26.86290.0033−13.1680
120.0087−9.39970.0103−7.8251
130.0039−35.34700.0050−21.7709
140.0023−70.94860.0034161.6506
Table A8. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order voltage harmonics for the IEEE 14-bus system.
Table A8. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order voltage harmonics for the IEEE 14-bus system.
BusV7hV9h
MagnitudePhaseMagnitudePhase
10.00104.13410.004254.3987
20.0010−5.09450.004455.7077
30.004461.23620.004566.6699
40.0020−80.69710.0039−78.4343
50.0060−158.36380.0024123.9482
60.001118.03560.0023−57.4554
70.005220.71740.0072−75.0046
80.001724.77580.004750.1152
90.001617.51930.004247.1212
100.004824.00060.0053−72.8899
110.004440.09470.0017−20.4970
120.010576.26870.013181.7285
130.006363.63260.006370.9736
140.005841.78600.0017−10.1490
Table A9. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order current harmonics for the IEEE 14-bus system.
Table A9. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order current harmonics for the IEEE 14-bus system.
BusI11hI13h
MagnitudePhaseMagnitudePhase
10.000045.00550.0040−26.2923
20.0003−109.02050.0048147.8820
30.0004−5.26980.00112.9953
40.0019128.73420.004153.8713
50.001016.33470.001732.1031
60.0014158.95900.0013158.3837
70.0020134.71390.002768.3087
80.0007−51.71500.0013−125.0988
90.0025−57.03380.0045−115.7826
100.0015132.57320.002174.3176
110.0003−30.30410.0002−61.6939
120.0007−18.94400.0007−18.7090
130.0005−42.87470.0004−79.9085
140.000421.19370.001021.3771
Table A10. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order voltage harmonics for the IEEE 14-bus system.
Table A10. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order voltage harmonics for the IEEE 14-bus system.
BusV11hV13h
MagnitudePhaseMagnitudePhase
10.000814.70250.00420.0000
20.000811.73250.0038−51.7604
30.000651.42350.00077.6872
40.0004−123.84580.0009166.9942
50.0001124.22320.0001131.4582
60.0002−64.73780.0004−45.8265
70.0007−99.82680.0005−140.0438
80.000924.75200.0029−33.3961
90.000828.79820.0020−16.5888
100.0005−97.99580.0003−132.1348
110.0001−15.98270.0003−18.9810
120.001369.72220.001561.6654
130.000962.00500.001044.9332
140.001279.70560.001878.7671

Appendix B.2

Table A11, Table A12, Table A13, Table A14, Table A15 and Table A16 present the fundamental and harmonic component values of current and voltage derived from the power flow and harmonic state analysis for the IEEE 30-bus system using the electromagnetic transient software PSCAD. These values are based on the current injection harmonics in randomly selected buses.
Table A11. Reference values of magnitude in p.u. and phase in degrees for the fundamental current and the 3rd- and 5th-order current harmonics for the IEEE 30-bus system.
Table A11. Reference values of magnitude in p.u. and phase in degrees for the fundamental current and the 3rd- and 5th-order current harmonics for the IEEE 30-bus system.
BusI1hI3hI5h
MagnitudePhaseMagnitudePhaseMagnitudePhase
12.42254.63030.0297101.79290.0277106.5383
20.3117−71.07120.0334101.64140.0339105.7678
30.0429114.53590.0212−56.69310.0201−19.8704
40.0914134.81750.0210−56.97810.0191−19.7931
50.6599176.91640.007794.65040.007390.4240
61.3085−4.01810.0446−86.31260.0646−63.5168
71.5177166.88490.055997.55240.055599.1196
80.4855178.14530.040893.27190.048994.4718
90.010078.08380.0007170.17960.000063.3482
100.2440−91.00780.0199−87.65240.0680−104.0201
110.1528−109.53620.007796.07480.009287.9318
120.1329128.18850.0257−66.12000.0016−61.6546
130.1036−111.28790.0129105.87290.008799.0378
140.0612150.47830.0048−88.18260.0030−44.4672
150.0819147.95320.0100−74.62880.0022−52.2500
160.0658144.69280.0161−64.14850.014623.7366
170.1221131.51750.0176135.18110.0238−152.6003
180.0320148.84930.0013177.87900.0017164.9190
190.0973144.71580.0046−100.57310.0010−80.6940
200.0265−40.67670.0037−33.60550.007650.4356
210.1989131.97940.0049−132.86600.0055151.1641
220.1321134.53460.0041−109.90210.0043130.7205
230.0344137.96970.0106−63.97120.0040−29.9911
240.0864150.33530.0053−133.67970.0140−134.4625
250.0002−105.99060.0000−98.64020.0000−119.1604
260.0415130.98710.0011161.91040.0017145.4195
270.010177.15640.0004173.34160.000077.4757
280.038778.59740.001674.37190.003786.3870
290.0254143.63330.0006163.87300.0009146.9362
300.1080152.98050.0025164.32740.0040143.5933
Table A12. Reference values of magnitude in p.u. and phase in degrees for the fundamental voltage and the 3rd- and 5th-order voltage harmonics for the IEEE 30-bus system.
Table A12. Reference values of magnitude in p.u. and phase in degrees for the fundamental voltage and the 3rd- and 5th-order voltage harmonics for the IEEE 30-bus system.
BusV1hV3hV5h
MagnitudePhaseMagnitudePhaseMagnitudePhase
11.06010.00000.00520.00000.00730.0000
21.0427−5.35120.0057−1.23330.0084−2.6467
31.0212−7.05930.01786.73740.024615.4302
41.0121−8.69740.01814.44000.02529.7507
51.0084−15.21660.0013−14.73330.0018−28.2391
61.0110−10.22070.0150−0.44750.02412.2046
70.9878−16.52000.0012−18.71180.0019−36.3926
81.0098−11.36160.0109−0.16090.01631.3840
91.0503−13.38010.0268−3.75030.0496−16.3213
101.0447−14.93230.0354−3.44050.0684−18.7517
111.0819−13.56060.0221−5.86400.0404−19.6643
121.0550−14.42370.04124.83730.0441−6.1853
131.0694−14.51650.03593.18530.0383−8.5845
141.0406−15.26930.04483.27480.0489−6.3038
151.0366−15.32010.04472.82810.0503−9.1232
161.0406−15.08960.04242.49620.0519−8.8895
171.0373−15.23100.0352−4.72720.0679−20.7700
181.0304−15.72520.04310.39660.0544−12.7432
191.0297−15.77450.0423−0.43530.0570−13.5611
201.0347−15.51630.0410−0.44600.0580−14.0665
211.0275−15.56460.0369−4.65230.0668−20.8442
221.0264−15.61080.0371−4.62550.0665−20.9705
231.0246−15.72660.04571.31090.0566−13.1335
241.0169−15.92750.0387−4.47810.0638−21.8784
251.0132−15.42300.0299−6.62770.0482−23.5341
260.9957−15.84260.0294−8.83840.0473−27.2579
271.0193−14.85020.0247−7.49290.0389−22.8831
281.0064−10.92480.0152−1.87670.0242−2.4239
290.9996−16.08120.0241−11.57980.0376−29.8528
300.9882−16.96390.0239−14.36490.0370−34.6621
Table A13. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order current harmonics for the IEEE 30-bus system.
Table A13. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order current harmonics for the IEEE 30-bus system.
BusI7hI9h
MagnitudePhaseMagnitudePhase
10.0141110.43830.0413113.4899
20.0163104.22810.0483112.3841
30.0196−26.80880.0306−40.0922
40.0196−27.22850.0311−43.7682
50.003381.17790.009181.5128
60.0330−53.70700.0733−51.3579
70.019487.84180.0444114.8477
80.015159.01120.0504122.2711
90.0000−113.36070.0000164.5245
100.0822−177.09480.0943135.7174
110.005224.25340.0029−10.7531
120.0040−28.45710.0097−64.0422
130.002576.95120.0045125.7667
140.0046−31.64950.0104−57.5931
150.0043−25.65930.0021−176.6715
160.0273−19.55030.0385−56.4748
170.0347154.79660.0462116.7126
180.0008113.29800.0004148.0787
190.0041−6.41420.0012102.6376
200.0107−12.72640.0175−53.5896
210.008346.67720.008338.7982
220.006926.15700.0080−29.1297
230.0051−25.68020.0110−47.8208
240.0144162.71080.0113132.1101
250.0000172.20810.0000155.8399
260.001179.53540.000566.4725
270.0000144.87360.0000147.7954
280.002365.03750.0076102.7021
290.000582.92380.0002118.4700
300.002276.40900.0010109.8880
Table A14. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order voltage harmonics for the IEEE 30-bus system.
Table A14. Reference values of magnitude in p.u. and phase in degrees for the 7th- and 9th-order voltage harmonics for the IEEE 30-bus system.
BusV7hV9h
MagnitudePhaseMagnitudePhase
10.00490.00000.01720.0000
20.0051−7.89350.0171−1.6754
30.017923.85680.058625.0809
40.015314.15170.050122.6895
50.0011−44.61490.0034−48.9590
60.0117−9.62840.032315.8900
70.0012−57.00780.0042−67.0073
80.0091−2.04140.02000.7528
90.0382−84.42560.0264−123.5510
100.0607−89.00560.0573−133.8570
110.0310−88.93320.0215−129.2516
120.0172−33.03420.037611.3306
130.0149−36.27150.03257.2377
140.0187−29.90110.042712.1779
150.0196−46.91060.02883.9094
160.0220−54.05630.02574.7365
170.0603−91.80580.0568−137.5212
180.0250−65.46570.0150−29.4747
190.0290−71.39680.0132−71.8265
200.0330−75.43170.0151−89.6589
210.0547−91.81770.0490−137.2982
220.0539−91.78080.0466−136.4792
230.0264−69.45140.0182−5.8386
240.0471−91.24770.0292−125.2941
250.0310−90.75130.0143−100.3610
260.0304−96.02230.0139−106.9738
270.0212−86.06160.0109−48.0127
280.0111−20.07520.029311.1500
290.0203−95.68080.0103−59.6021
300.0199−102.38930.0100−67.8574
Table A15. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order current harmonics for the IEEE 30-bus system.
Table A15. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order current harmonics for the IEEE 30-bus system.
BusI11hI13h
MagnitudePhaseMagnitudePhase
10.0883115.12300.0595115.7186
20.1114113.59600.0834116.7914
30.0331−6.53370.0285137.1271
40.0321−13.76710.046430.2879
50.018974.58160.012069.5983
60.1135−38.54790.0226−39.0112
70.0938121.94640.0649129.4856
80.1124130.85320.0589128.4795
90.0001159.36230.0000163.4049
100.1464123.58150.1146112.9497
110.0012−7.11670.0003−69.6160
120.0078−163.26420.0118−15.3862
130.0080137.94490.0067163.2837
140.0147−48.32040.0137−8.9153
150.0143−48.53220.0054−145.3856
160.0508−56.62060.0239−61.6100
170.0581114.22340.0332118.2853
180.0013−167.75960.0014−143.2664
190.0023−178.15610.0037−21.8539
200.0220−54.10450.0102−54.8897
210.0205−16.08480.017712.2930
220.0178−57.98620.0227−86.0469
230.0013−178.50320.01483.1012
240.0131−162.76410.0601−124.0382
250.0000−116.86350.0002−117.5470
260.0011141.67390.0027144.2810
270.0000−155.10660.0001−124.5232
280.0249110.80530.0257117.4348
290.0009155.76380.0014148.0514
300.0041144.40100.0061134.3415
Table A16. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order voltage harmonics for the IEEE 30-bus system.
Table A16. Reference values of magnitude in p.u. and phase in degrees for the 11th- and 13th-order voltage harmonics for the IEEE 30-bus system.
BusV11hV13h
MagnitudePhaseMagnitudePhase
10.04140.00000.03000.0000
20.0423−0.98040.03292.9704
30.128132.33110.077839.7848
40.111729.96430.080144.3757
50.0080−58.12720.0057−63.2138
60.079124.25930.062134.4277
70.0113−83.03250.0095−98.4306
80.05045.40610.053629.7745
90.0130−123.61120.0034170.5224
100.0615−144.66620.0378−151.1876
110.0106−130.51340.0027162.3888
120.078419.46350.074141.2042
130.067814.50410.064035.3749
140.092822.70130.088745.4470
150.072819.06590.071433.6934
160.061924.88060.047839.5820
170.0609−149.09680.0374−156.5337
180.043111.25050.047034.8823
190.02622.40010.032840.5601
200.0178−3.24160.022443.2083
210.0359−144.80970.0131−91.0183
220.0323−139.26600.0201−56.5909
230.04471.85680.079018.5544
240.0275−62.79140.1070−25.5891
250.0334−27.86840.0800−24.5293
260.0326−35.96360.0779−34.1237
270.0427−10.77820.0648−17.0999
280.078318.74120.072228.0978
290.0395−25.07530.0587−33.5944
300.0382−35.36650.0563−45.7243

References

  1. Ruuth, K.; Hilden, A.; Rekola, J.; Pakonen, P.; Verho, P. The Impact of LED Lighting Systems to the Power Quality and Recommendations for Installation Methods to Achieve the Expected Energy Efficiency. In Proceedings of the 25th International Conference on Electricity Distribution, Madrid, Spain, 3–6 June 2019. [Google Scholar]
  2. Schweppe, F.C.; Wildes, J. Power System Static-State Estimation Part I: Exact Model. IEEE Trans. Power Appar. Syst. 1970, PAS-89, 120–125. [Google Scholar] [CrossRef]
  3. Schweppe, F.C.; Rom, D.B. Power System Static-State Estimation Part II: Approximate Model. IEEE Trans. Power Appar. Syst. 1970, PAS-89, 125–130. [Google Scholar] [CrossRef]
  4. Schweppe, F.C. Power System Static-State Estimation Part III: Implementation. IEEE Trans. Power Appar. Syst. 1970, PAS-89, 130–135. [Google Scholar] [CrossRef]
  5. Heydt, G.T. Identification of Harmonic Sources by a State Estimation Technique. IEEE Trans. Power Deliv. 1989, 4, 569–576. [Google Scholar] [CrossRef]
  6. Beides, H.M.; Heydt, G.T. Dynamic state estimation of power system harmonics using Kalman filter methodology. IEEE Trans. Power Deliv. 1991, 6, 1663–1670. [Google Scholar] [CrossRef]
  7. Hartana, R.; Richards, G. Harmonic source monitoring and identification using neural networks. IEEE Trans. Power Syst. 1990, 5, 1098–1104. [Google Scholar] [CrossRef]
  8. Najjar, M.Y.; Heydt, G.T. A Hybrid Nonlinear-Least Squares Estimation of Harmonic Signal Levels in Power Systems. IEEE Trans. Power Deliv. 1991, 6, 282–288. [Google Scholar] [CrossRef]
  9. Mori, H.; Itou, K.; Uematsu, H.; Tsuzuki, S. An Artificial Neural-Net Based Method for Predicting Power System Voltage Harmonics. IEEE Trans. Power Deliv. 1992, 7, 402–409. [Google Scholar] [CrossRef]
  10. Lobos, T.; Kozina, T.; Koglin, H.-J. Power System Harmonics Estimation Using Linear Least Squares Method and SVD. IEE Proc. Gener. Transm. Distrib. 2001, 148, 567–572. [Google Scholar] [CrossRef]
  11. Yu, K.K.C.; Watson, N.R. Three-Phase Harmonic State Estimation Using SVD for Partially Observable Systems. In Proceedings of the 2004 International Conference on Power System Technology (PowerCon 2004); IEEE: Piscataway, NJ, USA, 2004; Volume 1, pp. 29–34. [Google Scholar]
  12. Negnevitsky, M.; Ringrose, M. Monitoring Multiple Harmonic Sources in Power Systems Using Neural Networks. In Proceedings of the 2005 IEEE Russia Power Tech; IEEE: Piscataway, NJ, USA, 2005; pp. 1–6. [Google Scholar]
  13. Arruda, E.F.; Kagan, N.; Ribeiro, P.F. Harmonic Distortion State Estimation Using an Evolutionary Strategy. IEEE Trans. Power Deliv. 2010, 25, 831–842. [Google Scholar] [CrossRef]
  14. Arruda, E.F.; Kagan, N.; Ribeiro, P.F. Three-Phase Harmonic Distortion State Estimation Algorithm Based on Evolutionary Strategies. Electr. Power Syst. Res. 2010, 80, 1024–1032. [Google Scholar] [CrossRef]
  15. Sepulchro, W.N.; Encarnação, L.F.; Brunoro, M. Harmonic Distortion and Power Flow State Estimation for Distribution Systems Based on Evolutionary Strategies. IEEE Lat. Am. Trans. 2015, 13, 3066–3071. [Google Scholar] [CrossRef]
  16. Santos, G.G.; Oliveira, T.L.; Oliveira, J.C.; Vieira, J.C.M. A Hybrid Method for Harmonic State Estimation in Partially Observable Systems. Int. Trans. Electr. Energy Syst. 2021, 31, e12763. [Google Scholar] [CrossRef]
  17. Warid, W.; Hizam, H.; Mariun, N.; Abdul Wahab, N.I. A novel quasi-oppositional modified Jaya algorithm for multi-objective optimal power flow solution. Appl. Soft Comput. 2018, 65, 360–373. [Google Scholar] [CrossRef]
  18. Rizwan, P.; Vikas, C.; Manasa, Y. Optimal Power Flow Analysis and Transmission System with Jaya Algorithm for Reducing Power Loss. J. Syst. Eng. Electron. 2024, 34. Available online: https://jseepublisher.com/wp-content/uploads/10-JSEE2104.pdf (accessed on 30 April 2024).
  19. Sepulchro, W.N.; Encarnação, L.F. Harmonic State and Power Flow Estimation in Distribution Systems Using Jaya Algorithm. In 2024 IEEE 18th International Conference on Compatibility, Power Electronics and Power Engineering (CPE-POWERENG); IEEE: Gdynia, Poland, 2024; pp. 1–6. [Google Scholar]
  20. Rechenberg, I. Cybernetic Solution Path of an Experimental Problem; Technical Report; Technical University of Berlin: Berlin, Germany, 1965. [Google Scholar]
  21. Schwefel, H.-P. Numerical Optimization of Computer Models; Wiley: Hoboken, NJ, USA, 1974. [Google Scholar]
  22. Rechenberg, I. Evolutionsstrategie: Optimierung Technischer Systeme Nach Prinzipien der Biologischen Evolution; Frommann-Holzboog: Stuttgart, Germany, 1971. [Google Scholar]
  23. Schwefel, H.-P. Numerische Optimierung von Computer-Modellen Mittels der Evolutionsstrategie. Ph.D. Thesis, Technische Universität Berlin, Berlin, Germany, 1981. [Google Scholar]
  24. Back, T.; Schwefel, H.-P. Evolutionary Computation: An Overview. In Proceedings of the IEEE International Conference on Evolutionary Computation, Nagoya, Japan, 20–22 May 1996. [Google Scholar]
  25. Eiben, A.E.; Smith, J.E. Introduction to Evolutionary Computing; Springer: New York, NY, USA, 2015. [Google Scholar]
  26. Rao, R. Jaya: A Simple and New Optimization Algorithm for Solving Constrained and Unconstrained Optimization Problems. Int. J. Ind. Eng. Comput. 2016, 7, 19–34. [Google Scholar]
  27. Rao, R.V. Jaya: An Advanced Optimization Algorithm and Its Engineering Applications; Springer: New York, NY, USA, 2019. [Google Scholar]
  28. Warid, W.; Hizam, H.; Mariun, N.; Abdulwahab, N.I. Optimization of Power Flow Using the Jaya Algorithm. Electr. Power Syst. Res. J. 2016, 140, 297–307. [Google Scholar]
  29. Gupta, S.; Kumar, N.; Srivastava, L.; Malik, H.; Pliego Marugán, A.; García Márquez, F.P. A Hybrid Jaya–Powell’s Pattern Search Algorithm for Multi-Objective Optimal Power Flow Incorporating Distributed Generation. Energies 2021, 14, 2831. [Google Scholar] [CrossRef]
  30. Grady, W.M.; Santoso, S. Understanding power system harmonics. IEEE Power Eng. Rev. 2001, 21, 8–11. [Google Scholar] [CrossRef]
  31. IEEE 14-Bus System. Available online: https://icseg.iti.illinois.edu/ieee-14-bus-system/ (accessed on 18 April 2024).
  32. IEEE 30-Bus System. Available online: https://icseg.iti.illinois.edu/ieee-30-bus-system/ (accessed on 18 April 2024).
  33. Manitoba Hydro International Ltd. IEEE 14 Bus System Technical Note. In Technical Documents of Manitoba Hydro International Ltd.; Revision 1; Manitoba Hydro International Ltd.: Winnipeg, MB, Canada, 2018; Available online: https://www.pscad.com/knowledge-base/download/ieee_14_bus_technical_note.pdf (accessed on 26 April 2024).
  34. Manitoba Hydro International Ltd. IEEE 30 Bus System Technical Note. In Technical Documents of Manitoba Hydro International Ltd.; Revision 1; Manitoba Hydro International Ltd.: Winnipeg, MB, Canada, 2018; Available online: https://www.pscad.com/knowledge-base/download/ieee_30_bus_technical_note.pdf (accessed on 26 April 2024).
  35. Gonzalez-Longatt, F.M. Power System Test Cases. Available online: http://fglongatt.org/OLD/Test_Case_IEEE_30.html (accessed on 25 February 2024).
Figure 1. Basic flowchart of the ES algorithm.
Figure 1. Basic flowchart of the ES algorithm.
Electronics 13 03559 g001
Figure 2. Flowchart of the basic Jaya method.
Figure 2. Flowchart of the basic Jaya method.
Electronics 13 03559 g002
Figure 3. Block diagram for the proposed methodology.
Figure 3. Block diagram for the proposed methodology.
Electronics 13 03559 g003
Figure 6. Flowchart of the meter allocation method presented.
Figure 6. Flowchart of the meter allocation method presented.
Electronics 13 03559 g006
Figure 7. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 14-bus system using the ES algorithm.
Figure 7. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 14-bus system using the ES algorithm.
Electronics 13 03559 g007
Figure 8. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 14-bus system using Jaya algorithm.
Figure 8. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 14-bus system using Jaya algorithm.
Electronics 13 03559 g008
Figure 9. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 30-bus system using the ES algorithm.
Figure 9. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 30-bus system using the ES algorithm.
Electronics 13 03559 g009
Figure 10. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 30-bus system using the Jaya algorithm.
Figure 10. Absolute voltage errors for each harmonic order and unmonitored bus for the IEEE 30-bus system using the Jaya algorithm.
Electronics 13 03559 g010
Figure 11. Error evolution curves for each harmonic order for the IEEE 14-bus system using the ES algorithm.
Figure 11. Error evolution curves for each harmonic order for the IEEE 14-bus system using the ES algorithm.
Electronics 13 03559 g011
Figure 12. Error evolution curves for each harmonic order for the IEEE 14-bus system using the Jaya algorithm.
Figure 12. Error evolution curves for each harmonic order for the IEEE 14-bus system using the Jaya algorithm.
Electronics 13 03559 g012
Figure 13. Error evolution curves for each harmonic order for the IEEE 30-bus system using ES algorithm.
Figure 13. Error evolution curves for each harmonic order for the IEEE 30-bus system using ES algorithm.
Electronics 13 03559 g013
Figure 14. Error evolution curves for each harmonic order for the IEEE 30-bus system using the Jaya algorithm.
Figure 14. Error evolution curves for each harmonic order for the IEEE 30-bus system using the Jaya algorithm.
Electronics 13 03559 g014
Figure 15. Comparison of THD errors between ES and Jaya methodologies for IEEE 14-bus system.
Figure 15. Comparison of THD errors between ES and Jaya methodologies for IEEE 14-bus system.
Electronics 13 03559 g015
Figure 16. Comparison of THD errors between ES and Jaya methodologies for IEEE 30-bus system.
Figure 16. Comparison of THD errors between ES and Jaya methodologies for IEEE 30-bus system.
Electronics 13 03559 g016
Figure 17. Comparison between the average times of the Jaya algorithm and the ES algorithm for the IEEE 14-bus system.
Figure 17. Comparison between the average times of the Jaya algorithm and the ES algorithm for the IEEE 14-bus system.
Electronics 13 03559 g017
Figure 18. Comparison between the average times of the Jaya algorithm and the ES algorithm for the IEEE 30-bus system.
Figure 18. Comparison between the average times of the Jaya algorithm and the ES algorithm for the IEEE 30-bus system.
Electronics 13 03559 g018
Figure 19. Comparison of the average total times between the Jaya algorithm and the ES algorithm for the IEEE 14-bus system and the IEEE 30-bus system.
Figure 19. Comparison of the average total times between the Jaya algorithm and the ES algorithm for the IEEE 14-bus system and the IEEE 30-bus system.
Electronics 13 03559 g019
Table 1. List of buses, number of connections, and neighboring buses for the IEEE 14-bus system.
Table 1. List of buses, number of connections, and neighboring buses for the IEEE 14-bus system.
BusConnectionsNeighbors
452, 3, 5, 7, 9
241, 3, 4, 5
541, 2, 4, 6
645, 11, 12, 13
944, 7, 10, 14
734, 8, 9
1336, 12, 14
122, 5
322, 4
1029, 11
1126, 10
1226, 13
1429, 13
817
Table 2. List of buses, number of connections, and neighboring buses for the IEEE 30-bus system.
Table 2. List of buses, number of connections, and neighboring buses for the IEEE 30-bus system.
BusConnectionsNeighbors
672, 4, 7, 8, 9, 10, 28
1066, 9, 17, 20, 21, 22
1254, 13, 14, 15, 16
241, 4, 5, 6
442, 3, 6, 12
15412, 14, 18, 23
27425, 28, 29, 30
936, 10, 11
22310, 21, 24
24322, 23, 25
25324, 26, 27
2836, 8, 27
122, 3
321, 4
522, 7
725, 6
826, 28
14212, 15
16212, 17
17210, 16
18215, 19
19218, 20
20210, 19
21210, 22
23215, 24
29227, 30
30227, 29
1119
13112
26125
Table 3. Parameters for Evolutionary Strategy methodology.
Table 3. Parameters for Evolutionary Strategy methodology.
ParameterValue for IEEE 14-BusValue for IEEE 30-Bus
Maximum number of generations20002000
Number of individuals in the initial population40150
Mutation operations per individual55
Recombination rate10%20%
Initial mutation step for voltage magnitudes11
Initial mutation step for voltage phase 2 π 3 2 π 3
Self-adaptation parameter β 55
Table 4. Parameters for Jaya.
Table 4. Parameters for Jaya.
ParameterValue for IEEE 14-BusValue for IEEE 30-Bus
Maximum number of generations20005000
Number of individuals in the initial population50150
Table 5. Absolute errors of estimated voltage (p.u.) for the IEEE 14-bus system using ES algorithm.
Table 5. Absolute errors of estimated voltage (p.u.) for the IEEE 14-bus system using ES algorithm.
BusHarmonic Order
135791113
20.0020.0010.0000.0000.0000.0000.000
30.0270.0100.0010.0010.0000.0000.000
50.0510.0030.0020.0040.0000.0000.000
70.0130.0000.0000.0020.0000.0000.000
90.0010.0000.0000.0000.0000.0000.000
110.0120.0000.0000.0000.0000.0000.000
120.0120.0000.0010.0050.0060.0000.000
130.0060.0000.0000.0000.0000.0000.000
Maximum Error0.0510.0100.0020.0050.0060.0000.000
Table 6. Absolute errors of estimated voltage (p.u.) for the IEEE 14-bus system using Jaya algorithm.
Table 6. Absolute errors of estimated voltage (p.u.) for the IEEE 14-bus system using Jaya algorithm.
BusHarmonic Order
135791113
20.0000.0000.0000.0000.0000.0000.000
30.0020.0020.0000.0000.0000.0000.000
50.0010.0010.0000.0000.0000.0000.000
70.0010.0010.0000.0000.0000.0000.000
90.0000.0000.0000.0000.0000.0000.000
110.0030.0010.0000.0000.0010.0000.000
120.0000.0010.0000.0000.0010.0000.000
130.0000.0020.0000.0000.0000.0000.000
Maximum Error0.0030.0020.0000.0000.0010.0000.000
Table 7. Absolute errors of estimated voltage (p.u.) for the IEEE 30-bus system using ES algorithm.
Table 7. Absolute errors of estimated voltage (p.u.) for the IEEE 30-bus system using ES algorithm.
BusHarmonic Order
135791113
10.0320.0010.0000.0030.0010.0020.002
20.0150.0010.0000.0000.0000.0020.002
40.0000.0000.0000.0000.0030.0010.001
70.0080.0000.0010.0010.0010.0000.001
80.0000.0000.0010.0000.0010.0010.000
90.0220.0000.0010.0000.0000.0010.001
100.0160.0000.0010.0010.0000.0020.000
130.0430.0000.0010.0010.0010.0020.003
140.0130.0020.0010.0010.0000.0020.002
150.0090.0020.0010.0000.0030.0000.004
160.0130.0000.0000.0010.0000.0010.005
190.0030.0000.0010.0010.0010.0010.000
220.0000.0000.0000.0000.0000.0000.000
230.0010.0030.0000.0040.0010.0010.000
250.0000.0000.0000.0010.0000.0020.002
280.0000.0020.0010.0000.0000.0060.001
290.0000.0020.0020.0000.0010.0030.001
300.0080.0010.0010.0000.0000.0010.002
Maximum Error0.0430.0030.0020.0040.0030.0060.005
Table 8. Absolute errors of estimated voltage (p.u.) for the IEEE 30-bus system using Jaya algorithm.
Table 8. Absolute errors of estimated voltage (p.u.) for the IEEE 30-bus system using Jaya algorithm.
BusHarmonic Order
135791113
10.0110.0010.0010.0000.0020.0070.003
20.0030.0010.0040.0000.0010.0070.002
40.0050.0000.0020.0010.0000.0020.001
70.0050.0000.0020.0010.0000.0040.001
80.0010.0010.0000.0000.0000.0020.003
90.0010.0010.0000.0020.0020.0010.000
100.0030.0000.0000.0000.0000.0000.003
130.0060.0000.0010.0010.0020.0040.002
140.0010.0000.0030.0000.0040.0020.005
150.0030.0000.0000.0010.0020.0010.003
160.0000.0000.0010.0010.0010.0010.000
170.0000.0000.0020.0010.0010.0010.001
180.0020.0010.0020.0010.0010.0010.001
190.0050.0000.0020.0010.0010.0010.001
200.0030.0030.0020.0010.0010.0020.006
210.0030.0010.0030.0000.0010.0050.002
Maximum Error0.0110.0030.0040.0020.0040.0070.006
Table 9. Reference and estimated values of THD and their absolute errors for the IEEE 14-bus system using the ES algorithm.
Table 9. Reference and estimated values of THD and their absolute errors for the IEEE 14-bus system using the ES algorithm.
BusReference THD (%)Estimated THD (%)Error (%)
22.4212.3780.043
31.8442.8330.990
54.5924.4800.112
73.5563.5590.003
92.7532.7710.017
112.4742.5110.037
122.6572.2620.395
132.2592.2670.008
Maximum Error 0.990
Medium Error 0.115
Table 10. Reference and estimated values of THD and their absolute errors for the IEEE 14-bus system using the Jaya algorithm.
Table 10. Reference and estimated values of THD and their absolute errors for the IEEE 14-bus system using the Jaya algorithm.
BusReference THD (%)Estimated THD (%)Error (%)
22.4212.4430.022
31.8442.0340.190
54.5924.4750.117
73.5563.4640.092
92.7532.7950.042
112.4742.5640.090
122.6572.6770.020
132.2592.4340.175
Maximum Error 0.190
Medium Error 0.053
Table 11. Reference and estimated values of THD and their absolute errors for the IEEE 30-bus system using the ES algorithm.
Table 11. Reference and estimated values of THD and their absolute errors for the IEEE 30-bus system using the ES algorithm.
BusReference THD (%)Estimated THD (%)Error (%)
15.1765.1560.020
25.5035.8600.357
414.85314.8620.009
71.5701.7240.153
1012.88013.0670.187
1310.55310.9410.388
1414.59214.6570.064
1512.26012.1550.105
1610.41810.7780.360
198.5878.5320.055
2210.81310.8220.010
2311.77911.5830.197
2510.77310.8650.092
2811.39210.9460.446
298.6758.7890.115
308.5098.5170.008
Maximum Error 0.446
Medium Error 0.095
Table 12. Reference and estimated values of THD and their absolute errors for the IEEE 30-bus system using the Jaya algorithm.
Table 12. Reference and estimated values of THD and their absolute errors for the IEEE 30-bus system using the Jaya algorithm.
BusReference THD (%)Estimated THD (%)Error (%)
15.1765.5320.355
25.5034.8940.609
414.85315.0570.204
71.5701.8860.316
1012.88012.9080.028
1310.55310.3550.198
1414.59214.7670.174
1512.26012.4460.186
1610.41810.6890.271
198.5878.6580.071
2210.81310.7830.030
2311.77911.7180.061
2510.77310.8820.109
2811.39211.4660.074
298.6758.4770.198
308.5098.3330.176
Maximum Error 0.609
Medium Error 0.111
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sepulchro, W.d.N.; Encarnação, L.F. Harmonic State Estimation in Power Systems Using the Jaya Algorithm. Electronics 2024, 13, 3559. https://doi.org/10.3390/electronics13173559

AMA Style

Sepulchro WdN, Encarnação LF. Harmonic State Estimation in Power Systems Using the Jaya Algorithm. Electronics. 2024; 13(17):3559. https://doi.org/10.3390/electronics13173559

Chicago/Turabian Style

Sepulchro, Walace do Nascimento, and Lucas Frizera Encarnação. 2024. "Harmonic State Estimation in Power Systems Using the Jaya Algorithm" Electronics 13, no. 17: 3559. https://doi.org/10.3390/electronics13173559

APA Style

Sepulchro, W. d. N., & Encarnação, L. F. (2024). Harmonic State Estimation in Power Systems Using the Jaya Algorithm. Electronics, 13(17), 3559. https://doi.org/10.3390/electronics13173559

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

Article Metrics

Back to TopTop