Economic Performance Optimization of a PV-BESS Power Generator : A Case Study La Reunion Island

This paper proposes an economic performance optimization strategy for a PV plant coupled with a battery energy storage system. The case study of La Reunion Island, a non-interconnected zone (NIZ) with a high level of renewable energy sources (RES), is considered. This last decade, to reach the ambitious target of electricity autonomy by 2030 set by the local authorities, local and national plans have been launched to promote renewable energy sources integration that led to a noticeable development of photovoltaic (PV) systems. To avoid a decrease of the grid reliability due to a large integration of intermittent energy sources into a non-interconnected grid, the authorities have introduced new regulatory rules for RES producers. The proposed optimization strategy relies on these new regulatory rules and takes into account the energy market data, the amount of PV production subject to penalties for imbalance, the batteries and the PV technological characteristics together with a PV production forecast model. Due to its high convergence rate to the true global minimum and its perfect suitability to practical engineering optimization problems, the recently developed Modified Cuckoo Search algorithm is used as optimization algorithm. The effectiveness and relevance of the proposed strategy are assessed on experimental data collected on a real PV power plant. An economical analysis demonstrates that the proposed optimization strategy is able to fulfill the new regulatory rules requirements while increasing the economic performance of the system.


Introduction
To date, reducing carbon emission has become a major concern.Among possible options, increasing shares of renewable energy sources (RES) such as solar, wind or biomass resources appear as a promising solution for a cleaner power generation.High shares of RES may then become a critical aspect of future energy systems.In this context, small islands that mainly rely on imported fossil fuels for energy production are likely to be pioneers in the development of decarbonized electricity production [1] [2].
In this study, the case of La Reunion Island, a non-interconnected zone (NIZ), is considered.Even if the territory has a high level of RES, its electricity production remains strongly based on imported fuels.In this context, local authorities have set the ambitious objective of reaching electricity autonomy by 2030.This last decade, to reach this target, local and national plans have been launched to promote RES integration [3] [4].Thus, supported by incentive mechanisms such as tax exemptions, direct subsidies or feed-in tariffs, photovoltaic (PV) systems have experienced a rapid and noticeable development [5].However, a large integration of intermittent sources into a non-interconnected grid raises critical technical issues due to the uncertainties of the energy production.The intermittency and unreliability of solar-generated power may reduce the network stability and lead to load shedding or to the interruption of electric service [6].To avoid such situations, the authorities have set a limit of 30% of intermittent sources in the instantaneous electricity production and have introduced new regulatory rules for RES development.Henceforth, RES producers have to declare to the grid operator, a day in advance, the power profile that will be injected to the grid.Then, if the power plants do not meet the submitted schedule for injected power, they face financial penalties.In La Reunion Island, if mismatches between actual and scheduled power injection exceed a given tolerance, RES producers are charged with imbalance penalties.In order to address the problem related to the intermittency of solar-generated electricity while reducing the amount of PV production subject to penalties for imbalance, energy storage systems (ESSs) appears as one of the most relevant option.Recently, several works related to the applicability, advantages and disadvantages of various ESS technologies for RES integration have been reported [7].As regards PV power plants coupling with ESS, several works dealing with technical issues and economic feasibility have been conducted [8] [9] [10] [11].However, from a regulatory point of view (incentive schemes and economic feasibility), nearly all works reported in the literature focus on the determination of the optimal sizing of the ESS [12] [13] [14] [15].
In the case of La Reunion, and according to our best knowledge, none work has been conducted to optimize the economical performance of existing hybrid photovoltaic-battery energy storage system (BESS) power generators, based on the latest regulatory rules.In this paper, an economical optimization of a hybrid PV-BESS power generator is developed.The proposed methodology relies on a metaheuristic optimization algorithm taking into account the energy market da-ta, the amount of PV-generated energy subject to penalties for imbalance, the PV and the batteries technological characteristics together with a PV production forecast model.To assess the effectiveness and relevance of the proposed strategy, the economic analyses are performed on data measured on a real power plant.Indeed, a one-year experimental data, collected from August 2013 to August 2014 on a 57 kWp PV farm coupled with a 78.5 kWh BESS, are considered.
The rest of this paper is organized as follows.The regulatory rules applied in La Reunion are presented in Section 2. Section 3 is dedicated to model design.In this section the PV production forecast model and the energy storage model are detailed.The performance of the proposed optimization strategy in terms of economical efficiency improvement is demonstrated in Section 4.

Regulatory Rules
In non-interconnected zones (NIZ) such as La Reunion, the large integration of intermittent sources raises critical technical issues issue related to the reliability of power supply.The reliability of an electrical grid can be defined by its ability to supply the aggregate electrical demand and energy requirements of the customers at all times, while withstanding sudden disturbances such as unanticipated loss of system elements (e.g.load or production fluctuations) [16] [17] [18].Currently the reliability of power supply in La Reunion is already lower than in Metropolitan France.In La Reunion the average duration of electricity not supplied is estimated at 4 h/year/consumer compared to 73 min in Metropolitan France [5].
In this context, and considering the rapid and important growth of PV systems in La Reunion the last decade, the authorities have recently decided to set up new regulatory rules to ensure the reliability of the power supply.Henceforth, producers have to declare a time generation profile that represents the day-ahead power injected by their plants, at a minute basis.If the mismatches between the actual injected power and the announced power exceed the admitted tolerance, financial penalties are applied.According to this regulatory framework, energy imbalance is calculated with minutely resolution, and the tolerance band is taken equal to ±5% of the installed PV power capacity ( ) The electricity tariff system relies on peak and off-peak hours.During peak hours, 7 PM to 9 PM, the electricity feed-in tariff is more attractive.However, during this time period, producers have to guarantee a constant power injection to the grid comprised between 20% and 70% of PV_peak P .The current electricity tariff system applied to PV power producers in La Reunion, including peak and off-peak feed-in tariffs, is summarized in Table 1.
Table 1.Summary of tariff system in La Reunion.Note that producers have the possibility to buy electricity from the grid.In some very specific cases, it could be interesting to buy electricity from the grid during off-peak, store the energy in an ESS, and sell it back during peak hours.
The producer's revenue is calculated each minute using the follow expression: where inj P denotes the power injected to the grid, out P the power extracted from the grid, and the financial penalties are calculated as follows: where bid P denotes the day-ahead schedule power profile to be injected to the grid.
Every time inj P is outside the tolerance band, the system is said to be in faulty condition and financial penalties are applied.
In this work, the economic effectiveness of the system is assessed using two criteria, which are the revenue and the daily fault rate (DFR).This last criterion defines the ratio of time where the system is in faulty condition each day: With fault τ the cumulated time of fault condition in minute.
To fulfill the requirements of the new regulatory rule, PV power producers have to announce a day in advance the power profile to be injected to the grid, which requires a PV production forecast model.Besides, regardless of the accuracy of the PV production forecast model, financial penalties due to imbalance are unavoidable.Therefore, to reduce financial penalties and above all take advantage of peak hours feed-in tariff, the use of ESS appears to be a relevant option.

PV Production Forecast Model
In the literature, a wide variety of parametric and non-parametric forecast models have been reported [19].Parametric models require a wide set of information about the PV power plant technology and its installation configuration.Nonparametric models are generally based on weather forecast models [20].These limitations make the reliability and the suitability for "on field" uses of parametric and non-parametric forecast models questionable (Table 2).
In this study, regarding practical purposes, the widely used persistence model is chosen to forecast the PV output power at a minute basis.This is a simple model based on the assumption that the PV production of today is the same as yesterday [21]: where ( ) PV _ m P t denotes respectively the measured and forecast PV output power at time t .
Even if this method does not take into account the intra-day variability of solar irradiance, it represents with a good accuracy the periodicity and seasonality of weather conditions (day/night and summer/winter cycles) [13].Besides, it has the merit of introducing no incremental cost, which cannot be underestimated for small cases applications.Obviously, the accuracy of the PV production forecast could be improved using more sophisticated models but at the price of increasing complexity and computational cost.

Energy Storage Model
Regarding optimization purposes and according to the considered time scale Subjected to constraints on power and capacity: ( ) ( )

Economic Performance Optimization
In this work, a minute dispatch strategy for a 57 kWp PV farm with 78.5 kWh BEES is implemented and an economical optimization of the dispatch strategy is proposed.BESS has two main applications: first, compensate PV production forecast errors during off-peak and thus reduce financial penalty due to imbalance.Second, inject power to the grid during peak hours and thus take advantage of the attractive feed-in tariff.
In this context, two parameters are introduced.The first one, denoted peak α , is related to the amount of energy dedicated to the evening peak hours ( ) , which have to be stored in the BESS meanwhile, and is written as a fraction of the estimated total energy produced by the PV plant ( ) Here . The second one, denoted forecast β , represents the fraction of the storage capacity that is dedicated to compensate power imbalance ( ) foreacast C due to forecast errors: To assess the effectiveness and relevance of the proposed strategy, the economic analyses are performed on a one-year experimental data, collected from August 2013 to August 2014 at La Reunion on a real PV power plant.Due to its high convergence rate to the true global minimum and its perfect suitability to practical engineering optimization problems, the recently developed Modified Cuckoo Search algorithm proposed by [22] is used as optimization algorithm.

Modified Cuckoo Search Algorithm
Modified Cuckoo Search (MCS) algorithm can be regarded as a modification of the standard Cuckoo Search (CS) proposed by Yang and Deb [23] to increase its convergence rate and makes it more suitable for practical engineering optimization problems.
CS algorithm relies on the observation of the reproduction strategy of cuckoos.
Cuckoos lay their eggs in the nests of other host birds that may be of different species.If the host bird discovers that the eggs are not its own it could either destroy the egg or abandon the nest [24].To apply this as an optimization tool, Yang and Deb [23] used three idealized rules: • Each cuckoo lays one egg (a set of solution) at a time and dumps it in a random nest; • The best nests, those containing the highest quality of eggs (solutions), will carry over to the next generation; • The number of nests is fixed, and there is a probability The main modification introduced by [22] consists in adding information exchange between the top eggs or the best solutions to increase the convergence rate of the method and extend its applicability to a wider range of applications.
Standard optimization benchmarking functions have been used to test the effects of this modification and it has been demonstrated that, in most cases, the MCS performed as well as, or better than, the standard CS, a particle swarm optimizer, and a differential evolution strategy.MCS and CS algorithms and their performance compared to other standard optimization strategies are fully detailed in [22].

PV/BESS Control Rules
The power injected to the grid is defined as the sum of the PV output power and the storage power exchanged with the AC bus: The BESS is used to adjust the PV power plant PV _ m P to maintain, as much as possible, the difference between the day-ahead announcement bid P and the actual power injected inj P to the grid within the tolerance band.If PV _ m P is above (respectively below) the tolerance limit, the BESS can be used, when it is possible, to store (respectively deliver) the imbalance power I P .Depending on whether the measured output PV is within, below or above the tolerance band, I P is defined as follows: Every time 0 I P ≠ financial penalties for imbalance are applied.Therefore, the storage charge/discharge process must be suitably controlled in order to reduce the DFR and thereby the amount of financial penalties, while ensuring that there is enough energy stored in the BESS for the evening peak.In this aim, a tolerance band control strategy is proposed and a specific control rules is designed for each zone: above, within and below the band.

Above the upper limit
When PV _ m P is above the tolerance band, the BESS is used to compensate the imbalance power and store the excess of energy whenever possible.Indeed, the imbalance power cannot always be compensated.Several conditions related to operational and technical limits of the BESS have to be verified (i.e state of charge, charge/discharge rate limits).In this case, the storage power is computed as follows: ( )

Below the lower limit
When PV _ m P is below the tolerance band, the imbalance power due to forecast error can be compensated using the energy stored in the BESS.However, this energy has to be manipulated very wisely in order to ensure that enough energy remains to guarantee peak hours.The storage power is calculated according to operational and technical limits of the BESS and subjected to constraints on _ MAX _ peak max max foreacast min , Within the tolerance band When PV _ m P is within the tolerance band, the power imbalance is equal to zero and there is not need to absorb/inject power from/to the grid through the BESS:

Economic Performance Improvement
The economical performance improvement relies on the estimation of two parameters peak α and forecast

β
. The effects of these parameters on the annual revenue and the average annual DFR are illustrated on Figure 1 and Figure 2.
As expected, while forecast β is increasing the DFR is decreasing.Indeed, forecast β is straightforwardly linked to the fraction of the storage capacity dedicated to compensate power imbalance due to forecast errors.However, it is important to highlight that decreasing the DFR and so the financial penalties do not necessary means increasing the revenue.In fact, while forecast β is increasing, the amount of energy that can be stored in the BESS for the evening peak hours decreases since foreacast C becomes closer to max C (Cf. Equation ( 8)).Indeed, Figure 1 illustrates that it could be more interesting to sell more energy during peak hours even if that means paying more penalties during off-peak due to forecast error.Which can be explained by the fact that peak hours feed-in tariff is more attractive than off-peak one.As illustrated on Figure 2, while peak α is increasing the DFR is increasing, whereas the revenue increasing to a maximum before decreasing.Which seems indicate that, for a fixed value of forecast β , there is an optimal amount of energy to store in the BESS for the evening peak hours.

Annual Optimization
In a first attempt, the economic performance improvement strategy consists on finding the optimal values of peak α and forecast β that maximizes the annual revenue.The optimization goal is to find the optimal set of parameter The optimization procedure, performed on experimental data collected from  The analysis of the forecast error reveals the same periodicity of 366 days, which means that the forecast error is linked to the season.Indeed, as illustrated on Figure 4, the forecast error is higher in summer than in winter.This can be mainly explained by the fact that the forecast model does not take into account the intra-day variability of solar irradiance and that in La Reunion Island the intra-day variability is higher in summer.In this context, when the intra-day variability increases the error modeling increases too, which leads to an increase of the DFR.In other worlds, the seasonality of the solar irradiance variability introduces a seasonal component into the forecast error that is transmitted to the DFR.
The analysis of the total energy produced by the PV plant each day ( ) PV _ _ total f E reveals a seasonal component that is due to the yearly solar irradiance variability (Cf. Figure 5).The amount of energy dedicated to the peak hours ( )

_ peak sto E
, which could have a strong influence on the revenue, is taken as a fraction of PV _ _ total f E using peak α .In this context, the seasonal compo-  .A thorough study of intra-day and yearly irradiance variability in La Reunion can be consulted in [25].
Regarding the seasonality of the solar irradiance, and since forecast β and peak α respectively influence the DFR and the amount of energy dedicated to the peak hours, it is likely that a seasonal-based or even a daily-based optimization of these parameters could increase the revenue.

Conclusions and Prospects
In this work, a minute dispatch strategy for a 57 kWp PV farm with 78.nomic performance of the system.The BESS is used during off-peak to compensate PV production forecast error and during peak hours to inject power to the grid.Therefore, two parameters have been introduced.The first one is related to amount of energy to store for the evening peak hours whereas the second one represents the fraction of the storage capacity that is dedicated to compensate power imbalance due to forecast errors.The optimization goal is to find the optimal value of these parameters that maximizing the revenue while taking into account the new regulatory rules constraints.The proposed optimization strategy takes into account the energy market data, the amount of PV production subject to penalties for imbalance, the batteries and the PV technological characteristics together with a PV production forecast model.The effectiveness and relevance of the proposed strategy have been assessed on experimental data collected on a real PV power plant.An economical analysis demonstrated that the proposed optimization strategy has been able to fulfill the new regulatory rules requirements while increasing the economic performance of the system.
Due to the seasonal behavior of solar radiation, it is likely that economical performance can be further increased using a seasonal-based optimization approach.Additional works are currently in progress to study if the revenue can be increased by taking into account the seasonal component contained in the forecast error and the total energy produced by the PV plant.
Peak hours (7 PM to 9 PM) [€ ct/kWh] Off-peak hours [€ ct/kWh] Selling price ( ) discovers an alien egg.If this happens, the host can either discard the egg or the nest, and builds a new nest in a new location (new random solutions at new locations).

Figure 2 .
Figure 2. Effects of peak α on the revenue and the DFR (with forecast 0.2 β = and calcu-

3 .
August 31 st 2013 to September 1 st 2014, leads to the optimal set of parameter to an annual revenue of 25 449.20 euros and an average DFR of 13.12%.The revenue and the DFR for each day are presented in FigureIt can be noticed that the DFR seems to contain a periodic component.The analysis of the DFR in the frequency domain reveals that the frequency component with the higher magnitude is located at 0.002732 day −1 , which corresponds to a periodicity of 366 days.Moreover, a thorough study of the DFR reveals a strong and significant correlation between DFR and the forecast error, with a Pearson's correlation coefficient of 82% (p-value < 0.0001).

Figure 3 .
Figure 3. Daily revenue and DFR obtain over one year.

Figure 5 .
Figure 5. Energy produced by the PV plant.
5 kWh BEES has been simulated, and an economical optimization of the dispatch strategy has been developed.This strategy has been designed to fulfill the requirements of the new regulatory rules set in La Reunion while optimizing the eco-C.Damour et al.
bid P Scheduled profile to be injected to the grid [kW] inj P Power injected to the grid [kW] peak α Parameter to be estimated forecast β Parameter to be estimated max C Maximum usable storage capacity [kWh] min SOC Min.energy storage level [% max C ] max SOC Max.energy storage level [% max C ] c η Efficiency of storage in charge d η Efficiency of storage in discharge s C Electricity selling price fault τ Cumulated time of faulty condition [minute]

Table 3 .
Storage system parameters values.
on the static characteristics of the battery (Cf.Table3) and neglects the transient dynamics of the process.In the sequel, powers are considered negative (respectively positive) during the discharge (respectively charge) phase.In this context, the battery state of charge (SOC) at time t is computed by: