Control Strategies to Smooth Short-Term Power Fluctuations in Large Photovoltaic Plants Using Battery Storage Systems

The variations in irradiance produced by changes in cloud cover can cause rapid fluctuations in the power generated by large photovoltaic (PV) plants. As the PV power share in the grid increases, such fluctuations may adversely affect power quality and reliability. Thus, energy storage systems (ESS) are necessary in order to smooth power fluctuations below the maximum allowable. This article first proposes a new control strategy (step-control), to improve the results in relation to two state-of-the-art strategies, ramp-rate control and moving average. It also presents a method to quantify the storage capacity requirements according to the three different smoothing strategies and for different PV plant sizes. Finally, simulations shows that, although the moving-average (MA) strategy requires the smallest capacity, it presents more losses (2–3 times more) and produces a much higher number of cycles over the ESS (around 10 times more), making it unsuitable with storage technologies as lithium-ion. The step-control shown as a better option in scenery with exigent ramp restrictions (around 2%/min) and distributed generation against the ramp-rate control in all ESS key aspects: 20% less of capacity, up to 30% less of losses and a 40% less of ageing. All the simulations were based on real PV production data, taken every 5 s in the course of one year (2012) from a number of systems with power outputs ranging from 550 kW to 40 MW. OPEN ACCESS Energies 2014, 7 6594


Introduction
The concern shown by a number of transmission system operators (TSO) over the increasing presence of large photovoltaic (PV) plants is leading to constant changes in grid codes, directed at facilitating PV penetration yet without detracting from the safety and quality of supply.The main reason is the fluctuating nature of PV generation due to variations in irradiance as a result of changes in cloud cover.Specifically, for very fast fluctuations, the TSO has an extremely limited response capacity and, if the fluctuations exceed the permitted limits, then there is a risk of a complete power system failure.In fact, some countries are currently working on grid codes that take this problem into account [1][2][3].This risk is particularly significant for high PV penetration rates, such as islands.For example, the Puerto Rico Electric Power Authority [1] requires the PV plant fluctuations to be less than 10%/min.In other countries such as Mexico, the regulations target greater restrictions, of around 1%/min-5%/min [2].The speed and magnitude of these fluctuations are well-known and have been quantified not only with irradiance data [4][5][6][7][8][9] but also with power output measurements [10][11][12].Likewise, the researchers [5,13,14] have demonstrated that the larger the PV system is, the lower the PV fluctuations are: typically, for this 1 min time window, a 1 MWp PV plant records fluctuations of up to 90%, a 9.5 MWp plant easily exceeds 70% [14], whilst a 45.6 MWp fluctuates by up to 33% [15].Therefore, as these new regulations gradually come into force, it will be necessary to equip PV plants with an energy storage system (ESS) in order to smooth out any fluctuations to within the maximum limits.
Logically, the installation of an ESS has a major impact on the energy/economic balance of the PV system, playing a key role in the viability of the future PV systems due to their high costs and reduced shelf life.As a result, parameters such as energy capacity, losses and the cycling degradation of the ESS, take on particular importance.Any reduction in both the ESS capacity required and charge/discharge cycles will have a positive impact on reducing the investment required to install and maintain the ESS.Consequently, the control strategy selected in order to smoothing fluctuations will be a crucial decision.Nevertheless, given the maximum fluctuation limitation, there is a range of control strategies to reduce the fluctuations to under this limit.At present, the two most often proposed strategies in the literature are ramp-rate control [16][17][18][19][20][21] and the moving-average (MA) [22][23][24][25].The main advantage of MA is that, if the system is equipped with an ideal converter and battery, by definition of the mean value, the value of the energy in the ESS at the beginning and end of any given day should be the same.In real practice, at the end of the day, the battery is discharged to a value equal to the energy lost in the charging/discharging processes.On the other hand, the key point of the ramp-rate control is that only acts when the fluctuation exceeds the maximum allowable ramp-rate value (i.e., lower degradation).There are also other strategies available such as the constant production strategy [26,27], however, this goes beyond smoothing out the fluctuations in short periods of time and requires a much larger ESS.Even so, no individual analysis is made for each of these two strategies as to the impact of them on the actual ESS requirements: energy capacity; charge/discharge cycles; cycling degradation; losses in the power electronic converter (PEC); losses in the ESS, etc.Only Marcos et al. [15] propose a method for sizing the ESS (with regard to power output and energy), but exclusively for ramp control.Accordingly, few comparisons have been made between the performance of the strategies ramp-rate control and the moving average working in the same PV plant-ESS system.The studies conducted are primarily focused on programming, monitoring and implementing the strategy [22,[28][29][30].
Likewise, these two strategies do not take into account a well-known smoothing effect over the fluctuations derived from the dispersion of PV plants: in a scenario of distributed generation, short-term power fluctuations of a group of PV plants are reduced linearly by the square root of the number of plants [5,31,32].With the goal of take advantage of this natural behavior, we first propose a novel control strategy (step-control) to improve the results in relation to the state-of-the-art strategies, ramp-rate control and MA, in terms of the performance and cyclability of the ESS.Secondly, we describes the method in order to calculate, for any PV plant size and maximum allowable ramp-rate, the energy storage requirements for both MA and step-control strategies.Finally, we have compared these three control methods in order to fit the grid quality demands.Each control method is evaluated regarding the quality of produced signal and the demand on the ESS.All the simulations were based on actual PV production data, taken every 5 s in the course of one year (2012) from a number of systems with power outputs ranging from 550 kW to 40 MW.The experimental ESS comprised a lithium-ion battery connected at the AC side through a two-way DC/AC PEC.The results shown can be used as a guide to select the best control strategy, and the optimum design and selection of the ESS, making it possible to assess the impact of the ESS on the PV plant energy/economic balance.

Database
Thanks to extensive monitoring, 5 s synchronized records of the output power of all the inverters from the Amaraleja (South Portugal, 38.1900 N, 7.2108 W) PV plant are available from May 2010.This plant, the property of Acciona Energí a S.A. (Sarriguren, Spain), occupies an area of 250 Ha and includes 2520 solar trackers with a rated output of 17.7-18.8kWp, up to a total peak power of 45.6 MWp (inverter power, P*, is 38.5 MW).Each tracker is tilted at an angle of 45° in relation to the ground with vertical tracking in an east-west direction (azimuthal).The plant is divided into 70 units, each comprising 36 tracking systems connected to a 550 kW DC/AC inverter.The minimum and maximum distances between the units are 220 m and 2.5 km, respectively.From this work, data was taken not only of the entire PV plant but also of five sections with P* between 0.55 kW and 11.5 MW (Figure 1), making it possible to study the relationship between the storage requirements and the size of the PV power plant.The simulations presented in this paper are based on data recorded over the course of one year (2012).

Power Fluctuations with No Energy Storage
Power fluctuation ΔPΔt(t) is defined at moment t for a given sampling period, Δt, as the difference between two consecutive samples of the power generated PPV(t), normalized to the rated power of inverter P* of the plant in question, Equation (1), which is: Therefore, given a power output time series PPV(t) as recorded in our experiment, it is possible to calculate an associated time series of fluctuations ΔPΔt(t); resulting from applying Equation (1) with a given time window Δt and moving it at the same rate as the monitoring resolution (5 s in our case).Having obtained vector ΔPΔt(t), it is then possible to determine the amount of time the fluctuations exceed a given permitted maximum ramp value in the course of one year.This result can be seen in Figure 2, showing the frequency of occurrence of fluctuations greater than a certain ramp rmax, in relation to the total production time for one whole year (2012).As was to be expected, significant fluctuations decrease as rmax increases.For example, for a restriction such as [2], where rmax is equal to 2%/min, even for the 38.5 MW PV system, more than 10% of the total production time exceeds the maximum ramp condition.This example shows that a battery would be required in order to comply with any power ramp limits imposed.

Figure 2.
Frequency over one year (2012) of PV power fluctuations calculated in 1-min time window, ∆P1min(t), are higher than a given ramp r (%/min).The frequency value is given in relative terms to the total production time (4380 h).

A Generic Strategy for Smoothing Fluctuations through Energy Storage
For our simulations, a lithium ion battery-based storage system was selected.Logically, the methodology put forward is equally valid for any other type of storage system (such as supercapacitors, flow batteries, fly wheels, etc.). Figure 3 shows the model used to simulate a generic strategy to smooth out fluctuations.For each sampling period, the algorithm reduces the fluctuations of PPV(t) to the maximum permitted level rmax and the result is the power injected into the grid, PG(t).The difference between PG(t) and PPV(t) is the power setpoint to either be provided by (PESS > 0) or to be absorbed by the ESS (PESS < 0).The next step is to apply the performance of the PEC ηPEC and that of the lithium ion battery, ηBAT, in order to obtain the battery power value PBAT required.
Figure 3. Generic model for the evaluation of the fluctuation smoothing strategies at a PV plant with energy storage.As can be seen, the performance of the energy storage system (ESS) has been divided into two: on the one hand, the performance of the power converter (ηPEC) and, on the other hand, that of the lithium ion battery (ηBAT).
To calculate the PEC losses, a real curve obtained from the Ingecon Sun Powermax ® 1 MW [33] commercial converter was used, Figure 4.This same converter is connected to a 1 MW transformer on the AC side, whose performance curve is also shown in Figure 4, together with the total for the PEC (Ingecon Sun + Transformer).With regard to the lithium ion battery losses, data was used from the manufacturer's catalogue, specifically the 1 MW/560 kW• h Intensium® Max battery made by SAFT (Bagnolet, France), whose round-trip efficiency can be considered constant and equal to 95% [34].The integral of PBAT is the energy EBAT to be provided by the battery in order to smooth the fluctuations.We would emphasize that this paper does not seek to achieve the accurate modelling of the performance of the ESS, but rather to make a comparison of the different strategies.

Ramp Control Strategy
Given a maximum allowable ramp condition for the power injected rmax, at each simulation step t the inverter tries to inject all the power to the grid, that is PG(t) = PPV(t).The control is enabled if the difference with the previous sample PG(t − Δt) exceeds the maximum ramp condition, that is, Equation (2): where rmax is expressed in units of power by time.This control has been improved by a state of charge (SOC) control proposed by [15] which prevents the natural tendency of the battery to discharge due to the effect of successive fluctuations (Figure 5).In this way, it is possible to maintain EBAT around the reference value EBAT,ref.The control action is applied prior to the ramp-rate limiter to ensure that condition rmax is met at all times.For our simulations, we will work with K = 6, a value which allows a good trade-off between speed and stability.Lower values might increase the risk of a complete discharge of CBAT, while higher values would result in extensive ESS usage.K values between 2 and 8 are recommended.As it was explained in [15], EBAT,ref should be half of the ESS capacity, CBAT.Note that the SOC control action is also smoothed by the ramp-limiter in order to guarantee that power fluctuations are always below rmax.
That previous work also proposed a model to determine CBAT (Figure 6): basically, the worst fluctuation for a PV plant can be modelled as the response of a first order system with a time constant τ in relation to a 90% irradiance fluctuation G(t) (global irradiance to diffuse irradiance).Furthermore, τ is dependent on the size of the PV plant (P*).Then, the energy required EBAT,ramp to smooth the worst fluctuation with the ramp-rate control strategy is equal to the integral of the difference between PPV(t) and PG,ramp(t), in other words, Equation (3): where P* is expressed in (kW), rmax in (%/s), τ in (s), and EBAT,ramp in (kW• h).Furthermore, as the sign of the first fluctuation (positive or negative) is unknown, battery CBAT,ramp must be sized for double this power, Equation ( 4): where P* is expressed in (kW), rmax in (%/s), τ in (s), and CBAT,ramp in (kW• h).
Figure 7a shows, for 15 February 2012, the evolution of the power generated PPV of Section B (1100 kW) and the simulated power output that would be injected to grid PG through this algorithm for rmax = 2%/min, a similar restriction to that of [2].According to Equation ( 4), CBAT needed to smooth the worst fluctuation is equal to 746 kW• h, so EBAT,ref = 373 kW• h. Figure 7b plots the corresponding PBAT for that specific day, whilst Figure 7c plots the integral over time EBAT.It can be seen how the SOC control operates correctly, returning the battery energy to its initial state.For this day, the required battery power is PBAT,MAX = 937 kW (or PBAT,MAX = 0.85• P*) and the required effective storage capacity is CBAT = EBAT,MAX − EBAT,MIN = 433 kW• h (or 23 min of capacity, equivalent to 0.39 h of PV plant production at P*).The total energy that has passed through the battery is 1600 kW• h, equivalent to 54% of daily production, whilst the total losses (battery plus PEC) are 3.87% of production.

Moving-Average Strategy
Given a power time series PPV(t) as recorded in our experiment, the power smoothed to be injected into the grid PG(t) is calculated as the mean production value in a time window with a duration of T, in other words, Equation ( 5): The greater the time for window T, the greater the smoothing of the fluctuations at PG(t).The key advantage of this strategy is that, if the system is equipped with an ideal converter and battery (ηBAT = ηPEC = 1), by definition of the mean value, the value of EBAT at the beginning and end of any given day should be the same.Therefore, there is no need for any type of SOC control to prevent the continuous battery discharge.However, if ηBAT and ηPEC <1, then this advantage disappears and, at the end of the day, the battery is discharged to a value equal to the energy lost in the charging/discharging processes in the ESS.This same phenomenon was already observed in [35].This paper proposes a simple solution to this problem, consisting in offsetting the mean value of the losses, also for a time window T (Figure 8).In the same manner as was done in ramp-rate control, it would be important to obtain a general expression to determine the storage capacity needed by the MA strategy.Figure 9 yet again resolves this problem.The analytical expression of function PG,MA is equal to Equation ( 6): Substituting T = 5400/r and resolving Equation ( 6), we obtain the general expression of PG,MA, Equation (7): The final term of this expression is responsible for the phase lag existing between PG,MA and PG,ramp which can be seen in Figure 9.The area between both curves AMA,ramp can approximately be calculated by Equation ( 8 where P* is expressed in (kW), rmax in (%/min), τ in (s), and AMA,ramp in (kW• h).Therefore, the MA strategy requires a battery equal to Equation ( 9  Figure 10 shows the example of applying this strategy once again for 15 February 2012, and with T = 2700 s.For this time window T, the maximum daily fluctuation in 1 min, is reduced to 1.95%/min.For this strategy, the required battery power is PBAT,MAX = 919 kW (or PBAT,MAX = 0.83• P*) and the required effective battery capacity is CBAT = EBAT,MAX − EBAT,MIN = 378 kW• h (or 20 min of capacity, equivalent to 0.34 h of PV plant production at P*).We would emphasise that the SOC control tends to return the battery power to its original state, thanks to the offsetting of losses.Likewise, although PG is slightly smoother than in the case of ramp-rate control, the ESS endures a few more cycles.The total energy that has passed through the battery is 2052 kW• h, equivalent to 70% of daily production, whilst the total losses in the battery and PEC are 5.2% of production.We should emphasize that this strategy must be operative even on clear days, because the reduction in fluctuations must be ensured in all cases.As will be shown later, this will cause an excessive cycling in the ESS. Figure 11 shows the example of applying this strategy during a clear day.In the absence of methods of predicting fluctuations, the MA must be applied since sunrise, preparing for a potential fluctuation that does not happen.As a result, the ESS suffered one charge-discharge cycle unnecessarily (Figure 11c).Up to now, the T value required to limit the fluctuations to below the rmax limit is not known.This paper has provided an answer to this question.To do so, this strategy has been simulated for all the 5 s data over one year (2012) and for all the PV plant sections (0.55-38.5 MW).For a given value of rmax successive iterations were made, increasing the duration of T (with a 60 s interval) until the condition |ΔPΔt(t)| < rmax was met. Figure 12 shows the relationship between the resulting time window T and the maximum allowable ramp rmax, with the PV plant power output as a parameter.As can be seen, time window T is solely dependent on rmax, and not on the size of the PV plant, P*.That is the reason why in Figure 12

Step-Rate Control Strategy
Logically, the specifications required for the ESS are the direct consequence of the standards and regulations to be applied.For example, the reference [1] requires the smoothing of fluctuations of a very short duration, exceeding 10%/min.Translated to another time scale, variations of 100% in 10 min are permitted, being the typical fluctuations of a PV system with no energy storage provisions (for ranges of 1-40 MW) as seen in [15].For this reason, based on the conclusions of this study, very small storage times are required, of around 6 min.For other scenarios, such as [2], storage becomes more important: for a maximum ramp of 2%/min (in other words 20% in 10 min) the required storage time is increased by up to 45 min [15].For these more severe constraints, it makes sense to try and seek strategies that optimize the ESS system.We shall now go on to propose an innovative step-rate control strategy based on strict compliance with the maximum ramp constraint rmax for the defined time window (for example 10 min).Hence, the proposed strategy attenuates the fluctuations in that particular time window and higher.Below that magnitude (high frequencies) we will take advantage of the geographical dispersion of a group of PV plants which strongly smoothens the fastest fluctuations.Previous studies [8,31,32] have well proved that N PV plants dispersed and separated at least by a few kilometers (6 km is enough) smoothens out the high frequencies (below 10 min) proportionally with √N.We understand N to be sufficiently large to produce the necessary smoothing effect in order to mitigate the effect of these steps.Therefore, the step-control strategy is coherent in a scenario with a number N of dispersed multi-MW PV-ESS plants, in which the step-control is implemented in each one, in the knowledge the TSO will see a reduction in the step effect simply due to geographical dispersion.Therefore, step control makes sense provided that the TSO permits strict compliance with rmax in the given time window.A detailed example of how this strategy works is shown in Figure 13, and it will be explained below.The algorithm of this strategy is as follows (Figure 14): the first step is to decide whether the present evolution of the PV generation is positive or negative.Then, a check is made to determine whether the ramp condition is met in a previous time window with a duration of n times the sampling time Δt (5 s in our case), in other words.This entails strict compliance with the ramp condition for times that are equal to or higher than n•Δt, but not below this value.For example, for grid regulations evaluating the fluctuations every 10 min, this would mean that the fluctuations would need to be below rmax = 20%/10 min.For this case, it would be necessary to choose n•Δt = 10 min (n = 120 as Δt is equal to 5 s). Figure 13 shows the ramp control response PG,ramp and that of the step-rate PG,step for a negative fluctuation recorded on 15 February 2012 at 14:25:10 am.Before the significant negative fluctuation (during the first instant in Figure 13), short and small fluctuations (high-frequency) can be seen.These fluctuations are lower than rmax (2%/min) evaluated in 10 min (20%), so step-control does not act.However, when a negative fluctuation of around 90% takes place in PPV, step-control acts smoothing this fluctuation in PG up to a 20% value.From an energetic point of view, thanks to the ESS, the power injected by the PV plant PG,step evolved from 100% to 10% in 40 min, similarly to ramp-rate control, PG,ramp.The only difference between PG,ramp and PG,step are the steps, fast fluctuations below 10 min.This high frequency will be strongly smoothed by geographical dispersion as it has been well proved in [8,31,32], and then, going unnoticed for the TSO.
The area within PG,ramp and PG,step corresponds to the necessary battery reduction.Logically, the higher the n•Δt value, the less battery required.To prevent the constant battery discharge, this strategy also requires the same SOC control proposed for the ramp strategy (Figure 5).
Once more, determining the capacity needed by this strategy can be easy solved thanks to the worst fluctuation model (Figure 15).The energy saving for a power plant P* using the step-rate strategy instead of the ramp-rate, Estep-ramp, (kW• h) is the integral of the difference between PG,ramp(t) and PG,step(t), or, in other words, the total area of each of the triangles Ntri of Figure 15, in other words, Equation ( 11 where P* is expressed in (kW), and n•Δt in (s).For example, for P* = 1.1 MW (τ = 6.14 s) and rmax = 20%/10 min, Equation (3) gives EBAT,MAX = 373 kW• h.From Equation (11), with n•Δt = 600 s the saving Estep-ramp is equal to 82.5 kW• h, that is 22%.As can be observed, this percentage is independent of rmax and the PV plant size (or, in other words, τ).By way of example, Figure 16 again shows the results of simulating this strategy for 15 February and with rmax = 20%/10 min (2%/min).For this strategy and according to Equation (11), CBAT needed to smooth the worst fluctuation is equal to 582 kW• h, so EBAT,ref = 291 kW• h.For that given day, the required battery power is PBAT,MAX = 879 kW (or PBAT,MAX = 0.79• P*) and the required effective battery capacity is CBAT = EBAT,MAX − EBAT,MIN = 312 kW• h (or 17 min of storage equivalent to 0.28 h).The total energy that has passed through the battery is 419 kW• h, equivalent to 32% of daily production, whilst the total losses in the battery and PEC are 2.3%.Obviously, PG is more fluctuating than in the case of step-rate control, but the stress in the ESS and capacity needed are lower.The steps produced in PG by this strategy can be better seen in Figure 13, which is precisely an enlarged version of Figure 16a.

A Comparison of the Smoothing Strategies
These three strategies were simulated for 5 s data for one year (2012) and for all the Amaraleja plant sections (from 0.55 MW to 38.5 MW), with a ramp restriction of rmax = 2%/min (20%/10 min for the step-control), a similar restriction imposed by [2].For the average moving strategy, this restriction requires a minimum T of 2700 s, according to Equation (10).For the step-rate control, n = 120 (120 × 5 s = 10 min).It should be pointed out that the SOC was calculated on the effective storage capacity CBAT, thereby allowing 100% variations in the SOC (effective SOC).As a result of this, in the case of the ramp control and the step rate control, the reference SOC is 50% whilst, for the moving average, the battery starts and ends each day fully discharged (SOC = 0%).The indices of merit selected to compare the strategies are the effective storage time tbat, the storage system losses, the degradation due to cycling and the quality of the wave injected into the grid.

Effective Storage Time, tbat
The difference in the storage capacity required for the different strategies depends on the following factors: rated power output of the PV plant P* and maximum allowable ramp rmax.The first factor affects the dynamics of the power generated PPV, the greater the size of the PV plant, the greater the filtering in the power fluctuations [36].Figure 17 shows the storage time required, tBAT, for each strategy, based on the size of the PV plant.This parameter was calculated after annual simulation, as the difference between the maximum and minimum annual energy value in the battery, EBAT,MAX − EBAT,MIN divided between the rated power output P* of the PV system.Likewise, solid lines in Figure 17 shows the theoretical storage time calculated by Equations ( 4), ( 9) and (11) for each strategy.4), ( 9) and (11).
The step-rate control requires an average 20% less storage than the ramp-rate control (the greater the value of parameter n•Δt, the less battery required).Clearly, the MA strategy requires the smallest battery capacity.However, this strategy does not benefit from the fluctuation smoothing due to the size of the PV plant, as shown in Equation ( 3), unlike the ramp-rate and step-rate control strategies [15].

RAMP-RATE MOVING-AVERAGE STEP-CONTROL
The best way to demonstrate this is Figure 18.It shows the relationship CBAT,ramp/CBAT,MA, constructed from Equations ( 4) and ( 9) for different values of rmax and P*.As rmax and P* increases, CBAT,ramp/CBAT,MA decreases.In any case, for restrictions around 2%/min, approximately half an hour of storage is sufficient to smooth the fluctuations.

Losses in the Storage System
Figure 19a shows the total energy loss in the ESS (PEC and battery) for rmax = 2%/min, expressed as a fraction of the annual production, based on the PV plant size.As can be seen, the losses for the moving average strategy are considerably higher than those for the other two strategies (2-3 times more).The reason lies in the fact that this strategy must act every day, regardless of whether there are any fluctuations or not, whilst the rate control strategies only act in the event of fluctuations greater than rmax.Figure 19b shows again the losses but for a less stringent ramp, rmax = 10%/min, as grid code in [1].Moving average losses are considerably higher than for ramp-rate control.Note that in this new perspective, step-control would not proceed, because it allows fluctuations of 100%/10 min.In any case, it should be highlighted that the losses are small for all strategies (<1.4% in the worst case).

Stress in the Storage System
Figure 20 shows the histogram for the annual effective SOC of the storage system for each of the strategies simulated in Section B (1.1 MW).As can be seen, the SOC distribution for the ramp-rate and step-rate strategies is very similar, being derived from the same algorithm: the control guarantees at all times that the effective SOC is maintained at around 50%, only moving away from this reference value when significant fluctuations occur (charges and discharges to 50%).These situations are extremely rare, consistent with the annual distribution of fluctuations already seen in [14].However, the SOC distribution for the MA is clearly bimodal, and is quite different from the rate-control strategies.Each day the battery undergoes a complete charge-discharge cycle (100% charge and discharge) with the maximum battery energy achieved at around midday, with full discharge at the end of the day, remaining in this state for the entire night.This explains the frequent SOC values close to 0% and 100%.However, Figure 20 is insufficient to determine the extent of the storage system degradation due to cycling.This is largely dependent on the number of Ncicl cycles and the depth of discharge (DOD).However, considering an annual effective SOC profile, such as those obtained in our simulations, there are methods to determine both parameters.One of the most-used algorithms and which provides the best results is Rainflow-counting [37].Initially developed to calculate mechanical fatigue [36], this method has also proved to be equally valid for calculating the ESS ageing due to cycling [38][39][40][41][42].In this way, any complex charging and discharging series can be broken down into a series of simple sub-cycles with a given DOD.The methodology followed in this document is similar to that shown in [39].
Figure 21 shows the result of applying the Rainflow-counting algorithm to the annual effective SOC for the various strategies and Section B (1.1 MW), for rmax = 2%/min and 10%/min.The figure shows the number of cycles Ncicl occurring (y axis) and the DOD (x axis).For the moving average strategy the ESS performs more cycles and with depths of discharge which are practically double those of the other strategies.(a) (b) Once these curves have been obtained, the calculation of the annual degradation due to cycling is immediate.Assuming again that our ESS is a lithium ion battery, with a life cycle curve similar to that shown in Figure 22 [34].This shows the number of maximum cycles Nmax for a given DOD.It is possible to calculate the annual degradation due to cycling CBAT,loss as the sum of each individual degradation to a given DOD, in other words the relationship between the number of annual cycles observed Ncycl, and the number of maximum possible cycles Ncycl,max, is given in Equation ( 12): The results obtained by Equation (12) show the degradation solely due to the effect of the ESS cycling.However, additional effects exist which are not going to be taken into account and which affect service life, such as operating temperature.The method applied here for lithium batteries is also valid for any other storage technology, simply by knowing the ESS life cycle curve (Figure 23).The results of applying Equation (12) to the values of Figure 21a are equal to 1%, 11.09% and 0.6% for ramp-rate, MA and step-rate, respectively.Given the fact that the manufacturer does not permit DOD values greater than 80%, the data in Figure 21a have been re-scaled to correspond to a DOD of between 0% and 80%.
This calculation methodology can be repeated for all the sections of the Amareleja plant observing the cycling degradation based on P* and the strategy considered (Figure 23a).The ramp-rate and step-rate control strategies show a very low cycling degradation (0.5%-2%) which hardly affects the useful life of the battery, unlike the moving average strategy with a degradation of around 10 times greater.As can be seen, the step-control cycling degradation is in any case lower than that of the ramp-rate.Once more, as it has been done with the losses analysis, cycling degradation has been repeated for a less stringent restriction, rmax = 10%/min (Figure 23b).In this situation, capacity requirements CBAT, are lower, therefore cycling degradation is bigger than for 2%/min.Anyway, the difference between the strategies is similar to the ones showed for rmax = 2%/min.

RAMP-RATE MOVING-AVERAGE
It is remarkable that CBAT,loss for the MA strategy is around 10 times bigger than ramp-rate control.It could be argued that this comparison is not fair, because necessary CBAT for the MA is, roughly speaking, the half respect ramp-rate control.Hence, CBAT,loss has been recalculated considering the same capacity CBAT for all the strategies and equal to the ramp-rate control capacity (Figure 24).Now, CBAT,loss for the MA considerable decreases, but it still presents a degradation close to four times more than ramp-rate control (3.6% versus 1% for rmax = 2%/min and 1.1 kW).This difference grows closer to 5-7 times with respect to the step-control.Recently, in order to avoid the daily MA cycling, Perez et al. [43] proposed to base the MA on the daily clear sky output ratio.Even so, it needs to resolve practical points such as control of the SOC (for example, how a midday negative fluctuation would be covered if the battery is discharged) or how to avoid constant discharge due to losses in the conversion elements.Likewise, it points out how forecasts can, to a large extent, help anticipate fluctuations and optimize the storage required and, in fact, the reference [43] discussed the improvement obtained with a perfect forecast, although how to face a situation in which the forecast fails has not been properly resolved.

Quality of the Signal Injected into the Grid
Figure 25 shows PG frequency spectrum for all the data for one year for Section B (1.1 MW) and for each of the strategies proposed in this article.The methodology followed is identical to that presented in [44,45].Basically the method consists in applying the discrete Fourier transform (DFT) to variable PG over an entire year, previously normalized (using P*) computed by a fast Fourier transform (FFT) algorithm.As can be observed, the three strategies analyzed offer a similar cut-off filtering frequency of around 0.2 mHz (1.4 h), smoothing higher dynamics, particularly the MA strategy.However, it can be observed how the effect of the step-rate control filtering is reduced after frequencies close to 10 min (1.2 mHz), allowing faster dynamics to pass through and, therefore, injecting poorer quality power into the grid.However, this faster dynamics are easily reduced when we take in account the production of a PV fleet geographically dispersed, as it was demonstrated in [5,31,32].

Conclusions
This paper has quantified the storage requirements needed to smooth short-term PV power fluctuations based on different control strategies and for a wide range of PV system power outputs.The one-year simulation of the strategies has demonstrated that the MA strategy requires the least storage capacity.However, it would be necessary to analyze in each particular case whether, firstly, this saving in size offsets the considerable increase in 2-3 times more of losses.Secondly, whether the storage technology is compatible with greater cycling (about 10 times more) and with SOC values close to zero for long periods of time.For example, for the case of lithium ion batteries, this strategy would have an extremely negative impact on the battery useful life, although this impact may be lower for ultra-capacitors or flow batteries.The new proposed strategy, namely the step-rate control, provides the most efficient use and with less ESS cycling in situations with several PV plants geographically dispersed and exigent ramp restrictions.Under less restrictive conditions, the ramp-rate control would be preferable against MA.Likewise, for a given grid restriction, an empirical expression is put forward to make it possible to determine the energy capacity required for each strategy.
We would highlight the lower losses achieved for the proposed system, particularly for the rate-control strategies (ageing due to cycling of around 1%/year-2%/year).This result draws attention to the potential use of the ESS for ancillary services (frequency regulation or time shifting), in an endeavor to maximize the value of its installation.To do so, a fluctuation prediction tool would be necessary, and this would need to be accurate enough to determine the critical times, in order to ensure the smart management of the energy stored.

Figure 1 .
Figure 1.Field distribution of the Amaraleja photovoltaic (PV) plant sections considered in this paper.

Figure 5 .
Figure 5. Ramp-rate control model modified with additional state of charge (SOC) control.Note that the SOC control action is also smoothed by the ramp-limiter in order to guarantee that power fluctuations are always below rmax.

Figure 6 .
Figure 6.Worst fluctuation model.The blue line represents the PPV(t) response to an irradiance fluctuation (yellow line) and the red one is the power injected to the grid PG with a ramp-rate control.The difference between PG and PPV is PBAT, the maximum difference corresponds to PBAT,MAX and the defined integral of PBAT corresponds to EBAT,MAX.

Figure 7 .
Figure 7. (a) Evolution of the generated power, PPV(t) on 15 February 2012 by Section B and the simulated power that would have been injected to the grid PG(t) if a battery had been available to limit fluctuations to rmax of 2%/min (0.166%/5s); (b) battery power, PBAT; and (c) battery energy, EBAT.The simulation has been made based on the model shown in Figure 5.

Figure 8 .
Figure 8. Model of the moving-average (MA) strategy, offsetting the ESS losses.

Figure 9 .
Figure 9. Model to calculate the storage capacity needed CBAT,MA for the MA strategy.

Figure 10 .
Figure 10.(a) Evolution of the generated power, PPV(t) on 15 February 2012 by Section B (1.1 MW) and the simulated power which would be injected to the grid PG(t) under the moving average strategy (T = 2700 s); (b) battery power, PBAT; and (c) battery energy, EBAT.The simulation has been made based on the model shown in Figure 8.

Figure 11 .
Figure 11.(a) Evolution of the generated power, PPV(t) during a clear sky day, 10 August 2011 by Section B (1.1 MW) and the simulated power which would be injected to the grid PG(t) under the moving average strategy (T = 2700 s); (b) battery power, PBAT; and (c) battery energy, EBAT.
all the curves are superimposed.The evolution of these values, calls for making a fit based on function T = m• rmax −1 , resulting in Equation (10): rmax is given in (%/min), and m = 5400 (s• %/min).The goodness of fit (R 2 = 0.99) confirms the validity of the Equation(10).

Figure 12 .
Figure 12.Relationship between the size of the power plant P*, the maximum allowable ramp, rmax and window T, Equation (10).

Figure 13 .
Figure 13.Evolution of the generated power, PPV(t) on 15 February 2012 by Section B (1.1 MW), the simulated power which would be injected to the grid PG(t) in the case of the availability of a battery which limits fluctuations to rmax of 20%/10min using the step-control strategy (n = 120, n•Δt = 600 s) and the ramp rate control.The reduction of the storage required compared to the ramp control is evident and corresponds to the area between PG,ramp and PG,step.

Figure 14 .
Figure 14.Flowchart of the step-rate control.

Figure 15 .
Figure 15.Model to calculate the storage capacity required for the worst fluctuation and step-rate strategy via EBAT,ramp.

Figure 16 .
Figure 16.(a) Evolution of the generated power, PPV(t) on 15 February 2012 and the simulated power which would be injected to the grid PG(t) under the step-rate strategy (n = 120, n•Δt = 600 s), limiting fluctuations to rmax of 20%/10 min; (b) battery power, PBAT; and (c) battery energy, EBAT.The simulation has been made based on the model shown in Figure 5.

Figure 17 .
Figure 17.Storage time required, tBAT, for each strategy, based on the size of the PV plant with a ramp restriction of rmax = 2%/min.Solid lines show the theoretical storage time calculated by Equations (4), (9) and(11).

Figure 18 .
Figure 18.Relation CBAT,ramp/CBAT,MA based on the PV plant size P* and the maximum permitted ramp rmax.

Figure 19 .
Figure 19.Losses at the PEC and battery, based on the size of the PV plant considered and the fluctuation smoothing strategy, with a ramp restriction of: (a) rmax = 2%/min; and (b) rmax = 10%/min.

Figure 20 .
Figure 20.Histogram for the effective SOC in the storage system resulting from the simulation of the ramp-rate, MA and step-rate strategies for rmax = 2%/min.Note that the y-axis is logarithmically scaled.

Figure 21 .
Figure 21.The number of observed cycles at a given depth of discharge (DOD) for the strategies presented in this paper, for a maximum ramp limitation of: (a) rmax = 2%/min in Section B (1.1 MW); and (b) rmax = 10%/min.Note that the y axis is logarithmically scaled.

Figure 23 .
Figure 23.Annual degradation of the ESS due to cycling for: (a) rmax = 2%/min; and (b) rmax = 10%/min for all strategies and based on the PV plant size.

Figure 24 .
Figure 24.Annual degradation of the ESS due to cycling for: (a) rmax = 2%/min; and (b) rmax = 10%/min for all strategies and based on the PV plant size assuming that all the strategies work with the CBAT required for ramp-rate control.

Figure 25 .
Figure 25.Spectrum for annual PG resulting from the simulation of the three control strategies for the data for PPV of Section B (P* = 1.1 MW).