Abstract
Identifying photovoltaic (PV) parameters accurately and reliably can be conducive to the effective use of solar energy. The grey wolf optimizer (GWO) that was proposed recently is an effective nature-inspired method and has become an effective way to solve PV parameter identification. However, determining PV parameters is typically regarded as a multimodal optimization, which is a challenging optimization problem; thus, the original GWO still has the problem of insufficient accuracy and reliability when identifying PV parameters. In this study, an enhanced grey wolf optimizer with fusion strategies (EGWOFS) is proposed to overcome these shortcomings. First, a modified multiple learning backtracking search algorithm (MMLBSA) is designed to ameliorate the global exploration potential of the original GWO. Second, a dynamic spiral updating position strategy (DSUPS) is constructed to promote the performance of local exploitation. Finally, the proposed EGWOFS is verified by two groups of test data, which include three types of PV test models and experimental data extracted from the manufacturer’s data sheet. Experiments show that the overall performance of the proposed EGWOFS achieves competitive or better results in terms of accuracy and reliability for most test models.
Keywords
Introduction
To mitigate problems associated with the lack of natural resources and to protect the environment, the development of renewable and clean energy has become the focus of current fields of research. Solar energy is generally recognized as a promising resource that can be applicable in the real world because it can be switched to daily electricity through photovoltaic (PV) systems [1, 2].
However, PV systems are easily affected by harsh outdoor environments. Constructing accurate PV models is conducive to promoting the efficiency of converting solar energy into electrical energy [3, 4, 5]. Currently, although some effective PV mathematic models, including double and single diodes, have been proposed [3, 5, 6, 7], the accuracy of these models is restricted by PV parameters, which need to be obtained through identification technology. Therefore, developing an effective parameter identification technology is important.
In addition, the identification of PV parameters is typically regarded as a multimodal problem [3, 7, 8], which is a challenging and difficult optimization problem. According to [7, 8, 9, 10], the traditional deterministic technique and meta-heuristic algorithm are two primary methods for solving PV parameter identification. Although deterministic technology has excellent local search ability on the entire, its disadvantage still may affect its application in PV parameter identification, which requires better proficiency in global exploration [11]. However, with the development of nature-inspired computing [12], many meta-heuristic optimization algorithms that are inspired by nature have appeared and include neural dynamic models [13], genetic algorithms [14], spiral dynamics algorithms [15], water drop algorithms [16], harmony search algorithms [17], and gravitational search algorithms [18]. Because this type of algorithm is not constrained by the form of the problem, it is suitable for complex application problems, such as reducing memory demand [19], renewable energy systems [20], servo systems [21, 22, 23], 3D mesh simplification [24], online scheduling [25], and smart grids [26, 27]. Considering these advantages, many scholars are committed to the development and improvement of meta-heuristic algorithms [28, 29, 30, 31, 32, 33], particularly the application of meta-heuristic algorithms in the parameter identification problem of the PV model. These algorithms include particle swarm optimization (PSO) [34], differential evolution (DE) [35], backtracking search algorithm (BSA) [36], whale optimization algorithm (WOA) [37], Jaya [38], grasshopper optimization algorithm (GOA) [39], equilibrium optimizer (EO) [40], sine cosine algorithm (SCA) [41], and multi-verse optimizer (MVO) [42]. Many other meta-heuristic algorithms have been recommended for PV parameter identification [43], and these meta-heuristics have achieved good results in terms of accuracy and reliability for PV parameter identification. However, most of the existing meta-heuristics still have difficulty determining the global optimal parameters [38]. Thus, there is still a lot of room to develop a potential meta-heuristic algorithm to identify PV parameters effectively.
Recently, a new meta-heuristic method called the grey wolf optimizer (GWO) was proposed in [44]. Due to its efficiency and few adjustment parameters, GWO has been widely used in many fields, including path planning [45], multidimensional knapsack problems [46], feature selection [47], image processing [48], economic load dispatch [49], and other optimization problems [50, 51, 52]. In addition, GWO has been selected to manage PV parameter identification [7], and promising results have been obtained.
However, as a recent meta-heuristic algorithm, GWO still has the potential to enhance accuracy and reliability in the problem of PV parameter identification. Therefore, it is worth considering how to improve GWO to solve the problem of PV parameter identification more effectively. Specifically, the search process of GWO is guided by three levels of wolves and performs exploitation well, but the resulting global exploration ability still has room for improvement during position updating [53, 54]. In addition, although GWO achieves good exploitation ability, it still faces the problem of insufficient reliability when solving PV parameter identification. Therefore, we focus on the exploration and exploitation mechanism of GWO to improve its PV parameter identification ability. Overall, no algorithm can be universally considered the best optimization algorithm according to [55], but developing optimization algorithms to solve the problem of PV parameter identification more effectively should be encouraged.
Therefore, an enhanced grey wolf optimizer with fusion strategies (EGWOFS) is constructed in this study to obtain the PV parameters more effectively. In EGWOFS, first, to improve the insufficient exploration ability of the GWO algorithm when solving PV parameter identification, a modified multiple learning backtracking search algorithm (MMLBSA) is constructed to enhance the exploration performance of GWO, and then, the accuracy of the identified PV parameters can also be improved. Second, although the global exploration of the GWO algorithm has been improved, the reliability of GWO improved by MMLBSA is still insufficient at obtaining the PV parameters; thus, we design a dynamic spiral updating position strategy (DSUPS) to improve the reliability of GWO by promoting the exploitation performance. Finally, to measure the effectiveness of EGWOFS, experiments are performed out with two groups of data, including three types of benchmark PV models and measured data extracted from the manufacturer’s data sheet. Results comprehensively demonstrate that the overall performance of the proposed EGWOFS achieved better or competitive results in terms of accuracy and reliability for most PV models.
The four contributions of this study are as follows:
EGWOFS is designed to better identify PV parameters in terms of accuracy and reliability; MMLBSA is constructed to promote the exploration performance of GWO; DSUPS is proposed to ameliorate the exploitation performance; The performance of EGWOFS is tested with comprehensive experiments.
The remainder of this paper is organized as follows. Section 2 introduces the solar cell model, PV module model and their parameter optimization. Section 3 describes the fundamentals of GWO. Section 4 describes the proposed EGWOFS. Section 5 describes the results and analysis. Section 6 provides conclusions of this study.
According to the literature [8, 36, 56, 57, 58, 59], there are three primary PV models: the double diode model (DDM), PV module model, and single diode model (SDM). Their parameter identification problems have been transformed into an objective function that can be optimized by algorithms. The three models are illustrated in the following subsections.
Solar cell model
According to [8, 36, 56, 57, 58, 59], two models including SDM and DDM, are described in the following two subsections.
SDM
The SDM has five unknown parameters, and its current output can be handled by the following formulae:
where output current
The DDM has seven unknown parameters, and its output current is described by the following formulae:
where the diffusion current
The current output of the single diode PV module is calculated using the following formula [8, 36, 56, 57, 58, 59]:
where the number of both parallel and series solar cells are represented by two variables
Objective function
Identifying PV parameters needs to be converted into a problem that can be optimized by an algorithm. Therefore, Eq. (14), which has been widely used in the literature [8, 36, 56, 57, 58, 59], is considered to be the objective function in this study. The calculation result of the objective function primarily depends on the error functions shown in Eq. (11), Eq. (12), and Eq. (13):
where Eq. (11), Eq. (12), and Eq. (13) are the error functions of the SDM, DDM, and single diode PV module model (
where
In addition to the above objective function, parameter variables are an important part of the optimization problem.
For SDM, five parameter variables need to be obtained by the meta-heuristic algorithm in Eq. (11):
For DDM, seven parameter variables need to be obtained by the meta-heuristic algorithm in Eq. (12):
For the PV module, five parameter variables need to be obtained by the meta-heuristic algorithm in Eq. (13):
Constraints
Constraints are another core component of the optimization problem. The constraint intervals of the specific parameter variables for the three models are listed in Table 1 and were used in [7, 60, 61]. More details are described below.
Search interval of identified parameters on three types of PV model
Search interval of identified parameters on three types of PV model
Social hierarchy
According to the literature [44], the ideology of the grey wolf optimizer (GWO) is primarily dependent on the daily capturing behaviour of grey wolves. The grey wolf population is generally delimited into four types: the first three levels grey wolves
Encircling prey
The behavior of encircling prey is primarily represented by the following model [44]:
where
The flowchart of PV parameter identification based on EGWOFS.
At this stage, the population depends on three levels of grey wolves to obtain the prey, and the positions of the remaining grey wolves are renewed via the prey position obtained by
where the vector
In GWO, according to [44], the process of exploitation is achieved by the change of some variables, and the value of vector
Prey search (exploitation)
Based on [44], although GWO primarily uses three levels of grey wolves to guide other grey wolves to seek prey, it is still insufficient for the exploration effect of GWO. Therefore, the deficiency of GWO is enhanced by adjusting the condition
Enhanced grey wolf optimizer with fusion strategies (EGWOFS)
Because the original GWO algorithm has deficiencies when solving PV parameters, an enhanced grey wolf optimizer with fusion strategies (EGWOFS) is proposed to more effectively identify PV parameters. A flowchart of PV parameter identification based on EGWOFS is shown in Fig. 1. In EGWOFS, we improved the process in two stages: global exploration (GE) and local exploitation (LE), where GE is achieved by reconstructing a modified multiple learning backtracking search algorithm (MMLBSA), and LE is realized by designing a dynamic spiral updating position strategy (DSUPS). Next, we examine how to better and effectively integrate the MLBSA and GWO, and fully utilize their respective advantages to compensate for each other’s shortcomings.
MMLBSA
In GWO, the search process primarily depends on three levels of grey wolves
However, solving this problem is more difficult than it appears. Directly combining MLBSA with GWO cannot fully utilize the exploration proficiency of MLBSA and may reduce the accuracy of PV parameter identification. Therefore, the proposed motivation is to flexibly integrate both MLBSA and GWO to better promote the exploration performance of GWO. To solve this problem, we primarily make the following modifications to the original MLBSA.
First, according to MLBSA [36], the historical population, current population, and best individual of the current population are all considered in Eq. (28), which is important for the exploration and convergence performance of MLBSA. However, the information of the current worst individual is still ignored. Although the current best individual can promote convergence speed, it may not help with extensive exploration under the leadership of the best individual. Therefore, we reconstruct Eq. (28) into Eq. (29) and call it the modified multiple learning strategy (MMLS). Here, the historical population, current population, and best and worst grey wolves of the current population are all considered during the exploration. Next, the greedy selection strategy used in MLBSA is selected to renew the population, and the updated population is sorted in ascending order so that we can obtain the first three best individuals to update the original first three levels of grey wolves.
Second, according to MLBSA [36], the current population has better individuals than the historical population, which means that the information diversity of the current population can affect the evolutionary trend of the population. Although MMLS has been considered to ameliorate the convergence speed and exploration performance of GWO, the quality of the current population also must be considered. Therefore, the EMBCLS of MLBSA is also employed to update the current best grey wolf. In addition, to minimize the impact on the accuracy and reliability of parameter identification and reduce unnecessary computational cost, we design diversified constraints shown in Eq. (31) to control the execution times of EMBCLS. Next, we perform the elite mechanism used in MLBSA to update the worst individual of the current population once there is a better solution generated by Eq. (32). Then, the updated population is sorted in ascending order so that we can obtain the first three best individuals to update the original first three levels of grey wolves.
According to the above description, the related models of the MMLBSA are as follows.
where
After these improvements, the MMLBSA can better assist GWO at improving its global exploration capability and obtain an accurate solution when solving PV parameter identification. However, this alone is insufficient: if the search mechanism of GWO is directly used to update the population obtained by MMLBSA, it may affect the reliability of the solution. Therefore, continuing to effectively combine MMLBSA and the search mechanism of GWO to improve the reliability of the solution is more difficult than it seems. To meet this challenge, we propose a dynamic spiral position update strategy (DSPUS) to strengthen the local search around the first three grey wolves. In DSPUS, there are three primary steps as follows.
First, inspired by the spiral position update strategy Eq. (36) of the whale optimization algorithm (WOA) [62], which has good flexibility and local search ability, the spiral position update strategy may be considered a promising method to further enhance the local search of GWO. However, if we directly employ Eq. (36) instead of updating mechanisms Eq. (22), Eq. (23), and Eq. (24), the three updating mechanisms of the first three wolves will become the same. This result means that the search of GWO will be guided only by the first-best wolf instead of the first three best wolves, which may easily affect its exploration performance and may also change the basic characteristics of GWO. Therefore, we modify Eq. (36) into the new Eq. (40), Eq. (42), and Eq. (44) to update the first three grey wolves.
Second, if we only employ Eq. (50) as the local search mechanism, this process does not fully consider the exploitation capabilities of Eq. (40), Eq. (42), and Eq. (44). Inspired by the hierarchy of GWO and the group search of improved JAYA (IJAYA) [8]. We design different updating mechanisms for different types of wolf individuals. Eq. (51), which is based on the idea of the DE mutation operator [63], is designed to update individual
Third, we change the range of
Finally, adjust the boundary of the new individual generated by any one of Eq. (51), Eq. (53), and Eq. (55), and calculate the fitness of new individual. When the new individual is better than the current individual, the new individual is accepted. Then, the new population is sorted to obtain the first three best individuals, and the original first three grey wolves are updated with them.
Under the close cooperation of the above steps, the reliability of the solution can be improved.
where
Related tests and results are performed and analyzed to demonstrate the effectiveness of the proposed EGWOFS in this section. First, the parameter settings of related algorithms are introduced, and two popular PV models are employed to test the generalization ability of PV parameter identification for the proposed EGWOFS. Then, various statistical indicators and visualizations are used to show the effect in Subsection 5.1. Second, Subsection 5.2 compares the proposed EGWOFS with more PV models and established algorithms reported in the literature to verify its advanced performance in PV parameter identification. Third, Subsection 5.3 verifies the effectiveness of the two introduced mechanisms of the proposed EGWOFS. Fourth, EGWOFS is applied to measured data extracted from the manufacturer’s data sheet at different irradiance conditions to test its practical application ability and reliability in Subsection 5.4. The test data of the SDM, DDM and PV module model are obtained from [60, 61, 68, 69] and are widely employed to prove the proficiency of meta-heuristic algorithms in solving the PV parameter identification problem [7, 36, 38, 70], and the measured data of the Multi-crystalline KC200GT are obtained from [71].
Test on two popular PV models
Experimental settings
Set parameters for seven algorithms
Set parameters for seven algorithms
Two types of PV models are first used to verify the accuracy and reliability of the proposed EGWOFS in solving PV parameter identification within an acceptable convergence speed in Section 5.1.
For the evaluation methods, we first use the average (Mean), minimum (Min), and maximum (Max) of the RMSE values and curve characteristics to test the accuracy. Then, we use both boxplot and standard deviation (SD) to verify the method’s reliability. The statistics of four indicators (‘Mean’, ‘Min’, ‘SD’ and ‘Max’) are dependent on the best RMSE gained through each algorithm with 30 independent runs. The experimental results in boldface indicate the best effects gained by related algorithms for a certain indicator. In addition, the convergence speed is confirmed by a curve comparison of each algorithm in the convergence graphs. Finally, Friedman’s rank test is performed to demonstrate the overall ranking performance of EGWOFS.
Moreover, to provide a fair comparison of related algorithms, the number of evaluations of both the proposed EGWOFS and comparison algorithms, including GWO, MVO [72], PSO [73], WOA, Salp Swarm Algorithm (SSA) [74], and EO [75], are all set to 48,000. The population size of each algorithm is 20. The other parameter settings primarily depend on the original literature or the suggestions of the relevant literature. In Table 2, for EGWOFS,
Statistical results of seven algorithms on SDM
Statistical results of seven algorithms on SDM
The 
In Table 3, the proposed EGWOFS achieves the best results in all four indicators for SDM. In particular, six comparison algorithms are much worse than the proposed EGWOFS in two indicators ‘Max’ and ‘Std’. The primary reason for these results is that the six comparison algorithms can easily fall into the local optimum prematurely and lack of the adequate exploitation capacity during the search for five optimal PV parameters, which can easily lead to poor accuracy of the maximum value and unstable solutions. While the proposed EGWOFS has a better global exploration ability for optimal SDM parameters, primarily due to the effectiveness of the MMLBSA mechanism in assisting EGWOFS in globally searching SDM parameters.
The curve feature of Fig. 2 indicates that the data gained by EGWOFS are generally coincident with the real data on SDM, which also verifies that the PV parameters identified by the proposed EGWOFS have high accuracy on SDM, primarily due to the effectiveness of MMLBSA mechanism in improving the global exploration capability of EGWOFS.
The boxplot in Fig. 3 shows that the smallest fluctuation range of the best RMSE is achieved by EGWOFS, which means that EGWOFS has the best reliability among the 7 algorithms on SDM, due to its better exploitation ability to search SDM parameters.
The fluctuation comparison of best RMSE in 30 runs of seven algorithms on SDM.
The convergence effect of seven algorithms on SDM.
Figure 4 shows that EGWOFS achieves the best convergence speed among the seven algorithms for SDM, primarily due to the EMBCLS employed by EGWOFS, in which the search process is performed under the guidance of elite individuals to improve the convergence speed. In addition, we find that EGWOFS, PSO and SSA have a better ability to jump out of the local optimum than the other four algorithms on SDM, which is primarily due to their better global search ability.
Therefore, compared with the original GWO algorithm, the global search ability and convergence speed of the proposed EGWOFS have been markedly improved, which also demonstrates the effectiveness of the improved mechanism.
In Table 4, the proposed EGWOFS achieves the best results in all four indicators for the DDM. The values obtained by the six comparison algorithms are much worse than those obtained by the proposed EGWOFS on the three indicators ‘Max’, ‘Mean’ and ‘Std’. Compared with the SDM, there are more parameters to be identified by the algorithms on the DDM, and the six comparison algorithms have insufficient global exploration and local exploitation in solving the parameter identification of the DDM, which easily leads to their insufficient accuracy and stability. While the introduced MMLBSA and DSUPS can enhance the ability of the EGWOFS to globally explore and locally exploit DDM parameters, respectively. So the proposed EGWOFS has a better ability to jump out of the local optimum on the DDM.
Statistical results of seven algorithms on DDM
Statistical results of seven algorithms on DDM
The 
Figure 5 shows that EGWOFS has produced satisfactory data that are generally coincident with real data on the DDM, which means that the proposed EGWOFS achieves high accuracy in PV parameter identification on the DDM. Because the MMLBSA is conducive to improving the global search ability of the proposed EGWOFS, the accuracy of the solution can be improved.
The boxplot of Fig. 6 shows that the smallest fluctuation range of the best RMSE values is achieved by EGWOFS among the 7 algorithms, which means that the proposed EGWOFS has the best reliability among the 7 algorithms on the DDM. Because the DSUPS mechanism improves the ability of the EGWOFS to locally search the optimal solution of the DDM, the reliability of the solution is better enhanced.
The fluctuation comparison of best RMSE in 30 runs of seven algorithms on DDM.
The convergence performance of each algorithm on DDM.
Figure 7 shows that EGWOFS gains better convergence precision than the other six algorithms on the DDM, while the WOA, EO and PSO achieve faster convergence speeds than EGWOFS, SSA, GWO and MVO on the DDM, and the convergence speeds among WOA, EO and PSO are relatively close. The primary reason is that the DDM needs to identify more parameters than SDM, which makes it more difficult for the algorithms to search the optimal solution. In addition, Fig. 7 shows that the proposed EGWOFS has the best ability to jump out of the local optimum. Therefore, the proposed EGWOFS can obtain the best convergence precision but requires a longer convergence time to globally search the optimal solution than the WOA, EO, and PSO on the DDM.
According to the above analysis, the proposed EGWOFS achieves the best convergence accuracy within the accepted convergence speed on the DDM.
To further verify the advantage of the EGWOFS in terms of the overall ranking, we perform a well-known statistical hypothesis test called Friedman’s rank test. In addition, based on the previous two test cases, the specific Friedman ranking test results of EGWOFS and six comparison algorithms on two PV models are shown in Table 5.
Ranking result obtained by the Friedman test of seven algorithms on two PV models
Ranking result obtained by the Friedman test of seven algorithms on two PV models
Comparison with established algorithms on SDM
