Prospect Prediction of Terminal Clean Power Consumption in China via LSSVM Algorithm Based on Improved Evolutionary Game Theory

In recent years, China’s terminal clean power replacement construction has experienced rapid development, and China’s installed photovoltaic and wind energy capacity has soared to become the highest in the world. Precise and effective prediction of the scale of terminal clean power replacement can not only help make reasonable adjustments to the proportion of clean power capacity, but also promote the reduction of carbon emissions and enhance environmental benefits. In order to predict the prospects of China’s terminal clean energy consumption, first of all, the main factors affecting the clean power of the terminal are screened by using the grey revelance theory. Then, an evolutionary game theory (EGT) optimized least squares support vector machine (LSSVM) machine intelligence algorithm and an adaptive differential evolution (ADE) algorithm are applied in the example analysis, and empirical analysis shows that this model has a strong generalization ability, and that the prediction result is better than other models. Finally, we use the EGT–ADE–LSSVM combined model to predict China’s terminal clean energy consumption from 2019 to 2030, which showed that the prospect of China’s terminal clean power consumption is close to forty thousand billion KWh.


Introduction
With the comprehensive deepening of China's power industry, the "Thirteenth Five-Year Plan" outline has addressed the major challenges of industrial transformation and put forward a plan for institutional reforms and upgrades in the next 5 years. The "13th Five-Year Plan" outline clearly proposes to promote revolution in four aspects, including energy consumption, supply, technology, and systems. Meanwhile, by 2020, the proportion of non-fossil energy in primary energy consumption must reach 15% [1]. The world's largest wind and photovoltaic market is in China and clean energy based on photovoltaic and wind power will continue to develop rapidly [2]. In 2019, the Chinese government actively developed photovoltaic and wind power-related construction without subsidy and parity in on-grid projects, fully implemented the competitive allocation mechanism for wind power and photovoltaic power plant projects, and established a new mechanism for renewable energy power consumption. All related policies promote the high-quality development of renewable energy sources. The clean energy power generation installations currently used in China mainly include wind energy, photovoltaic and biomass energy. According to "China Energy Development Report 2019" [3], from the perspective of power generation, in 2019, China's installed clean energy capacity reached 728 million kilowatts, with an increase of 12%. Among them, the installed hydropower capacity was 352 million kilowatts, the installed wind power From the perspective of terminal consumption, in 2019, China's total power consumption was 6.8 trillion KWh, of which wind power consumption was 395.7 billion KWh, and wind energy idle generation was 27.7 billion KWh, with a decrease of 14.2 billion KWh. The wind idle rate was 7%, with a decrease of 5%. At the same time, China's photovoltaic consumption was 183 billion KWh, and the idle amount of photovoltaic power was 5.49 billion KWh, which was 1.8 billion KWh lower than the previous year. China's hydropower consumption was 138.2 billion KWh, and the idle amount of water power was about 69.1 billion KWh. The national average water idle rate was 5%, with a decrease of 1.8%. China's biomass power consumption is 283.3 billion KWh. The proportion of terminal clean power consumption is shown in Figure 2. In order to accurately predict the scale of future terminal clean electricity consumption in China and provide planning and decision-making assistance for the Chinese government's energy strategy reform, this paper presents a combined least squares support vector machine (LSSVM) model based on evolutionary game theory and adaptive differential evolution. The algorithm is fully trained and From the perspective of terminal consumption, in 2019, China's total power consumption was 6.8 trillion KWh, of which wind power consumption was 395.7 billion KWh, and wind energy idle generation was 27.7 billion KWh, with a decrease of 14.2 billion KWh. The wind idle rate was 7%, with a decrease of 5%. At the same time, China's photovoltaic consumption was 183 billion KWh, and the idle amount of photovoltaic power was 5.49 billion KWh, which was 1.8 billion KWh lower than the previous year. China's hydropower consumption was 138.2 billion KWh, and the idle amount of water power was about 69.1 billion KWh. The national average water idle rate was 5%, with a decrease of 1.8%. China's biomass power consumption is 283.3 billion KWh. The proportion of terminal clean power consumption is shown in Figure 2.
Energies 2020, 13, x FOR PEER REVIEW 2 of 17 hydropower capacity was 352 million kilowatts, the installed wind power capacity was 184 million kilowatts, the installed capacity of photovoltaic power generation was 174 million kilowatts, and the installed capacity of biomass power generation was 17.81 million kilowatts. These values increased by 2.5%, 12.4%, 34%, and 20.7%, respectively, which can be shown in Figure 1. From the perspective of terminal consumption, in 2019, China's total power consumption was 6.8 trillion KWh, of which wind power consumption was 395.7 billion KWh, and wind energy idle generation was 27.7 billion KWh, with a decrease of 14.2 billion KWh. The wind idle rate was 7%, with a decrease of 5%. At the same time, China's photovoltaic consumption was 183 billion KWh, and the idle amount of photovoltaic power was 5.49 billion KWh, which was 1.8 billion KWh lower than the previous year. China's hydropower consumption was 138.2 billion KWh, and the idle amount of water power was about 69.1 billion KWh. The national average water idle rate was 5%, with a decrease of 1.8%. China's biomass power consumption is 283.3 billion KWh. The proportion of terminal clean power consumption is shown in Figure 2. In order to accurately predict the scale of future terminal clean electricity consumption in China and provide planning and decision-making assistance for the Chinese government's energy strategy reform, this paper presents a combined least squares support vector machine (LSSVM) model based on evolutionary game theory and adaptive differential evolution. The algorithm is fully trained and In order to accurately predict the scale of future terminal clean electricity consumption in China and provide planning and decision-making assistance for the Chinese government's energy strategy reform, this paper presents a combined least squares support vector machine (LSSVM) model based on evolutionary game theory and adaptive differential evolution. The algorithm is fully trained and has a good prediction accuracy performance. The rest of the article will be structured as follows: Section 2 reviews the literature from related fields of study in recent years to further clarify the Energies 2020, 13,2065 3 of 17 meaning and research method of the paper; Section 3 introduces the mathematical principles of LSSVM, adaptive differential evolution (ADE) and evolutionary game theory (EGT), and the overall algorithm programming process; in Section 4, the proposed EGT-ADE-LSSVM model is fully trained to reach the level of accurate prediction, and applied to China's terminal clean electricity through a comparison with the prediction results of a traditional backpropagation neural network (BP), LSSVM and ADE-LSSVM algorithms. An empirical analysis is carried out to prove that this paper proposes an intelligent model with a better prediction performance than other models. Finally, the model is used to predict China's terminal clean electricity consumption prospects for 2019-2030. Section 5 offers a conclusion.

Literature Review
Firstly, the prediction of terminal clean power consumption is a multi-dimensional complex system forecasting problem. At present, many scholars have used traditional prediction methods to study energy consumption and other issues, such as grey prediction theory [4,5], multiple regression [6][7][8][9], fuzzy judgment [10,11], etc. However, with the development of machine intelligence algorithms, the prediction accuracy of traditional prediction methods does not meet the needs of development [12]. Compared with traditional prediction methods, intelligent machine algorithms are more and more widely used in the prediction of energy planning. Among them, the LSSVM is the most frequently applied intelligent algorithm [13]. The LSSVM is an improved standard support vector machine (SVM) algorithm based on the least squares function and equality constraints [14]. Compared with standard support vector machines, it has a better prediction performance [15]. Wu et al. [16] used the Cloud-Based Evolutionary Algorithm (CBEA) to optimize LSSVM parameters, including penalty coefficient C and the radial basis kernel function parameter σ, and also verified that the model, after optimizing parameters, is better at predicting the power generation of wind turbines. Because the LSSVM uses equation constraints, Tan et al. [17] minimized the value. The multiplication system is defined as the objective function, thereby simplifying the prediction algorithm process and accelerating the prediction speed of wind turbines. De et al. [18] used the gray system (GM) model and SVM to realize the future of China's carbon emissions. It was found that the combined model works better. Ahmad et al. [19] analyzed in detail the application of the LSSVM and the packet data processing method in different prediction scenarios. Moreover, they proposed a combined model to predict the energy consumption of buildings with convincing results. De Giorgi et al. [20] used the LSSVM model to predict the photovoltaic power generation situation, and found that the combined prediction model has an optimized performance and successfully predicted 24h photovoltaic power generation.
Through our literature review, we found that, with different combinations of algorithms, a large improvement in the prediction effect of a single machine can be achieved. Liu et al. [21] combined the other algorithms to optimize the LSSVM algorithm. The prognostics health management (PHM) algorithm was used to predict the seasonal factor. Then, the combined PHM-LSSVM algorithm was used to model the equipment failure rate and good experimental results were obtained. Ying [22], in order to optimize the LSSVM, used the ant colony algorithm to realize the assessment of a college cyber security field. Karaboga et al. [23] studied the ant colony, bird swarm and other biological behavior mechanisms, and found that the artificial bee colony was surrounded by bee colony intelligence. They used the artificial bees group (ABG) algorithm to optimize multivariate functions. After comparison, it was found that the ABG algorithm has better results than traditional algorithms, such as the particle swarm optimization (PSO) algorithm. Zhao et al. [24] adopted the latest natural heuristic optimization method, namely the whale optimization algorithm (WOA) to search for the optimal values of two parameters of LSSVM; the algorithm can accurately predict future carbon emission trends, reflecting the applicational value of the combined algorithm. Yang et al. [25] described the basic principles of the PSO in detail, and various improvements to and applications of the PSO were introduced. A process for obtaining multiple prediction models using the PSO algorithm was proposed. Trelea et al. [26] analyzed the performance of the PSO algorithm. The application of the PSO algorithm to the optimization of LSSVM parameters laid an important foundation. Liu et al. [27] proposed an improved gravity search algorithm (AC-GSA) to optimize LSSVM parameters in order to artificially improve the GSA's performance. They used a new model to predict the thermal efficiency of a 600mw supercritical unit. The results show that the AC-GSA-LSSVM model is very effective for load forecasting. Gorjaei et al. [28] applied the LSSVM model to the prediction of two-phase flow in a wellhead packer and optimized two parameters of the LSSVM algorithm by using particle swarm optimization (PSO). The PSO-LSSVM model is very consistent with the actual measurement rate. The results show that the PSO-LSSVM model has good regression accuracy and generalization ability. Zhang [29] proposed a hybrid model that combines fuzzy clustering (FC), LSSVM, and the wolf pack algorithm (WPA), and tested the data on two cases. Test results show that the model has a high prediction accuracy and stability. Huan et al. [30] proposed a prediction model based on integrated empirical mode decomposition (EEMD) and LSSVM. The results show that the EEMD-LSSVM model is better than the wavelet denoising least squares support vector machine (WD-LSSVM) and the traditional LSSVM. Xue [31] optimized the LSSVM with an improved particle swarm optimization algorithm (IMPSO-LSSVM) and proposed a prediction model for the combined concrete compressive strength.
In recent years, evolutionary game theory (EGT) has been extensively applied in parameter optimization algorithms. In EGT, each participant is an anonymous and random game. Gintis [32] proposed EGT, which utilized the strategic behavior of game theory modeling based on a problem, and described the modeling process of evolutionary game theory in detail. Myerson [33] explained the core points of EGT, and put forward concerns about the process of building the model, laying the foundation for the application of later EGT models. The authors in [33,34] clearly defined the application direction of EGT. Zhang et al. [35] used the EGT model to model the problem of freeway accidents, and the number of accidental escapes was obtained. Therefore, EGT has been widely used in combination parameter selection, and evolutionary game theory is increasingly being introduced into intelligent algorithms [36]. Li et al. [37] used evolutionary stabilization strategies to mathematically coordinate the decision-making between private institutions and humanitarian organizations. The results showed that trust had a significant positive impact on promoting coordination. Selten [38] defined and classified asymmetric animal populations and modeled their evolution using large-scale data analysis. The evolution results of evolutionary games were obtained. Marzband et al. [39] applied evolutionary game theory modeling to integrated energy management systems, and considered a price bidding strategy based on the Nash equilibrium mechanism to clarify the supplier's response strategy after the game. Similarly, Xiao [40] proposed a new algorithm based on evolutionary game theory to solve the dynamic layout of virtual machines in energy centers and optimize the power consumption of energy centers. This modeling solution can reach theoretical optimality. Cai [41] applied a multi-person asymmetric evolutionary game model optimized by system dynamics to the governance of corporate environmental pollution problems. Simulation results showed that the gold factor is closely related to the limit of fines. Zhu et al. [42] modeled the relationship between Chinese local governments through evolutionary game theory and production enterprises, and analyzed their respective costs and benefits.

LSSVM
In order to solve the convex quadratic programming problem in the support vector machine algorithm, the square error term is introduced to transform it into a set of linear equations. This method transforms the inequality constraint problem in the classic support vector machine algorithm into an equality constraint problem, which is helpful in reducing the difficulty of solving [42].
In a set of training samples, D = (x k , y k ) x k ∈ R n , y k ∈ R, k = 1, 2, . . . , N , x k is the input vector of the n-dimensional integral energy of each system, y k is the output vector of the LSSVM model, n is the number of training samples.
The LSSVM model can be expressed in the feature space as follows: where ϕ(x) is a nonlinear mapping function, ω is a weight vector in a feature space and b is the amount of offset. According to the principle of structural risk minimization, the function estimation of LSSVM can be described as: where γ is the punishing parameters for the error, e k indicates the error variable, e k ∈ R.
The constraint can be expressed as follows: In order to extract the original spatial features of each system, this paper constructs a nonlinear mapping function in the original space to the highly dimensional space. According to (2), a Lagrangian equation can be defined as follows.
where α k is the Lagrange multiplier, α k ∈ R. According to KarushKuhnTucker (KKT) conditions, (4) makes the following optimizations: It can be known from Equation (5) that the value of each component α k is non-zero and is proportional to the error of the data sampling point. Therefore, the above problem can be transformed into a (n + 1) dimensional linear equation problem. By using the method of eliminating the weight vector w of the feature space and error variables e k , we get the matrix equation: According to the Mercer condition, there are mapping ϕ and kernel functions K(x).
We get α after calculating the combination of (6) and (7), and the regression function of the comprehensive prediction data of energy system is: Therefore, we get the output function of the LSSVM model: Energies 2020, 13, 2065 6 of 17

Adaptive Differential Evolution Algorithm
Adaptive differential evolution (ADE) is a heuristic random search algorithm that simulates the evolutionary differences of populations. The performance of the ADE algorithm is highly dependent on mutation strategies and controlling parameters, including population size NP, scaling factor H and cross constants CR. The ADE algorithm has four steps, including initialization, mutation operation, cross operation and selection operation.

Initialization
Let the population size be NP; if the dimension of the feasible solution space is D, then use X L to indicate L, the generation population.
Each individual is made up of D dimensional parameters, which can be expressed as:

Mutation Operation
Individuals in the parent population x L i can mutate to individuals v L+1 i through mutation operations. Two different groups were randomly selected in the population and their vectors were scaled to perform vector synthesis with the individual to be mutated: where H is the scaling factor, whose value range is from 0 to 1 [0,1]. x L r 1 is the randomly selected basis vector to be perturbed.

Cross Operation
The generated mutant individuals are crossed with the individuals of the original population to form new crossed individuals. The binomial crossed operation is: where i = 1, . . . NP; j = 1, . . . D. rand b is a random number from 0 to 1 in [0,1]. CR is the cross constant from 0 to 1 of [0,1].

Selection Operation
The selection operation follows the greedy genetic model of survival of the fittest, so individual offspring are always better than their parents or at least equal to their parents-that is, only when the fitness value of the individual of the new vector u i is greater than the fitness value of the target vector individual will the new vector individuall u i be accepted; otherwise, x i still remains in the next generation population, and continues to perform mutation and crossover operations as target vectors in the next iteration. The selection operation is as follows: where f (x) is the optimized objective function. Due to the fact that the population size range is fixed, the mutation factor and cross-factor are automatically adapted according to Equation (15) during the iterative process of evolution, and the optimal results are finally generated, as shown in Equations (16) and (17). During the continuous iteration of the program, the population tends to the optimal solution during evolution, the population favors the optimal solution during evolution, and the best mutation factors and cross-factors are allocated: where ∂ 1 and ∂ 2 are the probability of regulating CR i and H i , respectively. Because CR i and H i affect each other, the optimal values of these two parameters can be obtained after multiple simulation experiments. This article sets these two parameters as follows: where normrnd() is a normal random distribution datum with a mean of 0.5 and a deviation of 0.3.

Evolutionary Game Theory
The core idea of evolutionary game theory is that people are not completely rational, but rather that they establish a corresponding decision set through a trial and error mechanism. Therefore, the idea of evolutionary game theory is introduced in the adaptive differential evolution algorithm in this paper. The steps of evolutionary game theory are as follows.
Step 1: Set h(s, x ) as the utility function when the opponent uses strategy s, taking x as the strategy. Then, the stable conditions of evolution are as follows [41]: There are two populations, x and x , and their respective utility functions are: By playing against its own parent generation individuals, each population calculates its own average individual utility to determine its own fitness. The average fitness after the second game is calculated as follows: Then, the corresponding fitness of the individual H(x) can be calculated: where p i is the probability that the fitness value of an individual is occupied in the population.
Assume that half of the information comes from a population and the other half comes from the elite solution set. The probability of each individual is p i , so the probability that the generated children belong to the first population is as follows.
Energies 2020, 13, 2065 8 of 17 The flow of the EGT algorithm proposed in this paper is shown in Figure 3: Assume that half of the information comes from a population and the other half comes from the elite solution set. The probability of each individual is i p , so the probability that the generated children belong to the first population is as follows.
The flow of the EGT algorithm proposed in this paper is shown in Figure 3:

Combined EGT-ADE-LSSVM Model
In Section 3.3, it is clear that EGT can solve the shortcomings of the ADE algorithm well. In the ADE algorithm, disadvantaged individuals will be avoided after evolution, but these disadvantaged individuals will pass through N . After the second evolution, it is also possible to develop in the direction of the optimal solution. Meanwhile, the original superior individual may also proceed in the direction of the inferior solution, which is more in line with evolution in reality. Therefore, EGT can control the ADE algorithm to avoid developing in the direction of a local optimal solution. Therefore, the EGT-ADE algorithm can search for the penalty coefficient c and the radial basis kernel function global value  . Due to the fact that the initial parameter penalty coefficient c and the radial basis kernel function  of the traditional LSSVM model greatly affect the performance of the prediction model, this article uses the ADE algorithm improved by EGT to perform global search optimization on the parameters in order to increase the accuracy of the LSSVM algorithm.
During the operation of the combined EGT-ADE-LSSVM model, the optimal solution of parameters is always maintained to achieve accurate prediction of terminal clean power consumption in China. The algorithm principle of the clean power consumption prediction model in China based on EGT-ADE-LSSVM is divided into the following steps: (1) Collect experimental data and initialize all relevant parameters of the EGT-ADE-LSSVM model; (2) Operate the EGT-ADE algorithm to search for the penalty coefficient c and global value of radial basis kernel function  ; (3) Output the penalty coefficient c and the radial basis kernel function  corresponding to the global optimal fitness value into the LSSVM model; (4) Build training and testing sample sets of the LSSVM model and train multiple training sets; (5) Output the results of the testing sample sets.

Combined EGT-ADE-LSSVM Model
In Section 3.3, it is clear that EGT can solve the shortcomings of the ADE algorithm well. In the ADE algorithm, disadvantaged individuals will be avoided after evolution, but these disadvantaged individuals will pass through N. After the second evolution, it is also possible to develop in the direction of the optimal solution. Meanwhile, the original superior individual may also proceed in the direction of the inferior solution, which is more in line with evolution in reality. Therefore, EGT can control the ADE algorithm to avoid developing in the direction of a local optimal solution. Therefore, the EGT-ADE algorithm can search for the penalty coefficient c and the radial basis kernel function global value σ. Due to the fact that the initial parameter penalty coefficient c and the radial basis kernel function σ of the traditional LSSVM model greatly affect the performance of the prediction model, this article uses the ADE algorithm improved by EGT to perform global search optimization on the parameters in order to increase the accuracy of the LSSVM algorithm.
During the operation of the combined EGT-ADE-LSSVM model, the optimal solution of parameters is always maintained to achieve accurate prediction of terminal clean power consumption in China. The algorithm principle of the clean power consumption prediction model in China based on EGT-ADE-LSSVM is divided into the following steps: (1) Collect experimental data and initialize all relevant parameters of the EGT-ADE-LSSVM model; (2) Operate the EGT-ADE algorithm to search for the penalty coefficient c and global value of radial basis kernel function σ; (3) Output the penalty coefficient c and the radial basis kernel function σ corresponding to the global optimal fitness value into the LSSVM model; (4) Build training and testing sample sets of the LSSVM model and train multiple training sets; (5) Output the results of the testing sample sets.
The flow of the EGT-ADE-LSSVM algorithm is shown in Figure 4:

Influencing Factors Filter Model input
Using data from the "Energy Consumption" section of the China Statistical Yearbook and the statistical data published on the website of the China Energy Administration, this article selects gross domestic product(GDP), total population, urban population, coal price, GDP of primary industry, GDP of secondary industry, GDP of tertiary industry, net investment in energy industry, net investment in photovoltaic industry, net investment in wind energy industry, net investment in biomass energy industry, net installed wind power capacity, net installed photovoltaic capacity, fixed assets for biomass power generation, net installed capacity and net transmission and transformation investment and distribution network as the influencing factors for China's terminal clean energy consumption. Grey relevance theory is applied to eliminate some factors with a degree less than 0.4.

Grey Relevance Theory
The following are the relevant calculation steps for gray correlation theory: Step 1: Determine the relevant analysis series Determine the reference and the comparison series that reflect the behavior of the system series. The reference sequence refers to a data sequence that reflects the behavioral characteristics of the system. A series of data, as part of the factors that affect the behavior of the system, is called a

Influencing Factors Filter Model Input
Using data from the "Energy Consumption" section of the China Statistical Yearbook and the statistical data published on the website of the China Energy Administration, this article selects gross domestic product(GDP), total population, urban population, coal price, GDP of primary industry, GDP of secondary industry, GDP of tertiary industry, net investment in energy industry, net investment in photovoltaic industry, net investment in wind energy industry, net investment in biomass energy industry, net installed wind power capacity, net installed photovoltaic capacity, fixed assets for biomass power generation, net installed capacity and net transmission and transformation investment and distribution network as the influencing factors for China's terminal clean energy consumption. Grey relevance theory is applied to eliminate some factors with a degree less than 0.4.

Grey Relevance Theory
The following are the relevant calculation steps for gray correlation theory: Step 1: Determine the relevant analysis series Determine the reference and the comparison series that reflect the behavior of the system series. The reference sequence refers to a data sequence that reflects the behavioral characteristics of the system. A series of data, as part of the factors that affect the behavior of the system, is called a comparison series. Assume that the reference number is Y = {y(k)|k = 1, 2, . . . n }, and compare the series as follows: Step 2: Dimensionless variables Because the data of each factor column in the system may have different dimensions, it is inconvenient to compare these data, and it is difficult to draw the correct conclusion from this comparison. Therefore, in gray correlation theory, the data are usually dimensionless: Step 3: Calculate the correlation coefficient Calculate the correlation coefficient of x 0 (k) and x i (k): ρ ∈ (0, ∞) is the resolution coefficient. The smaller ρ, the greater the resolution. Under normal circumstances, the value range of ρ is (0, 1), the specific value depends on the situation. When ρ ≤ 0.5463, the resolution is the best, so we take ρ = 0.4.
Step 4: Calculate the correlation Because the concept of the correlation coefficient refers to the degree of correlation between the comparison sequence and the reference sequence at each moment (i.e., each point on the curve), the correlation coefficient has a lot of numbers, and too much information is not conducive to the overall comparison, so the correlation coefficient at each time (that is, each point on the curve) needs to be changed into a value, this value is then the average of the numbers, and is expressed as the number of correlations between the comparison sequence and the reference sequence. The correlation Equation r i is as follows: Step 5: Relevance ranking Sort the relevance based on size: if r 1 < r 2 , the reference sequence y is similar to the comparison sequence X 2 . The correlation coefficient between the X i (k) sequence and the Y(k) sequence are first calculated, and then the average of the correlation coefficients for each type is calculated, which is called the correlation degree between Y(k) and X i (k).
In this paper, grey relevance theory is used to select the main influencing factors that have a large correlation with the clean power consumption of the terminal. The set of influencing factors of terminal clean power consumption is shown in Table 1:

Data Normalization
In order to eliminate dimensional differences between data of different measures, all the corrected data are normalized. In the data set of M indicators, the Z-score data normalization method is used to correspondingly normalize the N groups of data.

Z mn =
x mn − x n σ , (m = 1, 2, . . . , j; n = 1, 2, . . . , j) Z mn is the normalized data, x mn is the correction data, x n is the mean of x mn and σ is the standard deviation of x mn .
For negative indices, this paper uses the inverse exponential positive processing method to obtain positive data: where X mn is the normalized data, max is the maximum value of the normalized data and min is the minimum value of the normalized data. Some details of the processed data are shown in Table 2.

Prediction Result Test of EGT-ADE-LSSVM Model
The data from 1990-2010 were input into the training set of the EGT-ADE-LSSVM model, and the data from 2011 to 2018 were used as the test set. The main parameters of the EGT-ADE-LSSVM model are shown in Table 3. In order to predict the accuracy of the model proposed in this paper, we simultaneously inputted sample data into three prediction models: a backpropagation (BP) neural network, LSSVM, and ADE-LSSVM, and used relative error (ER) to calculate the degree of prediction accuracy. The prediction results obtained by each prediction model are shown in Figure 5: whereŷ i is the predicted value and y i is the actual value.

Prediction Result Test of EGT-ADE-LSSVM Model
The data from 1990-2010 were input into the training set of the EGT-ADE-LSSVM model, and the data from 2011 to 2018 were used as the test set. The main parameters of the EGT-ADE-LSSVM model are shown in Table 3. In order to predict the accuracy of the model proposed in this paper, we simultaneously inputted sample data into three prediction models: a backpropagation (BP) neural network, LSSVM, and ADE-LSSVM, and used relative error (ER) to calculate the degree of prediction accuracy. The prediction results obtained by each prediction model are shown in Figure 5: where ˆi y is the predicted value and i y is the actual value.  The prediction results of different prediction models are compared as shown in Figure 6.
Energies 2020, 13, x FOR PEER REVIEW 13 of 17 The prediction results of different prediction models are compared as shown in Figure 6. As can be seen from Figures 5 and 6, this paper proposes a reliable EGT-ADE-LSSVM model that has a good prediction effect.
Through comparative analysis, we found that: (1) Parameter optimization is of great significance to machine intelligence algorithms, and prediction algorithms without parameter optimization have almost no effect on accurate predictions. Compared with the single model prediction, the fitting effect of the joint prediction model is obviously better. The parameter-optimized model has better prediction capabilities; (2) Different parameter optimization methods have different effects. The EGT-ADE-LSSVM model is a complex multi-combination prediction model that can achieve the complementary advantages of different algorithms; (3) Through the EGT-ADE algorithm continuously optimizing the parameters of the LSSVM algorithm, the model proposed in this article has more reasonable parameters than the ADE-LSSVM prediction algorithm.
Through the above analysis, the empirical analysis proves that the intelligent prediction model of the EGT-ADE-LSSVM method proposed in this paper is practical and effective in predicting the total consumption of China's terminal clean energy.

Model Output
We applied grey relevance theory to predict the main influencing factors from 2019 to 2030, which were applied as the input data for the EGT-ADE-LSSVM prediction model. Finally, we calculated China's total terminal clean energy consumption from 2019 to 2030, as shown in Figure 7:  As can be seen from Figures 5 and 6, this paper proposes a reliable EGT-ADE-LSSVM model that has a good prediction effect.
Through comparative analysis, we found that: (1) Parameter optimization is of great significance to machine intelligence algorithms, and prediction algorithms without parameter optimization have almost no effect on accurate predictions.
Compared with the single model prediction, the fitting effect of the joint prediction model is obviously better. The parameter-optimized model has better prediction capabilities; (2) Different parameter optimization methods have different effects. The EGT-ADE-LSSVM model is a complex multi-combination prediction model that can achieve the complementary advantages of different algorithms; (3) Through the EGT-ADE algorithm continuously optimizing the parameters of the LSSVM algorithm, the model proposed in this article has more reasonable parameters than the ADE-LSSVM prediction algorithm.
Through the above analysis, the empirical analysis proves that the intelligent prediction model of the EGT-ADE-LSSVM method proposed in this paper is practical and effective in predicting the total consumption of China's terminal clean energy.

Model Output
We applied grey relevance theory to predict the main influencing factors from 2019 to 2030, which were applied as the input data for the EGT-ADE-LSSVM prediction model. Finally, we calculated China's total terminal clean energy consumption from 2019 to 2030, as shown in Figure 7: Energies 2020, 13, x FOR PEER REVIEW 13 of 17 The prediction results of different prediction models are compared as shown in Figure 6. As can be seen from Figures 5 and 6, this paper proposes a reliable EGT-ADE-LSSVM model that has a good prediction effect.
Through comparative analysis, we found that: (1) Parameter optimization is of great significance to machine intelligence algorithms, and prediction algorithms without parameter optimization have almost no effect on accurate predictions. Compared with the single model prediction, the fitting effect of the joint prediction model is obviously better. The parameter-optimized model has better prediction capabilities; (2) Different parameter optimization methods have different effects. The EGT-ADE-LSSVM model is a complex multi-combination prediction model that can achieve the complementary advantages of different algorithms; (3) Through the EGT-ADE algorithm continuously optimizing the parameters of the LSSVM algorithm, the model proposed in this article has more reasonable parameters than the ADE-LSSVM prediction algorithm.
Through the above analysis, the empirical analysis proves that the intelligent prediction model of the EGT-ADE-LSSVM method proposed in this paper is practical and effective in predicting the total consumption of China's terminal clean energy.

Model Output
We applied grey relevance theory to predict the main influencing factors from 2019 to 2030, which were applied as the input data for the EGT-ADE-LSSVM prediction model. Finally, we calculated China's total terminal clean energy consumption from 2019 to 2030, as shown in Figure 7:   Table A2 in Appendix A section, combined with the terminal clean power consumption forecasted in this paper, which is shown in Table A1 in Appendix A section, the terminal clean power consumption trend of China in the next 10 years is finally obtained, as shown in Figure 8.
It can be seen from the trend in Figure 8 that the proportion of China's terminal clean power consumption has gradually increased, and it is estimated that, by 2030, this proportion will exceed 38%.
Energies 2020, 13, x FOR PEER REVIEW 14 of 17 According to the planned value of the entire power consumption of Chinese society proposed in the 13th Five-Year Plan of China, which is shown in Table A2 in Appendix A section, combined with the terminal clean power consumption forecasted in this paper, which is shown in Table A1 in Appendix A section, the terminal clean power consumption trend of China in the next 10 years is finally obtained, as shown in Figure 8.
It can be seen from the trend in Figure 8 that the proportion of China's terminal clean power consumption has gradually increased, and it is estimated that, by 2030, this proportion will exceed 38%.

Conclusions
This paper first defines terminal clean power consumption, and then proposes an innovative model to realize the extraction of influencing factors and forecast the prospect of China's terminal clean power consumption. The idea of evolutionary game theory (EGT) is applied to the differential evolution algorithm (ADE) to achieve the two parameters of LSSVM. We used the grey relevance theory to screen and obtain the set of main influencing factors that affect the terminal clean power consumption, and then normalized the original set. Then, we used the EGT theory to improve traditional ADE to form a new differential evolution algorithm and optimize the two parameters of the LSSVM. Using the data from 1995-2009 as the training set and the data from 2010-2018 as the test set, an innovative EGT-ADE-LSSVM forecasting model is used to predict China's terminal clean power consumption in 2010-2018. In order to verify the validity of this model, the same samples are inputted in the BP, LSSVM and ADE-LSSVM algorithms, etc. The prediction results are considerable, and verify that the EGT-ADE-LSSVM model has a high prediction accuracy and practicality. Finally, China's terminal clean power consumption and consumption proportion trend is predicted from 2019 to 2030 by using the prediction model and results proposed in this paper indicate that China's terminal clean energy consumption will reach 38.575 billion kWh in 2030 and the consumption proportion will reach 48.77%, indicating that China has great potential in the consumption of terminal clean power, and can implement the development goals of non-fossil energy production and consumption revolution.
Author Contributions: In this paper, S.Y. collected the relevant data and carried out the empirical analysis, results comparison and wrote the first manuscript. X.Z. and S.P. made improvements to the English and adjusted the article's structure. All authors have approved the submitted manuscript. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.

Conclusions
This paper first defines terminal clean power consumption, and then proposes an innovative model to realize the extraction of influencing factors and forecast the prospect of China's terminal clean power consumption. The idea of evolutionary game theory (EGT) is applied to the differential evolution algorithm (ADE) to achieve the two parameters of LSSVM. We used the grey relevance theory to screen and obtain the set of main influencing factors that affect the terminal clean power consumption, and then normalized the original set. Then, we used the EGT theory to improve traditional ADE to form a new differential evolution algorithm and optimize the two parameters of the LSSVM. Using the data from 1995-2009 as the training set and the data from 2010-2018 as the test set, an innovative EGT-ADE-LSSVM forecasting model is used to predict China's terminal clean power consumption in 2010-2018. In order to verify the validity of this model, the same samples are inputted in the BP, LSSVM and ADE-LSSVM algorithms, etc. The prediction results are considerable, and verify that the EGT-ADE-LSSVM model has a high prediction accuracy and practicality. Finally, China's terminal clean power consumption and consumption proportion trend is predicted from 2019 to 2030 by using the prediction model and results proposed in this paper indicate that China's terminal clean energy consumption will reach 38.575 billion kWh in 2030 and the consumption proportion will reach 48.77%, indicating that China has great potential in the consumption of terminal clean power, and can implement the development goals of non-fossil energy production and consumption revolution.
Author Contributions: In this paper, S.Y. collected the relevant data and carried out the empirical analysis, results comparison and wrote the first manuscript. X.Z. and S.P. made improvements to the English and adjusted the article's structure. All authors have read and agreed to the published version of the manuscript.