Open access peer-reviewed chapter

Static and Dynamic Photovoltaic Cell/Module Parameters Identification

Written By

Sid-Ali Blaifi and Bilal Taghezouit

Submitted: April 18th, 2019 Reviewed: August 30th, 2019 Published: April 1st, 2020

DOI: 10.5772/intechopen.89449

Chapter metrics overview

749 Chapter Downloads

View Full Metrics

Abstract

The accurate parameters extraction is an important step to obtain a robust PV outputs forecasting for static or dynamic modes. For these aims, several approaches have been proposed for photovoltaic (PV) cell modeling including electrical circuit-based model, empirical models, and non-parametrical models. Moreover, numerous parameter extraction methods have been introduced in the literature depending on the proposed model and the operating mode. These methods can be classified into two main approaches including automatic numerical and analytical approaches. These approaches are commonly applied in the static mode, whereas they can be employed for dynamic parameters extraction. In this chapter, as a first stage, the static parameters extraction for both single and double diodes models is exposed wherein Genetic Algorithm and outdoor measurements are considered for fixed irradiation and temperature. In the second stage, a dynamic parameters extraction is carried out using Levenberg-Marquardt algorithm, where 1 day profile outdoor measurement is considered. After that, the robustness of the proposed approaches is evaluated and the parameters obtained by the static method and that given by the dynamic technique are compared. The test is carried out using 3 days with different weather conditions profiles. The obtained results show that the parameters extraction by dynamic techniques gives satisfactory performances in terms of agreement with the real data.

Keywords

  • photovoltaic module
  • static parameters extraction
  • dynamic parameters extraction
  • empirical model
  • electrical model

1. Introduction

The increasing development of PV technologies brought out their potential to provide the energy abundance across the world. Hence, they have been interested by several research groups in the purpose to improve their behavior and extend their life-time. Meanwhile, giving an accurate forecasting of the PV outputs behavior has been always a real issue related to their nonlinearity. Two modes can characterize the PV module in terms of modeling, the first one is the static mode wherein the obtained model is characterized and validated for fixed weather conditions (irradiation and temperature), while the second is the dynamic mode where the validation is carried out using variable weather conditions. In this context, several models of the PV cell/module have been introduced in the literature based mainly on several approaches including electrical, empirical, and non-parametrical modeling. For the non-parametrical models, two approaches are introduced namely, Artificial Neural Network (ANN) and Neuron-Fuzzy based models. The accuracy and the robustness of these approaches rely strongly on the richness of the training dataset in terms of scenarios. Neuron-Fuzzy techniques has been used to predict cell short-circuit current and open-circuit voltage for the static representation [1]. The ANN-based technique is introduced to model the PV array power for the embedded systems implementation [2]. This technique has been tested for dynamic mode dealing [3]. The ANN technique is also used for the prediction of the PV cell/module voltage directed to amorphous silicon PV technology wherein the obtained network has been tested by real dynamic data [4].

In the empirical or the analytical approaches, several models have been proposed to estimate accurately both static and dynamic modes of PV cell/modules. An analytical model is proposed by [5] based on manufacturer characteristic. This model provides acceptable results for both static and dynamic working. Another model has been introduced by Sandia National Laboratory [6], which is widely employed for PV cell/model forecasting especially for the large-scale arrays. Other empirical models have been proposed in order to estimate the PV array power under uniform shading [7, 8].

In the electrical approaches, two models widely prevail owing to their simplicity based on equivalent circuits namely: single (SDM) and double (DDM) diodes-based models. These models can deal with both static and dynamic modes with an acceptable accuracy. Besides, several enhancements have been introduced in these models to minimize parameters number and give more simplicity [9, 10]. The expressions of both photo-generation and diode saturation currents have been improved to give more accuracy in the dynamic working [11].

All models possess unknown parameters, which should be identified according to the module used in the practice. For this end, numerous approaches have been introduced in the literature including analytical and automatic numerical methods. The analytical methods usually rely on specific points on the I-V curve and on some value given by manufacturer. However, a significant error can be engendered if one of more of selected points is incorrect [12].

Owing to their adequate results, automatic numerical methods are prevail in PV models parameters extraction either through the use of deterministic algorithms such as: Newton model modified with Levenberg [13], Levenberg-Marquardt [14], Simulated Annealing algorithm (SA) [15], Pattern Search (PS) [16], Nelder-Mead Simplex algorithm (NMS) [17], and hybrid Nelder-Mead and modified particle swarm optimization [18] or by introducing metaheuristic algorithms such as: Genetic Algorithms (GA) [19], Particle Swarm Optimization (PSO) [20, 21, 22, 23, 24, 25, 26], Cuckoo Search (CS) [27], Artificial Bee Colony (ABC) [26, 28], and Artificial Bee Swarm (ABS) [29]. Moreover, other algorithms have been introduced like (FPA) [30, 31], hybrid Bee Pollinator Flower Pollination Algorithm (BPFPA) [31, 32], Harmony Search (HS) [33], Artificial Fish Swarm Algorithm (AFSA) [34], and other algorithms. The majority of the aforementioned algorithms have been applied for static parameters extraction.

Numerical algorithms have been also applied for dynamic parameters extraction wherein the identification process is carried out using variables weather conditions [35, 36].

In this chapter, modeling and parameters extraction of PV cell/module are detailed. Where, comparison study among three models by applying static and dynamic identification using out-door measurement.

Advertisement

2. PV cell/module modeling

The PV cell presents outputs variation, which depends on weather conditions namely, irradiation and temperature. As illustrated in Figure 1, for load variation from open circuit to short circuit, the PV cell shows nonlinear characteristic that possess a maximum point of power. For an optimal working, the load should be adapted at this point. In this section, three PV cell models will be employed and improved using automatic parameters extraction namely: the empirical Sandia model and both single and double diodes electrical models.

Figure 1.

PV characteristic for different irradiation and temperature. (a) Current versus voltage for different irradiation; (b) Power versus voltage for different irradiation; (c) Current versus voltage for different temperature; (d) Power versus voltage for different temperature.

2.1 Sandia model

This empirical model given by SANDIA National Laboratories provides relatively accurate dynamic forecast for PV cell/module by describing the thermal, the electrical, and the optical characteristics. Also, this model can be destined for any technology and can be adapted with any scale of PV arrays. Furthermore, its simplicity can qualify it to be used for real-time online prediction. Expressions (1)(4) describe the variation of Impp, Vmpp, and Pmpp, respectively.

Imp=Imp_STCC0Ee+C1Ee21+αmpTTSTCE1
Vmp=Vmp_STC+C2Ns×δT×lnEe+C3NsδTlnEe2+βmpTTSTCE2
δT=n×k×T+273.15qE3
Pmp=Imp×VmpE4

where, C0–3 are empirical parameters to be identified, Imp_STC, Vmp_STC are the current and the voltage in the maximum power point under standard test condition, Ee is the effective irradiation, K is the Boltzmann constant, q is the electron charge, δ(T) is the thermal voltage, αmp and βmp are, respectively, the current and the voltage temperature coefficient [36].

2.2 Single diode based model (SDM)

This physical model is based on the electrical approach illustrated in Figure 2 wherein the PV cell is composed of: a photo-generation current source and a diode while joule losses are represented by two serial and parallel resistors.

Figure 2.

Single diode equivalent circuit.

From this electrical representation, expression (5) and (6) can be obtained to describe the evolution of both current and voltage. The output current is expressed as a sum of the photo-generation current Iph, the diode current Id, and the shunt current Ish.

I=IphIdIshE5
Id=I0expqV+RsInKT1E6
Ish=V+RsIRshE7

where, K is the constant of Boltzmann, q is the electron charge, T is the cell temperature, n is the diode ideality factor, and I0 is the current saturation due to diffusion and recombination.

After the substitution of Eqs. (6) and (7) in (5), the following expression is obtained:

I=IphI0expqV+RsInKT1V+RsIRshE8

Rs, Rsh, and n are parameters to be identified in the static study and can be adjusted in the dynamic study.

Diode saturation current I0 is expressed in Eq. (9) function of the cell temperature and the energy band-gap [11].

I0=Isc×T3×expq×EgKTexpq×VocnKTr1×Tr3×expq.EgrKTrE9

Eq. (10) describes the evolution of the energy band-gap Eg as function of the cell temperature.

Eg=Eg0α×T2β+TE10

where Eg0 and Egr are the energy band-gap of the silicon at 0°C and at the reference temperature Tr, respectively, αand βare constants of the material.

The photo-generation current is given by Eq. (11) as a function of the irradiation and the cell temperature.

Iph=G1000Isc+µTTrE11

where G is the input irradiation, T is the cell temperature, Isc is the module short-circuit current, and μ is the coefficient temperature/short-circuit current (given by the manufacturer) [11].

Eg0, α, β, μ, Isc, and Voc are parameters to be identified in the dynamic study.

2.3 Double diode-based model

From the electrical representation illustrated in Figure 3, the PV cell can be represented by a source of current that represents the photo-generation, two diodes and both parallel and serial resistances representing the loss of energy inside the cell.

Figure 3.

Double diodes equivalent circuit.

After applying nodes law, the output current is expressed as sum of: photo-generation current Iph, shunt current Ish and the diodes currents Id1 and Id2 (Eqs. (12)(15)).

I=IphId1Id2IshE12
Id1=I01expqV+RsIn1KT1E13
Id2=I02expqV+RsIn2KT1E14
Ish=V+RsIRshE15

In which I01–2 are currents saturation of the two diodes that resulted from diffusion and recombination, n1–2 are ideally factors.

By substituting Eqs. (13)(15) in (12), final description of the output current versus the voltage is obtained which is expressed in Eq. (16) [11].

I=IphI01expqV+RsInK1T1I02expqV+RsIn2KT1V+RsIRshE16

n1–2 and Rs, Rsh are parameters which will be identified in the static study and they can be adjusted in the dynamic study.

Eqs. (17) and (18) express the evolution saturation currents of the diodes I01–2 versus energy band-gap Eg and cell temperature [11].

I01=Isc×T3×expq×Eg1KTexpq×Vocn1KTr1×Tr3×expq×EgrKTrE17
I02=Isc×T3×expq×Eg2KTexpq×Vocn2KTr1×Tr3×expq×EgrKTrE18

The photo-generation current is represented by the same expression of the single diode model (Eq. (11)).

Parameters Eg01–2, α1–2, β1–2, μ, Isc and Voc will be identified in the dynamic study.

Advertisement

3. Static parameters extraction of PV module

For fixed irradiation and temperature, a static parameters extraction will be done to extract five parameters in SDM and seven parameters in DDM. A numerical stochastic optimization algorithm is used in this identification. This algorithm namely, Genetic Algorithm (GA), is employed to minimize the cost function given in Eq. (19) which expresses the root mean square error (RMSE) between the measured PV module I(v) characteristic and that given by the models. For this and, outdoor static measurements have been carried out using the peak measuring device tracer (PVPM 2540C), whose characteristics are illustrated in Table 1. This device has been programed to provide both I(V) and P(V) curves of 101 samples per 1 min.

RMSE=1Ni=1NfVmIP2E19

where Iis the simulated current, Vmis the measured voltage, N is the number of sample in I(V) characteristics. The error between the measured and simulated I(V) characteristics for the aforementioned models are expressed in Eqs. (20) and (21).

fVmIP=IphI0expqVm+RsInKT1Vm+RsIRshImE20
fVmIP=IphI01expqVm+RsInK1T1I02expqVm+RsIn2KT1Vm+RsIRshImE21

3.1 Genetic Algorithm

The Genetic Algorithm (GA) is a stochastic algorithm imitated from the biological genetic process used to find an approximate solution for optimization problems. Like in the natural concept, the chromosome is the holder of the genes that the child can probably get from his parents. By analogy, these genes represent the variables (parameters) of the function to be minimized. Five steps can characterize the GA namely, generation of initial population, evaluation of fitness, selection, crossover and mutation [37, 38].

3.1.1 Initial population

The process starts by the generation the initial population of N chromosome coded in binary. Each vector chromosome is formed of group of parameters in which its length M is given in Eq. (22) wherein n is the number of parameters and Nb is the length of the sub-string (number of bits) of each parameter as shown in Figure 4. The length of the integer part given by the vector Conv (Eq. (23)) is used to limit the research domain in which, Pi (Eq. (24)) is the parameter value in decimal code [37, 38].

Figure 4.

Chromosome string.

M=n×NbE22
Conv=2ni2ni120212niNb+1E23
Pi=a0a1...aNb1×ConvTE24

3.1.2 Fitness

In this stage, the parameters values that have been randomly generated and decoded in decimal base will be substituted in the cost function to be optimized. The fitness is the solution of the parameters in the RMSE (x) function calculated in Eq. (19). Its value is mathematically expressed in Eq. (23) [37, 38].

Fitnessx=11+RMSExE25

3.1.3 Selection

The chromosomes that will participate as parents to generate a new child are chosen in this step. Any chromosome in the generated population can be chosen however, the individual that presents a good fitness have a high probability. The technique used for the chromosome choice is the roulette wheel illustrated in Figure 5, wherein the selection probability Ps, expressed in Eq. (24), is calculated, consists of a cumulative sum of the fitness of each chromosome orderly relative to the sum of all fitness. After that, the process generates a random drawing probability Pr. Hence, the first chromosome corresponds to Pr < Ps is chosen for the next steps (crossover and Mutation).

Figure 5.

Roulette wheel example.

Psi=k=1iFitnesskj=1NFitnessjE26

where, k is an integer counter that varies from 1 to the current chromosome, and j is an integer counter that varies from 1 to the population size N [37, 38].

3.1.4 Crossover

After selecting the chromosomes, the algorithm gives birth to new children by performing a crossover between each two chromosomes. For this end, a drawing probability Pr is generated and compared with the crossover probability Pc (usually high probability). Hence, the parents chromosomes that corresponds to Pr < Pc will be chosen for child generation. If not, the same chromosomes are kept. As illustrated in Figure 6, the crossover by point is used wherein the bits after the point randomly chosen are swapped [37, 38].

Figure 6.

Crossover process.

3.1.5 Mutation

In this step, the algorithm introduces a change in some characters of the selected chromosomes in order to expand the search space if the initial population does not fall in the optimal solution. In the binary coding, the selected bit will change from 1 to 0 and vice versa as described in Figure 7. The mutation has low probability Pm in which, it will be affected for characters that correspond to Pr < Pm in which, Pr is the drawing probability (randomly generated) [37, 38].

Figure 7.

Mutation process.

Our system is formed of mono-crystalline PV module SANYO technology with the characteristics listed in Table 2, peak measuring tracer and the necessary sensors. The experimental platform is illustrated in both Figure 8a and b.

ApplicationDC voltageDC currentTemperatureIrradianceMeasuring period single measurementI-V curve samples
PV modules and small strings25/50/100/250 V2/5/10/40 A−40°C to +120°C with Pt10000–1300 (W/m2) (standard-sensor)0.02–2 (s)101

Table 1.

PVPM2540C characteristics.

DescriptionSANYO mono-crystalline
Cell number96
Cell typeMono-crystalline
Cell size156 × 156 mm
PV module dimension1319 × 894 × 35 mm
Nominal power180 W
Open circuit voltage Voc66.4 V
Short circuit current Isc3.65 A
Voltage Vmpp54 V
Current Impp3.33 A
Nominal operating temperature NOCT45 ± 2°C
Temperature coefficient (Pmax)−0.33%/°C
Temperature coefficient (Isc)1.10 mA/°C
Temperature coefficient (Voc)−0.173 V/°C

Table 2.

PV modules characteristics.

Figure 8.

(a) Pictures for I-V experimental platform and (b) synoptic scheme for I-V measurement.

After running of GA for 1000 cycles with the parameters listed in Table 3 for both SDM and DDM using outdoor measurement of the systems, wherein extracted parameters are listed in Table 4. Figure 9a and b illustrate the agreement between the measured and simulated I(V) and P(V) characteristics for SDM model whose obtained parameters are summarized in Table 5. Figure 10a and b show the agreement between the measured and simulated I(V) and P(V) characteristics for DDM model whose extracted parameters are summarized in Table 6.

GA parametersValue
Number of cycle1000
Population length500
Crossover probability0.7
Mutation probability0.2

Table 3.

The GA parameters.

The electrical parameterIph [A]I0 [A]nRs [Ω]Rsh [KΩ]
The identified value3.019549591e−0051.874*960.32738.1514

Table 4.

The obtained parameters for SDM.

Figure 9.

Matching between measured and simulated characteristics for SDM. (a) I(V) and (b) P(V).

The electrical parameterIph [A]I01 [A]I02 [A]n1n2Rs [Ω]Rsh [KΩ]
The identified value3.02896.1035e−0053.8147–0061.3658*961.9179*960.10175.992

Table 5.

The obtained parameters for DDM.

Figure 10.

Matching between measured and simulated characteristics for DDM. (a) I(V) and (b) P(V).

EΔE
NSNLZEPLPS
NSPSPSPLPSPS
NLPSPLPLPLPS
ZENLNLZEPLPL
PLNSNLNLNLNS
PSNSNSNLNSNS

Table 6.

The inference table.

Some parameters will be identified again in the dynamic study including the parameters involved in I0 and Iph equations, while the remaining will be adjusted to give more accuracy under variable weather conditions.

Advertisement

4. Dynamic parameters extraction of PV module

In this section, dynamic parameters identification will be described wherein the process is done by using 1 day profile of measurement. This allows to improve the nominal values given by the manufacturer, which can cause a significant error due to operating conditions and the consumed lifetime. Moreover, parameters obtained by static method can be adjusted by dynamic identification. For this end, automatic parameters adjustment using Levenberg-Marquardt optimization algorithm is employed.

As illustrated in Figure 11, the main idea is to take both PV module model and the MPPT as a single system with three outputs namely, Impp, Vmpp, and Pmpp. These outputs will be compared with 1 day profile of outdoor measurements. The process consists in minimizing the error between the model outputs and the real data. The whole system has been implemented in Matlab/Simulink tool.

Figure 11.

Dynamic parameters extraction process.

For Sandia model, the process is carried out without the use of MPPT considering that this model has been established to the dynamic forecasting.

4.1 The MPPT used

An Accurate fuzzy logic MPPT algorithm is employed in our system (for SDM and DDM) in order to get satisfactory results in terms of precision and accuracy. The algorithm is used to control a DC/DC boost converter for the purpose to keep the PV module working at the maximum point of power. Mamdani inference model is used with two inputs namely, the error Eand the variation of the error ΔE. The calculation of these attributes is expressed in Eqs. (27)(30).

ΔPpvn=PpvnPpvn1E27
ΔVpvn=VpvnVpvn1E28
En=ΔPpvnΔVpvnE29
ΔEn=EnEn1E30

Three steps can characterize the fuzzy algorithm; the first one is the fuzzification process that consists on the conversion of the numerical inputs values (E and ΔE) into linguistic values by the substitution in the membership functions. The second step is the inference process, which is considered as the main stage in the fuzzy algorithm wherein the relation between the inputs and the output is done. The third step is the defuzzification where the process converts the linguistic decision into numerical output. Figure 12 describes briefly the fuzzy processing steps [39].

Figure 12.

Fuzzy logic algorithm steps.

For both inputs and output, five trapezoidal and triangular membership functions have been employed namely: NS (negative strong), NL (negative low), ZE (zero), PL (positive low) and PS (positive strong). The center of gravity based method is used for the defuzzification to provide the control duty cycle after applying the Mamdani inference model given in Table 6. Figure 13a,b and d describes the used membership functions [39].

Figure 13.

(a) Error membership; (b) variation of error membership; and (c) duty cycle membership.

4.2 Simulation study

The Levenberg-Marquardt algorithm is implemented using 1 day profile of outdoor real measurement of dynamic PV outputs (Pmpp, Impp and Vmpp). The process consists in minimizing the error between simulated outputs of both SDM and DDM and 8 h of real data (09:00 am–05:00 pm). The peak measuring device tracer (PVPM 2540C) has been programmed to provide 1 sample per minute. Table 7 lists the lower and upper limits search of the extracted parameters. The extracted parameters using the dynamic method are summarized in Table 8. The inputs measurement of the irradiation and the temperature are illustrated in Figure 14a and b, respectively. Satisfactory results have been obtained in terms of matching between the real data and the simulated outputs Pmpp, Impp and Vmpp for SDM, DDM and Sandia as shown in Figure 15ac, respectively.

SDM parametersBoundariesDDM parametersBoundariesSandia parametersBoundaries
n[0, 2]*96n1, n2[0, 2]*96C0[0, 2]
Rs [Ω][0, 1]Rs [Ω][0, 1]C1[−1, 1]
Rsh [Ω][0, 104]Rsh [Ω][0, 104]C2[−10, 10]
Eg0[0, 1]Eg012[0, 2]C3[−10, 50]
A[0, 1]α12[0, 1]αImp[°C−1][0, 1]
B[0, 104]β12[0, 104]βVmp[V/°C][−1, 0]
μ[0, 1]μ[0, 1]n
Isc[3, 3.7]Isc[3, 3.7]
Voc[60, 66.8]Voc[60, 66.8]

Table 7.

Lower and upper boundaries selected for each model parameters.

SDM parametersValuesDDM parametersValuesSandia parametersValues
n105.73/96n1, n290.73/96; 73.39/96C01.058
Rs [Ω]0.82495Rs [Ω]0.3219C10.020
Rsh [Ω]8.371 × 103Rsh [Ω]4.9664 × 103C2−0.341
Eg0 [ev]1.4525Eg01–2 [ev]1.649; 1.31C3−9.997
α6.56 × 10−4α1–20.0018; 0.0132αImp[°C−1]2.53 × 10−14
β126.11β1–2694.84; 1020.76βVmp[V/°C]−0.203
μ0.0121μ0.0112n1.221
Isc3.671Isc3.629
Voc66.208Voc65.527

Table 8.

Dynamic extracted parameters.

Figure 14.

Eight hours profile (09:00 am–05:00 pm). (a) Irradiation and (b) temperature.

Figure 15.

Eight hours profile (09:00 am–05:00 pm) of: measured output VS obtained output for SANYO mono-crystalline. (a) Measured Pmpp vs obtained Pmpp from: SDM, DDM and Sandia, (b) measured Impp vs obtained Impp from: SDM, DDM and Sandia and (c) measured Vmpp VS obtained Vmpp from: SDM, DDM and Sandia.

Advertisement

5. Experimental validation

In this section, a validation with an unseen data is carried out to test and compare the effectiveness of the proposed enhancement. The three developed models will be compared with real measurement profile (09:00 am–05:00 pm) of irradiation and temperature for different weather conditions. Wherein, the SDM and DDM models using the developed parameters are compared with the former nominal parameters listed in Table 9, Sandia model and the real data. Three different skies of real measurement have been used for this validation namely, clear day, semi-cloudy day and cloudy day. The matching in the power (Pmpp) between and the real data and SDM and DDM with nominal parameters, SDM and DDM with the new parameters and Sandia model is illustrated in Figure 16ac for clear day, semi-cloudy day and cloudy day respectively. Besides, the agreement in the voltage and the current (Vmpp and Impp) for these models with the real data is shown in Figure 17ac and Figure 18ac, respectively.

ParametersValue
Eg01.16 eV
α4.73 × 10−4
β1000
μ1.10 mA/°C
Isc3.65 A
Voc66.4 V

Table 9.

Nominal parameters.

Figure 16.

Eight hours profile (09:00 am–05:00 pm) of: measured Pmpp vs obtained Pmpp from: SDM, SDM with nominal parameters, DDM, DDM with nominal parameters and Sandia. (a) clear day, (b) semi-cloudy day and (c) cloudy day.

Figure 17.

Eight hours profile (09:00 am–05:00 pm) of: measured Impp vs obtained Impp from: SDM, SDM with nominal parameters, DDM, DDM with nominal parameters and Sandia. (a) clear day, (b) semi-cloudy day and (c) cloudy day.

Figure 18.

Eight hours profile (09:00 am–05:00 pm) of: measured Vmpp vs obtained Vmpp from: SDM, SDM with nominal parameters, DDM, DDM with nominal parameters and Sandia. (a) clear day, (b) semi-cloudy day and (c) cloudy day.

It is clearly found that an improved agreement has been shown by models with new parameters compared to that given by the nominal parameters and the static method.

For more clarity, the hourly power efficiency given by the presented models and the real data has been calculated. It consists on the average of the power during 1 h versus the optimal PV module power (Eq. (31)) [40].

Hourly efficiency=k=1NPMPPk/NPPV_optimalE31

in which, N is the number samples per hour.

Figure 19ac show the bar-graph of the hourly power efficiency of the proposed models for the three weather conditions, namely the clear day, semi-cloudy day and cloudy day, respectively. The enhanced models present higher hourly power efficiency versus models with the former parameters and those given by the static technique. Furthermore, root mean square error (RMS) and the mean absolute error (MAE) between the real data and the studied models are calculated by Eqs. (19) and (32) to show the enhancement of the proposed method.

Figure 19.

Eight hours profile (09:00 am–05:00 pm) of: measured hourly efficiency vs obtained hourly efficiency from: SDM, SDM with nominal parameters, DDM, DDM with nominal parameters and Sandia. (a) Clear day, (b) semi-cloudy day, and (c) cloudy day.

MAE=1Ni=1NSimulatedoutputMeasuredoutputMeasuredoutputE32

where N is the number of samples [41].

The aforementioned results show clearly that the extracted parameters of the PV module using dynamic techniques present more accuracy compared with the static method and the parameters given by the manufacturer. Indeed, the parameters obtained by the static method are clearly improved for variable weather conditions (irradiation and temperature), which is confirmed using different skies. Table 10 summarizes the calculated RMS and MAE errors values which show obviously that the developed models present advantages comparing with real outdoor data of different weather conditions.

DayWeatherError [%]CurrentVoltage
SDMSDMnpDDMDDMnpSandiaSDMSDMnpDDMDDMnpSandia
22/01/2018Clear dayRMS0.721.120.690.930.681.482.211.632.341.49
MAE4.315.054.245.124.381.274.071.773.701.17
29/01/2018Semi-cloudyRMS0.220.450.240.320.233.309.106.537.871.13
MAE3.145.402.964.352.851.505.342.374.760.57
28/01/2018Cloudy dayRMS0.430.480.470.310.529.058.7712.212.31.83
MAE4.386.625.023.935.602.064.634.056.180.77
DayWeatherError [%]Power
SDMSDMnpDDMDDMnpSandia
22/01/2018Clear dayRMS3.658.493.454.343.79
MAE4.246.793.744.444.51
29/01/2018Semi-cloudyRMS2.484.302.484.300.82
MAE5.1411.74.5810.32.33
28/01/2018Cloudy dayRMS5.354.825.515.242.79
MAE1.371.561.551.575.63

Table 10.

Calculated RMS (%) and MAE (%).

Advertisement

6. Conclusion

In this chapter, both dynamic and static parameters identification methods have been highlighted and compared with real measurement. The SDM and DDM nominal parameters involved in I0 and Iph equations have been developed by dynamic method. This improved result has been compared with that given by the static technique and Sandia model versus out-door real data for different skies (clear day, semi-cloudy day and cloudy day). It was found that SDM and DDM based on the parameters extracted by dynamic method give satisfactory accuracy, which is confirmed by some calculated indicator such as: the hourly efficiency and both root mean square error (RMS) and the mean absolute error (MAE). This allows to solve modeling problems of PV module that apply for several applications such as fault detection.

References

  1. 1. AbdulHadi M, Al-Ibrahim AM, Virk GS. Neuro-fuzzy-based solar cell model. IEEE Transactions on Energy Conversion. September 2004;19(3):619-624
  2. 2. Mekki H, Mellit A, Salhi H, Khaled B. Modeling and simulation of photovoltaic panel based on artificial neural networks and VHDL-language. In: Electronics, Circuits and Systems, 2007. ICECS 2007. 14th IEEE International Conference.
  3. 3. Mellit A, Saglam S, Kalogirou SA. Artificial neural network-based model for estimating the produced power of a photovoltaic module. Renewable Energy. 2013;60:71-78
  4. 4. Askarzadeh A. Voltage prediction of a photovoltaic module using artificial neural networks. International Transactions on Electrical Energy Systems. 2014;24:1715-1725
  5. 5. Ortiz-Rivera EI, Peng FZ. Analytical Model for a Photovoltaic Module Using the Electrical Characteristics provided by the Manufacturer Data Sheet Power Electronics Specialists Conference. PESC '05. IEEE 36th. 2005
  6. 6. King DL, Boyson WE, Kratochvil JA. Photovoltaic array performance model. Sandia Report No. SAND 2004-3535. 2004. Available from:http://www.sandia.gov/pv/docs/PDF/King%20SAND.pdf
  7. 7. Zhou W, Yang H, Fang Z. A novel model for photovoltaic array performance prediction. Applied Energy. 2007;84:1187-1198
  8. 8. Deline C, Dobos A, Janzou S, Meydbray J, Donovan M. A simplified model of uniform shading in large photovoltaic arrays. Solar Energy. 2013;96:274-282
  9. 9. Chenni R, Makhlouf M, Kerbache T, Bouzid A. A detailed modeling method for photovoltaic cells. Energy. 2007;32:1724-1730
  10. 10. Rosell JI, Ibanez M. Modelling power output in photovoltaic modules for outdoor operating conditions. Energy Conversion and Management. 2006;47:2424-2430
  11. 11. Castaner L, Silvestre S. Front matter. In: Castaner RL, Silvestre S, editors. Modeling Photovoltaic Systems Using PSpice. Chichester, England: Wiley; 2003. p. 376
  12. 12. Ishaque K, Salam Z, Taheri H, Shamsudin A. A critical evaluation of EA computational methods for photovoltaic cell parameter extraction based on two diode model. Solar Energy. 2011;85:1768-1779
  13. 13. Easwarakhanthan T, Bottin J, Bouhouch I, Boutrit C. Nonlinear minimization algorithm for determining the solar cell parameters with microcomputers. International Journal of Solar Energy. 1986;4:1-12
  14. 14. Dkhichi F, Oukarfi B, Fakkar A, Belbounaguia N. Parameter identification of solar cell model using LevenbergeMarquardt algorithm combined with simulated annealing. Solar Energy. 2014;110:781-788
  15. 15. El-Naggar KM, AlRashidi MR, AlHajri MF, Al-Othman AK. Simulated annealing algorithm for photovoltaic parameters identification. Solar Energy. 2012;86:266-274
  16. 16. AlHajri MF, El-Naggar KM, AlRashidi MR, Al-Othman AK. Optimal extraction of solar cell parameters using pattern search. Renewable Energy. 2012;44:238-245
  17. 17. Chen Z, Wu L, Lin P, Wu Y, Cheng S. Parameters identification of photovoltaic models using hybrid adaptive Nelder-Mead simplex algorithm based on eagle strategy. Applied Energy. 2016;182:47-57
  18. 18. Hamid NFA, Rahim NA, Selvaraj J. Solar cell parameters identification using hybrid Nelder-Mead and modified particle swarm optimization. Journal of Renewable and Sustainable Energy. 2016;8. DOI: 10.1063/1.4941791
  19. 19. Jervase JA, Bourdoucen H, Al-Lawati A. Solar cell parameter extraction using genetic algorithms. Measurement Science and Technology. 2001;12(11):1922
  20. 20. Huang W, Jiang C, Xue L, Song D. Extracting solar cell model parameters based on chaos particle swarm algorithm. In: 2011 International Conference on Electric Information and Control Engineering. Vol. 2011. pp. 398-402
  21. 21. Meiying Y, Xiaodong W, Yousheng X. Parameter extraction of solar cells using particle swarm optimization. Journal of Applied Physics. 2009;105:094502-094508
  22. 22. Hachana O, Hemsas KE, Tina GM, Ventura C. Comparison of different metaheuristic algorithms for parameter identification of photovoltaic cell/module. Journal of Renewable and Sustainable Energy. 2013;5:053122
  23. 23. Ma J, Man KL, Guan SU, Ting TO, Wong PWH. Parameter estimation of photovoltaic model via parallel particle swarm optimization algorithm. International Journal of Energy Research. 2016;40:343e352
  24. 24. Jordehi AR. Time varying acceleration coefficients particle swarm optimization (TVACPSO): A new optimization algorithm for estimating parameters of PV cells and modules. Energy Conversion and Management. 2016;129:262-274
  25. 25. Muralidharan R. Parameter extraction of solar photovoltaic cells and modules using current-voltage characteristics. International Journal of Ambient Energy (2016). Europe-China cooperation and competition in the Green Industries; 2017;38(5):509-513. DOI: 10.1080/01430750.2016.1144525
  26. 26. Oliva D, Cuevas E, Pajares G. Parameter identification of solar cells using artificial bee colony optimization. Energy. 2014;72:93-102
  27. 27. Ma J, Ting TO, Man KL, Zhang N, Guan S-U, Wong PWH. Parameter estimation of photovoltaic models via Cuckoo search, Hindawi Publish. Corporat. Journal of Applied Mathematics. 2013;362619. 8 pages
  28. 28. Jamadi M, Merrikh-Bayat F, Bigdeli M. Very accurate parameter estimation of single- and double-diode solar cell models using a modified artificial bee colony algorithm. International Journal of Energy and Environmental Engineering. 2016;7:13-25
  29. 29. Askarzadeh A, Rezazadeh A. Artificial bee swarm optimization algorithm for parameters identification of solar cell models. Applied Energy. 2013;102:943-949
  30. 30. Alam DF, Yousri DA, Eteiba MB. Flower pollination algorithm based solar PV parameter estimation. Energy Conversion and Management. 2015;101:410-422
  31. 31. Benkercha R, Moulahoum S, Colak I, Taghezouit B. PV module parameters extraction with maximum power point estimation based on flower pollination algorithm. In: 2016 IEEE International Power Electronics and Motion Control Conference, PEMC. 2016. pp. 442-449
  32. 32. Ram JP, Babu TS, Dragicevic T, Rajasekar N. A new hybrid bee pollinator flower pollination algorithm for solar PV parameter estimation. Energy Conversion and Management. 2017;135:463-476
  33. 33. Askarzadeh A, Rezazadeh A. Parameter identification for solar cell models using harmony search-based algorithms. Solar Energy. 2012;86:3241-3249
  34. 34. Han W, Wang H-H, Chen L. Parameters identification for photovoltaic module based on an improved artificial fish swarm algorithm. Scientific World Journal. 2014;859239. 12 pages
  35. 35. Kichou S, Silvestre S, Guglielminotti L, et al. Comparison of two PV array models for the simulation of PV systems using five different algorithms for the parameters identification. Renewable Energy. 2016;99:270-279
  36. 36. Blaifi S-a, Moulahoum S, Taghezouit B, Saim A. An enhanced dynamic modeling of PV module using Levenberg-Marquardt algorithm. Renewable Energy. 2019;135:745-760
  37. 37. Blaifi S, Moulahoum S, Colak I, Merrouche W. An enhanced dynamic modelof battery using genetic algorithm suitable for photovoltaic applications. Applied Energy. 2016;169:888-898
  38. 38. Blaifi S, Moulahoum S, Colak I, Merrouche W. Monitoring and enhanced dynamic modeling of battery by genetic algorithm using LabVIEW applied in photovoltaic system. Electrical Engineering. June 2018;100(2):1021-1038
  39. 39. Blaifi S, Moulahoum S, Benkercha R, Taghezouit B, Saim A. M5P model tree based fast fuzzy maximum power point tracker. Solar Energy. 2018;163:405-442
  40. 40. Benyoucef AS et al. Artificial bee colony based algorithm for maximum power point tracking (MPPT) for PV systems operating under partial shaded conditions. Applied Soft Computing. 2015;32:38-48
  41. 41. Blaifi S, Moulahoum S, Colak I, Merrouche W. An enhanced dynamic model of battery using genetic algorithm suitable for photovoltaic applications. Applied Energy. 2016;169:888-898

Written By

Sid-Ali Blaifi and Bilal Taghezouit

Submitted: April 18th, 2019 Reviewed: August 30th, 2019 Published: April 1st, 2020