Research Article
Economic assessment of WECS for water pumping systems in the North Region of Cameroon
^{1}
Department of Renewable Energy, National Advanced School of Engineering of Maroua, University of Maroua, P.O. Box 46 Maroua, Cameroon
^{2}
Department of Physics, Higher Teachers' Training College, University of Maroua, P.O. Box 46 Maroua, Cameroon
^{3}
Laboratoire Matériaux Optiques, Photonique et Systèmes, LMOPS, EA 4423, Université de Lorraine, Metz, 57070, France
^{*} email: kidmokaoga@gmail.com
Received:
20
December
2020
Received in final form:
15
March
2021
Accepted:
16
March
2021
This paper analyses the potential utilization of wind electric pumping for water distribution in offgrid locations of the North Region of Cameroon (NRoC), using ground measured data as well as longterm satellitederived data. Furthermore, this paper puts emphasis on statistical indexes of accuracy for the comparison of measured (2007–2012) and longterm satellitederived (2005–2020) data. The outcome of this study clearly indicates that longterm satellitederived data, obtained through the Prediction of Worldwide Renewable Energy Resources, can be considered as a viable alternative to missing sitespecific data from ground stations, mainly in developing countries. The accuracy of satellitebased wind resource is deemed sufficient to provide a reasonable assessment in the initial phase of wind project planning, before in situ measurements with high accuracy are available. To model wind speeds characteristics, the energy pattern factor method (EPF) has been utilized as one of the reliable Weibull distribution methods for the assessment of wind energy potential at desired turbines heights in the selected locations. Five 20kW pitchcontrolled wind turbines (WT) with a hub height of 30 m, are considered to evaluate the power output and energy produced. Volumetric flow rates and costs of water produced are estimated. The results showed that, out of the eight chosen locations, Figuil shows the best combination of capacity factors (CF) and costs of energy (COE), while the site of Poli displays the worst values of CF and COE, no matter which of the five WT is employed. Thus, selecting WT for low wind speeds regimes, should require to combine location wind resource and WT characteristics such as lower cutin wind speed (1.5 m/s) as well as lower rated wind speed (8.0 m/s) in order to successfully and affordably implement the wind technology in the NRoC.
© D. Kaoga Kidmo et al., published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The sustainable development of a country is closely linked to the utilization of renewable energy resources (RES) at its disposal [1]. RES provide not only an important source of employment, wealth and profits to the country, but also play a critical role in building a sustainable and modern society [2–4]. Cameroon is endowed with conventional sources of energy such as natural gas and oil fuels, however associated with the emissions of both air pollutants, greenhouse gases and global warming. These conventional sources of energy share the fuel sources for electrical energy production (natural gas and oil fuels) in the country and other energy needs of industries and populace [5]. Moreover, fossil fuel sources, which are gradually depleting, need to be preserved and utilized efficiently [6]. As the global energy demand increases and fossil fuels become a gradually expensive commodity to humankind, the focus has been shifted towards renewable energy for a sustainable development [7].
RES, such as biomass, geothermal, hydro, solar, wind are inexhaustible, clean, free and offer many benefits such as nonpolluting and almost inexhaustible characteristics in contrast to fossil fuel [8]. To further add to RES advantages, they can be harnessed on a local basis for applications, such as decentralized production of electricity in remote or rural areas, with low density of population [9].
Amongst RES, wind energy is one of the world's most abundant and permanent energy source, with a highlevel of utilization. In addition, the potential of wind energy is credited, in the best of cases, with forty times the global electricity demand [10]. Over the last twenty years, the yearly growth rate of global wind energy has achieved an average growth rate of 22% [11]. Wind energy conversion systems (WECS) are considered as a mature and costcompetitive technology for electricity production [12], which is produced locally [13], but its dissemination remains marginal in most African countries. WECS necessitate to correctly characterize the wind resource over a location in addition to its geographic distribution [14]. As result, wind power is poised to play a key role in meeting global growing demand for electricity, especially in developing countries such as Cameroon, where it can be harnessed at an affordable cost on a local basis for water pumping applications.
The NRoC lays in a low wind speed regime and as a result the region is faced with numerous challenges in developing wind energy [15]. Measured wind data are mostly available at the ground measurements station at the Garoua main airport, while for the seven other considered sites located beyond a radius of at least 50 km of the main measuring wind station, there are no ground measurements stations. Therefore, longterm satellitederived data [16], obtained through the POWER (Prediction of Worldwide Renewable Energy Resources), are deemed suitable and a viable alternative to missing sitespecific data from ground stations [17–19]. These satellitederived data are at daily level, at a 0.5 × 0.5degree resolution. Additional scientific literature related to the use of longterm satellitederived data and their accuracy in comparison with in situ measurements are provided as referenced [20–23].
The aim of this paper is to acknowledge the use of satellitebased wind resource as appropriate data to provide a reasonable assessment, before higheraccuracy in situ measurements are accessible. Furthermore, this work explores wind power production using electric systems for water pumping in the NRoC. Measured wind speeds at 10 m height above ground level (agl), in hourly timeseries format, from January 2007 to January 2012 were collected through a cupgenerator anemometer for Garoua. When measured ground level data remained nonexistent, satellitederived data for the sites of Bashéo, Beka, Figuil, Pitoa, Poli, ReyBouba and Touboro, were obtained from the NASA Langley Research Center (LaRC) POWER Project funded through the NASA Earth Science/Applied Science Program [24]. The accuracy of satellitederived data for Garoua (the only site with ground measurements) was assessed using statistical indexes of accuracy, namely coefficient of determination (R^{2}), index of agreement (IOA), mean bias error (MBE), root mean square error (RMSE) and relative root mean square error (RRMSE).
Kidmo et al. [5] summarized most recent studies performed across Cameroon to investigate wind speeds and power characteristics, wind energy and costs of WECS. Recent and relevant research studies related to Garoua are detailed as followed. Kidmo et al. [25] performed a statistical analysis of wind speed distribution using six Weibull Methods for wind power assessment. Kidmo et al. [26] evaluated wind energy potential for small scale water pumping. Less attention was so far paid to the combination of location wind resource and wind turbines (WT) characteristics such as cutin wind speed as well as rated wind speed. This paper highlights the use of low wind speeds WT to take full advantages of local wind resource and WT characteristics. Capacity factors (CF), costs of energy generated (XAF/kWh) and water produced (XAF/m^{3}) are estimated using the Present Value Cost (PVC) method of energy produced per year, considering the frequency at which WT produce power.
To model wind speeds characteristics, the energy pattern factor method (EPF) [27,28] has been utilized as one of the reliable Weibull distribution methods for the assessment of wind energy potential at desired turbines heights in the selected locations. Five 20kW pitchcontrolled WT with a hub height of 30 m, are considered to evaluate the power output and energy produced. Volumetric flow rates and costs of water produced are estimated using a 50 m pumping head, for the sake of simplicity. The results showed that, out of the eight chosen locations, Figuil shows the best combination of capacity factors (CF) and costs of energy (COE), while the site of Poli displays the worst values of CF and COE, no matter which of the five WT is employed. Thus, selecting WECS for low wind speeds regimes, should require to combine location wind resource and WT characteristics such as lower cutin wind speed (1.5 m/s) as well as lower rated wind speed (8.0 m/s) in order to successfully and affordably implement the wind technology in the NRoC.
2 Methodology
2.1 Study area
The North Region as shown in Figure 1, is one of the ten regions of Cameroon, Located between 7.30 − 10.00 N latitude and 12.30 − 15.00 E longitude. It covers an area of 66 090 km^{2} and is bordered by two regions of the country, the Far North Region to the north and the Adamawa Region to the south, Chad to the east, Central African Republic to the southeast and Nigeria to the west. The North region accounts for one of the driest regions of Cameroon under the influence of the harmattan winds in the dry season which lasts eight months, temperatures rise at their highest (40–45 °C) and there is no rainfall. During the rainy season which lasts four months, average rainfall is between 900 and 1500 mm and torrential rains are observed as well as lower temperatures compared to dry season. The Benoue depression constitutes the primary land feature of the North region, with an elevation in the range of zero to 200 m. Land elevations are in the range of 200–500 and 500–1000 m, respectively in the north and south of the region [29]. The selected sites, namely, Bashéo, Beka, Figuil, Garoua, Pitoa, Poli, ReyBouba and Touboro, as shown in Figure 1, have been chosen based on relevant landscape information they provide in order to fully assess the potential of wind energy in the NRoC.
Fig. 1
Map showing the eight selected sites of the North Region of Cameroon. 
2.2 Wind data description and source
Two sources of data have been utilized in this work, measurements from a ground station for the site of Garoua and longterm satellitederived data for the seven other sites. Wind speed measurements for the locality of Garoua, at 10 m height agl, in hourly timeseries format, from January 2007 to January 2012 have been collected through a cupgenerator anemometer at the Garoua, main meteorological station [25]. For the seven other selected sites of the north region, measured ground level data are nonexistent. Therefore, longterm daily satellitederived data, from January 2005 to January 2020, obtained through the POWER, have been utilized. These data were obtained from the NASA Langley Research Center (LaRC) POWER Project funded through the NASA Earth Science/Applied Science Program [24]. Table 1 provides geographical coordinates of the eight sites, as well as mean ambient temperature and measurement period.
Geographical data for eight selected sites of the North Region.
2.3 Mean wind speed and standard deviation
Monthly mean wind speed v_{m} and standard deviation σ of wind speed data are calculated as Equations (1) and (2) [25]:(1) (2)Where:
v_{m} = mean wind speed [m/s]
σ = standard deviation of the mean wind speed [m/s]
v_{i}= wind speed [m/s]
N = number of wind speed data.
After the statistical analysis of wind data, monthly mean wind speeds and standard deviations for Garoua using measured and satellitederived data, are summarized in Table 2, while Figures 2 and 3 show, respectively monthly mean wind speeds and standard deviations using satellitederived data for for Bashéo, Beka, Figuil, Garoua, Pitoa, Poli, ReyBouba and Touboro. The magnitude of monthly mean measured wind speeds and standard deviations for Garoua at 10 m height agl, lay respectively within the range of 1.24 to 2.71 m/s and 0.73 to 1.39 m/s, while corresponding values using satellitederived data range from 1.67 to 3.38 m/s and from 0.47 to 1.00 m/s, in that order.
Monthly mean wind speeds and standard deviations for Garoua using measured and satellitederived data.
Fig. 2
Monthly mean wind speeds using satellitederived data for Bashéo, Beka, Figuil, Garoua, Pitoa, Poli, ReyBouba and Touboro. 
Fig. 3
Monthly standard deviations using satellitederived data for Bashéo, Beka, Figuil, Garoua, Pitoa, Poli, ReyBouba and Touboro. 
2.4 Weibull probability density function
The Weibull probability density function (PDF), which is a special case of a generalized twoparameter Gamma distribution, has been extensively utilized in scientific literature for wind speed forecasting and wind energy potential assessment [30–32]. Weibull PDF can be characterized by its PDF f (V) using Equation (3) and cumulative distribution function (CDF), F (V) using Equation (4) [33]:(3) (4)Where:
f (v) = probability of observing wind speed v;
v = wind speed [m/s];
C = Weibull scale parameter [m/s];
k = Weibull shape parameter.
Weibull parameters k and C are typically obtained using wellestablished estimation methods [25,34]. In this paper, Weibull parameters are computed using the energy pattern factor method (EPF). The first step is to define the power density expressed using Equation (5) [35].(5)Where:
ρ = air density at the site;
The air density (in kilograms per cubic meter) at a given location is calculated as the mass of a quantity of air (in kg) divided by its volume (in cubic meter). It depends on elevation and temperature above sea level and can be computed [36] using Equation (6).(6)Where:
Z = elevation (m);
T = temperature at the considered site (°K).
The second step of the EPF method is to determine the energy pattern factor (E_{pf}), as Equation (7) [37]:(7)Where:
v_{m} = mean wind speed [m/s];
= mean of wind speed cubes.
Once the energy pattern factor is calculated by using Equation (7), the Weibull shape parameter is estimated from Equation (8) [38]:(8)
where k is the dimensionless shape parameter illustrating how peaked the wind distribution is, while the scale parameter C indicates how ‘windy’ the wind location under consideration is [39].
The scale parameter is obtained using Equation (9):(9)
Where the standard gamma function is given by [40]:(10)
The gamma function can be approximated by [25,41]:(11)
2.5 Statistical indexes of accuracy
In order to evaluate the performance of the considered satellitederived data of Garoua, the following statistical indexes of accuracy are utilized:
1. Mean Bias Error (MBE) of satellitederived values using Equation (12) [42,43]:(12)where Gr_{i} denotes cumulative frequency distribution (CFD) of measured values at wind speed v_{i}, in time step i, Sat_{i} the respective satellitederived value, and N the number of nonzero wind speed data points. MBE provides basically the difference between CFD of average satellitederived and average measured values. Satellitederived values are overestimated if MBE > 0, while they are underestimated if MBE < 0.
2. Root mean square error (RMSE) using Equation (13) [44,45]:(13)
RMSE provides the deviation between satellitederived and measured values. Successful forecasts correspond to low values of RMSE, while higher indicate deviations [43]. RMSE should be as close to zero as possible.
3. Relative root mean square error (RRMSE) using Equation (14) [30,46]:(14)
RRMSE is calculated by dividing RMSE to the average of CFD of measured values. Ratings of satellitederived values' accuracy can be defined as [47]:
4. Coefficient of determination (R^{2}) using Equation (15) [30,48]:(15)
R^{2} determines the linear relationship between CFD of average satellitederived and average measured values. A higher R^{2} represents a better fit using satellitederived data and the highest value it can get is 1.
Where: Gr_{i} is the i^{th} CFD of ground (measured) wind speeds, Sat_{i} is the i^{th} CFD of satellitederived wind speeds, is the mean value of Sat_{i}, N the number of nonzero wind speed data points.
5. Index of Agreement (IOA) using Equation (16) [30,49]:(16)Where: denotes the average of Gr_{i}, is the mean value of Sat_{i}, N the number of nonzero wind speed data points. IOA estimates the accuracy of satellitederived wind speeds in predicting ground (measured) wind speeds. IOA values range from zero to one. IOA values above 0.5 imply efficiency in using satellitederived data [43]
Excellent:  RRMSE < 10%; 
Good:  10% < RRMSE < 20%; 
Fair:  20% < RRMSE < 30%; 
Poor:  RRMSE > 30%. 
2.6 Extrapolation of wind speeds at different hub heights
In this study, five WT, using a 30 m hub height have been considered. Wind speeds data obtained at 10 m height agl, are therefore adjusted to the relevant WT hub height. The vertical extrapolation of scale and shape parameters at an elevation (z), is defined by the relationships in Equations (17) and (18) [50,51]:(17) (18)Where z and z_{10} are in meters and the power law exponent n is given by Equation (19):(19)Where, the scale C_{10} and shape k_{10} parameters are determined at 10 m height agl.
2.7 Mean wind power density and energy density
Wind power density (P (v)), in watts per square meter, is estimated using Equation (20) [52].(20)
The mean wind power density (p_{D}) based on the Weibull probability density function can be calculated using Equation (21) [33].(21)The mean energy density (E_{D}) over a period of time T is expressed as Equation (22).(22)Where:
ρ = air density at the site;
A = swept area of the rotor blades [m^{2}].
2.8 Wind turbine and electric pumping systems
Five WT, from different manufacturers are considered, as illustrated by Table 3. The five WT are represented by WT_{1}, WT_{2}, WT_{3}, WT_{4}, and WT_{5} to avoid using registered names, trademarks etc. These WT are designed for low wind power density regime, for locations falling into Class 2 or lesser, based on PNL wind power classification scheme[50]. Although a 50 m pumping head is considered to calculate volumetric flow rates, a different pumping head can be considered, since the volumetric flow rate of water is inversely proportional to the pumping head.
Characteristics of the selected wind turbines.
2.9 Power curve model and capacity factor
The typical power curve of a 20kW pitchcontrolled WT is shown in Figure 4. Four different zones are observed as illustrated. For WS in the range of zero to cutin wind speed (V_{C}), there is no output yielded. From the cutin to rated WS (V_{C} to V_{R}), the power increases with the wind speed (WS). The power output curve is assumed to show a quadratic power shape, based on the selected power curve model in this work. From the rated WS (V_{R}) to the cutout WS (V_{F}), the WT yields a constant output at the rated power (P_{R}), regardless of WS variations. For WS higher than V_{F}, there is no output yielded.
Fig. 4
Power Curve of a 20kW pitchcontrolled WT, with technical specifications of WT_{3}. 
2.10 Power curve model and capacity factor
The power curve model of a WT can be modelled by four parameters: the cutin wind speed (v_{c}), the rated wind speed (v_{R}), the cutoff wind speed (v_{F}) and the rated electrical power (P_{eR}). In this study, the power output curve model of a WT can be approximated with the parabolic law, using the combination of Equation (23) [54]:(23)
The average power output (P_{e,ave}) from the turbine, which is related to the total energy production can be determined from Weibull distribution as Equation (24):(24)
The capacity factor CF is defined as the ratio of the average power output (P_{e,ave}) to the rated electrical power (P_{eR}) of the WT. The capacity factor CF can thus be computed as Equation (25) [55]:(25)
2.11 Water pumping capacity
The net hydraulic power output (P_{out}) required to deliver a volume of water V_{w} (m^{3}) can be expressed using Equation (26) [56].(26)Where: ρ_{w} = water density [kg/m^{3}];
g = acceleration due to gravity [m/s^{2}];
H = pump head [m]; η = system efficiency.
The volumetric flow rate of water is therefore evaluated using Equation (27).(27)
When taking into account the efficiency of the pump (η_{PUMP} = 62% for the submersible electric pump model) the water pumping capacity rate (F_{w}) can be expressed as Equation (28).(28)
2.12 Costs analysis
In this work, the adopted method to evaluate the costs of energy (COE) and costs of water (COW) produced is the present value of costs (PVC) of energy produced per year using Equation (29) [57].(29)
Where the following assumptions are made to estimate the cost of energy produced by the considered WT:

I is the investment cost, which includes the WT price in addition to 20% for civil works and other connections;

Average specific WT cost per kW is USD 1175, for WT rated power between 20 and 200 kW [53]

n is the useful lifetime of WT in years (20 yrs); i_{o} is the nominal interest rate (16%);

S is the scrap value (10% of the WT price); i is the inflation rate (3.6%);

C_{om} is the operation and maintenance costs (7.5% of the investment cost).
The discount rate (r ) is determined using Equation (30) [58]:(30)
The Wind availability (A) is determined using Weibull CDF of wind speeds at which WT produce energy. The total energy output (E_{WT}) over the WT lifetime (in kilowatthour) is computed as Equation (31).(31)
The cost of energy (COE) per unit kWh using the PVC method can be estimated using Equation (32).(32)
The annual volume of water V_{w} (m^{3}/year) produced by the WT is determined from Equation (33).(33)
The cost of water (COW) per unit m^{3} using the PVC method can be estimated using Equation (34).(34)
3 Results and discussion
3.1 Measured vs satellite data
Average monthly Weibull CDF values helped to perform the statistical comparison between measured and satellitederived data. Table 4 shows statistical indicators for the accuracy of satellitederived wind speed for the site of Garoua, which is the only site with measured ground level data. With the exception of May, June and July, satellitederived values are slightly overestimated (MBE > 0) for the rest of the year. The yearly average MBE value of 0.021 indicates a slight overestimation of satellitederived data. MBE values are in all cases close to zero, which demonstrate a match between locally measured and longterm satellitederived data. RMSE values range between 0.142 in February and 0.002 in May. These values, close to zero, show successful forecasts. RRMSE values, from April to November, are less than 10%, denoting an excellent accuracy of satellitederived data. RRMSE values for January, February, March and December are good, with values ranging from 11.4% to 15.5%. Overall, the accuracy of satellitederived data for yearly average is excellent, with a RRMSE value of 6.1%. R^{2} values vary from 0.847 in May to 0.696 in December. R^{2} value of the yearly average is 0.809. R^{2} values are generally higher and represents a better fit using satellitederived data. IOA values range between 0.993 and 1.000 and this implies efficiency in using satellitederived data. Therefore, satellitederived data can accurately represent measured data and when measured ground level data are missing, satellitederived data can be suitable to assess the potential of wind energy. Similar studies on the use of satellitederived data are found in the literature. Schmidt et al. [23] performed a comparison between in situ measurements of wind speed from Wave Glider (WG) and satellite/reanalysis products, including the wind speed error per wind speed category, RMSE, bias, and correlation coefficients. Accuracy values are within the same ranges as the present work. Barthelmie and Pryor [21] worked on satellite sampling of offshore wind speeds to represent wind speed distributions. Their study have stressed the difficulties to obtain through in situ methods wind speeds over the oceans. Remmers et al. [22] observed a very strong positive linear relationship between satellitederived wind speed data and the in situ measurements. Ayompe and Duffy [59] assessed the energy generation potential of photovoltaic systems in Cameroon using satellitederived solar radiation datasets due to the lack of a reliable network of surface observation stations for collecting weather data in the country. These studies concluded that the accuracy of satellitederived data varies according to empirical model functions utilized. Besides, the accuracy of satellitebased wind resource is sufficient to provide a reasonable assessment in the initial phase of wind project planning, before higheraccuracy in situ measurements are available.
Statistical indicators for the accuracy of satellitederived wind speed for Garoua.
3.2 Wind characteristics
It is observed that the probability of observing higher wind speeds is highest for Figuil. In addition, the probability of observing lower wind speeds is highest for Poli. For each site, the average air density is based on daily average temperature and location values of elevation. Air density values in the North Region vary between 1.10 and 1.13 kg/m^{3}, while the corresponding elevations are in the range of 291.79–651.63 m. At 10 m height agl, annual average wind speeds range from 2.24 m/s for Poli to 3.16 m/s for Figuil, corresponding wind power densities vary from 8.43 W/m^{ 2 } for Poli to 28.89 W/m^{ 2 } for Figuil. With annual average wind power densities less than 100 W/m^{ 2 }, the wind resource in the North Region can be ranked as wind power class 1, based on the scheme proposed by Battelle—Pacific Northwest Laboratory (PNL) [60]. Class 1 areas are considered unsuitable for large scale wind power development. It can be further observed that, for each of the eight sites, the peak of the PDFs, which indicates the most frequent wind speed, skews towards the higher values of wind speed. Table 5 proposes statistical Weibull analysis of wind speed data for 30 m elevation (tower of WT) for the eight sites. For each site, the average air density was adjusted to reflect the 30 mhub height tower. Annual average Weibull parameters were extrapolated to the 30 m height agl, as well as average wind speeds, power densities and energy densities.
Table 6 presents for WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5}, the annual average frequency for wind speeds greater or equal to WT's cutin wind speed for each of the eight selectd sites. It's observed that the probability of WT to produce power is the highest using WT_{1}, while this availability is the lowest using WT_{5}.
Annual wind characteristics at WT's hub height for the eight sites.
Annual average frequency for wind speeds greater or equal to WT cutin wind speed for different sites.
3.3 Wind turbine and WEPS
Figure 5 shows a typical wind electric pumping system. In order to highlight WECS characteristics needed for low wind speed locations, the five WT, with a rated capacity of 20 kW and a 30 m tower hub height, were considered for uniformity in the comparison. These WT are represented by WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5,} to avoid the use of trade names. Cutin and rated wind speeds are different for each of the five WT. WT_{1} has the lowest cutin and rated wind speeds, while WT_{5} has the highest cutin and rated wind speeds as illustrated by Table 3. Five 20kW pitchcontrolled WT, whose power curves are illustrated by Figure 6, combine Figuil wind characteristics with WT characteristics. The WECS are equipped with a submersible electric pump model to produce energy from wind and pump water using the produced energy for each of the eight selected sites. For the purpose of this work, a total dynamic head of 50 m was adopted as the average in the NRoC. Any other pumping head could be used since the volumetric flow rate of water is inversely proportional to pumping head.
Fig. 5
Typical wind electric pumping systems (WEPS). 
Fig. 6
Power curves of the five selected 20kW pitchcontrolled WT, which combine Figuil wind characteristics. 
3.4 Cost of energy
Table 7 provides the values of capacity factors (CF) and costs of energy (COE), using WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5,} for the eight chosen sites. The site of Figuil shows the best combination of CF and cost of energy, no matter the WT used, followed by the sites of Bashéo and Pitoa. The site of Poli displays the worst CF and COE. For the site of Figuil, the CF is equal to 26.52% and the COE is 49.05 XAF/kWh using WT_{1}. WT_{5} shows the worst performance, with CF equal to 7.11% and COE equal to 420.17 XAF/kWh. WT_{5} displays the poorest performance amongst WT.
CF and COE using WT_{1}, WT_{2}, WT_{3}, WT_{4}, and WT_{5}, for the eight sites.
3.5 Flow rate capacity
Figure 7 presents annual average flow capacity (m^{4}/h) for WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5} for the eight selected sites. Figuil shows the best flow rates capacity, 1076, 695, 421, 239 and 125.6 m^{4}/h, respectively using WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5}. Poli displays the worst flow rates capacity, 228, 105, 37, 9 and 1.4 m^{4}/h. For the site of Figuil, the annual average flow capacity are 613.76, 412.67, 233.55 and 191.28 m^{4}/h, using WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5}, in that order. WT_{5} displays the poorest performance amongst WT.
Fig. 7
Annual average flow capacity (m^{4}/h) for WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5} for the selected sites of the North Region. 
3.6 Cost of water
Table 8 displays the costs of water (COW) and flow rate (m^{3}/day), using WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5}, for the selected sites. It is observed that Figuil presents the best combination of COW and flow rate energy, while Poli shows the worst values, independently of WT used. For a total dynamic head of 50 m, the COW for Figuil are 4.01, 9.63, 26.20, 81.57 and 294.23 XAF/m^{3}, respectively using WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5}, while corresponding values of flow rate energy are 1515, 631, 232, 74 and 21 m^{3}/day. WT_{5}, in all cases, shows the poorest performance.
COW and volumetric flow rate using WT_{1}, WT_{2}, WT_{3}, WT_{4}, and WT_{5}, for the eight sites.
4 Conclusion
The economic assessment of wind energy conversions systems, for electric water pumping systems in the NRoC, has been explored, using measured and satellitederived wind data. The accuracy of satellitederived wind speeds data was assessed using MBE, RMSE, RRMSE, R^{ 2 } and IOA indicators. Satellitebased wind resources provided a reasonable accuracy and were deemed a viable alternative to missing sitespecific data from ground stations in Bashéo, Beka, Figuil, Pitoa, Poli, ReyBouba and Touboro. The energy pattern factor method (EPF) modeled wind speeds characteristics for the assessment of wind energy potential at desired turbines heights in the selected locations. Five 20kW pitchcontrolled WT of 30 m hub height, assessed the power output and energy produced. The assessment of the wind resource at 10 m height agl showed average annual power densities in the range of 8.43–28.89 W/m^{ 2 }. At 10 m height agl, this wind potential suggests that the NRoC falls under class 1 using the scheme proposed by BattellePacific Northwest Laboratory (PNL). Therefore, the wind resource is deemed unsuitable for large scale WT applications. Hence, the need to use low wind speed WT such as WT_{1} for electric water pumping applications. Based on its higher performance in comparison to the four other WT, WT_{1}, with the lowest cutin (1.5 m/s) and rated (8.0 m/s) wind speeds, showed the highest capacity factors, thus the lowest costs of energy and water produced. As a result, cutin and rated wind speeds have strong impact on the capacity factor, while cutout wind speed has negligible impact on it. Choosing WT for low wind speeds regime, require to combine location wind resource and WT characteristics such as cutin and rated wind speeds in order to successfully and affordably implement the wind technology in the NRoC.
Nomenclature
k_{10} : Shape parameter at 10 m height agl [−]
k_{z} : Shape parameter at z meters height agl [−]
σ : Standard deviation of the mean wind speed [m/s]
n : Useful lifetime of WT in years (20 yrs) [year]
F_{w} : Water pumping capacity rate [m^{4}/h]
v_{mG} : Measured wind speed [m/s]
σ_{G} : Standard deviation of the mean measured WS [m/s]
v_{mS} : Satellitederived wind speed [m/s]
σ_{S} : Standard deviation of the mean satellite derived WS [m/s]
f (v): Probability of observing wind speed v
g : Acceleration due to gravity [m/s^{2}]
C : Weibull scale parameter [m/s]
k : Weibull shape parameter [−]
E_{pf} : Energy pattern factor [−]
RRMSE : Relative root mean square error [%]
R^{2} : Coefficient of determination [−]
Gr_{i} : i^{th} CFD of ground (measured) wind speeds
Sat_{i} : i^{th} CFD of satellitederived wind speeds
C_{om} : Operation and maintenance costs [%]
N : Number of nonzero wind speed data points
C_{10} : Scale parameter at 10 m height agl [m/s]
C_{z} : Scale parameter at z meters height agl [m/s]
E_{WT} : Total energy output over the WT lifetime [kWh]
V_{w} : Annual volume of water [m^{3}/year]
v_{c} : Cutin wind speed [m/s]
v_{R} : Rated wind speed [m/s]
v_{F} : Cutoff wind speed [m/s]
P_{eR} : Rated electrical power [kW]
P_{e,ave} : Average power output [kW]
Q_{w} : Volumetric flow rate (m^{3}⁄day)
ρ_{w} : Water density [kg/m^{3}]
PVC : Present value of costs [XAF]
η_{PUMP} : Efficiency of the pump [−]
P_{out} : Net hydraulic power output [kW]
i_{o} : Nominal interest rate [%]
COE : Cost of energy [XAF/ kWh]
COW : Cost of water [XAF/m^{3}]
RMSE : Root mean square error [−]
Acknowledgments
The authors gratefully acknowledge the NASA Langley Research Center (LaRC) POWER Project funded through the NASA Earth Science/Applied Science Program for providing wind and temperature data used in this project.
References
 K.O. Odeku, An analysis of renewable energy policy and law promoting socioeconomic and sustainable development in rural South Africa, J. Hum. Ecol. 40 , 53–62 (2012) [Google Scholar]
 U. Lehr, A. Mönnig, R. Missaoui, S. Marrouki, G. Ben Salem, Employment from renewable energy and energy efficiency in Tunisia − new insights, new results, Energy Procedia 93 , 223–228 (2016) [Google Scholar]
 M.S. Salehudin, D.K. Prasad, P. Osmond, Renewable energy potential for energy efficient resort development in Malaysia, Solar 2011, the 49th AuSES Annual Conference , 2011, no. December, pp. 1– 10 [Google Scholar]
 M.O. Oseni, An analysis of the power sector performance in Nigeria, Renew. Sustain. Energy Rev. 15 , 4765–4774 (2011) [Google Scholar]
 D.K. Kidmo, K. Deli, B. Bogno, Status of renewable energy in Cameroon, Renew. Energy Environ. Sustain. 6 , 1–11 (2021) [Google Scholar]
 O.S. Ohunakin, Assessment of wind energy resources for electricity generation using WECS in NorthCentral region, Nigeria, Renew. Sustain. Energy Rev. 15 , 1968–1976 (2011) [Google Scholar]
 P.K. Halder, N. Paul, M.U.H. Joardder, M. Sarker, Energy scarcity and potential of renewable energy in Bangladesh, Renew. Sustain. Energy Rev. 51 , 1636–1649 (2015) [Google Scholar]
 K.S. Elie Bertrand, K. Abraham, M. Lucien, Sustainable energy through wind speed and power density analysis in Ambam, South Region of Cameroon, Front. Energy Res. 8 , 1–9 (2020) [Google Scholar]
 A. Fashina, M. Mundu, O. Akiyode, L. Abdullah, D. Sanni, L. Ounyesiga, The drivers and barriers of renewable energy applications and development in Uganda: a review, Clean Technol. 1 , 9–39 (2018) [Google Scholar]
 D.A. Salam, A. Riadh, Wind energy state of the art: present and future technology advancements, Renew. Energy Environ. Sustain. 5 , 1–8 (2020) [Google Scholar]
 D.K. Kidmo, B. Bogno, K. Deli, D. Goron, Seasonal wind characteristics and prospects of wind energy conversion systems for water production in the far North Region of Cameroon, Smart Grid Renew. Energy 11 , 127–164 (2020) [Google Scholar]
 IRENA, Global renewables outlook: energy transformation 2050 (Abu Dhabi, 2020). www.irena.org/publications [Google Scholar]
 M. Harries, Disseminating wind pumps in rural Kenya − meeting rural water needs using locally manufactured wind pumps, Energy Policy 30 , 1087–1094 (2002) [Google Scholar]
 C. Fant, B. Gunturu, A. Schlosser, Characterizing wind power resource reliability in southern Africa, Appl. Energy 161 , 565–573 (2016) [Google Scholar]
 J. Kenfack, J. Voufo, P.S. Ngohe Ekam, J.K. Lewetchou, U. Nzotcha, Overcoming local constraints when developing renewable energy systems for the electrification of remote areas in Africa, Renew. Energy Environ. Sustain. 5 , 1–10 (2020) [Google Scholar]
 NASA, POWER Single Point Data Access (2021). https://power.larc.nasa.gov/dataaccessviewer/ (accessed March 14, 2021) [Google Scholar]
 J.W. White, G. Hoogenboom, P.W. Wilkens, P.W. Stackhouse Jr., J.M. Hoel, Evaluation of satellitebased, modeledderived daily solar radiation data for the continental United States, Cimatol. Water Manag. 103 , 1242–1251 (2011) [Google Scholar]
 J. Bai, X. Chen, A. Dobermann, H. Yang, K.G. Cassman, F. Zhang, Evaluation of nasa satelliteand modelderived weather data for simulation of maize yield potential in China, Agron. J. 102 , 9–16 (2010) [Google Scholar]
 J.W. White, G. Hoogenboom, P.W. Stackhouse, J.M. Hoell, Evaluation of NASA satellite and assimilation modelderived longterm daily temperature data over the continental US, Agric. For. Meteorol. 148 , 1574–1584 (2008) [Google Scholar]
 C.B. Hasager et al., Remote sensing of environment offshore wind climatology based on synergetic use of Envisat ASAR, ASCAT and QuikSCAT, Remote Sens. Environ. 156 , 247–263 (2015) [Google Scholar]
 R.J. Barthelmie, S.C. Pryor, Can satellite sampling of offshore wind speeds realistically represent wind speed distributions? J. Appl. Meteorol. 42 , 83–94 (2003) [Google Scholar]
 T. Remmers, F. Cawkwell, C. Desmond, J. Murphy, E. Politi, The potential of advanced scatterometer (ASCAT) 12.5 km coastal observations for offshore wind farm site selection in Irish waters, 1–16 (2019) [Google Scholar]
 K.M. Schmidt, S. Swaart, C. Reason, S.A. Nicholson, Evaluation of satellite and reanalysis wind products with in situ wave glider wind observations in the Southern Ocean, J. Atmos. Ocean. Technol. 34 , 2551–2568 (2017) [Google Scholar]
 NASA, NASA Surface meteorology and Solar Energy: RETScreen Data (2020). https://eosweb.larc.nasa.gov/cgibin/sse/sse (accessed Dec. 07, 2020) [Google Scholar]
 D.K. Kidmo, R. Danwe, Y.S. Doka, N. Djongyang, Statistical analysis of wind speed distribution based on six Weibull Methods for wind power evaluation in Garoua, Cameroon, J. Renew. Energies 18 , 105–125 (2015) [Google Scholar]
 D.K. Kidmo, N. Djongyang, S.Y. Doka, D. Raidandi, Assessment of wind energy potential for small scale water pumping systems in the north region of Cameroon, Int. J. Basic Appl. Sci. 3 , 38–46 (2014) [Google Scholar]
 S.a. Akdaǧ, A. Dinler, A new method to estimate Weibull parameters for wind energy applications, Energy Convers. Manag. 50 , 1761–1766 (2009) [Google Scholar]
 N. Izadyar, H. Chyuan, W.T. Chong, K.Y. Leong, Resource assessment of the renewable energy potential for a remote area: a review, Renew. Sustain. Energy Rev. 62 , 908–923 (2016) [Google Scholar]
 North Region (Cameroon), From Wikipedia, Free Encycl., 2020. https://en.wikipedia.org/wiki/North%7B_%7DRegion%7B_%7D(Cameroon)%7B#%7DClimate (accessed Oct. 10, 2020) [Google Scholar]
 K. Mohammadi et al., Assessing different parameters estimation methods of Weibull distribution to compute wind power density, Energy Convers. Manag. 108 , 322–335 (2016) [Google Scholar]
 K. Mohammadi, O. Alavi, A. Mostafaeipour, N. Goudarzi, M. Jalilvand, Assessing different parameters estimation methods of Weibull distribution to compute wind power density, Energy Convers. Manag. 108 , 322–335 (2016) [Google Scholar]
 P.A. Costa Rocha, R.C. de Sousa, C.F. de Andrade, M.E.V. da Silva, Comparison of seven numerical methods for determining Weibull parameters for wind energy generation in the northeast region of Brazil, Appl. Energy 89 , 395–400 (2012) [Google Scholar]
 B. Safari, J. Gasore, A statistical investigation of wind characteristics and wind energy potential based on the Weibull and Rayleigh models in Rwanda, Renew. Energy 35 , 2874–2880 (2010) [Google Scholar]
 P.A. Costa Rocha, R.C. de Sousa, C.F. de Andrade, M.E.V. da Silva, Comparison of seven numerical methods for determining Weibull parameters for wind energy generation in the northeast region of Brazil, Appl. Energy 89 , 395–400 (2012) [Google Scholar]
 S.F. Khahro, K. Tabbassum, A.M. Soomro, L. Dong, X. Liao, Evaluation of wind power production prospective and Weibull parameter estimation methods for Babaurband, Sindh Pakistan, Energy Convers. Manag. 78 , 956–967 (2014) [Google Scholar]
 B. Bailey, S. McDonald, D. Bernadett, M. Markus, K. Elsholz, Wind resource assessment handbook: fundamentals for conducting a successful monitoring program, Music Educ. J. 41 , 65 (1997) [Google Scholar]
 A. Azad, M. Rasul, T. Yusaf, Statistical diagnosis of the best Weibull methods for wind power assessment for agricultural applications, Energies 7, 3056–3085 (2014) [Google Scholar]
 L. Bilir, M. İmir, Y. Devrim, A. Albostan, Seasonal and yearly wind speed distribution and wind power density analysis based on Weibull distribution function Science Direct Seasonal and yearly wind speed distribution and distribution function, Int. J. Hydrogen Energy 40 , 15301–15310 (2016) [Google Scholar]
 I. Fyrippis, P.J. Axaopoulos, G. Panayiotou, Wind energy potential assessment in Naxos Island, Greece, Appl. Energy 87 , 577–586 (2010) [Google Scholar]
 L. Olatomiwa, S. Mekhilef, O.S. Ohunakin, Hybrid renewable power supply for rural health clinics (RHC) in six geopolitical zones of Nigeria, Sustain. Energy Technol. Assess. 13 , 1–12 (2016) [Google Scholar]
 J.F. Manwell, J.G. McGowan, A.L. Rogers, Wind energy explained, Theory, Design and Application (BAffins Lane, Chichester, 2002) [Google Scholar]
 A. De Meij et al., Wind energy resource mapping of Palestine, Renew. Sustain. Energy Rev. 56 , 551–562 (2016) [Google Scholar]
 C. Stathopoulos, A. Kaperoni, G. Galanis, G. Kallos, Wind power prediction based on numerical and statistical models, J. Wind Eng. Ind. Aerodyn. 112 , 25–38 (2013) [Google Scholar]
 T.P. Chang, Wind speed and power density analyses based on mixture Weibull and maximum entropy distributions, Int. J. Appl. Sci. Eng. 8 , 39–46 (2010) [Google Scholar]
 D.K. Kidmo, R. Danwe, N. Djongyang, S.Y. Doka, Comparison of five numerical methods for estimating Weibull parameters for wind energy applications in the district of Kousseri, Cameroon, Asian J. Nat. Appl. Sci. 3 , 72–87 (2014) [Google Scholar]
 C. Neme, Statistical analysis of wind speed profile: a case study from Iasi Region, Romania, Int. J. Energy Eng. 3 , 261–268 (2013) [Google Scholar]
 K. Mohammadi, A. Mostafaeipour, O. Alavi, N. Goudarzi, M. Jalilvand, Assessing different parameters estimation methods of Weibull distribution to compute wind power density, Energy Convers. Manag. 108, 322–335 (2016) [Google Scholar]
 T.R.R. Ayodele, A.A.a. Jimoh, J.L.L. Munda, J.T.T. Agee, Wind distribution and capacity factor estimation for wind turbines in the coastal region of South Africa, Energy Convers. Manag. 64 , 614–625 (2012) [Google Scholar]
 G. Gualtieri, S. Secci, Methods to extrapolate wind resource to the turbine hub height based on power law: a 1h wind speed vs. Weibull distribution extrapolation comparison, Renew. Energy 43 , 183–200 (2012) [Google Scholar]
 A. Mostafaeipour, M. Jadidi, K. Mohammadi, A. Sedaghat, An analysis of wind energy potential and economic evaluation in Zahedan, Iran, Renew. Sustain. Energy Rev. 30 , 641–650 (2014) [Google Scholar]
 O.S.S. Ohunakin, M.S.S. Adaramola, O.M.M. Oyewola, Wind energy evaluation for electricity generation using WECS in seven selected locations in Nigeria, Appl. Energy 88 , 3197–3206 (2011) [Google Scholar]
 A.N. Celik, A. Albani, M.Z. Ibrahim, A statistical analysis of wind power density based on the Weibull and Rayleigh models at the southern region of Turkey, Renew. Energy 29 , 593–604 (2004) [Google Scholar]
 S. Diaf, G. Notton, D. Diaf, Technical and economic assessment of wind farm power generation at Adrar in Southern Algeria, Energy Procedia 42 , 53–62 (2013) [Google Scholar]
 M.S. Adaramola, O.M. Oyewola, O.S. Ohunakin, O.O. Akinnawonu, Performance evaluation of wind turbines for energy generation in Niger Delta, Nigeria, Sustain. Energy Technol. Assess. 6 , 75–85 (2014) [Google Scholar]
 E.K. Akpinar, S. Akpinar, An assessment on seasonal analysis of wind energy characteristics and wind turbine characteristics, Energy Convers. Manag. 46 , 1848–1867 (2005) [Google Scholar]
 S.S. Paul, S.O. Oyedepo, M.S. Adaramola, Economic assessment of water pumping systems using wind energy conversions in the southern part of Nigeria, ENERGY Explor. Exploit. 30 , 1–18 (2012) [Google Scholar]
 S. Rehman, T.O. Halawani, M. Mohandes, Wind power cost assessment at twenty locations in the kingdom of Saudi Arabia, Renew. Energy 28 , 573–583 (2003) [Google Scholar]
 S. Mathew, Wind energy: fundamentals, resource analysis and economics, Berlin Heidelberg: SpringerVerlag, 2007 [Google Scholar]
 L.M. Ayompe, A. Duffy, An assessment of the energy generation potential of photovoltaic systems in Cameroon using satellitederived solar radiation datasets, Sustain. Energy Technol. Assess. 7, 257–264 (2013) [Google Scholar]
 A. Ilinca, E. McCarthy, J.L. Chaumel, J.L. Rétiveau, Wind potential assessment of Quebec Province, Renew. Energy 28 , 1881–1897 (2003) [Google Scholar]
Cite this article as: Dieudonné Kaoga Kidmo, Bachirou Bogno, Kodji Deli, Michel Aillerie, Bello Pierre Ngoussandou, Economic assessment of WECS for water pumping systems in the North Region of Cameroon, Renew. Energy Environ. Sustain. 6, 6 (2021)
All Tables
Monthly mean wind speeds and standard deviations for Garoua using measured and satellitederived data.
Statistical indicators for the accuracy of satellitederived wind speed for Garoua.
Annual average frequency for wind speeds greater or equal to WT cutin wind speed for different sites.
CF and COE using WT_{1}, WT_{2}, WT_{3}, WT_{4}, and WT_{5}, for the eight sites.
COW and volumetric flow rate using WT_{1}, WT_{2}, WT_{3}, WT_{4}, and WT_{5}, for the eight sites.
All Figures
Fig. 1
Map showing the eight selected sites of the North Region of Cameroon. 

In the text 
Fig. 2
Monthly mean wind speeds using satellitederived data for Bashéo, Beka, Figuil, Garoua, Pitoa, Poli, ReyBouba and Touboro. 

In the text 
Fig. 3
Monthly standard deviations using satellitederived data for Bashéo, Beka, Figuil, Garoua, Pitoa, Poli, ReyBouba and Touboro. 

In the text 
Fig. 4
Power Curve of a 20kW pitchcontrolled WT, with technical specifications of WT_{3}. 

In the text 
Fig. 5
Typical wind electric pumping systems (WEPS). 

In the text 
Fig. 6
Power curves of the five selected 20kW pitchcontrolled WT, which combine Figuil wind characteristics. 

In the text 
Fig. 7
Annual average flow capacity (m^{4}/h) for WT_{1}, WT_{2}, WT_{3}, WT_{4} and WT_{5} for the selected sites of the North Region. 

In the text 