1. Introduction
There is an increasing need for reliable technical information to inform or act as a reference for investors in hydropower development in tropical run-off rivers especially in developing countries where there is a scarcity of time series data on river flow rates and a lack of published feasibility studies on plants. Hydropower is a widely deployed renewable energy technology, reliable, and in most cases, considered an inexpensive renewable energy source, with a levelized cost of electricity (LCOE) ranging from 0.05–0.10 US
$/kWh [
1]. This is reinforced by prior analyses (Panwar et al.; 2011; Chu and Majumdar, 2012; Jiang et al., 2018, as cited in [
2] indicating that hydropower is an affordable, reliable and the most deployed renewable energy technology.
Renewable energy are energy sources which are refilled at a higher rate compared to the rate of consumption and produce lower emissions in comparison to fossil fuels. Examples of renewable energy are solar energy, wind energy, geothermal energy, hydropower, ocean energy and bioenergy [
3].
In Africa, hydropower is the most developed renewable energy source, with a total installed capacity of 53 GW or approximately 63% of the total deployed renewable energy. However, Africa needs additional generation as a catalyst for economic development, while simultaneously controlling greenhouse emissions. As of this writing, the continent has developed only 11% of its hydropower potential [
4].
Hydropower plants are classified based on different factors such as water inflow regulation available, size of installed plant capacity, and type of load being supplied. Regarding water inflow regulation, hydropower power plants are classified into three categories: (1) Reservoir plants, which have a dam and reservoir that stores water and regulates water inflow; (2) “Run-of-river” (RoR) plants, which have minimal reservoir capacity—typically a small diversion dam and catchment—and the water supply is mainly determined by the current river flow; and (3) pumped storage plants, which have both upper and lower reservoirs. During low-power demand, pumped storage uses electricity to pump water from the lower to the upper reservoir, and uses water from the upper reservoir for electricity generation during high demand [
1].
Hydropower power can also be classified as large or small depending on the plant’s generation capacity. Small hydropower plants have a generation capacity up to 10 MW but specific definitions vary between countries [
5]. There are three other commonly accepted classes of hydropower plants: mini (100–1000 kW), micro (10–100 kW) and pico (<10 kW). Plants of these sizes are sometimes classified independently or included in the small hydropower plant category.
In Rwanda, the total installed capacity, by rated power, for all energy sources is 238 MW, of which renewable energy accounts for 55%. Total capacity, in descending order is: 44% hydropower, 25% thermal (primarily marine diesel), 8% methane gas (harvested from Lake Kivu), 8% imported and shared resources with neighboring countries, 6% peat-fired thermal plants and 5% solar photovoltaic [
6]. Small-scale hydropower includes 44 plants [
6], all of which are RoR. The country is characterized by a tropical climate with dry season river flows that are 50% or less than wet season flows [
7]. Rwanda is commonly characterized as having four seasons: A short rainy season (from September to November), a long rainy season (between March to May), a short dry season (between December to February) and a long dry season (from June to August). Recent climatic data indicates that rainy seasons are becoming shorter but with higher intensity rainfall events [
8]. This change is likely to increase the variability in river flows [
9].
Sebeya River is located in the western part of Rwanda and is the largest catchment in the region. The Sebeya originates in high elevation mountains along the edge of Nile–Congo divide in the Rutsiro District. The river runs through Gishwati national park and discharges into lake Kivu with a length of 110 km [
10]. The whole region of East-Central Africa has experienced an average temperature increase of 0.29 ˚C from 1985 to 2015, and Rwanda in particular has experienced large fluctuations in annual rainfall over the 55-year period from 1961 to 2016 [
11].
The plants which are the focus of this study—small RoR hydropower plants—have minimal reservoir capacity and use immediate river flows for electricity generation [
1]. As a result, these plants are sensitive to seasonal weather variations that impact river flow. In particular, this study deals with sizing a hydropower plant in these highly variable flow conditions, proposing a methodology for analysis and suggests general conclusions on the optimal design flow rate and turbine type for cost-effective electricity generation. This analysis is of interest to future investors in the run-of-river hydropower plants in similar climatic conditions.
Hydro turbines are classified into two types: impulse and reaction turbines. In reaction turbines, water flows over runner blades and energy production results from the pressure caused by changes in direction of the moving water. Reaction turbine types include Francis, Kaplan and propeller turbines. In impulse turbines runner buckets redirect a jet of water impacting a runner’s curved buckets, creating a change in momentum and resulting force [
12]. Impulse turbine types include Pelton, crossflow and Turgo turbines.
In turbine selection, there is a set of requirements and specifications which include either site conditions (head and flow rate) and/or output power requirements [
13]. In general, turbine selection is based on the available head and flow conditions in comparison to standard charts during initial sizing, while detailed design includes additional analysis to consider technical, social, environmental, and economic factors [
14]. The work performed here is focused on the initial design phase, when a developer is primarily interested in the initial sizing to estimate the cost of civil works, permitting, land-access costs, and similar preparatory steps, prior to detailed design.
The previous literature identified several technical parameters considered in the sizing of hydropower plants and predicting energy generation. In the design process, sizing of a hydropower plant size is primarily determined by head, reservoir size, minimum downstream flow rate, and river seasonal inflow [
11]. A study by Fairuz et al. [
15] on factors affecting mini hydropower production efficiency concluded that humidity and rainfall have a significant effect on the power generation in mini hydropower plants and can be used to predict energy generation. A study by Singh et al. [
16] concluded that head, discharge, and generator parameters affect the operation of RoR hydropower plants. A study by Aggidis et al. [
17] on the costs of small-scale hydropower production concluded that the cost of manufacturing mini hydropower turbines is determined by the hydraulic characteristics of the hydro resource available and changes with turbine size and type.
It also defined the levelized cost of electricity (LCOE) and identified parameters considered in computing LCOE of different renewable energy technologies. LCOE is a measure used in estimating the cost of generating electricity over a lifetime of any generation technology. It can also be used to determine the cost-effectiveness of electricity generation of any technology without taking in account the assumptions of price at which generated electricity is sold to both the end user and the grid. Again, from the LCOE one can determine whether a technology is profitable or not. If the LCOE is lower than the price at which the generated electricity is sold, it is profitable and when the LCOE is greater than the price at which the generated electricity is sold it is unprofitable [
18]. Different studies have identified different important parameters to be considered in calculating the LCOE of different technologies. A study by Patro et al. [
19] on clean development projects in India found that in recent years many developers of small hydropower plants conduct life cycle costing studies without ignoring any costs which were used before operating and maintaining the plant, replacement, certified emission and others. A review of the LCOE of solar photovoltaics by Branker et al. [
20] on how to get the correct LCOE of solar photovoltaics concluded that due to the highly variable LCOE inputs of solar technology, to get correct LCOE outcome, almost all costs associated with solar electricity generation should be included. A study by Ouyang et al. [
21] on the LCOE of renewable energies and required subsidies in China found the accurate cost estimations of renewable energies were urgently needed to accelerate renewable energy development. The study recommended that the existing feed–in–tariff (FIT) be improved and adjusted based on the LCOE and to provide subsidies in the short-term to reduce the high-cost of renewable energy.
Additionally, previous work has identified that low-quality pre- and full-feasibility studies contributed to project failure [
22], and a study by Fairuz recommended the use of machine learning to improve predictions of energy generation [
15]. In Rwanda, there is no similar published work, and, unfortunately, there is little time series data for river flow (discharge rates). Therefore, to systematically size mini-hydropower plants for Rwandan rivers, such as the Sebeya, it is necessary to both estimate the river flows and to use those flows to estimate the optimal sizing of the power plant. Thus, the Sebeya flows were predicted using a neural network (scaled conjugate gradient or SGC method) by inputting both secondary data (using precipitation and flow rate) and primary data (collected flow rate at the Sebeya river).
2. Materials and Methods
2.1. Flow Rate Estimation
To size a small RoR hydropower plant for any river, a time series flow rate of the river is needed. However, the only available historical consistent flow rate of the Sebeya from a Rwanda Water Resource Board was from 2017. By 2019, historical flow rate data was inconsistent or missing for the Sebeya river. There the report from 2017 was the best method for predicting, given that the lack of data is a real issue for other sites. Machine learning is being used in different approaches to help forecast energy generation, and improve the operation, design, and maintenance of hydropower plants. Among other applications, it has been implemented in the optimal dispatch of hydropower plants based on head and river flow [
23], to assess hydropower reliability due to climate change [
24,
25], reservoir water level forecasts [
26], and to help mitigate silt erosion in the hydro-mechanical components of the plant [
27]. Yang et al. used a deep learning method to predict hydropower generation at small plants based on historical precipitation data provided by weather stations with 93% accuracy, [
28]. However, access to reliable and consistent data often remains a challenge in such studies [
29]. This study used the Sebeya 2019 precipitation from the Rwanda Meteorology Agency, available Sebeya flow rates from 2017 (provided by the Rwanda Water Resources Board) since there were no historical consistent flow rates of the Sebeya for 2019, and Sebeya field visit measurements conducted in 2019 (19 April 2019, 17 May 2019, 29 June 2019 and 20 August 2019). Using these data, machine learning was used to predict the Sebeya flow for the observed weather conditions in 2019, a timeframe where near-continuous meteorological data was available for the watershed area of the river. The methodology used is shown in
Figure 1.
The machine learning method estimated river flow data for 2019, using daily precipitation in the Sebeya river watershed during 2019 after being trained on the available Sebeya flow rate data from 1 May 2017 to 31 December 2017. The training data represented the longest period of near-continuous flow rate data available—200 days of data, spanning wet and dry seasons, with coincident records for both precipitation and river flow. In contrast, only precipitation was available for 2019, the year when the study was conducted. Therefore, data from 2017 was utilized to train a neural network (NNET), and the trained model took as its input precipitation, and predicted the river flow.
Training used the scaled conjugate gradient (SCG) method. SCG is a supervised learning algorithm that minimizes the error between its results and a target variable, using a defined error function [
30]. The advantage of this method in comparison with others, like Adam and scaled gradient descent, is that it does not need a learning rate and therefore is less dependent on user inputs, making it effective for large-scale problems [
31].
According to the US Geological Survey [
32] and US rivers, basin lag time—the difference between peak rainfall and peak discharge of the river—varies from 0–105 h. To accommodate this lag, the model uses the 10 previous days of precipitation, approximately two times the estimated basin lag time.
The Sebeya river flows are characterized by high flow rate events shortly after heavy rain in the watershed of the river. Therefore, special precautions were taken to assure that predicted river flow rates—in practice, 1000 Monte Carlo iterations of the model—reflected the nature and magnitude of these events.
First, error was calculated using mean absolute error (MAE), rather than root mean squared error, as MAE assigns less weight to the large errors. This allows the fitted model to better retain points that deviated substantially from the norm, i.e., to better model high rainfall events.
Figure 2 provides an example. The middle plot compares the training data to the predicted data for one iteration of the model. From index 25 to 80, the river flow is nearly constant as this is the dry season, with few rainfall events. Later indices indicate the start of the rainy season, characterized by periodic rain events and a higher variation in river flow rate.
The training and error functions were selected to retain the extreme flow events, which are observed on the river following major rainstorms. One example of a predicted river flow rate is shown for 2019 in
Figure 3.
From
Figure 2, during the period with high variability in the rainy season river flow (index 80 onwards), the model output does not capture high variations in the river flow seen in the training data, i.e., from actual flow measurements. This is characteristic of machine learning models for this type of application. To add this high variability back into the model results, the predicted flow was adjusted by adding in the deviation between the neural network (NNET) prediction and the training data observations of the actual flow. The method used was as follows: NNET outputs were binned at 0.1 m
3/s intervals; an error distribution between the NNET prediction and the training data was calculated for each bin; the NNET prediction for 2019 was adjusted by randomly adding a deviation drawn from the error distribution of the correspondent NNET bin. Adding these adjustments further restores extreme river flow events to the smoothed NNET results.
Each run of the trained NNET results in a different prediction of the river flow rate. This process was performed 1000 times, creating a family of 1000 plausible daily average river flow rate predictions for the Sebeya in 2019, see
Figure 4.
As noted earlier, limited river flow rate data was available, and the data available had substantial periods with no, or suspect, flow readings. Therefore, to provide a limited ‘spot check’ of the model, the first author measured the river flow at a power plant located on the river on four dates: 19 April 2019, 17 May 2019, 29 June 2019, and 20 August 2019, resulting in measurements of 4.24 m
3/s, 3.3 m
3/s, 1.96 m
3/s, 1.61 m
3/s, respectively. A simple measurement process was performed: Timing the movement of a float on the surface of the river, measuring the cross-sectional depth of the river, and applying appropriate hydraulic calculations. We compared these data to the distribution of results for the measurement days, as shown in
Figure 5.
Two measurements were well within the inner quartile of the predicted flows, while the other two measurements were within the outer quartile and within the range of the outlier data, respectively. It should be noted that Rwanda was in drought conditions in 2017, while 2019 was considered a rainy year. While limited in scope, these measurements indicate that the model produces values, including outliers, which are representative of the observed conditions.
2.2. Power Production Estimation
The analysis utilized here develops a first-order simulation of the turbine performance, using the 1000 Monte Carlo iterations predicting the river flow rate.
Turbine produced power, is:
where
is the acceleration of gravity (m
3/s),
is the flow rate of the river in m
3/s used to select the turbine and design the power plant (‘designed flow rate’),
is the working head of the river net of friction losses in the penstock and head gates
,
is the efficiency of the turbine at the designed flow rate, and
is the water density
. For this study, we assumed the plant was designed such that the efficiency is maximized at the designed flow rate, i.e.,
is the maximum efficiency of the turbine, typically at or near full load.
For the Sebeya River, the flow rate () varies between 1.3 m3/s and 5.5 m3/s under non-flood conditions. In this study all plants were designed for flows below 3.5 m3/s. Water above this flow rate would be directed around the turbine via spillway gates and not contribute to power production.
For conditions away from the designed flow rate, Equation (1) can be considered as a general function of the river flow rate and efficiency, since all other parameters remain fixed after the plant is designed. Further, efficiency is also a function of the deviation of the current flow rate from the designed flow rate, resulting in:
where
and
are arbitrary, typically discretized functions of flow rate. Flow rate
is determined by the analysis in the previous section for all times considered in the study, and, in general,
is the result of empirical testing of the selected turbine; see below. In practice the analysis utilizes normalized calculation throughputs, where performance of the hydropower plant is expressed as fractional power output calculated at the given flow rate, relative to the designed flow rate. Fractional flow rate is:
where
, is the fraction of designed flow, in percent,
is the designed flow rate, and
is the predicted flow rate of the river.
Turbine efficiency for each turbine type was estimated by digitizing a generic turbine efficiency chart [
33] that provided the efficiency at flow rates relative to the designed flow rate of the turbine. Data is valid from 5–100% of the designed flow rate, and digitization was conducted at an interval of 5%, see
Figure 6.
Energy produced per turbine type was constrained to 100% of the designed flow rate, resulting in:
where
is energy produced relative to the designed conditions during one time step
,
is the time step of the simulation—1 day in this analysis—and
and
are calculated for each time step. Using the predicted flow rate, turbine net efficiency for the simulated time period is:
where
is the net turbine efficiency relative to the power produced, across all Monte Carlo simulations, for a given turbine, sized for a selected designed flow rate. The number of periods used to simulate
,
, is the product of the length of the simulated river flow (262 days) and the number of Monte Carlo iterations, 1000. The designed flow rate of the powerplant,
, was simulated for 1 m
3/s to 3.5 m
3/s. In practice, any duration of river flow rates, at any available time step, could be substituted providing peak and trough flow events, and are statistically represented in the simulation.
To compute the total energy produced for the purpose of computing the levelized cost of energy (LCOE), power in the water at the designed flow rate, as per Equation (1), is multiplied by net efficiency,
, and annualized to one year:
where
is the electricity generated in one simulated year.
2.3. Levelized Cost of Electricity (LCOE) Generated
A simplified cost model was utilized for plant sizing and costs, to provide general information on the tradeoffs between plant sizes, as is appropriate for the early planning stages simulated here.
The head is assumed for reasonable locations on the Sebeya river, and the type of turbine being used; a typical maximum gross head is 80 m. Costs can be broken into two components. The first component is an up-front investment in non-reoccurring engineering costs, permitting, legal fees, loan applications and processing fees, supervision cost, studies, etc. This value is uncertain, but for practical implementations in Rwanda, a reasonable range is between 250,000–500,000 USD for a small hydropower plant. Lacking any information about the distribution of these costs, we assumed uniform distribution.
The second component scales with the size of the plant. Using confidential data from two hydropower developers operating in Rwanda, the study assumed a required investment uniformly distributed between 2500–3000/kW USD. Using the renewable energy technologies cost analysis series report [
34] and estimates from hydropower owners in Rwanda, the study assumed the operation and maintenance (OM) cost per
$ capital cost uniformly distributed between 0.025–0.03
$/
$ capital investment.
To assume plant operating uptime and lifetime, owners of existing power plants were asked to estimate both parameters; resulting in an operating uptime of 0.85 and a plant design life of 18 years. Similarly, to assume a discounted rate, several banks in Rwanda were interviewed, resulting in a discounted rate of 12%.
Given the above assumptions, the LCOE for renewable energy is calculated as;
where
is the average levelized cost of electricity generated,
is investment expenditures in the year t,
is operation and maintenance expenditures in the year t,
is fuel expenditures in the year t,
is the electricity generated in year t,
is the discount and
is the economic life of the system.
2.4. Return on Investment (ROI)
Return on investors is a measure of business project profitablity.
where ROI is the return on investment,
is net gain and
is cost of investment. From return on investors, one can tell whether a business project is profitable or not [
35]. A positive return on an investment means a profitable business while a negative return on an investment means an unprofitable business.
4. Conclusions
This study illustrates a methodology for analyzing RoR hydropower plants in a tropical climate where rivers are characterized by highly variable flows, which are undocumented or poorly documented. This method allows for the analysis of RoR plants by using realistic flow distributions, rather than relying on reservoir assumptions that smooth river flow rates or single-point estimates based upon assumed ‘average flows’. Indeed, assuming average flow—a common reservoir assumption—will often over-size a small hydropower plant, producing an unexpectedly high LCOE in operational conditions. Additionally, an oversized plant may wear more quickly than one that is properly size, an issue often seen in the field in Rwanda.
The study found that a relatively simple neural network model, trained across two complete hydrological seasons, produced a similar flow behavior to the available river flow rate measurements. While not ideal, this type of extrapolation is valuable for early planning decisions, when flow data is not available, and common average flow assumptions are not valid for RoR plants.
Finally, using a Monte Carlo analysis approach provides a template to capture the uncertainty in inputs for first-order analyses of small power plants. The approach requires reasonable up-front engineering effort and easily captures variability in values for many key variables. The LCOE results coupled with the uncertainty estimates provides the analyst with both a specific target LCOE for decision-making, and guidance on the confidence to place in that number.
Finally, absent any similar analyses, investors are advised to design RoR plants at, or just below, the ‘knee’ in the flow exceedance curve. The Sebeya river is typical of tropical conditions, and plants designed substantially above that knee produced a significantly higher LCOE. This observation also illustrates that time series of daily river flow—and corrections to those flows for expected withdrawals over the life of the plant—are a critical input to investment decisions. Jurisdictions interested in small hydropower development are therefore encouraged to invest in river monitoring as a catalyst to development.