Abstract

Kenya is experiencing a fast increase in grid-connected intermittent renewable energy sources (RESs) to meet its increased power demand, and at the same time be able to fulfill its Paris Agreement obligations of abating greenhouse gas emissions. For instance, Kenya has 102 MW of grid-tied solar power and 410 MW of grid-tied wind power. However, these sources are very intermittent with low predictability. Thus, after their installation and integration into the grid, they impose a new challenge for the secure, reliable, and economic operation of the system. To mitigate these and to ensure proper planning of the system operations, accurate and faster prediction of the generation output of the wind energy resources and optimal design and sizing of storage for the large-scale wind energy integration into the grid are of paramount importance. Artificial intelligence (AI) and metaheuristic techniques have proven to be efficient and robust in offering solutions to complex nonlinear prediction and optimization problems. Therefore, this study aims to utilize backpropagation neural network (BPNN) algorithm to conduct hourly prediction of the generation output of Lake Turkana Wind Power Plant (LTWPP), a 310 MW plant connected to the Kenyan power grid, and optimally size its battery energy storage system (BESS) using genetic algorithm (GA) to guarantee its dispatchability. The historical weather data, namely wind speed, ambient temperature, relative humidity, wind direction, and generation output from LTWPP, are employed in the training, testing, and validation of the neural network. LTWPP and BESS are modelled in MATLAB R2016a software. Thereafter, the developed BPNN and GA algorithms are applied to the modelled systems to predict the wind output and optimize the storage system, respectively. BESS optimization with neural prediction reduces the BESS capacity and investment costs by 59.82%, while the overall dispatchability of LTWPP is increased from 73.36% to 90.14%, hence enabling the farm to meet its allowable loss of power supply probability (LPSP) index of 0.1 while guaranteeing its dispatchability.

1. Introduction

The utilization of renewable energy sources (RESs) has experienced remarkable growth in electrical grids worldwide. This is attributed largely to critical factors, namely a fast and unprecedented surge in volatile fuel prices, the scarcity of available primary energy resources utilized in conventional power plants, and the need to abate greenhouse gas emissions.

Wind power constitutes the renewable energy technology (RET) that has exhibited swift growth in comparison to other RETs being investigated currently. The capacity of wind worldwide is far larger than the world’s total energy consumption and the potential for development is substantial [1]. Globally, total capacities of about 651 GW of wind have been installed, with a yearly production of 60.4 GW in 2019 alone [2]. By 2023, it is foreseen that RES will be able to supply above 70% of the increase in electricity generation worldwide, with wind and solar taking the lead positions [3].

In Kenya, for example, at the beginning of the last decade, wind energy constituted only 0.24% of the total share of electricity generation. This has grown tremendously to constitute 11.74% by the end of the same decade, with still many more planned or incepted wind energy projects such as 100 MW Aeolus Ngong’ Wind and 60 MW Aperture Green Ngong’, among many others mentioned in [1]. The least-cost generation simulation included wind as a candidate project in Kenya for system optimization [4]. Still, there is a small level of competency and expertise in utilizing wind energy for large-scale power production in Kenya. Nevertheless, interest and awareness are gradually growing. The latest large investment in wind power is the 310 MW Lake Turkana Wind farm, the largest in Africa. LTWPP was anticipated not only to reduce the overall cost of energy but to reduce the greenhouse gas emission, hence placing it as a competitive energy source.

However, these merits and development of wind power come with several challenges while integrating it with the existing grids. These include erratic wind flow patterns/profiles which might introduce harmonics into the grid, interfere with the voltage levels, affect grid stability, and bring about unit commitment and scheduling problems [57]. To take care of these shortcomings, AI techniques such as artificial neural network (ANN) play a critical role in the development of faster and cheaper predictions for short-term, intermediate, and long-term predictions.

Recently, several wind impact researches have been conducted in various parts of the world in an attempt to address those challenges. The outcomes of these researches are normally linked to different aspects of wind power, notably its randomly varying nature, generator technologies, and sizing of the storage of the wind energy resources (WERs). Very few studies of that nature have been carried out in Kenya despite the increasing WER development projects. In January 2021, Clir Renewables, a Canadian software company, signed a contract with LTWPP to optimize power generation at the power station with the overall objective of improving the efficiency of the power plant, improving power output (reducing unnecessary turbine downtime), and grid stability [8].

In addition, it is worth noting that very few studies have integrated both the wind energy prediction and optimal storage sizing. In [9], the simulation performed demonstrated the efficiency of predicting the RES using a high-order neural network trained with an extended Kalman filter (EKF) and optimally sizing the grid-connected storage system utilizing both the GA and clonal selection algorithms to reduce CO2 emissions. Thus, there is a need to conduct studies to efficiently and precisely analyze and predict the fluctuating weather patterns associated with such renewables and be able to perform optimal storage designs that will make the grid robust and reliable in its operations. This study has explored a different prediction algorithm, that is, BPNN trained with Levenberg–Marquardt (LM) optimization and optimally sizing its BESS using GA to guarantee the dispatchability of LTWPP.

Despite the fact that Kenya has the largest wind power plant in Africa and is experiencing a fast growth of intermittent renewable sources, it does not have any BESS installations. This study therefore attempts to solve an existing problem that at the moment there is no available documentation that shows there is previous research that has been done on the same. The problem is that of the intermittent power output of WER at LTWPP which requires to be minimized so that the plant output can be dispatchable just like the other plants (geothermal, heavy fuel oil, and hydros) connected in the same grid. This intermittency has caused serious problems to the spinning generators having to undergo excessive hunting. This hunting has led to extensive bearing damages to the other generators called upon to either increase or decrease their outputs during intermittency periods. By researching the LTWPP output characteristic and arriving at the approximate contribution of this to the Kenyan grid, the authors were able to perform optimal sizing and costing of BESS required to take care of this intermittency. This study, therefore, contributes the following: First, it provides a unique research opportunity to analyze the incorporation of a large BESS system into LTWPP to achieve guaranteed and improved dispatchability. Second, it puts forward an optimized sizing strategy based on BPNN prediction, for calculating the optimal size and cost of BESS which meets all constraints. Finally, the study validates the results obtained and demonstrates the efficiency of this strategy via simulations of different scenarios.

The rest of the article is organized as follows: Section 2 presents a background overview of the literature review and mathematical modelling of the system components. The methodology used and problem formulation are described in detail in Section 3. Section 4 entails the performance analysis of BESS, while Section 5 presents a brief description of the case study. Section 6 provides the simulation results obtained and the accompanying discussions. Finally, the last section gives the conclusion of the study.

2. Literature Review

2.1. Background Overview

A simple time series-based approach for wind power prediction was first developed in 1984 by Brown et al. [10] by utilizing the utility's power curve. Thereafter, systematic investigations have been conducted in the field of predicting wind power or speed produced by WER and this has immensely contributed to the development of numerous techniques as well as reliable and effective tools which have been used with different success rates in various wind farms as proposed and reviewed in [11]. As a result, these accurate prediction tools have been found to yield better unit commitment and scheduling plans, achieve desired control of wind turbines [12], minimize the overall costs of operations, and enhance reliability linked with the increased penetration of wind power into the existing electrical grids [1315]. The most widely used approaches have been summarized in Table 1. Among all the techniques used, AI techniques, especially neural networks, have emerged to be more accurate and reliable while contrasted with traditional statistical models [12].

Using historical data and ANN modelling, the authors in [17] presented a tool that would be used by operators of RES systems to realize improved management and monitoring of their systems. The development of the prognostic tool was able to give an adequate forecasts 8 hours ahead of the absolute minimum, maximum, and mean hourly wind power. For short-term and long-term wind farm prediction, the ANN prediction models developed in [18] using the data mining approach were found to outperform the rest of the models.

Using the input parameters such as generation hours, relative humidity, and mean wind speed, the neural network model developed in [19] offers a reliable indicator for the output power from wind farms. The authors in [20] used ANN technique to approximate the wind speeds at a particular location utilizing the speeds of wind at a strong correlation location among neighboring locations. A wide comparison based on time horizons, performance analysis, and statistical distribution of normalized errors is conducted between five different types of ANN models, ARMA models, and ANFIS models in [21]. ANN models showed better performance. Local dynamic recurrent neural network models employed in [22] for long-term prediction of wind power and wind speed outshined the static rivals with reference to the improvement gained and forecast errors over the persistent method. The authors in [23] predicted wind speed automatically employing a hybrid neural network approach, comprising a radial basis function neural network (RBFNN) and a self-organizing map (SOM). The results indicated that the suggested method yielded improved output with small error rates.

An in-depth comparison study of three kinds of classical ANN approaches—adaptive linear element (ADALINE), radial basis function (RBF), and feed forward back propagation (FFBP)—is conducted in [24]. Results indicated that different model structures, inputs, and learning rates had a direct effect on the prediction accuracy. An advanced study was performed by the authors in [25] who proposed Bayesian model averaging (MBA) technique to merge the forecasts received from diverse ANN models, namely RBF, BPNN, and ADALINE neural networks. Utilizing their posterior model probabilities, BMA approach weights the individual forecasts. The better performing predictions acquire higher weights than the worse. This approached proved to be effective as shown by the simulation results.

Other various promising approaches for enhancing the accuracy of the prediction models includes: hybridizing different forecasting or numerical weather prediction (NWP) models as discussed in [2632], filtering of systematic errors emanating from NWPs using Kalman filtering [33], proper selection of input parameters [34], or the combination of any of the abovementioned promising approaches [15, 35]. Nonetheless, hybrid techniques may not at all times outperform the single techniques for all forecasting time horizons investigated as shown by the studies in [27].

Due to its ability to construct complex nonlinear relationships via training, and its ease of construction, ANN has been found to be generally a good choice for wind energy prediction. Most forecasting studies using ANN have offered the best performance compared to other techniques [12]. A more cost-effective option may be to utilize some optimal storage at the wind farm to mitigate the abovementioned problems as depicted by many ongoing studies and projects around the globe [36, 37]. Many researches have employed metaheuristic techniques such as GA, PSO, ABC, or a hybrid combination of AI and evolutionary computation such as ANN-PSO to optimally design and operate RES-based hybrid energy systems.

For instance, in [38], the optimization of an off-grid power system that consisted of the PV, WTG, diesel generator, and energy storage system (ESS) was investigated using GA. The proposed algorithm was found effective to aid the distribution network operators to reduce the total system cost that was related to the operation of a microgrid system. Authors in [39] performed simulations employing GA and the rule-based approach to optimize the cost of operation, while in [40], the authors used the firefly algorithm and considered the BESS’s depth of discharge while modelling the BESS operation cost. The cost of operation was well reduced as indicated by the results.

Various strategies have been suggested for obtaining the optimal size of BESS. For instance, a life cycle planning methodology of BESS in an off-grid was put forward in [41]. Using decomposition coordination algorithm, the optimal allocation of distribution energy resource (DER) capacity was carried out to match the demand growth while considering dynamic factors, namely demand growth and components’ uncertainties. Authors in [42] used an incremental cost approach to obtain the optimal values of BESS to realize the least running cost for an islanded DC microgrid, while in [43], the authors employed the grasshopper optimization algorithm to size the BESS while considering the energy cost and power supply probability efficiency.

In [44], the study found out that the grid-connected BESS was cheaper to operate than the stand-alone one, hence deserving consideration for efficient dispatch of wind power. Owing to the uncertainty of solar PV and wind, [45] asserted that capacity sizing was vital to completely meet the load demand. They formulated PSO algorithm to determine the optimal size for hybrid wind-PV with BESS while considering the unreliability in solar and wind energy generation.

To improve reliability, reduce cost, and to determine the optimal sizes of RES and BESS for a grid-connected microgrid system, the study in [46] used two constraint-based iterative search algorithms. Authors in [47] used ANN to validate the predicted load model and utilized GA to solve a chance-constrained model that handled the variability linked to RES. The results from the study showed that the storage system plays an essential role in renewable energy integration and in the reduction of environmental pollution.

From the background overview, it could be noted that the integration of RES into the existing large grids is very current and is gaining momentum globally, Kenya included. When a sector is new, there is always room for further studies. WER prediction offers one solution of efficiently integrating intermittent sources into the power grids while BESS acts as another solution as highlighted in the current studies reviewed in this section. Unlike these studies which have majorly employed these solutions independently, this study is unique in that it has employed the strength of both solutions to yield a better and improved way of guaranteeing dispatchability of LTWPP whereby the predicted output of the wind farm is utilized to optimize the BESS.

2.2. Mathematical Modelling
2.2.1. The Wind Turbine Generator (WTG) Model

The power output of WTG is estimated [38] according to the following equation:where is the density of air ( 1.225 kg/m3); V is the free wind speed (m/s); A is the area swept by the rotor (m2); is the power coefficient (dimensionless), are generator efficiency (%) and gear or bearing efficiency (%), respectively.

From the WTG profile in Figure 1, the wind power is given by the following equation [44]:where represents the wind power produced at time t, represents the wind speed at any hour t, represents the cut-in wind speed, represents the rated wind speed, and represents the cut-out wind speed.

2.2.2. BESS Modelling

Storing the surplus energy available, , produced by the WER is the first requirement that should be addressed in optimal BESS sizing [48] as shown in the following equation:where represents the demanded power at any given time t, represents the predicted output power from WTG at any given time t. The positive value of indicates that the wind farm produces surplus power, while the negative sign indicates a deficit in demand that the wind farm cannot meet.

The strategy is first to store the available energy from the WTG during hours of higher wind output and discharge the BESS when the output power is less than the demand until the BESS reaches a given minimum state of charge (SOC). The variable SOC(t) at any particular hour is the amount of energy stored in the BESS at that particular instant and is formulated as follows:where is the SOC at time , represents the charge or discharge power at time t, represents the BESS capacity, and represents the incremental time for the optimization.

The BESS in this study is modelled based on the SOC and depth of discharge (DOD) limits. Increase in DOD increases the cost of BESS due to more power discharged by the BESS. Cost is a very key factor in the development of BESS. Thus, a cheap BESS would be the most economically feasible. The minimum SOC is given as follows:

The charging and discharging energies at any given time are given in equations (6) and (7), respectively [49].where and are the discharge and charge efficacy of the BESS (%), respectively, while and are the discharged and charged powers, respectively.

The following assumptions are made in the optimal sizing of BESS:(i)Temperature affects the lifespan of lithium BESS. LTWPP site has a mean temperature 28°C, hence a temperature correction factor of 0.964 will be considered [50].(ii)As reported in [51], “sudden deaths” of batteries occur when its nominal capacity reaches 80% of its initial value and this usually marks the end of its life. It is for this reason that the aging correction factor of 20% has been selected since the aging process is accelerated after this point.

3. Methodology

3.1. Wind Data Analysis and Preprocessing of the Parameters in LTWPP Farm

First, the data which were obtained from Meteoblue [52] were preprocessed before being applied as a forecast model input. This was for eliminating any sharp peaks or errors in the data and to ascertain that there were not any missing data points. The wind speeds at a hub height of 80 m were then converted to the LTWPP WTG hub height of 44 m using the following expression [53]:where denotes the speed of wind at height , denotes the average speed of wind at , and denotes the friction coefficient or the Hellman’s exponent and is always dependent on the topography at a given site. For our case, α was taken as 1/7 since it is an open land.

3.2. Wind Power Output Prediction Using BPNN

BPNN has two stages in its learning algorithm. First, the training input parameters are fed to the network source nodes and are propagated from one hidden layer to another until the output parameter is produced at the output layer. Second, if the value of this parameter deviates from the anticipated output, an error is evaluated, which is propagated backwards as the network weights are modified. The common structure of BPNN is depicted in Figure 2.

The demerit of utilizing BPNN is its long training duration with very many iterations. To tackle this problem, Levenberg–Marquardt (LM) optimization was employed to train BPNN. Training was performed in batch model. First, the data were normalized. The database was then constructed from the historical real data collected from the LTWPP since it commenced its operation (October 2018 to September 2021).

Data for the first 2 years were used to train and validate the network while the remaining ones for a year, which was used to test the network. The optimum architecture of the BPNN, that is the choice of the optimal number of hidden layers, was performed using a combination of trial-and-error method and the proposed method in [24], where the optimal number of hidden layers (h) is given by the whole number near to log (T), where T denotes the total number of vectors used for training. The training, validation, and testing process of BPNN is illustrated in Figure 3.

The proposed optimal sizing strategy flowchart based on BPNN prediction for evaluating the BESS size and cost is shown in Figure 4.

3.3. Optimal Sizing of Battery Energy Storage System
3.3.1. Selection of BESS

Batteries have so far demonstrated to be an economically feasible energy storage technology. Previously, low round trip efficiency and high costs hampered the mass deployment of BESS. The high efficiency and reduction in cost of lithium-ion batteries led to the surge of BESS usage in the last few years for both large-scale grid-level deployments and small-scale behind-the-meter installations [54]. Table 2 shows the distinguishing characteristics of the various battery technologies currently in the market. It can be clearly seen that Li-ion ranks the best due to the superior characteristics it possesses in comparison to other battery types. The critical review done by [55] indicates that Li-ion batteries as the most suitable and with high potential for RES integration into the power grid. However, their cost still needs to be minimized if they are to be embraced fully. Therefore, this study proposes to use lithium-ion batteries.

3.3.2. BESS Sizing

The daily BESS average energy need, (MWh/day), depicted below is obtained from , evaluated from equation (3).

Based on the considerations and assumptions made in Section 2, the required BESS capacity (MWh) is given by [50];where is the correction factor for the effects of battery aging and temperature degradations and D represents the days of autonomy.

By applying GA, the BESS size obtained in (10) is further minimized with used as the maximum boundary limit. GA is a heuristic evolutionary algorithm used for hybrid search and optimization problems. It mimics the Charles Darwin’s theory of natural selection [56].

Two of the most notable merits of GA over other optimization algorithms are: parallelism, which helps in false peak reduction [47], and its capability in handling complex problems. It can effectively handle different kinds of optimization, whether the fitness function is linear or nonlinear, or continuous or discontinuous. Nevertheless, GA has few drawbacks. For the fitness function formulation, the choice of critical parameters, namely, crossover and mutation rates, and the selection criteria of the new population must be done carefully, otherwise it will be difficult for the algorithm to converge. GA still remains one of the most widely used techniques in modern nonlinear optimization despite the abovementioned disadvantages [57]. The GA optimization algorithm flowchart is shown in Figure 5.

3.4. Problem Formulation

Stable, efficient, and economical operation of the grid incorporated with WER is largely dependent on the accurate forecasting of the wind power output and optimally sized storage to take care of the uncertainty caused by the variable and intermittent nature of wind. There is a need to optimize WER when it is available. To improve the prediction accuracy of our system, other parameters that affect the wind speed and the average power production by the WTG such as relative humidity, wind direction, and temperature are added as inputs to the backpropagation neural network (BPNN) with Levenberg–Marquardt (LM) optimization technique. The wind power output from the forecast is used to calculate the excess surplus power in equation (3), which is the first key requirement in finding the optimal BESS size.

3.4.1. The Objective Functions

The proposed study aims at sizing an optimal BESS by minimizing the total cost of the lithium BESS investment, when incorporated into grid-connected LTWPP farm, the following equation represents an objective function (OF) to be minimized;where is the storage capacity of BESS in MWh, is the per-unit cost of the lithium BESS(US$/MWh), is the capital cost of the lithium BESS(US$/MWh), is the number of BESS replacements, and represents the operation and maintenance cost (US$/MWh). The capital cost consists of the buying, power conversion, and installation costs.

CRF is given equation (13) and represents the capital recovery factor. It enables the calculation of annual payments, which are spread equally over a stipulated time based on the initial payment where n is the BESS lifespan and i is the discounted interest rate.

3.4.2. Constraints

The objective function given by equation (11) is subject to:

(1) Balance Constraint. The total demanded power by the users must be equal to the total produced power by the WER and BESS, as depicted in the following equation:

(2) BESS Constraints. The power for charging and discharging BESS must be restricted according to the following equations:

Moreover, the energy stored must be within and at all times as depicted in the following equation:where and represents the lowest and highest BESS energy rating, respectively.

The BESS operates within and as shown in the following equation:

(3) WTG Power Limits. The forecasted wind power output is also limited by its rated power as presented in the following equation:(4) Reliability Constraint.with representing the maximum allowed LPSP.

4. The Case Study

LTWPP, the largest wind farm in Africa, is located in Sarima (2.49°N 36.8°E), between the foot slopes of Mt Kulal and the southeastern end of Lake Turkana, in Loiyangalani District, Marsabit County, Kenya. The area, often labelled as “The windiest place on Earth,” has distinct geographical conditions whereby the variations of daily temperature produce strong wind streams between the desert hinterland (with steep temperature variations) and Lake Turkana (with relatively constant temperature). It comprises of 365 Vestas V52, doubly fed asynchronous WTGs, each with a capacity of 850 kW, which generate a total of 310.25 MW of electricity. The technical specifications of the WTGs are given in Table 3. The power generated from the WTGs is transmitted to the wind farm internal overhead electric wires (collection grid) with the voltage of 33 kV via 0.69 kV/33 kV transformers to a sectionalized 33 kV substation that is located within the premises. Through a 33/220 (400) kV substation, the wind farm is AC-connected to the grid via a 400 kV high-voltage double transmission line [60]. Figure 6 shows the topology of the LTWPP incorporated with BESS.

5. Performance Analysis of Battery Storage

5.1. Reliability Model

In sizing BESS, the power reliability analysis is very crucial. When the energy generated by WTG plus the energy stored in the BESS fails to meet the load demand at hour t, the system experiences loss of power supply (LPS) [49] and is expressed in equation (22) as;where is the efficiency of the inverter.

Using LPS, the loss of power supply probability (LPSP) is evaluated as follows [9]:

The key aim in sizing BESS is to ensure that which is 0.1 in this study.

5.2. Power Shed

When , it necessitates load shedding in order for the RES without storage to match the dispatch curve. Alternative sources of energy such as diesel generators are usually assumed to cater for this deficit power. Thus, the revenue loss due to power shed, , is given by the following equation [50]:where is the cost of diesel generators and is estimated to be US$ 0.27/KWh of electricity.

6. Simulation Results and Discussion

Figure 7 is a representation of the histogram of various wind speeds at LTWPP from 01/09/2018 to 15/09/2021. The average value of the collected hourly wind speeds was between 10.5 and 11.5 m/s. Weibull parameters for the above data were found to be scale factor, c = 10.9096 and shape factor, k = 4.1525. The wind rose of LTWPP with 16 directions is represented in Figure 8. It can be deduced that the dominant wind which averages at 11 m/s blows from ESE direction. The data obtained from Meteoblue [52] were found to be good for usage since there were no missing data for the period under study as shown in the relative humidity, ambient temperature, wind speed data, and the calculated output power in Figure 9. In addition, the data were within the normal expected range and the overall pattern and mean are the same as the ones found in the literature [61].

The BPNN training performance results are tabulated in Table 4. The best performance MAE, MSE, and RMSE were found to be , respectively.

The regression and performance plots for the BPNN training are depicted in Figures 10 and 11(a), respectively. Results on Figure 11(b) indicates the effectiveness of the proposed algorithm in forecasting the output wind power of LTWPP since the predicted power output curve is almost the same as the actual power output curve. This neural forecasted output is used in the optimal sizing of the BESS.

The average per hour daily load demand to be supplied by LTWPP is given in Figure 12, while the hourly data for the 12 months used in this study is shown in Figure 13 for the period between October 2018 and September 2019. As it stands now, there is no specific amount of demand that is allocated to LTWPP. Instead, it is supposed to supply all power that it generates at any time. However, this poses a challenge to the system in terms of economical dispatch scheduling of other sources since either the geothermal or diesel generators have to be in standby mode to cater for any shortage unmet by LTWPP. In addition, the sudden ramping up or down of these standby sources causes mechanical stresses, hence reducing their longevity.

Therefore, with the incorporation of BESS to LTWPP, it can be able to supply a specific allocated amount to the grid. Any fluctuations of the wind are taken care of by the BESS, hence guaranteeing its dispatchability. This study assigns LTWPP to supply 17% of the total demand on the grid. This arises from the fact that LTWPP contributes an average of approximately 17% of the total demand on the Kenyan grid.

The variations in the predicted supply and demand allocated to LTWPP for a period of 1 year and for 1  day (11 April, 2019) are depicted in Figure 14. It is from these curves that the excess/deficit power that needs to be stored/discharged in the BESS is obtained. From the daily curve, it could be clearly seen that the power demanded in the very early hours of the day when people are still asleep is less than that produced by the LTWPP. However, when it is most needed during the day, especially between 0900 and 1600 hrs, enough is not produced. Thus, from the curves, we can deduce that without storage, we cannot guarantee the dispatchability of the power supplied to meet the load demanded. The system revenue losses resulting from power shedding amount to a sum of US$ 67,179,299.73 annually with a LPSP ratio of 17.44%.

The size and total lithium BESS cost required are evaluated using equation (10), and the results are tabulated in Table 5. By applying GA, the BESS size obtained in (10) is further minimized with used as the maximum boundary limit. The GA parameters that provide the optimum performance are depicted in Table 6. Table 5 also gives the results obtained after optimization.

Optimization of the BESS yielded remarkable changes in the investment costs of the system. The size is reduced to 19.65% of the BESS size evaluated using equation (10). The BESS size thus evaluated is approximately 14.31% of the rating of the LTWPP. These results, based on the proposed methodology, are found to be consistent and, in some instances, better compared to the ones got from previous studies. In a feasibility study conducted in [62], the authors suggested a BESS capacity of approximately 24% of the capacity of PV-wind system. The BESS size obtained utilizing optimization was 6–10% according to [63], while in [64] the battery size of 15–25% of hybrid RES size is proposed to be adequate for dispatching power effectively. And according to the study carried out in [65] using fuzzy and neural controllers, the BESS size was estimated to be 30–34%. Hence, a BESS of smaller size about 14.31% is proved to be sufficient in this current study.

Gains realized and energy dispatched before and after incorporating BESS with the LTWPP are contrasted. BESS aided to improve the power dispatched by 14.47%, therefore making it possible for the LTWPP to meet the planned dispatch curve for most of the times, thus decreasing the LPSP to 0.0986, which is below the reliability index, LPSP of 0.1. Losses due to power spillage are reduced considerably and hence, the energy saved is used to meet the demand when wind is unavailable. It is worth noting that the power delivered matches the load curve after the incorporation of BESS as shown in Figure 15. Moreover, the feed-in tariff rate for wind power in Kenya is fixed (at US$ 0.11). Thus, any gains earned are as a result of the reduction in power losses through spillage and shedding.

For the period under simulation, the total power discharged to the grid from the BESS is 73406.16 MW. This presents a contribution of 20.52% of the total power supplied by the LTWPP-BESS system. Plots in Figure 16 show the variation of SOC of the BESS with time for a period of 2.5 days. Figure 17 depicts the charging and discharging BESS power. During charging, the BESS power is positive while it becomes negative during discharging. The SOC varies within 98–19.6% limits. This demonstrates the effective usage and performance of the BESS.

The lifespan of any BESS largely depends upon the number of discharge/charge cycles it performs. This was determined to be approximately 300–350 cycles/year. Consequently, the lifetime of lithium-ion BESS in this study is calculated to be approximately 10 years since the total life cycles of the lithium-ion BESS is 3500 cycles. This result shows a better lifetime operation compared to those performed elsewhere. For instance, according to [66], the battery life is estimated to be 2–10 years using a battery dispatch strategy developed to maximize its lifetime. Therefore, this sizing strategy yields an optimal usage of BESS that ensures better dispatchability of LTWPP. To further demonstrate how effective the proposed methodology is, a contrast based on investment cost among diverse possible solutions is shown in Table 7.

Case 1: (a) and (b) represent LTWPP with no BESS. The second cases represent LTWPP with no forecasting but with BESS. Case 3 represents neural forecasted LTWPP with unoptimized BESS while the final cases represent neural forecasted LTWPP with BESS optimization using both GA and PSO. For the cases of LTWPP with no forecasting, the assumption made was to take the previous year data (October 2017 to September 2018) and use it as the estimate for the period October 2018 to September 2019. From the results in Table 3, various analyses were conducted as illustrated in Figures 1822.

It can be noted that BESS optimization is very critical as depicted in Figure 18. Optimization reduces the BESS capacity for non-forecasted LTWPP and the neural predicted LTWPP by 65.28% and 80.35%, respectively. This has a direct effect on the investment cost as illustrated in Figure 21. The BESS calculated without optimization and without any prior forecasting of wind power exhibits the highest investment cost, while the neural forecasted LTWPP with BESS optimization exhibits the lowest cost.

Wind energy prediction also plays a very critical role. First, the dispatchability of LTWPP to meet its scheduled demand is found to be generally low for all cases without wind power forecasting as shown in Figure 20. In addition, for Case 1, the neural prediction amounted to a reduction of revenue losses incurred due to load shedding by 34.55% and improved on the dispatchability of the wind farm by 9.2%, without considering any storage as illustrated in Figures 19 and 20, respectively. Furthermore, comparing Case 2a with Cases 4a and 4b, it is found that BESS optimization with neural prediction reduces the BESS capacity and investment costs by 59.82%, while improving the dispatchability of LTWPP to 90.14% from 72.97% (Figure 21).

As demonstrated in Figure 22, neural forecasting and BESS incorporation has a direct effect on the dispatchability of LTWPP. It could be seen that the highest dispatchability (94.06%) occurs when we have the prediction of wind together with the large unoptimized BESS capacity. However, the investment cost for this case is very high. Therefore, the optimum dispatch (90.14%), although lower, meets the allowable reliability index of 0.1 and realizes a well-minimized investment cost. This way, by managing the surplus energy through optimizing its storage, the dispatchability of the neural predicted WTG system is improved by 16.78%.

7. Conclusion

Kenya is set for a remarkable growth in the renewable energy sector. This study has suggested an optimized sizing strategy based on BPNN prediction. To effectively minimize the cost of the BESS, the optimization uses a genetic algorithm. Based on different scenarios, performance metrics were analyzed to explore the viability of the LTWPP-BESS system, a grid-connected RES. Due to uncertainties in wind power, BPNN prediction employed improved reliability of the system by improving its dispatchability by 9.72%. A better performance, that is dispatchability improvement of 16.78%, was obtained by utilizing both the predicted output of LTWPP and optimized storage. Moreover, this study found that BESS optimization with neural prediction reduces the BESS capacity and investment costs by 59.82%. LPSP is reduced from 17.44% to 9.86%, while overall dispatchability of LTWPP is increased from 73.36% to 90.14%, hence enabling the farm to meet its LPSP index while guaranteeing dispatchability. Simulation results indicated the sizing methodology based on the BPNN forecasting and dispatch strategy to be effective and efficient since congruous and better results were achieved, in comparison to other previous studies and different scenarios investigated herein.

However, the study did not attain 100% dispatchability, implying that there is still room for further improvement. This could be explored further in future comparative studies using other AI and metaheuristic technique combinations such as ANN-fuzzy and PSO, ANFIS and hybrid ABC-PSO, among many others, which have not been covered in this study. Investigations on the impact of climate change on WER prediction is also of great interest. This proposed research did not include such and could also be included in advancing this study. Inclusion of other key factors that affect the BESS, such as temperature and capacity fade in the constraint, which were excluded in this study, could also be considered. Lastly, future studies should test the designed LTWPP-BESS system on the Kenyan power grid or the equivalent standard IEEE system to evaluate its impact on various power qualities and stability issues.

Data Availability

The data are available upon request from the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest with regard to the publication of this article.

Acknowledgments

The authors acknowledge the support provided by the African Union Commission in funding this work. They also thank Meteoblue for providing LTWPP site data that greatly aided in conducting this research.