Carbon Reduction Potential of Resource-Dependent Regions Based on Simulated Annealing Programming Algorithm

In recent years, developing countries, especially resource-dependent regions, have been facing the paradox of ensuring both emissions reduction and economic development. Thus, there is a strong political desire to forecast carbon emissions reduction potential and the best way to achieve it. This study constructs a methodology to assess carbon reduction potential in a resource-dependent region. The Simulated Annealing Programming algorithm and the Genetic algorithm were introduced to create a prediction model and an optimized regional carbon intensity model, respectively. Shanxi Province in China, a typical resource-dependent area, is selected for the empirical study. Regional statistical data are collected from 1990 to 2015. The results show that the carbon intensity of Shanxi Province could drop 18.78% by 2020. This potential exceeds the 18% expectation of the Chinese Government in its ‘13th Five-Year Work Plan’ for Controlling Greenhouse Gas Emissions. Moreover, the carbon intensity of the province could be further reduced by 0.97 t per 10,000 yuan GDP. The study suggests that the carbon emissions of a resource-dependent region can be reduced in the following ways; promoting economic restructuring, upgrading coal supply-side reform, perfecting the self-regulation of coal prices, accelerating the technical innovation of the coal industry, and establishing a flexible mechanism for reducing emissions.


Introduction
The growing concentrations of greenhouse gases (GHGs), especially in resource-dependent regions, has been one of the major global challenges as it brings about environmental degradation and natural disasters threatening human safety and health [1].To reduce its influence and avoid more dangerous long-term effects, the Intergovernmental Panel on Climate Change (IPCC) called for limiting the increase in average global temperature to no more than 2 • C. To this end, we may need to see a reduction in carbon dioxide (CO 2 ) by at least 50% until 2050, which means that any future emissions leeway could become extremely constrained according to Pan et al. [2].Consequently, there is a strong political desire to forecast carbon emissions reduction potential and the best means of achieving it [3].Accordingly, we need to identify the factors influencing carbon reduction, understand in depth the mechanism of their effects, create a series of predictive models, and optimize carbon reduction potential.
The existing research on carbon reduction potential mainly covers the following; carbon emissions measurements [4][5][6][7], carbon emissions impact factors and mechanism analysis [8][9][10][11][12][13][14], scenario analysis, and unrelated.However, this assumption is impractical.Wang and Yang [39] and Zhang et al. [40] study the factors affecting carbon emissions in different levels in China.The results demonstrate that the driving factors of carbon emissions affected one another.To separate the interactions, a partial least squares regression and a principal component analysis are used in some studies [41,42].However, the clustering process brings new problems.The loss of the original information can be inconvenient to subsequent policy development.Therefore, a certain methodology system needs to be constructed that not only takes into account the interactions among factors but also retains the information of the original indicators.
This study contributes to the literature in the following ways: (1) The Simulated Annealing Programming (SAP) algorithm is introduced into carbon reduction research to assess carbon reduction potential.As a result, a new methodological system is constructed that covers gaps in existing studies.
(2) In the case example of Shanxi province in China, this study analyzes the carbon reduction potential of a resource-dependent region.The results provide evidence for policy makers in such regions.
(3) This study predicts the feasibility of the carbon reduction goals in the Chinese '13th Five-Year Work Plan' and, further, quantitatively calculates the potential of carbon reduction.
The remainder of this paper is organized as follows.Section 2 introduces the prediction method for carbon intensity based on the SAP.Section 3 presents the empirical analysis and results, and finally Section 4 discusses targeted suggestions for the policy-making process as relevant to carbon reduction.

Materials and Methods
As the carbon emissions system is a typical composite system composed of a natural and an artificial system, it has the characteristics of openness, being away from an equilibrium state, non-linearity, and the existence of random fluctuation.In order to study the mutual influence of each impact factor in the composite system, to adapt to environmental changes due to economic development and changes, and, ultimately, to be a reasonable state, this study uses complex system theory and methods, including selecting carbon intensity as the system dependent variable and the influencing factors as the independent variables.To solve regional carbon emissions reduction potential, we construct a model to capture the dynamic quantitative relationship between carbon intensity and the system elements and we then examine the prediction and optimization of carbon intensity.

Prediction Method of Carbon Intensity Based on the Simulated Annealing Programming
Due to the complexity of carbon emissions reduction, it is very difficult to establish a dynamic model that describes the evolution of the emissions reduction system.The traditional method is limited by the rationality of the system structure, the complexity of the parameter calculations, and the accuracy of the final conclusion.The SAP is a stochastic search algorithm that is suitable for solving large-scale combination optimization problems [43].It can overcome nonlinear, multi-dimensional, and complex system models and also determine the difficult structure of the equation [44].The simulated annealing programming has been applied in the fields of production scheduling, control engineering, machine learning, neural networks, and signal processing, among others.This study introduces it into the field of carbon emissions prediction.

Simulated Annealing Programming Algorithm
The algorithm is derived from the enlightenment of 'annealing' physical phenomena in thermodynamics.The lower the temperature is, the lower the energy of the object is; when the temperature is low enough, the liquid begins to condense and crystallize.In the crystalline state, the system's energy state lowers to a minimum [45].
The SAP algorithm is based on the simulation of the solid annealing process.Using the Metropolis acceptance criterion, a set of parameter control algorithm processes, namely the cooling schedule, is established to give the algorithm an approximate optimal solution in polynomial time [46].The physical image and statistical properties of the solid annealing process are the physical background of the SAP.The Metropolis acceptance criterion guarantees the algorithm will escape the local optimal 'trap', and the reasonable choice of the cooling schedule is the prerequisite for the application of the algorithm [47].The flow chart of the Simulated Annealing Programming Algorithm is shown in Figure 1.background of the SAP.The Metropolis acceptance criterion guarantees the algorithm will escape the local optimal 'trap', and the reasonable choice of the cooling schedule is the prerequisite for the application of the algorithm [47].The flow chart of the Simulated Annealing Programming Algorithm is shown in Figure 1.It is calculated by the initial solution and the initial value of the temperature control parameters; the current solution repeats the 'new solution → calculate the objective function difference → judge whether to accept → accept or discard' iteration and gradually attenuates the T value.The current solution at the end of the algorithm is the approximate optimal solution.The annealing process is controlled by the cooling schedule parameters, including the initial value of the control parameter and its attenuation factor, the number of iterations per value (called the length of a Markov chain), and the stop condition.
The SAP has the following characteristics [48]: (1) to avoid the loss of the local optimal solution, we accept the worse solution with a certain probability.(2) To improve the reliability of the optimal solution, the SAP algorithm slowly reduces the control parameters, improving the acceptance criterion until the control parameters tend to zero.(3) The objective function requirements are fewer It is calculated by the initial solution and the initial value of the temperature control parameters; the current solution repeats the 'new solution → calculate the objective function difference → judge whether to accept → accept or discard' iteration and gradually attenuates the T value.The current solution at the end of the algorithm is the approximate optimal solution.The annealing process is controlled by the cooling schedule parameters, including the initial value of the control parameter and its attenuation factor, the number of iterations per value (called the length of a Markov chain), and the stop condition.
The SAP has the following characteristics [48]: (1) to avoid the loss of the local optimal solution, we accept the worse solution with a certain probability.(2) To improve the reliability of the optimal solution, the SAP algorithm slowly reduces the control parameters, improving the acceptance criterion until the control parameters tend to zero.(3) The objective function requirements are fewer and are not subject to continuous micro-constraints, under which the definition domain can be arbitrarily set.
(4) To improve the quality and operation speed of the solution, the introduction of implicit parallelism is employed to solve the nonlinear problem.

Theoretical Model and Variable Selection
The theoretical model to solve the nonlinear model of regional carbon intensity and the key factors is as follows: where CI is the value of regional carbon intensity, X it is the t-th year of the i-th key factor, and f (x) is a non-linear function.The key influencing factors are entered into the SAP algorithm in order to find the quantitative relationship between regional carbon intensity and each influencing factor.Five influencing factors of carbon intensity, including economic development, industrial structure, technological advances, energy price, and social investment, are selected as independent variables for their significance.The mechanism of action is described as follows.
(1) Economic Development.According to Nag and Parikh [49] and Ang [50], there is an inverted U-shaped relationship between economic development and carbon intensity.That is, carbon intensity increases with growing GDP up to a threshold level beyond which carbon intensity drops with higher GDP.The development of the economy not only increases the total amount of carbon emissions, but also improves energy efficiency and reduces carbon intensity levels.Thus, the economy development provides the basic support for the development of technology, laying a strong foundation for the reduction of carbon emission.This study uses per capita GDP indicators to characterize the level of economic development.
(2) Industrial Structure.Yang and Liu [51] believe that the industrial structure is an important factor in differences in carbon emissions levels; Gao [52] believes that the industrial structure, the level of industrialization, and its openness to the public have a significant impact on carbon intensity.This study uses the proportion of tertiary industry to GDP to characterize the impact of industrial structure on carbon intensity.(3) Technological Advances.It was generally considered that technological advances have a positive effect on CO 2 emissions reduction [53][54][55].The effect is mainly manifested in two ways.First, technological advances can improve the efficiency of mechanical equipment, increase the use of artificial proficiency, and increase output in simple ways, directly reducing carbon intensity.Secondly, technological advances can improve unit energy output, increase the efficiency of the use of social resources, and thereby reduce the social needs of products, product prices, and overall investment, forcing the elimination of excess capacity and the adjustment of the industrial structure, evenly reducing carbon intensity.In this study, total productivity indicators are used to characterize the impact of energy technology.(4) Energy Price.Greening et al. [56] believed that energy prices had a significant impact on carbon intensity in direct or indirect ways.Chen and Tong [55] argued that the direct impact of energy prices was low.As an indirect impact, energy prices have an effect mainly by adjusting the energy consumption structure, thereby indirectly reducing carbon intensity.As a direct impact, a rise in prices can encourage companies to reduce carbon intensity, reduce their quantity of use, and utilize clean alternative energy sources to improve efficiency.This study uses raw materials, fuel, and the power purchase price index as indicators to characterize energy prices.(5) Social Investment.Research shows that the expansion of investment levels on carbon emissions has bidirectional effects on carbon intensity [57,58].On the one hand, the expansion of investment can increase total production scale, thus increasing carbon intensity; on the other hand, the investment increase can drive economic development, which plays a positive role in promoting technological innovation.Former studies show that a modest increase in investment can promote an increase in carbon intensity, but excessive growth is counterproductive.This study uses total societal fixed asset investment as the indicator to characterize the level of social investment.
Five indicators such as per capita GDP, the proportion of tertiary industry to GDP, total productivity, raw materials, fuel, power purchase price index, and total societal fixed asset investment were used to characterize the level of economic development, industrial structure, technological advances, energy price, and social investment respectively.

Prediction Method of Carbon Emissions Reduction Potential Based on the Genetic Algorithm
A genetic algorithm (GA) is a global optimization search algorithm proposed by a group search strategy and the information exchange between the individuals in the group when the adaptive system is established and gradually developed [59].Its main characteristic is that it does not depend on gradient information and realizes the evolution of the group by iteration under the premise of using the initial population and genetic operation.The result avoids the local optimum, the search scope is wider, and carbon intensity can be optimized under a multi-target constraint scenario.
According to national economic planning and regional economic and social situations, high growth scenarios, benchmark scenarios (planning scenarios), and low growth scenarios (three economic growth model) all set different policy scenarios, and the global parallel optimization of the GA based on the natural genetic mechanism can be used to optimize carbon intensity under different policy scenarios.The specific calculation process is shown in Figure 2 and described as follows.
study uses total societal fixed asset investment as the indicator to characterize the level of social investment.
Five indicators such as per capita GDP, the proportion of tertiary industry to GDP, total productivity, raw materials, fuel, power purchase price index, and total societal fixed asset investment were used to characterize the level of economic development, industrial structure, technological advances, energy price, and social investment respectively.

Prediction Method of Carbon Emissions Reduction Potential Based on the Genetic Algorithm
A genetic algorithm (GA) is a global optimization search algorithm proposed by a group search strategy and the information exchange between the individuals in the group when the adaptive system is established and gradually developed [59].Its main characteristic is that it does not depend on gradient information and realizes the evolution of the group by iteration under the premise of using the initial population and genetic operation.The result avoids the local optimum, the search scope is wider, and carbon intensity can be optimized under a multi-target constraint scenario.
According to national economic planning and regional economic and social situations, high growth scenarios, benchmark scenarios (planning scenarios), and low growth scenarios (three economic growth model) all set different policy scenarios, and the global parallel optimization of the GA based on the natural genetic mechanism can be used to optimize carbon intensity under different policy scenarios.The specific calculation process is shown in Figure 2 and described as follows.(i) Coding.The problem is described as a string through the fitness calculation, wherein each solution corresponds to a fitness value.The project code string is (X 1t , X 2t , . . ., X nt ).(ii) Forming the initial group.The initial population is generated by random method, indicating some things such as a random chromosome, the fixed number of groups, and the project is set to 100.(iii) Calculating the fitness.The fitness is an indicator of measuring the solution (chromosome) as good or bad and relates to the objective function of the algorithm model.The project fitness calculation function is fitness = f (X 1t , X 2t , . . ., X nt ); the corresponding objective function is CI optimization = Min(fitness).
(iv) Choice.This is the process of entering the next generation directly from the old group, choosing the finest individual.(v) Cross.This is the process of crossing and converting some parts of chromosomes (strings).(vi) Mutation.This is to randomly rewrite one or several genes representing a chromosome in a string to represent another trait so that the chromosome is updated.(vii) Termination.This should repeat (iii), (iv), (v), and (vi).When the optimal population reaches a certain requirement or the evolution of the algebra reaches a set value, it terminates.This study intends to stop when the evolutionary algebra reaches 20 to take the optimal value.
Substituting the key factor control range in different scenarios into the multi-objective optimization program of the above-mentioned GA, we get the CI optimization and input this value into Equation (1); namely, we obtain the regional carbon emissions reduction potential.In theory, if several factors can be reasonably regulated and balanced, the regional carbon intensity can reach the optimal value (lowest); that is, there is a certain emissions reduction leeway.At the same time, the value of the independent variables corresponding to the optimal carbon intensity represents the policy path that enables the regional carbon emissions reduction potential.

Results
As there are both theory and models as discussed above, it is necessary to explore the carbon reduction potential of resource-dependent regions in a typical region with a resource-based economy.Shanxi, a representative resource-dependent province, is chosen for investigation through the SAP model.

Study Area
Shanxi Province is located in the middle of the Yellow River Basin, which has a population of 36.64 million and covers an area of 156,700 square kilometers.Shanxi is the largest coal production and transportation province as well as an energy and heavy chemical industry base in China.The total area of coal, which accounts for about 40% of the province, is 64,800 square kilometers.Relying on coal resources, the Shanxi economy has achieved rapid development, but, at the same time, it has also created the phenomenon of a 'coal economy'.The volatility of coal prices has led to swings in regional economies (as shown in Figure 3).(iv) Choice.This is the process of entering the next generation directly from the old group, choosing the finest individual.(iv) Cross.This is the process of crossing and converting some parts of chromosomes (strings).(v) Mutation.This is to randomly rewrite one or several genes representing a chromosome in a string to represent another trait so that the chromosome is updated.(vi) Termination.This should repeat (iii), (iv), (v), and (vi).When the optimal population reaches a certain requirement or the evolution of the algebra reaches a set value, it terminates.This study intends to stop when the evolutionary algebra reaches 20 to take the optimal value.
Substituting the key factor control range in different scenarios into the multi-objective optimization program of the above-mentioned GA, we get the CIoptimization and input this value into Equation (1); namely, we obtain the regional carbon emissions reduction potential.In theory, if several factors can be reasonably regulated and balanced, the regional carbon intensity can reach the optimal value (lowest); that is, there is a certain emissions reduction leeway.At the same time, the value of the independent variables corresponding to the optimal carbon intensity represents the policy path that enables the regional carbon emissions reduction potential.

Results
As there are both theory and models as discussed above, it is necessary to explore the carbon reduction potential of resource-dependent regions in a typical region with a resource-based economy.Shanxi, a representative resource-dependent province, is chosen for investigation through the SAP model.

Study Area
Shanxi Province is located in the middle of the Yellow River Basin, which has a population of 36.64 million and covers an area of 156,700 square kilometers.Shanxi is the largest coal production and transportation province as well as an energy and heavy chemical industry base in China.The total area of coal, which accounts for about 40% of the province, is 64,800 square kilometers.Relying on coal resources, the Shanxi economy has achieved rapid development, but, at the same time, it has also created the phenomenon of a 'coal economy'.The volatility of coal prices has led to swings in regional economies (as shown in Figure 3).At the same time, even more unsettling is that coal resources have brought about environmental pollution and carbon emissions, which has led to Shanxi Province becoming one of the most carbon- At the same time, even more unsettling is that coal resources have brought about environmental pollution and carbon emissions, which has led to Shanxi Province becoming one of the most carbon-dense regions in China.A development model characterized as long-term and highly dependent on coal resources has raised a series of problems and contradictions by creating regional resource waste, capital outflows, the weakening of institutions, and a decrease in innovation capacity.All of this seriously restricts the development of the local society and economics.Thus, the study of carbon emissions reduction potential in Shanxi Province is a great reference point for the theory and practice of carbon reduction in resource-dependent areas.

Prediction Model of Carbon Intensity for Shanxi Province
In this section, the prediction model of carbon intensity for Shanxi Province was deduced through the SAP algorithm.The valuable values of the SAP algorithm were determined and the data required were collected.More importantly, a fit test was conducted to verify the validity of the prediction model.

Variable Setting and Data Resources
The carbon intensity in Province was selected as the dependent variable, expressed as CI; the influence factors on carbon intensity were selected as the independent variables.X 0 , the per capita GDP, is an indicator of economic development; X 1 , the proportion of the tertiary industry in GDP, is an indicator of industrial structure; X 2 , the total factor productivity, is an indicator of energy technology; X 3 , the raw material, fuel, and power purchase price index, is an indicator of energy prices; and X 4 , the total societal investment in fixed assets, is an indicator of social investment.The values of CI and X i are shown as Table 1.Considering the availability of data and the difference in the statistical caliber of each year, this study chooses time-series data from 1990 to 2015.The original data of variables in Table 1 were collected from the Shanxi Statistical Yearbook (1991-2016) and the China Energy Statistical Yearbook (1991-2016) and were standardized to eliminate the effect of dimensions.The detailed CI calculation and the carbon emissions coefficients of fossil fuels are provided in Appendix A. The detailed total factor productivity calculations and raw data are provided in Appendix B.

Prediction Model of Carbon Intensity
The data in Table 1 are input into the SAP model given in Section 2.1.The parameters are defined as: initial temperature T 0 = 200, termination temperature T min = 0, and temperature drop coefficient α = 0.95.The resulting fit coefficient R 2 is 0.91, and the fitted equation is given as follows: ĈI = x 4 +99.96 57.04x 2 (2)

Model Fit Test
In this section, goodness of fit tests are performed to verify the superiority of the SAP, and the results are shown in Table 2.As we can see in Table 2 and Equation ( 2), the maximum absolute residual of the SAP model is 1.809 in 1995, the minimum absolute residual is 0.002 in 2004, and the fitting degree is 0.91, which is greater than 0.9.Thus, the advantages and feasibility of using the SAP model to construct the prediction model of carbon intensity can be verified.The fitting chart of the carbon intensity prediction model is shown as Figure 4. Obviously, the trend line of the prediction value is similar to that of actual value in Figure 4.

Analysis on Reachability of Carbon Reduction Targets of 13th Five-Year Work Plan in Shanxi Province
As already stated, according to the '13th Five-Year Work Plan for Controlling Greenhouse Gas Emissions', Shanxi Province has had a target of carbon emissions reduction that translates to a carbon intensity reduction of 18% by 2020.We now estimate whether the carbon emissions reduction target is achievable using the following formula.Q denotes the reachability of the carbon reduction targets of '13th Five-Year Work Plan' in Shanxi Province.2), the result is CI ^2020 = 3.955 t/10,000 yuan.Table 1 shows the actual carbon intensity in 2015 as 4.87 (t/100,000 yuan); when input it into Equation (3), the result is Q = 18.78%.
As we forecast, the carbon intensity of Shanxi Province in 2020 would be about 18.78% less than that in 2015.This result indicates that the target of carbon reduction, 18% declared in the new fiveyear (i.e., the '13th-Five-Year Work Plan') will be achieved.

Estimation of Carbon Emissions Reduction Potential in Shanxi Province during the 13th Five-Year Work Plan
The above assessment shows that the carbon emissions reduction targets in the '13th Five-Year Work Plan' for Shanxi Province can be achieved.However, there appears to be room for additional reduction; we now use the GA to solve this.After taking into account the energy development trend in Shanxi Province, adopting the suggestions of authorities involved and targets set in the '13th Five-

Analysis on Reachability of Carbon Reduction Targets of 13th Five-Year Work Plan in Shanxi Province
As already stated, according to the '13th Five-Year Work Plan for Controlling Greenhouse Gas Emissions', Shanxi Province has had a target of carbon emissions reduction that translates to a carbon intensity reduction of 18% by 2020.We now estimate whether the carbon emissions reduction target is achievable using the following formula.Q denotes the reachability of the carbon reduction targets of '13th Five-Year Work Plan' in Shanxi Province.
According to the '13th Five-Year Work Plan' of Shanxi Province, the Shanxi Province Development and Reform Commission '13th Five-Year Work Plan' and development program, the Shanxi Province energy development planning, and other '13th Five-Year' control targets, X 0 to X 4 predictions of Shanxi Province are as follows.X 0 = 16,154.7 yuan/person, X 1 = 38.25%,X 2 = 1.15,X 3 = 440.51,X 4 = 561.78 billion yuan; input them into Equation ( 2), the result is ĈI 2020 = 3.955 t/10,000 yuan.Table 1 shows the actual carbon intensity in 2015 as 4.87 (t/100,000 yuan); when input it into Equation (3), the result is Q = 18.78%.
As we forecast, the carbon intensity of Shanxi Province in 2020 would be about 18.78% less than that in 2015.This result indicates that the target of carbon reduction, 18% declared in the new five-year (i.e., the '13th-Five-Year Work Plan') will be achieved.

Estimation of Carbon Emissions Reduction Potential in Shanxi Province during the 13th Five-Year Work Plan
The above assessment shows that the carbon emissions reduction targets in the '13th Five-Year Work Plan' for Shanxi Province can be achieved.However, there appears to be room for additional reduction; we now use the GA to solve this.After taking into account the energy development trend in Shanxi Province, adopting the suggestions of authorities involved and targets set in the '13th Five-Year Work Plan', and holding expert hearings, we reached the expected value of carbon intensity in 2020 in Shanxi Province as follows in Table 3.After the expected value range of Table 3 is input into the GA optimization program (Figure 1), the results are shown in Table 4.As can be seen in Table 4, according to the SAP model of 2020, the carbon intensity prediction value is ĈI = 3.955 t/10,000 yuan; after further optimization of influencing factors input into the GA program, we find a better carbon intensity value, namely, ĈI optimization = 2.985 t/10,000 yuan; this means that for each additional 10,000 yuan of GDP, emissions can be reduced by 0.97 t.

Discussion
This study follows the streams of system analysis, system prediction, system optimization and system decision-making.As such it combines the key influencing factors on carbon emissions found in the literature of economic development, industrial structure, energy technology, energy price, and social investment, and then uses a per capita GDP, tertiary industry as a proportion of GDP, total factor productivity, raw materials, fuel, power purchase price index, and total society fixed assets investment as characterization indicators.
(1) The study introduces a SAP algorithm to solve the quantitative relationship between carbon intensity and five influencing factors in Shanxi Province, and establishes a prediction model for carbon intensity.According to the model's test, the fit degree is strong and the residual values are low.The conclusion is that the SAP algorithm is superior for the prediction of carbon emissions.(2) According to the model's forecast, in Shanxi Province, at the end of the '13th Five-Year Work plan', a carbon intensity of 3.955 t/10,000 yuan will be achieved, which is 18.78% lower than the carbon intensity level in 2015; thus, as promulgated by China in its '13th Five-Year Work Plan for Controlling Greenhouse Gas Emissions', the Shanxi Province carbon intensity reduction task of 18% can be achieved.(3) In Shanxi Province, through the comprehensive evaluation of the policy direction, controlling objectives during the '13th Five-Year Work Plan' and optimizing the influencing factors, the GA can be used to further lower the carbon intensity and thereby further reduce the carbon intensity in Shanxi Province by 0.97 t/10,000 yuan; this means that each additional 10,000 yuan of GDP can achieve a further 0.97 t reduction of emissions in Shanxi Province.
According to the above analysis, we suggest that several measurements are identified in resource-dependent regions to enable carbon reduction potential.
(1) As a big coal region, Shanxi Province has a responsibility to the country; thus, when it implements a carbon emissions reduction policy, several factors should be taken into account, including its actual economic development, population density, and technological advancement, which all contribute to the establishment of a proper low-carbon industrial system.With a low-carbon industrial system, Shanxi not only can solve its outstanding problems, but can also promote equilibrium between coal supply and demand in China, which will promote continuous and scientific development throughout its regional economy, society, and energy and emissions reduction.(2) In order to provide a big impetus for the province's coal supply-side structural reform, which will promote its transformation and the upgrade of its economic structure, we should focus on the province's coal industry by decreasing production capacity, inventories, leverage, and costs, thereby making up for any shortage.In accordance with the law, a number of coalmines have been closed.According to regulations, the reorganization and consolidation of a number of coalmines, reductions, and the withdrawal of the need for excess coal capacity can generate an orderly exit through a market mechanism.In 2016, Shanxi closed 25 coal mines, reducing coal production capacity by 23.25 million tons.By 2020, the province will reduce excess coal capacity by more than 100 million tons.(3) The coal price formation mechanism needs improvement.A coal price mechanism that can correctly reflect the market supply and demand, the scarcity of resources, and the cost of environmental damage should be explored and established.The role of industry associations should be fostered and a coal pricing self-discipline mechanism should be established.To improve such a price self-discipline mechanism, pilot programs can be established through the Yang Coal Group and the Shanxi Coal Group.In addition, a coal price supervision mechanism could be established and promoted.Further, some reforms can be extended, including reform of the electricity price, the power trading system, the electricity plan, and electricity side sales.To that end, we can cultivate the main electricity sales market, establish and improve the electricity market trading mechanism, and also expand the field and scope of direct supply and raise the electricity market scale in and out of the province.We should also innovate the coal trading mechanism, striving to become a national coal-trading pilot in the near future.(4) The pace of scientific and technological innovation in the coal industry should be increased.We should establish the Shanxi Coal Clean Utilization Investment Fund.Among other projects, we should focus on supporting coal and electricity integration, the modern coal chemical industry, coal bed methane (gas) extraction and utilization, and carbon trading and carbon emissions reduction.Moreover, to emphasize the clean and efficient use of coal, we should implement a group of major technology innovated coal-based and low-carbon projects.In addition, among other new energy industries, we should foster the development of wind power, photovoltaic power generation, and biomass power generation, which will hasten the development and utilization of new energy industrialization.(5) In view of our regional economic development characteristics, a flexible carbon emissions reduction mechanism should be built.The Chinese government, over the years, has implemented an extensive administrative control mechanism.Under this system, local government administrations utilize mandatory means to meet the required targets for carbon emissions.
As it can directly help these regions to more readily attain these goals, the system is unlike any other, especially in terms of the control of industrial emissions reduction and the elimination of single enterprises.However, among other challenges, the system's shortcomings lie in issues such as information asymmetry, the high cost of environmental management, the inefficiency of government decision-making, and the inefficiency of government agencies.As for the implementation efficiency and the carbon emissions reduction effect, selecting a reasonable market-based policy tool is still the key to carbon emissions reduction; market-oriented regulatory policies are of two primary types, a price-based carbon tax regime and an aggregate control-based emissions trading scheme.The lack of market-oriented regulation means that enterprises still follow a cost-benefit principle in their investments in carbon emissions reduction, which generates environmental externalities.Compared with the above two policies, the effect of public participation regulation is more lasting, although it is slow.The most important characteristic of the control mechanism is that the driving force is endogenous.If the public's behavior is insufficient or irresponsible, the mechanism will fail.Only through the comprehensive use of these policies can a flexible mechanism of carbon emissions reduction be established, which will engender greater enthusiasm among the relevant subjects and ensure the realization of carbon emissions reduction targets.

Figure 1 .
Figure 1.The flow chart of the Simulated Annealing Programming Algorithm.

Figure 1 .
Figure 1.The flow chart of the Simulated Annealing Programming Algorithm.

Figure 2 .
Figure 2. The flow chart of the Genetic Algorithm.
(i) Coding.The problem is described as a string through the fitness calculation, wherein each solution corresponds to a fitness value.The project code string is (X1t, X2t, …, Xnt).(ii) Forming the initial group.The initial population is generated by random method, indicating some things such as a random chromosome, the fixed number of groups, and the project is set to 100.(iii) Calculating the fitness.The fitness is an indicator of measuring the solution (chromosome) as good or bad and relates to the objective function of the algorithm model.The project fitness calculation function is fitness = f(X1t, X2t, …, Xnt); the corresponding objective function is CIoptimization = Min(fitness).

Figure 2 .
Figure 2. The flow chart of the Genetic Algorithm.

Figure 3 .
Figure 3.An overview of economic development in Shanxi Province in China.Note: the data for GDP and industrial value-added are obtained from the Shanxi Statistical Yearbook (1979-2015).

Figure 3 .
Figure 3.An overview of economic development in Shanxi Province in China.Note: the data for GDP and industrial value-added are obtained from the Shanxi Statistical Yearbook (1979-2015).

Figure 4 .
Figure 4. Fitting chart of the carbon intensity prediction model.
According to the '13th Five-Year Work Plan' of Shanxi Province, the Shanxi Province Development and Reform Commission '13th Five-Year Work Plan' and development program, the Shanxi Province energy development planning, and other '13th Five-Year' control targets, X0 to X4 predictions of Shanxi Province are as follows.X0 = 16,154.7 yuan/person, X1 = 38.25%,X2 = 1.15,X3 = 440.51,X4 = 561.78 billion yuan; input them into Equation (

Figure 4 .
Figure 4. Fitting chart of the carbon intensity prediction model.

Table 1 .
The Carbon Intensity and its influence factors in Shanxi Province from 1990 to 2015.

Table 2 .
Goodness of Fit tests for the Simulated Annealing Programming model.

Table 3 .
Expected/predicted values of carbon intensity in 2020 in Shanxi Province.

Table 4 .
Comparison of carbon intensity prediction models before and after optimization.

Table A2 .
The Total Factor Productivity of Shanxi Province from 1990 to 2015.