Abstract

Correct extraction of the equivalent circuit model parameters of photovoltaic modules is of great significance for power prediction, fault diagnosis, and system optimization of photovoltaic power generation systems. Although there are many methods developed to extract the equivalent circuit model parameters of the photovoltaic module, it is still challenging to ensure the stability and operational efficiency of the extract method. In order to effectively extract the parameters of photovoltaic modules, this paper proposes a hybrid algorithm combining analytical methods and differential evolution algorithms for the extraction parameters of PV module. Firstly, the analytical method is applied to simplify the equivalent circuit model and improve the efficiency of the algorithm. Then, the adaptive algorithm is used to adjust the parameters of the differential evolution algorithm. Through the algorithm proposed in this paper, the parameters of the equivalent circuit model of the photovoltaic module can be extracted by the open-circuit voltage, short-circuit current, and maximum power point current and voltage provided by the manufacturer. The proposed method is applied to the extraction of the parameters of the dual-diode equivalent circuit model of different types of photovoltaic modules. The reliability and computational efficiency of the proposed algorithm are verified by comparison and analysis.

1. Introduction

With the emphasis on environmental and energy issues, photovoltaic power generation technology is widely used. In order to predict the power generation of photovoltaic power generation systems [15], improve the performance of photovoltaic modules, and perform fault diagnosis [610], extracting model parameters based on an equivalent circuit model of photovoltaic modules has become one of the research hotspots [1113].

The simulation models of solar cells currently widely used include single-diode models and two-diode models [14]. Traditional photovoltaic module modeling uses single-diode models. In recent years, with the deepening of research, researchers have found that the two-diode model takes into account the composite effect of the electrically neutral region and the space charge region, revealing the working mechanism of photovoltaic cells, making the two-diode model more accurate than the single-diode model [15]. However, the equivalent circuit model of the PV module based on the two-diode model is implicitly nonlinear. In equation, it is difficult to find the analytical solution by a traditional algebraic method. In recent years, researchers have proposed a large number of modeling methods of this problem. These methods can be divided into two categories: one is the numerical solution method, and the analytical solution to the model is obtained through some empirical formulas and mathematical techniques. It is shorter, because the method usually ignores some parameters or uses empirical formulas to obtain some parameters. If the parameters are improperly set in the application process, it will affect the accuracy of the modeling [1619]. The other type uses the intelligent optimization algorithm, such as the differential evolution (DE) algorithm, particle swarm optimization (PSO) algorithm, and genetic algorithm (GA). These methods do not approximate the parameters during use, which can improve the accuracy of the solution [2027]. However, these intelligent optimization algorithms are random optimization algorithms. If the algorithm parameters are set incorrectly, it will affect the stability and robustness of the solution. In addition, these algorithms take a long time to solve high-dimensional optimization problems. There are also two sources of data for PV module modeling. One is the data of a PV module obtained under experiment under certain conditions; the other is the standard test condition (STC) provided by the manufacturer: short-circuit current, open-circuit voltage, and voltage and current at the maximum power point. The first type of data comes from experiments. The established model is close to the actual working conditions, but the experimental workload is large, and the test process data will error due to the limitations of the test conditions, resulting in a large model error. The second is data. Standard test data from manufacturers are reliable data source, have small modeling workload, and are widely used in engineering, but this method has a small amount of data; if the modeling method is not appropriate, it will lead to large model error.

To find a fast, simple, and accurate method for modeling solar cells, Chin et al. [28] proposed a hybrid calculation method that combines numerical solutions with intelligent optimization algorithms to extract the parameters of the photovoltaic module circuit model using standard test data provided by the manufacturer. The solution method is used to obtain the photogenerated current in the solar cell model under standard test conditions, the two equivalent diode reverse saturation currents, and the parallel resistance. Then, the differential evolution algorithm is used to find the ideal factor and series resistance of the two diodes. This method combines two modeling methods to achieve fast and accurate modeling of solar cells. However, in order to maintain the accuracy of the modeling, the method does not make any approximations and assumptions in the numerical solution method, which leads to the complexity of the explicit expression of the parameters, which increases the difficulty of modeling and increases the calculation time. In addition, the evolutionary algorithm has the advantages of easy implementation and optimization for multipeak functions. However, in the evolutionary process of random search, the phenomenon of “premature maturity” is easy to occur, which leads to the local best, which leads to the decline of modeling accuracy.

Tao et al. [29] used the adaptive differential evolution algorithm to extract the parameters of the PV module based on the two-diode model, but it is modeled by the data obtained by the experiment, and the dimensionality of the solution is not caused by the dimension reduction processing of the model. It is higher, is more computational, and has time-consuming modeling. The adaptive differential evolution algorithm used in this paper only makes the scaling factor adaptively adjust with evolutionary algebra. Although this method can avoid the phenomenon of “premature,” the search efficiency and robustness of the algorithm are not improved much.

In order to maintain the advantages of the differential evolution algorithm and suppress the “premature maturity” phenomenon, researchers at home and abroad have improved the differential evolution algorithm. The improved method can be divided into two types, one is the improvement of the algorithm parameters and steps, such as Qin et al. [30] replacing the contraction factor of the algorithm with a Gaussian random number and randomly selecting an individual to perform mutation operations to improve the robustness and search efficiency of the algorithm. Vo-Duy et al. [31] in the variation operation, according to the average adaptation of the population to the optimal value, determine the scaling factor. The population adds partial perturbations in the later stages of evolution, maintaining population diversity and avoiding “premature maturity.” Another type is to add disturbance factors outside the population, maintain population diversity, and avoid falling into local parts, which is the most advantageous. Tvrdik and Krivy [32] used the best clustering method to select indivitwos involved in evolution, selected individual from different categories to form new populations of the iterative process, and randomly sorted the indivitwos in the cluster center. It searches for the global best. Ayala et al. [33] used the fusion of population information and manual intervention selection process to automatically identify the best search path. Qin and Suganthan [34] combine the above two improved methods, using improved crossover strategy and adaptive adjustment control parameters, to achieve global optimization. In order to improve the performance of the differential evolution algorithm, the researchers improved the traditional differential evolution algorithm in terms of parameter setting and mutation strategy and proposed many improved differential evolution algorithms [33, 3543]. In this paper, the parameter adaptive method proposed by Civicioglu and Besdok is used to adjust the parameters in the algorithm [33, 3543].

The aim of this paper is to develop a PV module parameter extraction algorithm that combines the analytical method with the improved differential evolution algorithm. The algorithm can accurately extract the PV module parameters based on the dual-diode circuit model and has higher computational efficiency. In the proposed algorithm, the equivalent circuit model is simplified by the analytical method to improve the efficiency of the algorithm. Then, the adaptive algorithm is used to adjust the parameters in the differential evolution algorithm to avoid the algorithm falling into local optimum and improve the convergence speed. Finally, using the above algorithm and the open-circuit voltage, short-circuit current, and maximum power point current and voltage of the photovoltaic module provided by the manufacturer, the equivalent circuit model parameters of the photovoltaic module can be extracted. The proposed method is applied to the extraction of dual-diode equivalent circuit model parameters of different types of photovoltaic modules. The comparison and analysis indicate that the proposed algorithm has high reliability and competitive computational efficiency.

The main contributions of this paper are as follows: (1)A hybrid algorithm combining analytical methods and differential evolution algorithms is developed for PV parameter estimation. In this method, the analytical method is applied to simplify the equivalent circuit model and improve the efficiency of the algorithm(2)The proposed algorithm employs an adaptive method to adjust the parameters in the differential evolution algorithm to avoid the algorithm getting into the local best and speed up the convergence of the algorithm(3)The proposed algorithm is applied to extract parameters of three different PV modules. The obtained results of the method proposed in this paper are compared with well-established algorithms to confirm its effectiveness

The rest of the paper is organized as follows: Section 2 introduces the photovoltaic module model based on two diodes. Section 3 describes PV module modeling and parameter optimization problem. Section 4 introduces the principles of the improved differential evolution algorithm. Section 5 presents experimental results and discussion. Section 6 concludes this paper.

2. The Photovoltaic Module Model Based on Two Diodes

The model of the photovoltaic module based on the two diodes is shown in Figure 1. Its mathematical model is as shown in equation (1). where is the output current (A) for the photovoltaic module; is the output voltage (V) for the photovoltaic module; is the photogenerated current (A); and are the two equivalent diode reverse saturation currents (A), respectively; and are the ideality factor of the two-diode model; is the equivalent series resistance of the photovoltaic cell (); is the equivalent parallel resistance of the photovoltaic cell (); and is the thermal temperature coefficient of the photovoltaic cell; its equation is where is the Boltzmann constant ( J/K), is the absolute temperature (K) for the photovoltaic cell, and is the electron charge ().

The short-circuit current , open-circuit voltage , and maximum power point voltage and current at different temperatures and different solar radiation intensities can be obtained by the parameters of the standard test conditions provided by the manufacturer and according to where , , , , , and are the short-circuit current, open-circuit voltage, maximum power point current, maximum power point voltage, solar radiation intensity, and temperature measured under standard test conditions; , , , and are the short-circuit current, open-circuit voltage, maximum power point current, and maximum power point voltage; and and are the radiation intensity and temperature of the photovoltaic module. The manufacturer does not provide the seven parameters , , , , , , and contained in equation (1), so the purpose of modeling the PV module was to estimate the above seven parameters from the data provided by the manufacturer. The seven parameters in (1) make the voltage-current curve of the PV module obtained by equation (1) consistent with the measurement curve provided by the manufacturer.

3. PV Module Modeling and Parameter Optimization

In this paper, using the test data provided by the manufacturer, combined with the numerical solution method and differential evolution algorithm, a method of modeling and parameter optimization of photovoltaic cells using only PV module test parameters is proposed.

First, the manufacturer provides the short-circuit status point data into formula (1): where and are the two-diode quality factors in the photovoltaic cell circuit models under standard test conditions, is the equivalent series resistance of the photovoltaic cell under test conditions, is the equivalent parallel resistance of the photovoltaic cell under the test conditions, and is the thermal temperature coefficient of the photovoltaic cell under test conditions; the mathematical expression is

Normally, , , then equation (1) can be simplified to

Similarly, it is available to substitute the open-circuit voltage state point (, 0) into equation (1):

Set

According to the analysis of Yahya-Khotbehsara and Shahhoseini [44], in general, is three to four times ; it can be expressed by the following formula:

Substituting equation (13) into equation (10) and because of the expression that can be obtained

Similarly, substituting the maximum power state point (, ) into equation (1) is available:

At this point, there are only three unknown parameters in the PV module model: , , and . Effectively reduce the dimension of the differential evolution algorithm to solve the problem. In order to use the differential evolution algorithm to solve the remaining unknown parameters, it is necessary to define a suitable objective function. By observing the curve of the PV module, it can be found that there is a “knee point,” that is, the maximum power point, and the position of the maximum power point is affected by three parameters: , , and , so the PV of the PV module can be optimized by optimizing three parameters. The gradient of the curve at the maximum power point is minimal. According to the above analysis, the output power of the maximum power state point is 0:

It can also be written as

Therefore, the objective function in the differential optimization algorithm is defined as where

When applying the differential evolution algorithm to search for, three parameters, in order to ensure that the solution is a feasible solution, the parameter changes must be constrained. According to the literature [11], the constraints are as follows:

When , or , or , or , or , where and are the minimum and maximum values of , respectively. If the parameter search range exceeds the above constraint range, the objective function will be set to . Figure 2 is the algorithm flow chart, where VTR is the preset value of the objective function.

4. Improve Differential Evolution Algorithm

The differential evolution algorithm is an intelligent optimization algorithm with simple programming and fast convergence. It has strong adaptability to multipeak function optimization problems. It has been optimized in mechanism, image recognition, target tracking, and so on. The field is widely used. With the deepening of human research, the DE algorithm is also applied to the extraction of PV module parameters. However, since the differential evolution algorithm is a random self-heuristic search algorithm, in high-dimensional and multipeak problems, the algorithm will appear “premature,” thus falling into the local best. The basic DE algorithm flow is as follows.

4.1. Initialization

The DE algorithm first randomly generates NP -dimensional decision vectors, where the th vector is expressed as where is the current algebra and is the dimension of the vector.

Each element in the vector is randomly generated according to where and are the upper and lower bounds of the search space, respectively, and is a random number that is uniformly distributed between 0 and 1.

4.2. Mutation Operation

Individuals are randomly selected in the current population, and the mutated individuals are generated by the scaling operation, thereby realizing the search for different regions in the feasible domain. The mutation formula is where , , and are random integers, and .

4.3. Cross Operation

In order to maintain the diversity of the population, the mutated individuals are exchanged with the original individual information to generate new individuals. The crossoperating operator can be expressed as where is a random number in , ensuring that at least one of the candidate individuals is from the information of the mutated individual, preventing the occurrence of invalid evolution. At the same time, the individual values are also constrained in the crossover operation. Avoid crossing the boundary and producing nonfeasible solutions.

4.4. Select Operation

Comparing the appropriateness of the original individual and the candidate individual, select the moderately good individual to enter the next iterative process, so as to optimize the decision variable. The selection operator is as in

Although the basic differential evolution algorithm has a fast convergence speed, its stability and robustness are poor. In order to maintain the fastness of the parameter optimization process and maintain good stability, literature [45] proposed an adaptive differential evolution algorithm, in which the scaling factor, crossover probability, and optimization strategy can be adaptively selected according to the objective function information to improve the performance of the algorithm. The scaling factor and crossfactor generation probability of each individual obey the Cauchy distribution and the normal distribution. The parameter success rate of the last iteration adjusts the control parameters of the probability distribution model so that the parameters are adaptively adjusted and the performance of the differential evolution algorithm is improved.

The specific method of parameter adaptive adjustment is as follows.

For an individual, its crossover probability follows a normal distribution:

In the formula, the effect of and is to limit to , the initial value of is 0.5, and it is updated after each iteration: where is a successful set of crossover probabilities, is the arithmetic mean operation, and is the weighting factor. The scaling factor is generated as follows:

In the formula, the initial value of is set to 0.5 and is updated as follows after each iteration: where is the Lehmer mean, calculated as

5. Experimental Results and Discussion

In order to effectively extract the parameters of the PV module, the parameter setting of the differential evolution algorithm is very important. In Reference [28], the scaling factor and the initial value of the crossover probability are set to 0.5 and 0.8, respectively. The choice of the crossover strategy is in the literature [28]. For the adaptive selection method, the search range of decision variables is as follows: , , and ; the size of the population is set to 30; and the maximum number of iterations is set to 1000.

In order to verify the proposed method, three different photovoltaic modules were selected, including monocrystalline silicon components, polycrystalline silicon components, and thin film components. The parameters were extracted and compared by different methods. The test data was obtained from the manufacturer’s instructions. The parameters of the three components provided by the manufacturer under STC are shown in Table 1. Under the MATLAB platform, the parameters in the diode model corresponding to the three components are extracted by the differential evolution (DE) algorithm [28] and the improved differential evolution (IDE) algorithm proposed in this paper, respectively. In order to make a fair comparison between the algorithm proposed in this paper and the algorithm in the literature [28], the two algorithms were run independently 30 times each to obtain their statistical results. Table 2 shows the statistical results of the experiment. The convergence curves of the two algorithms are shown in Figure 3. The parameters are shown in Tables 35. These tables show that the method used in this paper can accurately solve the parameters of monocrystalline silicon components, polycrystalline silicon components, and thin film components. It can be seen from Table 6 that the average absolute value of SP75 monocrystalline silicon modules extracted by this method under STC is absolute. The error is 0.015, the average absolute error of S25 polysilicon component under STC is 0.013, the average absolute error of ST36 thin film component under STC is 0.026, and the average absolute error is within 0.05. The calculation results of the above three components are, respectively, obtained. Substituting (1), using the Newton iteration method, can solve the characteristic curve of PV modules under different temperatures and different light intensities as shown in Figures 49. It can be seen from the figures that the parameters extracted by the literature [28] and the proposed method under the STC can well fit the output characteristic curves of the three PV modules, but for the method proposed in this paper, the accuracy of the fitting is higher than that of the literature [28].

Set the absolute error to where is the experimental data and is the data obtained by simulation; the mean absolute error (MAE) expression is where is the number of sampling points. The average absolute error of the photovoltaic modules in various environmental conditions is calculated according to equation (36) as shown in Table 6. It can be seen from the table that the average absolute error simulated by the parameters extracted by the proposed method under different environmental conditions is close to that of the literature [28], which proves the effectiveness of the proposed method.

With the CPU: Intel [email protected] GHz, memory: 12 GB computing environment, using the method proposed in this paper and the method of [28] to extract the two diodes of SP75, S25, and ST36 three photovoltaic modules under different environmental conditions, the time taken for the model parameters is shown in Table 7. It can be seen from the table that the time taken to extract parameters using this method in the same computing environment is smaller than that in the literature [28], indicating that the calculation speed can be improved after appropriate approximation.

Figures 1012 shows the absolute error comparison of the parameters of three different modules under STC using the literature [28] and the extraction method of this paper. It can be seen from the figures that for the SP75 monocrystalline silicon module and the S25 polysilicon module, the absolute error of the two methods is basically the same in the constant current region of the photovoltaic module, but in the constant voltage region, the accuracy of the proposed method is higher than that of the literature [28] method. However, for the ST36 thin film module, the method of [28] is higher than the method proposed in this paper, indicating that the conclusion of Yahya-Khotbehsara and Shahhoseini [44] for three to four times is not applicable to the thin film module.

6. Conclusion

This paper proposes a hybrid photovoltaic module parameter extraction algorithm that combines the analytical method with the improved differential evolution algorithm. It only needs four parameters that are the open-circuit voltage, short-circuit current, and maximum power point current and voltage of the photovoltaic module provided by the manufacturer; then, the parameters of the photovoltaic module based on the dual-diode equivalent circuit model are extracted. In the proposed algorithm, the equivalent circuit model is simplified by the analytical method to improve the efficiency of the algorithm. Then, the adaptive algorithm is used to adjust the parameters in the differential evolution algorithm to avoid the algorithm falling into local optimum and improve the convergence speed. In this paper, a comprehensive experimental test is carried out on the algorithm to study the performance of the algorithm on the parameter extraction of different types of photovoltaic modules. Compared with other recent methods, experimental and statistical analysis proves the superiority of the proposed algorithm in terms of accuracy, reliability, and computational efficiency.

Data Availability

The experimental data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

H.W. and Z.S. conceived the paper and designed and performed the simulations; H.W. wrote and revised the paper; H.W. and Z.S. finally analyzed the data.

Acknowledgments

This research was supported by the Special Financial Aid to Post-Doctor Research Fellow in Chongqing (grant Xm2014086); Science and Technology Research Project of the Chongqing Municipal Education Commission (grant KJ131321), and Science and Technology Research Project of Yangtze Normal University (grant 2013XJ2D004).