Comparative and Cost Analysis of a Novel Predictive Power Ramp Rate Control Method: A Case Study in a PV Power Plant in Puerto Rico

: One of the most important aspects that need to be addressed to increase solar energy penetration is the power ramp-rate control. In weak grids such as the one found in Puerto Rico, it is important to smooth power ﬂuctuations caused by the intermittence of passing clouds. In this work, a novel power ramp-rate control strategy is proposed. Additionally, a comparison with some of the most common power ramp-rate control methods is performed using a proposed model and real solar radiation data from the Coto Laurel photovoltaic power plant located in Ponce, Puerto Rico. The proposed model was validated using one-year real data from Coto Laurel. The power ramp-rate control methods were compared in real-time simulations using the OP5700 from Opal-RT Technologies considering power ramp rate ﬂuctuations, power ramp-rate violations, ﬂuctuations in the state-of-charge, among other indicators. Moreover, the proposed power ramp-rate control strategy, called predictive dynamic smoothing was explained and compared. Results indicate that the predictive dynamic smoothing produced a considerably reduced Levelized Cost of Storage compared to other power ramp-rate control methods and provided a higher lifetime expectancy for lithium batteries. O.F.R.-M.; software, J.D.V.-P., Y.V.G. and O.F.R.-M.; validation, J.D.V.-P., Y.V.G. and O.F.R.-M.; formal analysis, J.F.P.-M.; resources, F.A.; writing—original draft preparation, J.F.P.-M., Y.V.G. and O.F.R.-M.; writing—review and editing, J.F.P.-M. and F.A.; supervision, F.A.; funding acquisition, F.A.


Introduction
Electrical energy in Puerto Rico has been mainly produced using imported fossil fuels. For the fiscal year of 2020, petroleum-based power plants supplied almost half of the total power generation, while coal-based power plants supplied 29%, and renewable energies only supplied 2.5% [1]. The Puerto Rico Electric Power Authority (PREPA) is aimed to increase renewable energy penetration to 40% by 2025, 60% by 2040, and 100% by 2050 [2]. Particularly, solar energy is the renewable resource that has grown the most with an increase in supply from 0.3% in 2015 to 1.4% in 2020 [3]. Additionally, PREPA is aimed to add up 2740 MW and 1440 MW in solar power and batteries, respectively in the following years [4]. Due to the inherent stochastic nature of solar irradiance, factors such as frequency and voltage stability of a converter-dominated power system (CDPS) may be critically affected on cloudy days. To overcome this issue, some local power authorities across the world have imposed solar power plants to include power ramp-rate control (PRRC) methods in their power injection controllers. A maximum power ramp-rate (PRR) of 10% of the contracted power capacity per minute is allowed in Germany, whereas Hawaii, Ireland, and Denmark allow a maximum PRR of 2 MW/min, 30 MW/min, and 100 kW/s respectively [5,6].
The Puerto Rico Energy Bureau (PREB) has developed a regulation on renewable energy development-including its registration and operation-as a strategy for rebuilding and strengthening the electric power system. The Bureau's purpose is to increase the local renewable energy generation and provide a regulatory framework for the implementation of renewable energy across the island. One of the most important requirements described in the regulation is the Power Purchase and Operation Agreement (PPOA) or electric power contract. The PPOA is a document that establishes an energy contract between two parties (the seller and the buyer). This contract contains mutual obligations between the parts involved in the power sale, that governs the rights and duties of the parties [7]. One of the most important requirements for third-party photovoltaic (PV) power plants stipulated in the PPOA is the Minimum Technical Requirements for the Interconnection of Photovoltaic Installations (MTR). MTRs are mandatory requirements that the third-party seller must meet to help the main grid in case of contingencies. In the MTR, the ramp rate requirement is addressed, which sets the maximum rate of change for the system's output power. This maximum rate applies for both, increasing and decreasing the output power with a requirement of 10% of the contracted capacity. Additionally, the current MTR requires the third-party systems to be able to deliver 30% of their measured AC capacity for at least 25 min using their storage units [8].
To fulfill PRR requirements and smooth power fluctuations, large-scale solar power plants (LSSPP) rely on the use of energy storage systems (ESS). The state-of-charge (SOC), the output power P out , and the depth-of-discharge (DoD) of the ESS through a day may vary drastically depending on the PRRC method and the type of battery used. These factors directly impact battery degradation and the levelized cost of storage (LCOS) [9][10][11]. Thus, LSSPP may incur extra costs in ESS depending on the implemented PRRC method. The ramp-saturation (RS) is the simplest PRRC method [12,13]. In this method, the battery absorbs or injects power when the PRR at the point of interconnection (POI) exceeds the requirement. The Simple Moving Average (SMA) and the Exponential Moving Average (EMA) are also two popular methods due to their low computational cost and ease of implementation [14][15][16][17]. However, the RS, MA, and EMA methods require a large ESS storage capacity [13,18,19]. This implies an increase in costs and battery degradation. PRRC methods such as the Enhanced Linear Exponential Smoothing (ELES) [20][21][22][23], the First Order Low-Pass filter (FOLPF) [24][25][26][27], and the second-order LPF (SOLPF) have been proposed to reduce power fluctuations in the ESS [28]. A solution to avoid the use of ESS is the Active Power Curtailment (APC) method [5,[29][30][31][32][33][34]. This method moves the power converter away from the maximum power point to reduce the power injection to the main grid. Though, the APC is only valid for ramp-up events. For this purpose, the Forecasted APC (FAPC) method has been proposed in [35,36]. This method determines whether a ramp-down event will occur and reduces the injected power by moving the power reference away from the maximum power point. However, as both APC and FAPC are aimed to avoid the use of ESS, their reliability is highly affected by the accuracy of the forecast.
As each of the aforementioned PRRC methods has its own advantages and disadvantages, it becomes of high interest to study how they affect the operational costs and how they affect the battery degradation in LSSPP. Furthermore, it becomes important to study the performance of each PRRC method in an LSSPP in Puerto Rico, which is characterized to have a special solar radiation behavior and different grid characteristics compared to other territories. Other studies about the impact of PRRC have been made in Puerto Rico [5,37,38]. However, none of them have tested the impact of different PRRC methods nor their impact on battery degradation. This work is intended to propose a predictive PRRC method and to perform a study comparing the aforementioned PRRC methods applied to a case study in a highly detailed model of the Coto Laurel solar power plant. Coto Laurel is a solar farm in Ponce, a town located on the southern coast of Puerto Rico. The installed capacity is approximately 14 MW DC, distributed in fifteen (15) groups of solar panels, for a total of 55,416 panels. The contracted capacity at the POI is 10 MW. Figure 1 shows the satellite view of the solar farm, each of the fifteen groups is depicted with a white box. Moreover, each group was labeled with a black capital letter to simplify the discussion. The green boxes depict the location of the inverters and their respective transformers, whereas the orange box shows the farm's battery bank and its respective inverters and transformers to connect to the main grid. The model of the Coto Laurel solar power plant was developed, tested, and validated using an OP5700 simulator from OPAL-RT technologies. Furthermore, the proposed predictive PRRC method was compared against others PRRC methods to analyze the effect in the LCOS and the battery degradation.
The rest of this document is organized as follows: In Section 2, a brief overview of the most common PRRC methods is made. Then, in Section 3, the proposed predictive PRRC method is presented. Section 4 presents the methodology for the economic analysis of each of the PRRC methods. Section 5 presents the case study in Coto Laurel, Puerto Rico. Section 6 presents daily and yearly simulation results and a comparative analysis Finally, discussions and highlighting are presented in Section 7.

Overview of the Power Ramp-Rate Control Methods
A simplified scheme of an LSSPP with PRRC, similar to Coto Laurel, is shown in Figure 2. The DC voltage of each PV array is transformed to AC using three-phase voltagesource converters (VSC). The output of each converter is connected to PREPA using a central feeder transformer. Similarly, the ESS is composed of a set of battery banks. The voltage of each battery bank is transformed to AC using a VSC. The ESS has its own transformer to inject or absorb power from the AC bus. The PRRC block reads the power in the PV feeder and generates a power reference to the VSC so that the power injected into the main grid fulfills the PRR requirements. The output power in the ESS feeder can be positive (battery charge) or negative (battery discharge). As seen in Figure 3, when the PRR in the PV violates the requirement, the battery absorbs or injects energy to smooth the output power to the main grid. Depending on the PRRC method, the battery may achieve a higher or lower DoD, and the SOC at the end of the day may vary as well. These two factors drastically affect the battery degradation and the cost associated with ESS for PRRC.
Following, a brief description of each of the PRRC methods analyzed in this work is presented. To simplify notation, the power injected into the main grid is represented by: Furthermore, the PRR of the power injected by the PV array at time t = kT s is represented by: where T s is the sampling period. Additionally, the maximum allowed PRR at the POI by regulation is defined as PRR max .

Ramp Saturation (RS) Method
The RS is one of the most common and easy-to-implement PRRC methods. This method evaluates the PRR of the PV and absorbs or delivers energy from the ESS to saturate the PRR at the output. The RS method can be expressed as [12,13]: As PRR max decreases, the ESS VSC will require more energy to smooth power transitions at the output. Thus, the DoD of the ESS will increase because the SOC will reach lower values.

Simple Moving Average (SMA) Method
The SMA method is known for its ease of implementation and because it allows for adjusting PRR fluctuations. The mathematical expression for the SMA method is given by [14,15]: where W is the window size for calculating the average. For this method, as the window size W becomes larger, the power fluctuations at the output become softer. However, this implies more DoD and more computational cost.

Exponential Moving Average (EMA) Method
The EMA method is an extension of the SMA. However, this method takes priority on recent samples by using an exponential weighting function as follows [16,17]: where the exponential smoothing factor α is between 0 and 1. As α becomes closer to 1, more priority is given to older samples. As this method gives more priority to recent samples, it is more suitable for smoothing power peaks at the output while maintaining the mean value.

First Order Low-Pass Filter (FOLPF) Method
The FOLPF is used to suppress high frequencies from the input signal with an attenuation factor of −20 dB/dec after the cutoff frequency of the oscillations. In the case of PRRC, the cutoff frequency should be selected to attenuate oscillations in the PV power. Thus, it is important to analyze the worst-case power oscillation that can be observed during the day. The output equation for the FOLPF is given by [24][25][26][27]: where and Z −1 {·} is the inverse z-transform operator. The terms T s and T f are the sampling period and the time constant, respectively. As the ratio T s /T f becomes larger, the bandwidth of the FOLPF increases. It is important to note that, if T s /T f become larger than 2, the eigenvalues of G folpf (z) become located outside the unitary circle and the FOLPF will become unstable. Thus, the desired transient response must be selected in terms of the sampling period.

Second Order Low-Pass Filter (SOLPF) Method
The SOLPF is similar to the FOLPF. However, the attenuation factor for the SOLPF is −40 db/DB. This implies that the SOLPF is more selective below the cutoff frequency, compared to the FOLPF. This implies that high-frequency fluctuations in the PV power can be more attenuated, which helps to reduce PRR violations. The discrete transfer function for the SOLPF is given by [28]: where ω n is the natural frequency and ζ is the damping ratio. As the eigenvalue location of G solpf (z) is not as explicit as G folpf (z), it is important to analyze them before implementation.

Enhanced Linear Exponential Smoothing (ELES)
The ELES method performs two complementary exponential smoothing operations. Although this method has high mathematical complexity, it is known to have a better performance compared to the aforementioned methods [19]. The mathematical expression for the ELES is given by [20][21][22][23]: where P pv and P pv are auxiliary variables for computing the double complementary filter. Constant α is the smoothing factor, which is required to be selected between 0 and 1.

Proposed Predictive Dynamic Smoothing PRRC Method
The proposed Predictive Dynamic Smoothing (PDS) method uses predicted PV output power values to adjust the PRR considering future power fluctuations. For example, if the prediction indicates that there will be a ramp-down violation, the PDS method will charge the battery prior to that event. Thus, the DoD and the battery degradation will be reduced. For this work, it is assumed that the predicted PV output power values are already acquired. The prediction algorithms presented in [6,35,36,39,40] are examples with suitable results for the implementation of the proposed PDS method and their development is out of the scope of this work. In the results section, the PDS method was tested considering inaccuracies in the predicted PV output power.
The PDS uses a prediction window W to extrapolate the PRR considering the present output power P out [k] and the predicted PV output power P pv [k + W] at the end of the prediction window. To demonstrate the concept of the PDS method, and example with T s = 1 min, PRR max = (10%)P max /min, and W = 10 is presented in Figure 4, where P max is the maximum contracted power. The prediction horizon is located between the present time t ck and the final time t fk such that WT s = t fk − t ck . The control horizon moves with each sample instant. When a PRR violation event is detected in the prediction window, a line is projected between P out [k] and P pv [k + W]. This line is used to determine the PRR at the output (PRR out [k]) and the output power at the next time interval (P out [k + 1]). If P out [k] < P pv [k], the excess of energy is used to charge the energy storage system (ESS). Likewise, if P out [k] > P pv [k], the ESS is discharged. When no future PRR violations are detected, the output of the inverter is the output of the PV and the ESS will not be charged nor discharged. To guarantee that the PRR output of the PDS fulfills the regulations in the worst-case scenario, the prediction window W is selected such that: where T s is the sampling period in minutes. The predicted PRR of the PV array is given by: where i = 1, 2, . . . , W is the number of time periods in the future being forecasted. To detect if PRR pv [k + i] will be above PRR max at any time of the prediction window, a detection variable V[k] is defined by: where: The PRR of the power plant is given by: In this case, ∼ P out is used as an auxiliary variable for the estimated output power of the power plant. It is important to note that the PDS does not read the output power and that it only requires the current and predicted values of P pv . The output power of the ESS is then defined as P ess . A flowchart that describes the PDS algorithm is presented in Figure 5. The proposed PDS method can be used to reduce the ESS capacity requirements and battery degradation, reducing equipment costs and increasing the lifetime of the ESS. Simulation results are presented in the results section of this manuscript.

Methodology for the Levelized Cost of Storage Estimation
The Levelized Cost of Storage (LCOS) is an estimation of the costs associated with battery storage and is given in $/kWh/year [10,11]. The methodology to estimate the LCOS used in this work is based on the hybrid aging estimation made in [9]. This aging estimation defines the battery capacity fades associated with the storage time (calendar degradation) and with the number of complete electrical cycles of the battery (cycle degradation) for a Lithium Iron Phosphate (LFP) battery [41], which is the same battery type used in Coto Laurel: Constants α t = 3.087×10 −7 , α NC = 6.87×10 −5 , β t = 0.05146, β NC = 0.027, and n = 0.5 were calculated empirically. The time t is given in months and the temperature T is given in Kelvin. The number of normalized cycles NC is a compilation of the number of cycles per each DoD, which are calculated using the rainflow-counting (RFC) algorithm [9,42] and the Palmgren-Miner rule [43]. Then, the LCOS is calculated based on the estimated capacity degradation.

Battery Degradation Estimation for PRRC
The RFC algorithm counts the number of cycles per DoD. This computation is performed using the Matlab "rainflow" function, which receives the time vector with the SOC values and returns a vector N cyc with the number of cycles per DoD. The battery degradation parameter over the simulated time is obtained by using: where N max is the maximum number of cycles per each DoD before achieving the end-oflife (EOL) and is specified by the manufacturer. Figure 6 shows the curve of N max against DoD for the Intensium Max High Power 20P cells LFP type from SAFT batteries, which are the battery banks used in Coto Laurel. This curve is approximated by: Figure 6. Maximum number of cycles to achieve EOL for different DoD values for LFP batteries. The EOL capacity is defined by the manufacturer to be 70% of the initial capacity.
Once the overall battery degradation parameter D is obtained, it is required to transform it to a normalized number of cycles NC. Assuming a mean DoD of 80%, Figure 6 shows that the maximum number of cycles allowed before reaching the EOL is 10,000. Thus, the normalized number of cycles can be calculated using: Then, the NC value is evaluated in (19) to obtain the capacity fade in one year. To obtain the number of years to achieve EOL considering calendar capacity fade, it is necessary to solve the following equation: where RC is the residual capacity (70% in this case), T is the temperature and y EOL is the number of years before achieving EOL. For this work, the operating temperature is assumed constant at 299.95 • K or 26.8 • C.

Levelized Cost of Storage (LCOS) Estimation
The LCOS is the cost of using energy storage per year. In the case of ESS for PRRC, it is the sum of the initial battery investment and the operational costs divided by the lifetime utilization factor (LUF) [9,10]. The LUF is a time variable that considers capacity fading and AC-to-AC efficiency n y , which is 92% for the type of batteries used in Coto Laurel (see Section 5.1). The equation for the LCOS is defined by: where I o is the initial investment of the complete ESS in $/kWh, op cost is the operational cost assumed as 3.5% of the initial investment with a discount rate of 4%, C y = 100% × (1 − D) y is the remaining capacity at year y, and ∆t is the time step of the estimation in years (1 in this case) [9,10]. It can be seen from (24) that the LCOS is directly affected by the lifetime expectancy y EOL and the degradation parameter D, which are calculated using the procedure presented in the previous subsection.

Ramp Violations
Ramp violations are difficult to monetize because it varies depending on the established contract between the power company and the government. However, these violations are counted in this work to show the effectiveness of each PRRC method. The violations are counted as the cumulated measure in minutes when the PRR is above the requirements from PREPA's MTR.

Coto Laurel Case Study
In this section, the Coto Laurel power plant is presented. A proposed model was developed and validated using the OP5700 simulator from Opal-RT Technologies. For this validation, the simulated data were compared against measured data from the real Coto Laurel power plant.

General Description of Coto Laurel Solar Power Plant
Coto Laurel's PV arrays are distributed into two feeders F1 and F2, each feeder is composed of seven and eight groups, respectively. The inverters used in the Coto Laurel project for both solar panels and batteries are the "Ingecon Sun 880TL U X360 Outdoor". For this project, a total of 21 identical inverters were used, along with 3 MWh of LFP batteries (Saft Intensium Max 20P Energy Storage Container with 8ESSUs) connected to the feeder F3. The model of the LFP batteries was developed using lookup tables based on the manufacturer's datasheet. The technical characteristics of the ESS and the PV array are described in Tables 1 and 2, respectively. Additionally, the one-line diagram of the solar farm is depicted in Figure 7.

Environmental Characteristics of Coto Laurel Solar Power Plant
Several studies have shown that the geographic position of Puerto Rico (latitude 18 • ) is an excellent location for the implementation of photovoltaic plants since it receives approximately five hours of direct peak sunlight every single day of the year [44]. Additionally, Puerto Rico is in the tropical region, where no significant variations in solar irradiance can be observed throughout the year. In other words, the irradiance monthly variation is negligible as depicted in Figure 8. However, significant hourly variations can be observed throughout the day as shown in Figure 9. These hourly variations generate inconsistencies in the energy production and may generate stability issues in the main grid.

Proposed Model
The implemented model describes the dynamics of the PV cell injection, considering the voltage of each array and the irradiance and temperature variations. The inverter was modeled as a continuous ideal AC current source. Additionally, a Direct-Quadrature (DQ) transformation was implemented to control the grid's power injection. Finally, a phase-locked loop (PLL) was used to perform the synchronization with the grid.
The configuration of each set of arrays was carried out using the information provided by the Coto Laurel farm. The parameters of the PV panels were extracted from the corresponding solar panel's datasheet, using information available in the standard test conditions (STD), and nominal operating cell temperature (NOCT). Additionally, the irradiance profile was measured in a meteorological station located near the system's arrays. Moreover, the voltage of each set of arrays was obtained from the information provided by the solar farm. Using the farm's measured parameters, a model was created using Matlab/Simulink. This model is composed of solar panels, inverters, transformers, distribution lines, and batteries.
Coto Laurel's lithium-ion battery bank was modeled using an equivalent circuit model. This model is composed of a voltage source with a series resistor, which represents the battery's open-circuit voltage (OCV) and internal resistance, respectively. Moreover, two RC networks in series were used to represent the electrochemical and concentration polarization of the battery. The voltage source was modeled using a lookup table that describes the relationship between the OCV and the battery level, which can be used to calculate the dynamic values of the battery's internal resistance and the two RC networks.
Considering the nature of this study and the model's time response (e.g., minutes), a continuous ideal AC current source was used to model the inverters on the entire system. This model is suitable for the proposed analysis since the time response of the power electronic elements and the applied current control techniques is faster (e.g., milliseconds) than the exhibited by the model. Moreover, considering the previous assumptions, a PLL with a continuous ideal AC current source was used, since it offered an accurate approximation for the current system.
The distribution and internal transmission line connections for the different elements (PV arrays, inverters, and transformers) were modeled using the branch block of the Simulink's SimPower systems library. Moreover, a resistive three-phase branch was used to model the transmission lines available on the Coto Laurel farm, since this model provides accurate approximations for transmission lines shorter than 80 Km.

Model Validation
To test and validate the proposed model, a comparison between the active power calculated using the proposed model and the real data measured by the Coto Laurel's POI's meter was performed. To perform the simulation, the yearly average irradiance per minute was used as an input parameter. To improve accuracy and testing time, the simulation was carried out in the OP5700 simulator from Opal-RT Technologies. After running the simulation, the average power per minute was obtained, which closely matches the average power data measured at the Coto Laurel farm. Figure 10 shows the simulated mean power (blue), the measured mean power data (orange), and the absolute error (green) for each minute. The average error during the day is approximately 3.05%. It is important to highlight that at the beginning and at the end there are the highest errors of the model, this happens due to the transient that occurs in the simulation due to the change in irradiance from 0 to a certain value. However, the model demonstrates to be accurate during the rest of the day. Results from this study validate the accuracy of the model describing the dynamics of Coto Laurel.

Simulation Results
Each of the PRRC methods presented in Sections 2 and 3 was implemented on the model presented in Section 5 using the OP5700 real-time simulator from Opal-RT Technologies. Simulation results were obtained in two stages. First, three variations of each PRRC method were evaluated for one day. The proposed PDS method was also evaluated considering variations in the prediction window and the accuracy of the prediction. In this stage, the performance indicators were the SOC variations, the SOC at the end of the day, the PRR, the PRR violations, the number of normalized cycles, the battery degradation, and the variance of the output PRR. Second, the best variation of each PRRC method was tested during an entire year to evaluate the PRR violations, the LCOS, and the degradation of the ESS.

First Stage: One Day Evaluation
The solar radiation profile used for the one-day evaluation is shown in Figure 9. This profile corresponds to 18 May 2020. This day was selected since it combines rough and smooth variations with high radiation intervals.
The same radiation profile was introduced to the modeled PV arrays shown in Figure 7. If the SOC was above 100% or below 30%, the battery banks stopped absorbing or injecting energy, respectively. The initial SOC was set to 60% to reduce the probability of reaching the SOC limits. The sampling period of all the simulations was set to 1 min.

Ramp Saturation Method
The simulation results for the RS method are shown in Figure 11. It can be noted that, as PRR max increased, the SOC variations and the power at the ESS (P ess ) were reduced. A summary of the simulation results for the RS method is shown in Table 3.  Simulation results for the SMA are shown in Figure 12. It is noted that, as the window size is increased, the output PRR is reduced and the SOC variations increase. A summary of the simulation results for the RS method is shown in Table 4. A window size of W = 10 causes half of the degradation compared to W = 20. Furthermore, the variance of PRR out for W = 10 is more than three times less than W = 5.  Simulation results for the EMA method are presented in Figure 13. All the simulations were made using a window size of 30. Although output and ESS powers are similar for each α, it is noted that the SOC reached 100% at t = 13:12 for α = 0.113. Thus, the ESS stopped absorbing energy and PRR violations started to appear. A summary of the simulation results for the RS method is shown in Table 5. There are no violations with α = 0.123 since it has the least PRR variance in all three tests.  Simulation results for the FOLPF method are shown in Figure 14. It is noted that, although results were similar, there was a violation at t = 10:15 for T f = 450. A summary of the simulation results for the FOLPF method is shown in Table 6. The least degradation without having PRR violations occurred when T f = 550.  Simulation results for the SOLPF method are shown in Figure 15. It is noted that the SOC reached the upper limit twice during the day for ω n = 1/35. Furthermore, there were multiple violations for ω n = 1/35. A summary of the simulation results for the SOLPF method is shown in Table 7. The least degradation without having PRR violations occurred when ω n = 1/25.  Simulation results for the ELES method are shown in Figure 16. It is noted that, although results were similar, there was a violation at t = 9:50 for α = 0.07. A summary of the simulation results for the ELES method is shown in Table 8. The least degradation without having PRR violations occurred when α = 0.06.  The PDS was evaluated by varying the prediction window size W (5 and 10). However, since this method relies on the accuracy of the predictions, its robustness against prediction inaccuracies was also evaluated. For this purpose, the assumed prediction was contaminated with Gaussian noise as follows: where rand is a random scalar with standard normal distribution and µ(%) is the variable noise parameter in the percent of the contracted output power. This means that a µ = 20% incurred in an estimated power estimation error of approximately 2 MW. Simulation results for the PDS method with W = 5 are shown in Figure 17. It is noted that, compared to other methods, the SOC variations got reduced and there were no PRR violations. It is also noted that the SOC variations are critically affected by the µ parameter. A summary of these results is presented in Table 9. The PDS shows to be robust against an estimation noise of at least 20%. The battery degradation is also small compared to other methods.  Simulation results for the PDS method with W = 10 are shown in Figure 18. It is noted that this window size significantly reduced the PRR and SOC variations. It is concluded that, as the window size increases, the PDS is more robust against uncertainties in the predictions. It is also noted that the SOC variations are critically affected by the µ parameter. A summary of these results is presented in Table 10. Although none of the tests showed PRR violations and the PRR variance was significantly reduced, a PDS method with a W = 10 showed twice the battery degradation compared with a window size of W = 5. Figure 19. Comparison of the effect of selecting the window size using the PDS method. Figure 19 shows the main effects of selecting the window size.   19. Comparison of the effect of selecting the window size using the PDS method.

Economic Analysis through One Year
From the previous subsection, one parameter configuration was selected for each PRRC method to be tested using a whole year's data. These configurations were selected considering the minimum amount of degradation without having PRR violations. Additionally, the mean of the SOC and the variance of the output PRR were considered. The SOC at the end of the day is an indicator of SOC variations and is also used as a selection criterion. With these criteria in mind, the following PRRC methods were considered for this section: • RS with a PRR max of 10%. First, a battery degradation analysis was made to each of the PRRC methods. Then, the estimated cost per kWh per year was compared. Finally, the violations throughout the year were analyzed. The solar radiation data for this test was selected from 20 May 2019 to 19 May 2020. Results were summarized at the end of this section. It is important to mention that, as the SOC at the end of each day is not the same as the one at the beginning of the day, the mean SOC tended to drift during the year until it reached the maximum or minimum limits. To account for this issue, a proportional SOC controller was added at the output of each PRRC controller. The mathematical expression for this controller is dictated by: where P ess is the adjusted power reference for the ESS, P ess is the output of the PRRC method, SOC ref is the desired SOC at the end of the day, and K soc is the proportional control constant, which is selected small to not affect the dynamics of the PRRC.

Battery Degradation
The summarized number of cycles using the rainflow algorithm are depicted in Figure 20. The cycles are a histogram representing the number of cycles per DoD during an entire year for the Coto Laurel plant under different PRRC methods using the rainflow algorithm. The DoDs were accumulated in intervals of 10%. With this information, it was possible to determine the capacity degradation and the lifetime of the ESS by using  According to Figure 21, The PDS showed to be the PRRC method that offers the most estimated battery lifetime (118 months), followed by the RS (103 months), the SMA (95 months), the ELES and EMA (91 months), the FOLPF (79 months), and the SOLPF (74 months). This was because the PDS method produced fewer battery cycles for high DoD values. According to Figure 6, the battery life is exponentially affected as the number of cycles of high DoD values increase.

LCOS
The initial investment for the LCOS analysis was determined based on [9], which used as base various reports from international technology centers and specialized consultancies. According to this, the estimated initial investment for a 3 MWh LFP ESS is about 2.35M€ or 2.87M$ (USD as of May 2021 exchange rate). Figure 22 shows a comparison of the LCOS (in $/kWh/year) for each PRRC method. These results were obtained using (24). These results were consistent with those presented in Figure 21. Results from Figure 22 indicate that the LCOS is highly dependent on battery degradation, which directly affects the battery lifetime and the capacity degradation of the ESS. One of the main advantages of the proposed PDS is that the DoD is considerably reduced by charging or discharging the ESS before a PRR event arrives. This also helps to reduce cumulated SOC drifts throughout the year. In order to increase battery lifetime, it is also recommended to reduce the amplitude of the SOC fluctuations. Furthermore, it is important to consider the penalties for violating the PRR requirements. This must be considered along with the LCOS to decide which PRRC method is more convenient for a particular LSSPP.

Violations
It is difficult to assign a monetary value to the PRR violations since this is unique for each agreement. However, PRR violations throughout the year are presented in this work to serve as an indicator of the accuracy and performance of each PRRC method. Figure 23 shows the cumulated measurement of the minutes when the output PRR was greater than PRR max . It is noted that most of the PRR violations in Coto Laurel occurred during the summer, which is characterized for having the highest solar radiation profiles according to Figure 8. The main cause of these violations over the year was because the ESS reached its minimum (30%) or maximum capacity (100%). This implied that the ESS was unable to support the PRRC. However, the PRR violations were relatively small, considering that the RS was the PRRC method with the most violations throughout the year (84 min). The PDS and the SOLPF were the PRRC methods that did not show any PRR violations throughout the year.

Discussion
The model of the Coto Laurel solar power plant was developed, tested, and validated using an OP5700 simulator from OPAL-RT technologies. Simulation results throughout the year validated the accuracy of the proposed Coto Laurel model. With this model, a novel PRRC method was proposed. This proposed method and other PRRC methods were tested using daily and yearly data. Daily results indicate that there is a tradeoff between PRR regulation and DoD cycles. It is important to determine the storage capacity required to optimize PRR support and battery degradation. For this work, the marginal required storage capacity was used. However, less degradation could be appreciated when this capacity is increased. Furthermore, as the storage capacity is increased, the PRR violations are expected to decrease. This is because it is less probable to reach the minimum or maximum SOC values. When this event occurs, mainly during summer, the PRRC method is unable to support the PRR regulation.
Among the PRRC methods found in the literature, the RS showed to be the one with the least LCOS value. However, there are some days when the ESS may reach the maximum or minimum value, limiting the capacity of the RS method to support the PRRC. The ELES, SMA, and EMA showed similar LCOS values. However, the EMA showed the most violations of these three methods. Although the SOLPF did not show PRR violations, it is the PRRC method with the highest LCOS, followed by the FOLPF.
The proposed PDS method showed the least LCOS values among all the PRRC methods. This was because the PDS used predicted data to charge or discharge the ESS to reduce DoD. This helped to reduce high-valued DoD cycles, which considerably increased the battery life while regulating the PRR without violations. Moreover, the PDS demonstrated to be robust against solar prediction inaccuracies of at least 20% of the real value with a minimum prediction window of 5 min. At shorter sampling periods, the algorithm would have more resolution correcting sudden PV power drops. However, if the window size is not increased, more violations would be expected, and the degradation would decrease because the ESS will not be used frequently because their output will not be smoothed. It is important to note that the selection of the sampling period was made based on the regulations, which dictated the PRR limit in minutes.
It is important to consider that the predictive nature of the PDS may incur extra costs depending on the selected prediction method. There are plenty of solar prediction methods available in the literature that may be used for the PDS implementation [6,35,36,39,40]. However, it is important to consider the accuracy, the costs, and the effectiveness of these methods against radiation profiles found in places such as Puerto Rico. As future work, it is of interest to study the complete implementation of the PDS method in a real solar power plant using one of the solar prediction methods found in the literature.

Data Availability Statement:
The data available in this study is openly available in IEEE Dataport, doi:10.21227/kkj4-fh6.