Next Article in Journal
Acknowledgement to Reviewers of Applied Sciences in 2018
Next Article in Special Issue
Optimal Energy Routing Design in Energy Internet with Multiple Energy Routing Centers Using Artificial Neural Network-Based Reinforcement Learning Method
Previous Article in Journal
Laboratory and Field Experiments on the Effect of Vinyl Acetate Polymer-Reinforced Soil
Previous Article in Special Issue
Analysis and Impact Evaluation of Missing Data Imputation in Day-ahead PV Generation Forecasting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Applications of ANN to Solar Radiation Estimation and Forecasting for Energy Applications

1
Renewable Energy Department, University of Corsica, UMR CNRS 6134, Route des Sanguinaires, 20000 Ajaccio, France
2
Castelluccio Hospital, Radiotherapy Unit, BP 85, 20177 Ajaccio, France
3
Laboratory Physical and Mathematical Engineering for Energy, Environment and Building, University of Reunion Island, 15 Avenue René Cassin, BP, 97715 Saint-Denis CEDEX, France
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(1), 209; https://doi.org/10.3390/app9010209
Submission received: 4 December 2018 / Revised: 27 December 2018 / Accepted: 31 December 2018 / Published: 8 January 2019
(This article belongs to the Special Issue Applications of Artificial Neural Networks for Energy Systems)

Abstract

:
In solar energy, the knowledge of solar radiation is very important for the integration of energy systems in building or electrical networks. Global horizontal irradiation (GHI) data are rarely measured over the world, thus an artificial neural network (ANN) model was built to calculate this data from more available ones. For the estimation of 5-min GHI, the normalized root mean square error (nRMSE) of the 6-inputs model is 19.35%. As solar collectors are often tilted, a second ANN model was developed to transform GHI into global tilted irradiation (GTI), a difficult task due to the anisotropy of scattering phenomena in the atmosphere. The GTI calculation from GHI was realized with an nRMSE around 8% for the optimal configuration. These two models estimate solar data at time, t, from other data measured at the same time, t. For an optimal management of energy, the development of forecasting tools is crucial because it allows anticipation of the production/consumption balance; thus, ANN models were developed to forecast hourly direct normal (DNI) and GHI irradiations for a time horizon from one hour (h+1) to six hours (h+6). The forecasting of hourly solar irradiation from h+1 to h+6 using ANN was realized with an nRMSE from 22.57% for h+1 to 34.85% for h+6 for GHI and from 38.23% for h+1 to 61.88% for h+6 for DNI.

1. Introduction

Solar thermal or electrical systems require high quality solar radiation measurement instruments in order to accurately measure solar energy received on the plant. Poor quality data or too short data series can generate errors in plant design, performance, and production forecasting, negatively impacting return on investment.
Unfortunately, measures of solar radiation are sparse and inaccurate over the world; there are still large areas without any solar radiation observations [1]. Investments and maintenance costs for each measurement site are not negligible and even in industrialized countries, the national network often consists in a relatively small number of solar radiation stations [2]; and the measurement quality varies from a network to another, often by lack of maintenance and calibration. The measuring devices’ price is an important part of the process cost of collecting solar data, especially for non-profit institutions, such as schools or universities [2].
The amount of meteorological stations measuring solar irradiance through the world is difficult to count because various sources give different information [3]. Even so, only 1000 continental stations around the world are measuring solar radiation [4].
In these conditions, it is interesting to look for some relations between sparse solar radiation data and more measured meteorological data as temperature, humidity, and wind speed. Satellite observations are used for determining the solar irradiation on the ground, with a good accuracy, but the time step of the estimated data is relatively large (minimum hourly).
A bibliographical study showed [5] that artificial neural networks (ANN) were developed between meteorological parameters and solar irradiations, but generally only for averaged values of solar irradiations (on a monthly or annual basis); today, data with a short time granularity (minute, 5-min) must be known for interesting applications in solar energy. In the first part of this paper, 10 meteorological inputs (measured or calculated) are available and then 210 − 1 = 1023 combinations of input data are possible, from 10 ANN with only one input to one ANN with 10 inputs. All these combinations are tested and the best configurations are discussed.
The solar panels are rarely fixed horizontally and the global tilted irradiation (GTI) is rarely measured and must be estimated from global horizontal irradiation (GHI), more often available; using pyranometers with various inclinations is costly and their maintenance is constraining. Thus, it is useful to develop accurate methods for determining GTI from only GHI. This objective is difficult to reach with conventional physical relations because the sky anisotropy makes the modeling of the sky diffuse radiation difficult [6,7,8]. ANN methods generally realize the same conversion with an improved adequacy (generally, at an hourly scale), and they outperform [9,10,11,12,13] the traditional methods due to the inherent non-linearity in solar radiation data. In this work, an ANN model is developed and optimized for “tilting” 5-min time step solar data.
The stochastic and intermittent behavior of the solar resources poses numerous problems for the electricity grid operator [5] and limits the future development of the phovoltaic (PV) and concentrated solar power (CSP) plants. To improve the integration of such systems, the solution consists in introducing energy storages and developing smart grids as well as implementing production and consumption forecasting.
The forecasting of the output power of solar energy systems is required for good operation of the power grid and for optimal management of electrical flows [3]. It is essential to estimate the energy reserves, to schedule the power systems, to optimally manage the storage, and to trade in the electricity market [14,15,16]. Thus, predicted and anticipated events are easier to manage. Electricity must be produced by CSP and/or PV plants; the first ones convert the direct normal irradiation (DNI) into heat through focusing receivers and PV ones enable direct conversion of GHI into electricity through semiconductor devices [3,17]. The literature shows that the most efficient methods for a forecast at a short time horizon from one hour (h+1) to six hours (h+6) are time series analysis and artificial intelligence methods [18]. If a large literature exists about the GHI forecasting [3,15,17,18,19], this literature is poorer concerning DNI being more difficult to predict because its variations are deeper and more frequent [20,21]. ANN predictive models are implemented to forecast hourly GHI and DNI from h+1 to h+6.
Section 2 presents the data used in this paper (Bouzareah, Algeria for estimation purposes and Odeillo, France for forecasting purposes), the preprocessing used on this data, and the calculated error metrics to estimate the accuracy of ANN models.
Section 3 gives some information on the ANN implementation for estimation and forecasting.
In Section 4 and Section 5, the main results on the estimation of 5 min-GHI from other meteorological parameters and of the 5 min-GTI from GHI measurements are presented.
Section 6 shows the results of the ANN forecasting of hourly GHI and DNI for a time horizon from h+1 to h+6.

2. Meteorological Stations and Data

2.1. The Meteorological Stations

The two first studies were realized using meteorological data measured in the meteorological station belonging to the Renewable Energies Development Centre (CDER) located in Bouzareah near Algiers (latitude: 36.8° N; longitude: 3.17° E) at an altitude of 347 m. The site is characterized by a Mediterranean climate with dry and hot summers and damp and cool winters. The data were measured each second and stored each 5 min from April 2011 to April 2013 (24 months of 5-min data). The measured data and the calculated astronomical data (horizontal extraterrestrial irradiation (EHI), solar declination, δ, and zenith angle, θ z [22]), are presented in Table 1.
The tilted solar data were measured for a 36.8° tilt angle equal to the latitude of Bouzareah (optimal angle for a maximum annual irradiation). The data basis contains 12 5-min parameters, 9 measured and 3 calculated. Each data was previously verified in order to extract outliers.
The forecasting work was realized from GHI and DNI data provided by the PROMES laboratory (CNRS UPR 8521) located in the south of France in Odeillo (Pyrénées Orientales, France, 42°29 N, 2°01 E, 1550 m asl), the station is located in the mountains, at about 100 km from the Mediterranean sea and presents an often high nebulosity. The solar data are measured and stored with a 1 min time granularity. This meteorological station is in altitude, the climate is very perturbed, the rainfall continues to be present during the driest months, the variability of solar radiation is high, and thus its forecasting is more difficult to realize. Two years of hourly data were available i.e., 17 520 data, for both GHI and DNI.

2.2. Cleaning and Preprocessing

For Bouzareah, each 5-min data were first verified to extract outliers or missing data. Then, the data during which the sun rises or sets were deleted because the mask effect of the environment and the no-reliable response of pyranometers at a high zenith angle (cosine effect) introduced some errors. Thus, over the 2 years, 75674 validated 5-min data were available for each parameter.
For the Odeillo’s data, an automatic quality control used in the frame of the GEOSS project (Group on Earth Observation System of System) [23] was applied. Before introducing the solar data into the machine learning process, the data were cleaned and filtered.
For forecasting purposes, it is common to filter out the data to remove night hours and to conserve them only between sunrise and sunset. As for Bouzareah, the data near sunset and sunrise are sources of errors and a pre-processing operation was applied based on the solar elevation: Solar radiation data for which the solar elevation is lower than 10° were removed [15,24]. Two years of hourly data were used in this study. After cleaning and filtering, the total number of hourly data for each solar component (GHI and DNI) was 10559 (about 60% of the data were not used (2% for outliers’ data and 58% for sun elevation less than 10°)). These solar data were then transformed into stationary data by a method described in Section 6.

2.3. Statistical Index for Accuracy Evaluation

There are no well-defined error metrics standards, which makes the forecasting and estimation methods difficult to compare [25]. A benchmarking exercise was realized within the framework of the European Actions Weather Intelligence for Renewable Energies (WIRE) [26], with the objective to evaluate the state of the art concerning models’ performances for short term renewable energy forecasting. They concluded that: “More work using more test cases, data and models needs to be performed in order to achieve a global overview of all possible situations. Test cases located all over Europe, the US and other relevant countries should be considered, trying to represent most of the possible meteorological conditions”.
In this paper, these five error metrics were used:
-
The mean absolute error (MAE) defined by:
M A E = 1 N × i = 1 N | X ^ t i X t i |
X ^ t i being the forecasted outputs (or predicted values), X t i the observed data, and N the number of observations.
-
The root mean square error (RMSE), more sensitive to important forecast errors, and hence suitable for applications where small errors are more tolerable than larger ones, as in utility applications. It is probably the reliability factor that is the most widely used:
R M S E = 1 N × i = 1 N ( X ^ t i X t i ) 2
-
The mean bias error (MBE), mainly used to estimate the bias of the model:
M B E = 1 N × i = 1 N ( X ^ t i X t i )
These errors were then normalized, and the mean value of irradiation is generally used as the reference:
n R M S E = R M S E X ¯
n M A E = M A E X ¯
with X ¯ the average value of X calculated on the N data.

3. ANN Method and Implementation

3.1. General Description of ANN Structure

An ANN [27] is a modelling tool able to find complex relationships between inputs and outputs. It is considered as “intelligent” because it works as a human brain:
-
A neural network acquires knowledge through learning; and
-
a neural network’s knowledge is stored within inter-neuron connection strengths known as synaptic weights.
ANNs can represent both linear and non-linear relationships and have their ability to learn these relationships directly from the measured data. Among the various ANN methods, the multilayer perceptron (MLP) using feed-forward back-propagation is often used for empirical estimation in general [28], and in particular [19,29,30], for the estimation of solar radiation. An ANN has a parallel structure with an input layer receiving data, an output layer sending the computed data, and one or several hidden layers lying between the input and output layers as illustrated in Figure 1.
A neuron receives input from other neurons and/or an input data, which represents an external source. In the feed-forward propagation MLP configuration, this connection is unidirectional. Each input, x k , j , has an associated weight, w k , j (related to the j-th neuron among p of the k-th layers), which can be modified during the learning phase. The weighted sum, j = 1 p w k , j x k , j , is called the net input to unit j, often written net. The unit computes a function, f, of this weighted sum and is called the activation or transfer function; this function, f, produces an output, O, of a neuron if this sum exceeds a given threshold denoted biases. A bibliographical study conducted to use as activation transfer functions a sigmoid one for the hidden layers and a linear function for the output layer. For the jth neuron of the layer (k+1), O is given by [31]:
O = f ( n e t ) = f ( j = 1 p w k , j x k , j )
This output is then distributed to other neurons as inputs.

3.2. ANN Implementation for GHI and GTI Irradiation

Several steps are applied in view to find the optimized MLP:
-
Choice of the network size (number of hidden layers and hidden nodes per layer): Too small a number of hidden nodes does not allow good learning, but an oversized number increases the training time with a marginal improvement [32,33], need more data and the ANN can be over-trained. In accordance with the principle of parsimony and with the literature [34], only one hidden layer was used.
-
Determination of the optimal number of neurons in the hidden layer: It is realized in testing various configurations and calculating the adequacy. Some empirical rules exist, but their efficiency is not really proven: The number of hidden neurons equal to the inputs number [35], to 75% of it [36], to the square root of the product of the number of inputs and outputs [37]. Here, the number of hidden neurons was taken between 1 and the number of inputs +1. Each MLP architecture was trained 8 times per architecture in order to avoid random effects.
-
Learning (or training) process: It consists in modifying the weights until the gap between the actual and simulated outputs reaches a desired accuracy. The Levenberg–Marquardt learning algorithm (LM) was used as in most studies. Another preprocess called k-fold sampling was used with the dataset [38,39], this cross-validation is a statistical method used to estimate the skill of machine learning models. It is commonly used in applied machine learning to compare and select a model for a given predictive modeling problem, which has a lower bias than other methods [40,41]: It consists in dividing randomly the data set into a training data set (80%) and a test data set (20%), the training set and the test set are different for each fold; this process is repeated k times and the value of the reliability metrics given in this paper are the average value on the k-fold. Here, k is taken equal to 10. Thus, the results are independent of the set of data used for the training; using only one data set (with its own statistical particularities) can reduce the robustness of the conclusions.
Other information specific to each study will be given in the corresponding chapter.

4. Estimation of GHI from Other Meteorological Data

4.1. Method

In this study, the data used were measured in Bouzareah (see Section 2.1). Figure 2 shows the ANN structure with all the available inputs.
The sunshine duration is the time expressed in minutes during which the solar irradiance exceeds 120 W/m2. It is strongly linked with the solar radiation (Angstrom relation for example).
The number of inputs, 10 (7 measured, 3 calculated), makes the optimization of the MLP structure long and arduous: With 10 inputs, 210 − 1 = 1023 combinations of input data are possible.
The choice of the best inputs combination is a prerequisite stage because the parsimony is a basic principle in ANN elaboration, essential for its generalization. Some of the variables bring little information, sometimes no information at all, some of them are redundant, even worse they reduce the model performance. Moreover, an increase of the input number is accompanied by an increase of the hidden neurons and of the calculation time.
The Pearson’s correlation coefficient between each input parameter and the output is determined before using an exhaustive selection (testing 1023 architectures).

4.2. Relationship between Input and Output Data

The Pearson’s correlation coefficient (R) between input variables and GHI and between inputs variables themselves were calculated. When the absolute value of R is near 1, there is a high degree of linear correlation between the two variables; if R = 0, there is no linear correlation, but other relation types can exist.
Computing R between input variables allows an estimation of whether the inputs are redundant and interdependent. The first objective is to rank the statistical linear dependences between the inputs and output; for a large sample of data, the R threshold from which there is a significant link between parameters is very low. Table 2 shows the values of the Pearson correlation coefficient between meteorological variables.
We mainly see:
  • The only high correlation is between GHI and S (69%), the other are just weak correlations;
  • The ranking of inputs from the R point of view (excepted S) is H (18%), EHI (13%), θ z (13%), and p (10%); and
  • A high value of R between inputs data, θ z , δ, and EHI (between 34% and 95%); it was obvious because θ z and δ are used in the calculation of EHI.
This preliminary study allows to have an idea about the link between variables, but only linear ones. The results were not significant enough to avoid an exhaustive study for all combinations; the Pearson coefficient allows estimation of only the linear dependency between the data while the MLP is a non-linear model (sigmoid activation function), thus it is not surprising that the analysis of the Pearson coefficient is not sufficient to customize an MLP.

4.3. Results

Figure 3 shows the average, minimum, and maximum values of nRMSE and the standard deviation versus the number of inputs.
The minimum nRMSE was obtained for the 10 inputs model with a value of 18.65% compared to 73.91% for the worst configuration (2 inputs: WD and WS). Table 3 presents the two best configurations for the same number of inputs; as an example, for 9 inputs, 10 combinations were possible, and only the two best ones (from an nRMSE point of view) are reported in Table 3. The models were classified in descending order of performance (ranking).
Some models with a lower number of inputs can be better than models with a higher number of inputs. The declination, δ, appeared very rarely, WD and Pr never appeared, S was always present, and T, H, EHI, and p were often present; these 5 inputs had a relatively good R with GHI (but only 6% for temperature).
Without S (strongly linked with GHI), the best nRMSE dropped to 32.07% compared with 18.65% for the best configuration with 10 inputs. Table 4 shows the results for the best configurations without S input.
Figure 4 shows the average, minimum, and maximum nRMSE values and the standard deviation versus the number of input data without S input.
The ANN reliability can be considered as correct, particularly when S was an input. Pr and WD had a low influence on GHI estimation (low correlation with GHI). The nRMSE of the 6-inputs model (T, H, p, WS, S, EHI) was 19.35% compared with the nRMSE of 10-inputs model with 18.65%; this combination had a good performance with a minimum of inputs.
Estimated 5-min GHI is plotted versus the measured 5-min GHI for four architectures in Figure 5.
-
ANN structure with 10 inputs;
-
ANN structure with six inputs;
-
ANN structure with nine inputs without sunshine duration; and
-
best ANN structure with five inputs (without sunshine duration).
Few differences appeared between the ANN structure with 10 and six inputs in term of reliability and data dispersion. Without S, it appears clearly a more important spread of data compared with the results obtained with ANN structures with S in the inputs set. The performances of the best ANN structures without S for 5-min data were correct with an nMAE between 28.5% and 31% and an nRMSE between 32% and 35%.
The presence of some meteorological inputs in the “best” configurations seems sometimes surprising as WD and WS; it is difficult to understand the physical relations between GHI and other meteorological parameters. One of the major criticisms that could be levelled at the ANN model is that it is a black box model, allowing it to find some relations between data as often difficult to interpret, and ANN is a data driven method.
Even if some estimated data are far away from the real values, we can consider that the performance of this model is satisfying because determining GHI with a time granularity of 5-min from other meteorological data is a very complex task (high variability phenomenon and anisotropic aspect); keeping in mind that such a method is generally applied only for daily average values [5].
A bibliographical study [5] was realized on ANN methods used for such an estimation of GHI from exogenous meteorological data and this study showed that:
-
For the estimation of monthly mean values of daily GHI, the nRMSE was between 4.07% and 9.4%, but the process, on average, monthly, allows smoothing of the anisotropic effects and sometimes linear relationships are sufficient to link GHI with other parameters.
-
For the estimation of the daily GHI, nRMSE around 6% and nMAE around 5% were found. The time granularity was much higher than in our work.
Note that the determination coefficient (R2) between measured and estimated data was between 0.86 and 0.95 for the four graphs in Figure 5.
The mean bias error (MBE) was also computed for Figure 5 and was equal to −0.08 Wh/m2, for 10 and six inputs, −1.14 Wh/m2 for nine inputs without S, and −1.87 Wh/m2 for five inputs; thus, all the ANN models slightly underestimated GHI.

5. Estimation of Tilted Global Irradiation (GTI) from Horizontal Global Irradiation (GHI)

In sizing or simulation software for solar systems, the solar collector inclination is introduced as an input and the horizontal solar data (generally hourly) collected from several meteorological stations are “tilted”. The accuracy and quality of GTI used as an input in these software have an impact on the reliability of the results.
It is difficult to develop a simple model for converting GHI into GTI [6] because the radiation received by a tilted plane includes the radiation reflected by the ground and scattered by the sky; this last component is difficult to estimate; when the collector is inclined, it sees only a part of the sky; moreover, the sky diffuse radiation depends on the inclination or orientation of the collector, on the elevation and azimuth of the sun, but also on the sky state with complex anisotropic effects [7,8]. The larger the time-step is, the more this anisotropy decreases (time-averaging and compensating effect) and tends towards an isotropic distribution; the shorter the time-step is, the more it is difficult to realize this conversion with good accuracy. The conversion of GHI to GTI is a complex issue often dealt with in the scientific literature [7,42,43,44,45,46].

5.1. Method

As in Section 4, the data used here were measured in Bouzareah.
Four data, among them GHI, were used as input:
-
The declination representing the position of the Earth from the Sun depending on the day number;
-
the zenith angle characterizing the sun position, which influences the quantity and the quality of the sun radiation; when the sun is high in the sky (low zenith angle), the solar radiation is maximal (in clear skies). Moreover, as the optical path is minimal, the incident radiation is less absorbed;
-
the extraterrestrial irradiation, EHI, used as a reference; depending on sky conditions, several values of GTI correspond to the same GHI. In diffuse radiation models, the clearness or diffuse index are often used to characterize the sky. When the clearness index is high, then the sky is clear and GHI is mainly composed of BNI.
According to the rules described in Section 3.2, the number of hidden neurons in only one hidden layer will vary from one to eight. The ANN structure is shown in Figure 6.

5.2. Results

The five error metrics are presented in Table 5 (calculated on the basis of eight runs (each run corresponds to a different random weight initializing). The first column contains the number of neurons in the hidden layer.
The nRMSE mean values and its corresponding standard deviation are presented in Figure 7 as an error-bar graph. An improvement appears until reaching four hidden neurons, then the nRMSE becomes almost constant and no improvement was observed. The dash points define the 95% confidence interval of the prediction errors (calculated based on eight runs per architecture), the triangles and squares are the minimum and maximum observed errors, respectively. We observed the same trend for the variation of the nMAE. The best configuration is encircled in red.
We conclude that an ANN with one hidden layer of four neurons is the best model. Moreover, it appears that the use of the azimuth does not provide any improvement. Thus, we will retain an ANN with four inputs and one hidden layer of four neurons, which have an average nRMSE of 8.81%, the best simulation with this ANN structure conduced to an nRMSE of 8.27%.
To illustrate the good reliability of this optimized MLP, a period of seven days unknown to the network was plotted with measured and calculated data in Figure 8.
A good relationship is observed between the modelled and measured data whatever the state of the sky is (clear, partially cloudy, and cloudy) because the nRMSE was under 10%, which is a good value for an nRMSE for such a short time step (5-min).

6. Forecasting of Hourly Direct Normal (DNI) and Horizontal Global (GHI) Irradiation for a Time Horizon from 1 h to 6 h

This forecasting work was realized from global horizontal (GHI) and normal beam (BNI) data measured at Odeillo, Pyrénées Orientales, located in the south of France (42°29 N, 2°01 E, 1550 m asl).

6.1. Stationnarization of Solar Data

Machine learning methods are efficient tools for forecasting time series with a stationary behavior. An MLP is a stationary model, which must use stationary data as input. To make solar irradiation time series stationary and to separate the climatic effects and the seasonal effects, the solar data are generally transformed in unitless variables called “clearness index”, and denoted kt; kt is the ratio of the solar radiation on the earth, GHI, to that outside the atmosphere, EHI, and defined by Equation (6) [35]:
k t = G H I E H I
It is the clearness index series, kt, that induces randomness, caused by the diversity of atmospheric components (dusts, aerosols, clouds motion, and humidity) on the solar irradiation measured at the Earth‘s surface.
Numerous studies have showed that EHI can be efficiently replaced by the clear sky solar irradiation [22] taking into account the climatic conditions of the meteorological site; thus the clearness index is replaced by the clear sky index, kg,cs, defined by:
k g , c s = G H I G H I C S
with G H I C S the global horizontal solar irradiation in clear sky conditions.
For DNI, a similar index, kd,CS [47,48], is defined by:
k d , c s = D N I D N I C S
Various models of clear sky solar irradiations are available in the literature, which differ from each other mainly in the inputs needed by each model [49]. Solar irradiance models by clear sky denoted in the following clear sky models used meteorological variables (as ozone layer thickness, precipitable water, optical aerosol depth, etc.) and used solar geometry (solar elevation and air mass), using radiative transfer models to consider the absorption and diffusion effects of solar radiation into the atmosphere [50,51]. The most widely used clear sky model is the Solis model developed by Mueller et al. [52] and simplified by Ineichen [53], the European Solar Radiation Atlas (ESRA) model [54], and the Reference Evaluation on Solar Transmittance 2 (REST2) model [55].
Thus, the simplified Solis clear sky model [53] was used here. It allowed calculations of the GHICS, and DNICS. This clear sky model was validated for each month by comparison with experimental solar radiation data measured in clear sky conditions. For illustration purposes, experimental and modelled solar irradiances by clear sky are plotted in Figure 9 for one day in April and in September.

6.2. Choice of the Number of Input Data

The purpose of Section 6 is to predict the future hourly solar irradiation (at different time horizons) based on the past observed data, i.e., mathematically:
X ^ t + h = f ( X t , X t 1 , X t 2 , X t 3 , , X t n )
A variable, X, with the symbol, ^ , represents a forecasted data; without this symbol, X is a measured data. The solar data at future time step (t+h), X ^ t + h is forecasted from the observed data X measured at the times (t, t − 1…, tn); thus, the first objective consists of determining the value of n, i.e., the dimension of the input matrix; to do it, an auto mutual information method [56,57,58] was used. The auto mutual information is a property of the time series, it depends on each dataset and is characteristic of the degree of statistical dependence between X t + h and X t i with 0 ≤ in. It is a dimensionless quantity with units of bits. High mutual information indicates a large reduction in uncertainty about one random variable, X t + h , given knowledge of another ( X t , X t 1 , X t 2 , X t 3 , , ) . The auto-mutual information method showed that the number of inputs (value of n in Equation (9)) for predicting GHI is six and for DNI, it is seven.

6.3. Methods

A first forecasting method, a naïve model, easy to implement and requiring no training step, i.e., no historical data set, was used as a reference model to compare it with more sophisticated models in terms of accuracy. It allowed us to see the improvement due to the use of the ANN forecaster.
The persistence model, the simplest forecasting model, assumes that the future value is identical to the previous one (Equation (10)). The persistence forecast accuracy decreases significantly with the forecasting horizon [59].
X ^ t + h = X t   with   X = GHI   or   DNI
The smart persistence (SP) is an improved version of the persistence one taking into account the diurnal solar cycle: The clear sky solar radiation profile over the day was used [41]:
X ^ t + h = X t · X C S t + h X C S t   with   X = GHI   or   DNI
X C S t is GHI or DNI for a clear sky condition calculated at time t. This smart persistence model was applied in this paper and used mainly as a reference model.
For the ANN model, as explained in Section 6.1, the clear sky index will be forecasted because it is a stationary series. Once the clear sky index is forecasted, the value of the forecasted solar irradiation, ( G H I ^ t + h ,   o r   D N I ^ t + h ) is obtained by multiplying it by the calculated clear sky irradiation ( G H I C S t + h   o r   D N I C S t + h ) .

6.4. Results for GHI

Table 6 gives the values of the error metrics calculated on the test data set (RMSE, MAE, and MBE are given in Wh.m−2).
As the ranking of the model is always identical from a RMSE point of view or a MAE point of view, we only present in Figure 10 the results in terms of RMSE and nRMSE expressed in percentage.
The smart persistence, a naive model, was used as a reference. This model has a good RMSE and MAE for a time horizon, h+1, but its performances decrease rapidly with the time horizon. The gap in term of performances between ANN and SP increases with the time horizon.

6.5. Results for DNI

Table 7 gives the values of the performance metrics computed on the test data set (RMSE, MAE, and MBE are given in Wh.m−2) for DNI.
The results in terms of RMSE and nRMSE are presented in Figure 11 for BNI.
The DNI forecasting is more difficult and the models’ performances were less satisfying than with GHI particularly because DNI was more sensitive to meteorological conditions and because its variation was more rapid and of a greater magnitude than for GHI.
Some differences are noted in terms of ranking according to the metric used (nRMSE or nMAE), the nRMSE gives more importance to large gaps between predicted and measured data and generally the forecasting models were better compared in term of nRMSE than nMAE.

6.6. Comparison between GHI and DNI Forecasts

It is impossible to compare the performances of the models according to the solar component in terms of the absolute value of RMSE (or MAE) because the daily curve of GHI and DNI are very different in term of amplitude and form. Thus, we plotted in Figure 12 a comparison of the performances in term of nRMSE because it is the most common error metric used in the solar radiation prediction; in the case of Figure 11, as we compared two different solar radiations (GHI and DNI) with a different scale, the normalized value of RMSE seems to be the most adapted metric.
As previously underlined, GHI was forecasted with a better accuracy compared with DNI. It is probably due to the fact that in GHI, the two components, diffuse one and beam one, have compensating effects (when diffuse increases, beam decreases) and the variation rate of GHI is less rapid than for DNI. With SP and ANN methods, DNI is predicted with an nRMSE nearly twice as high than for GHI, but this difference was reduced when the forecast horizon increased and for (h+6), the accuracy for DNI prediction was the same than for GHI prediction.
Antonanzas et al. [50] reviewed the intra-day ahead forecast performances for PV production (using GHI as renewable resource) using different numerical prediction models. Various error metrics were used and calculated according to different definitions, moreover, the forecasting methods were applied in different meteorological stations; thus, it is very difficult to make a comparison of our results with the literature in these conditions.
A short bibliographical study [60] on DNI forecasting concludes that the DNI forecasting is obtained with a lower accuracy than for GHI forecasting and that only a small number of articles are written on the DNI forecasting at short time horizons as confirmed by Law et al. [20].

7. Conclusions

In this paper, the efficiency of the ANN method was shown for the estimation and the forecasting of solar irradiation.
Successively, several ANN models were developed:
-
An ANN model to determine the 5-min GHI from more available meteorological data (a maximum of seven measured meteorological data). The nRMSE of the optimized 6-inputs model was 19.35%.
-
An ANN model to compute the 5-min GTI from GHI with an accuracy around 8% for the optimal configuration.
-
At last, a forecasting ANN tool was developed to predict hourly DNI and GHI i for a time horizon from h+1 to h+6. The forecasting of hourly solar irradiation from h+1 to h+6 using ANN was realized with an nRMSE from 22.57% for h+1 to 34.85% for h+6 for GHI and an nRMSE from 38.23% for h+1 to 61.88% for h+6 for DNI.
ANN methods are very promising even if new artificial intelligence methods are in development, such as regression trees and random forest.

Author Contributions

Software: A.F. and J.L.D.; Methodology: C.V., G.N. and M.L.N.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Atwater, M.A.; Ball, J.T. A numerical solar radiation model based on standard meteorological observations. Sol. Energy 1978, 21, 163–170. [Google Scholar] [CrossRef]
  2. Hidalgo, F.G.; Martinez, R.F.; Vidal, E.F. Design of a Low-Cost Sensor for Solar Irradiance. Available online: http://oceanoptics.com (accessed on 25 November 2018).
  3. Paulescu, M.; Paulescu, E.; Gravila, P.; Badescu, V. Weather Modeling and Forecasting of PV Systems Operation, Green Energy and Technology; Ó Springer-Verlag: London, UK, 2013. [Google Scholar]
  4. World Radiation Data Center (WRDC). WRDC Online Archive, National Renewable Energy Laboratory, US Department of Energy. 2012. Available online: https://www.re3data.org (accessed on 3 May 2017).
  5. Dahmani, K.; Notton, G.; Voyant, C.; Dizene, R.; Nivet, M.L.; Paoli, C.; Tamas, W. Multilayer Perceptron Approach for estimating 5-min and hourly horizontal global radiation from exogenous meteorological data in locations without solar measurements. Renew. Energy 2016, 90, 267–282. [Google Scholar] [CrossRef]
  6. Behr, H.D. Solar radiation on tilted south-oriented surfaces: Validation of transfer-models. Sol. Energy 1997, 61, 399–413. [Google Scholar] [CrossRef]
  7. Robledo, L.; Soler, A. Modelling irradiance on inclined planes with an anisotropic model. Energy 1998, 23, 193–201. [Google Scholar] [CrossRef]
  8. Padovan, A.; Del Col, D. Measurement and modeling of solar irradiance components on horizontal and tilted planes. Sol. Energy 2010, 84, 2068–2084. [Google Scholar] [CrossRef]
  9. Hontoria, L.; Riesco, J.; Zufiria, P.; Aguilera, J. Improved generation of hourly solar radiation artificial series using neural networks. In Proceedings of the Engineering Applications of Neural Networks (EANN’99), Warsaw, Poland, 13–15 September 1999. [Google Scholar]
  10. Hontoria, L.; Aguilera, J.; Zufiria, P. Generation of hourly irradiation synthetic series using the neural network multilayer perceptron. Sol. Energy 2002, 72, 441–446. [Google Scholar] [CrossRef]
  11. Jiang, Y. Computation of monthly mean daily global solar radiation in China using artificial neural networks and comparison with other empirical models. Energy 2009, 34, 1276–1283. [Google Scholar] [CrossRef]
  12. Notton, G.; Paoli, C.; Ivanova, L.; Vasileva, S.; Nivet, M.L. Neural network approach to estimate 10-min solar global irradiation values on tilted plane. Renew. Energy 2013, 50, 576–584. [Google Scholar] [CrossRef]
  13. Elminir, H.K.; Azzam, Y.A.; Younes, F.I. Prediction of hourly and daily diffuse fraction using neural network, as compared to linear regression models. Energy 2007, 3, 1513–1523. [Google Scholar] [CrossRef]
  14. Heinemann, D.; Lorenz, E.; Girodo, M. Forecasting of solar Radiation. In Solar Energy Resource Management for Electricity Generation from Local Level to Global Scale; Nova Science Publishers: New York, NY, USA, 2006. [Google Scholar]
  15. Lauret, P.; Voyant, C.; Soubdhan, T.; David, M.; Poggi, P. A benchmarking of machine learning techniques for solar radiation forecasting in an insular context. Sol. Energy 2015, 112, 446–457. [Google Scholar] [CrossRef]
  16. McCandless, T.C.; Haupt, S.E.; Young, G.S. Short term solar radiation forecast using weather regime dependent artificial intelligence techniques. In Proceedings of the 2th Conference on Artificial and Computational Intelligence and its Applications to the Environmental Sciences, Atlanta, GA, USA, 2–6 February 2014. [Google Scholar]
  17. COST. Weather Intelligence for Renewable Energies (WIRE). Current State Report No.ES1002. 2012. Available online: www.cost.eu/actions/ES1002 (accessed on 8 May 2017).
  18. Diagne, M.; David, M.; Lauret, P.; Boland, J.; Schmutz, N. Review of solar irradiance forecasting methods and a proposition for small-scale insular grids. Renew. Sustain. Energy Rev. 2013, 27, 65–76. [Google Scholar] [CrossRef]
  19. Kalogirou, S.A.; Şencan, A. Artificial intelligence techniques in solar energy applications. In Solar Collectors and Panels, Theory and Applications; Manyala, R., Ed.; Intechopen Publisher: London, UK, 2010; ISBN 978-953-307-142-8. [Google Scholar]
  20. Law, E.W.; Prasad, A.A.; Kay, M.; Taylor, R.A. Direct normal irradiance forecasting and its application to concentrated solar thermal output forecasting—A review. Sol. Energy 2014, 108, 287–307. [Google Scholar] [CrossRef]
  21. Ghofrani, M.; Ghayekhloo, M.; Azimi, R. A novel soft computing framework for solar radiation forecasting. Appl. Soft Comput. 2016, 48, 207–216. [Google Scholar] [CrossRef]
  22. Iqbal, M. An Introduction to Solar Radiation; Academic Press: Don Mills, ON, Canada, 1983; ISBN 0-12-373752-4. [Google Scholar]
  23. Global Earth Observation System of Systems (GEOSS). Available online: www.earthobservations.org/geoss.php (accessed on 3 May 2017).
  24. David, M.; Ramahatana, F.; Trombe, P.J.; Lauret, P. Probabilistic forecasting of the solar irradiance with recursive ARMA and GARCH models. Sol. Energy 2016, 133, 55–72. [Google Scholar] [CrossRef] [Green Version]
  25. Notton, G.; Voyant, C. Forecasting of Intermittent Solar Energy Resource. In Advances in Renewable Energies and Power Technologies; Yahyaoui, I., Ed.; Elsevier Science: Amsterdam, The Netherlands, 2018; pp. 77–109. ISBN 978-012-8131855. [Google Scholar]
  26. Sperati, S.; Alessandrini, S.; Pinson, P.; Kariniotakis, G. The “Weather Intelligence for Renewable Energies” Benchmarking Exercise on Short-Term Forecasting of Wind and Solar Power Generation. Energies 2015, 8, 9594–9619. [Google Scholar] [CrossRef] [Green Version]
  27. Haykin, S. Neural Networks: A Comprehensive Foundation, 2nd ed.; Prentice-Hall: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  28. Yildiz, N. Layered feedforward neural network is relevant to empirical physical formula construction: A theoretical analysis and some simulation results. Phys. Lett. A 2005, 345, 69–87. [Google Scholar] [CrossRef]
  29. Mellit, A.; Benghanem, M.; Hadj Arab, A.; Guessoum, A. An adaptive artificial neural network model for sizing stand-alone photovoltaic systems: Application for isolated sites in Algeria. Renew. Energy 2005, 30, 1501–1524. [Google Scholar] [CrossRef]
  30. Mellit, A.; Pavan, A.M. A 24-h forecast of solar irradiance using artificial neural network: Application for performance prediction of a grid-connected PV plant at Trieste, Italy. Sol. Energy 2010, 84, 807–821. [Google Scholar] [CrossRef]
  31. Abraham, A. Artificial Neural Networks. Handbook for Measurement Systems Design; Sydenham, P., Thorn, R., Eds.; John Wiley and Sons Ltd.: London, UK, 2005; pp. 901–908. ISBN 0-470-02143-8. [Google Scholar]
  32. Krishnaiah, T.; Srinivasa Rao, S.; Madhumurthy, K.; Reddy, K.S. Neural network approach for modelling global solar radiation. J. Appl. Sci. Res. 2007, 3, 1105–1111. [Google Scholar]
  33. Alam, S.; Kaushik, S.C.; Garg, S.N. Assessment of diffuse solar energy under general sky condition using artificial neural network. Appl. Energy 2009, 86, 554–564. [Google Scholar] [CrossRef]
  34. Cybenko, G. Approximation by superposition of sigmoidal function. Math. Control Signals Syst. 1989, 2, 303–314. [Google Scholar] [CrossRef]
  35. Wierenga, B.; Kluytmans, J. Neural nets versus marketing models in time series analysis: A simulation study. In Proceedings of the 23th Annual Conference of the European Marketing Academy, Maastricht, The Netherlands, 17–20 May 1994; pp. 1139–1153. [Google Scholar]
  36. Venugopal, V.; Baets, W. Neural networks and statistical techniques in marketing research: A conceptual comparison. Mark. Intell. Plan. 1994, 12, 30–38. [Google Scholar] [CrossRef]
  37. Shepard, R.N. Neural nets for generalization and classification: Comment on Staddon and Reid. Psychol. Rev. 1990, 97, 579–580. [Google Scholar] [CrossRef] [PubMed]
  38. Wiens, T.S.; Dale, B.C.; Boyce, M.S.; Kershaw, G.P. Three way k-fold cross-validation of resource selection functions. Ecol. Model. 2008, 212, 244–255. [Google Scholar] [CrossRef]
  39. Kuhn, M.; Johnson, K. Applied Predictive Modelling; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  40. Cigizoglu, H.K.; Kişi, Ö. Flow prediction by three back propagation techniques using k-fold partitioning of neural network training data. Hydrol. Res. 2005, 36, 49–64. [Google Scholar] [CrossRef]
  41. Voyant, C.; Soubdhan, T.; Lauret, P.; David, M.; Muselli, M. Statistical parameters as a means to a priori assess the accuracy of solar forecasting models. Energy 2015, 90, 671–679. [Google Scholar] [CrossRef] [Green Version]
  42. Wenxian, L.; Wenfeng, G.; Shaoxuan, P.; Enrong, L. Ratios of global radiation on a tilted to horizontal surface for Yunnan province. China. Energy 1995, 20, 723–728. [Google Scholar] [CrossRef]
  43. Li, D.H.W.; Lau, C.C.S.; Lam, J.C. Predicting daylight illuminance on inclined surfaces using sky luminance data. Energy 2005, 30, 1649–1665. [Google Scholar] [CrossRef]
  44. Cheng, C.L.; Chan, C.Y.; Chen, C.L. An empirical approach to estimating monthly radiation on south-facing tilted planes for building application. Energy 2006, 31, 2940–2957. [Google Scholar] [CrossRef]
  45. De Rosa, A.; Ferraro, V.; Kaliakatsos, D.; Marinelli, V. Calculating diffuse illuminance on vertical surfaces in different sky conditions. Energy 2008, 33, 1703–1710. [Google Scholar] [CrossRef]
  46. Pandey, C.K.; Katiyar, A.K. A note on diffuse solar radiation on a tilted surface. Energy 2009, 34, 1764–1769. [Google Scholar] [CrossRef]
  47. Linares-Rodríguez, A.; Ruiz-Arias, J.A.; Pozo-Vázquez, D.; Tovar-Pescador, J. Generation of synthetic daily global solar radiation data based on ERA-Interim reanalysis and artificial neural networks. Energy 2011, 36, 5356–5365. [Google Scholar] [CrossRef]
  48. Kaur, A.; Nonnenmacher, L.; Pedro, H.T.C.; Coimbra, F.M. Benefits of solar forecasting for energy imbalance markets. Renew. Energy 2016, 86, 819–830. [Google Scholar] [CrossRef]
  49. Badescu, V.; Gueymard, C.A.; Cheval, S.; Oprea, C.; Baciu, M.; Dumitrescu, A.; Iacobescu, F.; Milos, J.; Rada, C. Computing global and diffuse solar hourly irradiation on clear sky. Review and testing of 54 models. Renew. Sustain. Energy Rev. 2012, 16, 1636–1656. [Google Scholar] [CrossRef]
  50. Antonanzas, J.; Osorio, N.; Escobar, R.; Urraca, R.; Martinez-de-Pison, F.J.; Antonanzas-Torres, F. Review of photovoltaic power forecasting. Sol. Energy 2016, 136, 78–111. [Google Scholar] [CrossRef]
  51. Ineichen, P. Validation of models that estimate the clear sky global and beam solar irradiance. Sol. Energy 2016, 132, 332–344. [Google Scholar] [CrossRef] [Green Version]
  52. Mueller, R.; Dagestad, K.; Ineichen, P.; Schroedter-Homscheidt, M.; Cros, S.; Dumortier, D. Rethinking satellite-based solar irradiance modeling: The SOLIS clear-sky module. Remote Sens. Environ. 2004, 91, 160–174. [Google Scholar] [CrossRef]
  53. Ineichen, P. A broadband simplified version of the Solis clear sky model. Sol. Energy 2008, 82, 758–762. [Google Scholar] [CrossRef]
  54. Rigollier, C.; Bauer, O.; Wald, L. On the clear sky model of the ESRA—European Solar Radiation Atlas—with respect to the Heliosat method. Sol. Energy 2000, 68, 33–48. [Google Scholar] [CrossRef]
  55. Gueymard, C.A. REST2: High-performance solar radiation model for cloudless-sky irradiance, illuminance, and photosynthetically active radiation—Validation with a benchmark dataset. Sol. Energy 2008, 82, 272–285. [Google Scholar] [CrossRef]
  56. Huang, D.; Chow, T.W.S. Effective feature selection scheme using mutual information. Neurocomputing 2005, 63, 325–343. [Google Scholar] [CrossRef]
  57. Jiang, A.H.; Huang, X.C.; Zhang, Z.H.; Li, J.; Zhang, Z.Y.; Hua, H.X. Mutual information algorithms. Mech. Syst. Signal Process. 2010, 24, 2947–2960. [Google Scholar] [CrossRef]
  58. Parviz, R.K.; Nasser, M.; Motlagh, M.R.J. Mutual Information Based Input Variable Selection Algorithm and Wavelet Neural Network for Time Series Prediction. In Proceedings of the International Conference on Artificial Neural Networks (ICANN 2008), Prague, Czech Republic, 3–6 September 2008; Kůrková, V., Neruda, R., Koutník, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2008; pp. 798–807. [Google Scholar]
  59. Huang, R.; Huang, T.; Gadh, R.; Li, N. Solar generation prediction using the ARMA model in a laboratory-level micro-grid. In Proceedings of the IEEE Third International Conference on Smart Grid Communications (SmartGridComm), Tainan, Taiwan, 5–8 November 2012; pp. 528–533. [Google Scholar]
  60. Benali, L.; Notton, G.; Fouilloy, A.; Voyant, C.; Dizene, R. Solar radiation forecasting using artificial neural network an random forest methods: Application to normal beam, horizontal diffuse and global components. Renew. Energy 2019, 132, 871–884. [Google Scholar] [CrossRef]
Figure 1. Architecture of an artificial neuron and a multi-layered neural network.
Figure 1. Architecture of an artificial neuron and a multi-layered neural network.
Applsci 09 00209 g001
Figure 2. ANN structure for GHI estimation.
Figure 2. ANN structure for GHI estimation.
Applsci 09 00209 g002
Figure 3. Average, minimum, and maximum values of the nRMSE and its standard deviation versus the number of inputs (Bouzareah site).
Figure 3. Average, minimum, and maximum values of the nRMSE and its standard deviation versus the number of inputs (Bouzareah site).
Applsci 09 00209 g003
Figure 4. Average, minimum, and maximum values of the nRMSE and its standard deviation versus the number of inputs without S as input (Bouzareah site).
Figure 4. Average, minimum, and maximum values of the nRMSE and its standard deviation versus the number of inputs without S as input (Bouzareah site).
Applsci 09 00209 g004
Figure 5. Estimated 5-min GHI versus measured 5-min GHI for various ANN architectures (Bouzareah site); (a) using 10 inputs; (b) using 6 inputs (with S); (c) using 9 inputs (without S); (d) using 5 inputs (without S).
Figure 5. Estimated 5-min GHI versus measured 5-min GHI for various ANN architectures (Bouzareah site); (a) using 10 inputs; (b) using 6 inputs (with S); (c) using 9 inputs (without S); (d) using 5 inputs (without S).
Applsci 09 00209 g005aApplsci 09 00209 g005b
Figure 6. ANN architectures for the estimation of ETI.
Figure 6. ANN architectures for the estimation of ETI.
Applsci 09 00209 g006
Figure 7. nRMSE evolution vs the number hidden neurons (Bouzareah site).
Figure 7. nRMSE evolution vs the number hidden neurons (Bouzareah site).
Applsci 09 00209 g007
Figure 8. Validation of the model for seven randomly chosen days (Bouzareah site).
Figure 8. Validation of the model for seven randomly chosen days (Bouzareah site).
Applsci 09 00209 g008
Figure 9. Experimental and modelled solar irradiance curves in clear sky conditions (hour in true solar time) (Odeillo site), (a) GHI, April; (b) DNI, April; (c) GHI, September; (d) DNI, September.
Figure 9. Experimental and modelled solar irradiance curves in clear sky conditions (hour in true solar time) (Odeillo site), (a) GHI, April; (b) DNI, April; (c) GHI, September; (d) DNI, September.
Applsci 09 00209 g009
Figure 10. Comparison of forecasting models for various horizons for hourly GHI (Odeillo site); (a) in term of nRMSE; (b) in term of RMSE.
Figure 10. Comparison of forecasting models for various horizons for hourly GHI (Odeillo site); (a) in term of nRMSE; (b) in term of RMSE.
Applsci 09 00209 g010
Figure 11. Comparison of forecasting models for various horizons for hourly DNI (Odeillo site); (a) in terms of nRMSE; (b) in term of RMSE.
Figure 11. Comparison of forecasting models for various horizons for hourly DNI (Odeillo site); (a) in terms of nRMSE; (b) in term of RMSE.
Applsci 09 00209 g011
Figure 12. Comparison of forecasting models for GHI and DNI in terms of nRMSE (Odeillo site).
Figure 12. Comparison of forecasting models for GHI and DNI in terms of nRMSE (Odeillo site).
Applsci 09 00209 g012
Table 1. Available meteorological data in Bouzareah.
Table 1. Available meteorological data in Bouzareah.
DataSymbolUnity
Measured data
Global horizontal solar irradiationGHIWh.m−2
Global tilted solar irradiation (36.8°)GTIWh.m−2
Ambient temperatureT°C
Relative humidityH%
Wind speedWSm.s−1
Wind directionWDdegree
PrecipitationPrmm
Sunshine durationSminutes
Atmospheric pressurePmbar
Calculated data
Extraterrestrial Horizontal irradiationEHIWh.m−2
Solar declinationδdegree
Zenith angle θzdegree
Table 2. Pearson correlation coefficients between meteorological variables for the Bouzareah site.
Table 2. Pearson correlation coefficients between meteorological variables for the Bouzareah site.
ParameterGHIδθZTHpPrWSWDSEHI
GHI1.000−0.052−0.1270.063−0.1800.095−0.063−0.0200.0160.6870.127
δ 1.000−0.3940.4310.170−0.243−0.012−0.003−0.1630.051−0.340
θz 1.000−0.060−0.048−0.068−0.012−0.017−0.0310.097−0.955
T 1.000−0.387−0.073−0.049−0.124−0.0920.054−0.065
H 1.0000.0890.0300.039−0.213−0.059−0.006
P 1.000−0.070−0.162−0.2290.112−0.090
Pr 1.0000.0130.0130.0150.006
WS 1.0000.1090.011−0.009
WD 1.0000.017−0.044
S 1.000−0.068
ETI 1.000
Table 3. Best configurations according to the number of inputs for the Bouzareah site (1 = present data, 0 = absent data).
Table 3. Best configurations according to the number of inputs for the Bouzareah site (1 = present data, 0 = absent data).
Nbr inputsRankδ θ z THpPrWSWDSETIMAE Wh/m2nMAE %RMSE Wh/m2nRMSE %MBE Wh/m2
10211111111115.08315.506.11618.65−0.084
9111111111105.08315.506.11518.650.011
9311111011115.11815.616.15218.76−0.034
8601111011115.19915.856.23419.01−0.018
8710111110115.20915.886.24519.04−0.012
71501111010115.28516.126.32319.28−0.038
72411111010105.35816.346.39919.51−0.032
61800111010115.30616.186.34519.35−0.081
63200111010115.41516.516.45719.69−0.021
510800011010115.72317.456.77520.660.018
511700110010115.75617.556.80920.76−0.015
419600110000116.06618.507.12821.74−0.007
422600011000116.21318.957.28022.20−0.018
336000010000117.16221.848.25725.18−0.003
336301000000117.24522.098.34325.44−0.005
237200000000117.38922.538.49125.890.003
237901000000108.08624.669.20928.080.062
1753010000000013.76241.9715.05845.920.062
1754000000000113.78042.0215.07645.97−0.004
Table 4. Best configurations according to the number of inputs (without S) (Bouzareah site). (1 = present data, 0 = absent data).
Table 4. Best configurations according to the number of inputs (without S) (Bouzareah site). (1 = present data, 0 = absent data).
Nbr inputsRankδ θ z THpPrWSWDEHIMAE Wh/m2nMAE %RMSE Wh/m2nRMSE %MBE Wh/m2
93821111111119.35428.5210.5232.07−1.142
83830111111119.42728.7510.5932.29−1.008
83841011111119.43528.7710.6032.32−1.058
73891011111019.53629.0810.7032.64−1.059
73910111110119.56629.1710.7332.73−1.010
63980111110019.66129.4610.8333.03−1.038
63991011110019.66129.4610.8333.031.002
54210111100019.80029.8810.9733.47−1.872
54260011100119.83830.0011.0133.58−1.013
446801111000010.07630.7311.2634.34−1.042
449600111000110.22631.1811.4134.80−0.929
358700011000111.93136.3813.1740.16−0.816
360501011000012.12036.9613.3740.750.918
263500010000112.41537.8613.6741.681.033
265001010000012.56638.3213.8342.161.063
175301000000013.76241.9715.0645.921.062
175400000000113.78042.0215.0745.97−0.904
Table 5. Average statistical parameters between measured and estimated global solar 5 min-irradiation on a 36.8° tilted plane for the station of Bouzareah. The bold line is the results for the best architecture.
Table 5. Average statistical parameters between measured and estimated global solar 5 min-irradiation on a 36.8° tilted plane for the station of Bouzareah. The bold line is the results for the best architecture.
Hidden Neurons NumberMAEnMAERMSEnRMSEMBE
Wh.m−2%Wh.m−2%Wh.m−2
15.9416.326.8918.93−0.97
22.958.123.8110.48−0.29
32.807.713.6510.04−0.50
42.416.633.218.810.31
52.567.043.479.54−0.80
62.757.563.7710.36−1.14
72.687.383.6810.12−1.09
82.707.433.6910.14−1.12
Table 6. Performance metrics (in Wh/m2 for RMSE, MAE, and MBE) for GHI (in bold the best predictor for each horizon and each error metric) (Odeillo site).
Table 6. Performance metrics (in Wh/m2 for RMSE, MAE, and MBE) for GHI (in bold the best predictor for each horizon and each error metric) (Odeillo site).
MetricModelh+1h+2h+3h+4h+5h+6
RMSESP97.72132.41157.08176.48193.06202.66
ANN101.79126.65141.95150.28154.82157.27
nRMSESP21.67%29.36%34.82%39.12%42.79%44.91%
ANN22.57%28.08%31.47%33.31%34.31%34.85%
MAESP56.9780.8298.70112.83124.64130.85
ANN72.8790.97106.83112.60117.56118.59
nMAESP12.63%17.92%21.88%25.01%27.62%28.99%
ANN16.16%20.17%23.68%24.96%26.05%26.28%
MBESP4.623.790.33−4.639.42−12.9
ANN−2.39−5.677.1210.1412.3214.7
Table 7. Performance metrics (in Wh/m2 for RMSE, MAE, and MBE) for DNI (in bold is the best predictor for each horizon and each error metric) (Odeillo site).
Table 7. Performance metrics (in Wh/m2 for RMSE, MAE, and MBE) for DNI (in bold is the best predictor for each horizon and each error metric) (Odeillo site).
MetricModelh+1h+2h+3h+4h+5h+6
RMSESP207.86287.64338.86378.19412.67434.51
ANN212.33270.10297.73321.59336.98344.00
nRMSESP37.42%51.77%60.98%68.04%74.23%78.16%
ANN38.23%48.62%53.58%57.85%60.61%61.88%
MAESP125.24187.12230.34266.01298.23317.71
ANN168.14223.27244.91274.60283.68299.16
nMAESP22.55%33.68%41.45%47.85%53.64%57.15%
ANN30.27%40.19%44.07%49.40%51.02%53.82%
MBESP2.282.101.280.63−0.58−2.66
ANN3.342.76−4.12−6.15−6.24−6.88

Share and Cite

MDPI and ACS Style

Notton, G.; Voyant, C.; Fouilloy, A.; Duchaud, J.L.; Nivet, M.L. Some Applications of ANN to Solar Radiation Estimation and Forecasting for Energy Applications. Appl. Sci. 2019, 9, 209. https://doi.org/10.3390/app9010209

AMA Style

Notton G, Voyant C, Fouilloy A, Duchaud JL, Nivet ML. Some Applications of ANN to Solar Radiation Estimation and Forecasting for Energy Applications. Applied Sciences. 2019; 9(1):209. https://doi.org/10.3390/app9010209

Chicago/Turabian Style

Notton, Gilles, Cyril Voyant, Alexis Fouilloy, Jean Laurent Duchaud, and Marie Laure Nivet. 2019. "Some Applications of ANN to Solar Radiation Estimation and Forecasting for Energy Applications" Applied Sciences 9, no. 1: 209. https://doi.org/10.3390/app9010209

APA Style

Notton, G., Voyant, C., Fouilloy, A., Duchaud, J. L., & Nivet, M. L. (2019). Some Applications of ANN to Solar Radiation Estimation and Forecasting for Energy Applications. Applied Sciences, 9(1), 209. https://doi.org/10.3390/app9010209

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop