Forecasting of Electrical Generation Using Prophet and Multiple Seasonality of Holt–Winters Models: A Case Study of Kuwait

Electrical generation forecasting is essential for management and policymakers due to the crucial data provided for resource planning. This research employs the Prophet model with single and multiple regressors to forecast the electricity generation in Kuwait from 2020 to 2030. In addition, multiple seasonality Holt–Winters models were utilized as a benchmark for comparative analysis. The accuracy, generalization, and robustness of the models were assessed based on different statistical performance metrics. The triple seasonality Holt–Winters model achieved superior performance compared with the other models with R2 = 0.9899 and MAPE = 1.76%, followed by the double seasonality Holt–Winters model with R2 = 0.9893 and MAPE = 1.83%. Moreover, the Prophet model with multiple regressors was the third-best performing model with R2 = 0.9743 and MAPE = 2.77%. The forecasted annual generation in the year 2030 resulted in 92,535,555 kWh according to the best performing model. The study provides an outlook on the mediumand long-term electrical generation. Furthermore, the impact of fuel cost is investigated based on the five forecasting models to provide an insight for Kuwait’s policymakers.


Introduction
Accurate electrical generation forecasting is essential for the management and policymakers of local and national power plants. Electric power generation is one of the complex processes where many parameters are involved and need to be optimized to deliver the electricity continuously and efficiently to the whole population [1]. On the other hand, power generation is faced with many challenges, such as environmental legislation, fluctuations in fuel prices, and the need to optimize the available resources. Moreover, different sources for power generation are used nowadays. These sources include nuclear energy, fossil fuel such as coal, natural gas, crude oil, other gasses, and renewable energy such as hydropower, solar, wind, municipal solid waste, and geothermal. Therefore, appropriate forecasting is a need for these power plants in order to support continued prosperity. Table 1. The growing population and per capita consumption in three decades [21]. One of the main points proposed [28] for fixing the high generation's need is to tackle the subsidized energy prices, which inspire greater electricity. Another viable idea is the use of dynamic pricing [29], which can change the pattern or reduce the quantity of people's electricity consumption. Increasing electricity tariffs is expected to be a significant contributor to change the consumption pattern of people. However, there are some difficulties in implementing these reforms due to their political nature [28]. Alajmi [30] used energy audit techniques and implemented non-retrofitting and retrofitting on a two-story educational facility that results in more than 50% annual total saving for governmental buildings. The Ministry of electrical and water, especially with the need to mitigate climate change, aims to reach a generation of 15% of renewable energy by 2035 [31]. Moreover, the government is encouraging people to use renewable energy and new technology, which might lead to decentralize the power systems and might also reduce the effect on government [22]. The introduction of new sustainable buildings and cities that implement renewable energy and efficiently use the electricity, especially with new residential cites in Kuwait, is feasible and efficient [22]. Digitalization through the use of smart meters and smart grid is also believed to provide more reduction and sustainable use of the electricity and are discussed by the Ministry of electrical and water and should come into action [22]. However, remedies for the increase in the generation and their consequences are out of this study's scope.
The continuous increase in electricity generation and the different plans of government necessitate carrying out more studies in forecasting. Studies are needed to estimate the total electricity generation for the coming years to assist the policymakers in achieving the proper decisions regarding the future. Several studies were carried out regarding long-term forecasting for Kuwait. However, most of them were about peak load forecasting [18,[32][33][34][35][36][37]. Alajmi and Phelan [32] used a bottom-up approach to create a baseline for the residential sector's end-use energy profile in Kuwait until 2040. Recently, Alhajeri et al. [35] investigate the effects of COVID-19 on electrical consumption from 1st March to 30th May. Atalla and Hunt [37] studied high residential electricity demand drivers in the Gulf Cooperation Council countries (GCC) using the time series model and suggested some recommendations to reduce the demand.
To the best of authors' knowledge, limited studies have investigated the Prophet model's use in forecasting the electricity generation. The originality of this study is that it explores the performance of prophet model with multi-regressors in long-term forecasting of electricity generation. The multi-regressors are the historical daily maximum and minimum load, temperature and population. In addition, comparative quantitative analyses with multi-seasonality Holt-Winters model is presented to assess the forecasting performance. Furthermore, the generalizability and robustness of the Prophet and Holt-Winters methods for forecasting long-term electricity generation are explored and presented. Various factors of the future are changing, such as population, weather temperatures, and peak loads, directly or indirectly related to electricity generation were taken into account in the forecasting. The study also provides and discusses the estimation of the expected future fuel consumption cost until 2030.
The paper is structured as follows. Section 2 presents a description of the Prophet and the multi-seasonality Holt-Winter models. Section 3 illustrates the results and analyses based on the performance indicators. In addition, Section 3 provides insight for policymakers on the implications of electrical generation on fuel costs.

Methodologies
The electrical generation's real data between January 2015 and July 2020 was attained from the Ministry of Electricity and Water (MEW) is used in this research. The data is plotted in Figure 1, showing the daily generation for approximately six years. A vivid seasonality of the data is observed with an annual maximum generation in the hot summer, especially June to August. There is an increasing trend of the generation, as noticed from Figure 1, due to the growth of population, the building of new cities, as well as economic development, and this is expected to continue in the upcoming years.
suggested some recommendations to reduce the demand.
To the best of authors' knowledge, limited studies have investigated the Prophet model's use in forecasting the electricity generation. The originality of this study is that it explores the performance of prophet model with multi-regressors in long-term forecasting of electricity generation. The multi-regressors are the historical daily maximum and minimum load, temperature and population. In addition, comparative quantitative analyses with multi-seasonality Holt-Winters model is presented to assess the forecasting performance. Furthermore, the generalizability and robustness of the Prophet and Holt-Winters methods for forecasting long-term electricity generation are explored and presented. Various factors of the future are changing, such as population, weather temperatures, and peak loads, directly or indirectly related to electricity generation were taken into account in the forecasting. The study also provides and discusses the estimation of the expected future fuel consumption cost until 2030.
The paper is structured as follows. Section 2 presents a description of the Prophet and the multi-seasonality Holt-Winter models. Section 3 illustrates the results and analyses based on the performance indicators. In addition, Section 3 provides insight for policymakers on the implications of electrical generation on fuel costs.

Methodologies
The electrical generation's real data between January 2015 and July 2020 was attained from the Ministry of Electricity and Water (MEW) is used in this research. The data is plotted in Figure 1, showing the daily generation for approximately six years. A vivid seasonality of the data is observed with an annual maximum generation in the hot summer, especially June to August. There is an increasing trend of the generation, as noticed from Figure 1, due to the growth of population, the building of new cities, as well as economic development, and this is expected to continue in the upcoming years.
The data consists of five years with daily observations for electrical generation from January 2015 to May 2020 with 1978 data points. Taking a closer look at Figure 1, a yearly seasonality occurring within July-August, can be observed. During these months, the high electricity demand is due to Kuwait's hot summer, with temperatures reaching 55 °C. The data also reveals a weekly seasonality with low demand on Fridays and Saturdays each week due to the factories' and industries' closure on weekends. Furthermore, the third seasonality corresponds to the year's meteorological seasons and will be addressed hereafter as quarters to avoid the confusion with methods' seasonality patterns.  The data consists of five years with daily observations for electrical generation from January 2015 to May 2020 with 1978 data points. Taking a closer look at Figure 1, a yearly seasonality occurring within July-August, can be observed. During these months, the high electricity demand is due to Kuwait's hot summer, with temperatures reaching 55 • C. The data also reveals a weekly seasonality with low demand on Fridays and Saturdays each week due to the factories' and industries' closure on weekends. Furthermore, the third seasonality corresponds to the year's meteorological seasons and will be addressed hereafter as quarters to avoid the confusion with methods' seasonality patterns.
In this study, long-term forecasting for electrical generation for the coming ten years is aimed. Therefore, a search for models that comprise data with time-series nature with different seasonality resulted in choosing the relatively new Prophet model from Facebook and the Holt-Winters model for this work.
The relatively new Prophet model possesses some features such as tunable parameters by an expert during the analysis, which allows for fine-tuning and the ability to reach robust forecasting. The model's automatic tuning does not require a full understanding of the underlying model, making it easier to use. Limited studies are reported to use this model and proved to achieve vigorous results [18,38]. Whereas the Holt-Winters model, with different variations, is used extensively for long-term forecasting and providing excellent results [39,40].
One of the Prophet model's points of superiority is its ability to deal with missing data and reach outstanding results. The Prophet model's flexibility for complex data can be achieved using an analyst-in-the-loop feature by integrating multiple-seasonality of different periods. The availability of single and multivariate forecasting in a straightforward manner adds a new strength for the Prophet model. On the other hand, Holt-Winters deals only with univariate forecasting. There exist different Holt-Winters model variations in the literature where it was generalized to include double and triple seasonalities [41,42]. The Holt-Winter model was then generalized and introduced in literature to include n-seasonalities and was therefore known as nHWT [38]. Initialization methods of nHWT were proven to have a crucial impact on the accuracy of the forecast [19].
In this research, the multivariate Prophet forecasting model is utilized to forecast the annual load generation needed in Kuwait for the upcoming 10 years and is compared with the univariate nHWT forecasting model with initialization methods. In the next sections, the implementation of the various Prophet and Holt-Winters forecasting methods is presented in detail.

Prophet Forecasting Method
A Prophet model is used with simple or complicated time-series data that includes single or multiple seasonality, holidays, and data trends. Multiple seasonality can include different patterns, such as days, weeks, months, and years. As reported by Taylor and Letham [17], the mathematical representation of the decomposed time series model: where g(t) denotes the data trend function, s(t) denotes the seasonality, and h(t) denotes holidays effect that can be added within specific points of the data and as an extra regressor. The error term, ε t , denotes any distinctive features of the data that are not fitted by the model. Prophet trend function, g(t), can be signified by a piecewise linear growth model or a saturating growth model. Since the electrical generation does not exhibit a saturating growth, a piecewise linear growth model is utilized as: where k represents the growth rate, δ represents adjustment rate, m represents the offset parameter, and γ is the trend changepoints, s j , and is set as −s j δ j , with a(t) defined as: The change points allow the analyst in the loop to adjust the resulting forecast based on a previous experience. Therefore, the trend of the forecast can be hence fine-tuned and results in an enhanced forecast. The seasonality function s(t) can be modelled by Fourier series to represent daily, weekly, and yearly seasonality as well as to incorporate more complex seasonality patterns by higher order Fourier series. The seasonality function is hence written as: where P is assumed to be 365.25 for yearly seasonality pattern. Furthermore, Prophet model allows adding extra regressors to enhance the forecast results. The holiday effects, as example, can be incorporated using h(t) function and define a list of holiday dates as a matrix of regressors Z(t) defined as: with D as the set of holiday dates, κ ∼ Normal 0, v 2 with v as the holiday smoothing parameter. In this work, additional regressors were utilized and defined as the temperature, the maximum load, the minimum load, and the population to result in a more reliable forecast that reflects the growth in population and Kuwait's average temperatures.

Holt-Winters Forecasting Model
Holt-Winters models are well-established models that have two seasonal variations: additive and multiplicative types. The additive method is appropriate for fitting a time series data with a constant seasonal variation. Whereas multiplicative method is appropriate for a time series with an increasing seasonal pattern relative to the data level. The generalized Holt-Winters model was reported in literature [19,41,43]. In this research, the generalized Holt-Winters model with multiple seasonality and initialization methods is used for forecasting the electrical generation of Kuwait until 2030.
The Holt-Winters model consists of forecast equation and smoothing equations of the level, trend, and seasonality of the time series. Assuming that S t represents the seasonality, T t to represent the trend and that I t corresponds seasonality, the additive trend multiplicative seasonality HW model can be defined as: where S t represents the level, T t as the trend, I i t to correspond to multiple seasonality andX t (k) is the k-step ahead forecast. The smoothing parameters are defined as α for the level smoothing, γ as the trend smoothing parameter, and δ i as the smoothing parameters of each seasonal pattern with cycle length of S i . The term ϕ k AR is an adjustment for the first autocorrelation error. There exist various initialization methods for the level, trend, and seasonality in order to enhance the forecast accuracy [19]. The trend is initialized as T 0 = 0 for an additive trend. Whereas the level equation is initialized by the moving average value of the dataset and is obtained by: Appl. Sci. 2020, 10, 8412 The seasonality is initialized by the method presented by Brockwell and Davis [39] and adopted by the National Institute of Standards (NIST) which depends on calculating the weights of the data series against the multiple seasonality pattern values. The proposed method is presented in the following steps: Step 1: Compute the yearly average as A i m for each seasonality of length s m and has a pattern of n q times in the dataset Step 2: Divide the observations by the yearly averages as: . .
Step 3: Write each seasonality as: Step 4: Write the seasonal indices as: The simulation parameters are presented in Table 2 with the values of level, trend and seasonality smoothing parameters as α, γ, δ i respectively. The seasonality smoothing parameters are presented as δ y , δ q and δ w that corresponds to the yearly, quarterly, and weekly seasonality smoothing parameter values.

Performance Indicators
In order to evaluate the performance of each model in terms of accuracy, various statistical metrics were adopted, such as root mean square error (RMSE), the mean absolute percentage error (MAPE), coefficient of determination (R 2 ), mean absolute error (MAE), and coefficient of variation of root mean square error (CVRMSE) that can be expressed as in the following equations: Appl. Sci. 2020, 10, 8412 where y,ŷ, andy represent the measured, predicted, and averaged values respectively.

Results and Discussion
Different models are employed to forecast the electrical generation for Kuwait up to 2030 using real data from the Kuwait Ministry of electrical and water. The simulations were carried out within Rstudio and Matlab software environments. The period between January 2015 to December 2018 was used for estimation and adjustment, while the period from January 2019-May 2020 was used for the validation, and the forecasting horizon is between June 2020 to December 2030.
Single (PSR) and multiple (PMR) regressors were used for Prophet model, and three methods of seasonality Holt-Winters models: single (HWSS), double (HWDS), and triple (HWTS). Each model's performance is assessed from different angles by accuracy indicators. The future electrical generation forecasting of Kuwait until the year 2030 using the five models is plotted in Figure 2a,b. The models are assessed in the following section and then discussed, and different fuel cost scenarios are investigated.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 19 where , , and represent the measured, predicted, and averaged values respectively.

Results and Discussion
Different models are employed to forecast the electrical generation for Kuwait up to 2030 using real data from the Kuwait Ministry of electrical and water. The simulations were carried out within Rstudio and Matlab software environments. The period between January 2015 to December 2018 was used for estimation and adjustment, while the period from January 2019-May 2020 was used for the validation, and the forecasting horizon is between June 2020 to December 2030.
Single (PSR) and multiple (PMR) regressors were used for Prophet model, and three methods of seasonality Holt-Winters models: single (HWSS), double (HWDS), and triple (HWTS). Each model's performance is assessed from different angles by accuracy indicators. The future electrical generation forecasting of Kuwait until the year 2030 using the five models is plotted in Figure 2a

Model Assessments
Several performance indicators were employed to investigate the reliability of the forecasting of the used models. These indicators are discussed in this section to assess the accuracy and generalization of the used models. Five accuracy performance statistical metrics were employed in this study. These metrics include MAPE, MAE, RMSE, CVRMSE, and R 2 . Each of the metrics shed light on one angle of the accuracy of the data. MAPE can be considered one of the utmost used tools for evaluating the accuracy of models. MAE reveals the difference between the estimated value and the real value using the absolute error. RMSE assesses the variability of model response regarding variance and sensitivity to large errors. CVRMSE standardizes the forecasted error and provides a unit-less metric that evaluates the variability of the errors between predicted and real values.  [2]. All models achieved a low MAPE value of less than 5%, which indicates the models' high accuracy. It is benchmarked [40] that When MAPE is less than 10%, it is considered a highly accurate model. Highly accurate models were then categorized into four levels, with the best model denoted by I level when MAPE is ≤1.2%. Other levels are as follows: II (1.2-2.8%), III (2.8-4.6%), and IV (4.6-10%) [2]. Three of the used models, namely the multiple Prophet model and the double and triple seasonality of Holt-Winters models, fall in the II level of this rating, whereas the other two, namely single Prophet and single seasonality of Holt-Winters model, fall in the III level. The triple seasonality Holt-Winters model outperforms all the other models in MAPE value. Superiority is also observed for both double and triple seasonality Holt-Winters models in MAPE values compared to other models. The MAPE values were 1.76% and 1.83% for triple and double seasonality Holt-Winters models, respectively. The coefficients of determination, R 2 , of the five models are plotted with excellent values in Figure 3, showing small differences among the forecasted and real data values. The best-fitting was for the Holt-Winters model's triple seasonality with R 2 = 0.9899, then to the Holt-Winters model's double seasonality R 2 = 0.9893 followed by multiple Prophet model with R 2 = 0.9743. The least values of R 2 were 0.9709 and 0.9641 for the single Prophet model and single seasonality of the Holt-Winters model, respectively. The MAE of the best model, which is the triple seasonality of Holt-Winters, was found to be 47.82, whereas the value of MAE for the least model, which was the single seasonality Holt-Winters model, was 127.82, which is approximately three times of the triple model. The RMSE value of the Holt-Winters model's triple seasonality was also the lowest with a value of 67.05 among the other models, followed by the double seasonality with a value of 75.09, whereas the single seasonality of the Holt-Winters model was the maximum with a value of 165.06. Similarly, the CVRMSE percentage of the Holt-Winters model's triple seasonality was the least compared to other models, followed by the double seasonality. Therefore, based on the statistical performance metrics, it is vividly clear that the triple seasonality of the Holt-Winters model has a superior fitting and better accuracy compared to other used models followed by the double seasonality of the Holt-Winters model.
A model's generalization is measured by a model's capability to estimate samples beyond the training zone. The five models' relative error variations are shown in Figure 4 with a maximum relative error value of less than 40% for all models. The triple seasonality of the Holt-Winters model values was less than 21%, and the maximum was for the single seasonality of the Holt-Winters model. The worst relative error was for the single seasonality of the Holt-Winters model.  The RMSE value of the Holt-Winters model's triple seasonality was also the lowest with a value of 67.05 among the other models, followed by the double seasonality with a value of 75.09, whereas the single seasonality of the Holt-Winters model was the maximum with a value of 165.06. Similarly, the CVRMSE percentage of the Holt-Winters model's triple seasonality was the least compared to other models, followed by the double seasonality. Therefore, based on the statistical performance metrics, it is vividly clear that the triple seasonality of the Holt-Winters model has a superior fitting and better accuracy compared to other used models followed by the double seasonality of the Holt-Winters model.
A model's generalization is measured by a model's capability to estimate samples beyond the training zone. The five models' relative error variations are shown in Figure 4 with a maximum relative error value of less than 40% for all models. The triple seasonality of the Holt-Winters model values was less than 21%, and the maximum was for the single seasonality of the Holt-Winters model. The worst relative error was for the single seasonality of the Holt-Winters model. The RMSE value of the Holt-Winters model's triple seasonality was also the lowest with a value of 67.05 among the other models, followed by the double seasonality with a value of 75.09, whereas the single seasonality of the Holt-Winters model was the maximum with a value of 165.06. Similarly, the CVRMSE percentage of the Holt-Winters model's triple seasonality was the least compared to other models, followed by the double seasonality. Therefore, based on the statistical performance metrics, it is vividly clear that the triple seasonality of the Holt-Winters model has a superior fitting and better accuracy compared to other used models followed by the double seasonality of the Holt-Winters model.
A model's generalization is measured by a model's capability to estimate samples beyond the training zone. The five models' relative error variations are shown in Figure 4 with a maximum relative error value of less than 40% for all models. The triple seasonality of the Holt-Winters model values was less than 21%, and the maximum was for the single seasonality of the Holt-Winters model. The worst relative error was for the single seasonality of the Holt-Winters model.  The models' robustness is assessed by incorporating different Gaussian white noise intensities into the real tested data similar to relevant previous studies [44][45][46]. The data were divided into a training dataset and several testing sets with each set corresponding to a different noisy intensity. Four different intensities of 20%, 40%, 60%, and 80% were utilized to investigate the robustness of the five models. The models' robustness was evaluated through the accuracy of the coefficient of determination R 2 by studying the random distributed Gaussian white noise on the training dataset. Table 4 shows the variance of R 2 affected by different noise intensities. The effects of different noise intensities are also plotted in Figure 5 to show the effects of noise on the five models' overall trends. The general trends showed an association between the increase in noise intensity and a reduction in the determination coefficient's value as expected. High robustness is noticed by the low variance of R 2 of the triple seasonality of the Holt-Winters model with the least value of 0.9816 at 80% noise intensity. The Holt-Winters model's single seasonality showed the least robustness among the other models with a minimum value of 0.9552. The other models showed excellent robustness across the different noise intensities. These results drive to the same previous conclusion that the triple seasonality then double seasonality of the Holt-Winters model outperformed the other models in terms of accuracy, generalization, and robustness.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 19 The models' robustness is assessed by incorporating different Gaussian white noise intensities into the real tested data similar to relevant previous studies [44][45][46]. The data were divided into a training dataset and several testing sets with each set corresponding to a different noisy intensity. Four different intensities of 20%, 40%, 60%, and 80% were utilized to investigate the robustness of the five models. The models' robustness was evaluated through the accuracy of the coefficient of determination R 2 by studying the random distributed Gaussian white noise on the training dataset. Table 4 shows the variance of R 2 affected by different noise intensities. The effects of different noise intensities are also plotted in Figure 5 to show the effects of noise on the five models' overall trends. The general trends showed an association between the increase in noise intensity and a reduction in the determination coefficient's value as expected. High robustness is noticed by the low variance of R 2 of the triple seasonality of the Holt-Winters model with the least value of 0.9816 at 80% noise intensity. The Holt-Winters model's single seasonality showed the least robustness among the other models with a minimum value of 0.9552. The other models showed excellent robustness across the different noise intensities. These results drive to the same previous conclusion that the triple seasonality then double seasonality of the Holt-Winters model outperformed the other models in terms of accuracy, generalization, and robustness.

Future Generation Forecasting
Long-term electricity forecasting studies can include the peak load forecasting as well as electricity generation forecasting as in the current study [18]. Each of the two forecasting approaches contributes to providing different information for the experts and reveals a comprehensive overview of the electricity grid. One of the common misconceptions is that the generated electricity's annual peak occurs on the same day as the annual peak load. However, this may not be true for all datasets and case studies due to several factors. These factors can be attributed to the nature of the annual peak load date, resulting from the high consumer consumption at a specific time due to temperature, the time during the day, and consumption pattern. The annual peak generation is associated with these factors during the whole day and not at a specific time. As a result, the annual peak load dates are not identical with the annual peak generation dates. Table 5 provides a comparison between the peak load dates, from [18], and the peak generation of this study in Kuwait between 2015-2020. It is evident that forecasting the peak load and the peak generation would be advantageous to policymakers to prepare for the expected loads in terms of logistics needed for the electrical generation. Moreover, it would help in making decisions for expanding the electricity grid and to assess the need for boosting from the current interconnected GCC power grid. The long-term electric generation forecasting for the year 2020-2030 of the five models, and the actual data are tabulated in Table A1. The forecasted data of the annual electrical generation of the five models are plotted in Figure 6.

Future Generation Forecasting
Long-term electricity forecasting studies can include the peak load forecasting as well as electricity generation forecasting as in the current study [18]. Each of the two forecasting approaches contributes to providing different information for the experts and reveals a comprehensive overview of the electricity grid. One of the common misconceptions is that the generated electricity's annual peak occurs on the same day as the annual peak load. However, this may not be true for all datasets and case studies due to several factors. These factors can be attributed to the nature of the annual peak load date, resulting from the high consumer consumption at a specific time due to temperature, the time during the day, and consumption pattern. The annual peak generation is associated with these factors during the whole day and not at a specific time. As a result, the annual peak load dates are not identical with the annual peak generation dates. Table 5 provides a comparison between the peak load dates, from [18], and the peak generation of this study in Kuwait between 2015-2020. It is evident that forecasting the peak load and the peak generation would be advantageous to policymakers to prepare for the expected loads in terms of logistics needed for the electrical generation. Moreover, it would help in making decisions for expanding the electricity grid and to assess the need for boosting from the current interconnected GCC power grid. The long-term electric generation forecasting for the year 2020-2030 of the five models, and the actual data are tabulated in Table A1. The forecasted data of the annual electrical generation of the five models are plotted in Figure 6. Figure 6. The long-term annual electrical generation forecasting of the five models. Figure 6. The long-term annual electrical generation forecasting of the five models.
According to the forecast of the triple seasonality of the Holt-Winters model, the annual electric generation in Kuwait will reach about 92,535,555 kWh by the year 2030. This forecasted value is around 23% higher than the actual generation in 2019. The other models are higher in the forecasted annual electrical generation, which reaches up to 102,262,507 kWh in single seasonality of Holt-Winters models, which is higher by 36% than the generation of 2019. The management needs to assess the needed facility or equipment needed for such a generation to be established and prepared.
The forecasted monthly consumptions of the coming two years for different models are depicted in Table A2 and Figure 7. In addition to the study's primary aim, analyzing the monthly electrical generation can provide mid-term information for the logistics needed for this generation. The highest generation is clearly in two months, July and August, where Kuwait has its harshest weather temperature and humidity. Moreover, the higher electrical generation period of each year is between May and September. The maximum monthly electrical generation in 2021 is July for all the models, whereas August is the maximum in 2022 for the PSR, PMR, and HWSS models, and July 2022 is the maximum for the HWDS and HWTS models. Such information provides management and policymakers with vital tools to reach the right decision for the proper generation of electricity and planning for major maintenance or shutdowns.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 19 According to the forecast of the triple seasonality of the Holt-Winters model, the annual electric generation in Kuwait will reach about 92,535,555 kWh by the year 2030. This forecasted value is around 23% higher than the actual generation in 2019. The other models are higher in the forecasted annual electrical generation, which reaches up to 102,262,507 kWh in single seasonality of Holt-Winters models, which is higher by 36% than the generation of 2019. The management needs to assess the needed facility or equipment needed for such a generation to be established and prepared.
The forecasted monthly consumptions of the coming two years for different models are depicted in Table A2 and Figure 7. In addition to the study's primary aim, analyzing the monthly electrical generation can provide mid-term information for the logistics needed for this generation. The highest generation is clearly in two months, July and August, where Kuwait has its harshest weather temperature and humidity. Moreover, the higher electrical generation period of each year is between May and September. The maximum monthly electrical generation in 2021 is July for all the models, whereas August is the maximum in 2022 for the PSR, PMR, and HWSS models, and July 2022 is the maximum for the HWDS and HWTS models. Such information provides management and policymakers with vital tools to reach the right decision for the proper generation of electricity and planning for major maintenance or shutdowns. The monthly electrical generation introduces useful information to determine a map of electrical demand and associated maintenance program. In contrast, the long-term annual electrical generation forecasting help in estimating the fossil fuel power plants' fuel budget. Therefore, addressing fuel costs evaluation constitutes a positive addition in the current research and reflects the primary goal of electrical generation prediction for all-electric power entities.

Implications of Elecrical Generation on Fuels Cost
The long-term annual electrical generation forecasting can be used to estimate the fuel cost of the future till 2030. Since Kuwait is dependent on oil export as the primary source for the country's revenue, the electricity generation directly relates to the country's economy, especially with the depletion of oil due to burning fossil fuel as a fuel in electricity generation. Kuwait's power plants use four different types of fossil fuel: natural gas, gas oil, crude oil, and heavy fuel oil (HFO). The rate of increase in electricity consumption has significant effects on the stability of the country's economy. Therefore, the long-term forecasting of electrical generation is essential for policymakers to keep Kuwait's prosperous future. The sharp increase in electrical generation can deplete the oil The monthly electrical generation introduces useful information to determine a map of electrical demand and associated maintenance program. In contrast, the long-term annual electrical generation forecasting help in estimating the fossil fuel power plants' fuel budget. Therefore, addressing fuel costs evaluation constitutes a positive addition in the current research and reflects the primary goal of electrical generation prediction for all-electric power entities.

Implications of Elecrical Generation on Fuels Cost
The long-term annual electrical generation forecasting can be used to estimate the fuel cost of the future till 2030. Since Kuwait is dependent on oil export as the primary source for the country's revenue, the electricity generation directly relates to the country's economy, especially with the depletion of oil due to burning fossil fuel as a fuel in electricity generation. Kuwait's power plants use four different types of fossil fuel: natural gas, gas oil, crude oil, and heavy fuel oil (HFO). The rate of increase in electricity consumption has significant effects on the stability of the country's economy. Therefore, the long-term forecasting of electrical generation is essential for policymakers to keep Kuwait's prosperous future. The sharp increase in electrical generation can deplete the oil resources capacity and is critical to Kuwait's economic future. The total cost of fuel consumed in the electrical generation between 2012 and 2019 is tabulated in Table 6 to illustrate the fuel cost change. Table 6. The consumption of fuels and the total cost of fuel consumption between years 2012-2019 [21][22][23][24].

Year
Natural It is clear from Table 6 that the total cost of fuel is fluctuating and associated with world fuel prices. Assessment of the power and desalination plants' fuel future requirements is vital to secure power generation and water supply stability. However, this may require thorough investigations to address all the parameters and be examined in the future. Instead, the analysis of possible simplified estimation of the total cost of fuel associated with the long-term electrical generation forecasting is introduced to address the uncertainties of fuel prices. The fuel cost represents the first source of expenses based on the life cycle assessment for the whole plant lifespan, which exceeds the capital investment cost. Fossil fuel constitutes about 99.6% of the total energy supply sources to the Kuwait power sector [21][22][23][24], and the high percentage level of dependency is typical in most oil countries.
There are two natural gas sources in Kuwait, local and imported liquefied natural gas (LNG). However, each type of fossil fuel has its own chemical and physical characteristics, and the price fluctuates according to several factors such as source, quality of the fuel, and demand-supply bases. The percentage of fuel consumptions map are varied year per year in power generation sectors and mainly governed by the strategies of Kuwait Petroleum Corporation (KPC), which is oriented economically. The fixed baseline of fuel price was selected throughout this work, as illustrated in Table 7, based on the average values extracted from the past five years data to evaluate the fuel budget for the next ten years using the five models' extracted outcomes. Single and triple seasonality Holt-Winters models are the two models that cover the upper and lower band of the estimated cost of the forecasted annual electricity generation for Kuwait from 2020-2030 as it can be observed in Figure 8. Estimating the fuel budget is not limited only to the economic criteria but also to technical factors, such as type of unit, heat rate, degradation, and fuel quality. For the sake of simplicity, the average degradation value for the unit heat rate was assumed 1% as per registered and proposed by the manufacturers at standard conditions. Figure 8  All models elucidate consistent trend through the gradual increase in fuel cost due to increasing power generation. The gap between the upper model (HWSS) and lower model (HWTS) outcomes, as it can be noticed in Figure 8, increases at the end of the tested period compared to the beginning because of the rise in expected power generation difference, which exceeds 9 million kWh total annual production. The moderate growth of the fuel budget might be attributed to the rising dependency on natural gas on other liquid fossil fuels. Relying more on natural gas as a significant fuel source is evident in the electricity value between them.
The primary contribution of fuel to the total fossil fuel budget comes from natural gas with LNG, followed by Heavy fuel oil, gas oil, and crude oil for the year 2020, while gas oil and crude oil have the same contribution by 5% in the year 2030 as illustrated in Figure 9. The fuel distribution estimation was determined based on the actual fuel scenarios over the last five years, which do not depend on fuel price only but govern by fuel availability and type of available units. The distribution pattern shows that the natural gas share increases over the examined period by 10.9%, whereas the liquid fossil fuel reduces by 8.1%. The present study has confirmed the advantage of securing more natural gas on fuel future map by cost-effectively All models elucidate consistent trend through the gradual increase in fuel cost due to increasing power generation. The gap between the upper model (HWSS) and lower model (HWTS) outcomes, as it can be noticed in Figure 8, increases at the end of the tested period compared to the beginning because of the rise in expected power generation difference, which exceeds 9 million kWh total annual production. The moderate growth of the fuel budget might be attributed to the rising dependency on natural gas on other liquid fossil fuels. Relying more on natural gas as a significant fuel source is evident in the electricity value between them.
The primary contribution of fuel to the total fossil fuel budget comes from natural gas with LNG, followed by Heavy fuel oil, gas oil, and crude oil for the year 2020, while gas oil and crude oil have the same contribution by 5% in the year 2030 as illustrated in Figure 9. All models elucidate consistent trend through the gradual increase in fuel cost due to increasing power generation. The gap between the upper model (HWSS) and lower model (HWTS) outcomes, as it can be noticed in Figure 8, increases at the end of the tested period compared to the beginning because of the rise in expected power generation difference, which exceeds 9 million kWh total annual production. The moderate growth of the fuel budget might be attributed to the rising dependency on natural gas on other liquid fossil fuels. Relying more on natural gas as a significant fuel source is evident in the electricity value between them.
The primary contribution of fuel to the total fossil fuel budget comes from natural gas with LNG, followed by Heavy fuel oil, gas oil, and crude oil for the year 2020, while gas oil and crude oil have the same contribution by 5% in the year 2030 as illustrated in Figure 9. The fuel distribution estimation was determined based on the actual fuel scenarios over the last five years, which do not depend on fuel price only but govern by fuel availability and type of available units. The distribution pattern shows that the natural gas share increases over the examined period by 10.9%, whereas the liquid fossil fuel reduces by 8.1%. The present study has confirmed the advantage of securing more natural gas on fuel future map by cost-effectively The fuel distribution estimation was determined based on the actual fuel scenarios over the last five years, which do not depend on fuel price only but govern by fuel availability and type of available units. The distribution pattern shows that the natural gas share increases over the examined period by 10.9%, whereas the liquid fossil fuel reduces by 8.1%. The present study has confirmed the advantage of securing more natural gas on fuel future map by cost-effectively selecting new energy systems. The MEW's future fuel budget can be reduced by securing more local natural gas and implementing large-capacity renewable energy projects.

Conclusions
Long-term forecasting of electrical generation in Kuwait was investigated using five models to predict the electrical generation until 2030. Different statistical performance metrics were employed to assess the accuracy, generalization, and robustness of the models. The Holt-Winters model's triple seasonality outperformed the other models in terms of the MAPE, RMSE, CVRMSE, MAE, and R 2 . According to the different performance indicators, the double seasonality of the Holt-Winters model was the second-best performing model with superior accuracy. The Prophet model, with multiple regressors, has shown comparable performance to the Holt-Winters double seasonality model. In contrast, the least performing models were the single regressor Prophet model followed by the single seasonality Holt-Winters model. The estimated annual electrical generation for the year 2030 of the Holt-Winters models' triple and single seasonality models were 92,535,555 and 102,262,507 kWh, respectively. The generation of electricity in Kuwait during the upcoming ten years was then presented and discussed. The results reveal that the higher electrical generation period occurs between May and September of the upcoming two years. The five models' long-term electrical generation forecast was then used to estimate the fuel's total cost to provide an overview for policymakers to support the execution of planning decisions.