Determination of a Suitable Solar Radiation Model for the Sites of Chad ()
1. Introduction
The energy sector of Chad, which is still weakly developed, is characterized by high consumption of wood fuels (wood and charcoal), which accounts for more than 90% of the total energy consumption of the country. The use of conventional energies (petroleum products and electricity) occupies a marginal part of the national energy balance. These energies, although crucial in the development of a modern economy, account for only about 10% of total energy consumption in the country. There is no interconnected network in the country. Chad, however, has significant energy potential such as hydrocarbons, biomass and renewable energies, including solar and wind energy, whose exploitation could have contributed to the development of the sector. Solar energy applications require, above all, knowledge of the global solar radiation of a site. Thus, a reliable estimation of global solar radiation for a site is fundamental. Doing this will allow adequate knowledge of how to channel such application for either electricity generation, water heating, or irrigation, to mention a few [1] [2] [3] . In a developing country, the data are recorded directly only in a few meteorological stations. In many stations, the spatial coverage of radio stations is insufficient. For example, in Chad, there is not only a lack of measuring equipment, but also of training technicians. In addition, lack of maintenance or calibration of solar radiation sensors gives erroneous measurements or missing data [4] . Researchers have developed a large number of methods to estimate global solar radiation due to lack of reliable data on solar radiation. These methods are based on an empirical model linking variables such as humidity [5] [6] , temperature [7] [8] [9] [10] , elevation [11] [12] , duration of the sun [13] [14] [15] [16] and latitude [17] [18] . It may be easier to exploit solar energy resources when the site under consideration is equipped with a pyranometer functioning regularly for several years. However, there is need to use approximate methods to predict solar radiation characteristics if local measures do not exist. Numerous models have been developed that connect global solar radiation to extraterrestrial solar radiation in order to estimate the amount of solar energy incident on a horizontal surface. Among these models, the Angstrom-Prescott model has been developed and employed by many researchers and with certain meteorological parameters [19] . In the estimation of global solar radiation, each of these factors contributes significantly.
This study was therefore focused to develop a mathematical model to estimate the solar radiation coming from extraterrestrial radiation, with meteorological and geographical data as governing parameters. The model was validated by comparing its results with experimentally measured data across the twelve sites of Chad. The impact of this is the fact that such model can now be used to analyse and make informed decisions on solar technology applications without recourse to several years of experimental measurements around the studied sites and across the region.
2. Data Base
Solar data employed for the study were obtained from the General directorate of the National Meteorology of Chad, in June 2014 covering the period of 63 years (i.e. 1950-2013). They are monthly data of the relative humidity, maximum temperature, minimal temperature, and sunshine duration. The geographical coordinates of the twelve stations of the National office of Meteorology (ONM) are given in Table 1.
3. Used Models
To calculate the global solar radiation, one has recourse to the ideal models. These models are in the form of empirical relations which connect the components of the solar radiation to the principal weather parameters and the astronomical parameters. The weather parameters are the ambient temperature, the relative humidity, the sunshine duration. Amongst the astronomical parameters one has the maximum duration of the day, the variation of the sun, the variation of the ground-sun distance and the solar radiation in the extraterrestrial radiation [20] . As part of this work, these include the Angstrom-Prescott model, the Allen model and the Sabbagh model.
3.1. Estimation of Extraterrestrial Radiation, H0
The monthly mean of the daily extraterrestrial solar radiation on a horizontal surface is determined according the following relation [21] [22] [23] [24] [25] :
(1)
Table 1. Geographical locations of the sites considered.
where δ and
are respectively the monthly mean of the daily solar declination and the sunshine hour angle defined by [26] [27] [28] :
(2)
(3)
where
is the solar constant (Isc =1367 W/m2), L is the location latitude, Dn is the number of the day in the year.
3.2. Global Solar Radiation on a Horizontal Level with Angstrom-Prescott Model
Angstrom has been the first to propose an ideal model (linear model) to estimate the horizontal global solar radiation with in entry the data over the sunshine duration [29] [30] [31] . Prescott and Page [32] [33] have improved this model and considered the horizontal global solar radiation of the sunshine duration according to the relation (4).
(4)
The parameters a, bare respectively defined by the Equations (5) and (6).
(5)
(6)
The possible maximum monthly mean of the daily sunshine duration is [34] :
(7)
3.3. Global Solar Radiation on a Horizontal Level with the Allen Model
Allen [35] [36] estimated the monthly mean of the global solar radiation as a function of H0, the monthly mean of the maximum temperature (TM), and the monthly mean of the minimum temperature (Tm) as:
(8)
Where Kr is defined as:
(9)
In the relation (9), Kra = 0.17 and P/P0 may be defined as:
(10)
where P and P0 are respectively the values of local and standard atmospheric pressure, and h is the altitude of the place in meters.
3.4. Global Solar Radiation on a Horizontal Level with the Sabbagh Model
While being based on data relating to several countries of the Gulf, in particular, the sites of Saudi Arabia, Sabbagh et al. developed two empirical relations binding the various weather parameters which affect the attenuation of the solar radiation, namely: sunshine duration, relative humidity, the maximum temperature, the altitude, the geographical situation (longitude, latitude) and its situation compared to the sea and a lake of water characterized by the characteristic factor of the zone, which is given by the following relation [37] :
, (11)
With:
, (12)
, (13)
RH and Tmax are respectively the monthly average per day of the sunshine duration, the relative humidity and the maximum average temperature of the considered month.
, climatic factor
n, number of the month considered
S, monthly average daily bright sunshine duration (h)
4. Statistical Test
In order to compare the data of the solar radiation provided by NASA (National Aeronautic and Space Administration) with those obtained from the various presented models, all the different models were implemented by creating a code using MATLAB and Excel. From each of these programs and for each studied site, on the one hand we drew up in the same graph, the values of whole-body radiation by NASA and those calculated, and in the other hand we drew up the relative error. The presented models in paragraph 3 permitted to evaluate the calculated global radiation Hi,c in order to be compared to the measured global radiation Hi,m.
Several statistical indicators used in the literature [20] [21] [37] - [42] were employed to determine the degree of accuracy of the estimated values as compared to measured values. These indicators included: RMSE, MBE, MPE and NSE.
RMSE (Root Mean Square Errors) is a measure of the variation of calculated values; it gives information about the performance of the model and is always positive values. The model is best when its RMSE value is the smallest. It is defined by the relation:
(14)
n is the number of the month.
The mean bias error (MBE) is the mean inclination error giving information on the performance of the long-term model. To this end, a negative value refers to underestimation, while positive value refers to an overestimation. It is given by Equation (15):
(15)
The MPE (Mean Percentage Error) is defined by the relation:
(16)
For this indicator, for a given model, an error expressed as a percentage between −10% and +10% is acceptable.
The NSE (Nash-Sutcliffe Equation) represents a measure of the precision of the model results. The NSE is defined by the relation:
(17)
where:
is the mean measured global radiation.
5. Results and Discussions
To estimate the overall solar radiation using the Angstrom-Prescott model, apart from the geographical coordinates of the site, the average daily insolation duration was considered.
Meteorological data such as relative humidity, maximum temperature and average daily sun exposure measured the overall solar radiation using the Sabbagh model. Maximum temperature and minimum temperature were considered in the Allen model.
Figure 1 presents the map of Chad with the three climatic zones. The different sites that make up these areas can be observed.
Figure 2 presents the results of the monthly extraterrestrial radiation for the twelve sites in the three climatic zones of Chad. It can be noted in this table that the smallest value of the extraterrestrial radiation of 8.112 kWh/m2 is found in Faya-Largeau. The highest value of 9.364 kWh/m2 is found at Moundou.
Figure 3 presents the values of the monthly global solar radiation in the form of histograms obtained through the three appropriate models of estimation for each studied site.
Tables 2-5 present the total results obtained through the relations of Angstrom-Prescott, Allen, Sabbagh equation and Sabbagh model for the twelve sites of the three climatic zones of Chad.
Figure 1. Presentation of the map of Chad with the three climatic zones.
Figure 2. Monthly extraterrestrial solar radiation.
In Table 2 the values of the various parameters of Angstrom-Prescott model for the twelve sites of the three climatic zones of Chad are presented. Thus based on Angstrom-Prescott model calculation, one can conclude that Am-timan is the suitable place for the exploitation of solar energy compared to the other sites, because of his highest monthly global solar radiation of 5.917 kWh/m². The smallest value of the monthly global solar radiation calculated through Angstrom-Prescott model is 5.457 kWh/m² (obtained for the site of Faya-Largeau).
Table 2. Angstrom-Prescott model for the twelve sites.
Table 3. Allen model for the twelve sites.
The Angstrom-Prescott model is the most suitable for the calculation of global solar radiation based on the relative error with NASA data especially for the sites: Bongor (1.590%), Pala (1.482%), Am-timan (−0.887%) and Mongo (−1.829%).
Table 4. Sabbagh equations model for the twelve sites.
Table 5. Sabbagh model for the twelve sites.
Table 3 presents the estimation of the global solar radiation through the relation of Allen shows that it is Bongor which has the best solar potential with a value of 6.085 kWh/m², compared to the other sites of Chad. On the other hand, one notices that in term of relative error with measured data, the Allen model is more adapted for the calculation of the global solar radiation of N’Djamena, Pala and Sarh, respectively with a mean relative error of 0.078%, 0.844% and 0.878%. Moreover, Abeche is the site where the temperature is the highest with a value of 37.42˚C while the lowest temperature is 20.54˚C (in the site of Am-timan).
Table 4 presents the obtained values of the Sabbagh equations for the twelve sites of Chad. Thus, we note in this table the variation of the coefficients which is of the order of
to Pala and
to Abeche.
Table 5 presents the obtained values of the parameters of Sabbagh model for the twelve sites of Chad. Abeche is the site where the sunshine duration is the highest (11.62 h); the smallest value is observed in Sarh (7.17 h).The highest value of the duration of the day is noted in Moundou (11.62 h) while the lowest value is observed in Faya-Largeau (11.20 h). For the relative humidity, it is in Sarh which has the highest value (59.31%); the smallest value is observed in Faya-Largeau (20.37%). In addition, the highest value of the maximum temperature is observed in Abeche (37.42˚C) and the lowest value is observed in Pala (33.92˚C). The highest valueof sunshine rate is observed in Abeche (1.02) and the lowest value is observed in Sarh (0.62). The estimation of the global solar radiation for all sites through Sabbagh model shows that the highest value is observed in Abeche (6.354 kWh/m2) and the lowest value is observed in Pala (5.285 kWh/m2). The various results show indeed that the three zones of Chad do not have the same climatic conditions.
The Sabbagh model is the most suitable for calculating global radiation based on the relative error with NASA data especially for the sites: Faya-Largeau (−6.766%), Abeche (3.704%), N’Djamena (−0.348%), Ati (0.356%) and Moussoro (3.400%).
It can also be noted that all three models can only be applied if weather data and geographic parameters are available for a given site. The reliability of these models is the correct estimate of global solar radiation without going through direct and diffuse radiation. Moreover, it is observed that the solar radiation is affected by the meteorological parameters because the decrease of the parameters such as the temperature and the relative humidity leads to the reduction of the solar radiation. For example, the month of August seems the most unfavorable because the more it rains, the more radiation decreases.
Table 6 compares the three models used for the twelve studied sites. The comparison between the obtained results through the three models of calculation and the measured data shows that the Angstrom-Prescott model gives the best estimation of the global solar radiation for the sites of Bongor (MPE (%) = 1.5902; RMSE = 0.3368; MBE = 0.0159; NSE = 0.6999), Pala (MPE (%) = 1.4822; RMSE = 0.3338; MBE = 0.0148; NSE = 0.7241), Am-timan (MPE (%) = −0.8875;
(a) (b) (c)
Table 6. Statistics for the validation of the selected models (kWh/m2).
RMSE = 0.4256; MBE = −0.0089; NSE = 0.291) and Mongo (MPE (%) = −1.8291; RMSE = 0.5609; MBE = −0.0183; NSE = −0.0373).
The Allen model is appropriate for Moundou (MPE (%) = −0.8437; RMSE = 0.3962; MBE = −0.0084; NSE = 0.7181), Sarh (MPE (%) = −0.8779; RMSE = 0.5114; MBE = −0.0088; NSE = 0.5245) and Bokoro (MPE (%) = 1.1811; RMSE = 0.799; MBE = 0.0118; NSE = −0.9464). For the cities of Faya−Largeau (MPE(%)= −6.7662; RMSE = 0.513; MBE = −0.0677; NSE = 0.6468), Abeche (MPE(%) = 3.7042; RMSE = 0.5418; MBE = 0.037; NSE = 0.0899), N’Djamena (MPE (%) = −0.3477; RMSE = 0.3842; MBE = −0.0035; NSE = 0.6927), Ati (MPE(%) = 0.3564; RMSE = 0.2314; MBE = 0.034; NSE = 0.5697) and Moussoro (MPE(%) = 3.4003; RMSE = 0.3879; MBE = −0.0677; NSE = 0.6468), it is the Sabbagh model which is adapted for the calculation of the global solar radiation.
We can justify that one model is more suitable than another if the statistical values tend towards zero.
The Comparison between the measureddataandthe estimated values of the monthly global solar radiation is presented in Table 7.
6. Conclusions
In this work, the most adapted mathematical model of estimating the global solar radiation has been determined for twelve sites of Chad. The main results show that Abeche, a site in the Sahelian zone, has an radiation of 6.354 kWh/m2, while Ati has solar radiation of 5.523 kWh/m2. In the Saharan zone to the north, Faya-Largeau has the best solar radiation potential around 6.013 kWh/m2. In southern of Chad, the site with highest solar radiation is Bongor (with 5.839 kWh/m2), but Sarh is the site with the lowest solar radiation of 5.751 kWh/m2. Three models, including the Angstrom-Prescott, Sabbagh and Allen models have been used to estimate the global solar radiation of each city, and these methods have been compared to the various statistical tests used to choose the appropriate model.
In terms of theoretical contributions, three equations were chosen:
Acknowledgements
The authors would like to thank the persons in charge of the National Meteorology for Chad, that of Bongor as well as the personnel who deal with the collection and the processing weather data on these sites, to have placed at our disposal the data used in this work.