Open access peer-reviewed chapter

Static and Dynamic Photovoltaic Cell/Module Parameters Identification

By Sid-Ali Blaifi and Bilal Taghezouit

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

DOI: 10.5772/intechopen.89449

Downloaded: 57

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.

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.

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 I is the simulated current, Vm is 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.

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 E and 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.

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 (%).

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.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sid-Ali Blaifi and Bilal Taghezouit (April 1st 2020). Static and Dynamic Photovoltaic Cell/Module Parameters Identification, Advanced Statistical Modeling, Forecasting, and Fault Detection in Renewable Energy Systems, Fouzi Harrou and Ying Sun, IntechOpen, DOI: 10.5772/intechopen.89449. Available from:

chapter statistics

57total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Advanced Statistical Modeling, Forecasting, and Fault Detection in Renewable Energy Systems

Edited by Fouzi Harrou

Next chapter

Performance Analysis of a 20 MW Grid-Connected Photovoltaic Installation in Adrar, South of Algeria

By Nouar Aoun

Related Book

First chapter

Potential of the Solar Energy on Mars

By Dragos Ronald Rugescu and Radu Dan Rugescu

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us