Power Output Fluctuations in Large PV Plants
Power Output Fluctuations in Large PV Plants
net/publication/316820035
CITATIONS READS
6 215
4 authors, including:
Miguel García
Universidad Pública de Navarra
37 PUBLICATIONS 942 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Eduardo Lorenzo Pigueiras on 06 October 2017.
Abstract. The power generated by large plants located in Spain and with output powers
photovoltaic (PV) plants is of an intermittent nature, ranging from 1 MWp to 9.5 MWp (a total of 20
due to cloud passage, which can create considerable MWp). The research focussed on studying the
fluctuations in output. Therefore, if the penetration of irradiance and power fluctuations for time intervals of
PV power in the overall electricity system is high less than 10 minutes, paying particular attention to the
enough, these fluctuations could have a serious effect characterisation of the magnitude and frequency of
on the quality and safety of the power supply. This occurrence of these fluctuations. The study also
article aims to characterise these fluctuations and to analysed the influence of the PV plant size and the
look at the smoothing effect of the PV plant size and number of PV plants grouped in a multi-site system,
that of the geographical dispersion of the PV plants. on the magnitude of the fluctuations. Finally, a PV
plant was modelled, making it possible to simulate the
Key words: dynamics of the power generated by a PV plant based
on irradiance data.
large PV plants; grid-connected; power fluctuations;
PV plant model; geographical dispersion. 2. Experimental set-up
1s 10
2
-2
Rada 1780 1400 8,7 10
-3
10
Castejón 2640 2.000 11,8
-4
10
Milagro 9500 7243 52
-5
10
Socuéllamos 2600 1975 18
-6
10
-100 -80 -60 -40 -20 -1 +1 +20 +40 +60 +80 +100
Total 19906 15348 - Irradiance Fluctuations (%)
∆P∆t (t ) =
[P (t + ∆t ) − P (t ) ] (1)
irradiance fluctuations recorded at the various sites for
an entire year in relation to the different sampling
P* times ∆t considered. The value of the maximum
fluctuation is extremely useful for the power utility
This definition can also be applied to irradiance operator, making it possible to check the robustness of
fluctuations, ∆G∆t, normalizing by G* = 1000 W/m2. a grid. The values observed for ∆t greater than 2
seconds are over 80%. Therefore, although the
80
100
Max-48 kWp
60 Max-143 kWp
80 Max-958 kWp
40 Max-990 kWp
Max-Arguedas Max-1438 kWp
60
Max-1780 kWp
Max-Castejon 20 Max-2640 kWp
Max-Cintruenigo Max-9500 kWp
40
Max-Milagro 0
Max-Rada 1 2 5 10 20 40 60 120 300 600
20 Max-Sesma Sample period (s)
Max-Socuellamos Fig. 5: Maximum power fluctuations recorded at each PV
0
1 2 5 10 20 40 60 120 300 600 plant vs. ∆t.
Sample Time (s)
Fig. 3: Maximum irradiance fluctuation recorded at each PV The relationship between thee smoothing effect
plant site vs. ∆t.
and PV plant size is embodied in the following
With regard to the power fluctuations produced by proposed empirical Eq. (2), which determines the
a PV plant, the magnitude of these fluctuations is peak power fluctuation value of a PV plant of any size
influenced not only by ∆t, but also by the PV plant S [Ha] and for any sampling time ∆t:
size. It is therefore to be expected that, the larger the
Max(∆P∆t) =90%·(1-e-0.24·∆t)·S –c (2)
plant size, the lower the power fluctuations in relation
to the irradiance fluctuations. This reasoning is based Smoothing coefficient c for ∆t= 1s is 0.5, thus the
on the fact that a cloud will always take longer to smoothing effect on the fluctuations can be described
cover an entire PV plant than an irradiance quite accurately by a 1/√S law. However, for ∆t=600s,
measurement module. Thus, the shorter the sampling coefficient c is equal to 0.02 which supports the idea
time, the greater the smoothing effect. This fact can that, for these sampling times, the plant size is not a
clearly be observed in Fig. 4, showing the annual smoothing factor for the fluctuations. A more detailed
power fluctuation distributions for the powers and supplementary account of the results shown in
recorded at each PV plant for ∆t = 20s. As was to be this and the previous section can be found in [8].
expected, the filtering effect on the fluctuations
increases in line with the increased plant size whilst, 5. Fluctuations in a multi-site PV plant
on the other hand, it decreases as ∆t increases.
system
Likewise, this smoothing effect can be clearly seen in
Fig. 5, showing the value of the maximum power
We shall now go on to quantify the smoothing of PV
fluctuation recorded over a year at each plant in
power fluctuations by geographical dispersion.
relation to ∆t. However, for ∆t= 600s, the smoothing
Thanks to the synchronised monitoring made, it is
effect is non-existent, and fluctuations of 90% can be
possible to determine, second by second, the
reached, regardless of the plant size.
combined power generated by a group of PV plants,
10
2 by monitoring plant behaviour, in the same way as the
48 kWp 10
2
system operator.
1 143 kWp
10 958 kWp Fig. 6 shows the normalised distributions for the
0
990 kWp
1438 kWp
fluctuations recorded in 2009 for ∆t = 60s
Relative frequency (%)
10 0
1780 kWp 10
-10 -1 +1 +10 corresponding to just one plant (Arguedas,
2640 kWp
10
-1
9500 kWp P1* = 0.775 MW), to a group of 3 plants (Castejón,
48 kWp
Milagro and Sesma, P3* = 10.043 MW, COMB3) and
-2
10 a further group of 6 plants (all located in southern
-3
Navarre, P6* = 13,373 MW, COMB6). It can clearly
10
9500 kWp be observed how the magnitude and frequency of the
10
-4 fluctuations decreases in line with an increase in the
number of plants grouped together, N. It is important
-100 -80 -60 -40 -20 -1 +1 +20 +40 +60 +80 +100
to emphasise the fact that the smoothing effect
Power Fluctuations (%) extends as far as the long sampling times (∆t = 600 s)
in contrast to the effect observed for plant size.
Fig. 4: Distributions of the power fluctuations, ∆P∆t(t),
recorded at each PV plant over an entire year, for ∆t = 20 s.
0 0
Relative frequency (%) 10 10
-2
-2 Arguedas 10
10
COMB3 COMB6
-4
10
-4
10 1s
60s 20s
600s
-6
-100 -80 -60 -40 -20 -1 +1 +20 +40 +60 +80 +100 10
-100 -80 -60 -40 -20 -1 +1 +20 +40 +60 +80 +100
Power Fluctuations (%) Power Fluctuations (%)
Fig. 6: Distributions of the power fluctuations ∆P∆t,N, Fig. 8: Distribution of the power fluctuations ∆P∆t,2,
recorded in the course of an entire year (2009) at the recorded over an entire year (2009) at the combination of
Arguedas PV plant and the combination of N = 3 (COMB3) Arguedas–Castejón PV plants (distance = 6 km) and
and N = 6 (COMB6) for multi-site PV plants, all in Navarre, Arguedas–Socuéllamos (distance = 345 km) for ∆t = 1, 20,
for ∆t = 60s. 60, and 600 s.
Fig. 7 below shows the value of the peak power The data shown in Fig. 7 was used to obtain an
fluctuation recorded versus the sampling time ∆t empirical Eq. (3) which provides the value of the peak
based on the number of plants grouped together, N, power fluctuation for a group of N PV plants of a
for all plants located in Navarre. It is possible to certain size S, for any ∆t of less than 10 minutes.
observe a considerable reduction in the magnitude of
the fluctuations, irrespective of the ∆t. Furthermore, it 99th(∆P∆t,N)= 99th(∆P600,1)·(1-e-0.24·∆t)·S -c·N–a; (3)
was found that for a given N, no plant combination a,c> 0
had a greater smoothing effect than another. In other
words, the location of each PV plant in relation to the This equation describes the smoothing of the
others within the group, had no effect on the fluctuations resulting from the plant size and also the
magnitude of the smoothing effect. Therefore, the geographical dispersion. Whenever coefficients a are
fluctuations of a group of PV plants will be smoothed greater than c, as in Table 4, the smoothing due to
provided that there is sufficient distance from one dispersion is greater than that due to plant size.
plant to another. This leads on to the question of what
is the minimum distance between plants in order to ∆t c a
ensure that the plant fluctuations remain independent 1s 0.49 25.55
of each other and to guarantee smoothing. For this, 5s 0.29 64.09
20s 0.11 88.2
Fig. 8 shows the distribution of fluctuations for the
60s 0.05 91.6
entire year of 2009 for the Arguedas-Castejón group, 600s 0.02 94.07
separated by a distance of 6 kms, and the distribution
for the Arguedas-Socuéllamos pair, 345 kms apart. As Table 4: Estimated parameters a and c for empirical
can be seen, there are hardly any differences between Eq. (3).
the fluctuations of any two-plant combinations.
Therefore, a distance of 6 kms, being the shortest Based on Eq. (3), Fig. 9 shows the fluctuations of
distance between any plant and our observatories a total PV power of 100 MW in relation to the number
(located at Arguedas and Castejón) is sufficient to of plants amongst how this power is distributed and
guarantee the independence of plant power sampling time ∆t. As can be seen, after approximately
fluctuations for time intervals of less than 10 minutes 100 plants (1 MW each) the fluctuations are smoothed
or, to put it another way, to guarantee smoothing by considerably, and it is debatable whether the power
geographical dispersion. needs to be distributed still further, given the fact that
this would have little effect on the smoothing
achieved. All these conclusions, together with some
additional analyses can be found in [9].
-3
maximum daily fluctuations, for ∆t equal to 600s at
10
Milagro PV plant. As can be seen, the frequency
-4 distributions are very similar. A more detailed
10 GN
description of this model can be found in [10].
PN 0.99 MWp
-5
10
PN 9.5 MWp 40
Real
-6
10 -6 -4 -2 35 Simulated
10 10 10
f (Hz)
30
Relative Frequency (%)
performance, with a downward trend f -0.7. However, Fig. 12: Histograms of the real maximum daily power
after a certain frequency fc, which differed from one fluctuation observed in the course of a year at Milagro PV
plant to another, a change in trend was observed, plant and the simulated fluctuations for ∆t =600 s.
implying the increased smoothing of the higher
frequencies. As was to be expected, the greatest 7. Conclusions
smoothing was obtained for the plant with the largest
size. The mathematical expression of the new trend is The quality and safety of the electricity supply can be
form f -1.7 and is independent of the PV plant seriously affected by large PV plant power
considered. The cut-off point for functions f -0.7 fluctuations caused by the variability in irradiance due
and f -1.7 can be identified as cut-off frequency fc for to cloud passage. This article has presented the study
the systems in question. Furthermore, the difference of these fluctuations, obtaining the following results:
between the slope for the two functions is the unit,
which means that the corresponding low pass filter is - A unique monitoring system, in place since April
first order, with a pole at cut-off frequency fc. This 2008, and which makes synchronous recordings with
exercise was extended to the other PV plants, a one second resolution, of the power delivered to the
obtaining the corresponding value for fc. These grid by seven PV plants, with powers ranging from 1
experimental values were generalised to obtain MWp to 9.5 MWp and with a total combined power
of almost 20 MWp.
Acknowledgements