Next Article in Journal
A Novel Co-Phase Power Supply System for Electrified Railway Based on V Type Connection Traction Transformer
Next Article in Special Issue
Residential Photovoltaic Profitability with Storage under the New Spanish Regulation: A Multi-Scenario Analysis
Previous Article in Journal
Effect of the Reactant Transportation on Performance of a Planar Solid Oxide Fuel Cell
Previous Article in Special Issue
A Genetic Algorithm Approach as a Self-Learning and Optimization Tool for PV Power Simulation and Digital Twinning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neural Network Approach for Global Solar Irradiance Prediction at Extremely Short-Time-Intervals Using Particle Swarm Optimization Algorithm

by
Ahmed Aljanad
1,2,*,
Nadia M. L. Tan
2,*,
Vassilios G. Agelidis
2,3 and
Hussain Shareef
4
1
Department of Science, Technology, Engineering and Mathematics, American University of Afghanistan, Darul Aman 4001, Kabul, Afghanistan
2
Institute of Power Engineering, Universiti Tenaga Nasional, Kajang 43000, Selangor, Malaysia
3
Department of Electrical Engineering, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
4
Department of Electrical Engineering, United Arab Emirates University, Al Ain 15551, Abu Dhabi, United Arab Emirates
*
Authors to whom correspondence should be addressed.
Submission received: 24 December 2020 / Revised: 21 January 2021 / Accepted: 22 January 2021 / Published: 23 February 2021

Abstract

:
Hourly global solar irradiance (GSR) data are required for sizing, planning, and modeling of solar photovoltaic farms. However, operating and controlling such farms exposed to varying environmental conditions, such as fast passing clouds, necessitates GSR data to be available for very short time intervals. Classical backpropagation neural networks do not perform satisfactorily when predicting parameters within short intervals. This paper proposes a hybrid backpropagation neural networks based on particle swarm optimization. The particle swarm algorithm is used as an optimization algorithm within the backpropagation neural networks to optimize the number of hidden layers and neurons used and its learning rate. The proposed model can be used as a reliable model in predicting changes in the solar irradiance during short time interval in tropical regions such as Malaysia and other regions. Actual global solar irradiance data of 5-s and 1-min intervals, recorded by weather stations, are applied to train and test the proposed algorithm. Moreover, to ensure the adaptability and robustness of the proposed technique, two different cases are evaluated using 1-day and 3-days profiles, for two different time intervals of 1-min and 5-s each. A set of statistical error indices have been introduced to evaluate the performance of the proposed algorithm. From the results obtained, the 3-days profile’s performance evaluation of the BPNN-PSO are 1.7078 of RMSE, 0.7537 of MAE, 0.0292 of MSE, and 31.4348 of MAPE (%), at 5-s time interval, where the obtained results of 1-min interval are 0.6566 of RMSE, 0.2754 of MAE, 0.0043 of MSE, and 1.4732 of MAPE (%). The results revealed that proposed model outperformed the standalone backpropagation neural networks method in predicting global solar irradiance values for extremely short-time intervals. In addition to that, the proposed model exhibited high level of predictability compared to other existing models.

1. Introduction

Sustainable energy sources such as photovoltaics (PV) generation is becoming increasingly important in the current times due to the depletion of natural resources, increasing energy demand, high cost of new fossil-fuel generation and transmission power system infrastructure, and worsening greenhouse gas effects. However, PV output power is not dispatchable in terms of supply and demand because of inherent intermittency in solar irradiance. Be it large-scale or nanogrid PV generation, energy storage devices such as batteries and ultracapacitors are required to manage energy and transient power demand, respectively [1]. Therefore, prediction of solar irradiance is essential in order to correctly size a solar PV power plant and energy storage system and to optimize the corresponding energy management system.
There are many predictive and artificial intelligence (AI) techniques that are proposed to predict solar irradiance and eliminate the limitations of measuring instrument, which feature high efficiency, flexibility, consistency, reliability, and so on. To obtain GSR data, different models are proposed by researchers, such as empirical [2], physical [3], and machine learning [4] models. Among them, machine learning models, presenting high accuracy and wider applicability, are widely studied worldwide [5]. Most of the prediction of GSR using the two methods have been for daily, hourly, or half-hourly time intervals. There has been no extensive studies on short-time, few minutes or few seconds, interval GSR prediction, which is advantageous for real-time compensation by energy management system (EMS) of PV generation and with battery or ultracapacitor.
In classical modelling methods, GSR are structured as standalone models. Artificial neural network is one of the more common AI technique that is used for prediction of GSR [6,7,8]. They include a variety of models such as multilayer perception neural networks (MLP-NNs) [9], radial basis function neural networks (RBF-NNs) [10], generalized regression neural networks (GR-NNs) [11], extreme learning machines [12], deep learning neural networks [13], convolutional neural networks [14], and so on. For example, authors of [15] compared MLP-NNs, RBF-NNs, and GR-NNs for daily global solar radiation estimation in three climate zones of China, illustrating that MLP-NNs and RBF-NNs performed much better than GR-NNs. Authors of [16] used latitude, longitude, and altitude geographical parameters and meteorological parameters like relative humidity, air pressure, clearness index (Kt), and average temperature (Ta) as inputs of ANN models to predict the GSR of stations in Zimbabwe. The results showed that Kt is the most effective feature to predict the GSR. Authors of [17] evaluated the performance of GSR prediction over other ANN techniques and architectures such as multilayer perceptron neural network, a neural network using wavelet transform, Elman Neural Network, and radial basis neural network. The performance of GSR prediction using other AI techniques such as support vector machine (SVM) [18], adaptive neuro-fuzzy inference system (ANFIS) [19], kernel-based nonlinear extension of Arps decline model and gradient boosting with categorical features support [20], and support vector regression [21] have also been investigated. Empirical methods are also widely discussed and evaluated as a predictive model of GSR in many countries [22]. Authors of [23], developed an empirical model for GSR in southwest Turkey where two models, cubic and quadratic are found to be the appropriate models. The drawbacks and limitations of these standalone and classical models are the insufficient capacity to capture and ascertain the deterministic features of meteorological data. The utilization of the temperature as a crucial input key features in ANN model is suggested by some researchers to have a high predictive performance [24].
Hybrid models are implemented as data-driven models to enhance and increase the capability of GSR prediction. Ibrahim and Tamer proposed a hybrid RF technique-based firefly algorithm (FFA) for predicting hourly GSR [25]. FFA is used to optimize the number of trees and leaves of a random forest (RF) technique. The RF technique is compared with different standalone and hybrid prediction models like ANN and ANN based FFA, respectively. The Mycielski-Makrov hybrid model is proposed to predict solar irradiance by using historical solar irradiance data without needing any other parameters [26]. The results of the model achieved better accuracy when it is compared with standalone Mycielski model. In [27], the solar irradiance performance and efficiency of hybrid model using auto-regressive average moving (ARMA) and feed-forward time-delay neural network showed improvement in prediction accuracy compared with classical ARMA model. In a similar manner, a number of independent hybrid day-of-the-year based models such as ANFIS with chaotic firefly algorithm and ANFIS with whale optimization algorithm with simulated annealing and roulette wheel selection have been employed in estimating GSR in China, where the variability of climate is wide-ranging [28]. The hybrid methods would increase the estimation accuracy of GSR despite the wide-ranging climate variability. Improvement in GSR prediction is also performed by hybridizing extreme-learning-machine and neural-grouping genetic algorithms into a single model, where optimal feature selection is achieved [29].
This study aims to overcome the drawback of standalone AI models by examining the adaptability, efficiency and accuracy of the hybrid backpropagation neural network and particle swarm optimization (BPNN-PSO) method for prediction of solar irradiance in very short-time intervals and fast changing climate conditions as there has been no extensive studies on prediction of GSR in those conditions in a tropical country. This paper also designs and develops a high-efficiency prediction model through the use of the actual meteorological data from Kajang, Malaysia. The meteorological data include pressure, humidity, temperature, wind speed, wind direction, and diffuse and direct irradiances.
The proposed PSO is used inside the BPNN to improve the prediction performance by reducing the error with actual data and improving the convergence rate by minimizing the objective function. The PSO optimizes hidden layers, neurons, and learning rate of BPNN architecture. The best minimized combination value of these three parameters in BPNN is the main objective function. The accuracy of the developed hybrid BPNN-PSO model is compared and evaluated with other existing model using reliable statistical indicators, such as root mean square error (RMSE), mean squared error (MSE), mean absolute error (MAE), and mean absolute percentage error (MAPE).
This paper is organized as follows: first, a literature review is presented. In Section 2, the BPNN algorithm is presented. In Section 3, data preparation for training and testing the PSO algorithm working concepts, and the performance evaluation of the model is discussed. In Section 4, the simulations results are analyzed. Finally, the achievements of the proposed method are highlighted in the conclusion.

2. Theory of Backpropagation Neural Network Structure

The ANN models nonlinear systems and make every effort to create the function of the human brain in a simulation environment through distributed processing. The ANN consists of neurons, which are simple connected elements that have the advantage in solving complex nonlinear relationship between system input and output [30,31].
Figure 1 shows the structure of the BPNN model employed. The model consists of seven meteorological inputs that are obtained from Kajang, Malaysia, which include pressure (P), temperature (T), humidity (H), wind speed (WS), wind direction (WD), irradiance direct (IDR), and irradiance diffuse (IDF), a few hidden layers and one GSR as output parameter. In order to predict GSR, this paper employs the Levenberg-Marquardt backpropagation algorithm to train the multi-layer perceptron of the ANN model in Matlab. The Levenberg-Marquardt algorithm has been selected due to its minimal localization error as well as its efficiency and speed [32]. The training of ANN using Levenberg-Marquardt back propagation algorithm involves three phases: (a) The feedforward phase, (b) the computation and backpropagation of the associated error, and (c) the adjustment of the weights. All inputs and outputs of the BPNN are expressed by the following equations:
I n p u t = T 1 T 2 T n H 1 H 2 T n P 1 P 2 P n W S 1 W S 2 W S n W D 1 W D 2 W D n I D R 1 I D R 2 I D R n   I D F 1 I D F 2 I D F n
O u t p u t = G S R 1 G S R 2 G S R n  
where GSR is the global solar irradiance of the meteorological datasets.
During the feed-forward phase of ANN training, data samples consisting of seven input values (P, T, H, WS, WD, IDR, and IDF) and one output or target value G S R are presented into the ANN. The explanation of BPNN working progress is summarized in the following steps:
Step 1: The weight and bias are randomly initialized.
Step 2: The input layer relays the input signals to the hidden nodes. Then, the variables in the hidden nodes are calculated by using [33],
n e t = Z i n _ j = P w 1 j + T w 1 j + H w 1 j + W S w 1 j + W D w 1 j + I D R w 1 j + I D F w 1 j
where Z i n the input of the hidden nodes, and j is the number of hidden nodes that is calculated by using the PSO and w is the weight factor.
Step 3: The hidden layer is calculated by using the sigmoid function and it is given by,
Z n e t = 2 1 + e n e t 1
For input pattern p, the i-th input layer node holds xp,i. Net input to j-th node in the hidden layer is,
n e t j =   i = 0 n w j , i x p , i + θ j , i
where, wj,i is the weight from the input layer to the hidden layer, θj,i represents the bias from the input layer to the hidden layer. Output of j-th node in the hidden layer is,
x p , j =   Z j i = 0 n w j , i x p , i + θ j , i
Net input to k-th node in the output layer is,
n e t k =   j w k , j x p , k + θ k , j
where, wk,i, θk,j, are the weight and bias from the hidden layer to the output layer, respectively. Output of k-th node in the output layer is,
O p , k =   Z k j w k , j x p , k + θ k , j
Step 4: In the BPNN training, after executing the feed-forward phase, the next phases are to compute the backpropagation of the associated error and then adjust the weights. During training, the GSR is compared with its target value in the sample data to determine the associated error (e). Based on this error, factor δ k that is used to distribute the error at the output layer back to all hidden nodes is given by,
δ k = Z k 1 Z k T k Z k
where T k is the true output (GSR). The error in the hidden layer is calculated as
δ j = Z j 1 Z j δ k w k , j
Step 5: BPNN in this phase, updates error and biases. Weights are updated using the following equations
Δ w k , j = α δ k S j
w k , j = w k , i + Δ w k , j
Δ w j , i = α δ k x p
w j , i = w j , i +   Δ w j , i
where α is the learning rate which can be assigned values between 0 and 1.
Biases are updated using the following equations
Δ θ k , j = α δ k
θ k , j = θ k , j + Δ θ k , j
Δ θ j , i = α δ j
θ j , i = θ j , i +   Δ θ j , i
The adjustment of the weights from the input to the hidden layer is based on factor δ j in (10) and the activation of the input features. Using (4) to (18), the BPNN process is repeated for all the data set samples to achieve one epoch. The training process will continue until the error goal or the predefined epoch is achieved. After the training process, the ANN can be utilized to generate the reference GSR with new input data.

3. Methodology Proposed Hybrid BPNN-PSO for Predicting GSR

The BPNN architecture is implemented at the process of training the GSR due to its capability to minimize output error by optimizing the input weight values of the output layer. However, three significant parameters in BPNN architecture namely hidden layers (HL), number of neurons (Ne), and learning rate (LR) are normally set based on trial and error method, which is time consuming. In order to further enhance the performance of the GSR prediction and overcome the time-consuming process of trial and error method, the PSO algorithm is utilized to find the optimal best combination of the dimension array. In addition, PSO is used to give reliable GSR prediction with less variance, less error, and the best fitting for the prediction function. Figure 2 shows the schematic diagram for prediction GSR using BPNN based PSO. The process is categorized into three phases which are clearly illustrated as follows:
Phase I: This phase begins with the collection of seven input variables at 5-s and 1-min time intervals. This study has included two profiles for the reason of improving the prediction performance by capturing the nonlinear association of patterns between different meteorological parameters, such as temperature, pressure, humidity, wind speed, wind direction, irradiance direct, and irradiance diffuse. After that, the sample data is moved through the normalization method.
Phase II and III: In these phases, 10-fold cross validation method is used for data pre-processing and division into training and testing observations. Then, the BPNN model is used for GSR prediction where the optimal number of hidden layer neurons and learning rate of BPNN model are computed based on PSO algorithm to enhance the accuracy of GSR.

3.1. Particle Swarm Optimization Algorithm

PSO is a heuristic optimization technique, which is based on population inspired by bird flocking and fish schooling for solving nonlinear problems with both discrete and continuous variables. The PSO algorithm is robust and easy to implement with global exploration capability in various applications [34]. In PSO, the potential solution to the problem being considered is randomly generated though a population individual particles also known as “swarm”. Each particle will move at arbitrary velocity across a dimensional search space to find two locations. The swarm will keep two locations. The first location is the best position in the current iteration also known as the local best, and the second location is the best point found in all previous iterations also known as the global best. The velocity and position factors are updated as [35]:
V i d t + 1 = w V i d t + c 1 r 1 P i d t X i d t + c 2 r 2 P t d t X i d t
  X i d t + 1 = X i d t + V i d t + 1
where   c 1 is the social rate, and   c 2 is the cognitive rate. r 1   and   r 2   denote the randomness in the interval (0,1), V is the velocity factor of agent i at iteration d , t is the present iteration, w is the inertia factor, and X is the position factor.

3.2. Performance Evaluation

To investigate the performance and validate the accuracy of each model, four statistical index errors are implemented. In addition, a previous benchmark study at different sites are compared with the proposed hybrid model. There are various statistical indices that have been used by other researchers [36]. However, MSE, MAE, and MAPE are the most common statistical error indexes [37]. Therefore, in this study, the performance evolution of the proposed techniques are investigated through the following indexes:
M A E =   i = 1 n 1 n   x   G S R A i G S R P i   W   m 2
R M S E = 1 n i = 1 n G S R A i G S R P i 2   W   m 2
M S E = 1 n i = 1 n G S R A i G S R P i 2   ( W   m 2 ) 2
M A P E = 1 n i = 1 n G S R A i G S R P i G S R P i %
where error = G S R A i G S R P i , G S R A i   is the state of charge of the actual data, and G S R P i is the state of charge of the predicted data, and n is the number of samples.

3.3. Data Preparation and Model Execution

The prediction of the GSR for one-day and three-day profiles begins with the collection of seven input variables at 5-s and 1-min time intervals. This study has included two profiles for the reason of improving the prediction performance by capturing the nonlinear association of patterns between different meteorological parameters, such as temperature, pressure, humidity, wind speed, wind direction, irradiance direct, and irradiance diffuse. At the first stage, the data have been collected from the TNBR—Solar Resource Monitoring Station, located in Kajang, Malaysia, with different time intervals. The data is collected from 1 March 2013 to 15 February 2014 using a high sampling data logger at a sampling rate of 5 and 30 s for 1, 5, 30, and 60 min interval. This measurement was taken as part of a Seeding Fund Project TNBR/SF140/2010 entitled “Development of Solar Research Facility for studies of Grid Connection of Utility Scale Solar Power Plant”. The Solar Resource Measurement was performed as follow: (a) GSR was measured by a CMP11 pyranometer, (b) Diffuse irradiance was measured by a CMP11 pyranometer and shaded by a shading ball attached, and (c) Direct irradiance was measured by a CHP1 pyrheliometer. The irradiance and temperature were measured using a solar pyranometer sensor and a temperature sensor, respectively. All measurements were performed instantaneously for every 5-s and 1-min time intervals and the measurements were recorded using the DT80 data logger.
The available meteorological data features would guarantee a high efficiency prediction of the solar irradiance during the short time intervals (5-s and 1-min) without any new features. Moreover, to increase and improve the accuracy of the output GSR parameter, all other possible data features that may have a direct effect on the output parameter performance are included in the predictive models. All the saved input datasets are normalized before being trained by BPNN model to increase the robustness and efficiency of the system. Besides, it enhances the convergence rate of the predicted GSR. In this study, the data normalization range is from 0 to 1 and is determined as:
N o r m = D a t a D a t a min D a t a max D a t a min × 0.9 + 0.1
where the Datamax and Datamin are the maximum and minimum value of the BPNN trained dataset. The testing dataset is also normalized with the same range limit. The trained dataset has been divided to 70% for training and 30% for testing, where 10-fold cross validation method is implemented in all input arrays.
Figure 3 shows the flowchart of the proposed hybrid BPNN-PSO implementation for predicting the PV GSR in detail. The implementation procedure of the proposed hybrid technique in determining the optimum number of hidden neurons (Ne), hidden layers (HL), and learning rate (LR).
Table 1 presents the parameters used in the initialization of the algorithm. Then the initial local and global best positions, Pbest and Gbest are randomly generated. In order to train the ANN and evaluate the fitness value, the number of learning rate, hidden neurons, and layers should be selected. The objective function MAE is evaluated based on (21) and the new Pbest and Gbest are updated. The velocity and position are computed and updated using (19) and (20). If the position is rejected, other combinations have to be identified. Training of ANN and evaluation of MAE should be repeated until the maximum population size is reached or when Pbest is lower than Gbest. If the local best value is less than global best value, the previous is updated as the global best value. The process is iterated until the stop criterion is reached so as to obtain the optimal hidden layers and hidden neurons.

4. Discussion and Results

4.1. Objective Function Performance of PV Solar Irradiance

Figure 4 and Figure 5 depict solar irradiance convergence results using BPNN with the assistance of PSO algorithm, which are evaluated under different time intervals of 5-s and 1-min for both 1-day and 3-days profiles. The performance evaluation is represented by statistical index error or objective function (MAE), which compares the simulation convergence results over 100 iterations. The optimization parameters of population size and the iteration numbers are standardized with 4 and 100 maximum iterations, respectively. In Figure 4, the proposed hybrid algorithm BPNN-PSO shows less error compared to actual GSR data, where the minimum objective function MAE values obtained from 3-Days profile for 5-s and 1-min time intervals are 0.7537, and 0.2754, respectively. Figure 5 presents the one-day profile with better convergence rate over the 3-days profile for both time intervals of 5-s and 1-min with 0.1000 and 0.0956, respectively.
Other statistical index errors like MSE, RMSE, and MAPE are used to give more extensive analysis of the output performance for both profiles in 1 day and 3 days as shown in Table 2. Under the nonlinearity behaviors of GSR, the proposed predictive hybrid technique is investigated with the aforementioned analytical index error tools, where all other statistical index errors like RMSE, MSE, and MAPE, show the superiority of the proposed BPNN-PSO technique in 3-days profile with 1.7078, 0.0292, and 31.4348 for 5-s, and 0.6566, 0.0043, and 1.4732 for 1-min time intervals, respectively. Likewise, the 1-day profile shows the superior performance of the proposed hybrid model with 0.1911 (RMSE), 0.0004 (MSE), and 0.7484 (MAPE) for 5-s, and 0.2032, 0.0004, and 1.1271 for 1-min time intervals.

4.2. PV Solar Irradiance Optimal Parameters

The optimal parameters HL, Ne, and LR, acquired after implementing the heuristic optimization algorithms PSO for the 1-day and 3-days profiles are shown in Table 3. In the 3-days profile, the BPNN-PSO algorithm attains hidden layers of 1 (7 Ne) and 2 (14 and 9 Ne) after 22 and 35 iterations for 5-s and 1-min time intervals. In the 1-day profile, the proposed BPNN-PSO achieved the optimal value of hidden layers of 2 (8 and 2 Ne) and (9 and 11 Ne) after 76 and 91 iterations for 5-s and 1-min time intervals, respectively. In contrast, the best learning rate values of 3-days and 1-day profiles are 0.1295, 0.7373, 0.5946, and 0.6481 during both time intervals of 5-s and 1-min, respectively.

4.3. PV Solar Irradiance Prediction

The input parameters considered during the training process of BPNN and PSO are the seven input as shown in Figure 2. The selection of the input features is set to seven after conducting a set of trial-and-error tests to the predictive models, the increase of the input features from one to seven increases the performance of the proposed predictive model BPNN-PSO. Any further increase in the number of input features after seven did not show further increase in the performance of the BPNN and PSO.
Figure 6 and Figure 7 present the predicted solar irradiance using the proposed hybrid BPNN-PSO algorithm and compare the predicted GSR with the reference or actual data. In the figures, the red line represents the actual global solar irradiance data obtained from Kajang, Malaysia, while the blue line is the solar irradiance prediction of the proposed hybrid BPNN-PSO algorithm. It is clearly observed that the performance of the BPNN-PSO over-classed other techniques in both 3-days and 1-day profiles, In Figure 6a and Figure 7a, the predicted solar irradiance of BPNN-PSO is almost aligned with the actual data in both profiles. The time domain response agrees well with the different statistical index errors in Table 2, which proves the superior performance of the proposed technique over the other conventional test techniques under the two time intervals of 5-s and 1-min, respectively. From Figure 7a, it is also noticeable that the results obtained from the proposed BPNN-PSO technique is robust and able to track the fast variation of the actual environmental data. The absolute error between predicted and actual values of 3-Days profile at 5-s and 1-min intervals is represented with the enlarged visual box to provide more clarity on the GSR prediction results, as depicted in Figure 6b and Figure 7b.
Moreover, the 1-day profile is also tested to investigate the ability and adaptability of predictive models under low data set and fast nonlinear environmental change of the solar irradiance. The 1-day profile are used to train the BPNN-PSO with the meteorological data obtained from Kajang, Malaysia, which were recorded on 22 February 2014. The predicted GSR of the proposed BPNN-PSO technique is compared at 5-s and 1-min time intervals. The BPNN-PSO model has a good alignment with the actual data, as depicted in Figure 8a and Figure 9a. The absolute error between predicted values of GSR using the proposed model with the target test values is presented in both 5-s and 1-min time intervals, respectively, as shown in Figure 8b and Figure 9b. The maximum error is approximately 8%, which is negligible.

4.4. Performance Comparison Using Regression Coefficient

The regression coefficient (R) is used as an indicator of the predictive model’s training process performance. The predicted data is displayed in black circle, while the blue line represents the reference value (Actual target). The regression coefficient results are very close to unity, which validates the accuracy of the model. The regression values of 3-Days profile with 5-s and 1-min time interval are 0.99951 and 0.99993, as shown in Figure 10. Moreover, Figure 11 shows that the 1-day profile have regression coefficients of 0.99999 at 5-s and 1-min intervals, respectively.

5. Model Validation with the Existing Methods

Table 4 presents the GSR prediction results of the proposed method that are benchmarked with existing results. The proposed hybrid model is compared with various empirical techniques, artificial intelligence, and hybrid artificial intelligence. The prediction accuracy of the proposed and existing methods is investigated through the four statistical index errors RMSE, MSE, MAE, and MAPE, respectively. The results obtained from the proposed study show greater performance in terms of predictability and improved forecasting capability as compared with the rest of the other methods.

6. Conclusions

This paper has presented a hybrid prediction model using BPNN based PSO for the enhancement the GSR prediction performance. Two profiles, 3-days and 1-day have been investigated during training and validation process. The main contribution of this paper is developing a robust and consistent BPNN-PSO model for prediction of global solar irradiance at tropical country like Malaysia in extremely short-time intervals. Secondly, the implementation of PSO algorithm has significantly enhanced the classical BPNN architecture, by finding the optimal values of the architecture parameters, namely, hidden layers, neurons, and learning rate. The performance results of the proposed BPNN-PSO model have been compared with other widely used neural network models. Statistical error indicators RMSE, MAE, MSE, and MAPE have been used for performance and precision evaluation of all models. When the proposed model is used to predict solar irradiance based on the dataset of one region in Malaysia, the developed model has shown remarkable prediction improvements, proving that the model is superior to other techniques in terms of reliability, adaptability, and accurate correlation in GSR prediction of fast, short-time intervals, and nonlinear nature. From the results obtained, the 3-days profile performance evaluation of the BPNN-PSO are 1.7078 of RMSE, 0.7537 of MAE, 0.0292 of MSE, and 31.4348 of MAPE (%), at 5-s time interval, whereas the obtained results of 1-min interval are 0.6566 of RMSE, 0.2754 of MAE, 0.0043 of MSE, and 1.4732 of MAPE (%). In contrast, the 1-day profile performance evaluation of the proposed method are 0.1911 of RMSE, 0.1000 of MAE, 0.0004 of MSE, and 0.7484 of MAPE (%), at 5-s time interval, whereas the obtained results of 1-min interval are 0.2032 of RMSE, 0.0956 of MAE, 0.0004 of MSE, and 1.1271 of MAPE (%). Even with the high accuracy performance of the model in different time intervals, the high performance is restricted to the availability of the aforementioned meteorological parameters. Moreover, the execution of the model needs to be extended to include spatial database using the proposed model for GSR prediction in short-time intervals. In addition, the optimization and calibration of the model could be proposed for the future work to make the model adaptable in different world regions.

Author Contributions

Conceptualization, A.A. and V.G.A.; methodology, A.A.; simulation, A.A.; validation, A.A. and N.M.L.T.; resources, H.S.; writing—original draft preparation, A.A., N.M.L.T., and V.G.A.; writing—review and editing, A.A., N.M.L.T., and V.G.A.; supervision, N.M.L.T.; project administration, N.M.L.T.; funding acquisition, N.M.L.T. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to acknowledge the financial support received from Universiti Tenaga Nasional, Malaysia, BOLD Research Grant 2020 (BOLD 2020), Project Code: RJO10517844/110.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

GSRGlobal solar irradiance
BPNNBackpropagation neural network
PSOParticle swarm optimization
PVPhotovoltaics
ANNArtificial neural network
ANFISAdaptive neural-fuzzy inference system
MLFFNNMultilayer feedforward neural network
MLP-NNMultilayer perception neural network
MAEMean absolute error
GR-NNGeneralized regression neural network
RBF-NNRadial basis function neural network
RMSERoot mean square error

References

  1. Jung, S.; Yoon, Y.T. Optimal Operating Schedule for Energy Storage System: Focusing on Efficient Energy Management for Microgrid. Processes 2019, 7, 80. [Google Scholar] [CrossRef] [Green Version]
  2. Zhou, Y.; Wang, D.; Liu, Y.; Liu, J. Diffuse solar radiation models for different climate zones in China: Model evaluation and general model development. Energy Convers. Manag. 2019, 185, 518–536. [Google Scholar] [CrossRef]
  3. Verbois, H.; Huva, R.; Rusydi, A.; Walsh, W. Solar irradiance forecasting in the tropics using numerical weather prediction and statistical learning. Sol. Energy 2018, 162, 265–277. [Google Scholar] [CrossRef]
  4. Wang, Z.; Tian, C.; Zhu, Q.; Huang, M. Hourly Solar Radiation Forecasting Using a Volterra-Least Squares Support Vector Ma-chine Model Combined with Signal Decomposition. Energies 2018, 11, 68. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, Y.; Zhou, Y.; Chen, Y.; Wang, D.; Wang, Y.; Zhu, Y. Comparison of support vector machine and copula-based nonlinear quan-tile regression for estimating the daily diffuse solar radiation: A case study in China. Renew. Energy 2020, 146, 1101–1112. [Google Scholar] [CrossRef]
  6. El Mghouchi, Y.; Chham, E.; Zemmouri, E.; El Bouardi, A. Assessment of different combinations of meteorological parameters for predicting daily global solar radiation using artificial neural networks. Build. Environ. 2019, 149, 607–622. [Google Scholar] [CrossRef]
  7. Marzouq, M.; Bounoua, Z.; El Fadili, H.; Mechaqrane, A.; Zenkouar, K.; Lakhliai, Z. New daily global solar irradiation estimation model based on automatic selection of input parameters using evolutionary artificial neural networks. J. Clean. Prod. 2019, 209, 1105–1118. [Google Scholar] [CrossRef]
  8. Wang, F.; Xuan, Z.; Zhen, Z.; Li, K.; Wang, T.; Shi, M. A day-ahead PV power forecasting method based on LSTM-RNN model and time correlation modification under partial daily pattern prediction framework. Energy Convers. Manag. 2020, 212, 112766. [Google Scholar] [CrossRef]
  9. Jadidi, A.; Menezes, R.J.A.; De Souza, N.; Lima, A.D.C. A Hybrid GA–MLPNN Model for One-Hour-Ahead Forecasting of the Global Horizontal Irradiance in Elizabeth City, North Carolina. Energies 2018, 11, 2641. [Google Scholar] [CrossRef] [Green Version]
  10. Jiang, H. A novel approach for forecasting global horizontal irradiance based on sparse quadratic RBF neural network. Energy Convers. Manag. 2017, 152, 266–280. [Google Scholar] [CrossRef]
  11. Lotfinejad, M.M.; Hafezi, R.; Khanali, M.; Hosseini, S.S.; Mehrpooya, M.; Shamshirband, S. A Comparative Assessment of Predicting Daily Solar Radiation Using Bat Neural Network (BNN), Generalized Regression Neural Network (GRNN), and Neuro-Fuzzy (NF) System: A Case Study. Energies 2018, 11, 1188. [Google Scholar] [CrossRef] [Green Version]
  12. Ghimire, S.; Deo, R.C.; Downs, N.J.; Raj, N. Self-adaptive differential evolutionary extreme learning machines for long-term solar radiation prediction with remotely-sensed MODIS satellite and Reanalysis atmospheric products in solar-rich cities. Remote Sens. Environ. 2018, 212, 176–198. [Google Scholar] [CrossRef]
  13. Ghimire, D.; Raj, M. Deep Learning Neural Networks Trained with MODIS Satellite-Derived Predictors for Long-Term Global Solar Radiation Prediction. Energies 2019, 12, 2407. [Google Scholar] [CrossRef] [Green Version]
  14. Dong, N.; Chang, J.-F.; Wu, A.-G.; Gao, Z.-K. A novel convolutional neural network framework based solar irradiance prediction method. Int. J. Electr. Power Energy Syst. 2020, 114, 105411. [Google Scholar] [CrossRef]
  15. Wang, L.; Kisi, O.; Zounemat-Kermani, M.; Salazar, G.A.; Zhu, Z.; Gong, W. Solar radiation prediction using different techniques: Model evaluation and comparison. Renew. Sustain. Energy Rev. 2016, 61, 384–397. [Google Scholar] [CrossRef]
  16. Elsheikh, A.H.; Sharshir, S.W.; Elaziz, M.A.; Kabeel, A.E.; Guilan, W.; Zhang, H. Modeling of solar energy systems using artificial neural network: A comprehensive review. Sol. Energy 2019, 180, 622–639. [Google Scholar] [CrossRef]
  17. Meenal, R.; Selvakumar, A.I. Assessment of SVM, empirical and ANN based solar radiation prediction models with most influencing input parameters. Renew. Energy 2018, 121, 324–343. [Google Scholar] [CrossRef]
  18. Mohammadi, K.; Shamshirband, S.; Tong, C.W.; Alam, K.A.; Petković, D. Potential of adaptive neuro-fuzzy system for prediction of daily global solar radiation by day of the year. Energy Convers. Manag. 2015, 93, 406–413. [Google Scholar] [CrossRef]
  19. Wu, L.; Huang, G.; Fan, J.; Zhang, F.; Wang, X.; Zeng, W. Potential of kernel-based nonlinear extension of Arps decline model and gradient boosting with categorical features support for predicting daily global solar radiation in humid regions. Energy Convers. Manag. 2019, 183, 280–295. [Google Scholar] [CrossRef]
  20. Khosravi, A.; Koury, R.; Machado, L.; Pabon, J. Prediction of hourly solar radiation in Abu Musa Island using machine learning algorithms. J. Clean. Prod. 2018, 176, 63–75. [Google Scholar] [CrossRef]
  21. Fan, J.; Wu, L.; Zhang, F.; Cai, H.; Zeng, W.; Wang, X.; Zou, H. Empirical and machine learning models for predicting daily global solar radiation from sunshine duration: A review and case study in China. Renew. Sustain. Energy Rev. 2019, 100, 186–212. [Google Scholar] [CrossRef]
  22. Bayrakçı, H.C.; Demircan, C.; Keçebaş, A. The development of empirical models for estimating global solar radiation on horizontal surface: A case study. Renew. Sustain. Energy Rev. 2018, 81, 2771–2782. [Google Scholar] [CrossRef]
  23. Wang, H.; Liu, Y.; Zhou, B.; Li, C.; Cao, G.; Voropai, N.; Barakhtenko, E. Taxonomy research of artificial intelligence for deterministic solar power forecasting. Energy Convers. Manag. 2020, 214, 112909. [Google Scholar] [CrossRef]
  24. Ibrahim, I.A.; Khatib, T. A novel hybrid model for hourly global solar radiation prediction using random forests technique and firefly algorithm. Energy Convers. Manag. 2017, 138, 413–425. [Google Scholar] [CrossRef]
  25. Hocaoglu, F.O.; Serttaş, F. A novel hybrid (Mycielski-Markov) model for hourly solar radiation forecasting. Renew. Energy 2017, 108, 635–643. [Google Scholar] [CrossRef]
  26. Ji, W.; Chee, K.C. Prediction of hourly solar radiation using a novel hybrid model of ARMA and TDNN. Sol. Energy 2011, 85, 808–817. [Google Scholar] [CrossRef]
  27. Zang, H.; Cheng, L.; Ding, T.; Cheung, K.W.; Wang, M.; Wei, Z.; Sun, G. Estimation and validation of daily global solar radiation by day of the year-based models for different climates in China. Renew. Energy 2019, 135, 984–1003. [Google Scholar] [CrossRef]
  28. Aybar-Ruiz, A.; Jiménez-Fernández, S.; Cornejo-Bueno, L.; Casanova-Mateo, C.; Sanz-Justo, J.; Salvador-González, P.; Salcedo-Sanz, S. A novel Grouping Genetic Algorithm–Extreme Learning Machine approach for global solar radiation prediction from numerical weather models inputs. Sol. Energy 2016, 132, 129–142. [Google Scholar] [CrossRef]
  29. Qazi, A.; Fayaz, H.; Wadi, A.; Raj, R.G.; Rahim, N.; Khan, W.A. The artificial neural network for solar radiation prediction and designing solar systems: A systematic literature review. J. Clean. Prod. 2015, 104, 1–12. [Google Scholar] [CrossRef]
  30. Alsina, E.F.; Bortolini, M.; Gamberi, M.; Regattieri, A. Artificial neural network optimisation for monthly average daily global solar radiation prediction. Energy Convers. Manag. 2016, 120, 320–329. [Google Scholar] [CrossRef]
  31. Premalatha, N.; Arasu, A.V. Prediction of solar radiation for solar systems by using ANN models with different back propagation algorithms. J. Appl. Res. Technol. 2016, 14, 206–214. [Google Scholar] [CrossRef] [Green Version]
  32. Yadav, A.K.; Chandel, S. Solar radiation prediction using Artificial Neural Network techniques: A review. Renew. Sustain. Energy Rev. 2014, 33, 772–781. [Google Scholar] [CrossRef]
  33. Mohandes, M.A. Modeling global solar radiation using Particle Swarm Optimization (PSO). Sol. Energy 2012, 86, 3137–3145. [Google Scholar] [CrossRef]
  34. Ghazvinian, H.; Mousavi, S.-F.; Karami, H.; Farzin, S.; Ehteram, M.; Hossain, S.; Fai, C.M.; Bin Hashim, H.; Singh, V.P.; Ros, F.C.; et al. Integrated support vector regression and an improved particle swarm optimization-based model for solar radiation prediction. PLoS ONE 2019, 14, e0217634. [Google Scholar] [CrossRef]
  35. Ehteram, M.; Ahmed, A.N.; Chow, M.F.; Afan, H.A.; El-Shafie, A. Accuracy Enhancement for Zone Mapping of a Solar Radiation Forecasting Based Multi-Objective Model for Better Management of the Generation of Renewable Energy. Energies 2019, 12, 2730. [Google Scholar] [CrossRef] [Green Version]
  36. Malvoni, M.; Hatziargyriou, N. One-day ahead PV power forecasts using 3D Wavelet Decomposition. In Proceedings of the 2019 International Conference on Smart Energy Systems and Technologies (SEST), Porto, Portugal, 9–11 September 2019; pp. 1–6. [Google Scholar]
  37. Zou, L.; Wang, L.; Xia, L.; Zou, L.; Hu, B.; Zhu, H. Prediction and comparison of solar radiation using improved empirical models and Adaptive Neuro-Fuzzy Inference Systems. Renew. Energy 2017, 106, 343–353. [Google Scholar] [CrossRef]
  38. Ghimire, S.; Deo, R.C.; Downs, N.J.; Raj, N. Global solar radiation prediction by ANN integrated with European Centre for medium range weather forecast fields in solar rich cities of Queensland Australia. J. Clean. Prod. 2019, 216, 288–310. [Google Scholar] [CrossRef]
  39. Çelik, Ö.; Teke, A.; Yıldırım, H.B. The optimized artificial neural network model with Levenberg–Marquardt algorithm for global solar radiation estimation in Eastern Mediterranean Region of Turkey. J. Clean. Prod. 2016, 116, 1–12. [Google Scholar] [CrossRef]
  40. Shamshirband, S.; Mohammadi, K.; Chen, H.-L.; Samy, G.N.; Petković, D.; Ma, C. Daily global solar radiation prediction from air temperatures using kernel extreme learning machine: A case study for Iran. J. Atmospheric Solar-Terrestrial Phys. 2015, 134, 109–117. [Google Scholar] [CrossRef]
  41. Quej, V.H.; Almorox, J.; Ibrakhimov, M.; Saito, L. Empirical models for estimating daily global solar radiation in Yucatán Peninsula, Mexico. Energy Convers. Manag. 2016, 110, 448–456. [Google Scholar] [CrossRef]
  42. Anis, S.; Jamil, B.; Ansari, A.; Bellos, E. Generalized models for estimation of global solar radiation based on sunshine duration and detailed comparison with the existing: A case study for India. Sustain. Energy Technol. Assess. 2019, 31, 179–198. [Google Scholar] [CrossRef]
  43. Olatomiwa, L.; Mekhilef, S.; Shamshirband, S.; Mohammadi, K.; Petković, D.; Sudheer, C. A support vector machine–firefly algorithm-based model for global solar radiation prediction. Sol. Energy 2015, 115, 632–644. [Google Scholar] [CrossRef]
  44. Khosravi, A.; Nunes, R.; Assad, M.E.H.; Machado, L. Comparison of artificial intelligence methods in estimation of daily global solar radiation. J. Clean. Prod. 2018, 194, 342–358. [Google Scholar] [CrossRef]
  45. Halabi, L.M.; Mekhilef, S.; Hossain, M. Performance evaluation of hybrid adaptive neuro-fuzzy inference system models for predicting monthly global solar radiation. Appl. Energy 2018, 213, 247–261. [Google Scholar] [CrossRef]
Figure 1. Backpropagation neural networks (BPNN) structure with seven inputs, one output, and multi hidden nodes.
Figure 1. Backpropagation neural networks (BPNN) structure with seven inputs, one output, and multi hidden nodes.
Energies 14 01213 g001
Figure 2. The structure of the proposed backpropagation neural network based on Particle Swarm Optimization Algorithm (PSO) algorithm.
Figure 2. The structure of the proposed backpropagation neural network based on Particle Swarm Optimization Algorithm (PSO) algorithm.
Energies 14 01213 g002
Figure 3. Flowchart of the proposed hybrid backpropagation neural network based on particle swarm optimization algorithm (BPNN-PSO) model.
Figure 3. Flowchart of the proposed hybrid backpropagation neural network based on particle swarm optimization algorithm (BPNN-PSO) model.
Energies 14 01213 g003
Figure 4. The convergence performance curves of solar irradiance prediction for 3-days, with 1-min and 5-s time interval.
Figure 4. The convergence performance curves of solar irradiance prediction for 3-days, with 1-min and 5-s time interval.
Energies 14 01213 g004
Figure 5. The convergence performance curves of solar irradiance prediction for 1-Day, with 1-min and 5-s time interval.
Figure 5. The convergence performance curves of solar irradiance prediction for 1-Day, with 1-min and 5-s time interval.
Energies 14 01213 g005
Figure 6. The photovoltaics solar prediction results using 3-days profile with 5-s time interval. (a) Solar irradiance prediction. (b) Solar irradiance prediction error.
Figure 6. The photovoltaics solar prediction results using 3-days profile with 5-s time interval. (a) Solar irradiance prediction. (b) Solar irradiance prediction error.
Energies 14 01213 g006
Figure 7. The photovoltaics prediction results using 3-days profile with 1-min time interval. (a) Solar irradiance prediction. (b) Solar irradiance prediction error.
Figure 7. The photovoltaics prediction results using 3-days profile with 1-min time interval. (a) Solar irradiance prediction. (b) Solar irradiance prediction error.
Energies 14 01213 g007
Figure 8. The prediction results using 1-day profile with 5-s time interval. (a) Photovoltaics (PV) solar irradiance prediction. (b) Photovoltaics (PV) solar irradiance prediction error.
Figure 8. The prediction results using 1-day profile with 5-s time interval. (a) Photovoltaics (PV) solar irradiance prediction. (b) Photovoltaics (PV) solar irradiance prediction error.
Energies 14 01213 g008
Figure 9. The prediction results using 1-day profile with 1-min time interval. (a) Photovoltaics (PV) solar irradiance prediction. (b) Photovoltaics (PV) solar irradiance prediction error.
Figure 9. The prediction results using 1-day profile with 1-min time interval. (a) Photovoltaics (PV) solar irradiance prediction. (b) Photovoltaics (PV) solar irradiance prediction error.
Energies 14 01213 g009
Figure 10. Regression performance of 3-days profile (a) BPNN-PSO, with 5-s time interval. (b) BPNN-PSO, with 1-min time interval.
Figure 10. Regression performance of 3-days profile (a) BPNN-PSO, with 5-s time interval. (b) BPNN-PSO, with 1-min time interval.
Energies 14 01213 g010
Figure 11. Regression performance of 1-day profile. (a) BPNN-PSO with 5-s time interval. (b) BPNN-PSO with 1-min time interval.
Figure 11. Regression performance of 1-day profile. (a) BPNN-PSO with 5-s time interval. (b) BPNN-PSO with 1-min time interval.
Energies 14 01213 g011
Table 1. PSO algorithm initialization parameters.
Table 1. PSO algorithm initialization parameters.
ParametersSymbolValue
Population SizeN10
No. of DimensionsD5
No. of IterationsT500
Maximum WeightWmax0.9
Minimum WeightWmin0.4
Acceleration Coefficientc1, c22
Table 2. Performance comparison using statistical index errors of the 1-day and 3-days profiles.
Table 2. Performance comparison using statistical index errors of the 1-day and 3-days profiles.
Type of ProfileTime IntervalsRMSE × 10−2
(W m−2)
MAE × 10−2
(W m−2)
MSE × 10−2MAPE [%]
3-Days Profile5-s1.70780.75370.029231.4348
1 min0.65660.27540.00431.4732
1-Day Profile5-s0.19110.10000.00040.7484
1 min0.20320.09560.00041.1271
Table 3. Optimal parameters of 1 and 3-Days profile at 5-s and 1-min time intervals.
Table 3. Optimal parameters of 1 and 3-Days profile at 5-s and 1-min time intervals.
ProfileTime IntervalHidden LayersNeuron 1Neuron 2Neuron 3Learning Rate Number of Iterations
3-Days5-s17000.129522
1-min214900.737335
1-Day5-s28200.594676
1-min291100.648191
Table 4. Comparison of statistical index errors between proposed method and other existing methods.
Table 4. Comparison of statistical index errors between proposed method and other existing methods.
ModelTime IntervalStatistical Error IndexesStudy Location
RMSE (MJ/m2/day)MAE (MJ/m2/day)MSEMAPE (%)
ANN, Genetic Programming (GP) [38]1 h1.613,
2.142
1.146,
1.629
Australia
ANN [39]1 h3.288Turkey (Mersin)
SVR [40]1 h2.5243Iran
Empirical [41]1 h2.52216.078Mexico (Yucatan Peninsula/Calakmut)
RF-FFA [24]1 h18.97976.3826Malaysia
Generalized Models [42]1 h1.79251.3800India
ANFIS
[18]
1 h1.04824.6402Iran
SVM-FFA [43]1 h0.69886.1768Nigeria
MLFFNN [44]1 h0.3214 (kWh/m2/day)0.25310.1033 (kWh/m2/day)3.316Iran
ANFIS,
ANFIS-PSO,
ANFIS-GA,
ANFIS-DE [45]
1 h0.3712,
0.3121,
0.3285,
0.3765
Malaysia
BPNN-PSO (Proposed Model)5 s, 1 min0.1911, 0.20320.1000, 0.09560.0004, 0.00040.7484, 1.1271Malaysia
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aljanad, A.; Tan, N.M.L.; Agelidis, V.G.; Shareef, H. Neural Network Approach for Global Solar Irradiance Prediction at Extremely Short-Time-Intervals Using Particle Swarm Optimization Algorithm. Energies 2021, 14, 1213. https://0-doi-org.brum.beds.ac.uk/10.3390/en14041213

AMA Style

Aljanad A, Tan NML, Agelidis VG, Shareef H. Neural Network Approach for Global Solar Irradiance Prediction at Extremely Short-Time-Intervals Using Particle Swarm Optimization Algorithm. Energies. 2021; 14(4):1213. https://0-doi-org.brum.beds.ac.uk/10.3390/en14041213

Chicago/Turabian Style

Aljanad, Ahmed, Nadia M. L. Tan, Vassilios G. Agelidis, and Hussain Shareef. 2021. "Neural Network Approach for Global Solar Irradiance Prediction at Extremely Short-Time-Intervals Using Particle Swarm Optimization Algorithm" Energies 14, no. 4: 1213. https://0-doi-org.brum.beds.ac.uk/10.3390/en14041213

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