Scenario Analysis of Natural Gas Consumption in China Based on Wavelet Neural Network Optimized by Particle Swarm Optimization Algorithm

Natural gas consumption has increased with an average annual growth rate of about 10% between 2012 and 2017. Total natural gas consumption accounted for 6.4% of consumed primary energy resources in 2016, up from 5.4% in 2012, making China the world’s third-largest gas user. Therefore, accurately predicting natural gas consumption has become very important for market participants to organize indigenous production, foreign supply contracts and infrastructures in a better way. This paper first presents the main factors affecting China’s natural gas consumption, and then proposes a hybrid forecasting model by combining the particle swarm optimization algorithm and wavelet neural network (PSO-WNN). In PSO-WNN model, the initial weights and wavelet parameters are optimized using PSO algorithm and updated through a dynamic learning rate to improve the training speed, forecasting precision and reduce fluctuation of WNN. The experimental results show the superiority of the proposed model compared with ANN and WNN based models. Then, this study conducts the scenario analysis of the natural gas consumption from 2017 to 2025 in China based on three scenarios, namely low scenario, reference scenario and high scenario, and the results illustrate that the China’s natural gas consumption is going to be 342.70, 358.27, 366.42 million tce (“standard” tons coal equivalent) in 2020, and 407.01, 437.95, 461.38 million tce in 2025 under the low, reference and high scenarios, respectively. Finally, this paper provides some policy suggestions on natural gas exploration and development, infrastructure construction and technical innovations to promote a sustainable development of China’s natural gas industry.


Introduction
With fast growing urbanization and economic development, worldwide energy consumption has increased by 30% in the last 25 years [1].Compared with other fossil fuels such as coal and oil, natural gas is more efficient, clean and cheap, and thus many countries including China have been trying to increase the consumption of natural gas to improve the environment quality.In China, natural gas consumption has increased with an average annual growth rate of about 10% between 2012 and 2017, and the total natural gas consumption accounted for 6.4% of consumed primary energy resources in 2016, up from 5.4% in 2012, making China the world's third-largest gas user.Moreover, in the 13th Five-Year Plan of China, the development and natural gas industry has been adopted as an important task of the government.In such circumstances, forecasting natural gas consumption is a very important aspect of any country's energy policy and planning.Improper estimation of the natural gas consumption may subject final consumers to economic losses and lead to mismanagement of supplies and infrastructures [2].Therefore, accurately forecasting of natural gas consumption is of great significance for all the market participants to organize indigenous production, foreign supply contracts and infrastructures in a better way.
In the past few years, many efforts have been made, and many models have been proposed for different forecasting problems in the energy area.The proposed models can be generally classified into the following three categories: (1) statistical models; (2) artificial intelligence forecasting models; and (3) hybrid forecasting models.As for the statistical models, the widely used models mainly include auto-regressive moving average (ARMA), auto-regressive integrated moving average (ARIMA), grey forecasting model (GM) and generalized autoregressive conditional heteroscedasticity (GARCH).For example, Xu et al. established a new model with improved GM-ARMA based on HP Filter to forecast the final energy consumption in Guangdong Province, China [3].Sen et al. used ARIMA based models for forecasting energy consumption and greenhouse emissions, and obtained satisfied results [4].Zeng and Li presented a self-adapting intelligent grey model for predicting the natural gas demand in China during 2015-2020 [5].Wang and Wu forecasted the energy market volatility using both univariate and multivariate GARCH-class models [6].
The most often used artificial intelligence forecasting models include artificial neural network (ANN), extreme learning machine (ELM), support vector machine (SVM) and least squares support vector machine (LSSVM).For example, Wang et al. proposed a novel hybrid electricity price forecasting model by combining a two-layer decomposition technique and back propagation (BP) neural network, and the results show the efficiency of the proposed model [7].Geng et al. established a predictive modeling method based on the ELM integrated fuzzy C-Means integrating analytic hierarchy process (FAHP-ELM) for energy saving and forecasting problems [8].Ahmad provided a review on applications of ANN and SVM for building electrical energy consumption forecasting [9].Barman et al. proposed a hybrid forecasting model by integrating the grasshopper optimization algorithm and SVM model for the electricity load forecasting under local climatic conditions in Assam, India [10].Niu and Dai established a new short-term load predicting method by combining the empirical mode decomposition, gray relational analysis, and modified particle swarm optimization with the LSSVM model.The former two modules were applied to process the original load series, and the latter was used to predict the preprocessed subsequences [11].Among the above various forecasting methods, artificial neural networks (ANNs) have become the most used for energy related forecasting problems due to their attractive properties such as fault tolerance, strong ability in pattern recognition and distributed associative memory [12].Therefore, in the past several years, many different ANN-based forecasting models have been established for different forecasting problems [13,14].However, many studies have illustrated that ANN-based models trained using gradient-based approaches such as back propagation or its variants have some weaknesses when applied to various complex fields such as power load forecasting [12,14].
Based on the above considerations, to further improve the forecasting performance of the ANN-based models, many researchers have proposed various hybrid forecasting models based on ANNs.For example, Fard and Akbari-Zadeh [15] investigated the short-term load forecasting problem using a hybrid method based on wavelet, ANN and ARIMA model, and the experiments illustrated its good performance.Moreover, in standard ANN-based models, the initial weight values are usually determined randomly which may result in low robustness and convergence speed.To overcome the disadvantages of standard ANNs, some researchers have combined the optimization algorithms such as genetic algorithms (GA) and particle swarm optimization algorithms (PSO) with ANNs for constructing the forecasting models in which the optimization algorithm is used for training the connection weights of the network or optimizing design of the network structure [16].For instance, Muralitharan et al. proposed a novel neural network based optimization approach for energy demand prediction, and the proposed neural network includes the following two steps: (1) the conventional neural network approach is employed to find the required energy demand prediction at the consumer end; and (2) neural network based genetic algorithm (NNGA) and neural network based particle swarm optimization (NNPSO) approaches are implemented where the weights of the neural network are automatically adjusted.The experimental result reveals that the proposed NNGA approach performs better for short term load forecasting and proposed NNPSO is more suitable for the long term energy prediction [17].Ren et al. established a modified elman neural network (ENN) with a new learning rate scheme, and the experiments demonstrated the effectiveness of the proposed scheme from the aspects of convergence speed and consumption time with some popular schemes such as the original ENN, and PSO-ENN which uses PSO algorithm to search the best structure of ENN [18].
Based on aforementioned studies, the main novelties and contributions of this study can be denoted in the following two respects: (1) This study establishes a novel ANN-based hybrid model by combining the PSO algorithm, wavelet analysis and a modified learning rate scheme for the scenario analysis of natural gas consumption in China.Even through some of the previous studies have proposed some ANN-based hybrid models, such as ENN with a novel learning rate scheme [18], GA-ANN [19], differential evolution (DE) algorithm optimized WNN [20], and DE-ELM [21], and have proven the efficiency of the proposed models, the ANN-based hybrid model combining PSO algorithm and wavelet neural network with a modified learning rate scheme needs to be paid more attention.(2) This study conducts the scenario analysis of natural gas consumption from 2017 to 2025 in China based on three scenarios: low scenario, reference scenario and high scenario.
The reminder of the paper is organized as follows.Section 2 introduces the main methodology adopted in this paper.Section 3 introduces the combination strategies and proposes the PSO-WNN hybrid forecasting model.Section 4 first presents the main factors affecting China's natural gas consumption and then tests the performance of the proposed forecasting model.Section 5 first establishes three different scenarios according to three different development situations, and then analyzes the natural gas consumption under the three scenarios and provides some suggestions.Section 6 concludes this paper.

Methodology
The research methods used in this study include particle swarm optimization algorithm and wavelet neural network.Brief descriptions of those methods are stated as follows.

Particle Swarm Optimization (PSO)
Particle swarm optimization (PSO) is an efficient intelligent algorithm for solving the optimized problems [22].Similar to other swarm-based evolutionary computations, this algorithm also adopts the population concept and evolving iteration to achieve the optimization purpose.At the beginning of the algorithm, all the designed variables are collected to be a parameter vector which is called the particle or the individual from the viewpoint of PSO.Each particle also represents a candidate solution of the optimized problem, and a pre-specified number of particles further constitute a population.In the PSO algorithm, the particle movement is fully dominated by two particular particles: the individual best denoted by P best and the global best denoted by G best .The velocity updating formula of the algorithm utilizes these two particles' information to produce a new particle velocity.After that, next particle position is then derived by the position updating formula using the particle velocity.Based on the evolutionary mechanisms, all particles guided by the global best eventually converge to some system optimum of the optimized problem by successively executing a certain number of iterations.
In the PSO algorithm, let X = [x 1 , x 2 , ..., x n ] be a particle, where x i , i = 1, 2, ..., n, is the ith system variable of the optimized problem and n is the number of variables.To find the solution of the optimized problem, the following two formulas should be used to guide the particle movement in the iteration process of PSO algorithm. (1) where Equations ( 1) and ( 2) represent the velocity formula and position formula, respectively; x ij , pbest ij , and gbest j are the jth position component of the ith particle, the ith individual best, and the global best, respectively; v ij is the jth velocity component of the ith particle; ω is the inertia weight which balances the global and local search; c 1 and c 2 are two positive constants; r 1 and r 2 are two uniformly random numbers chosen from the interval [0,1]; and k is the iterative counter.
The pseudocode of PSO algorithm is listed in Algorithm 1.
Algorithm 1 Pseudocode of PSO algorithm.

Begin
Create the particle swarm and initialize particles repeat For all

End if
Update the particles using the following two equations

End for
Until the termination condition is satisfied Output the best particle position End

Wavelet Neural Network (WNN)
The wavelet theory, proposed by Morlet in 1980, appears to be more effective than Fourier transform in studying particularly non-stationary time series [23,24].Fourier transform gives frequency domain representation for temporal processes of signal or function that is only available in the time domain.Fourier analysis provides frequency information that can only be extracted for the complete duration of the signal, and it provides no information about the local variations in time and requires the signal to be stationary.Fourier transforms decompose stationary signals into linear combinations of sine and cosine waves, while the continuous wavelet transforms decompose non-stationary signals into linear combinations of the wavelets.
The continuous wavelet transform begins with: where W f (a, b) are the wavelet coefficients; a is scaling function used to stretch or compress mother wavelet ψ(t) that is related to the frequency of the signal; b is the translation function used to shift mother wavelet ψ(t) to a time domain of the signal; and f (t) is the input signal.The mother wavelet ψ(t) is defined as follows.
In this study, the Morlet function (Equation ( 5)) is taken as the mother wavelet in the analysis process.
In terms of the theory of wavelet and the learning ability of neural network, wavelet and the neural networks can be integrated into the wavelet neural network (WNN) to improve the approximation ability, especially on the catastrophe points, of the neural networks.In WNN, wavelet basis functions are used as node activation functions, WNN introduces the wavelet decomposition property into a general neural network, and combines the advantage of time-frequency location of the wavelet transform and self-learning capability of artificial neural networks.Thus, it is approximate and robust.
The WNN in this paper is designed as a three-layer structure with an input layer, a hidden layer, and an output layer.Each layer has one or more nodes.Figure 1 shows the schematic diagram of the proposed three-layer WNN.

Hybrid PSO-WNN Forecasting model
To avoid falling into the local optimum, the weight values of WNN and the parameters of wavelet, a and b, are optimized using PSO algorithm, see the Figure 2. Let X = (x 1 , ..., x m ) and Y = (y 1 , ..., y n ) be the input and output vectors of WNN, and l,W ik and W kj (i = 1, 2, ..., m, k = 1, 2, ..., l,j = 1, 2, ..., n) be the total number of hidden neurons, weight values of the input layer and output layer, respectively.Moreover, suppose that the hybrid forecasting model has a number of P studying samples which are denoted by I = [X 1 , ..., X P ], O = [Y 1 , ..., Y P ], where X p = x 1p , ..., x mp , Y p = y 1p , ..., y np , p (= 1, ..., P) represents the p th input sample.Let Z k = z k1 , ..., z kp be the output of the k th hidden neuron.In the PSO algorithm, the four parameters w ik, w kj and a k , b k are coded as one particle, and then the detailed steps of this hybrid PSO-WNN forecasting model can be described as follows.
Step 1: Initialization.Initialize the neural network structures including the number of layers, the number of nodes in each layer, the weight values of WNN, w ik, and w kj , and the parameters of wavelet, a and b.Initialize the parameters of PSO algorithm including the two positive constants c 1 and c 2 , the maximum particle velocity, and the total number of particles denoted by popsize.
Step 2: Population generation.Code the four neural network parameters including w ik, w kj and a k , b k as the particle position vector, i.e., where pos(i) represents the position of the ith particle.Then, randomly generate a number of popsize particles as an initial population denoted by pop(0) using both random and artificial ways.The position of each particle represents one feasible solution.
Step 3: Fitness calculation.Calculate the difference between network's output and real values and take this difference as the fitness function of PSO algorithm, i.e., then, calculate the fitness of each particle according to Equation (7), and take the smallest fitness value as the initial values of P best and G best .
Step 4: Update the P best .Compare the fitness of each particle ( f (x i )) with P best .If f (x i ) < P best , then replace P best by f (x i ), i.e., P best = f (x i ).
Step 5: Update the G best .Compare the fitness of each particle ( f Step 6: Update the position and velocity vectors of particles using the Equations (1) and (2).
Step 7: Stop the PSO algorithm if one of the following is reached: (1) the current iteration number has reached the maximum number of generations; or (2) the fitness value of the particles remains constant for 50 iterations.Then, output the best particle which includes the best combination of w ik, w kj and a k , b k , and go to Steps 8-10 for training the WNN using the training sample.Otherwise, return to Step 3.
Step 8: Calculate the output vector.Input the p th sample, obtain the output of the k th hidden neuron through the following Equation (8).
where ψ is the Morlet wavelet function, . For the p th sample, calculate the output of the j neuron in the output layer through the following way: where w j0 is the threshold value of the j neuron in the output layer.
Step 10: Train the network until a set of w ik, w kj and a k , b k that satisfies where ξ is the pre-specified error and O p = o 1p , ..., o np is the real value vector related to input sample X p , is found.In the training stage of WNN, a dynamic learning rate mechanism is adopted to update the parameter values to improve the convergence speed of WNN, and the steps of this dynamic learning rate mechanism is presented as follows.
Step 1: In the training process of the network, if the total error denoted by ER t+1 of (t + 1) th iteration is larger than that of t th iteration, and the difference between them is larger than the pre-specified value ς (specified as 3% in this study), then this updating process is ignored and modify the learning rate by η t+1 = (1 − α)η t .
Step 3: If ER t+1 > ER t , while the increasing rate is less than ς, update w ik, w kj and a k , b k and keep the current learning rate value.
Step 4: The above process of updating the learning rate can be summarized as follows.
According to the entropy function theory, for the minimization function, entropy function has faster convergence than mean square error.Thus, in the above updating process of the learning rate value, we set the error function ER as follows:

Performance Test of the Hybrid PSO-WNN Forecasting Model
In this section, the real data collected from 1995 to 2016 is adopted to test the performance of the proposed hybrid PSO-WNN forecasting model.Before the performance test process, we firstly analyze the factors affecting natural gas consumption in China.

Affecting Factors of Natural Gas Consumption in China
Natural gas consumption is generally affected by several factors such as GDP, gas price, economic structure, per-capita gas consumption, and proportion in the energy consumed [25].China's energy system is huge and complex with many uncertainties due to the driving forces of energy requirements.Thus, scientific studies on issues surrounding energy demand and consumption forecasting are challenging to conduct.Even more complicated are the factors affecting the energy consumption in China.The detailed description of key factors is presented as follows. 1.
Economic growth.Previous studies show that China's economic growth is the Granger cause of energy consumption growth.Energy demand and consumption tends to grow in line with GDP, although typically at a lower rate.Natural gas, which not only plays the role of important energy in the social development, but also an important raw chemical material in the industrial production, is the important original force of the economic development.Therefore, economic growth reflects the consumption of natural gas.However, due to the huge population of China, the GDP may not be related to an individual's income.Thus, in this study, per capita GDP is used to measure the economic growth.

2.
Total amount of gas production (TP).The energy production has close relation with consumption.Natural gas, which is a non-renewable resource, should be produced according to the demand.Therefore, natural gas production has a positive or negative effect on consumption.In recent years, the increasing proportion of natural gas in the energy production structure promotes the growth of consumption.

Household consumption level (HCL)
. With increasing of household consumption level, the public environmental awareness and quality of life improve greatly, which promote the wide use of natural gas in China.

4.
Population with access to gas (PAG).The population with access to gas is an important index to reflect the number of gas users.The population with access to gas increases with the development of gas infrastructure, which will drive more gas consumption.

5.
Urbanization ratio (UR).Urban and rural residents behave differently in terms of gas use and consumption level.Urban residents have access to better gas services because of better gas supply infrastructure such as urban gas pipelines.It is obviously that urban population growth stimulates the gas consumption.However, rural residents are less privileged in this aspect since they could not afford commercial gas consumption due to their income level.At current stage, China's rural residents contribute little to the total gas consumption.Therefore, urbanization ratio is considered as an important gas consumption factor.6.
Gas price (GP).According to accepted economic theory, price is the primary factor that affects consumption and the important lever for the balance of supply and demand.Price has a significant effect on consumer behavior.When the gas price is too high, users will use some other energy sources such as electricity or coal.In China, even though the gas pricing mechanisms has been changed from the previous cost-plus pricing method to the net market value pricing method, the current gas pricing process is still regulated by the government, and thus cannot fully reflect the scarcity of natural gas [26].Therefore, the gas price is not considered as an affecting factor in this study.

Correlation Analysis between the Affecting Factors and Natural Gas Consumption
Correlation analysis measures the relationship between two variables, and its resulting value called the "correlation coefficient" shows if changes in one variable (independent variable) will results in changes in the other one (dependent variable), which helps us understand an indicator's predictive abilities [27,28].The direction of the dependent variable's change depends on the sign of the coefficient.If the coefficient is a positive number, then the dependent variable will move in the same direction as the independent variable; if the coefficient is negative, then the dependent variable will move in the opposite direction of the independent variable.The correlation coefficient can range [−1, 1].A low correlation coefficient suggests that the relationship between two variables is weak or nonexistent.A high correlation coefficient indicates that the dependent variable will usually change when the independent variable changes.There are several correlation coefficients, and the most common of these is the Pearson correlation coefficient which is adopted in this paper.Table 1 shows the general decision range of two variables' correlation.According to the actual data for 1995-2016, which were obtained from China Statistical Yearbook [29], the Pearson correlation coefficient of affecting factors listed in Section 4.1 and natural gas consumption are calculated using the SPSS 19.0, as shown in Table 2. Table 2 indicates that there is a significant or strong correlation between natural gas consumption and affecting factors including per capita GDP (unit: RMB yuan), total amount of gas production (unit: million tce ("standard" tons coal equivalent)), household consumption level (unit: RMB yuan), population with access to gas (unit: ten thousand people), urbanization ratio (unit: percent), and, therefore, the combination of these five affecting factors has strong predictive ability for forecasting the natural gas consumption in China.

Models Comparison
The data of natural gas consumption and the above five affecting factors between 1995 and 2016 are collected from China Statistical Yearbook [29], as shown in Table 3.The per capita GDP, total amount of gas production, household consumption level, population with access to gas and the urbanization ratio are taken as the inputs of the forecasting models, and the amount of natural gas consumption is adopted as the output.programmed in MATLAB language.The data are divided into two sets, in which the first set including the first 18 years of data (from 1995 to 2012) is taken as the training sample, and the second set including the remaining data (from 2013 to 2016) is adopted as the test sample.This paper adopts the mean absolute percentage error (MAPE) to test the efficiency of the proposed hybrid PSO-WNN model.The computational formulas of MAPE is provided as follows: where P is the number of samples, and Y p and O p are the forecast and real values of natural gas consumption.The fitting values and its relative errors of ANN, WNN and PSO-WNN models are shown in Table 4.Moreover, to well present the difference between these three forecasting models, the performance of ANN, WNN and PSO-WNN models are shown in Figure 3.In Table 4 and Figure 3, it is obviously that the PSO-WNN model preforms better than ANN and WNN based models in both simulation and forecasting stages, which illustrates the good forecasting performance of PSO-WNN model.The reason that the PSO-WNN model outperform ANN and WNN based models can be interpreted by the following three aspects: (1) by combining the wavelet analysis and neural network, the WNN obtains strong function approximation ability, especially on the catastrophe points; (2) by adjusting the wavelet parameters and applying a dynamic learning rate mechanism for updating the connection weight values, WNN can effectively make up the disadvantage of falling into local optimum of traditional ANNs; and (3) by optimizing the connection weight values and wavelet parameters using PSO algorithm, the convergence efficiency and forecasting precision of WNN are effectively improved.

Scenario Analysis of Natural Gas Consumption in China during 2017-2025
Natural gas consumption is influenced by many factors which contain uncertainty that will be passed on to the forecast.To reduce this uncertainty, scenario analysis can be applied to investigate different outcomes depending on alterations in the underlying assumptions.Scenario analysis is a process of analyzing possible future events by considering alternative possible outcomes [30,31].Thus, scenario analysis, which is a main method of forecasting, does not try to show one exact picture of the future.Instead, it presents consciously several alternative future developments.
In this paper, three different scenarios are established and categorized as high, reference and low scenarios.Parameters such as per capita GDP, total amount of gas production, household consumption level, population with access to gas, and urbanization ratio are chosen to reflect the storylines of the different scenarios.The increasing rates of these five affecting factors and their corresponding values are shown in Tables 5 and 6.Based on the above three scenario settings, this study takes the five affecting factors and the natural gas consumption as the input and output of PSO-WNN forecasting model, respectively.Based on this proposed model, the natural gas consumption in China from 2017 to 2025 has been analyzed under low, reference and high scenarios, and the results are listed in Table 6 and Figure 4. From Figure 4, the following three observations can be obtained.

1.
Based on the results of scenario analysis, the China's natural gas consumption is going to be 342.70,358.27, 366.42 million tce ("standard" tons coal equivalent) in 2020, and 407.01, 437.95, 461.38 million tce in 2025 under the low, reference and high scenarios, respectively.2.
The natural gas consumption in the high, reference and low scenarios have a similar increasing trend, while, over time, the gap of natural gas consumption between high and low cases becomes larger and larger.

3.
In all three scenarios, natural gas consumption increases relatively rapid from 2017 to 2020, while, after 2020, it increases relatively slow and has the trend to be relatively stable after 2025, which may be interpreted as, in the first four years (2017-2020), the five affecting factors have direct and important influences on the natural gas consumption, while, afterwards, many other factors such as policy and international energy environment, which are not considered in this study, will have influences on the forecasting results.

Conclusions
Accurately forecasting natural gas consumption becomes very crucial for all market participants to organize indigenous production, foreign supply contracts and infrastructures in a better way.This paper first presents the main factors affecting China's natural gas consumption, and then establishes a hybrid forecasting model by combining the Particle Swarm Optimization algorithm and Wavelet Neural Network.Finally, this study conducts scenario analysis of natural gas consumption from 2017 to 2025 in China based on low, reference and high scenarios.Based on the results, the following conclusions and strategic enlightenments are obtained: 1.
The combination of five affecting factors, namely per capita GDP, total amount of gas production, household consumption level, population with access to gas and urbanization ratio, obtained using the correlation analysis, has significant or strong predictive ability for natural gas consumption in China.

2.
The PSO-WNN model successfully predicts gas consumption as reflected by a MAPE value of 2.32% for prediction, and outperforms others.Based on the experiment shown in Figure 4 The proposed hybrid PSO-WNN model has satisfactory performance for natural gas consumption forecasting, which can be adopted as an important tool for government or company energy related decision-making processes.

Figure 1 .
Figure 1.Basic structure of the WNN model.

Table 1 .
The general decision range of two variables' correlation.

Table 2 .
Correlation degree between the affecting factors and natural gas consumption.

Table 3 .
Consumption of China's natural gas and values of affecting factors between 1995 and 2016.

Table 4 .
Performance comparison of ANN, WNN and PSO-WNN models.

Table 5 .
The increasing rate of each affecting factor in Low, Reference and High scenarios.

Table 6 .
Results for each scenario under different settings of affecting factors.
, it is obviously that WNN model performs better than ANN model, which can be explained as, by combining the wavelet and neural work, WNN obtains strong function approximation ability, especially on the catastrophe points.Moreover, by adjusting the wavelet parameters and applying a dynamic learning rate mechanism for updating the connection weight values, WNN can effectively avoid falling into the local optimum.PSO-WNN model outperforms WNN model, which can be interpreted as the optimization of network weights and wavelet parameters using PSO algorithm effectively improves the forecasting precision and reduces fluctuation of WNN model.3. Natural gas consumption in China will keep a relatively rapid growing tendency.According to the results of scenario analysis, the China's natural gas consumption is going to be 342.70,358.27, 366.42 million tce ("standard" tons coal equivalent) in 2020, and 407.01, 437.95, 461.38 million tce in 2025 under the low, reference and high scenarios, respectively.To satisfy the increasing demand of natural gas consumption, the Chinese government should take some constructive measures: (a) The Chinese government should promote new exploration and development of natural gas reserves, especially in the southwest and northwest regions of China.The current estimation of natural gas resources indicates a URR of 22 trillion cubic meters, mainly distributed in southwest and northwest China; (b) The Chinese government should promote infrastructure construction such as enlargement of natural gas pipeline networks; (c) The Chinese government should accept that natural gas imports are unavoidable in future, thus seeking new suppliers and forging new relationships and collaborations in the gas sector are essential.