International Journal of Photoenergy
Volume 2015 (2015), Article ID 824832, 10 pages
http://dx.doi.org/10.1155/2015/824832

Research Article

Accurate Maximum Power Tracking in Photovoltaic Systems Affected by Partial Shading

Pierluigi Guerriero, Fabio Di Napoli, Vincenzo d’Alessandro, and Santolo Daliento

Department of Electrical Engineering and Information Technology (DIETI), University of Naples Federico II, Via Claudio 21, 80125 Naples, Italy

Received 28 November 2014; Accepted 19 January 2015

Academic Editor: Ying Dai

Copyright © 2015 Pierluigi Guerriero et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A maximum power tracking algorithm exploiting operating point information gained on individual solar panels is presented. The proposed algorithm recognizes the presence of multiple local maxima in the power voltage curve of a shaded solar field and evaluates the coordinated of the absolute maximum. The effectiveness of the proposed approach is evidenced by means of circuit level simulation and experimental results. Experiments evidenced that, in comparison with a standard perturb and observe algorithm, we achieve faster convergence in normal operating conditions (when the solar field is uniformly illuminated) and we accurately locate the absolute maximum power point in partial shading conditions, thus avoiding the convergence on local maxima.

1. Introduction

In the last few years photovoltaic (PV) installations have become familiar in the landscape of our countries. Their diffusion was initially encouraged by even aggressive, apparently uneconomical, feed-in tariffs; nevertheless the merit of those policies is now fully appreciable and consists in the impressive prices fall for all components of a PV plant, so much so that less than 2 /Wp are now required as total investment. According to [1] such low prices led the cost of PV energy to be comparable with conventional energy sources and facilitated self-sustainment of the PV market. However, some limiting factors, which are inherently embedded in the PV technology, still persist. One of them is the large dependence on shadows, as they cause dramatic degradation of energy production. This issue prevents PV diffusion in the urban context where neighbor buildings can likely affect daily irradiance actually available on rooftop sites. The effect of shadowing is the deformation of the power-voltage curve exhibited by affected subfield. As will be shown in the next sections, indeed, the power-voltage (P-V) curve of locally shaded strings (with the term string we indicate a group of series connected solar panels) is characterized by the presence of multiple local maxima; the consequence is that maximum power point tracking (MPPT) algorithms often stabilizes the PV system at an output power that is lower than the maximum achievable. This problem can be completely avoided by adopting distributed conversion schemes exploiting either microinverter [2–4] or solar panel optimizers [5–8]. Both of them lead each solar panel working at its own MPP; the drawback is that a complex power electronic circuit is brought on board solar panels, thus inducing a not negligible fault probability. As an alternative some centralized inverters allow a periodic scanning of the whole power-voltage curve which univocally individuates the absolute maximum on the curve; in this case a trade-off exists between the scan frequency and energy production because each time a scan is performed the system is driven far from the MPP; thus, when the string is uniformly irradiated the scan is definitively detrimental [9].

In this paper an approach that first identifies all maxima present in the power voltage curve and then drives the operating point of the solar system to the absolute maximum is described.

The procedure is based on the availability of a distributed sensor network devised to the measurement of the actual short circuit current () and open circuit voltage () for each panel forming the solar field. The algorithm for the tracking of the absolute maximum is based on the information achieved by the sensors; thus it is hereafter denoted as information based MPPT. The sensors allow on demand disconnection and by passing of solar panels from the string; electrical parameters are measured during the disconnection time (few milliseconds); thus they are not affected by the operating point of others solar panels belonging to the same string.

The paper is organized as follows. Section 2 reports the analytical model of the proposed approach. In Section 3 the distributed sensor network providing individual panels data is presented. In Section 4 numerical simulation results evidencing the effectiveness of the information based MPPT algorithm are reported. Section 5 shows experimental result concerning the comparison between the proposed approach and the behavior of a standard algorithm. Conclusions are drawn in Section 6.

2. First-Order Algorithm

The idea underlying the proposed MPPT algorithm assumes that it is possible to reconstruct an approximate version of the actual current-voltage (I-V) curve of a partially shaded string, formed by series connected solar panels, once the couple and is known for each panel. The algorithm estimates the voltage corresponding to the absolute maximum power point () and set this value as starting point for a standard perturb and observe (P&O) tracking procedure [10–13]. Reconstruction of the string I-V curve is performed by superposing a simplified form of the single panels I-V curves which are approximated by trapezoidal shapes. Namely, once both and of a given panel are known, the coordinates of the maximum power point are calculated according to the following relations:

where the coefficients and are taken from panel datasheets. Then the trapezoidal I-V curve of the solar panel is drawn as shown in Figure 1, where the actual I-V curve is also reported for comparison. It is worth noting that a critical point of this procedure is the strong temperature dependence of aforementioned coefficients; this implies that the working temperature of each solar panel is needed; actually and depend (weakly) on the irradiance as well. Both issues will be discussed in detail shortly.

Figure 1: Current-voltage characteristic of a PV panel (solid grey line) compared to both rectangular (dotted black line) and trapezoidal (solid black line) approximation. , , , and of the panel are also shown.

As an example of string characteristic reconstruction Figure 2 shows the hypothetical case where a string is formed by three solar panels, each of them, protected with one bypass power diode, exhibiting a different short circuit current (i.e., they are subjected to different irradiation levels). In particular Figure 2(a) shows the “actual” I-V and P-V curves evaluated by means of standard circuit simulations assuming a one diode model for the solar panel [14]; as it was expected three local maxima are observable in the P-V curve. Figure 2(b) reports the corresponding approximated curves; as can be seen all voltages corresponding to the maxima are almost coincident for the two cases.

Figure 2: (a) Simulated I-V (solid line) and P-V (dashed line) curve. (b) Trapezoidal approximation of both I-V (solid line) and P-V (dashed line) curve. Cycles indicate the local maxima.

The previous qualitative example can be easily generalized. Let us consider a string made by series connected solar panels and assume that the string is subjected to different irradiance levels. It is important to note that the number of local maxima appearing in the P-V curve coincides with the number of irradiance levels; thus we expect the presence of local maxima. In other words, each group of panels exhibiting the same contributes to the same local maximum. The voltage , corresponding to the th local MPP, can be computed as follows:

where is the open-circuit voltage exhibited by all panels belonging to group , which is defined as the panel group sharing the th irradiance level, is the number of panels forming group , and is defined as follows:

Equations (2) and (3) can be interpreted by referring again to the simplified case shown in Figure 2. Let us evaluate, for example, the voltage of the second maximum. As can be seen is almost equal to the sum of the voltages corresponding to the MPP of the first two panels because the third panel is bypassed. This condition is taken into account in (3) by setting ; indeed, the fact that the short circuit current of the third panel is lower than the others means that the third panel is bypassed.

Equation (2), however, does not account for the voltage drop across conducting bypass diodes [15, 16]. In the previous example, indeed, the actual is lower than the one given by (2) because the voltage drop across panel 3, coinciding with the forward voltage across the conducting bypass diode (about 0.7 V), must be considered. Thus a more detailed model for is

where is the forward voltage drop of bypass diodes. Equation (4) simply states that all modules that do not contribute to with their negatively affect with their .

Above equations can be conveniently expressed, for computational purposes, in a matrix form as follows:

where are -by- matrices defined according to (3).

With the aim to evaluate the power corresponding to each the respective current is also needed. The th can be straightforward written as

so that the power values of the maxima are expressed in a matrix form as

where is a -by- matrix defined by

The coefficients and depend on temperature and irradiance; they can be analytically determined as a function of both the open circuit voltage and the short circuit current as follows.

Let us consider the standard equation for a solar panel formed by series connected solar cells

By multiplying both members of (9) by the voltage the power supplied by the solar cell is achieved:

By definition is the voltage where (10) keeps its maximum; thus it can be calculated by deriving (10):

this derivative must be zero when :

By neglecting 1 in the square brackets (12) can be rewritten as

Equation (13)  can be further simplified by considering that as

By applying the logarithm to both sides, (14) becomes

In (15) we recognize the standard expression of the open circuit voltage

which allows the desired relation between and where the dependence on temperature is shown:

Equation (17) is graphically represented in Figure 3.

(a)

(b)

Figure 3: (a) MPP voltage against open-circuit voltage and (b) the percentage of with respect to the . The grey area highlights the typical range.

It is interesting to note that in the typical operating range for (0.6 V-0.7 V) is a fraction almost constant of ; from the measurement of a reliable estimation of can be achieved.

A similar procedure can be adopted to obtain as a function of .

By neglecting 1 in (9) we can write

hence

that, is zero when . Consider

which is the required relation between and . By numerically solving (20) the current at the MPP can be calculated from the measurement of .

In our approach the measurement of and is also exploited to gain the operating temperature of the solar panel and the irradiance. As it is well known, indeed, the short circuit current is a linear function of the irradiance while is an almost linear function (decreasing) of the temperature [17]. Respective coefficients are given in the data sheets accompanying the solar panels; hence the comparison of measured data with nominal values of and allows straightforward evaluation of actual illumination and temperature conditions.

3. The Disconnection Circuit

The algorithm discussed above requires the measurement of both and . These values are usually not accessible in series connected solar panel. In our system we exploit a devised electronic circuit which temporary (few milliseconds) disconnects each solar panel from the string, performs the measurements of own panel parameters, and transmits measured data via wireless communication. The system is schematically represented in Figure 4.

Figure 4: Schematic block diagram of the disconnection circuit.

As can be seen the disconnecting element is a power MOSFET which is series connected to the panel. When a measurement command is sent by the central coordinator the MOSFET is set in the OFF state and the solar panel behave as an isolated one. During the time that the panel is disconnected all parameters that are measured by the measurement unit are not affected by the operating point of the remaining string; hence the actual values of and can be measured. It is worth noting that during the disconnection time the string current is not interrupted thanks to the additional bypass diode which gives a free path for the string current. Just as an example Figure 5 shows an over day and collection for a solar panel embedded in a ten-panel string.

Figure 5: Individual panel short-circuit current (grey solid line) and open-circuit voltage (black solid line) against the time during the whole day.

A network controller (named PV Unit) collects individual panel operating parameters (open-circuit voltage, short-circuit current) measured by the distributed sensor network [18, 19]. The first stage of the information based MPPT controller receives panels data from the PV Unit and estimates the voltage () corresponding to the absolute maximum power point; then, the second stage, relying on an optimized perturb and observe (P&O) approach [10, 12], tracks the actual MPP voltage.

4. Numerical Simulations

The performances of the algorithm were tested by means of circuit level numerical simulations performed in the PSpice environment. A solar string formed by ten 50 Wp solar panels was designed. For reliable simulations each solar panel was divided into two subpanels, each protected by a devised bypass diode. The parameters of the solar panels were taken from the datasheets of the model ET-solar ET-M54050 [20], the same used for the experiments described in the next section.

Figure 6 shows the results we achieved when the string was subjected to four different irradiation levels, so as to have four local maxima in the P-V curve.

Figure 6: Current-voltage (a) and power-voltage (b) characteristics. The actual local maxima (crosses) and the maxima estimated by means of information based MPPT algorithm (circles) are also reported.

In the figure the crosses indicate the actual position of each maximum while circles indicate the position estimated by the tracking algorithm. As can be seen the match is almost perfect; hence, once the voltage corresponding to the absolute power maximum has been calculated the subsequent P&O algorithm can be effectively driven to track that maximum.

As it is well known the maximum power point shows a significant decrease when the temperature increases. This effect can be almost entirely attributed to the dependence on temperature of the open circuit voltage

where the term is given by

Thus, the key parameters affecting the temperature dependence of PV modules are the reverse saturation current of the p-n junction and recombination lifetime . Direct measurement of these parameters for analytical determining the temperature dependence of the MPP is quite complicated [21–24]; thus, in order to verify that the direct measurement of single panel and always leads to a close estimation of the actual , independently of operating conditions, numerical simulations were performed in a wide range of temperature (between 10°C and 80°C) and irradiance levels (715 different shading shapes obtained by shading up to 4 panels in the string with 4 different shading levels, corresponding to 100 W/m², 200 W/m², 300 W/m², 500 W/m², and 800 W/m²). Results were collected in plots like the one shown, as an example, in Figure 7. The figure depicts the relative voltage error in local maxima estimation corresponding to all the 715 shading cases at .

Figure 7: Relative voltage error in local maxima estimation. Cycles represents the estimation error with respect to the actual numerical result corresponding to all the 715 shading cases at .

The figure reports the percentage error on the estimation of as a function of the actual for all simulated cases. Two important comments can be made. First, the absolute error is always less than 2.5%, which is a very small value. Second, the error is more likely an underestimation more than an overestimation. The results achieved at higher temperatures, which are likely reached during normal operation of solar panels, do not differ from Figure 7. They are summarized in Figure 8 where the percentage error of the estimation and the corresponding error bar are reported for each simulated temperature.

Figure 8: The percentage error on estimation and the corresponding error bar for each considered temperature value.

5. Experiments

Experiments were conducted with the aim of both experimentally verifying the correct behavior of the whole system and comparing the proposed algorithm with the operation of a commercial inverter. Figure 9 shows the solar string exploited for the experiments.

Figure 9: Ten-panel PV string. The desired shading profile is obtained superposing different kind of optical plastic film on selected PV panels.

It is composed by ten solar modules whose parameters were described in the previous section. As can be seen an uneven illumination frame was induced by covering some solar cells with semitransparent dyed sheets.

Two cases were considered. In the first case the string was kept under uniform illumination. The resulting P-V curve we measured is shown in Figure 10(a).

(a)

(b)

Figure 10: Sunny case. (a) Experimental power-voltage characteristic (solid line) of the 10-panel test string under full irradiation. The operating points of the string (triangles) by adopting a classic P&O algorithm are shown. (b) String power against the time.

The evolution of the output voltage imposed by the tracking algorithm of the commercial inverter adopting a conventional P&O algorithm is shown in Figure 10(b).

The algorithm assumes that the MPP lies at about 80% of ; this value is assumed as starting point for the “hill climbing.” At the end the operating point stabilizes at around 180 V, that is, the actual voltage of the MPP.

Figure 11 shows the behavior of the tracking algorithm on the same P-V curve.

(a)

(b)

Figure 11: Sunny case. (a) Experimental power-voltage characteristic (solid line) of the 10-panel test string under full irradiation. The operating points of the string (triangles) by adopting the information based MPPT algorithm are shown. (b) String power against the time.

In this case the starting voltage was estimated from the knowledge of and and was immediately achieved. The two algorithms are compared in Figure 12 which reports the output power achieved by each of them when uniform illumination conditions are assumed.

Figure 12: Comparison of string power behavior against the time corresponding to both P&O (dashed line) and information based MPPT (solid line) algorithms.

As can be inferred the proposed information based MPPT keeps the MPP faster.

A more interesting case is the second one we considered, where illumination was not uniform, so as to have the P-V curve shown in Figure 13(a).

(a)

(b)

Figure 13: Partial-shading case. (a) Experimental power-voltage characteristic (solid line) of the 10-panel test string under partial-shading conditions. The operating points of the string (triangles) by adopting a classical P&O algorithm are shown. (b) String power against the time.

Also in this case we compared the performances of the standard algorithm with the proposed one. Figure 13(b) shows the output voltage given by the standard one.

Again, the algorithm starts from and tries, as starting point, 80% of that value; then the voltage is increased to stabilize at about 200 V. By comparing with the P-V curve we observe that the algorithm tracked the local maximum on the right side of the PV curve.

Figure 14 shows the operation of the information based MPPT.

(a)

(b)

Figure 14: Partial-shading case. (a) Experimental power-voltage characteristic (solid line) of the 10-panel test string under partial-shading conditions. The operating points of the string (triangles) by adopting the information based MPPT algorithm are shown. (b) String power against the time.

In this case the starting point is a slight underestimation of the actual absolute maximum which is immediately reached by the system as shown in Figure 14(b).

Figure 15 shows the comparison performed on the above P-V curves.

Figure 15: Algorithms comparison. The tracking performance corresponding to the commercial inverter (a) and the information based MPPT (b) are compared. The operating points (triangle) over the P-V curves (solid lines) are reported.

The most important observation is that the information based MPPT correctly recognizes the absolute maximum thus giving a power gain of about 75 W over the standard system.

6. Conclusions

Local information about the operation of each panel embedded in a PV string can strongly improve the overall efficiency of a PV system. In this paper an information based MPPT algorithm, exploiting the information provided by individual panel sensor network, has been presented.

Experiments comparing the operation of a standard P&O algorithm with the proposed approach demonstrate the capability of the present approach to locate the true MPP, thus avoiding power losses due to tracking failures. Moreover, a faster convergence has been reached with respect to a classic P&O algorithm.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

G. Masson, S. Orlandi, M. Rekinger et al., Global Market Outlook for Photovoltaics 2014–2018, European Photovoltaic Industry Association, 2014, http://www.epia.org/.

D. Meneses, F. Blaabjerg, Ó. García, and J. A. Cobos, “Review and comparison of step-up transformerless topologies for photovoltaic AC-module application,” IEEE Transactions on Power Electronics, vol. 28, no. 6, pp. 2649–2663, 2013.

S. B. Kjaer, J. K. Pedersen, and F. Blaabjerg, “A review of single-phase grid-connected inverters for photovoltaic modules,” IEEE Transactions on Industry Applications, vol. 41, no. 5, pp. 1292–1306, 2005.

M. Coppola, S. Daliento, P. Guerriero, D. Lauria, and E. Napoli, “On the design and the control of a coupled-inductors boost dc-ac converter for an individual PV panel,” in Proceedings of the 21st International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM '12), pp. 1154–1159, June 2012.

N. Femia, G. Lisi, G. Petrone, G. Spagnuolo, and M. Vitelli, “Distributed maximum power point tracking of photovoltaic arrays: novel approach and system analysis,” IEEE Transactions on Industrial Electronics, vol. 55, no. 7, pp. 2610–2621, 2008.

G. Petrone, G. Spagnuolo, and M. Vitelli, “An analog technique for distributed MPPT PV applications,” IEEE Transactions on Industrial Electronics, vol. 59, no. 12, pp. 4713–4722, 2012.

B. V. P. Chong and L. Zhang, “Controller design for integrated PV-converter modules under partial shading conditions,” Solar Energy, vol. 92, pp. 123–138, 2013.

S.-M. Chen, T.-J. Liang, and K.-R. Hu, “Design, analysis, and implementation of solar power optimizer for DC distribution system,” IEEE Transactions on Power Electronics, vol. 28, no. 4, pp. 1764–1772, 2013.

M. Boztepe, F. Guinjoan, G. Velasco-Quesada, S. Silvestre, A. Chouder, and E. Karatepe, “Global MPPT scheme for photovoltaic string inverters based on restricted voltage window search algorithm,” IEEE Transactions on Industrial Electronics, vol. 61, no. 7, pp. 3302–3312, 2014.

D. P. Hohm and M. E. Ropp, “Comparative study of maximum power point tracking algorithms,” Progress in Photovoltaics: Research and Applications, vol. 11, no. 1, pp. 47–62, 2003.

M. A. G. de Brito, L. Galotto, L. P. Sampaio, G. de Azevedo Melo, and C. A. Canesin, “Evaluation of the main MPPT techniques for photovoltaic applications,” IEEE Transactions on Industrial Electronics, vol. 60, no. 3, pp. 1156–1167, 2013.

K. Ishaque and Z. Salam, “A review of maximum power point tracking techniques of PV system for uniform insolation and partial shading condition,” Renewable & Sustainable Energy Reviews, vol. 19, pp. 475–488, 2013.

N. Onat, “Recent developments in maximum power point tracking technologies for photovoltaic systems,” International Journal of Photoenergy, vol. 2010, Article ID 245316, 11 pages, 2010.

J. Cubas, S. Pindado, and A. Sanz-Andrés, “Accurate simulation of MPPT methods performance when applied to commercial photovoltaic panels,” The Scientific World Journal. In press

V. D’Alessandro, P. Guerriero, and S. Daliento, “A simple bipolar transistor-based bypass approach for photovoltaic modules,” IEEE Journal of Photovoltaics, vol. 4, no. 1, pp. 505–513, 2014.

S. Daliento, L. Mele, P. Spirito, R. Carta, and L. Merlin, “Experimental study on power consumption in lifetime engineered power diodes,” IEEE Transactions on Electron Devices, vol. 56, no. 11, pp. 2819–2824, 2009.

S. Daliento, V. d'Alessandro, and P. Guerriero, Dispositivi e sistemi fotovoltaici, Edises, 2014.

M. Gargiulo, P. Guerriero, S. Daliento et al., “A novel wireless self-powered microcontroller-based monitoring circuit for photovoltaic panels in grid-connected systems,” in Proceedings of the International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM '10), pp. 164–168, June 2010.

P. Guerriero, V. d'Alessandro, L. Petrazzuoli, G. Vallone, and S. Daliento, “Effective real-time performance monitoring and diagnostics of individual panels in PV plants,” in Proceedings of the 4th International Conference on Clean Electrical Power (ICCEP '13), pp. 14–19, June 2013.

“Datasheet of ET-solar PV module ET-M54050,” http://www.etsolar.com.ar/

S. Bellone, G. D. Licciardo, S. Daliento, and L. Mele, “Experimental measurements of majority and minority carrier lifetime profile in SI epilayers by the use of an improved OCVD method,” IEEE Electron Device Letters, vol. 26, no. 7, pp. 501–503, 2005.

S. Daliento, A. Sanseverino, and P. Spirito, “An improved model for extraction of strongly spatial dependent lifetimes with the ac lifetime profiling technique,” IEEE Transactions on Electron Devices, vol. 46, no. 8, pp. 1808–1810, 1999.

S. Daliento, L. Mele, E. Bobeico, L. Lancellotti, and P. Morvillo, “Analytical modelling and minority current measurements for the determination of the emitter surface recombination velocity in silicon solar cells,” Solar Energy Materials & Solar Cells, vol. 91, no. 8, pp. 707–713, 2007.

S. Daliento and L. Mele, “Approximate closed-form analytical solution for minority carrier transport in opaque heavily doped regions under illuminated conditions,” IEEE Transactions on Electron Devices, vol. 53, no. 11, pp. 2837–2839, 2006.