Next Article in Journal
Marine Energy Environmental Permitting and Compliance Costs
Next Article in Special Issue
Challenges in All-DC Offshore Wind Power Plants
Previous Article in Journal
An Insight to the Degradation Behaviour of the Parallel Connected Lithium-Ion Battery Cells
Previous Article in Special Issue
European and Indian Grid Codes for Utility Scale Hybrid Power Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Energy Arbitrage Using Energy Storage Systems: A Wind Park’s Perspective

1
Department of Energy Technology, Aalborg University, 9220 Aalborg, Denmark
2
Department of Mechanical and Manufacturing Engineering, Aalborg University, 9220 Aalborg, Denmark
3
Lead-Performance Optimization, Blue Power Partners, 9000 Aalborg, Denmark
*
Author to whom correspondence should be addressed.
Energies 2021, 14(16), 4718; https://doi.org/10.3390/en14164718
Submission received: 14 June 2021 / Revised: 20 July 2021 / Accepted: 27 July 2021 / Published: 4 August 2021
(This article belongs to the Collection Women's Research in Wind and Ocean Energy)

Abstract

:
With the growing application of green energy, the importance of effectively handling the volatile nature of these energy sources is also growing in order to ensure economic and operational viability. Accordingly, the main contribution of this work is to evaluate the revenue potential for wind parks with integrated storage systems in the day-ahead electricity markets using genetic algorithm. It is achieved by the concept of flexible charging–discharging of the Energy Storage System (ESS), taking advantage of the widespread electricity prices that are predicted using a feedforward-neural-network-based forecasting algorithm. In addition, the reactive power restrictions posed by grid code that are to be followed by the wind park are also considered as one of the constraints. Moreover, the profit obtained with a Battery Energy Storage System (BESS) is compared with that of a Thermal Energy Storage System (TESS). The proposed method gave more profitable results when utilizing BESS for energy arbitrage in day-ahead electricity markets than with TESS. Moreover, the availability of ESS at wind park has reduced the wind power curtailment.

1. Introduction

Increasing the penetration of Renewable Energy Sources (RES) is of high importance, as many countries have already set their goals for increasing RES share in the energy mix [1,2]. The technical and economical analysis of both standalone and grid-connected RES is presented in [3,4]. Climatic conditions highly influence the wind farm output that often causes voltage and frequency deviations in the connected power grid. Using the governor actions to solve these instability issues due to either shortage or surplus wind power would be detrimental to speed-control systems [5]. For example, in Denmark, 10% of hourly small variations are most common for wind farms, but up to 30% of hourly fluctuations still occur once in 1.1 years [6]. Energy storage systems not only enable the effective integration of widespread RES into power systems, but also provide various services including load shifting, energy management, frequency regulation, wind power smoothing, or energy arbitrage to the power grid. Among various Energy Storage System (ESS) technologies, Battery Energy Storage System (BESS) turned out to be a suitable choice for integration of wind farms and for maximizing profit [7]. There has been a large interest towards investing in BESS for a decade; for example, Tesla completed installation of a 129-MWh, grid-scale, lithium-ion battery connected to Neoen’s Hornsdale Wind Farm near Jamestown, South Australia in 2017 [8]. Further, ABB won an order from DONG Energy, the leading Denmark-based sustainable energy group, to deliver a 2-MegaWatt BESS to support the integration and transmission of power from the Burbo Bank offshore wind farm near Liverpool, UK [9].

Study of Literature: Battery Energy Storage System (BESS) and Thermal Energy Storage System (TESS)

The Energy Storage System (ESS) plays a key role, even in the reserve market; however, cross border and intraday trade have lessened the potential for ESS in the Nordic intraday market [10], resulting in a lower return on investment compared to the day-ahead market. The most significant revenue stream for BESS is energy arbitrage in short-term electricity markets, such as the day-ahead electricity market. In general, the arbitrage is a transaction that generates revenue by using price differences over a period, where the energy is stored during the low-demand periods or low energy prices and dispatches during high-demand periods or high energy prices [11]. For storage units, real-time arbitrage is a major source of revenue, but the uncertain nature of prices pose difficulties in designing good strategies. The difficulties arise from the fact that the future electricity prices are unknown, highly stochastic, and difficult to forecast as well [12]. In general, the demand dynamics dictate the price dynamics in a cyclical pattern, which can be daily, weekly, seasonal, and with some unpredictable factors. Application of appropriate optimization methods would maximize the returns from the combination of wind farms and ESS.
The general solutions that can be employed for handling the intermittency of the wind power systems include flexibility technologies such as demand response technologies [13], energy storage technologies [14], etc. Several research works demonstrated that the effective coordination of Wind and ESS can handle the future fluctuations in wind power generation [15]. In [16], BESS is used to reduce the shortfalls of load, and reserves are reduced along with a decrease in wind power curtailment in day-ahead scheduling. It is assumed that the BESS is a system asset operated by the system operator. A comprehensive analysis of the fallouts of grid-scale ESS in a real thermal power system under a centralized spot market is carried out in [17]. The ESS is employed to provide energy arbitrage, and primary and secondary reserves. A novel wind dispatch control scheme is proposed in [18], where the major decision lies in discharging the BESS during peak price times and storing it at off-peak periods in agreement with the Australian National Electricity Market. Further, the price arbitrage-based feasibility study is carried out for BESS in relation to a RES park in [19]. In the literature, several tools have been developed to analyze the economic viability of the wind-storage systems [20]. The storage technologies used for integrating renewable generation are not only in the form of electrical energy, but also chemical and thermal energy [21]. Thermal energy storage systems, especially with latent heat storage methods, have also been proven to be a competing energy dispatcher, due to their large-scale, being environmentally friendly, and long life [22]. In TESS technologies, latent heat storage is preferred over sensible due to its high storage capacities. The molten-salt-based phase change material (PCM) is used for energy storage in the form of latent heat at high temperatures due to its recognition as the best heat transfer fluid (HTF) in Thermal Energy Storage System (TESS) [23] and the cost of molten salt being low. In countries such as Denmark, where more than 60% of heating in Danish houses is supplied by district heating and, in addition, high wind power penetration, a considerable amount of wind energy can be used to support DH system in low electricity price periods [24]. For half of the concentrated solar power (CSP) throughout the world, TESS is used as a viable storage [25]. Moreover, thermal energy present in the solar irradiance can be stored in the TESS and used to meet the heat load [26].
The present work prioritizes the use of battery and thermal energy storage systems for maximizing the profit of the wind farm operator obtained by selling the wind energy through energy arbitrage. The problem is formulated as a maximization function that measures the economic benefits obtained from the dispatched power from the wind farm against the cost of ESS. Both BESS and TESS types are considered and their effect on the wind park’s profit is evaluated. Considering the point of common connection (PCC) node prices, the ESS is scheduled for either charging/discharging. According to the grid codes of various countries, for maintaining power factor at the PCC as a function of voltage, the wind farms not only supply active power but also reactive power into the grid [27]. Accordingly, the reactive power compensation at the PCC is considered one of the constraints, which is a regulatory imposition by the Transmission System Operator (TSO) to be followed by the wind park operator for injecting its power production into the grid. The modeled wind farm is used to find the reactive power set points and passed onto optimization routine in the form of a constraint. This is an important aspect, as in the process of maintaining the reactive power at the PPC and within wind farm, there can be a possibility to adjust the real power output of the wind farm, which has a direct effect on energy sold to the grid, thereby on the hourly profit. However, the wind turbine dynamics are not considered in this work, as the time scale is hourly precision. The dynamics, including electrical converter controls and turbine mechanical controls, is a study of a seconds or less time range but this is out of the scope of present work. The advantage of this work is that both the performances of BESS and TESS are compared with respect to the wind park’s revenue. The novelty of this work lies in considering the energy storage arbitrage along with grid code fulfillment in terms of reactive power response for finding a trade-off solution to maximize the revenue generated by selling wind energy to the grid. In addition, the neural network model proposed in this work for forecasting the electricity price is compared with the methods that are proposed in literature [28,29] and were revealed to be better.
Figure 1 shows the generic wind park model with the ESS. The ESS (either battery or thermal) is installed at the wind park site and connected to the main grid through a charge controller/inverter. The energy capacity of ESS describes the energy stored in the battery for a given hour and the power that can be supplied or stored under the rated charge/discharge time interval. The selected power and energy rating are given in the wind park modeling section. The details of the BESS and TESS mathematical models are described in the following subsections.

2. Battery Energy Storage Systems

There are several BESS models reported in the literature. The charging/discharging process for both the long-term and short-term applications is illustrated in [30]. The present work considers the long-term operation (e.g., over hours/days) for BESS/TESS charging and discharging processes, which is relevant for energy arbitrage applications. Besides this, a numerical model in [31] is adapted for charge/discharge for BESS modeling, as given in Equation (1).The equation gives the state of charge (SoC) at any time step k with respect to charge ( P c h , B E S S ) and discharge ( P d c h , B E S S ) instances corresponding to given charge ( η c h )/discharge ( η d c h ) efficiencies. The SoC is having upper limit ( S O C m a x ) and lower limit ( S O C m i n ) and the SoC for any time interval ( Δ t ) are limited by charge ( r c h )/discharge ( r d c h ) rates.
S o C ( k + 1 ) =   S o C ( k )   +   P c h , k η c h , B E S S Δ t E B E S S , k S o C ( k + 1 ) =   S o C ( k )     P d c h , k η d c h , B E S S Δ t E B E S S , k S o C m i n     S o C ( k )     S o C m a x r c h     S o C ( k + 1 )     S o C ( k )     r d c h

3. Thermal Energy Storage Systems

Thermal energy storage systems are considered for wind parks due to many factors including cost effectiveness, large-scale potential, and environmental friendliness. The most widespread configuration of TESS is the two-tank sensible heat system with Solar salt as Heat Transfer Fluid (HTF) and a storage medium due to its high thermal efficiency and strong operability [32]. The molten-salt-based TESS system is one of the best HTFs due to its low cost, which makes it economically attractive. The present work considers a thermal storage model proposed by Stiesdal, currently under development [33], which is basically a two-tank model with a compressor connected to a generator with a governor system. End-use TESS stores electricity from off-peak periods in the form of heat using hot or cold storage tanks in underground aquifers, water or ice tanks, molten salts, or other storage materials. Heat losses would be there, which can be considered in cost calculation. In the present model, the associated delays are due to the ramping time of the generator dynamic system. The state of energy (SoE) that is stored at any time step k is derived from the total thermal energy ( E c a p ) and current energy ( E k ) available in the TESS. The binary variable ( u k ) is used to avoid simultaneous charge/discharge and ( l k ) is accounting for the internal losses within the thermal storage tank. are given in Equation (2):
E ( k + 1 ) =   ( 1     u k l k ) E ( k )   +   η c h , B E S S P c h , k Δ t     P d c h , k η d c h , B E S S Δ t S o E ( k ) =   E ( k ) E c a p
The power capacity ( P c a p ) of both the thermal and battery storage systems is used in formulating the charge and discharge constraints, which are given in Equation (3):
P c h , k     u k P c a p P d c h , k     ( 1     u k ) P c a p
Unlike BESS, the thermal storage systems have a ramp constraint, expressed in terms of SoE, which is given in Equation (3). Further, the constraint on the minimum energy ( E m i n ) and maximum energy ( E m a x ) that TESS can handle is given in Equation (4):
R k Δ t     S o E ( k + 1 )     S o E ( k )     R k Δ t S o E m i n     S o E ( k )     S o E m a x E m i n     E k     E m a x
In the present work, the ramp rate ( R k ) is seen as rate at which the dispatch of power takes place from TESS to grid.

4. Problem Statement and Optimization Formulation

The optimization formulation was solved using a Genetic Algorithm (GA) due to its parallel optimization capability, thereby giving a group of solutions at once. In addition, the conventional optimization solvers might end up in local minima and oscillation effects, thereby failing to arrive at global optima. The Darwinian principle of natural evolution was adapted in GA involving a population, which finds the solution for a given problem [34]. In comparison with other global search algorithms, GA proves to be superior over global convergence. This method solves the problem by considering an appropriate selection of mutation and cross-over operators, thereby updating the individual population in each iteration [35].
The objective is to maximize the revenue of the wind park, as shown in Equation (5), using intelligent charging/discharging based on predicted electricity prices for relevant constraints on storage systems including their SoC limits, ramp rates for thermal storage, charge/discharge rates and reactive power limits (( Q m i n ), ( Q m a x )). The power management strategy for ESS ( P E S S , i ) is based on given wind production profiles ( P W , i ) and forecasted electricity prices ( E p r , i ) by the maximization optimization of the revenue.
M a x i = 1 24 E p r , i ( P W , i + P E S S , i ) C o n s t r a i n t s E q u a t i o n ( 1 )     B E S S E q u a t i o n s ( 2 )     ( 4 )     T E S S Q m i n     Q i , W     Q m a x
From the above formulation, the revenue obtained by selling energy from both wind and ESS ( E E S S ) to the grid is obtained. The net profit is obtained from the difference of revenue and the costs that is given in Equation (6). Since the focus of this paper is on energy storage facility, the corresponding capital costs of both the wind farm and converters are assumed to be constant. This equation describes the costs of both the TESS and BESS, whereas the cost of storing energy is comparatively less than that of the cost of operating for the TESS system. This is due to the fact that the TESS consists of a steam turbine along with a generator for generating electricity from thermal form. So, the operating cost of TESS is assumed to be high compared with BESS.
P r o f i t   =   α P s o l d     β P E S S     γ E E S S
In this study, it is considered that the predicted values of wind speed profiles are already available and electricity prices are forecasted from historical data using neural networks. The real data for a period of two years (2018–2019) corresponding to wind speed profiles of the wind farm under study and electricity nodal pricing is provided by our project partners. The forecasting is carried out only for electricity prices, with price being the main decision maker in the profit maximization.
The simulation model is as shown in Figure 2, illustrating the neural-network-based forecasting model for day-ahead electricity prices. The input for the model contains hour of the day, day of the week, and time of year, which are used deterministically, and this dataset is divided into a training set to train the neural network for day-ahead price forecasting. The other dataset is the test set to validate the accuracy of the forecast model. The model accuracy for the given sample periods is computed with the Mean Absolute Error (MAE) and Mean Absolute Percent Error (MAPE) metrics. The forecasting toolbox available in Matlab is used to carry out the prediction algorithm.
Once the day-ahead electricity prices are forecasted, the GA-based optimization is carried out for two cases, (1) BESS and (2) TESS, and the working flowchart is as shown in Figure 3. The optimization formulation given in Equation (5) is solved once for BESS and once for TESS. Within the GA optimization routine, the control goes to the DIgSILENT (power flow program) routine to simulate the total wind farm model and obtain the values of reactive power to be inserted in the Matlab routine. In the profit maximization, along with constraints of BESS and TESS, meeting the reactive power requirements is also formulated as one of the constraints.

5. DIgSILENT Model of Wind Park

The DIgSILENT model of the wind farm illustrating the basic load flow results is shown in Figure 4. This model is one of the examples within the software, and the original wind farm under study is modeled in the same way. However, it is not shown here due to data protection policy. Parameters such as copper losses and short circuit voltage are calculated [36].
A fully rated converter-based wind turbine generator of 3.45 MW model is used to model the wind turbines. The maximum real power output of the park is around 197 MW. The reactive response to be seen at the PCC, which is a 345 kV connection that is governed by the Vestas Power Plant Controller (PPC) utilizing the capabilities of the turbines. As ordered by the TSO, the park controller sends the set points for real and reactive powers either directly to wind turbine units or through a cluster controller that is commonly used in large wind parks [37]. When reactive power supply is insufficient, the PPC will operate the capacitor bank installed at the wind park site, thereby avoiding drawing reactive power from the grid, which would otherwise have led to grid code violations [38]. The bank is sized into two steps of 12.6 MVar each-total 25.2 MVar compensation. The real data of wind speed profiles are given as input to the DIgSILENT models of wind turbines and the load flow simulations are carried out for determining the total wind park production at the PCC for a month, as shown in Figure 5.
The active power injection of the wind turbines mainly depends on the wind speed. The real power production at PCC and real power losses within the wind park with respect to wind sweep speeds are shown in Figure 6.

Reactive Power Capability

The concept of capability curve is used to determine the active and reactive powers that can be delivered by a wind farm. The real power–reactive power (PQ) diagrams depict the performance of the individual turbines or wind farm as a whole based on the limits of active and reactive powers. It is possible to control the reactive power at PCC with the park controller on-site. The reactive power capability of the wind park without and with park controller are analyzed and the results are represented as Case 1 and Case 2, respectively. The reference curve for reactive power capability (red color) of the wind park for a given reactive power requirement (blue color) without a park controller system is shown in Figure 7. The active power is swept from the starting value to maximum using two cases, i.e., Q m i n and ( Q m a x ) correspond to reactive power for underexcited and overexcited cases, respectively. This is formulated as a constraint in the optimization problem given in Equation (5). The DIgSILENT output is given in Table 1. In Figure 7, the blue curve (requirement of reactive power) is exceeding the red curve (capability of the wind park) for both underexcited (left side of origin) and overexcited (right side of origin) cases.
It can be seen from both Figure 7 and the DIgSILENT output data in Table 1 that the wind park is unable to meet the required reactive power in few operating conditions. For overexcited case, the reactive power requirement is 58.99 MVar, whereas the wind park capability is only 41.37 MVar. Likewise, during underexcited cases, the wind park could not meet the required reactive power. Therefore, the park controller is made responsible for meeting the reactive power requirement either by switching the capacitor bank or modifying the active power output of the wind park.
The modified reactive power capability of the wind park by park controller in accordance with the reference (blue color) is as shown in Figure 8. The capacitor bank is switched on by PPC in order to meet the reactive power requirement. The park controller present in the wind farm switches ON the capacitor bank of 25.2 MVar in two steps depending on the reactive power requirement at the PCC. It can be observed from the DIgSILENT report that the capability is increased to 82.58 Mvar for underexcited and 83.14 Mvar for overexcited cases compared with the results obtained without park controller. As can be observed from Figure 8, the blue curve (requirement) is well within the red curve (capability of the wind park) for both underexcited (left side of origin) and overexcited (right side of origin) cases. The DIgSILENT output is given in Table 2. The PPC finds new active and reactive power set points, thereby operating the capacitor bank in order to meet the reactive power requirements. It can be seen from the results that there is a fair margin between capability and requirements, enabling the wind park to accommodate any future rise in reactive power requirement.
Table 3 presents the simulation results for PQ-analysis carried on the wind park before and after PPC operation. It can be observed from the PQ-analysis results, that the maximum and minimum voltage at shunt terminal are above 1 p.u. in Case 2 and are higher than in Case 1 because the capacitor bank is switched ON, thereby injecting reactive power. Further, in Case 1, the transformer should be loaded to 100%, whereas it is only loaded to 99.9%, showing that the real power production of wind park is also reduced to meet the reactive power requirement. This reduction in active power is significantly lower for both cases, but may be significant during fault disturbances within or outside the wind farm. This is not of main interest for this work. The DIgSILENT model is used to simulate the system and obtain the operating points for the later optimization problem that will be solved in Matlab, where a cosimulation between Matlab and DIgSILENT is established.

6. Simulations Results

The performance of the profit maximization algorithm is validated for various seasons where electricity prices are widespread, and the comparative analysis is shown for both BESS and TESS. The power and energy ratings, SoC/SoE limits, and ramp rates for both the storage models are given in Table 4. Considering the hourly time-steps, the response time of BESS is only from 20 ms-s [39], where BESS characteristics show no delay in terms of ramp rate. In addition, the amortized cost parameters are taken from [6] for the present work. Considering the BESS costs, the investment cost is $29 million for a battery cost of 116 $ /kWh. This can be recovered in a span of 3–4 years with the proposed energy arbitrage solution for the wind farm.
The ability of deep learning regarding the input–output relationship for a given input data makes neural networks preferable over linear regression techniques. Figure 9 shows the electricity price results obtained from the forecasting algorithm. The residuals plot shows the estimated error between actual data and modeled data. The forecasting algorithm is based on Neural Network (NN)-based regression. The simple linear regression uses an approach where a line is used to best fit the given data for minimizing the error. However, when the data is scattered, this approach leads to either overfitting or underfitting. Even though the overfitting problem is also an issue in neural networks, in this work, it is handled by decreasing the complexity of the model through strategically reducing the number of neurons. There are many methods to resolve this problem, but this solution came out to be feasible for handling the present formulation. A three-layer feedforward neural network is used for forecasting the day-ahead electricity prices of hourly interval. The NN model consists of an input layer, hidden layer with 10 neurons, and output layer connected by weight functions. These weight functions are adjusted to better fit the model, thereby reducing the error. The model considers inputs including historical hourly prices, temperature, and holidays. The predictor matrix contains information about hour of the day, day in a week, previous day’s average price, and price from the same hour on the previous day. Then, this whole data is split into training for building the model and test data for forecasting the prices using estimated parameters from the training data set. The MAPE for a day is as given in Equation (7):
M A P E   =   1 24 i = 1 24 P i , a c t u a l P i , m o d e l P ¯ a c t u a l P ¯ a c t u a l   =   1 24 i = 1 24 P i , a c t u a l
It can be observed from Figure 9 that the model is able to predict the electricity prices using the historical data, where the model is able to track the actual data. There exists large MAE during the large peaks for certain periods, this is due to the fact that sometimes the NN model could not capture the sudden uncertainty in the electricity prices, which has no relation with the historical data. This uncertainty in electricity price can be caused by the demand, as electricity price is sensitive to load changes. The average MAPE is 6.89%, MAE is 5.27 MWh, and Daily Peak MAPE is 6.33% for the considered data set. The MAE provides information about the deviation in the estimated model and predicted model—the less the deviation, the better the predicted model. The average MAPE result obtained from the proposed NN model is 6.89%. Table 5 shows the comparison of the results from the proposed method with the results from the literature. It is to be noted that the MAPE obtained from the proposed method is lower, which shows the effectiveness of the proposed method. Figure 10 shows the daily forecast results corresponding to various seasons. The MAPE came out to be higher in the fall season, i.e., 7.05%, as shown in Figure 10c; it is worth noticing that the prediction is quite good during morning and evenings. However, the inaccuracy exists, especially during higher price ranges caused due to the uncertainty in the power demand. For the winter season, as shown in Figure 10d, the ANN forecast price is almost same as the actual price from morning till evening, this is due to the fact that the training data has better accuracy that could lead to increased efficiency in the forecasting.
The predicted electricity prices and the simulated wind production profiles are given as known variables and the ESS charge/discharge strategies are decision variables in the GA optimization routine. As a first step, GA initializes the population and creates the next generation, including elite children with better fitness function, cross-over children from a pair of parents, and mutated children with random changes of same parent. The cross-over generation extracts the best features from the different populations and mutation helps in adding diversity—this way, GA generates superior individuals with a better fitness function that increases the maximum likelihood. For these heuristic algorithms, the initial population plays a significant role in finding the feasible solution. The convergence of objective/fitness function for four days corresponding to four seasons is shown in Figure 11, Figure 12, Figure 13 and Figure 14. The selection function provides information about the selection of offspring from the cross-over and mutation of parents/individuals from the population of 200. The GA tool in Matlab is a minimization algorithm, however, the present study is a maximization problem. Accordingly, the objective function is given as a negative function, as is the negative value for the fitness function. The day in July could capture more profit compared to other seasonal days, which is a positive outcome from energy arbitrage using BESS. Further, the number of iterations that are needed by GA to meet the stopping criteria are different.
The results obtained from GA are shown in Figure 15a,b for the SoC (BESS: min—20% and max—90%) and SoE (TESS: min—0% and max—95%) of both the BESS and TESS, respectively. The depth of discharge is different for the two storage systems, which can be observed from charge and discharge profiles as shown in Figure 16 and Figure 17. The thermal storage system discharges fully to zero, whereas the battery storage system’s minimum discharge limit is 20% of the total capacity. Internal losses are comparatively higher in the thermal storage system than in BESS, as can be seen from the negative axes of Figure 16 and Figure 17. Even though the TESS could charge up to 180 MW, due to more internal heat losses, it could only discharge 150 MW at any given time interval.
The optimal fitness function values, which is the profit that is obtained without ESS and with BESS/TESS, for the year 2019 are given in Table 6. It is to be noticed that the profit obtained from BESS is more than that in the TESS case. Nevertheless, if the objective is a trade-off between being ecologically friendly and maximizing the profit along with leveled costs, then TESS could be a better solution in spite of P2X losses.
Without the introduction of ESS, the wind farm is subjected to curtailment due to wind speed prediction errors, reactive power restrictions, etc. In this study, the curtailment of wind power for the year 2019 is compared for three cases (No ESS, BESS, and TESS), given in Table 7. It can be observed that the curtailment is less in the presence of BESS than with TESS case. This is due to the fact that there exists operational delays, including ramp rate charge/discharge for TESS, which is negligible in BESS.

7. Discussion and Conclusions

Grid-scalable energy storage systems, including BESS and TESS, are utilized to study their impact on wind park owner’s revenue through selling power to the grid, taking advantage of the price arbitrage. The thermal storage system behavior is different to the BESS system, which is due to the ramp rate constraints and the higher operating losses. Even though the depth of discharge of TESS is more than that of BESS, due to internal heat losses in the conversion of heat to electricity, BESS gives more profit than TESS. In this work, the energy from the TESS is used to produce electrical energy, thereby selling to the grid. Whereas, there can be other option of meeting heat demand with the TESS. The main conclusions can be summarized as follows.
  • The present work deals with a real-world project, proposing an optimal strategy that helps the wind park owner to invest in large-scale ESS in order to prudently sell wind energy into the grid along with adhering to grid code requirements, thereby maximizing the profit.
  • The proposed forecasting algorithm using feedforward neural network obtained better results compared to the results from the literature. The overfitting problem of the NN is handled by strategically reducing the model complexity by decreasing the number of hidden neurons.
  • Profit maximization by GA mainly depends upon the accuracy of the electricity price and wind power forecasting, which are inputs to the optimization problem. In this work, the electricity price is considered as a major factor, so NN-based regression analysis is performed to have better prediction values. Furthermore, the mean absolute error is only 6.89% with the proposed NN model, which shows the accuracy of the predicted model regarding the profit.
  • The utilization of wind energy (reduced curtailment) is much better in BESS than in TESS because of the ramp rate constraint and also significant power to heat and heat-to-power losses. However, with respect to the cost associated with operation and maintenance, TESS would be in high competition to become the first choice for a wind park operator.
  • The reactive power support is provided by the wind farm, whereas the ESS is only responsible for real power charge/discharge with respect to electricity prices.
  • The integrated system either with the BESS or the TESS gives better revenue to wind park owners when compared to the case where there is no ESS facility.
  • The BESS of large capacity will be a better option for short-term or long-term energy arbitrage, i.e., in day-ahead electricity markets. However, the cost is a major factor in this decision.
  • Further, balancing markets and reserves could also be possible marketplaces for energy arbitrage with ESS.

Author Contributions

Data curation, J.A.; Formal analysis, P.P.; Methodology, P.P.; Project administration, J.A.; Software, P.P.; Supervision, B.B.-J.; Validation, P.P.; Writing—original draft, P.P.; Writing—review & editing, P.P., B.B.-J. and B.V.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding and the APC was funded by Electric power systems and Microgrids section in the Department of Energy, Aalborg University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The work is carried out as a part of EU project 886074:SST Energilagring, Vindmølleindustrien, EU, Regionalfonden.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
BESSBattery Energy Storage Systems
CSPConcentrated solar power
EESEnergy storage systems
GAGenetic algorithm
HTFHeat transfer fluid
MAEMean absolute error
MAPEMean absolute percent error
MWMegawatt
NNNeural Network
PCCPoint of common coupling
PCMPhase change material
PPCPower plant controller
PQReal power-reactive power
RESRenewable energy sources
SoCState of Charge
SoEState of Energy
TESSThermal energy storage systems
TSOTransmission system operator
Nomenclature
Δ t Time interval
P c h , k BESS charging power at time k
P d c h , k BESS discharging power at time k
S o C ( k ) BESS State of Charge (SoC) at time k
S o C m i n Minimum limit on SoC of BESS
S o C m a x Maximum limit on SoC of BESS
r c h Allowable charging rate
r d c h Allowable discharge rates
η c h , B E S S BESS charging efficiency
η d c h , B E S S BESS discharging efficiency
u k A binary variable used to avoid simultaneous charge and discharge
l k Internal losses occurring due to self-discharge in each time step
E c a p Energy capacity of the TESS system
S o E ( k ) SoE of TESS at time k
S o E m i n Minimum limit on SoE of the TESS system
S o E m a x maximum limit on SoE of the TESS system
R k Ramp rate
E p r , i Forecasted electricity market price
P W , i Wind production of the wind park
P E S S , i Available power within ESS
Q i , W Reactive power capability of the wind park
Q m i n Minimum reactive power capability of the wind park at its maximum active power output
Q m a x Maximum reactive power capability of the wind park at its maximum active power output
α Wind energy unit price in dollar per kilowatthour
β Amortized capital costs in dollar per kilowatt
γ Amortized capital costs in dollar per kilowatthour
P s o l d Wind power sold to grid
P B E S S Power rating of BESS
E B E S S Energy capacity of BESS
P i , a c t u a l Actual price for i t h hour
P i , m o d e l Forecasted price from NN model for i t h hour
P ¯ a c t u a l Average actual price for that day

References

  1. Beaudin, M.; Zareipour, H.; Schellenberglabe, A.; Rosehart, W. Energy for sustainable development energy storage for mitigating the variability of renewable electricity sources: An updated review. Energy Sustain. Dev. 2010, 14, 302–314. [Google Scholar] [CrossRef]
  2. Mehrjerdi, H. Modeling, integration, and optimal selection of the turbine technology in the hybrid wind-photovoltaic renewable energy system design. Energy Convers. Manag. 2020, 205, 112350. [Google Scholar] [CrossRef]
  3. Mazzeo, D.; Matera, N.; De Luca, P.; Baglivo, C.; Congedo, P.M.; Oliveti, G. Worldwide geographical mapping and optimization of stand-alone and grid-connected hybrid renewable system techno-economic performance across Köppen-Geiger climates. Appl. Energy 2020, 276, 115507. [Google Scholar] [CrossRef]
  4. Mazzeo, D.; Matera, N.; De Luca, P.; Baglivo, C.; Congedo, P.M.; Oliveti, G. A literature review and statistical analysis of photovoltaic-wind hybrid renewable system research by considering the most relevant 550 articles: An upgradable matrix literature database. J. Clean. Prod. 2021, 295, 126070. [Google Scholar] [CrossRef]
  5. Luo, C.; Ooi, B.-T. Frequency deviation of thermal power plants due to wind farms. IEEE Trans. Energy Convers. 2006, 21, 708–716. [Google Scholar] [CrossRef]
  6. Lipman, T.; Ramos, R.; Kammen, D. An Assessment of Battery and Hydrogen Energy Storage Systems Integrated with Wind Energy Sources in California; Technical Report; University of California: Berkeley, CA, USA, 2005. [Google Scholar]
  7. Divya, K.C.; Ostergaard, J. Battery energy storage technology for power system—An overview. Electr. Power Syst. Res. 2009, 79, 511–520. [Google Scholar] [CrossRef]
  8. ‘Undeniable Success’: South Australia’s 129 MWh Tesla Battery. Available online: https://www.energy-storage.news/news/undeniable-success-south-australias-129mwh-tesla-battery#:~:text=Hornsdale%20Power%20Reserve%20(HPR)%2C,into%20operation%20in%20November%202017 (accessed on 14 July 2020).
  9. ABB to Provide Innovative Energy Storage Solution for UK Offshore Wind Farm. Available online: https://new.abb.com/news/detail/46554/abb-to-provide-innovative-energy-storage-solution-for-uk-offshore-wind-farm (accessed on 14 July 2020).
  10. Zakeri, B.; Syri, S. Value of energy storage in the Nordic Power market-benefits from price arbitrage and ancillary services. In Proceedings of the 13th International Conference on the European Energy Market (EEM), Porto, Portugal, 6–9 June 2016; Volume 12, pp. 1–5. [Google Scholar]
  11. Hakim, M.; Danley, D.B.; Caplan, D.R. Energy Arbitrage by Load Shifting. U.S. Patent 20090326729, 31 December 2009. Available online: http://patentimages.storage.googleapis.com/pdfs/US20090326729.pdf (accessed on 14 July 2020).
  12. Kim, T.T.; Poor, H.V. Scheduling power consumption with price uncertainty. IEEE Trans. Smart Grid 2011, 2, 519–527. [Google Scholar] [CrossRef]
  13. Jiang, Y.; Xu, J.; Sun, Y.; Wei, C.; Wang, J.; Ke, D.; Li, X.; Yang, J.; Peng, X.; Tang, B. Day-ahead stochastic economic dispatch of wind integrated power system considering demand response of residential hybrid energy system. Appl. Energy 2011, 190, 1126–1137. [Google Scholar] [CrossRef]
  14. Sullivan, P.; Short, W.; Blair, N. Modelling the benefits of storage technologies to wind power. In Proceedings of the American Wind Energy Association Wind Power Conference, Golden, CO, USA, 1–4 June 2008; Curran: Red Hook, NY, USA, 2008; pp. 1–4. [Google Scholar]
  15. Chowdhury, M.; Haque, M.; Aktarujjaman, M.; Negnevitsky, M.; Gargoom, A. Grid integration impacts and energy storage systems for wind energy applications—A review. In Proceedings of the IEEE Power and Energy Society General Meeting, Detroit, MI, USA, 24–28 July 2011; pp. 1–8. [Google Scholar]
  16. Li, N.; Uckun, C.; Constantinescu, E.M.; Birge, J.R.; Hedman, K.W.; Botterud, A. Flexible operation of batteries in power system scheduling with renewable energy. IEEE Trans. Sustain. Energy 2016, 7, 685–696. [Google Scholar] [CrossRef]
  17. Suazo-Martínez, C.; Pereira-Bonvallet, E.; Palma-Behnke, R.; Zhang, X. Impacts of Energy Storage on Short Term Operation Planning under Centralized Spot Markets. IEEE Trans. Smart Grid 2014, 5, 1110–1118. [Google Scholar] [CrossRef]
  18. Khatamianfar, A.; Khalid, M.; Savkin, A.V.; Agelidis, V.G. Improving Wind Farm Dispatch in the Australian Electricity Market with Battery Energy Storage Using Model Predictive Control. IEEE Trans. Sustain. Energy 2013, 4, 745–755. [Google Scholar] [CrossRef]
  19. Daggett, A.; Qadrdan, M.; Jenkins, N. Feasibility of a battery storage system for a renewable energy park operating with price arbitrage. In Proceedings of the IEEE PES Innovative Smart Grid Technologies Conference Europe (ISGT-Europe), Turin, Italy, 26–29 September 2017; pp. 1–6. [Google Scholar]
  20. Dicorato, M.; Forte, G.; Pisani, M.; Trovato, M. Planning and operating combined wind storage system in electricity market. IEEE Trans. Sustain. Energy 2012, 3, 209–217. [Google Scholar] [CrossRef]
  21. Amrouche, S.O.; Rekioua, D.; Rekioua, T.; Bacha, S. Overview of energy storage in renewable energy systems. Int. J. Hydrog. Energy 2016, 41, 20914–20927. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Ding, W. Increasing Wind Power Consumption & Absorbability, Solving Grid Connection Issue by Use of Thermal Energy Storage Technology. In Proceedings of the Chinese Society for Electrical Engineering Annual Meeting, Beijing, China, 21–24 November 2012. [Google Scholar]
  23. Hasnain, S.M. Review on sustainable thermal energy storage technologies, Part I: Heat storage materials and techniques. Energy Convers. Manag. 1998, 39, 1127–1138. [Google Scholar] [CrossRef]
  24. Wu, Y.; Zhang, X.; Wang, H.; Sun, J. Molten salt heat storage and supply technology based on heating using abandoned wind power PV power or off-peak power. Sino-Glob. Energy 2017, 22, 93–99. [Google Scholar]
  25. Tian, Y.; Zhao, C.Y. A review of solar collectors and thermal energy storage in solar thermal applications. Appl. Energy 2013, 104, 538–553. [Google Scholar] [CrossRef] [Green Version]
  26. Zhao, Z.; Arif, M.T.; Oo, A.M.T. Solar thermal energy with molten-salt storage for residential heating application. Energy Procedia 2017, 110, 243–249. [Google Scholar] [CrossRef]
  27. On-Netz, E. Grid Code. High and Extra High Voltage. 2006. Available online: http://www.pvupscale.org/IMG/pdf/D4_2_DE_annex_A-3_EON_HV_grid__connection_requirements_ENENARHS2006de.pdf (accessed on 14 July 2020).
  28. Mandal, P.; Senjyu, T.; Urasaki, N.; Funabashi, T.; Srivastava, A.K. A Novel Approach to Forecast Electricity Price for PJM Using Neural Network and Similar Days Method. IEEE Trans. Power Syst. 2007, 22, 2058–2065. [Google Scholar] [CrossRef]
  29. Rodriguez, C.P.; Anders, G.J. Energy price forecasting in the Ontario competitive power system market. IEEE Trans. Power Syst. 2004, 19, 366–374. [Google Scholar] [CrossRef]
  30. Zhan, C.J.; Wu, X.G.; Kromlidis, S.; Ramachandaramurthy, V.K.; Barnes, M.; Jenkins, N.; Ruddell, A.J. Two electrical models of the lead-acid battery used in a dynamic voltage restorer. IEE Proc. Gener. Transm. Distrib. 2003, 150, 175–182. [Google Scholar] [CrossRef]
  31. Cervone, A.; Carbone, G.; Santini, E.; Teodori, S. Optimization of the battery size for PV systems under regulatory rules using a markov-chains approach. Renew. Energy 2016, 85, 657–665. [Google Scholar] [CrossRef]
  32. Pacheco, J.E.; Reilly, H.E.; Kolb, G.J.; Tyner, C.E. Summary of the Solar Two Test and Evaluation Program; Sandia National Laboratories: Albuquerque, NM, USA, 2000. [Google Scholar]
  33. Stiesdal Storage Technologies GridScale Battery. Available online: https://www.stiesdal.com/material/2019/02/Stiesdal-GridScale-01.01.19.pdf (accessed on 1 May 2020).
  34. Goldberg, D.E.; Holland, J.H. Genetic algorithms and machine learning. Mach. Learn. 1988, 3, 95–99. [Google Scholar] [CrossRef]
  35. Riechmann, T. Genetic Algorithms and Economic Evolution; Physica Fachbereich Wirtschaftswiss, University of Heidelberg: Heidelberg, Germany, 2002. [Google Scholar]
  36. Lou, J.; Bhobe, A.; Shu, Y.; Yu, J. Analytical calculation of transformer parameters by S-parameters, Solving Grid Connection Issue by Use of Thermal Energy Storage Technology. In Proceedings of the IEEE International Symposium on Electromagnetic Compatibility and IEEE Asia-Pacific Symposium on Electromagnetic Compatibility, (EMC/APEMC), Singapore, 14–18 May 2018; IEEE: Piscataway Township, NJ, USA, 2018; pp. 1310–1313. [Google Scholar] [CrossRef]
  37. Diaz-Dorado, E.; Carrillo, C.; Cidras, J. Coordinated reactive compensation in a Wind Park. In Proceedings of the 9th International Conference on Electrical Power Quality and Utilisation, Barcelona, Spain, 9–11 October 2007; pp. 1–6. [Google Scholar]
  38. National Grid Electricity Transmission. 2015. Available online: https://www.nationalgrid.com/uk/electricity-transmission/ (accessed on 1 September 2020).
  39. Behabtu, H.A.; Messagie, M.; Coosemans, T.; Berecibar, M.; Fante, K.A.; Kebede, A.A.; Van Mierlo, J. A review of energy storage technologies’ application potentials in renewable energy sources grid integration. Sustainability 2020, 12, 10511. [Google Scholar] [CrossRef]
  40. RTO Unit Commitment Test System; Staff Report; Federal Energy Regulatory Commission: Washington, DC, USA, 2020.
Figure 1. Wind park schematic diagram under study. BESS—Battery Energy Storage System; TESS—Thermal Energy Storage System; TSO—Transmission System Operator.
Figure 1. Wind park schematic diagram under study. BESS—Battery Energy Storage System; TESS—Thermal Energy Storage System; TSO—Transmission System Operator.
Energies 14 04718 g001
Figure 2. Schematic model of forecasting algorithm.
Figure 2. Schematic model of forecasting algorithm.
Energies 14 04718 g002
Figure 3. Flowchart showing the proposed method.
Figure 3. Flowchart showing the proposed method.
Energies 14 04718 g003
Figure 4. DIgSILENT model of Wind Park POC: Point of Conection, PCC: Point of common coupling.
Figure 4. DIgSILENT model of Wind Park POC: Point of Conection, PCC: Point of common coupling.
Energies 14 04718 g004
Figure 5. Simulated real power production at the PCC for a March month.
Figure 5. Simulated real power production at the PCC for a March month.
Energies 14 04718 g005
Figure 6. Simulated real power production and losses versus wind sweep speed.
Figure 6. Simulated real power production and losses versus wind sweep speed.
Energies 14 04718 g006
Figure 7. Reactive power capability of wind park without PPC.
Figure 7. Reactive power capability of wind park without PPC.
Energies 14 04718 g007
Figure 8. Modified reactive power capability of wind park after PPC set points.
Figure 8. Modified reactive power capability of wind park after PPC set points.
Energies 14 04718 g008
Figure 9. Forecast results of electricity prices.
Figure 9. Forecast results of electricity prices.
Energies 14 04718 g009
Figure 10. ANN results for a day during various seasons.
Figure 10. ANN results for a day during various seasons.
Energies 14 04718 g010
Figure 11. Selection and fitness functions from GA for Spring (26 April 2019).
Figure 11. Selection and fitness functions from GA for Spring (26 April 2019).
Energies 14 04718 g011
Figure 12. Selection and fitness functions from GA for summer (26 July 2019).
Figure 12. Selection and fitness functions from GA for summer (26 July 2019).
Energies 14 04718 g012
Figure 13. Selection and fitness functions from GA for fall (10 September 2019).
Figure 13. Selection and fitness functions from GA for fall (10 September 2019).
Energies 14 04718 g013
Figure 14. Selection and fitness functions from GA for winter (16 December 2019).
Figure 14. Selection and fitness functions from GA for winter (16 December 2019).
Energies 14 04718 g014
Figure 15. SoC/SoE results for both (a) BESS and (b) TESS.
Figure 15. SoC/SoE results for both (a) BESS and (b) TESS.
Energies 14 04718 g015
Figure 16. Charge (a), Discharge (b) profile for BESS.
Figure 16. Charge (a), Discharge (b) profile for BESS.
Energies 14 04718 g016
Figure 17. Charge (a), Discharge (b) profile for TESS.
Figure 17. Charge (a), Discharge (b) profile for TESS.
Energies 14 04718 g017
Table 1. Simulation results from the DIgSILENT model without PPC set points.
Table 1. Simulation results from the DIgSILENT model without PPC set points.
At Maximum Active PowerRequirement (Mvar)Capability (Mvar)Mismatch (Mvar)
Underexcited−58.99−41.5917.41
Overexcited58.9941.3717.62
Table 2. Simulation results from the DIgSILENT model after correction using PPC set points.
Table 2. Simulation results from the DIgSILENT model after correction using PPC set points.
At Maximum Active PowerRequirement (Mvar)Capability (Mvar)Margin (Mvar)
Underexcited−58.99−82.5823.5941
Overexcited58.9983.1124.15
Table 3. Simulated results for PQ-analysis before/after PPC operation.
Table 3. Simulated results for PQ-analysis before/after PPC operation.
DatasetSimulated Results
Before PPC Operation (Case 1)After PPC Operation (Case 2)
Max. cable loading53.315%55.022%
Max. transformer loading100.00%99.974%
Max. voltage inside power plant1.022 p.u.1.095 p.u.
Min. voltage inside power plant0.952 p.u.1.023 p.u.
Max. voltage at generation unit terminals1.001 p.u.1.071 p.u.
Min. voltage at generation unit terminals0.960 p.u.1.034 p.u.
Max. voltage at shunt terminal0.953 p.u.1.024 p.u.
Min. voltage at shunt terminal0.952 p.u.1.023 p.u.
Table 4. Parameters of storage systems.
Table 4. Parameters of storage systems.
ParametersBESSTESS
Power (MW)197197
Energy (MWh)250250
SoC limits (%)20–900–95
Round trip efficiency (%)75–9065–75
Ramp rateNo delay3 MW/min [40]
BESS capital Cost116 $/kWh100 $/100l
BESS maintenance Cost29 $/kWh/year1 $/100l
α 0.065 $/kWh0.065 $/kWh
β 0.0116 $/kW/h
γ 0.0099 $/kWh/h
Table 5. Comparison of proposed ANN results with literature.
Table 5. Comparison of proposed ANN results with literature.
ReferencesMethodsMAPE for a Day
[28]Selected 8 days for year 20107.81% (Average value)
[29]Selected 4, 8, and 12 neurons for hidden layer23.3% (Best value)
Proposed methodSelected 10 days for year 2019 with 10 hidden neurons6.89% (Average value)
Table 6. GA results: Profit obtained for cases without and with ESS.
Table 6. GA results: Profit obtained for cases without and with ESS.
Without ESS (Million $)BESS (Million $)TESS (Million $)
14.242.538.6
Table 7. Wind power curtailment without and with ESS.
Table 7. Wind power curtailment without and with ESS.
Without ESS (GW)BESS (GW)TESS (GW)
79.9443.5860.55
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ponnaganti, P.; Bak-Jensen, B.; Wæhrens, B.V.; Asmussen, J. Assessment of Energy Arbitrage Using Energy Storage Systems: A Wind Park’s Perspective. Energies 2021, 14, 4718. https://doi.org/10.3390/en14164718

AMA Style

Ponnaganti P, Bak-Jensen B, Wæhrens BV, Asmussen J. Assessment of Energy Arbitrage Using Energy Storage Systems: A Wind Park’s Perspective. Energies. 2021; 14(16):4718. https://doi.org/10.3390/en14164718

Chicago/Turabian Style

Ponnaganti, Pavani, Birgitte Bak-Jensen, Brian Vejrum Wæhrens, and Jesper Asmussen. 2021. "Assessment of Energy Arbitrage Using Energy Storage Systems: A Wind Park’s Perspective" Energies 14, no. 16: 4718. https://doi.org/10.3390/en14164718

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop