Next Article in Journal
Application of the Navigational Air-Sea Methane Exchange Flux Observation System in the Qiongdongnan Basin of the Northern South China Sea
Previous Article in Journal
MODELING of Rarefied Gas Flows Inside a Micro-Nozzle Based on the DSMC Method Coupled with a Modified Gas–Surface Interaction Model
Previous Article in Special Issue
A Study on the Vulnerability Cascade Propagation of Integrated Energy Systems in the Transportation Industry Based on the Petri Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting Electricity Consumption in the Kingdom of Saudi Arabia

1
Sadat Academy for Management Sciences, Cairo 2222, Egypt
2
Department of Economics & Management Sciences, NED University of Engineering & Technology, Karachi 75270, Pakistan
3
College of Business Administration, University of Business and Technology, Jeddah 21361, Saudi Arabia
4
Faculty of Management, University of Primorska, SI-6000 Koper-Capodistria, Slovenia
5
Faculty of High Asian Studies, Zagazig University, Zagazig 31527, Egypt
*
Authors to whom correspondence should be addressed.
Energies 2023, 16(1), 506; https://doi.org/10.3390/en16010506
Submission received: 27 September 2022 / Revised: 29 November 2022 / Accepted: 22 December 2022 / Published: 2 January 2023

Abstract

:
Forecasting energy consumption in Saudi Arabia for the period from 2020 until 2030 is investigated using a two-part composite model. The first part is the frontier, and the second part is the autoregressive integrated moving average (ARIMA) model that helps avoid the large disparity in predictions in previous studies, which is what this research seeks to achieve. The sample of the study has a size of 30 observations, which are the actual consumption values in the period from 1990 to 2019. The philosophy of this installation is to reuse the residuals to extract the remaining values. Therefore, it becomes white noise and the extracted values are added to increase prediction accuracy. The residuals were calculated and the ARIMA (0, 1, 0) model with a constant was developed both of the residual sum of squares and the root means square errors, which were compared in both cases. The results demonstrate that prediction accuracy using complex models is better than prediction accuracy using single polynomial models or randomly singular models by an increase in the accuracy of the estimated consumption and an improvement of 18.5% as a result of the synthesizing process, which estimates the value of electricity consumption in 2030 to be 575 TWh, compared to the results of previous studies, which were 365, 442, and 633 TWh.

1. Introduction

The Kingdom of Saudi Arabia (KSA) is currently undergoing significant development in all fields, especially in the production of electric energy and electric energy consumption. The government is formulating future plans for the advancement of this sector due to its importance in targeted sustainable development and the mandate to generate electricity through alternative energy sources, as contained in Vision 2030. This is represented in the establishment of a number of economic development projects that pave the way to the goals of the KSA’s Vision 2030. This includes goals such as raising the competitiveness of renewable energy and managing carbon emissions in KSA.
KSA is the 14th largest consumer of electricity in the world. Its consumption is similar to that of the most densely populated countries, for example, Mexico, whose population in 2019 reached 127.5 million, compared to 34.2 million for KSA. It is also on par with the most advanced economies, for example, Italy, whose gross domestic product (GDP) for the year 2019 amounted to USD 2151.4 billion, compared to USD 704.0 billion for KSA. The electric energy consumption in KSA is affected by many economic variables, including economic growth, population growth, and income of individuals, and by energy prices, their fluctuations, and their direct and strong connection with the current successive global economic shocks [1,2].
The research problem and its challenges are summed up in the existence of a large variation in the predicted values of consumption in previous studies.
From this point, the importance of predicting electricity consumption until 2030 with greater accuracy using a composite model that helps avoid the large disparity in predictions in previous studies lies in clarifying the predicted electricity consumption for decision makers, in a way that enables them to develop and activate plans to keep pace with increasing energy demand—as proven by the results of this research—by using more renewable energy sources.

1.1. Prior Studies

In this part, we review the most important previous studies from two perspectives. First, the predicted values of electric energy consumption in KSA for the year 2030 are discussed. Second, a number of models used in prediction are reviewed.

1.1.1. The First Perspective: Predicting Consumption Values in 2030

The first study was conducted by Somani and Gracie in 2022 [1] on “Projecting Saudi sectoral electricity demand in 2030 using a computable general equilibrium model”. Their study predicted that the value of the total consumption would be 365.4 TWh in 2030. The second study was conducted by Al-Harbi and Shala in 2019 [3] on Saudi Arabia’s electricity: Energy supply and demand future challenges, in which they predicted electric energy consumption in KSA for the year 2030 to be 442 TWh. The third study was conducted by Al-Mulla in 2014 [4] on Gulf Cooperation Council (GCC) countries 2040 energy scenario for electricity generation and water desalination, in which he predicted that KSA’s electric energy consumption would reach 633.34 TWh in 2030. From the above, it is clear that there is a discrepancy in the estimates of future consumption, which justifies the need for further studies to accurately estimate consumption, especially in light of current global economic challenges.

1.1.2. The Second Perspective: Models Used in Forecasting

Liang and Liang 2017 [5] used a hybrid of the gray model and the logistic model to predict China’s electricity consumption from 2016 to 2020. It is noted that these are both distinct models. Muhammed and Podger 2005 [6] used the logistic model after modifying it by making the saturation level a function of population, electricity price, and gross national product (GNP). The modified model was used to predict electricity consumption in New Zealand in the period (2000–2020). Ogungbemi et al. 2017 [7] used the Harvey Model and a first order autoregressive model to predict the industrial electricity consumption in Nigeria from 2015 to 2029. Gharib 2022 [8] used a set of specific and stochastic models to predict the consumption of solar PV energy in China in the period 2019–2030. Al-Samman and Ahmed 2021 [9] used ten polynomial models to predict the consumption of solar PV energy in China in the period (2019–2030).
The research is based on the following hypotheses:
1.
There will be a steady and continuous increase in KSA’s electric energy consumption until 2030.
2.
It is possible to synthesize the polynomial models and the ARIMA models.
3.
Prediction accuracy using the compound models is better than prediction accuracy using a single polynomial or a single stochastic model.
The rest of the study structure consists of the following sections: Section 2 presents previous studies. Section 3 develops the general equation of the polynomial models in the following steps: calculating the parameters of the models, comparing models and selecting the best one, testing the significance of the second order polynomial model, prediction using the second order polynomial model in the sample period from 1990 to 2019, modeling residuals using the ARIMA model, calculating residuals of the residuals and making sure that they become white noise, calculating the autocorrelation function (ACF) for residuals of the residuals, testing ACF parameters one by one, synthesizing the second order polynomial outputs with the outputs of the autoregressive integrated moving average (ARIMA) model in the sample period 1990–2019, comparing the polynomial residuals with the compound model residuals in the sample period, and predicting electric energy consumption in KSA in the period from 2020 until the end of 2030. Section 4 explains and discusses the results. The final Section 5 concludes with recommendations.

2. Materials and Methods

The inductive approach, with its descriptive and analytical parts, was used [10,11], where the electric energy consumption in KSA in the period from 1990 to 2019 was described. Then, these data were analyzed, and polynomial models were created to predict consumption. Various models were evaluated and the best of them selected; then, a stochastic model was developed for the residuals of the selected polynomial model, and a compound model was created consisting of a polynomial model and a stochastic model. The significance of the selected polynomial model was tested using the F-statistical test and it passed successfully; the significance of the selected ARIMA model was tested using the Z-statistical test to ensure that its residuals became white noise.

2.1. The General Equation of the Polynomial Models Used

The general equation of these models is [12]:
y = a 0 + a 1 t + a 2 t 2 + + a n t n
where: y is the dependent variable and expresses the value of the electric energy consumption (the unit is terawatt-hours); t is the independent variable and expresses time (the unit is the year); a 0 , a 1 , a 2 , …. are the model parameters; n is the order of the polynomial and will be from 1 to 10.
We enter these parameters ( a 0 , a 1 , a 2 ). into a program called PHP, which calculates the values, and then put the parameter values in Excel to obtain the estimated y ^ .

2.2. Calculating the Parameters of the Models

The actual electricity consumption data Table A1 will be used [13] in calculating the model parameters by applying the least squares method and using Excel and Hypertext Preprocessor (PHP) programs [14].
The results are shown in the Appendix, where y ^ is the predicted consumption.
In these formulas, numbers in 17 decimal digits are used for accuracy because an important part of this research is related to the study of residuals, which are small quantities.

2.3. Comparing Models

The comparison between the ten models is made on the basis of the coefficient of determination R2, which is a strong measure of the quality of the model’s fit for the sample data. It also expresses the percentage of data interpreted by the model.
The coefficient of determination R 2 is calculated as follows: [15,16,17]:
R 2 = S S F S S T = 1 S S R S S T
where
S S T = S S R + S S F
SST is the sum of the squared deviations of the actual consumption values from their mean. SSR is the sum of the square differences between the actual electricity consumption and estimated electricity consumption. SSF is the sum of the squared deviations of the estimated electricity consumption values from their mean:
S S T = t = 1 T ( y t y ¯ ) 2
S S R = t = 1 T ( y t y ^ t ) 2
S S F = t = 1 T ( y ^ t y ¯ ) 2        
T is the size of the study sample (the number of actual electricity consumption values used in the research).
y t is the actual electricity consumption.
y ^ t is the estimated electricity consumption.
y ¯ is the average of the actual electricity consumption values.
y ^ ¯ is the average of the estimated electricity consumption values.
From Figure 1 and Table A2, there is an improvement in the value of the coefficient of determination caused by increasing the polynomial order, which becomes limited after the second order polynomial. Since the lower the polynomial order, the better it is, due to the simplicity of the model, the second order polynomial model was selected.

2.4. Testing the Significance of the Second Order Polynomial Model

In this part, the “F-statistical test” is performed [16,18] to determine the significance of the selected model as a whole.
The general equation of the model is:
y t = a + b   t + c   t 2
The test aims to determine if there is a relationship between the dependent variable yt and a subset of the variables t, t2.
The hypotheses are the null hypothesis H0 and the alternative hypothesis H1, and are as follows:
H0: a = b = c = 0
H1: at least one of a, …, c ≠ 0
Rejecting the null hypothesis H0 results in at least one of the variables t and t2 contributing significantly to the model.
Test procedures include an analysis of variance (the mean of sum of squared errors) by dividing the sum of the total squared errors (SST) into two parts. The first part is the sum of squared errors resulting from the model (regression process) SSF, and the second is the sum of the squared errors resulting due to the residual SSR, i.e.,
SST = SSF + SSR
where
S S F = i = 1 30 ( y ^ i y ^ ¯ ) 2 = 241915.7399
S S R = i = 1 30 ( y i y ^ i ) 2 = 2213.57356
The test statistic of the significance of the model is:
F 0 = S S F / k S S R / ( n p )
where k is the number of independent variables in the model (it is equal to 2 in the case under study), n is the number of observations (equal to 30 in the case under study), and p is the number of parameters in the model (equal to 3 in the case under study). In the case under study, it will be:
F 0 = 241915.7399 / 2 2213.57356 / ( 30 3 ) = 1475.380149 .
By comparing the value of   F 0 to the value
F α , k , n p
taken from the general table of the probability distribution “F” with a level of significance α = 0.05 and degrees of freedom for the numerator k = 2, and degrees of freedom for the denominator np = 27, i.e.,
F 0.05 , 2 , 27 = 3.35
By comparison, it was found that
F 0 > F 0.05 , 2 , 27
Thus, the null hypothesis H0 was rejected with a confidence level of 95% and the alternative hypothesis H1 was accepted. That is, there is a significant relationship between the dependent variable and a subset of the independent variables; therefore, we accept the model.
The model is significant with a confidence level of 95%, and it reconciles the data to an excellent degree. The conclusion is that the model can be relied upon in the prediction process.

2.5. Prediction Using the Second Order Polynomial Model in the Sample Period 1990–2019

Table A3 shows estimated electricity consumption by the second order polynomial model (calculating the actual and estimated consumption in the sample period 1990–2019) after concluding that the model as a whole is significant, with a confidence level of 95%, that it reconciles the data to an excellent degree, and that the model can be relied upon in the prediction process. The residuals, which will be inputs to the random model, can be calculated in Table A4.

2.6. Residues of the Second Order Polynomial Model

Table A4 and Figure 2 present the residues of the second order polynomial model by calculating the difference between the actual consumption and estimated consumption E = y y ^ within the period 1999–2019.
Figure 2 shows that this period (T = 26) (2015) has the biggest difference between actual and estimated consumption.

2.7. Modeling Residuals Using the ARIMA Model

By entering the residuals into the SPSS program and using the EXPERT MODELER, we obtained the optimal model for the residuals, which was ARIMA (0, 1, 0) [15,19,20] with a constant of (−0.616). This is shown in Table A5; by entering the residuals into the SPSS program and using the EXPERT MODELER, we obtained the optimal model for the residuals, which was ARIMA (0, 1, 0) [15,19,20] with a constant of (−0.616). This is a step to calculate the residuals of the residuals.

2.8. Calculating Residuals of the Residuals and Making Sure They Become White Noise

The calculation of residuals of the residuals shown in Table A6 makes sure that it becomes white noise; this is very important for the next step of calculating the autocorrelation function for residuals of the residuals.

2.9. Calculating Autocorrelation Function (ACF) [16,20,21] for Residuals of Residuals

Table A7 shows the calculation of the autocorrelation function. The autocorrelation analysis helps detect patterns and check for randomness. This is especially important when intending to use the (ARIMA) model for forecasting because it helps to determine its parameters.

2.10. Testing ACF Coefficients One by One

If the time series is white noise, then the distribution of the correlation coefficients for the sample follows the normal distribution with a mean of zero and a variance of 1/T, where T is the number of time periods; in the case under study, T = 30. This can be expressed as follows [19,20,21]:
r k ~ N ( 0 , 1 T )
Thus, the null hypothesis:
H 0 :     r k = 0
and can be tested using the parametric statistical test:
Z 0 = r k 1 T = r k T
If it is:
| Z 0 | < Z α 2
where Z α 2 is the upper percentage point α/2 in the standard normal distribution. This indicates that the coefficients of the ACF are statistically zero. That is, we accept the null hypothesis at the confidence level α = 0.05.
Note that Z α 2 = Z 0.025 = 1.96 is from the standard normal distribution table.
From Table A8, it is clear that:
| Z 0 | < Z α 2
This is for all rk coefficients. Thus, the null hypothesis can be accepted; all coefficients are statistically zero at the level of significance α = 0.05. Therefore, the residual series of the stochastic model is white noise, which means that the model is sufficient and significant with a confidence level of 95%.

2.11. Synthesizing the Second Order Polynomial Outputs with ARIMA Period 1990–2019

Table A9 presents a compound model. It consists of a polynomial model and a stochastic model. It examines the process of superposition between the second order parametric model and the ARIMA model.

2.12. Comparing Polynomial Residuals with Compound Model Residuals (Sample Period 1990–2019)

To show the improvement that occurred in the results due to the synthesizing process of the second order polynomial model and the ARIMA model, the residuals of the compound model were calculated, and we compared the sum of the residuals square and the square root of the mean of residuals of the two models shown in Table A10.
The residual sum of squares of the polynomial model RSS1 = 2213.57356.
The residual sum of squares of the compound model RSS2 = 1471.04359.
The amount of improvement resulting from the synthesizing process = RSS1 − RSS2 = 742.52997.
Improvement percentage = (RSS1 − RSS2)/RSS1 × 100 = 33.5444%.
Root-mean-square of residuals of the polynomial model RMSE1 = 8.58987.
Root-mean-square of residuals of the compound model RMSE2 = 7.00248.
Difference = RMSE1-RMSE2 = 1.58739.
Improvement percentage = (RMSE1 − RMSE2)/RMSE1 × 100 = 18.47979%.
As shown above, the synthesizing process of the two models led to a significant improvement in the accuracy of the results.
Figure 3 shows the synthesizing process of the two models that led to a significant improvement in the accuracy of the results in line y-(yhat + Ehat) compared to line E.

2.13. Predicting Electricity Consumption in KSA from 2020 to 2030 Using Compound Model

Table 1 shows the last step in predicting the electricity consumption values for the period 2020–2030 by using the compound model.

3. Explanation of the Results

The polynomials of the first order to the tenth order were experimented, where their parameters were calculated by entering the actual electricity consumption data, available in the study sample, into the PHP program on the Internet based on the method of least squares. The coefficient of determination R2 was calculated for each polynomial.
The coefficient of determination is a measure of how well the model fits the sample data; that is, the extent of convergence between the estimated electricity consumption and actual electricity consumption in the sample period. It also expresses the percentage of the data explained by the model from the total data in the sample.
By comparing the values of the coefficient of determination of the polynomials, it was found that the improvement in the value of the coefficient of determination was little for the polynomials from the third order to the tenth order. Considering that the simpler the model, the better, the second order polynomial was selected as a model for the sample data.
A graph was made showing the relationship between the polynomial and the value of determining factor. From the graph, it was found that the improvement in accuracy became imperceptible after the second order polynomial; this compendium could reduce the number of selected orders, but this was not possible before this experiment.
The significance of the model was tested using the F-statistical test at the level of significance α = 0.05, where the test was successfully passed. The residuals of the model were calculated, showing the difference between the actual electricity consumption and estimated electricity consumption in the sample period.
The ARIMA model was developed for these residuals by entering them into the SPSS statistical program and using the Expert Modeler available in the software. The ARIMA (0, 1, 0) model was produced with a constant, which is the best model for expressing residuals. The estimated residuals were calculated in the sample period using the ARIMA model, and then the residuals of the residuals were calculated, revealing the difference between the estimated residuals and the actual residuals. To make sure that the residuals of residuals became white noise, i.e., that they became statistically zero, the coefficients of the ACF were calculated and the Z-test was applied at the level of significance α = 0.05; the test was passed successfully. This means that the residuals estimated by the ARIMA model are all the remaining meaningful values in the actual residuals, and they can then be added to the electricity consumption estimated by the second order polynomial model in order to increase its accuracy.
In order to obtain the outputs of the compound model, the outputs of the ARIMA model were added to the outputs of the second order polynomial model. In order to find out the extent of improvement in the results obtained from the synthesizing process, the sum of the root-mean-square error in the sample period was calculated for each of the second order polynomial model and the compound model; a decrease of 33.5% was found, the root mean square deviation was calculated in both cases, and a decrease of 18.5% was found. All this means an improvement in the prediction accuracy using the compound model. The compound model was used in calculating the prediction for the period 2020–2030, where the value of the predicted electricity consumption in the year 2030 was 575 TWh. This result can be compared with the results of previous studies, mentioned in the current study, which were 365, 442, and 633 TWh. It should be noted that the actual electricity consumption for the year 2020 was 341 TWh, while the estimated electricity consumption in the compound model was 375 TWh, with an increase of about 10% from the actual. This can be understood as a result of the impact of the COVID-19 pandemic on actual electricity consumption in that year. As for the year 2021, the actual electricity consumption value had not yet been published at the time of preparation of this research.
There will be a steady and continuous increase in the KSA’s electricity consumption until 2030. It is also possible to make a synthesis of the polynomial models and the ARIMA models.
Prediction accuracy using compound models solves the research problem of a large discrepancy in the predicted values of consumption in previous studies. It also proves the validity of the research hypotheses, which are:
-
There will be a steady and continuous increase in the KSA’s electric energy consumption until 2030.
-
It is possible to synthesize the polynomial models and the ARIMA models.
-
Prediction accuracy using the compound models is better than prediction accuracy using a single polynomial or a single stochastic model.

4. Discussion

Energy is the mainstay of comprehensive and sustainable development in all societies and the artery of development in various economic and social fields, in addition to being one of the most important pillars of national security. Comprehensive sustainable development plans are linked to the state’s ability to provide the energy resources needed for these plans during trade-offs between state regulation, and deregulation and liberalization of energy markets [22,23].
In the case of KSA’s economy, the Ministry of Energy is working on diversifying the energy mix used in the production of electricity by increasing the share of gas and renewable energy sources in it in a more efficient and less costly manner. This is done by replacing liquid fuels with natural gas in addition to renewable energy sources. This makes the process of predicting electric energy consumption a very important component of the future planning of any economy, especially in light of current challenges and external shocks in the global economy from the impact of the COVID-19 pandemic and the Russia–Ukraine war. Therefore, the importance of this research is to achieve accurate forecasting of electric energy consumption by making an overlay between polynomial models and ARIMA models, and to confirm that the accuracy of the prediction using complex models is better than prediction accuracy using single polynomial or singular stochastic models. From this point of view, the current research provides an accurate vision of what will be the level of energy consumption in KSA up to 2030, contrary to what has been presented by previous studies, which gives a clearer picture for decision makers to make the best decision regarding the generation of electrical energy through a variety of renewable energy sources. This is in line with the sustainable development goals in KSA’s 2030 vision.
KSA’s goal of increasing and diversifying the use of renewable energy sources in generating electric energy by nearly 50% in 2030 on the one hand, and the existence of a large discrepancy in the predicted values in most of the studies conducted on this subject on the other, prompted the researchers to perform this study in order to predict KSA’s electric energy consumption for the period from 2020 to 2030. The current study was based on a sample size of 30 observations representing KSA’s actual electric energy consumption in the period from 1990 to 2019 [13]. A composite model was constructed consisting of two parts. The first part is parametric and the second part is the ARIMA model. The limits were tested from the first degree to the tenth degree, where their parameters were calculated by entering the actual consumption data, available in the study sample, into the PHP program, which is based on the method of least squares. An ARIMA model was built for these residuals by entering them into the SPSS statistical program, and with the help of the Expert Modeler available in the program, the ARIMA (0, 1, 0) model was produced with a fixed term, which is the best model for expressing the residuals.
In order to obtain the outputs of the composite model, the outputs of the ARIMA model were added to the outputs of the parametric model of the second degree. In order to find out the extent of improvement in the results resulting from the fitting process, the sum of the squared errors RSS was calculated in the sample period for each of the parametric models of the second degree and the composite model, and a decrease of 33.5% was found. All of this means an improvement in prediction accuracy using the composite model.
The composite model was used to calculate the forecast for the period 2020–2030, where the predicted consumption value in the year 2030 is 575 terawatt-hours (TWh). This result can be compared to the results of previous studies that were mentioned in this study, where they were 365, 442, and 633 TWh. It should be noted that the actual consumption for the year 2020 was 341 TWh, while the estimated consumption in the composite model was 375 TWh, an increase of about 10% from the actual. This can be understood as a result of the impact of the COVID-19 pandemic on the actual consumption in that year [24,25,26]. As for the year 2021, the actual consumption value had not yet been published by the date of preparation of this research.
Although the actual electricity consumption value for the year 2020, which was 341 TWh [27], was available, it was not used in creating the prediction models and was left to compare with the predicted values.
Therefore, the research problem is presented in the presence of a large discrepancy in the predicted values of electricity consumption in previous studies, taking these values as an indicator for sustainable development plans in the generation of electric energy and comparing them with the target for 2030. This is done in order for decision makers to be able to do what is required to keep pace with the increase in electric energy consumption by generating it through a variety of renewable energy sources and reducing carbon emissions in line with the goals of sustainable development and KSA’s vision 2030. One of the potential sources can be solar energy [28], as well as more efficient use of energy in sustainable economic development [29].
The aim of the present research is to shed light on polynomial predictive models, indicate the degree of confidence in their outputs using appropriate statistical methods, predict the electric energy consumption values using the compound model representing the period from 2020 until the end of 2030, and create a synthesis of the selected best ‘polynomial’ and its ARIMA model in order to increase prediction accuracy. It also aims to draw a future picture of electric energy consumption in KSA.
Finally, it is useful to mention that in the literature, several other forecasting approaches have been developed, such as those based on statistical and artificial neural network models, projections based on artificial intelligence [30,31,32,33], machine learning models [34,35,36], and those using novel hybrid models [37,38,39,40]. These are challenging issues for KSA, as well as other countries in the Middle East and Asian region, and are thus among the issues deserving of research in future.

5. Conclusions

The research problem was represented in forecasting the consumption of electricity consumption in the KSA in the period of the years 2020 to 2030, depending on a sample size of 30 observations, representing the actual consumption of electrical energy by Saudi Arabia in the period from 1990 to 2019. Although the actual consumption value for the year 2020 was available, which was 341 Twh, it was not used in building forecast models and was left to compare with the forecast values. The motivation for this research came from the presence of a large discrepancy in the estimated values presented in previous research.
The importance of the current research pertains to KSA’s goal of increasing and diversifying the use of renewable energy sources in generating electric energy by nearly 50% by 2030 on the one hand, and the presence of a large discrepancy in the predicted values in most of the studies conducted on this subject, on the other hand. This prompted the researchers to perform this study to predict KSA’s electric energy consumption for the period from 2020 to 2030, based on a sample size of 30 observations representing the actual electric energy consumption by KSA in the period from 1990 to 2019. Accordingly, a two-part compound model was developed. The first part is the polynomial model and the second part is the ARIMA model. Polynomials have the advantage that they can be used to express a wide range of mathematical models.
An improvement in the prediction accuracy was observed using the compound model. The compound model was used in calculating the prediction for the period 2020–2030, where the value of the predicted electricity consumption in the year 2030 was 575 TWh. This result can be compared with the results of previous studies that were mentioned in this study, where they were 365, 442, and 633 TWh. It should be noted that the actual electricity consumption for the year 2020 was 341 TWh, while the estimated electricity consumption in the compound model was 375 TWh, with an increase of about 10% from the actual.
This study recommends using the composite model as a tool that has proven its effectiveness and accuracy through the findings of this research, in order to provide accurate data and results for decision makers in KSA to determine the current and future consumption of electricity in line with the goals of sustainable development and KSA’s Vision 2030 to achieve economic benefit for KSA.

Author Contributions

Conceptualization, M.S.E.F., F.A. and Š.B.; methodology, F.D. and M.M.G.; software, F.A.; validation, M.S.E.F., F.A. and Š.B.; formal analysis, F.D. and M.M.G.; investigation, M.S.E.F. and F.A.; resources, F.A. and Š.B.; data curation, F.D. and M.M.G.; writing—original draft preparation, M.S.E.F., F.A. and F.D.; writing—review and editing, F.A. and Š.B.; visualization, M.S.E.F. and M.M.G.; supervision, F.A. and Š.B.; project administration, M.S.E.F. and M.M.G.; funding acquisition, Š.B. All authors have read and agreed to the published version of the manuscript.

Funding

This study has not received any external funding.

Data Availability Statement

Data can be provided on request.

Acknowledgments

Authors would like to thank the AGM Institute of Research & Development (AGM-IRD) Pakistan and Development Policy Practice and Research Center (DPPRC) Pakistan for their technical support in completing this study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Contains all the tables and formulas that were used in the model, including: Calculating the Parameters of the Models a; The relationship between the polynomial order and the coefficient of determination; Modeling residuals using the ARIMA model; Calculating residuals of the residuals and making sure they become white noise; Calculating autocorrelation function; Calculation of the significance test parameter of the stochastic model of the second order polynomial residuals; and a table presenting the compound model, which was created by combining a polynomial model and a stochastic model.
Table A1. Electricity Consumption Time Series (Research Sample).
Table A1. Electricity Consumption Time Series (Research Sample).
YearPeriod
T
Actual Consumption
TWh
YearPeriod
T
Actual Consumption
TWh
YearPeriod
T
Actual Consumption
TWh
1990179.9200011138.7201021240.1
1991285200112146.1201122250.1
1992393.5200213154.9201223271.1
19934102.7200314166.6201324284.1
19945106200415173.4201425311.8
19956109.9200516191.1201526338.5
19967112.2200617196.3201626345.6
19978119.1200718204.4201728355.2
19989126.2200819204.2201829359.2
199910131200920217.3201930357.4
Data source: reference [3].
Table A1 will be used [3] in calculating the model parameters by applying the least squares method and using Excel and Hypertext Preprocessor (PHP) programs [15].
Table A2. The relationship between the polynomial order and the coefficient of determination.
Table A2. The relationship between the polynomial order and the coefficient of determination.
Polynomial Order, nCoefficient of Determination, R2
10.954939
20.990933
30.990954
40.994186
50.996316
60.998177
70.998177
80.998873
90.998929
100.998929
Source: Authors’ calculations.
Table A3. Estimated electricity consumption by the second order polynomial.
Table A3. Estimated electricity consumption by the second order polynomial.
Period
T
Actual Consumption
Y
Estimated Consumption
y Hat
Period
T
Actual Consumption
Y
Estimated Consumption
y Hat
Period
T
Actual Consumption
Y
Estimated Consumption
y Hat
179.985.9985911138.7139.236421240.1243.6417
28589.0198412146.1147.374422250.1256.8965
393.592.5527613154.9156.024123271.1270.6629
4102.796.5973614166.6165.185424284.1284.941
5106101.153615173.4174.858525311.8299.7308
6109.9106.221616191.1185.043226338.5315.0322
7112.2111.801217196.3195.739527345.6330.8453
8119.1117.892518204.4206.947628355.2347.1701
9126.2124.495519204.2218.667329359.2364.0066
10131131.610120217.3230.898730357.4381.3547
Source: Table 1, reference [3], and authors’ calculations.
Table A4. Residues of the second order polynomial model.
Table A4. Residues of the second order polynomial model.
Period
T
Actual Consumption
Y
Estimated Consumption
by Second Order Polynomial
y ^
Residuals
E = y − y ^
Period
T
Actual Consumption
Y
Estimated Consumption
Second Order Polynomial y ^
Residuals
E = y − y ^
179.985.99858871−6.0985887116191.1185.04315916.05684093
28589.01983732−4.0198373217196.3195.73952780.56047215
393.592.552760610.9472393918204.4206.9475713−2.54757131
4102.796.597358576.1026414319204.2218.6672894−14.4672894
5106101.15363124.8463687820217.3230.8986823−13.5986823
6109.9106.22157853.6784214621240.1243.6417498−3.54174976
7112.2111.80120050.3987994622250.1256.8964919−6.79649194
8119.1117.89249721.2075027823271.1270.66290880.43709121
9126.2124.49546861.7045314224284.1284.9410003−0.84100032
10131131.6101146−0.6101146125311.8299.730766512.0692335
11138.7139.2364353−0.5364353226338.5315.032207423.4677926
12146.1147.3744307−1.2744307227345.6330.84532314.754677
13154.9156.0241008−1.1241007928355.2347.17011328.02988678
14166.6165.18544551.4145544729359.2364.0065781−4.80657814
15173.4174.858465−1.4584649630357.4381.3547177−23.9547177
Source: Table 1, reference [3], and authors’ calculations.
Table A5. The estimated residuals using the ARIMA model in the sample period 1990–2019.
Table A5. The estimated residuals using the ARIMA model in the sample period 1990–2019.
Period
T
Actual Consumption
Y
Estimated Consumption
y ^
Actual Residuals EEstimated Residuals E ^
179.985.99859−6.09859
28589.01984−4.01984−6.71432
393.592.552760.94724−4.63557
4102.796.597366.102640.33151
5106101.15364.846375.48691
6109.9106.22163.678424.23064
7112.2111.80120.398803.06269
8119.1117.89251.20750−0.21693
9126.2124.49551.704530.59177
10131131.6101−0.610111.08880
11138.7139.2364−0.53644−1.22584
12146.1147.3744−1.27443−1.15216
13154.9156.0241−1.12410−1.89016
14166.6165.18541.41455−1.73983
15173.4174.8585−1.458460.79883
16191.1185.04326.05684−2.07419
17196.3195.73950.560475.44111
18204.4206.9476−2.54757−0.05526
19204.2218.6673−14.46729−3.16330
20217.3230.8987−13.59868−15.08302
21240.1243.6417−3.54175−14.21441
22250.1256.8965−6.79649−4.15748
23271.1270.66290.43709−7.41222
24284.1284.941−0.84100−0.17864
25311.8299.730812.06923−1.45673
26338.5315.032223.4677911.45350
27345.6330.845314.7546822.85206
28355.2347.17018.0298914.13895
29359.2364.0066−4.806587.41416
30357.4381.3547−23.95472−5.42231
Source: Table 1, reference [3], and authors’ calculations.
Table A6. Calculation of the residuals of residuals in the sample period 1990–2019.
Table A6. Calculation of the residuals of residuals in the sample period 1990–2019.
Time
T
Actual Error
E
Estimated Error
E ^
Difference between Actual Error and Estimated Error
E     E ^
the Residuals of Residuals
1−6.09859
2−4.01984−6.714322.69448
30.94724−4.635575.58281
46.102640.331515.77113
54.846375.48691−0.64054
63.678424.23064−0.55222
70.398803.06269−2.66389
81.20750−0.216931.42443
91.704530.591771.11276
10−0.610111.08880−1.69892
11−0.53644−1.225840.68941
12−1.27443−1.15216−0.12227
13−1.12410−1.890160.76606
141.41455−1.739833.15438
15−1.458460.79883−2.25729
166.05684−2.074198.13103
170.560475.44111−4.88064
18−2.54757−0.05526−2.49231
19−14.46729−3.16330−11.30399
20−13.59868−15.083021.48434
21−3.54175−14.2144110.67266
22−6.79649−4.15748−2.63901
230.43709−7.412227.84931
24−0.84100−0.17864−0.66236
2512.06923−1.4567313.52596
2623.4677911.4535012.01429
2714.7546822.85206−8.09739
288.0298914.13895−6.10906
29−4.806587.41416−12.22074
30−23.95472−5.42231−18.53241
Source: Authors’ calculations.
Table A7. Coefficients of the two autocorrelation functions.
Table A7. Coefficients of the two autocorrelation functions.
LagCoefficients of the Autocorrelation Function
ACF
LaCoefficients of the Autocorrelation Function
ACF
10.258130.000
20.11214−0.123
3−0.08815−0.067
4−0.32216−0.116
5−0.036170.047
6−0.274180.014
7−0.174190.007
8−0.079200.034
9−0.041210.049
100.180220.110
110.136230.081
120.15624−0.021
Source: Authors’ calculations.
Table A8. Second order polynomial residuals.
Table A8. Second order polynomial residuals.
Lagrk Z 0 = r k   Ö 30 |Z0|
10.2581.4131.413
20.1120.6130.613
3−0.088−0.4820.482
4−0.322−1.7641.764
5−0.036−0.1970.197
6−0.274−1.5011.501
7−0.174−0.9530.953
8−0.079−0.4330.433
9−0.041−0.2250.225
100.1800.9860.986
110.1360.7450.745
120.1560.8540.854
130.0000.0000.000
14−0.123−0.6740.674
15−0.067−0.3670.367
16−0.116−0.6350.635
170.0470.2570.257
180.0140.0770.077
190.0070.0380.038
200.0340.1860.186
210.0490.2680.268
220.1100.6020.602
230.0810.4440.444
24−0.021−0.1150.115
Source: Authors’ calculations.
Table A9. Synthesis of the second order polynomial outputs with the ARIMA model outputs.
Table A9. Synthesis of the second order polynomial outputs with the ARIMA model outputs.
TYSecond Order Polynomial Outputs
y ^
EARIMA Model Outputs
E ^
Compound Model Outputs
  y ^ + E ^  
179.985.99859−6.09859 85.99859
28589.01984−4.01984−6.7143282.30552
393.592.552760.94724−4.6355787.91719
4102.796.597366.102640.3315196.92887
5106101.15364.846375.48691106.64051
6109.9106.22163.678424.23064110.45224
7112.2111.80120.39883.06269114.86389
8119.1117.89251.2075−0.21693117.67557
9126.2124.49551.704530.59177125.08727
10131131.6101−0.610111.0888132.6989
11138.7139.2364−0.53644−1.22584138.01056
12146.1147.3744−1.27443−1.15216146.22224
13154.9156.0241−1.1241−1.89016154.13394
14166.6165.18541.41455−1.73983163.44557
15173.4174.8585−1.458460.79883175.65733
16191.1185.04326.05684−2.07419182.96901
17196.3195.73950.560475.44111201.18061
18204.4206.9476−2.54757−0.05526206.89234
19204.2218.6673−14.4673−3.1633215.504
20217.3230.8987−13.5987−15.083215.81568
21240.1243.6417−3.54175−14.2144229.42729
22250.1256.8965−6.79649−4.15748252.73902
23271.1270.66290.43709−7.41222263.25068
24284.1284.941−0.841−0.17864284.76236
25311.8299.730812.06923−1.45673298.27407
26338.5315.032223.4677911.4535326.4857
27345.6330.845314.7546822.85206353.69736
28355.2347.17018.0298914.13895361.30905
29359.2364.0066−4.806587.41416371.42076
30357.4381.3547−23.9547−5.42231375.93239
Source: Table 1, reference [3], and authors’ calculations.
Table A10. Residues of the Compound Model.
Table A10. Residues of the Compound Model.
TY y ^ Boundary Residuals E = y y ^ E ^ y ^ + E ^ Residuals of the Compound Model
y ( y ^ + E ^ )
179.985.99859−6.09859 85.99859−6.09859
28589.01984−4.01984−6.7143282.305522.69448
393.592.552760.94724−4.6355787.917195.58281
4102.796.597366.102640.3315196.928875.77113
5106101.15364.846375.48691106.64051−0.64051
6109.9106.22163.678424.23064110.45224−0.55224
7112.2111.80120.39883.06269114.86389−2.66389
8119.1117.89251.2075−0.21693117.675571.42443
9126.2124.49551.704530.59177125.087271.11273
10131131.6101−0.610111.0888132.6989−1.6989
11138.7139.2364−0.53644−1.22584138.010560.68944
12146.1147.3744−1.27443−1.15216146.22224−0.12224
13154.9156.0241−1.1241−1.89016154.133940.76606
14166.6165.18541.41455−1.73983163.445573.15443
15173.4174.8585−1.458460.79883175.65733−2.25733
16191.1185.04326.05684−2.07419182.969018.13099
17196.3195.73950.560475.44111201.18061−4.88061
18204.4206.9476−2.54757−0.05526206.89234−2.49234
19204.2218.6673−14.46729−3.1633215.504−11.304
20217.3230.8987−13.59868−15.083215.815681.48432
21240.1243.6417−3.54175−14.2144229.4272910.67271
22250.1256.8965−6.79649−4.15748252.73902−2.63902
23271.1270.66290.43709−7.41222263.250687.84932
24284.1284.941−0.841−0.17864284.76236−0.66236
25311.8299.730812.06923−1.45673298.2740713.52593
26338.5315.032223.4677911.4535326.485712.0143
27345.6330.845314.7546822.85206353.69736−8.09736
28355.2347.17018.0298914.13895361.30905−6.10905
29359.2364.0066−4.806587.41416371.42076−12.22076
30357.4381.3547−23.95472−5.42231375.93239−18.53239
Source: Table 1, reference [3], and authors’ calculations.
Calculating the Parameters of the Models
The first order polynomial model:
y ^ = 41.19057471264368 + 10.184694104560624 t
The second order polynomial model:
y ^ = 83.48901477832513 + 2.253736592245352 t + 0.25583733910694423 t 2
The third order polynomial model:
y ^ = 84.83395730706076 + 1.7721439626681563 t + 0.2940459336732974 t 2 0.0008216902057280256 t 3
The fourth order polynomial model:
y ^ = 63.33479221927498 + 13.867819507823924 t 1.395524075804146 t 2 + 0.08300972446366839 t 3 0.0013521195914418776 t 4
The fifth order polynomial model:
y ^ = 86.16489832007073 4.099943959460825 t + 2.400591494097555 t 2 0.23516904422015494 t 3 + 0.010093454991176316 t 4 0.00014768483332410572 t 5
The sixth order polynomial model:
y ^ = 57.75731211317419 + 24.83996800501148 t 6.041637787034118 t 2 + 0.8068538789773372 t 3 0.051680440105507985 t 4 + 0.0015929760862161242 t 5 0.00001871678408107774 t 6
The seventh order polynomial model:
y ^ = 57.643512858955226 + 24.9819845514518 t 6.095377306021259 t 2 + 0.815862184315436 t 3 0.05245280907224429 t 4 + 0.0016282585012546156 t 5 0.00001953344097961959 t 6 + 0.00000000752679169163 t 7
The eighth order polynomial model:
y ^ = 90.76698650674774 23.644465498609527 t + 16.69659216991278 t 2 4.0949632466064285 t 3 + 0.5120294538898259 t 4 0.035190806450623516 t 5 + 0.0013459517043590157 t 6 0.00002680117056660032 t 7 + 0.00000021619917224429 t 8
The ninth order polynomial model:
y ^ = 77.18611694154131 0.8930437908518888 t + 3.950515271184289 t 2 0.7057407804364993 t 3 + 0.015718240582896285 t 4 + 0.007782396358829107 t 5 0.0009080646201480251 t 6 + 0.00004348308707838665 t 7 0.00000098162131905363 t 8 + 0.00000000858652681934 t 9
The tenth order polynomial model:
y ^ = 77.67766416789804 1.8113357174773226 t + 4.54824364199519 t 2 0.8957648439793898 t 3 + 0.04986512462338578 t 4 + 0.004043738455536497 t 5 0.0006497488448332843 t 6 + 0.00003217491257888407 t 7 0.00000067763036404404 t 8 + 0.00000000401026876840 t 9 0.00000000002952424549 t 10

References

  1. Soummane, S.; Ghersi, F. Projecting Saudi sectoral electricity demand in 2030 using a computable general equilibrium model. Energy Strat. Rev. 2022, 39, 100787. [Google Scholar] [CrossRef]
  2. Omran, K.M.; Zaalan, R.A.-I. Using Some Statistical Methods to Predict the Consumption of Electric Energy in the Kingdom of Saudi Arabia; Economic Sciences; University of Basra, College of Administration and Economics: Basra, Iraq, 2012; Volume 8, p. 29. [Google Scholar]
  3. Al Harbi, F.; Csala, D. Saudi Arabia’s Electricity: Energy Supply and Demand Future Challenges; Lancaster University: Lancaster, UK, 2019. [Google Scholar]
  4. Almulla, Y. Gulf Cooperation Council (GCC) Countries 2040 Energy Scenario for Electricity Generation and Water Desalination; KTH Industrial Engineering and Management: Stockholm, Sweden, 2014. [Google Scholar]
  5. Liang, J.; Liang, Y. Analysis and Modeling for China’s Electricity Demand Forecasting Based on a New Mathematical Hybrid Method. Information 2017, 8, 33. [Google Scholar] [CrossRef] [Green Version]
  6. Mohamed, Z.; Bodger, P.S. A variable asymptote logistic (VAL) model to forecast electricity consumption. Int. J. Comput. Appl. Technol. 2005, 22, 65–72. [Google Scholar] [CrossRef] [Green Version]
  7. Oluropo, O.E.; Archibong, E.J.; Affia, N.J. Time Series Analysis of Industrial Electricity Consumption in Nigeria Using Harvey Model and Autoregressive Model. Int. J. Energy Power Eng. 2017, 6, 40–46. [Google Scholar] [CrossRef] [Green Version]
  8. Ghareeb, M.M. Forecasting Solar PV Consumption in China by Compound Application of Specific and Random Models; Zagazig University: Zagazig, Egypt, 2022. [Google Scholar]
  9. Al-Samman, A.H.; Mohamed, M.; Ahmed, G. Predicting the Consumption of Solar Photovoltaic Energy in China Using ten Frontiers. Sci. J. Fac. Commer. Sect. 2021. [Google Scholar]
  10. Badawi, A.R. Scientific Research Methods, 3rd ed.; Publications Agency: Kuwait, 1977. [Google Scholar]
  11. Sarhan, M.; Al-Mahmoudi, A. Scientific Research Methods, 3rd ed.; Dar Al-Kutub: Sana’a, Yemen; Al-Azhar University: Cairo, Egypt, 2019. [Google Scholar]
  12. BP. Statistical Review of World Energy, 69th ed.; Microsoft Excel; BP: London, UK, 2020. [Google Scholar]
  13. Ali Shalaby, O. Principles of Mathematics for Traders; Part One; Faculty of Commerce; Zagazig University: Zagazig, Egypt, 2001; p. 86. [Google Scholar]
  14. Polynomial Regression. Available online: http://polynomialregression.drque.net/online.php (accessed on 22 December 2022).
  15. Montgomery, D.C.; Jennings, C.L.; Kulahci, M. Introduction to Time Series Analysis and Forecasting; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015. [Google Scholar]
  16. Al-Rahman, A.; Al-Salihi, I. Introduction to Statistics; Faculty of Commerce, Zagazig University: Zagazig, Egypt, 1997; pp. 176–177. [Google Scholar]
  17. Casmir, L.C. Commercial Statistics; Mustafa, M.G., Translator; International House for Cultural Investments: Cairo, Egypt, 2003. [Google Scholar]
  18. Chong, J. Powerful Forecasting with MS Excel; XL Pert Enterprise: Perth, WA, USA, 2010. [Google Scholar]
  19. Abdel, A.; Sabry, R. Applied Statistics with SPSS System; Methodology House for Publishing and Distribution: Amman, Jordan, 2015; p. 65. [Google Scholar]
  20. Mohamed, A.; Samie, A. Econometrics and Statistical Analysis of Economic Relationships—A Modern Introduction Using Windows SPSS; Faculty of Commerce, Zagazig University: Zagazig, Egypt, 2009. [Google Scholar]
  21. Borgonovo, E. Sensitivity Analysis: An Introduction for the Management Scientist; Springer International Publishing: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  22. Bojnec, Š.; Križaj, A. Electricity Markets during the Liberalization: The Case of a European Union Country. Energies 2021, 14, 4317. [Google Scholar] [CrossRef]
  23. Bahmanyar, A.; Estebsari, A.; Ernst, D. The impact of different COVID-19 containment measures on electricity consumption in Europe. Energy Res. Soc. Sci. 2020, 68, 101683. [Google Scholar] [CrossRef] [PubMed]
  24. Abulibdeh, A. Modeling Electricity Consumption Patterns during the COVID-19 Pandemic across Six Socioeconomic Sectors in the State of Qatar. Energy Strategy Rev. 2021, 38, 100733. [Google Scholar] [CrossRef]
  25. Jaramillo, M.; Carrión, D. An Adaptive Strategy for Medium-Term Electricity Consumption Forecasting for Highly Unpredictable Scenarios: Case Study Quito, Ecuador during the Two First Years of COVID-19. Energies 2022, 15, 8380. [Google Scholar] [CrossRef]
  26. BP. Statistical Review of World Energy, 70th ed.; BP: London, UK, 2021. [Google Scholar]
  27. Scaron; Tefan, B.; Drago, P. Can Solar Electricity Production and Use Help to Promote Economic Development in Africa? Afr. J. Bus. Manag. 2011, 5, 6949–6956. [Google Scholar]
  28. Bojnec, Š.; Papler, D. Economic efficiency, energy consumption and sustainable development. J. Bus. Econ. Manag. 2011, 12, 353–374. [Google Scholar] [CrossRef] [Green Version]
  29. da Silva, F.L.C.; da Costa, K.; Rodrigues, P.C.; Salas, R.; López-Gonzales, J.L. Statistical and Artificial Neural Networks Models for Electricity Consumption Forecasting in the Brazilian Industrial Sector. Energies 2022, 15, 588. [Google Scholar] [CrossRef]
  30. Jin, N.; Yang, F.; Mo, Y.; Zeng, Y.; Zhou, X.; Yan, K.; Ma, X. Highly Accurate Energy Consumption Forecasting Model Based on Parallel LSTM Neural Networks. Adv. Eng. Inform. 2022, 51, 101442. [Google Scholar] [CrossRef]
  31. Yang, S.; Deng, Z.; Li, X.; Zheng, C.; Xi, L.; Zhuang, J.; Zhang, Z.; Zhang, Z. A novel hybrid model based on STL decomposition and one-dimensional convolutional neural networks with positional encoding for significant wave height forecast. Renew. Energy 2021, 173, 531–543. [Google Scholar] [CrossRef]
  32. Gao, T.; Niu, D.; Ji, Z.; Sun, L. Mid-term electricity demand forecasting using improved variational mode decomposition and extreme learning machine optimized by sparrow search algorithm. Energy 2022, 261, 5328. [Google Scholar] [CrossRef]
  33. Shin, S.-Y.; Woo, H.-G. Energy Consumption Forecasting in Korea Using Machine Learning Algorithms. Energies 2022, 15, 4880. [Google Scholar] [CrossRef]
  34. Lee, M.H.L.; Ser, Y.C.; Selvachandran, G.; Thong, P.H.; Cuong, L.; Son, L.H.; Tuan, N.T.; Gerogiannis, V.C. A Comparative Study of Forecasting Electricity Consumption Using Machine Learning Models. Mathematics 2022, 10, 1329. [Google Scholar] [CrossRef]
  35. Kauffman, R.J.; Kim, K.; Lee, S.Y.T.; Hoang, A.P.; Ren, J. Combining machine-based and econometrics methods for policy analytics insights. Electron. Commer. Res. Appl. 2017, 25, 115–140. [Google Scholar] [CrossRef]
  36. Wang, Q.; Li, S.; Li, R. Forecasting energy demand in China and India: Using single-linear, hybrid-linear, and non-linear time series forecast techniques. Energy 2018, 161, 821–831. [Google Scholar] [CrossRef]
  37. Fan, G.F.; Wei, X.; Li, Y.T.; Hong, W.C. Forecasting electricity consumption using a novel hybrid model. Sustain. Cities Soc. 2020, 61, 102320. [Google Scholar] [CrossRef]
  38. Lee, J.; Cho, Y. National-scale electricity peak load forecasting: Traditional, machine learning, or hybrid model. Energy 2022, 239, 122366. [Google Scholar] [CrossRef]
  39. Dagar, V.; Ahmed, F.; Waheed, F.; Bojnec, Š.; Khan, M.K.; Shaikh, S. Testing the Pollution Haven Hypothesis with the Role of Foreign Direct Investments and Total Energy Consumption. Energies 2022, 15, 4046. [Google Scholar] [CrossRef]
  40. Kousar, S.; Sabir, S.A.; Ahmed, F.; Bojnec, Š. Climate Change, Exchange Rate, Twin Deficit, and Energy Inflation: Application of VAR Model. Energies 2022, 15, 7663. [Google Scholar] [CrossRef]
Figure 1. The relationship between the polynomial order (n) and the coefficient of determination R2 according to the data of Table A2. Source: Authors’ calculations.
Figure 1. The relationship between the polynomial order (n) and the coefficient of determination R2 according to the data of Table A2. Source: Authors’ calculations.
Energies 16 00506 g001
Figure 2. Second order polynomial residuals versus time intervals. Source: Table 1, reference [3], and authors’ calculations.
Figure 2. Second order polynomial residuals versus time intervals. Source: Table 1, reference [3], and authors’ calculations.
Energies 16 00506 g002
Figure 3. Comparing the residuals of the polynomial model E with the residuals of the compound model y-(yhat + Ehat). Source: Authors’ calculations.
Figure 3. Comparing the residuals of the polynomial model E with the residuals of the compound model y-(yhat + Ehat). Source: Authors’ calculations.
Energies 16 00506 g003
Table 1. Electricity consumption prediction values for the period 2020–2030.
Table 1. Electricity consumption prediction values for the period 2020–2030.
TYearEstimated Consumption by the Second Order PolynomialResiduals Estimated by the ARIMA ModelEstimated Consumption by the Compound Model
312020399.2145−24.5705374.64405
322021417.586−25.1862392.39983
332022436.4692−25.8019410.6673
342023455.864−26.4176429.44637
352024475.7705−27.0334448.73714
362025496.1887−27.6491468.53961
372026517.1186−28.2648488.85378
382027538.5601−28.8806509.67955
392028560.5133−29.4963531.01702
402029582.9782−30.112552.8662
412030605.9548−30.7277575.22707
Source: Authors’ calculations.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fahmy, M.S.E.; Ahmed, F.; Durani, F.; Bojnec, Š.; Ghareeb, M.M. Predicting Electricity Consumption in the Kingdom of Saudi Arabia. Energies 2023, 16, 506. https://doi.org/10.3390/en16010506

AMA Style

Fahmy MSE, Ahmed F, Durani F, Bojnec Š, Ghareeb MM. Predicting Electricity Consumption in the Kingdom of Saudi Arabia. Energies. 2023; 16(1):506. https://doi.org/10.3390/en16010506

Chicago/Turabian Style

Fahmy, Marwa Salah EIDin, Farhan Ahmed, Farah Durani, Štefan Bojnec, and Mona Mohamed Ghareeb. 2023. "Predicting Electricity Consumption in the Kingdom of Saudi Arabia" Energies 16, no. 1: 506. https://doi.org/10.3390/en16010506

APA Style

Fahmy, M. S. E., Ahmed, F., Durani, F., Bojnec, Š., & Ghareeb, M. M. (2023). Predicting Electricity Consumption in the Kingdom of Saudi Arabia. Energies, 16(1), 506. https://doi.org/10.3390/en16010506

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