Abstract
The accurate estimation of power signal parameters allows smart grids to optimize power delivery efficiency, improve equipment utilization, and control power flow among generation nodes and loads. However, practically it becomes a challenging task because of the presence of harmonic distortions. In this study, a parameter estimation of the power system harmonics is investigated through swarm intelligence–based optimization strength of the cuckoo search algorithm. The performance evaluation is conducted in detail for different generations and particle sizes and for different signal-to-noise ratios. The simulation results reveal that the cuckoo search optimization heuristic accurately estimates the amplitude and phase parameters of the power system harmonics and is robust against different signal-to-noise ratios.
1 Introduction
The control of the power systems and synchronization of grid-connected electrical devices require parameters of voltage and current signals, such as amplitude, frequency, and phase (; ). Accurate information on these signal parameters allows smart grids to optimize power delivery efficiency, improve equipment utilization, and control power flow among generation nodes and loads (Sun et al., 2019). However, practically the presence of harmonics and inter-harmonic distortions in power signals complicates the situation because the amplitude and phase estimations of harmonics is challenging for evaluating the quality of service characteristics in smart grids (). Thus, it is imperative to investigate the development of accurate, robust, and stable estimation methods to mitigate the adverse effects of harmonics in smart grid efficiency.
Different researchers have investigated the domain of estimating power system harmonics and proposed various techniques. For example, ) developed the ADALINE-based least mean square (LMS) algorithm and compared its performance with the recursive least square, while ) presented the novel smooth variable structure filters for voltage harmonics and voltage imbalance. ) introduced the filter bank–based ESPRIT approach for increased efficiency in harmonics detection and estimation. developed a hybrid estimation scheme by combining recursive least squares with the iterated extended Kalman filter, and ) presented the frequency shifting/filtering method. Xu et al. (Xu and Ding, 2018; Xu and Song, 2020) investigated the parameter estimation of power signals using iterative/recursive methods, such as multi-innovation stochastic gradient (Xu and Ding, 2017), separable multi-synchronous multi-innovation gradient (Xu et al., 2022), separable Newton recursive (Xu, 2022a), separable multi-innovation Newton iterative (Xu, 2022b), and hierarchical principle-based recursive least squares (Xu et al., 2021). Chaudhary et al. introduced fractional gradient-based estimation algorithms for power signals, such as fractional LMS () and innovative fractional order LMS (), while Zubair et al. (2018) presented the momentum term–based fractional scheme. , had applied evolutionary and swarming heuristics for parameter estimation of power signals.
The schemes based on swarm and evolutionary heuristics have established their significance through effective application in solving various challenging optimization tasks (; ; ; ; ) such as power system harmonics estimation (; ; ; ; ; Singh et al., 2016). Yang and Deb (2009), Yang and Deb (2014) introduced a metaheuristic inspired by the search mechanism of cuckoo called the cuckoo search optimization (CSO) algorithm. The CSO has been used extensively for effective solutions in a number of engineering and applied problems with many variants (; ), such as photovoltaic model (), social media sentiment analysis (), path planning (Song et al., 2020), power control in salt reactors (), damage detection infrastructures (Tran-Ngoc et al., 2019), chemoinformatics (), economic load dispatch (Yu et al., 2020), and many others. The effective performance of CSO in illustrated applications motivated us to explore parameter estimation power system harmonics through well-established optimization strength of the CSO. In this study, parameters of power system harmonics are estimated through CSO for different generation and particle size. Detailed and in-depth performance analyses are conducted to check the accuracy, diversity, and robustness of the CSO for harmonics estimation.
The remaining article is set as follows: the estimation model for power system harmonics along with the optimization procedure of CSO is described in Section 2. The results of simulation studies are elaborated in Section 3. The conclusions and future works are listed in Section 4.
2 Materials and methods
The system model for power harmonics signal is first introduced, then the proposed methodology for optimization of the fitness function for the estimation of the harmonic signal is presented in this section. The overall flow diagram of the proposed scheme is presented in Figure 1 by means of different process block structures.
FIGURE 1
2.1 System model
The general harmonic signal in terms of its constituent parameters () is defined asand the variables in Eq. 1 are defined as J represents harmonic order, is amplitude, is angular frequency represented as with as fundamental frequency, is phase, and d is AWGN. Then, by writing a discrete version of Eq. 1 after sampling with period p,, and by assuming , Eq. 1 is rewritten as
Expanding Eq. 2 through the fundamental trigonometric identitywith the assumption thatusing the assumptions given in Eq. 4 into Eq. 3 gives Eq. 5:
The harmonics estimation model with information-vector and parameter-vector is written as
The aim is to estimate the amplitude and phase parameters of the harmonics through minimizing the error-based criterion function defined as is the actual signal, is an estimated harmonic signal calculated through estimated parameter-vector by using the proposed CSO-based heuristic. The relationship between the intermediate variables and the actual parameters is given by
2.2 Optimization method: Cuckoo search optimization algorithm
Yang and Deb (2009), Yang and Deb (2014) introduced a metaheuristic inspired by the search mechanism of cuckoo called CSO with exhaustive applications in different fields of engineering design and optimization. The formulation of the CSO is based on three fundamental concepts: firstly, each cuckoo lays a single egg in a single instance of time, while dumping the respective egg in an arbitrary selected nest; secondly, the nests with the best fitness, i.e., having the best quality of eggs, proceed to the next generations; and thirdly, the number of host nests is fixed, and the probability of discovery of next by the host cuckoo is set between 0 and 1. The basis of these three fundamental steps in the CSO is proposed with the process flow structure as shown in Figure 1.
Proceeding by the actual behavior of the cuckoo, mathematical models were introduced by Yang et al.; by taking the new candidate solution xi(t+1) that represents an i-th cuckoo at t flight/iteration index and considering the Levy flight, we have the following expression iterative update of CSO ashere α represents the step size based on the scale of the optimization problem, and generally its value is set equal to unit value in most of the cases and the Levy flight is represented with a random walk procedure on the basis of Levy distribution ashere λ is a constant taken between 1 and 3, and the Levy flight has an infinite mean and variance. Further details on the mathematical terms, convergence proofs, and applications can be referred from citations mentioned therein: Yang and Deb (2009); Yang and Deb (2014); ; ; Tran-Ngoc et al. (2019); ; ; Song et al. (2020); Yu et al. (2020); ; . Inspired by the optimization performance of the CSO algorithm, we have implemented the CSO for optimization of parameters of the system models presented in Eqs. 1–8. The CSO algorithm was implemented in the presented research based on the routine available at the MATLAB Central File Exchange (Yang, 2022). The optimization strength of the CSO may be enhanced by integrating it with the kernel theory (; ; ; ; Sweis et al., 2022).
3 Results and discussion
Harmonics estimation of power systems is carried out in this research work by applying the CSO algorithm to two examples, and the results are given in a tabular form with the necessary discussion along with graphs. The simulations are conducted using MATLAB with a sampling frequency of 2 KHz. Three levels of additive white Gaussian noise d with 30 DB, 60 DB, and 90 DB are introduced in the system to check the heftiness of the proposed scheme. The investigation is carried out by considering four generation (G) sizes and three particle (P) sizes in the CSO. The considered values of G are 500, 1000, 1,500, and 2000, while the values of P are 50, 100, and 150.
Example 1The harmonic signal () considered in the first simulation study isTen parameters are taken under consideration in this problem. Five phase and five amplitude parameters that were taken are given below, while, , , , , and are the frequencies with values 50 Hz, 150 Hz, 250 Hz, 350 Hz, and 550 Hz, respectively:
Example 2The harmonic signal considered in the second simulation study isThe parameter-vector consisting of five amplitude and five phase parameters isThe results of Example 1 in terms of parameter estimates through the CSO along with the MSE value are given in Tables 1–4 for G = 500, 1,000, 1,500, and 2,000, respectively. While the respective results in the case of Example 2 are presented in Tables 5–8. The learning curves for Example 1 are given in Figure 2, while for Example 2, the convergence plots are provided in Figure 3. The results clearly indicate that the CSO gives a better accuracy for 90 dB SNR than it does for 60 dB and 30 dB. Moreover, the accuracy of the CSO for harmonics estimation increases by increasing the generation size.
TABLE 1
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.5028 | 0.4903 | 0.2059 | 0.1384 | 0.1010 | 1.3943 | 1.0096 | 0.7530 | 0.5487 | 0.5573 | 9.46E-04 |
| 100 | 1.4904 | 0.4991 | 0.1894 | 0.1581 | 0.0975 | 1.4002 | 1.0211 | 0.7449 | 0.6339 | 0.5726 | 9.70E-04 | |
| 150 | 1.5087 | 0.4904 | 0.2023 | 0.1463 | 0.0967 | 1.3955 | 1.0529 | 0.7615 | 0.6934 | 0.6441 | 8.64E-04 | |
| 60 | 50 | 1.5026 | 0.4985 | 0.1989 | 0.1496 | 0.0982 | 1.3965 | 1.0456 | 0.7820 | 0.6432 | 0.5329 | 1.02E-05 |
| 100 | 1.5003 | 0.5008 | 0.1978 | 0.1507 | 0.1007 | 1.3967 | 1.0410 | 0.7978 | 0.6495 | 0.5096 | 1.79E-05 | |
| 150 | 1.5005 | 0.5010 | 0.2018 | 0.1537 | 0.1005 | 1.3956 | 1.0455 | 0.7833 | 0.6221 | 0.5477 | 1.38E-05 | |
| 90 | 50 | 1.4988 | 0.5024 | 0.1988 | 0.1486 | 0.0999 | 1.3964 | 1.0468 | 0.7795 | 0.6266 | 0.5078 | 7.40E-06 |
| 100 | 1.4988 | 0.5009 | 0.1997 | 0.1483 | 0.1006 | 1.3954 | 1.0426 | 0.7934 | 0.6270 | 0.5264 | 7.11E-06 | |
| 150 | 1.4991 | 0.4975 | 0.2018 | 0.1504 | 0.1006 | 1.3965 | 1.0449 | 0.7793 | 0.6296 | 0.5074 | 8.01E-06 | |
| 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 0 |
Results of parameter estimates of Example 1 for G = 500.
TABLE 2
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.4952 | 0.4968 | 0.2042 | 0.1611 | 0.0934 | 1.4038 | 1.0503 | 0.7960 | 0.5386 | 0.5653 | 7.11E-04 |
| 100 | 1.4977 | 0.4896 | 0.1984 | 0.1546 | 0.1018 | 1.3969 | 1.0433 | 0.7826 | 0.6401 | 0.5742 | 7.17E-04 | |
| 150 | 1.5014 | 0.5028 | 0.2011 | 0.1480 | 0.0941 | 1.3938 | 1.0509 | 0.8318 | 0.6224 | 0.6175 | 7.12E-04 | |
| 60 | 50 | 1.5000 | 0.5000 | 0.1998 | 0.1501 | 0.1000 | 1.3961 | 1.0470 | 0.7862 | 0.6260 | 0.5211 | 7.07E-07 |
| 100 | 1.4999 | 0.5004 | 0.2002 | 0.1500 | 0.1001 | 1.3959 | 1.0468 | 0.7881 | 0.6290 | 0.5188 | 8.60E-07 | |
| 150 | 1.5001 | 0.4998 | 0.1999 | 0.1497 | 0.0997 | 1.3963 | 1.0468 | 0.7855 | 0.6267 | 0.5263 | 8.06E-07 | |
| 90 | 50 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7849 | 0.6280 | 0.5227 | 2.65E-09 |
| 100 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0471 | 0.7850 | 0.6280 | 0.5229 | 3.54E-09 | |
| 150 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1001 | 1.3960 | 1.0470 | 0.7850 | 0.6281 | 0.5228 | 5.67E-09 | |
| 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 0 |
Results of parameter estimates of Example 1 for G = 1,000.
TABLE 3
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.5062 | 0.5012 | 0.2063 | 0.1391 | 0.1015 | 1.3967 | 1.0482 | 0.7961 | 0.6084 | 0.4847 | 6.81E-04 |
| 100 | 1.4972 | 0.5052 | 0.1938 | 0.1573 | 0.1014 | 1.3949 | 1.0526 | 0.7824 | 0.6144 | 0.5761 | 6.08E-04 | |
| 150 | 1.4962 | 0.5003 | 0.1958 | 0.1505 | 0.0907 | 1.3955 | 1.0315 | 0.8229 | 0.6524 | 0.5335 | 6.34E-04 | |
| 60 | 50 | 1.5000 | 0.5001 | 0.1999 | 0.1501 | 0.0997 | 1.3961 | 1.0472 | 0.7848 | 0.6273 | 0.5270 | 6.83E-07 |
| 100 | 1.5000 | 0.5001 | 0.1997 | 0.1498 | 0.1001 | 1.3961 | 1.0467 | 0.7847 | 0.6292 | 0.5218 | 6.23E-07 | |
| 150 | 1.5000 | 0.5000 | 0.2000 | 0.1498 | 0.0998 | 1.3961 | 1.0467 | 0.7841 | 0.6268 | 0.5203 | 6.57E-07 | |
| 90 | 50 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6279 | 0.5230 | 6.68E-10 |
| 100 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 7.68E-10 | |
| 150 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5231 | 7.82E-10 | |
| 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 0 |
Results of parameter estimates of Example 1 for G = 1,500.
TABLE 4
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.5024 | 0.5020 | 0.2011 | 0.1462 | 0.0958 | 1.3909 | 1.0393 | 0.7911 | 0.6087 | 0.5671 | 6.10E-04 |
| 100 | 1.4909 | 0.4982 | 0.2001 | 0.1471 | 0.0913 | 1.3968 | 1.0594 | 0.7639 | 0.5954 | 0.6038 | 5.51E-04 | |
| 150 | 1.4930 | 0.5076 | 0.2050 | 0.1453 | 0.1043 | 1.3977 | 1.0479 | 0.8170 | 0.5725 | 0.5632 | 6.00E-04 | |
| 60 | 50 | 1.5001 | 0.4999 | 0.2000 | 0.1502 | 0.1001 | 1.3959 | 1.0468 | 0.7849 | 0.6286 | 0.5216 | 5.94E-07 |
| 100 | 1.5002 | 0.4999 | 0.2001 | 0.1498 | 0.0999 | 1.3960 | 1.0475 | 0.7848 | 0.6282 | 0.5219 | 5.74E-07 | |
| 150 | 1.5000 | 0.5003 | 0.1999 | 0.1498 | 0.0999 | 1.3959 | 1.0467 | 0.7848 | 0.6292 | 0.5194 | 6.15E-07 | |
| 90 | 50 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5229 | 6.59E-10 |
| 100 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 6.67E-10 | |
| 150 | 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 6.59E-10 | |
| 1.5000 | 0.5000 | 0.2000 | 0.1500 | 0.1000 | 1.3960 | 1.0470 | 0.7850 | 0.6280 | 0.5230 | 0 |
Results of parameter estimates of Example 1 for G = 2000.
TABLE 5
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.1894 | 0.8090 | 0.2030 | 0.1721 | 0.0846 | 1.3078 | 0.9609 | 0.8036 | 0.6879 | 0.6185 | 7.91E-04 |
| 100 | 1.1887 | 0.8069 | 0.1964 | 0.1717 | 0.1024 | 1.3088 | 0.9683 | 0.6965 | 0.7246 | 0.5261 | 9.62E-04 | |
| 150 | 1.1867 | 0.7975 | 0.1961 | 0.1844 | 0.1074 | 1.3148 | 0.9610 | 0.7408 | 0.6736 | 0.4660 | 8.57E-04 | |
| 60 | 50 | 1.2029 | 0.8000 | 0.2001 | 0.1829 | 0.1022 | 1.3100 | 0.9597 | 0.7827 | 0.6976 | 0.5053 | 1.34E-05 |
| 100 | 1.2018 | 0.7972 | 0.2015 | 0.1777 | 0.0971 | 1.3077 | 0.9612 | 0.7884 | 0.6951 | 0.5268 | 1.58E-05 | |
| 150 | 1.1984 | 0.7962 | 0.2008 | 0.1790 | 0.0979 | 1.3080 | 0.9620 | 0.7898 | 0.6939 | 0.5613 | 2.19E-05 | |
| 90 | 50 | 1.1987 | 0.7990 | 0.1996 | 0.1796 | 0.0997 | 1.3095 | 0.9583 | 0.7956 | 0.7191 | 0.5277 | 1.14E-05 |
| 100 | 1.1974 | 0.7991 | 0.2004 | 0.1790 | 0.0983 | 1.3075 | 0.9591 | 0.7817 | 0.7046 | 0.5416 | 1.01E-05 | |
| 150 | 1.2008 | 0.7975 | 0.2042 | 0.1780 | 0.1001 | 1.3100 | 0.9559 | 0.7866 | 0.6948 | 0.5428 | 2.01E-05 | |
| 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 0 |
Results of parameter estimates of Example 2 for G = 500.
TABLE 6
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.2096 | 0.8051 | 0.2019 | 0.1866 | 0.0936 | 1.3101 | 0.9696 | 0.8023 | 0.6930 | 0.5116 | 7.02E-04 |
| 100 | 1.1960 | 0.8067 | 0.1974 | 0.1822 | 0.1000 | 1.3079 | 0.9672 | 0.7963 | 0.7429 | 0.5082 | 7.32E-04 | |
| 150 | 1.1969 | 0.8024 | 0.1955 | 0.1762 | 0.1092 | 1.3089 | 0.9584 | 0.8046 | 0.6871 | 0.4465 | 6.98E-04 | |
| 60 | 50 | 1.1998 | 0.8001 | 0.1998 | 0.1802 | 0.0997 | 1.3092 | 0.9590 | 0.7865 | 0.6978 | 0.5212 | 7.53E-07 |
| 100 | 1.2000 | 0.8002 | 0.1998 | 0.1800 | 0.1005 | 1.3094 | 0.9586 | 0.7861 | 0.7009 | 0.5285 | 8.08E-07 | |
| 150 | 1.2001 | 0.8001 | 0.2000 | 0.1802 | 0.1001 | 1.3088 | 0.9593 | 0.7865 | 0.6980 | 0.5252 | 8.10E-07 | |
| 90 | 50 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5229 | 1.45E-09 |
| 100 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7848 | 0.6981 | 0.5231 | 2.16E-09 | |
| 150 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9591 | 0.7851 | 0.6983 | 0.5229 | 4.69E-09 | |
| 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 0 |
Results of parameter estimates of Example 2 for G = 1,000.
TABLE 7
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.1985 | 0.8005 | 0.2013 | 0.1755 | 0.0949 | 1.3134 | 0.9585 | 0.7586 | 0.6795 | 0.5023 | 6.40E-04 |
| 100 | 1.1909 | 0.8036 | 0.1962 | 0.1757 | 0.1137 | 1.3091 | 0.9548 | 0.7750 | 0.6668 | 0.6236 | 6.61E-04 | |
| 150 | 1.2013 | 0.8014 | 0.1918 | 0.1790 | 0.0984 | 1.3135 | 0.9679 | 0.8190 | 0.7242 | 0.6072 | 6.39E-04 | |
| 60 | 50 | 1.2001 | 0.8000 | 0.1999 | 0.1799 | 0.1000 | 1.3089 | 0.9589 | 0.7858 | 0.6976 | 0.5225 | 6.91E-07 |
| 100 | 1.1998 | 0.7997 | 0.2000 | 0.1802 | 0.0999 | 1.3091 | 0.9591 | 0.7847 | 0.6988 | 0.5215 | 6.82E-07 | |
| 150 | 1.1999 | 0.7999 | 0.2002 | 0.1800 | 0.1001 | 1.3091 | 0.9588 | 0.7866 | 0.6979 | 0.5224 | 6.11E-07 | |
| 90 | 50 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 6.71E-10 |
| 100 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 6.73E-10 | |
| 150 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 7.60E-10 | |
| 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 0 |
Results of parameter estimates of Example 2 for G = 1,500.
TABLE 8
| d | p | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 30 | 50 | 1.1992 | 0.7993 | 0.1921 | 0.1931 | 0.1024 | 1.3071 | 0.9623 | 0.8177 | 0.7112 | 0.5600 | 5.57E-04 |
| 100 | 1.1976 | 0.8012 | 0.2052 | 0.1809 | 0.1048 | 1.2983 | 0.9477 | 0.8055 | 0.7149 | 0.5158 | 5.55E-04 | |
| 150 | 1.2013 | 0.8022 | 0.2006 | 0.1794 | 0.1015 | 1.3068 | 0.9566 | 0.7844 | 0.7002 | 0.6340 | 6.02E-04 | |
| 60 | 50 | 1.2001 | 0.8000 | 0.2002 | 0.1802 | 0.0999 | 1.3092 | 0.9590 | 0.7843 | 0.6989 | 0.5230 | 6.11E-07 |
| 100 | 1.1999 | 0.8002 | 0.2000 | 0.1800 | 0.1001 | 1.3090 | 0.9591 | 0.7845 | 0.6993 | 0.5232 | 6.04E-07 | |
| 150 | 1.2001 | 0.8001 | 0.1998 | 0.1802 | 0.0998 | 1.3090 | 0.9594 | 0.7859 | 0.6976 | 0.5214 | 6.33E-07 | |
| 90 | 50 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 6.67E-10 |
| 100 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 6.73E-10 | |
| 150 | 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 6.45E-10 | |
| 1.2000 | 0.8000 | 0.2000 | 0.1800 | 0.1000 | 1.3090 | 0.9590 | 0.7850 | 0.6980 | 0.5230 | 0 |
Results of parameter estimates of Example 2 for G = 2000.
FIGURE 2
FIGURE 3
4 Conclusion
This study exploits a swarm intelligence–based cuckoo search optimization, CSO, heuristic for parameter estimation of power system harmonics. The CSO accurately estimates the amplitude and phase parameters associated with the first, third, fifth, and eleventh harmonic components. Simulation studies conducted on the mean square error–based evaluation metric indicate that the accuracy of the CSO increases with an increase in the generation size, while increasing particle size has not shown a significant rise in the accuracy level. Moreover, the CSO has shown a robust performance in estimating the parameters of power system harmonics for different scenarios of additive white Gaussian noise.
Future studies may investigate the applying of the CSO algorithm for real-time harmonics estimation and for other engineering optimization problems (; ; ; Yang et al., 2020; ).
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.
Author contributions
Conceptualization: NC and MZ; methodology: NM, NC, and MZ; software: NM; validation: ZK, AK, MZ, and NC; resources: AK and CC; writing—original draft preparation: NM; writing—review and editing: NC, ZK, AK, and MZ; supervision: CC and NC; project administration: AM and AA; funding acquisition: AM and AA.
Funding
The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia has funded this project, under grant no. (KEP-MSc: 122-135-1443).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, editors, and reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbdelAtyA. M.FoudaM. E.EltawilA. (2022). Parameter estimation of two spiking neuron models with meta-heuristic optimization algorithms. Front. Neuroinform.16, 771730. 10.3389/fninf.2022.771730
2
AltafF.ChangC. L.ChaudharyN. I.RajaM. A. Z.CheemaK. M.ShuC. M.et al (2022). Adaptive evolutionary computation for nonlinear hammerstein control autoregressive systems with key term separation principle. Mathematics10 (6), 1001. 10.3390/math10061001
3
ArqubO. A.Al-SmadiM. (2020). Numerical solutions of Riesz fractional diffusion and advection-dispersion equations in porous media using iterative reproducing kernel algorithm. J. Porous Media23 (8), 783–804. 10.1615/jpormedia.2020025011
4
ArqubO. A. (2016). Approximate solutions of DASs with nonclassical boundary conditions using novel reproducing kernel algorithm. Fundam. Inf.146 (3), 231–254. 10.3233/fi-2016-1384
5
ArqubO. A. (2020). Numerical simulation of time-fractional partial differential equations arising in fluid flows via reproducing Kernel method. Int. J. Numer. Methods Heat. Fluid Flow.30, 4711–4733. 10.1108/hff-10-2017-0394
6
ArqubO. A. (2018). Numerical solutions for the Robin time-fractional partial differential equations of heat and fluid flows based on the reproducing kernel algorithm. Int. J. Numer. Methods Heat. Fluid Flow.28, 828–856. 10.1108/hff-07-2016-0278
7
BaradaraniF.ZadehM. R. D.ZamaniM. A. (2014). A phase-angle estimation method for synchronization of grid-connected power-electronic converters. IEEE Trans. Power Deliv.30 (2), 827–835. 10.1109/tpwrd.2014.2362930
8
BeleiuH. G.MaierV.PavelS. G.BirouI.PicăC. S.DărabP. C. (2020). Harmonics consequences on drive systems with induction motor. Appl. Sci.10 (4), 1528. 10.3390/app10041528
9
ChaudharyN. I.LatifR.RajaM. A. Z.MachadoJ. T. (2020). An innovative fractional order LMS algorithm for power signal parameter estimation. Appl. Math. Model.83, 703–718. 10.1016/j.apm.2020.03.014
10
ChaudharyN. I.ZubairS.RajaM. A. Z. (2017). A new computing approach for power signal modeling using fractional adaptive algorithms. ISA Trans.68, 189–202. 10.1016/j.isatra.2017.03.011
11
ChaudharyN. I.RajaM. A. Z.HeY.KhanZ. A.Tenreiro MachadoJ. (2021). Design of multi innovation fractional LMS algorithm for parameter estimation of input nonlinear control autoregressive systems. Appl. Math. Model.93, 412–425. 10.1016/j.apm.2020.12.035
12
Cuong-LeT.MinhH. L.KhatirS.WahabM. A.TranM. T.MirjaliliS. (2021). A novel version of Cuckoo search algorithm for solving optimization problems. Expert Syst. Appl.186, 115669. 10.1016/j.eswa.2021.115669
13
DasS. R.RayP. K.SahooA. K.BalasubramanianK.ReddyG. S. (2020). Improvement of power quality in a three-phase system using an adaline-based multilevel inverter. Front. Energy Res.8, 23. 10.3389/fenrg.2020.00023
14
do Nascimento SepulchroW.EncarnaçãoL. F.BrunoroM. (2014). Harmonic state and power flow estimation in distribution systems using evolutionary strategy. J. Control Autom. Electr. Syst.25 (3), 358–367. 10.1007/s40313-014-0110-1
15
ElnadyA.Al-ShabiM.AdamA. A. (2020). Novel filters based operational scheme for five-level diode-clamped inverters in microgrid. Front. Energy Res.8, 11. 10.3389/fenrg.2020.00011
16
Elvira-OrtizD. A.Jaen-CuellarA. Y.Morinigo-SoteloD.Morales-VelazquezL.Osornio-RiosR. A.Romero-TroncosoR. d. J. (2020). Genetic algorithm methodology for the estimation of generated power and harmonic content in photovoltaic generation. Appl. Sci.10 (2), 542. 10.3390/app10020542
17
EnayatiJ.MoravejZ. (2017). Real‐time harmonics estimation in power systems using a novel hybrid algorithm. IET Gener. Transm. &. Distrib.11 (14), 3532–3538. 10.1049/iet-gtd.2017.0044
18
FrancescaG.BirattariM. (2016). Automatic design of robot swarms: Achievements and challenges. Front. Robot. AI3, 29. 10.3389/frobt.2016.00029
19
GudeS.JanaK. C. (2020). Parameter extraction of photovoltaic cell using an improved cuckoo search optimization. Sol. Energy204, 280–293. 10.1016/j.solener.2020.04.036
20
HousseinE. H.HosneyM. E.ElhosenyM.OlivaD.MohamedW. M.HassaballahM. (2020). Hybrid Harris hawks optimization with cuckoo search for drug design and discovery in chemoinformatics. Sci. Rep.10 (1), 14439–14522. 10.1038/s41598-020-71502-z
21
JafarpishehB.MadaniS. M.ShahrtashS. M. (2016). A new DFT-based phasor estimation algorithm using high-frequency modulation. IEEE Trans. Power Deliv.32 (6), 1–2423. 10.1109/tpwrd.2016.2629762
22
JanaG.MitraA.PanS.SuralS.ChattarajP. K. (2019). Modified particle swarm optimization algorithms for the generation of stable structures of carbon clusters, Cn (n= 3–6, 10). Front. Chem.7, 485. 10.3389/fchem.2019.00485
23
KabalciY.KockanatS.KabalciE. (2018). A modified ABC algorithm approach for power system harmonic estimation problems. Electr. power Syst. Res.154, 160–173. 10.1016/j.epsr.2017.08.019
24
KarahanO. (2021). Design of optimal fractional order fuzzy PID controller based on cuckoo search algorithm for core power control in molten salt reactors. Prog. Nucl. Energy139, 103868. 10.1016/j.pnucene.2021.103868
25
LiX.YinM. (2015). Modified cuckoo search algorithm with self adaptive parameter method. Inf. Sci.298, 80–97.10.1016/j.ins.2014.11.042
26
MalikN. A.ChangC. L.ChaudharyN. I.RajaM. A. Z.CheemaK. M.ShuC. M.et al (2022). Knacks of fractional order swarming intelligence for parameter estimation of harmonics in electrical systems. Mathematics10 (9), 1570. 10.3390/math10091570
27
MehmoodA.RajaM. A. Z.ZameerA.ChaudharyN. I. (2021). Design of backtracking search heuristics for parameter estimation of power signals. Neural comput. Appl.33 (5), 1479–1496. 10.1007/s00521-020-05029-9
28
MehmoodA.ChaudharyN. I.ZameerA.RajaM. A. Z. (2020). Novel computing paradigms for parameter estimation in power signal models. Neural comput. Appl.32 (10), 6253–6282. 10.1007/s00521-019-04133-9
29
MontoyaF. G.BañosR.AlcaydeA.Arrabal-CamposF. M. (2019). Analysis of power flow under non-sinusoidal conditions in the presence of harmonics and interharmonics using geometric algebra. Int. J. Electr. Power & Energy Syst.111, 486–492. 10.1016/j.ijepes.2019.04.032
30
PandeyA. C.RajpootD. S.SaraswatM. (2017). Twitter sentiment analysis using hybrid cuckoo search method. Inf. Process. Manag.53 (4), 764–779. 10.1016/j.ipm.2017.02.004
31
PhannilN.JettanasenC.NgaopitakkulA. (2018). Harmonics and reduction of energy consumption in lighting systems by using LED lamps. Energies11 (11), 3169. 10.3390/en11113169
32
RayP. K.SubudhiB. (2012). BFO optimized RLS algorithm for power system harmonics estimation. Appl. Soft Comput.12 (8), 1965–1977. 10.1016/j.asoc.2012.03.008
33
RayP. K.SubudhiB. (2015). Neuro-evolutionary approaches to power system harmonics estimation. Int. J. Electr. Power & Energy Syst.64, 212–220. 10.1016/j.ijepes.2014.07.035
34
RivasA. E. L.da SilvaN.AbrãoT. (2020). Adaptive current harmonic estimation under fault conditions for smart grid systems. Electr. Power Syst. Res.183, 106276. 10.1016/j.epsr.2020.106276
35
SabirZ.RajaM. A. Z.GuiraoJ. L.ShoaibM. (2020). A neuro-swarming intelligence-based computing for second order singular periodic non-linear boundary value problems. Front. Phys.8, 224. 10.3389/fphy.2020.00224
36
SantosE.KhosravyM.LimaM. A.CerqueiraA. S.DuqueC. A. (2020). ESPRIT associated with filter bank for power-line harmonics, sub-harmonics and inter-harmonics parameters estimation. Int. J. Electr. Power & Energy Syst.118, 105731. 10.1016/j.ijepes.2019.105731
37
ShuaiZ.ZhangJ.TangL.TengZ.WenH. (2018). Frequency shifting and filtering algorithm for power system harmonic estimation. IEEE Trans. Ind. Inf.15 (3), 1554–1565. 10.1109/tii.2018.2844191
38
SinghS. K.SinhaN.GoswamiA. K.SinhaN. (2016). Robust estimation of power system harmonics using a hybrid firefly based recursive least square algorithm. Int. J. Electr. Power & Energy Syst.80, 287–296. 10.1016/j.ijepes.2016.01.046
39
SongP. C.PanJ. S.ChuS. C. (2020). A parallel compact cuckoo search algorithm for three-dimensional path planning. Appl. Soft Comput.94, 106443. 10.1016/j.asoc.2020.106443
40
SunJ.AboutaniosE.SmithD. B.FletcherJ. E. (2019). Robust frequency, phase, and amplitude estimation in power systems considering harmonics. IEEE Trans. Power Deliv.35 (3), 1158–1168. 10.1109/tpwrd.2019.2936885
41
SweisH.ShawagfehN.ArqubO. A. (2022). Fractional crossover delay differential equations of Mittag-Leffler kernel: Existence, uniqueness, and numerical solutions using the Galerkin algorithm based on shifted Legendre polynomials. Results Phys.41, 105891. 10.1016/j.rinp.2022.105891
42
Tran-NgocH.KhatirS.De RoeckG.Bui-TienT.WahabM. A. (2019). An efficient artificial neural network for damage detection in bridges and beam-like structures by improving training parameters using cuckoo search algorithm. Eng. Struct.199, 109637. 10.1016/j.engstruct.2019.109637
43
XuL.ChenF.DingF.AlsaediA.HayatT. (2021). Hierarchical recursive signal modeling for multifrequency signals based on discrete measured data. Int. J. Adapt. Control Signal Process.35 (5), 676–693. 10.1002/acs.3221
44
XuL.DingF. (2018). Iterative parameter estimation for signal models based on measured data. Circuits Syst. Signal Process.37 (7), 3046–3069. 10.1007/s00034-017-0705-4
45
XuL.DingF. (2017). Recursive least squares and multi-innovation stochastic gradient parameter estimation methods for signal modeling. Circuits Syst. Signal Process.36 (4), 1735–1753. 10.1007/s00034-016-0378-4
46
XuL.DingF.ZhuQ. (2022). Separable synchronous multi-innovation gradient-based iterative signal modeling from on-line measurements. IEEE Trans. Instrum. Meas.71, 1–13. 10.1109/tim.2022.3154797
47
XuL. (2022). Separable multi-innovation Newton iterative modeling algorithm for multi-frequency signals based on the sliding measurement window. Circuits Syst. Signal Process.41 (2), 805–830. 10.1007/s00034-021-01801-x
48
XuL. (2022). Separable Newton recursive estimation method through system responses based on dynamically discrete measurements with increasing data length. Int. J. Control Autom. Syst.20 (2), 432–443. 10.1007/s12555-020-0619-y
49
XuL.SongG. (2020). A recursive parameter estimation algorithm for modeling signals with multi-frequencies. Circuits Syst. Signal Process.39 (8), 4198–4224. 10.1007/s00034-020-01356-3
50
YangX. S. (2022). MATLAB Central File Exchange. Avaliable at: https://www.mathworks.com/matlabcentral/fileexchange/29809-cuckoo-search-cs-algorithm September 30, 2022).
51
YangX. S.DebS. (2009). “Cuckoo search via Lévy flights,” in Proceedings of the 2009 World congress on nature & biologically inspired computing (NaBIC), Coimbatore, India, December 2009 (IEEE), 210–214. 10.1109/NABIC.2009.5393690
52
YangX. S.DebS. (2014). Cuckoo search: Recent advances and applications. Neural comput. Appl.24 (1), 169–174. 10.1007/s00521-013-1367-1
53
YangX.YangY.LiuY.DengZ. (2020). A reliability assessment approach for electric power systems considering wind power uncertainty. IEEE Access8, 12467–12478. 10.1109/access.2020.2966275
54
YuJ.KimC. H.RheeS. B. (2020). Clustering cuckoo search optimization for economic load dispatch problem. Neural comput. Appl.32 (22), 16951–16969. 10.1007/s00521-020-05036-w
55
ZubairS.ChaudharyN. I.KhanZ. A.WangW. (2018). Momentum fractional LMS for power signal parameter estimation. Signal Process.142, 441–449. 10.1016/j.sigpro.2017.08.009
Summary
Keywords
cuckoo search, swarm optimization, harmonics, parameter estimation, smart grid
Citation
Malik NA, Chang C-L, Chaudhary NI, Khan ZA, Zahoor Raja MA, Kiani AK, Milyani AH and Azhari AA (2022) Parameter estimation of harmonics arising in electrical instruments of smart grids using cuckoo search heuristics. Front. Energy Res. 10:1059132. doi: 10.3389/fenrg.2022.1059132
Received
30 September 2022
Accepted
25 October 2022
Published
15 November 2022
Volume
10 - 2022
Edited by
Khalid Mehmood Cheema, Fatima Jinnah Women University, Pakistan
Reviewed by
Asad Aziz, Southeast University, China
Arshad Rauf, Nanjing University of Aeronautics and Astronautics, China
Omar Abu Arqub, Al-Balqa Applied University, Jordan
Updates
Copyright
© 2022 Malik, Chang, Chaudhary, Khan, Zahoor Raja, Kiani, Milyani and Azhari.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Adiqa Kausar Kiani, adiqa@yuntech.edu.tw
This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.