Prediction of the Diffuse Solar Energy on Horizontal at Different Selected Locations ()
1. Introduction
Solar energy is considered one of the most important sources of renewable energy. Accurate knowledge of available solar energy with its direct and diffuse components in a particular place is of great importance in designing and sizing of solar energy conversion systems. A crucial input required in the simulation of buildings’ energy performance is the availability of detailed information on the magnitudes of diffuse and direct irradiance data. Moreover, configuration and sizing of solar energy systems (e.g. photovoltaic cells, solar-thermal collectors) necessitates reliable solar radiation measurements. However, concurrent measured data of global and diffuse irradiance on horizontal surface or direct normal solar irradiance are available only for a limited number of locations [1] [2] [3]. The measurement of global horizontal irradiance is rather simple and cost-effective. It can be, conceivably, an integral part of the sensory equipment of every building. Given global solar irradiation measurements on a horizontal surface (as the most widely available data), direct and diffuse solar radiation components can be obtained through various correlations [4]. The models are usually expressed in terms of first to fourth degree polynomial functions relating the diffuse fraction kd (ratio of the diffuse-to-global solar radiation) with the clearness index kt (ratio of the global-to-extraterrestrial solar radiation on horizontal surface), as well as to other variables such as solar altitude, air temperature, relative humidity. Although these models are typically derived following sound approaches, their performance appears to lessen once they are applied to regions other than those, which provided the initial data for model development [5].
Solar radiation data is the basic input for solar applications, such as photovoltaic, solar thermal systems and passive solar design. The data should be reliable and can be used at any time to design, optimize and evaluate the performance of solar technology at any site. However, it is not economically feasible to install solar radiation measuring instruments wherever possible. Therefore, the use of mathematical models to forecast the solar radiation in a given area has proved to be a viable option based on the measurement results of limited locations [6] [7] [8]. Unfortunately, in many developing countries, solar radiation measurement is not yet available because they cannot afford the equipment and measurement technology. Hence, it is vital to create methods for evaluating solar radiation based on more promptly accessible meteorological data. Within the design and execution examination of solar energy projects, particularly within the design and measure assurance of solar PV as a future alternative energy source, exact forecast of diffuse solar radiation (DSR) is essential. Careful thought of DSR can much better assess the productivity of the solar system [9] [10] [11] [12]. In addition, in several regions of the world, there is no or very little measurement of diffuse solar radiation. Because of their wide application in other places, they can measure total horizontal irradiance and other standard meteorological variables, such as sunshine duration, temperature and relative humidity. In the field of meteorology and agriculture, given global solar radiation data and some meteorological parameters, the diffusion component can be obtained through various correlations. Recently, extensive research has been conducted in many parts of the world to estimate DSR using the most widely available data. Many authors have proposed empirical formulas to use the clearness index (Kt: the ratio of global solar radiation to extraterrestrial solar radiation), or use the fraction of hours of sunshine to forecast the monthly average on a horizontal surface of DSR, or in combination. A model for predicting daily DSR using sunshine fraction (the ratio of sunshine duration to the maximum possible sunshine time), clarity index and haze factor is also proposed. Facts have proved that the model that uses both the clarity index and the sunshine score is the best choice for estimating DSR [13] [14] [15] [16] [17].
The utilization of solar energy will facilitate scale back the demand for standard energy. Therefore, owing to this, solar power is taken into account to be the proper resolution to the energy crisis facing the globe nowadays. During this method, radiation information should be established for places of interest that is sometimes a requirement for the institution and commission of star facilities. In any case, it is not possible to form careful observations of native environmental condition. This can be the case for several developing countries (such as Saudi Arabia). Though these countries could have high radiation potential, they lack comfortable radiation info results in fewer energy plans to be explored and enforced [18] [19] [20] [21] [22].
The quality of radiation is typically outlined in step with its composition, specifically incident and diffuse solar radiation. Among these elements, the number of diffuse radiation is usually unsure, because of additionally to location parameters; it is chiefly littered with several native geographic factors and climatically characteristics. Most accessible information bases are equipped with data on world radiation and need info on diffuse radiation. This can be because of the sometimes-higher value of putting in a meteorological workplace to examine elements. Therefore, empirical models are sometimes wont to measure diffuse radiation. The written kind proposes a range of various models to estimate the common monthly average radiation, that uses input factors (such as world radiation and daytime) and different climatically factors (such as humidness, pressure, precipitation, and temperature) to estimate. Among these, world radiation and sunshine amount are imperative factors used at intervals the advancement of experimental models for diffuse radiation [14] [23] [24] [25] [26].
The main objective of the present research is to correlate the monthly average daily diffuse fraction Hd = Hwith the monthly average daily clearness index Kt or the fraction of possible number of sunshine hours S/So for selected locations in Saudi Arabia. In addition, the ratio of Hd to H0 or the diffuse transmittance is correlated the ratio S/So for the considered locations. Therefore, the data available for the considered locations are then combined and used to develop correlations for all Saudi Arabia that can be used for estimation of Hd for any location of Saudi Arabia. The data were obtain; part from the Meteorological and Environmental Protection Agency (MEPA) and other part from the station of Al-Baha University in Saudi Arabia.
2. Data and Methodology
In the present research, the monthly average of daily global solar radiation H, diffuse solar radiationand the number of bright sunshine hours S available for three selected sites in Saudi Arabia (Al-Aqiq, Hail and Dammam) for the years from 1980 and 2019 are used. The geographical information of the selected locations are summarize in Table 1.
In the present research, the regression analysis is used for the proposed models, where the estimate and is the diffuse fraction (Kd) or diffuse coefficient (KD) and the predictors are sunshine ratio (S/So) and clearness index (Kt). Thus, the models can be defined for diffuse fraction and diffusion coefficient as follows [27] [28] [29] [30] [31]:
(1)
(2)
where H0, H, and Hd are the monthly mean daily extraterrestrial solar radiation, global solar radiation and diffuse solar radiation on a horizontal surface, respectively. Mathematically, sunshine ratio (S/So) and clearness index (Kt = H/Ho), S and So are the sunshine duration and maximum possible sunshine duration, respectively.
The monthly average daily extraterrestrial solar radiation on a horizontal surface is calculated from the following equation [32] - [38]:
(3)
where Hsc is the solar constant, nis the Julian day of the year, φis the location latitude, and δis declination angle, ωs is the sunset hour angle. δand ωsare mathematically defined as:
(4)
(5)
We can obtained the maximum possible sunshine duration (So) from
(
).
Table 1. The geographical information of the selected sites in the present research.
The correlations to which the measured data are fitted are as follow [39] - [44]:
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
In the present research, the results mentioned in the above section are used with the following correlations to express the dependence of diffuse radiation on various parameters in models of Equations (6)-(14) as follow [41] [42]:
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
(23)
where a,b, c and d are regression coefficients that depend on the site. The measured values of global solar radiation and diffuse solar radiation obtained with Epply Pyranometer, together with the corresponding sunshine duration values for the different selected locations in Saudi Arabia in the present research.
The accuracy of estimation of Hd is tested by calculating the mean bias error (MBE), root mean square error (RMSE) and the mean percentage error (MPE). Generally, low values of RMSE and MPE are desirable. Positive MBE shows overestimation while negative MBE indicates under estimation. The MBE, RMSE and MPE are defined as in the following equations [45] [46] [47] [48] [49]:
(24)
(25)
(26)
where
and
are the ith calculated and measured values of Hdand n is the number of observations.
3. Results and Discussion
The performance and evaluation of the statistical indicators mean bias error MBE, root mean square error RMSE, mean percentage error MPE%, correlation coefficients R, stander error S.E. and regression coefficients for the different selected location in Saudi Arabia during the period time from 1980 to 2019 are listed in Tables 2-4. From the analysis of the combined data for selected locations in the present research, relationships are obtained to express the diffuse radiation from various parameters. The obtained values of the regression coefficients, statistical indicators, correlation coefficient (R) and standard error estimate (S.E) of the models (15)-(23) are summarize in Tables 2-4. From these tables, we may noticed that, the values of regression constants a, b, c and d for all models are different values according to geographical information of site to another one during the period time in the present research, with the stander error estimate (S.E.). From the results of Tables 2-4, it is seen that the values of correlation coefficients are higher than 0.86 for all models with except in a few models. The Models (17), (20) and (23) are given the higher values of correlation coefficients 0.935, 0.962 and 0.974 in Al-Aqiq site respectively, while Models (16), (20) and (21) are given 0.957, 0.936 and 0.942 in Hail site respectively, also Models (17), (20) and (21) are given 0.957, 948 and 965 in Dammam site respectively. For statistical indicators, the values of MBE for all models are given some negative and other positive, this due to the correlations are overestimate and
Table 2. The values of statistical indicators MBE, RMSE, MPE%, R2, S.E. and regression coefficients for Al-Aqiq site in the present study.
Table 3. The values of statistical indicators MBE, RMSE, MPE%, R2, S.E. and regression coefficients for Hail site in the present study.
underestimate values of diffuse fraction Hd. The low values of the RMSE for all correlations indicate good agreement between the measured and calculated values of Hd. The negative values of MPE for all models show that for all locations, the proposed correlations slightly overestimate Hd. The results in this research are similar results in the previous study [24] [35]. Generally, the absolute values of MPE never reach 1.35%. This means that, very good agreement between the measured and calculated data and a good fitting exists between the monthly average daily values of the diffuse fraction (Hd/H)and (Kt) or (S/So) and between the diffuse transmittance (Hd/H0) and (S/So) for all selected locations in the present research.
Table 4. The values of statistical indicators MBE, RMSE, MPE%, R2, S.E. and regression coefficients for Dammam site in the present study.
The relationship between measured and calculated values of Hd by using different models for selected locations during the period time in the present research are show in Figures 1-3. From these figures, we indicate that the behavior of the measured values of the diffuse solar energy and the calculated values, which obtained from models 15 to 23 for all selected sites in the present study are nature shape with the path of the diffuse solar energy in the day. On the other meaning from these figures, the maximum values occur around in summer months, while the minimum occur in winter months for all location in the present research. The differences between measured values of diffuse solar energy and estimated by different models varies from 1.5% to 3.4% for all sites in this study. Generally, the results of these figures and the latter results which confirmed graphically for each month for all locations under present research. Therefore, the obtained results are show in Figures 1-3 for sites in this study. From the results of Figures 1-3, it is concluding that the second and third order correlations between Hd/Hand Kt do not significantly improve the accuracy of estimation of diffuse radiation see in Figure 1. On the other hand, a second order correlation between Hd/Hor Hd/H0 and S/So gives the desired accuracy of estimation of diffuse radiation for all considered locations see in Figure 2 and Figure 3 and the values of regression coefficients in Tables 2-4. The percentage error for a single month rarely exceeds ±7% for any of the locations in the present research. Furthermore, Models (15), (19) and (22) with the aid of the values of the regression constants given in Tables 2-4 are recommended for estimating the horizontal monthly average diffuse radiation at selected locations in the present study. The results in this research are similar results in the previous study [24] [35].
The empirical models correlations in the form of Equations (15), (19) and (22) for all Saudi Arabia are developed. They may then be used for estimating Hd for
Figure 1. The relationship between measured and calculated values of Hd by using different models at Al-Aqiq site in the present research.
Figure 2. The relationship between measured and calculated values of Hd by using different models at Hail site in the present research.
Figure 3. The relationship between measured and calculated values of Hd by using different models at Dammam site in the present research.
any location of Saudi Arabia. For this purpose, the measured data available from the selected sited in the present research are combined and analyzed (3 × 12 = 36 sets of values). Figure 4 and Figure 5 provide the variations of Hd/Hwith Kt and S/So, respectively. Figure 6 shows the variations of Hd/H0 and S/So using the data collected for the three selected locations in the present research. A linear correlation between Hd/Hand Ktfound to fit the measured data see in Figure 4. Furthermore, second order correlations between Hd/H and S/So see in Figure 5, and Hd/H0 and S/So are remarkable in Figure 6.
The following correlations have been obtained for all Saudi Arabia:
Model A:
,
(27)
Figure 4. The correlation of the diffuse fraction with the clearness index in the present research.
Figure 5. The correlation of the diffuse fraction with the sunshine hours in the present study.
Model B:
,
(28)
Model C:
,
(29)
The models from (27) to (29) are then employed for calculating the diffuse solar energy Hd for the selected locations in the present research. Figures from (7) to (9) present comparisons between the measured and calculated diffuse solar energy Hd for Al-Aqiq, Hail and Dammam locations respectively. From the results of Figures 7-10, it is obvious that there is a good agreement between the
Figure 6. The correlation of the diffuse transmittance with the sunshine hours in the present research.
Figure 7. The difference between measured and calculated values of Hd at Al-Aqiq site in the present study.
Figure 8. The difference between measured and calculated values of Hd at Hail site in the present study.
Figure 9. The difference between measured and calculated values of Hd at Dammam site in the present study.
measured and calculated data of Hd. The best estimate is obtained for Hail site see in Figure 8 and Hail site see in Figure 9, where the maximum percentage error is found to be ±9%. The maximum percentage errors are ±12% for Al-Aqiq site see in Figure 7. Therefore, comparisons between the measured and calculated annual averages of Hd indicated that the all Saudi Arabia correlations (Models A, B and C) could be used for estimating the annual averages of Hd with excellent accuracy. The results in this research are similar results in the previous study [24] [35]. The average annual differences between measured and calculated values of diffuse solar energy Hd on horizontal at selected sites in the present research are summarizes in Figure 10, the results obtained are a good
Figure 10. The average annual differences between measured and calculated values of diffuse solar energy on horizontal at selected sites in the present research.
agreement is clear.
4. Conclusion
The available measured data of global solar energy and diffuse solar energy for the three selected locations in Saudi Arabia during the period time from 1980 to 2019 are used to develop the model empirical correlations as a function of sunshine duration (S/So) and clearness index (Kt) for predicted horizontal diffuse radiation energy at these sites. The values of regression constants a, b, c and d for all models are different values according to geographical information of site to another one during the period time in the present research, with the stander error estimate (S.E.). The correlation coefficients are higher than 0.86 for all models except for a few models. Models (17), (20) and (23) are given the higher values of correlation coefficients 0.935, 0.962 and 0.974 in Al-Aqiq site respectively, while Models (16), (20) and (21) are given 0.957, 0.936 and 0.942 in Hail site respectively, also Models (17), (20) and (21) are given 0.957, 948 and 965 in Dammam site respectively. For statistical indicators, the values of MBE for all models are given some negative and other positive, this is due to the correlations are overestimated and underestimated values of diffuse fraction Hd. The low values of the RMSE for all correlations indicate good agreement between the measured and calculated values of Hd. The negative values of MPE for all models show that for all locations, the proposed correlations slightly overestimate Hd. The absolute values of MPE never reach 1.35%. This means that very good agreement between the measured and calculated data and a good fitting exists between the monthly average daily values of the diffuse fraction (Hd/H)and (Kt) or (S/So) and between the diffuse transmittance (Hd/H0) and (S/So) for all selected locations in the present research. The first order of correlation modeling between the diffuse solar fraction and the clearness index has been proposed. In addition, the second order of correlations modeling between the diffuse solar fraction or the diffuse transmittance and the maximum possible number of sunshine hours are found to fit the measured data. The comparisons between the measured and calculated annual averages of Hd indicated that all Saudi Arabia correlations (Models A, B and C) could be used for estimating the annual averages of Hd with excellent accuracy. The average annual differences between measured and calculated values of diffuse solar energy Hd on horizontal at selected sites in the present research are summarized. All Saudi Arabia modeling correlations (Models A, B and C) can be used for predicting of diffuse solar energy Hd during the period time in the present research at the present locations as well as at any Saudi Arabia site with a reasonable accuracy. Therefore, the proposed modeling correlations can be used to accurately for estimating the annual averages of horizontal diffuse solar energy, which helps in predicted of the long-term performance of the various solar energy devices.