1. Introduction
According to [
1], a fault can be defined as “an abnormal condition that may cause a reduction in, or loss of, the capability of a functional unit to perform a required function”. Faults in photovoltaic (PV) installations can significantly affect their energy yield, which is why the development of fault detection and diagnosis techniques has become an essential topic of research in recent years, Refs. [
2,
3,
4]. Research works that analyze failure mode in solar modules, e.g., Refs. [
5,
6,
7,
8,
9,
10,
11], have identified that the power reduction in PV modules under fault causes a deformation of the current-voltage curve, e.g., curves in [
12].
Recently, Ref. [
13] presented a performance analysis with 30 faulty modules of different brands. The study shows that, on average, the power loss after two years, with several faults, drops around
; also, the voltage variation at maximum power was, on average
. Similar research in India concluded that, on average, power degradation in PV modules is found to increase 1.4% per year over 25 years [
14]. But this degradation rate is not a static value; it depends on the life stage of the PV installations as is shown in [
15], the first 18 years around 0.1% per year; then the next period of 10 years, the rate raises around 1%/year, and the final life stage, the rate increases more than 1.2%/year.
It is important to highlight that failure modes at the module level affect the whole PV installation performance because the PV array is a symmetrical composition of PV modules as is demonstrated by Gokmen in [
16]. Eighteenth types of faults in solar modules have been described in [
17], their detectable effects sometimes overlap, so it is not easy to distinguish among them. Furthermore, the diagnoses of faults require combining knowledge of different domains (visual inspection, thermography, electrical, chemical, material analyses, so on), Refs. [
18,
19]. This means that the specific failure modes could have, for example, the same behavior in the electric domain.
The detection criterion used in this paper is based on the widely accepted five parameters model, which, according to [
20], can represent a single solar cell or several cells connected in series (PV module). This model is mathematically represented as follows,
where
I and
V represent the current and voltage of the module,
is the photo-current generated,
is the saturation current of the diode,
and
are the parasitic series and parallel resistances,
is the number of cells in series,
represent the ideality factor, and
is the thermal voltage.
An increase in the PV module series resistance can model several fault mechanisms, as it has been indicated in [
21,
22]. However, using other parameters to explain other failure modes is scarce in the literature. In this regard, one of the contributions of this work is to propose parameter ranges in the PV model that represent others failure modes not included in the typical models, e.g., only series resistance variation. This work will focus on the most common types of degrading faults listed next.
Series resistance degradation (SRD) faults. The increase of the parasitic resistance
could be expected by thermal cycling, solder corrosion, broken ribbons, and homogeneous soldering disconnection, so on. For example, Ref. [
23] indicates that normally the increment is
year. However, the variation of
could be up to ten times its original value, according to [
24].
Optical homogeneous degradation (OD) faults. The photocurrent
is affected by homogeneous glass corrosion, homogeneous contamination, loss of transparency, homogeneous corrosion of the anti-reflection coating. For instance, the yellowing and the browning are faults that cause power losses up to
[
25]. According to [
23] this parameter decreases its value between
year, for simulation purpose a window of 20 years were simulated.
Optical heterogeneous degradation (OHD) faults. This fault family is similar to the one presented before, but in this case, the level of affectation is not homogeneous in two ways, location and degree of degradation. According to [
17] the parasitic parallel resistance
decreases, as well as
. The variation of
was taken as equivalent to OD faults, and the variation of
was taken as the worse case reported in PID faults.
Potential induced degradation (PID) faults. This fault occurs when high potential voltages of the system cause a leakage current between the frame and the solar cells. According to [
26], the shunt resistance
decreases as well as the short-circuit current
, but for this work, we assume that the changes in
are caused by
because these variables are proportional
. In [
27] has proposed that both parameter decrease proportionally,
. Hence, the range of variation for these parameters are
and
. We call this variant as
.
Also, another PID fault variation, we call
, is found in [
28], and they show that the parameters
moves between
,
varies between
and
moves between
. These percentages were calculated by analyzing the maximum and minimum values reported in the graphical results.
Micro cracks (MC) faults. Micro crack can reduce energy production up to
of a module with only
of the affected area in one cell [
29]. For these faults, according to [
30],
and
decrease proportionally because both parameters are inversely proportional to the affected area; hence the range of variation is
.
Light induced degradation (LID) faults. According to [
31] LID is defined as the increment of the recombination current in the base in P-type silicon wafers. The effect is also observed in the open-circuit voltage [
17]. For a particular module, the saturation current
could be moved in a range of
times; this variation can produce changes of open-circuit voltage
in the range of
[
20]. These percentages were calculated by analyzing graphical results.
A second contribution is the proposal of new detection criteria based on the measurements of the module voltages. Due to its simplicity, this criterion is suitable for implementation in monitoring strategies such as the one presented in [
32]. In fact, several approaches use voltage measurements to detect and diagnose faults in the array, for example [
33,
34,
35,
36,
37,
38], but they require other variables like the string current, time, or even the ambient variables, to generate alarms adequately. Also, in those papers, the nature of faults analyzed is different because they focus mainly on faults such as diode short-circuit, open-circuits, partial shadows, and degradation faults represented only as an increase of the serial resistance. Furthermore, a similar indicator based on voltage is presented in [
39], which is called
, but this indicator requires knowing precisely the output voltage of the array and the theoretical voltage at its maximum power. This indicator is quite different from our approach since our criterion is calculated from the voltages of the modules in the PV string.
This work presents and analyzes new detection criteria to demonstrate that is capable to detect all the progressive faults only with the module voltage variable in multi-crystal PV modules. This numerical simulation study is divided into the following main sections:
Section 2 describes the theoretical framework used to derive the detection criterion.
Section 3 links specific fault conditions with changes in the variation of some parameters of the PV array. This section also explains the simulated experiment, circuits and software used, and the calibration process.
Section 4 and
Section 5 show the main results and the discussion. Finally, the main conclusions are highlighted in
Section 6.
2. Theoretical Framework: Detection Criterion Fundamentals
The explicit expression of (
1) for the module voltage is obtained using Lambert-
transformation, if we have an equation with the form
, the term
x could be obtained as
for any
. Therefore, the voltage is expressed as follows,
It is clear that if a failure affects the PV module performance, the parameters in (
1) and (
2) change their values and hence the voltage
V in (
2) also varies [
28]. If the initial values of all the parameters are known it is possible to identify the degradation condition of one PV module using (
2). Nevertheless, the aforementioned is seldom the case. Not only the initial values are unknown but also a PV installation is comprised of several PV modules connected in series and parallel and the relationships between the parameter variation in one PV module will affect the others.
In a typical photovoltaic plant, PV modules are connected using a configuration known as series-parallel (SP). In an SP configuration, the PV modules are first connected in series forming strings which are then connected in parallel as depicted in
Figure 1, forming a PV array of
modules. The voltages in the
j-string, formed by
m modules, are governed by Kirchhoff’s voltage law [
20], as follows:
where
is the differential voltage of the module
i in the string
j as in Equation (
2),
is the voltage of the blocking diode, and
is the operational voltage point of the array.
When all the modules are identical (
3) can be simplified as follows,
Now consider the case in which at least one module is degraded or faulty. Due to the SP configuration the sum of all the voltages in the string with the faulty PV module does not change, see Equation (
3). Therefore, the PV modules that are part of the string in which a faulty module is located will change their electrical operating point. We call the aforementioned elements
affected PV modules, as it can be seen in
Figure 2.
If we denote
as the number of modules with the same fault, (
3) can be rewritten as,
where
is the voltage of the faulty module (
degradeted) and
is the voltage of the affected modules.
Notice that the sum of the voltages at the faulty and at the affected modules in a given string should be equal to the sum of the voltages of non-faulty modules in a string parallel to it. Hence, the variation in the voltages of the affected modules and the faulty modules must behave oppositely, i.e., when the voltage at the faulty decreases the voltage at the affected modules increases. Moreover, the difference between the voltages at the affected and the faulty modules can be expressed as follows,
Given that the current trough the string,
I, is shared by the modules, (
6) can be rewritten as
where
and
are the power produced by the affected and faulty modules, and
I is the current of the
j-string (
j). Therefore (
6) represents a measure proportional to the power reduction in a PV module.
4. Results
The results of the SRD fault are presented in
Figure 3. The SRD fault was simulated by varying the resistance R4 up to ten times its original value.
Figure 3a shows how the current of the faulty string changes as the degradation process caused by resistor R4 continues. This means that the differential voltages
and
are separating as the resistor R4 is increasing; see
Figure 3b. The relative percentage difference calculated with Equation (
6) is presented in
Figure 3c. In this case, it is shown that there is a power loss by the faulty module concerning the non-faulty modules in the same string; the arrow shows the degradation process. The array’s output power is shown in
Figure 3d, and it goes down according to the increment of the serial resistance R4; also, the degradation process is indicated with the arrow.
Table 5 shows how the resistor value affects the voltage at the maximum power-point (
), the maximum power point (
), and the total power loss (TPL). The TPL indicates how much power is lost in the whole PV array with respect to an expected non-faulty value of 1600 W, as indicated in Equation (
8). Additionally, the table shows the value of voltage for the affected (
) and faulty (
) PV modules as well as the relative percentage difference (
) between them. The last column in the table shows the module power loss (MPL) by the faulty module with respect to its maximum power as reported in the datasheet (200 W).
The graphical results for the following six fault types are presented in the
Appendix A, each figure belongs to one fault type, and it contains the same four charts explained before for the SRD fault. In the
Appendix A, the
Figure A1 presents OD fault, and the results were obtained moving I4 values randomly.
Figure A2 presents OHD fault; this simulation changes the values for R8 and I4 randomly and separately.
Figure A3 presents PID
fault; here, the elements R8 and I4 move randomly but proportionally. Next,
Figure A4 presents PID
fault; for this simulation, R8, I4, and the parameter
in diode D4 move freely and randomly. Further, the results for MC fault are presented in
Figure A5, here I4, and the parameter
in the diode D4 move randomly but in a proportional way. Finally,
Figure A6 shows the results for LID fault when the parameter
moves in the diode D4. It is interesting to see that all fault types present similar patterns and behaviors in all the charts.
5. Discussion
Only the PV module under examination has one type of progressive fault in the previous circuit simulations. The rest of the elements, such as PV modules, wires, bypass diodes, were assumed to work in a non-faulty condition. The aforementioned does not mean that the results are specific to the selected module location; if the faulty module is located at any other place in the array, the obtained charts would be equivalent to the graphs presented in this work.
Based on the simulation results it is confirmed that when a progressive fault affects a photovoltaic module, its differential voltage changes according to the severity of the failure mode.
Figure 3b confirms that the progressive fault unbalances the voltage in the string; this means that as the progressive fault progress, the faulty module’s differential voltage decrease, and the differential voltage of the affected modules increase. Hence, the relative percentage difference (
) increases, and the output power of the array decreases; these facts are easily checked in
Figure 3c,d. For all the fault types analyzed, the same performance was observed in sub-charts in
Figure A1,
Figure A2,
Figure A3,
Figure A4,
Figure A5 and
Figure A6.
The degree of fault affectation in the module is calculated with the relative percentage difference (
6), and this indicator is equivalent to the percentage of power loss if the modules work at the same string current.
Table 5 presents the numerical results for the SRD fault. These results correspond to the evolution of the maximum power-point. In these maximum power points, it seems that the power loss indicator (
) is highly correlated with the power loss of the module calculated at standard test conditions, MPL. Actually, for this fault condition, the estimated power
and the MPL are highly correlated, the Pearson Coefficient [
43] is, in this case,
; therefore, it is possible to calculate a linear regression as,
Moreover, the obtained linear regression has a very low variability, given a determination coefficient
of
. This result is not unique for this failure mode, the correlation between
and MPL appears for all the fault types analyzed. Graphically, it can be appreciated in
Figure 4 and the parameters for the best fit line equation presented in
Table 6.
In this PV array model of 1.6 kWp, as the progressive fault increase, the
moves affecting all the modules’ power production; this means generating small losses in all the modules in the array. These small changes in power production caused by the progressive fault could be hidden because they are too small to be detected by the monitoring system. For the 1.6 kWp example, if
double its value in whatever PV module, the total power loss of the array drops about 0.86%, which could be negligible. However, this percentage represents one faulty module losing around 6.26% of its energy (See the third row in
Table 5). This hidden effect is even more drastic in large PV installations; let’s suppose that a faulty module has a constant power loss, and if the number of PV modules in the array is high enough, the total power loss tends to zero. However, the faulty module is still losing the same percentage of energy and may evolve into a more severe fault condition.
The correlation
can be used to detect abnormal behavior. Moreover, it is possible to develop a simple criterion for fault detection and location in SP arrays. The Equation (
10) is based on (
6), and it checks if the relative percentage difference between two voltages is more than a threshold value called
. Here
represents the maximum differential voltage of all the modules in the string
j, and
is the differential voltage of the analyzed module. It should be noticed that
is the higher voltage in the string, and it belongs to the module with less degradation. This criterion is simple to incorporate into other fault detection proposals that use differential voltages as input signals such as [
44,
45].
To define the fault threshold, it is crucial to know the environment in which the solar array is placed. For instance, soiling may not be a problem in tropical regions benefiting from rainfall cleaning. On the contrary deserts or dry places, the PV array could be affected by dust storms or air pollution that could reduce the general performance [
46,
47]. In [
48] it has been reported that power losses in outdoor conditions could be reduced between
and for laboratory conditions, it is possible to reduce it up to
. For instance, in tropical weather like Phitsanulok Tayland, it has been reported in [
49] a decrease in solar radiation of
when the dust deposition rate (DDR) is 425 mg/m
d in 60 days. On the contrary dry cities like Mexico City for also 60 days, the DDR reported is
to reduce the performance ratio up to
.
Detecting permanent faults such as PID, hotspots, or micro-cracks with online real-time methods is always a challenge. For instance, if a hotspot is considered mild, its temperature is just
C higher than the other parts of the cells; however, if the hot spot is considered severe, it presents a temperature higher by approximately 18
C. These facts mean that power production could be reduced between
and
[
50]. A similar analysis can be done with micro-crack; for instance, in a PV module formed with 60 cells, if just one cell has an inactive area of
, the power loss in the whole module is about
[
30]. Therefore, a rule to avoid false faults due to normal soiling or other temporal issues should use a
value between 5 and 10% as a threshold.