Forecasting Natural Gas Production and Consumption in United States-Evidence from SARIMA and SARIMAX Models

: Research on forecasting the seasonality and growth trend of natural gas (NG) production and consumption will help organize an analysis base for NG inspection and development, social issues, and allow industrials elements to operate effectively and reduce economic issues. In this situation, we handle a comparison structure on the application of different models in monthly NG production and consumption forecasting using the cross-correlation function and then analyze the association between exogenous variables. Moreover, the SARIMA-X model is tested for US monthly NG production and consumption prediction via the proposed method for the ﬁrst time in the literature review in this study. The performance of that model has been compared with SARIMA ( p , d , q ) * ( P , D , Q ) s . The results from RMSE and MAPE indicate that the superiority of the best model. By applying this method, the US monthly NG production and consumption is forecast until 2025. The success of the proposed method allows the use of seasonality patterns. If this seasonal approach continues, the United States’ NG production (16%) and consumption (24%) are expected to increase by 2025. The results of this study provide effective information for decision-makers on NG production and consumption to be credible and to determine energy planning and future sustainable energy policies.


Introduction
Global warming has a critical problem in countries all over the world, seriously threatening the development and survival of human beings especially in recent years. However, most countries have started to explore clean and low-greenhouse gas energy transitions and are beginning to decrease greenhouse gas emissions [1]. Likewise, most countries have improved a transition in the energy matrix to reduce greenhouse gas emissions without restricting economic growth [2] and environmental politics [3]. An energy outlook report has forecast that the global demand for natural gas will increase in the next 15 years [4]. The comparison of petroleum and coal, natural gas blaze absolutely, and remainder content were minimized. Therefore, it purposes low emission more than fossil fuels [5]. Moreover, both liquefied petroleum gas and natural gas could act as intermediate fuels to control environmental pollution [6].
Natural gas (NG) is the only non-renewable energy resource [7], and the most applied energy resource in certain societies [8]. Its dependency is increasing day by day in oil and gas developed countries [9][10][11][12]. However, according to the 2020 statistics energy review report by British Petroleum, natural gas production has increased by 132 (bcm) or 3.4%, and natural gas consumption has increased by 78 (bcm) or 2%, one of the strongest growth rates seen in 2018 [13].
Energy management is the most challenging issue as today's situation will greatly affect future outcomes. Managing energy sources is the most important duty of decisionmakers and lawmakers [14], along with policymakers and other stakeholders [11]. Moreover, energy management is the main concern in preparing energy policy for developing countries [15]. The most used energy management variables are natural gas, all of which are imported. Evaluating natural gas consumption such as annually, monthly, weekly, daily, and hourly, will lead to the management of plans such as transportation, distribution, and natural gas production [16,17], along with sampling frequency, such as yearly, daily, and hourly [18,19]. In daily studies, NG demand was predicted. A successful daily NG forecasting will lead to determining temperature, pipeline pressures, the viscosity of the gas [8,20], and exogenous variables [21].
Furthermore, advantages in evaluating NG consumption in the short term can include "Sales agreements and supply of gas" [22]. The layout of the reference pipelines and gas storage facility development empowers a more capable NG network optimization [23,24], and NG transmission network [25]. This ensures that supply and supply are in balance [8,18]. An advantageous short-term NG prediction estimation contributes to social contributions [1,25]. It grants industrial chain components to the event effectively with minimal economic loss [26,27]. The existing literature tending to the individual financial and commodity markets has significantly shifted through the recent global pandemic [28].
In recent years, different NG production and consumption prediction methods have been continuously applied. These forecasting methods can be divided into: time series methods and machine learning methods. Linear regression and time series-based models analyze the association between the 'output variable' and exogenous factors; commonly, the relationship is indicated linearly. This model is applied for NG forecasting [26]. Though linear regression models are surely applicable, they prove a weakness in indicating "nonlinear" relationships. It is insufficient to interpret via a non-linear and linear relationship such as NG consumption [29]. Time-series methods can be determined as natural relationships. They are widely used in natural gas consumption forecasts [10,[30][31][32][33][34]. Although machine learning models can be determined, their relationships are usually non-natural. They are widely applied in evaluating natural gas consumption [35][36][37]. Furthermore, NG production and consumption forecasts are necessary for the sustainable growth of all countries.
In this study, according to the available statistics report, the production and consumption for NG will increase in future years equal to the increasing growth of the United States. Studies of the production and consumption forecast in the review of the literature are performed yearly. Although, the aim of seasonal production and consumption for NG cannot be expected through annual studies. For this reason, this study focuses on the United States natural gas production and consumption forecasted monthly by focusing on the seasonal balance, natural gas supply/demand, peak shaving, residential, commercial, transportation, and industrial. For correlation between monthly production and consumption forecast, the more accurate SARIMAX (p, d, q) * (P, D, Q) s forecasting model was proposed as a new approach. This is the first used in the field of natural gas. To show the performance of the models, the basic SARIMA (p, d, q) * (P, D, Q) s model-a popular technique-is performed, and two methods have been compared with forecast success. Therefore, correlativity between the monthly NG production and consumption and seasonal exogenous variables was reliably predicted by SARIMAX (p, d, q) * (P, D, Q) s model to determine high forecast success. Finding the right model based on natural gas data is part of our contributions. Based on these exogenous factors, the superiority in order accuracy was achieved.

Material and Methods
Natural gas production and consumption of the United States (US) is forecast until 2025 on a monthly basis. For this study, a system process in the collection of data and dataset partitions, models building, comparative analysis of the models, prediction, and estimations steps are followed. The research flowchart is indicated in Figure 1.
Energies 2021, 14, x FOR PEER REVIEW 3 of 18 of our contributions. Based on these exogenous factors, the superiority in order accuracy was achieved.

Material and Methods
Natural gas production and consumption of the United States (US) is forecast until 2025 on a monthly basis. For this study, a system process in the collection of data and dataset partitions, models building, comparative analysis of the models, prediction, and estimations steps are followed. The research flowchart is indicated in Figure 1.

Data Measurement
United States natural gas data was taken from the Energy Information Administration. Then, we selected the variables-natural gas production and consumption (NGP and NGC). This information contains monthly natural gas and renewal energy between January 1983 to January 2021. The dataset is shown graphically in Figure 2. It indicates that the dataset is seasonal with a growth trend. The dataset was divided into training and test sets. The last 40 seasonal years in the training sets is used to organize an identical estimation of SARIMAX (p, d, q) * (P, D, Q) s models.

Data Measurement
United States natural gas data was taken from the Energy Information Administration. Then, we selected the variables-natural gas production and consumption (NGP and NGC). This information contains monthly natural gas and renewal energy between January 1983 to January 2021. The dataset is shown graphically in Figure 2. It indicates that the dataset is seasonal with a growth trend. The dataset was divided into training and test sets. The last 40 seasonal years in the training sets is used to organize an identical estimation of SARIMAX (p, d, q) * (P, D, Q) s models.  The Box-Jenkins model was used (Box et al. 1994). The ARIMA model is considered to be the future values of a linear function with past values and random errors [38]. One of the most common methods used in the time series analysis is the ARIMA model. The  The Box-Jenkins model was used (Box et al. 1994). The ARIMA model is considered to be the future values of a linear function with past values and random errors [38]. One of the most common methods used in the time series analysis is the ARIMA model. The ARIMA model consists of the terms autoregressive model (AR), integration (I), and moving average (MA) [39]. ARIMAX proposed the model by expanding the ARIMA model [40]. In this method, extra independent variables are applied in addition to the ARIMA method. The ARIMAX method can be specified as follows [41]: where ϕ(G) denotes the AR coefficients, γ(G) presents the MA coefficients, l i indicates the degree of lag, e t shows the error analysis in regression, a t Zero average with error terms time series.

SARIMA Models
The seasonal ARIMA (p, d, q) * (P, D, Q) s model is a development of ARIMA which was introduced to advance the performance of auto-regressive integrated moving average in modeling the seasonal series [42][43][44][45]. The general form [46,47] of the seasonal ARIMA method is expressed as given in the following Equation (2): Expression ϕ p (G) and γ q (G) indicate the order of the characteristic polynomial of non-seasonal autoregressive (AR) and non-seasonal moving average (MA) components.
The terms φ P G seasonal ω Q G seasonal and indicate the seasonal autoregressive (SAR) and seasonal moving average (SMA) polynomial. (1 − G) and 1 − G seasonal indicates the non-seasonal and seasonal time series are differencing components, respectively. The terms are d and D indicate the ordinary differenced terms of non-seasonal ARIMA, seasonal differenced terms of the seasonal ARIMA of the series [48,49]. In addition, X t indicates the observed value at the time t, e t stands for error terms of prediction, and s is the length of the seasonal pattern (e.g., s = 12, monthly series) and G is the backshift operator terms, respectively.
The SARIMA components can be written as follows:

Seasonal ARIMA with eXogenous Factors (SARIMAX)
Seasonal ARIMAX is an advancement of the seasonal ARIMA with external feature variables (X) called SARIMAX (p, d, q) * (P, D, Q), to improve its prediction and performance. The seasonal ARIMAX model can be suggested as follows [50]; where y k,t are the corresponding observations of the k th represented as the number of exogenous variables at the time t, and α k represents the correlation coefficient value of the k th exogenous (X) input variables. The seasonal ARIMA model with exogenous factors develops the capacity of the seasonal ARIMA model by integration of contextual information such as air temperature, humidity, average cold cover, wind speed, and other meteorological variables that affect natural gas demand and consumption [32]. Similarly, the indoor pollutant CO2 [51], mumps cases [52], and scarlet fever [53] were found in a study in the literature. In the current study, the SARIMAX method was organized to investigate natural gas production and consumption. The determining methods used to produce the model are given below.
Step 1: The data of time series is integration determines the order d(1,2) stationery.
Step 2: The selection of order AR and MA in the model parameters are determined for use in seasonal and non-seasonal ACF and PACF diagrams using determined values of (p, d, q) * (P, D, Q) s parameters estimated for the SARIMA and SARIMAX models. We assessed the relationship between pair sequences with robust autocorrelation. The cross-correlation function diagram was used to estimate the association between the NGP, NGC, and exogenous variables-also determined with covariates and lags selected for the best model. The Akaike information criterion and Schwarz Bayesian criterion values were used to assess the performance of the best model fit. Minimal AIC and BIC values indicate the best model fitting [54].

Analysis and Measures of Performance
To evaluate, the measure used to forecast the models in this study is the root mean square error (RMSE) and mean absolute percentage error (MAPE) [55][56][57], (Barnston 1992) shown in Equations (4) and (5), given below; where (E t ) and (F t ) represents the actual values and forecast values at t the time respectively, n is the estimate values. Table 1 presents summary statistics of natural gas production and consumption between January 1983 to January 2021. The US have maximally natural gas consumption at 3417.27, therefore, its minimum natural gas production is 3232.45. The maximal values for the natural gas market production and NGPL production are 3232.45 and 240.70, respectively. Similarly, residential, commercial, industrial, transport, and total natural gas consumption is 1037.19, 571.74, 952.25, 109.53, and 3417.27, respectively. Moreover, the minimal natural gas market production and NGPL production are 1284.80 and 59.96, respectively. Similarly, residential, commercial, industrial, transport, and total natural gas consumption is 99.78, 88.74, 435.47, 36.30, and 939.93, respectively. Looking at the standard deviation, natural gas consumption is more than natural gas production. Superior and inferior kurtosis can be found in natural gas production variables and natural gas consumption variables, which ensures the presence of a fat tail distribution. The natural gas production and consumption are almost positively skewed as the small p-values of the JB (Jarque-Bera test) reject the normality of the null hypothesis.  ) for the monthly dataset used in this studies for the natural gas production sector and consumption sector defined in this work. All the series in the table are in actual terms: natural gas production and consumption. JB refers to the statistic of the Jarque-Bera (JB) test normality (p > 0.05 presented the normality of series). Asterisk ***,** represents statistical significance at the 1% and 5% level, accordingly.

The Findings of the ADF, KPSS, and PP Unit Root Tests
The preliminary steps needed for time series methods is the inquiry of the integration of order into all variables. In the present study, we have used unit root tests, which are [58], ADF [59,60], and KPSS [61], and PP test [62,63] to extract information concerning the integration of order, and the unit root results are shown in Table 2. The unit root tests results suggest that all variables are non-stationarity by level form; however, we firstorder differenced [i.e., I(1)] them to make them stationary with Market Production, NGPL, Total NGC, Residential, Commercial, Industrial, and Transport been series are stationary statistical significance at 1%, 5%, and 10% levels. All stationary as represented by the ADF, KPSS, and PP test results.

Correlation Coefficient Matrix
The correlation coefficient matrix is shown in Table 3. There is a statistically significant positive correlation between Natural gas production and NG consumption. The correlation coefficient between market production and NGPL production in the United States is 0.979771; in the commercial sector and residential sector with superior consumption growth, it is 0.985891, and total Natural gas consumption and transportation sector is 0.931583, the results are presented in Table 3.

Application of SARIMA and SARIMAX Models on Natural Gas Production and Consumption
The proposed SARIMA and SARIMAX (X is exogenous variables) time-series methods are advanced for the monthly series of natural gas production and consumption with other exogenous factors (production sector and consumption sector). The ADF [59,60], KPSS [61], and PP [62,63] unit root tests for the series natural gas production and consumption data are tests for stationary (reported in Table 2). In the series of the monthly data after the first order integration process and seasonal integration, values of both D and d are 1. The ACF and PACF graphic of the natural gas production sector and consumption sector are presented in Figure 3. The autocorrelation function of lag values exceeds the critical values in the natural gas production and consumption sector. In the ACF of lag, values such as natural gas production and consumption sectors are cross-effective of the non-seasonal and seasonal autocorrelation. Moreover, the appropriate maximal values of the non-seasonal criteria (Q) and the seasonal criteria (q) are 1 and 2. Furthermore, the sample PACF values of lag that are significant are natural gas production and consumption sectors (Figure 3), so the appropriate maximal values of the non-seasonal criteria (p) and the seasonal criteria (P) are 2 and 4. Table 4 presents the best fit model based on significant parameters and coefficient statistics. The p-values are always p > 0.05, this model's criteria estimation appears statistically significant in the appropriate structure of seasonal ARIMA. The bestfitted model is NG production SARIMA (4,1,1) * (1,0,1) 12 and consumption SARIMA (2,1,1) * (5,0,0) 12 . Accordingly, the selected lowest values of AIC and BIC are presented in Table 4. The sufficiency of the seasonal ARIMA (4,1,1) * (1,0,1) 12 and SARIMA (2,1,1) * (5,0,0) 12 model is verified by the LB (Q) test statistics (0.01) for which the p-values (Production 0.90, Consumption 0.91) are presented in Table 4. The Ljung-Box test exceeds the p-value (0.05), but autocorrelation in the SARIMA model residual rejects the null hypothesis.  and d are 1. The ACF and PACF graphic of the natural gas production sector and consumption sector are presented in Figure 3. The autocorrelation function of lag values exceeds the critical values in the natural gas production and consumption sector. In the ACF of lag, values such as natural gas production and consumption sectors are cross-effective of the non-seasonal and seasonal autocorrelation. Moreover, the appropriate maximal values of the non-seasonal criteria (Q) and the seasonal criteria (q) are 1 and 2. Furthermore, the sample PACF values of lag that are significant are natural gas production and consumption sectors (Figure 3), so the appropriate maximal values of the non-seasonal criteria (p) and the seasonal criteria (P) are 2 and 4. Table 4 presents the best fit model based on significant parameters and coefficient statistics. The p-values are always p > 0.05, this model's criteria estimation appears statistically significant in the appropriate structure of seasonal ARIMA. The best-fitted model is NG production SARIMA (4,1,1) * (1,0,1) 12 and consumption SARIMA (2,1,1) * (5,0,0) 12 . Accordingly, the selected lowest values of AIC and BIC are presented in Table 4. The sufficiency of the seasonal ARIMA (4,1,1) * (1,0,1) 12 and SARIMA (2,1,1) * (5,0,0) 12 model is verified by the LB (Q) test statistics (0.01) for which the p-values (Production 0.90, Consumption 0.91) are presented in Table 4. The Ljung-Box test exceeds the p-value (0.05), but autocorrelation in the SARIMA model residual rejects the null hypothesis.  The cross-correlation function (CCF) is used to examine relationships between natural gas production and consumption with other exogenous factors NPGL production, residential, commercial, industrial, and transportation sectors. The CCF between the exogenous variables are monthly NPGL production, residential, commercial, industrial, and transportation is presented in Figure 4. Accepting only positive lags, all of the exogenous variables (NPGL production, residential, commercial, industrial, and transportation) are significantly associated with natural gas production and consumption of at least some of the lags. Similarly, Table 5 presents the best fit model based on significant parameters and coefficient statistics. The p-values are always p > 0.05, and this model's criteria estimation appears statistically significant in the appropriate structure of SARIMAX (X-exogenous variable). The best-fitted model is NG production SARIMAX (2,1,0) * (1,0,2) 12 and consumption SARIMAX (1,1,0) * (1,0,1) 12 . Accordingly, the selected lowest values of AIC and BIC are presented in Table 5. These results indicate that NPGL production, residential, commercial, industrial, and transportation sectors at integration I(1) affected natural gas production and consumption at integration I(1) after fitted the time series model. Natural gas production and consumption sectors at integration I(1) (coefficient =20.8187, 0.7572, 1.0017, 1.0670, 7.8810, p < 0.05) show a positive relationship with production and consumption, which was the best fitted SARIMAX model. The sufficiency of the SARIMAX (2,1,0) * (1,0,2) 12 and SARIMAX (1,1,0) * (1,0,1) 12 model is verified by the Ljung-Box (Q) test statistics (0.06, 0.04) for which the p-values (production 0.81, consumption 0.83) are presented in Table 5. The Ljung-Box test exceeds the p-value (0.05), but autocorrelation in the SARIMAX model standardized residual rejects the null hypothesis.     Authours compilation; Significance level codes; p < 0.001 ***, p < 0.01 ** indicate that the null hypothesis is rejected at 1%, and 5% levels, respectively. SE = standard error; AIC = Akaike information criterion; BIC = Bayesian information criterion.
The results show that diagnostic tests on the SARIMAX (2,1,0) * (1,0,2) 12 , SARIMAX (1,1,0) * (1,0,1) 12 (right side), SARIMA (4,1,1) * (1,0,1) 12 , and SARIMA (2,1,1) * (5,0,0) 12 (left side) are presented in Figures 5 and 6. The correlogram and standardized residual tests (Cryer and Chan 2008) prove that the residual is stationary time series. The Kernel density estimation (KDE) line of the Q-Q plots and histogram plots indicate that the line of the residual follows a normal distribution (Stoffer and Dhumway 2010). These results indicate that the diagnostic statistics test satisfies the inference of SARIMAX and SARIMA modeling of natural gas production and consumption.     The results were obtained for both sets after using SARIMA and SARIMAX. The received results are presented in Table 6. Furthermore, Figure 7 indicates the production and consumption forecasting performance of these models as a diaram for the test set. Furthermore, although SARIMA models show the best performance than SARIMAX in the training set data, this SARIMAX model shows better performance for the test data set in natural gas production and consumption. In addition, the best performance of SARIMAX is nearby the SARIMA model for the training set. Natural gas production and consumption on the SARIMAX model performance is 'best' and 'well' according to RMSE and MAPE measures for the training data set, respectively. The SARIMA models have good performance according to RMSE and MAPE measures for the test data set. The results were obtained for both sets after using SARIMA and SARIMAX. The received results are presented in Table 6. Furthermore, Figure 7 indicates the production and consumption forecasting performance of these models as a diaram for the test set. Furthermore, although SARIMA models show the best performance than SARIMAX in the training set data, this SARIMAX model shows better performance for the test data set in natural gas production and consumption. In addition, the best performance of SARI-MAX is nearby the SARIMA model for the training set. Natural gas production and consumption on the SARIMAX model performance is 'best' and 'well' according to RMSE and MAPE measures for the training data set, respectively. The SARIMA models have good performance according to RMSE and MAPE measures for the test data set  The SARIMA has not been as successful as the SARIMAX in forecasting the monthly NG production and consumption with the seasonal and trends. SARIMA (p, d, q) * (P, D, Q) s has six parameters. The various values of previous periods affect the forecast progress. The suggested approach is better evaluation of the SARIMA model. This main aim was identical to reach the SARIMA with exogenous variables.
To examine the forecast performance of the SARIMAX model, another comparison is made with four similar studies from the literature review. These first study authors [64] have achieved a prediction of the best performance with around 8% MAPE measure. In Figure 7. Forecast visualization of natural gas production and consumption with SARIMAX and SARIMA models with selected exogenous factors for the test set. The SARIMA has not been as successful as the SARIMAX in forecasting the monthly NG production and consumption with the seasonal and trends. SARIMA (p, d, q) * (P, D, Q) s has six parameters. The various values of previous periods affect the forecast progress. The suggested approach is better evaluation of the SARIMA model. This main aim was identical to reach the SARIMA with exogenous variables.
To examine the forecast performance of the SARIMAX model, another comparison is made with four similar studies from the literature review. These first study authors [64] have achieved a prediction of the best performance with around 8% MAPE measure. In that research, the monthly prediction values for only one year were tested data. Our studies are natural gas forecasting for the long-term in the test data set, and a more robust model was obtained. The second study was [65]. In this study, the researchers achieved a prediction best performance of around 5% MAPE measure. Furthermore, our studies are electricity demand for the long-term data were applied in the modeling set. The third study was [10]. In this study, the researchers achieved a prediction best performance with around Energies 2021, 14, 6021 13 of 17 9% MAPE measure. However, our studies are natural gas demand for the long-term data were applied in the tested set. The fourth study was [34]. In this study, the researchers achieved a prediction of best performance with around 0.357% MAPE measure. However, Long-term data were used in the training set. The comparison of the results is presented in Table 7. These results indicated that the SARIMAX is successful on seasonal prediction in a few points.  These best-proposed methods have been determined, while the result shows those compared with both models and from similar literature reviews. However, the United States' monthly natural gas production and consumption until 2025 is forecast by the SARIMAX model as covariates. The prediction values of natural gas production and consumption are presented in Table 8 and indicated in Figure 8. while the results were analyzed, it was evaluated that natural gas production and consumption will be increasing seasonally. The production and consumption increase was a regular condition in light of the industrial, residential, commercial, transportation, population growth, and economic developments in the current years. When the maximal production was conducted in December, the minimum production is in February. Similarly, when the maximal consumption was conducted in January, the minimum is in May. The yearly natural gas production and consumption in 2025 is expected to grow by 16% and 24% compared to 2020 if the seasonal trend continues. In a country that imports all of its natural gas, it is needful to modify the situation and apply different and mostly local energy sources such as nuclear power, coal, and renewable energy.  -3046  2754  2334  2175  2232  2554  2475  2267  2402  2510  2022  3279  3059  2795  2414  2270  2320  2610  2540  2351  2473  2572  2023  3272  3072  2833  2487  2355  2399  2661  2598  2427  2538  2629  2024  3266  3083  2867  2552  2432  2471  2707  2650  2496  2596  2681  2025  3561  3094  2898  3612  2501  2535  2749  3698  2558  2649  2727 Authors own compilation. Figure 8. Graph of the United States monthly natural gas production and consumption until 2025.

Conclusions and Policy Implications
The future prediction of energy, which is a critical component, is of signif developed and developing countries. Moreover, it will be feasible to determine ergy policy, plan energy investments, and take necessary precautions. In the s there are many benefits to evaluating natural gas production and consumption

Conclusions and Policy Implications
The future prediction of energy, which is a critical component, is of significance for developed and developing countries. Moreover, it will be feasible to determine future energy policy, plan energy investments, and take necessary precautions. In the short term, there are many benefits to evaluating natural gas production and consumption. In recent years, natural gas evaluations have been made using various methodologies. The most successful forecasting in monthly wise predictions studies can be obtained by improving the seasonal prediction methods. In this study, a new approach which is called the SARIMAX prediction method is recommended for forecasting the United States' monthly wise natural gas production and consumption.
The proposed methods focus on finding the most relevant past observation values by making identical estimations. Moreover, it performs tests for all applicable values on the seasonal coefficient. In this way, seasonality performance in the dataset is more successfully modeled. The success of the model has been detected using some comparisons. Firstly, this method is compared with SARIMA (p, d, q) * (P, D, Q) s and SARIMAX (p, d, q) * (P, D, Q) s along with the same dataset. The forecast accuracy measures as RMSE and MAPE results indicate the proposed models were higher than other models. The natural gas production and consumption time series dataset of this study were modeled and improved the forecast performance of the SARIMAX model with RMSE and MAPE values. Moreover, the conquest of the proposed method has been determined by comparing with similar articles related to monthly wise natural gas production and consumption. Finally, by use of the proposed method, monthly wise natural gas production and consumption for the United States was predicted until 2025.
The result shows seasonal fluctuation in the year will be determined with the monthly forecast model, it is feasible to analyze such a situation around seasonal balance, shortage of gas supply, and peak for decision-makers. This study organizes important outputs to the country's energy planning, decision-makers, and energy policies.
Suggestions for future policymakers are given below: (a) Establish a useful natural predict model that will play a key role in determining future gas purchase agreements, liquefied natural gas (LNG) and natural gas storage applications, consumption, production, and storage planning of renewable energy and other fossil energy sources, and finally, the country's energy policies; (b) improving on-demand flexibility and performance; (c) sustaining and enhancing infrastructure; (d) supporting energy sectors innovation; (e) accelerating the uses of sustainable energy, utilization, and carbon capture technologies; (f) invest in a long-term perspective and safe areas in nuclear energy; (g) policymakers need methods to measure and evaluate the current and future effects of energy use on human society, human