3Es System Optimization under Uncertainty Using Hybrid Intelligent Algorithm: A Fuzzy Chance-Constrained Programming Model

Harmonious development of 3Es (economy-energy-environment) system is the key to realize regional sustainable development.The structure and components of 3Es system are analyzed. Based on the analysis of causality diagram, GDP and industrial structure are selected as the target parameters of economy subsystem, energy consumption intensity is selected as the target parameter of energy subsystem, and the emissions of COD, ammonia nitrogen, SO 2 , and NOX and CO2 emission intensity are selected as the target parameters of environment system. Fixed assets investment of three industries, total energy consumption, and investment in environmental pollution control are selected as the decision variables. By regarding the parameters of 3Es system optimization as fuzzy numbers, a fuzzy chance-constrained goal programming (FCCGP) model is constructed, and a hybrid intelligent algorithm including fuzzy simulation and genetic algorithm is proposed for solving it.The results of empirical analysis on Shandong province of China show that the FCCGPmodel can reflect the inherent relationship and evolution law of 3Es system and provide the effective decision-making support for 3Es system optimization.


Introduction
How to balance the economic development, energy security, and environmental protection to achieve sustainable development has been one of the greatest challenges of the 21st century.Referred to as the 3Es (economy-energyenvironment) issues, they often involve simultaneous consideration of the speed and quality of economic development, energy consumption, and pollutants and GHGs (greenhouse gases) emissions.Therefore, it is imperative to optimize 3Es system in order to achieve coordinated development.The coordination is not only the basic assurance for the energy security and sustainable development of economy, but also an effective way to reduce environmental pollution.
There are some researches on the optimization of 3Es system, and a number of optimization models have been proposed.According to the modelling method, 3Es models can be divided into three categories, that is, top-down model, bottom-up model, and mixed model [1].CGE (computable general equilibrium) model is a representative top-down model, which takes economic model as the core technology, and focuses on the interactions between energy price, economic elasticity, energy consumption, and energy production [2][3][4].MARKAL (market allocation model), EFOM (energy flow optimization model), and LEAP (long range energy alternatives planning model) are widely used bottom-up models, and they apply optimization and simulation to predict and optimize the elements of 3Es.The most representative mixed models are NEMS (national energy modelling systems) and IIASA-WEC 3Es (IIASA-WEC energyeconomic-environment model), and they optimize the whole 3Es system through emulation and simulation.Because topdown models do not have to describe energy and environmental technologies in detail, many researchers constructed top-down models.Zhang and Li constructed TH-3EM (Tsinghua hybrid energy-economy-environment assessment model) [5], which consist of CSCGEM, IESOM, and SICGE and other modules.Zhang

Economy
Energy subsystem Environment subsystem Figure 1: Elements and composition structure of 3Es system.[6].DeCarolis et al. applied ESOMs to uncover knife-edge solutions, explore alternative system configurations, and suggest different ways to achieve policy objectives under conditions of deep uncertainty [7].Kambo et al. constructed a linear goal programming model for urban energy-economyenvironment interaction [8], which used energy demand, energy conservation, air pollution, economic cost of energy system, energy imports, and utilization capacity as objectives and constraints.Oliveira and Antunes constructed a multiple objective model of 3Es interactions based on the linear input-output structure of interindustry [9].Wei et al. applied the historical statistical data of Beijing city to construct a multiobjective goal programming model for Beijing's coordinated development of population, resources, environment, and economy [10].Hu et al. put forward an economy-energy-electricity-environment framework and summarized the energy usage management and integrated resource strategic planning of China [11].The above models are all the deterministic programming models with real number coefficients.Actually, because of the uncertainty of system evolution and forecast, it is better to see the parameters of 3Es system model as uncertain parameters.Liu et al. constructed a chance-constrained mixed-integer linear programming model with hybrid interval parameters for nonrenewable energy resources management under uncertainty, with an objective of maximizing economic return under constraints of resources availability and environmental regulations [12].Cai et al. developed an interactive decision support system based on an interval parameter linear programming model to aid decision makers in planning energy management systems [13].Fragnière and Haurie constructed a MARKAL model of Geneva [14], which regarded different policy scenes as a random event with some kind of possibility distribution, and carried out energy planning by random decision tree technology.Lin and Huang developed an interval parameter two-stage stochastic municipal planning model for supporting decisions of energy systems planning and GHG emission management at a municipal level [15].Borges and Antunes considered the right side of inputoutput model as triangular fuzzy numbers and constructed a fuzzy multiple objective decision support model for energyeconomy planning [16].Weber and Martinsen constructed an energy system optimization model based on a fuzzy linear program and achieved an optimal compromise between system sustainability and cost minimization [17].Zhang et al. adopted the fuzzy integral methods to offer an effective decision support for sustainable energy planning of Nanjing city [18].
In summary, 3Es system optimization is one of the significant ways to achieve the regional sustainable development.Extensive mathematical programming models have been put forward to demonstrate the relationships of 3Es system and obtain the optimal decision scheme.From the perspective of reality, uncertain models including interval, stochastic, and fuzzy programming models are better than deterministic models.Since plenty of statistical data is needed in the stochastic programming model to determine the probability distributions of variables, it is rarely used for 3Es optimization directly.Many interval and fuzzy programming models were constructed and applied for micro energy systems optimization, which were transformed into deterministic models with interaction parameters.Few studies constructed specifically the uncertain models of regional 3Es system considering targets of economy, energy, and environment subsystems simultaneously.In this paper, a fuzzy chance-constrained goal programming model (FCCGP) and its hybrid intelligent algorithm are put forward for optimization of 3Es system, and an empirical analysis on Shandong province is conducted for application of the above model and algorithm.
This paper is organized as follows: Section 2 analyzes the structure and components of 3Es system, selects the target parameters and decision variables, and constructs a FCCGP model for optimization of 3Es system.Section 3 provides a hybrid intelligent algorithm for solving the fuzzy model.Section 4 gives an empirical analysis on 3Es system optimization of Shandong province to testify the feasibility and validity of the FCCGP model and hybrid intelligent algorithm.Conclusions are presented in the last section.

FCCGP Model for Optimization of 3Es System
2.1.Economy-Energy-Environment System Analysis.3Es system consists of three subsystems, namely, energy subsystem, economy subsystem, and environment subsystem, which interrelate with and influence each other.The development of each subsystem not only depends on its own structure but also depends on the other subsystems.The elements and composition structure of 3Es system are shown in Figure 1.Energy subsystem is the driving force of economy development and provides the power sources for the economic subsystem.Once the energy supply is threatened and unreliable, it will directly affect the economic output and even lead to economic recession.Meanwhile, economic subsystem plays important positive roles on the energy subsystem: firstly, economic development stimulates energy demand The excessive exploitation and utilization of the energy subsystem can aggravate the environmental bearing capacity, and the degradation and destruction of the environment subsystem will make energy subsystem lose the foundation of exploitation.The development of economic subsystem, especially the industrial parts, accompanied by waste water, waste gas, solid waste, and GHGs, will cause environmental pollution, do harm to human health, and hinder the economic development.
Due to the fact that the impact value of energy investment and environmental pollution on the economic output is difficult to be estimated, the backward effects from energy of environment subsystem to economy subsystem are generally not considered into 3Es system optimization models.Referring to the literatures [8][9][10] and provincial development plan of China, the typical indicators are selected to construct the causality diagram shown in Figure 2. The main causal relationships are as follows: (1) Investments in fixed assets → + GDP → + energy consumption → + CO 2 emissions → + CO 2 emissions intensity.

Selection of Target Parameters and Decision Variables.
Each subsystem of 3Es system has the corresponding targets.
Generally, economy subsystem pursues high GDP growth rate and scientific industrial structure.GDP and three industry structures are selected as the target parameters of economic subsystem: where  0 and   are the GDP of the basic year and the target year and  is the average annual GDP growth rate.
Energy subsystem pursues lower growth rate of energy consumption and lower energy consumption intensity.Based on the target reduction ratio, the energy consumption intensity of the target year can be obtained by where  0 and   are the energy consumption intensities of the basic year and the target year and  is the reduction ratio of energy consumption intensity.Environment subsystem pursues fewer emissions of pollutants and GHGs.The pollutants and GHGs usually include chemical oxygen demand (COD), ammonia nitrogen, sulphur dioxide (SO 2 ), nitrogen oxides (NO X ), and carbon dioxide (CO 2 ).Based on the reduction ratios, they can be computed by the following formula: where  0 ,  0 ,  0 , and  0 are the emissions of COD, ammonia nitrogen, SO 2 , and NO X of the basic year,   ,   ,   , and   are the emissions of the target year,  0 and   are carbon emission intensities of the basic year and the target year, and   ,   ,   ,   , and   are the reduction ratios of the above pollutants and CO 2 .
According to the principles of controllability and exogenesis, the decision variables of the 3Es system include the investments in fixed assets of the three industries, total energy consumption, and total investment in environmental pollution control.
2.3.Formulation of FCCGP Model.3Es system is a complex system and contains multiple optimization objectives.Generally, the target of economy subsystem includes two aspects: GDP of the target year is not lower than the planning value; industrial structure achieves the planning target; that is, the GDP proportions of primary industry and secondary industry are not higher than their target values and that of tertiary industry is not lower than its target value.The target of energy subsystem is that the energy consumption intensity is not higher than the planning value.The target of environment subsystem includes two aspects: the emissions of pollutants including COD, ammonia nitrogen, SO 2 , and NO X are not higher than their planning values; CO 2 emissions intensity is not higher than its planning value.
Due to the small sample of historical data and the uncertainty of forecasting, many parameters of 3Es system can be seen as fuzzy parameters.In this paper, a fuzzy chance-constrained goal programming model (FCCGP) is constructed for optimizing 3Es system.According to the above targets of 3Es system, the constraints consist of two types: (a) the uncertain constraint that each subsystem achieves its targets holding with a certain confidence level; (b) technical constraints.The uncertain constraints are as follows.
(1) GDP of the target year equals the sum of the added values of three industries: where  is GDP of the target year,   is the fixed assets investment in the th industry,   is the added value of unit investment in the th industry, and   are fuzzy variables,  = 1, 2, 3.
(2) The fuzzy chance constraints that three industry structures of the target year should achieve the desired structure are satisfied with certain confidence levels: where Cr{A} express the credibility that fuzzy event  occurs,   is the target value of GDP proportion of the th industry, and   are the predetermined confidence levels by the decision-maker,  = 1, 2, 3.
(3) The constraint that GDP of the target year is not lower than the planning value is satisfied with a certain confidence level: where   is the target value of GDP and  4 is the predetermined confidence level.(4) The constraint that energy consumption intensity of the target year is not higher than the planning value is satisfied with a certain confidence level: where  4 is total energy consumption,  is the target value of energy consumption intensity, and  is the predetermined confidence level.(5) The constraints that the emissions of COD, ammonia nitrogen, SO 2 , NO X , and CO 2 are not higher than their planning values are satisfied with certain confidence levels: where  5 is total investment in environmental pollution control,   are the production coefficients of COD, ammonia nitrogen, SO 2 , and NO  , that is, the productions of COD, ammonia nitrogen, SO (1) Each decision variable should be within a certain range: where   − and   + are the lower limit and upper limit of the th decision variable.
(2) The sum of investments in fixed assets of three industries should be within a certain range: where  − and  + are the lower limit and upper limit of the fixed assets investment of three industries.(3) All the deviation variables are nonnegative: The objective of the model is to minimize the total deviations according to the priority structure and target levels set by the decision-maker: where   1 ,   2 ,   3 are the preemptive priority factors which express the relative importance of economy, energy, and environment targets, respectively.For example, if the priority order of them is economy ≻ energy ≻ environment, then

Solution Method of FCCGP Model
3.1.Fuzzy Simulation.Fuzzy simulation is a technology of sampling test for the fuzzy system model.If  = ( 1 ,  2 , . . .,   ) is a fuzzy vector defined on the credibility space (Θ, P(Θ), Cr),  is predetermined confidence level, and the algorithm of fuzzy simulation for max{ | Cr{(x, ) ≥ } ≥ } is as follows [19].
Step 2. Find the maximal value  such that () ≥  by employing bisection search method, where Step 3. Return  as the estimation of .

Hybrid Intelligent Algorithm.
A hybrid intelligent algorithm can be employed for solving the fuzzy chanceconstrained goal programming model, which integrates fuzzy simulation and genetic algorithm.The procedure of this algorithm is as follows [20].
Step 4. Compare all the objective values of  chromosomes, and rearrange them in descending order.Calculate the cumulative possibility of the sorted chromosome: where  is the evaluation function parameter, 0 <  < 1.
Step Step 6. Set a parameter   as the probability of crossover, and repeat the following process from  = 1 to : generating a random number  from the interval [0, 1], the chromosome x  is selected as a parent if  <   .Make the linear crossover operation on the parent, and output the new population x  = (  1 ,   2 , . . .,   3 ),  = 1, 2, . . ., .
Step 7. Set a parameter   as the probability of mutation, and repeat the following process from  = 1 to : generating a random number  from the interval [0, 1], the chromosome x  is selected as a parent if  <   .Make mutation operation on the parent, and output the new population x  = (  1 ,   2 , . . .,   5 ),  = 1, 2, . . ., , where all the chromosomes satisfy the constraints.
Step 8. Repeat the second to seventh steps for  cycles, and output the best chromosome as the optimal solution.

Empirical Analysis on Shandong Province
4.1.Data of Parameters.As a major economic province of China, Shandong province has achieved significant economic growth, and since 2007, its economic output has been ranked third in 31 provinces.Although many efforts have been made for energy saving and pollutants reduction, the energy efficiency level of Shandong province still falls behind the national average level [21,22], and the major pollutants' emissions are still the largest.In this paper, empirical analysis is conducted on Shandong province, which provides the modelling process of the above FCCGP model and effective decision-making supports for its 3Es system optimization.
The data of the fixed assets investments of three industries at 2000 price index, total energy consumptions, and investment in environmental pollution control from 2000 to 2014 are shown in Table 1.All the indicators show the upward trend; therefore, regression, grey forecasting, and other forecasting methods can be used to forecast the future values theoretically.
The forecast of the fixed assets investment in primary industry is taken as an example to illustrate the application of various methods.Using the annual sequence number as independent variable, the linear regression (LR) equation of the fixed assets investment in primary industry can be obtained as follows: Note: data in the parentheses denote t-test values.t-test and F value show that the regression equation is significant and can be used for short-term forecasting.Inputting  = 16 into the equation, the fixed assets investment in primary industry in 2015 is 117.9 billion Yuan.
By first-order accumulating operation on the raw data, the accumulated values of the fixed assets investment in primary industry can be represented as a row vector  (1) = {77.1,155.3, . . ., 1428.7}, which satisfies the grey exponential law basically.According to grey forecasting (GF) theory, the albino equation can be written as follows: (1)    +  (1) = , where  and  can be derived by using least squares estimation, and the calculation formulas are as follows: The results of  and  are −0.0295 and 75.6988, and the grey forecast model is Because the above formula is the accumulated forecast model, regressive operation should be implemented for restoring the original value.Grey forecast results and actual values of the fixed assets investment in primary industry from 2000 to 2014 are shown in Figure 3.
The results of residual test and posterior variance test are as follows.The maximum relative error is 5.10%, the average relative error is 1.89%, the forecast accuracy is 98.11%, the posterior error ratio is 0.1945, and the small error frequency is 1.The test results show that grey forecast model has high precision and can be used for short-term forecasting.Inputting  = 16 into the model, the fixed assets investment of the first industry in 2015 is 119.6 billion Yuan.
In addition to the above forecasting methods, exponential regression (ER), power function regression (PR), quadratic All the correlation coefficient square values are greater than 0.8, so the above regression equations are significant.Inputting  = 16 into them, the forecast results are 120.2,109.6,130.7, and 111.9, respectively.Then the maximum and minimum values of the above results can be taken as the lower and upper bounds of the variable  1 in 2015.In the same way, the forecast results of the other four decision variables in 2015 and their lower and upper bounds can be obtained as shown in Table 2.
The added values per unit investment of three industries from 2000 to 2014 are shown in Figure 4, and all three curves  show rising trends.By using regression and grey forecast methods, the forecast values in 2015 are shown in Table 3.
The production coefficients of four pollutants and CO 2 emission coefficients from 2000 to 2014 are shown in Table 4.
The forecast values of them by different methods in 2015 are shown in Table 5, where, due to the forecast significance, CO 2 emission coefficients are forecasted based on the samples of 2005-2014.
The removal coefficients of COD, ammonia nitrogen, SO 2 , and NO X from 2000 to 2014 are shown in Table 6, and the forecast values of them in 2015 are shown in Table 7.
All the forecast methods are applied according to their premises, so the above results are reasonable theoretically.Due to less effective forecast results, the above indicators can be seen as interval or fuzzy variables.In order to make full use of the forecast results, all the indicators are considered to obey the triangular fuzzy distributions.The triangular fuzzy variable  ( 1 ,  2 , and  3 ) is fully determined by the triplet ( 1 ,  2 , and  3 ) of crisp numbers with  1 <  2 <  3 , whose membership function is given by    According to the requirements of new normal economy that emphasize structural adjustment, the elimination of energy bottleneck and green low-carbon development, and the effects that Shandong province has achieved at the early time, the GDP proportion of tertiary industry, energy consumption intensity, and CO 2 emission intensity should be paid more attention in 2015.Assuming the predetermined confidence levels of all the targets are  1 = 0.9,  2 = 0.9,  3 = 0.95,  4 = 0.9,  = 0.95,  1 = 0.9,  2 = 0.9,  3 = 0.9,  4 = 0.9, and  5 = 0.95, the priority orders of three subsystems are economy ≻ energy ≻ environment, and the weighting factors of deviations corresponding to economy subsystem and environment subsystem are all equal to 1, and a FCCGP model can be constructed as follows: min { 1 ( where the lower and upper bounds of all the decision variables are shown in Table 2, and the values of triangular distribution parameters of fuzzy variables are shown in Table 8.

Solution and Discussions.
According to the above hybrid intelligent algorithm, the population size  = 50, evaluation function parameter  = 0.05, crossover probability   = 0.3, mutation probability   = 0.2, and generation cycles  = 3000.Running the hybrid intelligent algorithm, the optimal solutions of 3Es system optimization model of Shandong province are shown in Table 9.The results of economic priority scenario show that economic subsystem achieves its targets, but energy and environment subsystems do not achieve their targets.The unachieved indicators include the following: the energy consumption intensity is 54.49tce per million Yuan, the emissions of SO 2 and NO X are 1566900 t and 979302 t, and CO 2 emission intensity is 139.22 tons per million Yuan.
In economic priority scenario, there are greater differences between the optimized results and the targets in energy and environment subsystems.Considering the importance of environment subsystem for sustainable development, environment priority scenario is defined with the priority orders of three subsystems that are environment ≻ economy ≻ energy.Running the hybrid intelligent algorithm, the optimal solutions are shown in Table 10.The results show that environment and energy subsystem achieve their targets, but economic subsystem does not achieve its target.The unachieved indicators include the following: the optimized GDP is 4332.194 billion Yuan, and the industrial structure is 8 : 50 : 42.A scenario endowing the equivalent status for three subsystems is put forward to compare with the first two scenarios.The optimal solutions are shown in Table 11.The results show that all three subsystems do not achieve their targets.The unachieved indicators include the following: the optimized GDP is 4338.292 billion Yuan, the industrial structure is 8 : 50 : 42, the energy consumption intensity is 56.12 tce per million Yuan, and the emissions of SO 2 are 1340146 t.
It can be known from the comparison of three scenarios that the results of the scenario of equivalent status are not improved better than those of environment priority scenario, while environment priority scenario is better than economic priority scenario.Therefore, the optimization results of environment priority scenario can be applied to guide practice The optimization results show that Shandong province should increase the investments in fixed assets of three industries and environmental pollution control in 2015.The investments in fixed assets of the primary, secondary, and tertiary industries in 2015 can be 17%, 4.73%, and 7.59% higher than those in 2014, respectively.The total investments in fixed assets and environmental pollution control in 2015 can be 6.55% and 6.81% higher than those in 2014.Meanwhile, the total energy consumption in 2015 can be merely 2.63% more than that in 2014, which means the growth should be slowed.Practically, Shandong province is undergoing a critical period of structural adjustment at present.The secondary industry consumes considerable amount of energy and causes environmental pollutions and carbon emissions; therefore, it is reasonable to reduce the investment growth of the secondary industry.On the contrary, increasing the investment in the tertiary industry and speeding up its development can reduce energy consumption and environmental pollutions.In addition, the carrying capacity of energy resources in Shandong province has been continuously improving, so it is necessary to reduce the growth of energy consumption while ensuring economic development, which means that Shandong province should strengthen energy conservation, improve energy efficiency, and develop alternative energy sources.

Conclusions
In this study, we have analyzed the structure of 3Es system and selected the typical indicators to demonstrate the main causal relationships.The target parameters and decision variables are set and a FCCGP model is put forward for 3Es system optimization.In the FCCGP model, uncertain parameters are fuzzy variables, and fuzzy chance constraints are satisfied with certain confidence levels.The proposed model is essentially flexible, and any number and form of fuzzy goals and constraints can be contained into it.A hybrid intelligent algorithm including fuzzy simulation and genetic algorithm is put forward for solving the uncertain model.Depending upon the nature of the study region, the user can modify the model and algorithm inputs to get desired outputs.After the empirical analysis on Shandong province, we can conclude that the parameters of 3Es system optimization can be regarded as triangular fuzzy numbers, and FCCGP models with fuzzy variables can demonstrate the relations of interaction among three subsystems.The results of 3Es system optimization of Shandong province testify the feasibility of FCCGP models, which can be applied to guide the practice of 3Es system optimization so as to realize their harmonious development.According to the optimization results, Shandong province should obey the rule of environment priority.In this scenario, 3Es system can realize relatively more scientific development, although it may lower the economic development to a certain extent.
Limited by the definition of research object and availability of data, in this paper, we neither include imports and exports, household consumption and other contents into economic subsystem, nor take into consideration land, forest, mineral, water resources, population, science and technology, and other subsystems.This weakened the complex relationship between each subsystem in some certain extent and might affect the objectivity of the programming results, which should be improved in the future studies.

Figure 3 :
Figure 3: Comparison of grey forecast results and actual values.

Figure 4 :
Figure 4: The added values per unit investment of three industries from 2000 to 2014.
et al. applied IESOM to simulate energy consumptions, energy structures, and carbon emissions in the different condition of emission mitigation , , , and  are the target emissions of COD, ammonia nitrogen, SO 2 , and NO X ,  is the target CO 2 emission intensity,   ,   , and  are fuzzy variables,  = 1, 2, 3, 4, and  1 ,  2 , . ..,  5 are the predetermined confidence levels by the decision-maker.The technical constraints include the following.
2, and NO X per unit of the added value of secondary industry,  = 1, 2, 3, 4, respectively,   are the removal coefficients of COD, ammonia nitrogen, SO 2 , and NO X , that is, the removals of COD, ammonia nitrogen, SO 2 , and NO X per unit of investment in environmental pollution control,  = 1, 2, 3, 4, respectively,  is CO 2 emission coefficient, that is, CO 2 emissions per unit of energy consumption, and  1 ≫  2 ≫  3 .  and   are the weighting factors corresponding to the deviations of economy subsystem and environment subsystem,  = 1, 2, 3, 4,  = 1, 2, . . ., 5, respectively, and they are all positive.

Table 1 :
Data of the fixed assets investments, energy consumptions, and investment in environmental pollution control of Shandong province from 2000 to 2014.
5. Select chromosomes by spinning the roulette wheel; that is, repeat the following process  times: generating a random number  in (0,   ), the rth chromosome is selected if  −1 <  ≤   .Output the new population x  = ( Sources: Shandong Statistics Yearbook 2001-2015 and China Energy Statistical Yearbook 2001-2015.Note: -indicates the lack of data.

Table 2 :
Forecast results of decision variables in 2015.

Table 3 :
Forecast results of the added values per unit investment of three industries in 2015 (unit: Yuan/Yuan).

Table 4 :
The production coefficients of pollutants and CO 2 emission coefficients from 2000 to 2014.
Note: -indicates the lack of data.

Table 5 :
Forecast results of the production coefficients of four pollutants and CO 2 emissions coefficients in 2015.: * indicates that the value does not conform to the reality and should be discarded. Note
Note: -indicates the lack of data.

Table 7 :
Forecast results of the removal coefficients of pollutants in 2015 (unit: t/10 8 Yuan).*indicates that the value does not conform to the reality and × indicates the corresponding method does not pass significant test.where 1 ,  2 , and  3 present the minimum, average, and maximum of all the effective forecast results, respectively.The values of the parameters in triangular fuzzy variables of FCCGP model are shown in Table8. Note:

Table 8 :
The values of the parameters in triangular fuzzy variables of FCCGP model.

Table 9 :
Solution results of FCCGP model of Shandong province in economic priority scenario.

Table 10 :
Solution results of FCCGP model of Shandong province in environment priority scenario.

Table 11 :
Solution results of FCCGP model of Shandong province in the scenario of equivalent status.optimization of Shandong province.In 2015, 13.07, 171.62, and 200.8 ten billion Yuan can be invested into the primary, secondary, and tertiary industries, respectively, total energy consumption should be controlled below 370.801 million tce, and 130.76 ten billion Yuan can be invested into environmental pollution control.