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

Authors: Jiekun Song; Kaixin Zhang; Zijian Cao
Journal: Scientific Programming
 (2016)

---

## Abstract

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, SO2, 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 FCCGP model can reflect the inherent relationship and evolution law of 3Es system and provide the effective decision-making support for 3Es system optimization.


---

## Body

## 1. 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-energy-environment) 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–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 energy-economic-environment model), and they optimize the whole 3Es system through emulation and simulation. Because top-down 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 et al. applied IESOM to simulate energy consumptions, energy structures, and carbon emissions in the different condition of emission mitigation [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-economy-environment 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 input-output model as triangular fuzzy numbers and constructed a fuzzy multiple objective decision support model for energy-economy 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: Section2 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.


## 2. 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 Figure1. 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 and provides the development space for energy subsystem; secondly, economic growth provides the material foundation of human, material, and financial resources for the development of energy subsystem; thirdly, economic growth drives the technological investment and progress in energy subsystem and consequently promotes the development of energy subsystem. Energy subsystem and environmental subsystem both conflict and coordinate with each other. 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.Figure 1
Elements and composition structure of 3Es system.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–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 → + CO2 emissions → + CO2 emissions intensity.(2)
GDP→ + GDP of secondary industry → + pollutants production → + pollutants emissions.(3)
GDP proportion of tertiary→ − energy consumption → + energy consumption intensity.(4)
Investment in environmental pollution control→ + pollutants removals → − pollutants emissions.Figure 2
Causality diagram of 3Es system in Shandong province.


### 2.2. 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:(1)Gt=G01+dt,where G0 and Gt are the GDP of the basic year and the target year and d 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(2)et=e01-r,where e0 and et are the energy consumption intensities of the basic year and the target year and r 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 (SO2), nitrogen oxides (NOX), and carbon dioxide (CO2). Based on the reduction ratios, they can be computed by the following formula:(3)Dt=D01-rD,At=A01-rA,St=S01-rS,Nt=N01-rN,Ct=C01-rC,where D0, A0, S0, and N0 are the emissions of COD, ammonia nitrogen, SO2, and NOX of the basic year, Dt, At, St, and Nt are the emissions of the target year, C0 and Ct are carbon emission intensities of the basic year and the target year, and rD, rA, rS, rN, and rC are the reduction ratios of the above pollutants and CO2.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, SO2, and NOX are not higher than their planning values; CO2 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:(4)G=∑i=13Ii×xi,where G is GDP of the target year, xi is the fixed assets investment in the ith industry, Ii is the added value of unit investment in the ith industry, and Ii are fuzzy variables, i=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:(5)CrI1×x1G-g1≤d1+≥α1,CrI2×x2G-g2≤d2+≥α2,Crg3-I3×x3G≤d3-≥α3,where CrA express the credibility that fuzzy event A occurs, gi is the target value of GDP proportion of the ith industry, and αi are the predetermined confidence levels by the decision-maker, i=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:(6)CrGt-G≤d4-≥α4,where Gt 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:(7)Crx4G-e≤d5+≥β,where x4 is total energy consumption, e is the target value of energy consumption intensity, and β is the predetermined confidence level.(5)
The constraints that the emissions of COD, ammonia nitrogen, SO2, NOX, and CO2 are not higher than their planning values are satisfied with certain confidence levels:(8)CrP1×I2×x2-E1×x5-D≤d6+≥γ1,CrP2×I2×x2-E2×x5-A≤d7+≥γ2,CrP3×I2×x2-E3×x5-S≤d8+≥γ3,CrP4×I2×x2-E4×x5-N≤d9+≥γ4,Crc×x4G-C≤d10+≥γ5,where x5 is total investment in environmental pollution control, Pi are the production coefficients of COD, ammonia nitrogen, SO2, and NOX, that is, the productions of COD, ammonia nitrogen, SO2, and NOX per unit of the added value of secondary industry, i=1,2,3,4, respectively, Ei are the removal coefficients of COD, ammonia nitrogen, SO2, and NOX, that is, the removals of COD, ammonia nitrogen, SO2, and NOX per unit of investment in environmental pollution control, i=1,2,3,4, respectively, c is CO2 emission coefficient, that is, CO2 emissions per unit of energy consumption, D, A, S, and N are the target emissions of COD, ammonia nitrogen, SO2, and NOX, C is the target CO2 emission intensity, Pi, Ei, and c are fuzzy variables, i=1,2,3,4, and γ1,γ2,…,γ5 are the predetermined confidence levels by the decision-maker.The technical constraints include the following.(1)
Each decision variable should be within a certain range:(9)xi-≤xi≤xi+,i=1,2,…,5,where xi- and xi+ are the lower limit and upper limit of the ith decision variable.(2)
The sum of investments in fixed assets of three industries should be within a certain range:(10)I-≤∑i=13xi≤I+,where I- and I+ are the lower limit and upper limit of the fixed assets investment of three industries.(3)
All the deviation variables are nonnegative:(11)d1+,d2+,d3-,d4-,d5+,d6+,d7+,d8+,d9+,d10+≥0.The objective of the model is to minimize the total deviations according to the priority structure and target levels set by the decision-maker:(12)minPi1w1d1++w2d2++w3d3-+w4d4-,Pi2d5+,Pi3u1d6++u2d7++u3d8++u4d9++u5d10+,where Pi1,Pi2,Pi3are 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 i1=1, i2=2, i3=3, and P1≫P2≫P3. wi and uk are the weighting factors corresponding to the deviations of economy subsystem and environment subsystem, i=1,2,3,4,  k=1,2,…,5, respectively, and they are all positive.


## 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 Figure1. 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 and provides the development space for energy subsystem; secondly, economic growth provides the material foundation of human, material, and financial resources for the development of energy subsystem; thirdly, economic growth drives the technological investment and progress in energy subsystem and consequently promotes the development of energy subsystem. Energy subsystem and environmental subsystem both conflict and coordinate with each other. 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.Figure 1
Elements and composition structure of 3Es system.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–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 → + CO2 emissions → + CO2 emissions intensity.(2)
GDP→ + GDP of secondary industry → + pollutants production → + pollutants emissions.(3)
GDP proportion of tertiary→ − energy consumption → + energy consumption intensity.(4)
Investment in environmental pollution control→ + pollutants removals → − pollutants emissions.Figure 2
Causality diagram of 3Es system in Shandong province.


## 2.2. 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:(1)Gt=G01+dt,where G0 and Gt are the GDP of the basic year and the target year and d 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(2)et=e01-r,where e0 and et are the energy consumption intensities of the basic year and the target year and r 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 (SO2), nitrogen oxides (NOX), and carbon dioxide (CO2). Based on the reduction ratios, they can be computed by the following formula:(3)Dt=D01-rD,At=A01-rA,St=S01-rS,Nt=N01-rN,Ct=C01-rC,where D0, A0, S0, and N0 are the emissions of COD, ammonia nitrogen, SO2, and NOX of the basic year, Dt, At, St, and Nt are the emissions of the target year, C0 and Ct are carbon emission intensities of the basic year and the target year, and rD, rA, rS, rN, and rC are the reduction ratios of the above pollutants and CO2.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, SO2, and NOX are not higher than their planning values; CO2 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:(4)G=∑i=13Ii×xi,where G is GDP of the target year, xi is the fixed assets investment in the ith industry, Ii is the added value of unit investment in the ith industry, and Ii are fuzzy variables, i=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:(5)CrI1×x1G-g1≤d1+≥α1,CrI2×x2G-g2≤d2+≥α2,Crg3-I3×x3G≤d3-≥α3,where CrA express the credibility that fuzzy event A occurs, gi is the target value of GDP proportion of the ith industry, and αi are the predetermined confidence levels by the decision-maker, i=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:(6)CrGt-G≤d4-≥α4,where Gt 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:(7)Crx4G-e≤d5+≥β,where x4 is total energy consumption, e is the target value of energy consumption intensity, and β is the predetermined confidence level.(5)
The constraints that the emissions of COD, ammonia nitrogen, SO2, NOX, and CO2 are not higher than their planning values are satisfied with certain confidence levels:(8)CrP1×I2×x2-E1×x5-D≤d6+≥γ1,CrP2×I2×x2-E2×x5-A≤d7+≥γ2,CrP3×I2×x2-E3×x5-S≤d8+≥γ3,CrP4×I2×x2-E4×x5-N≤d9+≥γ4,Crc×x4G-C≤d10+≥γ5,where x5 is total investment in environmental pollution control, Pi are the production coefficients of COD, ammonia nitrogen, SO2, and NOX, that is, the productions of COD, ammonia nitrogen, SO2, and NOX per unit of the added value of secondary industry, i=1,2,3,4, respectively, Ei are the removal coefficients of COD, ammonia nitrogen, SO2, and NOX, that is, the removals of COD, ammonia nitrogen, SO2, and NOX per unit of investment in environmental pollution control, i=1,2,3,4, respectively, c is CO2 emission coefficient, that is, CO2 emissions per unit of energy consumption, D, A, S, and N are the target emissions of COD, ammonia nitrogen, SO2, and NOX, C is the target CO2 emission intensity, Pi, Ei, and c are fuzzy variables, i=1,2,3,4, and γ1,γ2,…,γ5 are the predetermined confidence levels by the decision-maker.The technical constraints include the following.(1)
Each decision variable should be within a certain range:(9)xi-≤xi≤xi+,i=1,2,…,5,where xi- and xi+ are the lower limit and upper limit of the ith decision variable.(2)
The sum of investments in fixed assets of three industries should be within a certain range:(10)I-≤∑i=13xi≤I+,where I- and I+ are the lower limit and upper limit of the fixed assets investment of three industries.(3)
All the deviation variables are nonnegative:(11)d1+,d2+,d3-,d4-,d5+,d6+,d7+,d8+,d9+,d10+≥0.The objective of the model is to minimize the total deviations according to the priority structure and target levels set by the decision-maker:(12)minPi1w1d1++w2d2++w3d3-+w4d4-,Pi2d5+,Pi3u1d6++u2d7++u3d8++u4d9++u5d10+,where Pi1,Pi2,Pi3are 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 i1=1, i2=2, i3=3, and P1≫P2≫P3. wi and uk are the weighting factors corresponding to the deviations of economy subsystem and environment subsystem, i=1,2,3,4,  k=1,2,…,5, respectively, and they are all positive.


## 3. Solution Method of FCCGP Model

### 3.1. Fuzzy Simulation

Fuzzy simulation is a technology of sampling test for the fuzzy system model. Ifξ=fξ1,ξ2,…,ξn is a fuzzy vector defined on the credibility space (Θ,P(Θ),Cr), α is predetermined confidence level, and the algorithm of fuzzy simulation for max⁡f-∣Crfx,ξ≥f-≥α} is as follows [19].Step 1.
Generateθk from Θ such that Crθk≥ε/2 (k=1,2,…,N), where ε is a sufficiently small number.Step 2.
Find the maximal valuer such that Lr≥α by employing bisection search method, where(13)Lr=12max1≤k≤Nvk∣fx,ξθk≥r+min1≤k≤N1-vk∣fx,ξθk<r.Step 3.
Returnr as the estimation of f¯.


### 3.2. Hybrid Intelligent Algorithm

A hybrid intelligent algorithm can be employed for solving the fuzzy chance-constrained goal programming model, which integrates fuzzy simulation and genetic algorithm. The procedure of this algorithm is as follows [20].Step 1.
Initialize a population withm chromosomes randomly, and the rth chromosome xr=(x1r,x2r,…,x5r) satisfies the constraints as follows:(14)xi-≤xir≤xi+,i=1,2,…,5;r=1,2,…,m,I-≤∑i=13xir≤I+.Step 2.
According to the credibility distributions of all the fuzzy variables includingIi, Pj, Ej, and c,i=1,2,3,j=1,2,3,4, fuzzy simulation is applied to generate the data of uncertain functions (e1,e2,…,e10) like(15)e1=max⁡r1∣CrI1×x1G≤r1≥α1,e2=max⁡r2∣CrI2×x2G≤r2≥α2,e3=max⁡r3∣CrI3×x3G≥r3≥α3,e4=max⁡Crr4∣G≥r4≥α4,e5=max⁡r5∣Crx4G≤r5≥β,e6=max⁡r6∣CrP1×I2×x2-E1×x5≤r6≥γ1,e7=max⁡r7∣CrP2×I2×x2-E2×x5≤r7≥γ2,e8=max⁡r8∣CrP3×I2×x2-E3×x5≤r8≥γ3,e9=max⁡r9∣CrP4×I2×x2-E4×x5≤r9≥γ4,e10=max⁡r10∣Crc×x4G≤r10≥γ5.Step 3.
Calculate the values of deviation variables of each chromosome:d1+=e1-g1∨0, d2+=e2-g2∨0, d3-=g3-e3∨0, d4-=Gt-e4∨0, d5+=e5-e∨0, d6+=e6-D∨0, d7+=e7-A∨0, d8+=e8-S∨0, d9+=e9-N∨0, and d10+=e10-C∨0, and calculate their weighted sum(16)Pi1w1d1++w2d2++w3d3-+w4d4-+Pi2d5++Pi3u1d6++u2d7++u3d8++u4d9++u5d10+as the objective value.Step 4.
Compare all the objective values ofm chromosomes, and rearrange them in descending order. Calculate the cumulative possibility of the sorted chromosome:(17)pr=∑t=1ra1-at-1,r=1,2,…,m,where a is the evaluation function parameter, 0<a<1.Step 5.
Select chromosomes by spinning the roulette wheel; that is, repeat the following processm times: generating a random number p in (0, pm), therth chromosome is selected if pr-1<p≤pr. Output the new population xf=(x1f,x2f,…,x5f), f=1,2,…,m.Step 6.
Set a parameterpj as the probability of crossover, and repeat the following process from f=1 to m: generating a random number q from the interval [0,1], the chromosome xf is selected as a parent if q<pj. Make the linear crossover operation on the parent, and output the new population xc=(x1c,x2c,…,x3c), c=1,2,…,m.Step 7.
Set a parameterpb as the probability of mutation, and repeat the following process from c=1 to m: generating a random number s from the interval [0,1], the chromosome xc is selected as a parent if s<pb. Make mutation operation on the parent, and output the new population xb=(x1b,x2b,…,x5b), b=1,2,…,m, where all the chromosomes satisfy the constraints.Step 8.
Repeat the second to seventh steps forT cycles, and output the best chromosome as the optimal solution.


## 3.1. Fuzzy Simulation

Fuzzy simulation is a technology of sampling test for the fuzzy system model. Ifξ=fξ1,ξ2,…,ξn is a fuzzy vector defined on the credibility space (Θ,P(Θ),Cr), α is predetermined confidence level, and the algorithm of fuzzy simulation for max⁡f-∣Crfx,ξ≥f-≥α} is as follows [19].Step 1.
Generateθk from Θ such that Crθk≥ε/2 (k=1,2,…,N), where ε is a sufficiently small number.Step 2.
Find the maximal valuer such that Lr≥α by employing bisection search method, where(13)Lr=12max1≤k≤Nvk∣fx,ξθk≥r+min1≤k≤N1-vk∣fx,ξθk<r.Step 3.
Returnr as the estimation of f¯.


## 3.2. Hybrid Intelligent Algorithm

A hybrid intelligent algorithm can be employed for solving the fuzzy chance-constrained goal programming model, which integrates fuzzy simulation and genetic algorithm. The procedure of this algorithm is as follows [20].Step 1.
Initialize a population withm chromosomes randomly, and the rth chromosome xr=(x1r,x2r,…,x5r) satisfies the constraints as follows:(14)xi-≤xir≤xi+,i=1,2,…,5;r=1,2,…,m,I-≤∑i=13xir≤I+.Step 2.
According to the credibility distributions of all the fuzzy variables includingIi, Pj, Ej, and c,i=1,2,3,j=1,2,3,4, fuzzy simulation is applied to generate the data of uncertain functions (e1,e2,…,e10) like(15)e1=max⁡r1∣CrI1×x1G≤r1≥α1,e2=max⁡r2∣CrI2×x2G≤r2≥α2,e3=max⁡r3∣CrI3×x3G≥r3≥α3,e4=max⁡Crr4∣G≥r4≥α4,e5=max⁡r5∣Crx4G≤r5≥β,e6=max⁡r6∣CrP1×I2×x2-E1×x5≤r6≥γ1,e7=max⁡r7∣CrP2×I2×x2-E2×x5≤r7≥γ2,e8=max⁡r8∣CrP3×I2×x2-E3×x5≤r8≥γ3,e9=max⁡r9∣CrP4×I2×x2-E4×x5≤r9≥γ4,e10=max⁡r10∣Crc×x4G≤r10≥γ5.Step 3.
Calculate the values of deviation variables of each chromosome:d1+=e1-g1∨0, d2+=e2-g2∨0, d3-=g3-e3∨0, d4-=Gt-e4∨0, d5+=e5-e∨0, d6+=e6-D∨0, d7+=e7-A∨0, d8+=e8-S∨0, d9+=e9-N∨0, and d10+=e10-C∨0, and calculate their weighted sum(16)Pi1w1d1++w2d2++w3d3-+w4d4-+Pi2d5++Pi3u1d6++u2d7++u3d8++u4d9++u5d10+as the objective value.Step 4.
Compare all the objective values ofm chromosomes, and rearrange them in descending order. Calculate the cumulative possibility of the sorted chromosome:(17)pr=∑t=1ra1-at-1,r=1,2,…,m,where a is the evaluation function parameter, 0<a<1.Step 5.
Select chromosomes by spinning the roulette wheel; that is, repeat the following processm times: generating a random number p in (0, pm), therth chromosome is selected if pr-1<p≤pr. Output the new population xf=(x1f,x2f,…,x5f), f=1,2,…,m.Step 6.
Set a parameterpj as the probability of crossover, and repeat the following process from f=1 to m: generating a random number q from the interval [0,1], the chromosome xf is selected as a parent if q<pj. Make the linear crossover operation on the parent, and output the new population xc=(x1c,x2c,…,x3c), c=1,2,…,m.Step 7.
Set a parameterpb as the probability of mutation, and repeat the following process from c=1 to m: generating a random number s from the interval [0,1], the chromosome xc is selected as a parent if s<pb. Make mutation operation on the parent, and output the new population xb=(x1b,x2b,…,x5b), b=1,2,…,m, where all the chromosomes satisfy the constraints.Step 8.
Repeat the second to seventh steps forT cycles, and output the best chromosome as the optimal solution.


## 4. 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 Table1. All the indicators show the upward trend; therefore, regression, grey forecasting, and other forecasting methods can be used to forecast the future values theoretically.Table 1
Data of the fixed assets investments, energy consumptions, and investment in environmental pollution control of Shandong province from 2000 to 2014.
Year Fixed assets investment (108 Yuan) Total energy consumption (104 tce) Investment in environmental pollution control (108 Yuan) Primary industry Secondary industry Tertiary industry 2000 77.1 1176.7 1288.8 11361.9 76.2 2001 78.2 1193.1 1306.9 9955.0 89.8 2002 79.1 1206.3 1321.3 14599.0 144.6 2003 81.4 1241.2 1359.6 16624.8 156.4 2004 87.4 1333.1 1460.2 19623.7 191.9 2005 89.9 1371.1 1501.8 24162.0 238.8 2006 91.5 1395.8 1528.8 26759.3 258.1 2007 95.2 1452.2 1590.6 29176.6 320.8 2008 102.5 1564.0 1713.1 30570.0 432.2 2009 99.4 1515.5 1660.0 32420.2 459.5 2010 103.0 1570.7 1720.4 34807.8 483.9 2011 110.0 1678.1 1838.1 37132.0 614.1 2012 110.9 1692.2 1853.5 38899.3 739.1 2013 111.4 1698.8 1860.7 35358.0 848.0 2014 111.7 1703.9 1866.3 36511.0 823.8 Sources:Shandong Statistics Yearbook 2001–2015 and China Energy Statistical Yearbook 2001–2015. Note: — indicates the lack of data.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:(18)FI1=2.8271t+72.629560.781821.5116R2=0.9727,F=426.7498,sig.=1.51e-11.Note: data in the parentheses denotet-test values.t-test andF value show that the regression equation is significant and can be used for short-term forecasting. Inputting t=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 vectorX1=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:(19)dx1dt+ax1=b,where a and b can be derived by using least squares estimation, and the calculation formulas are as follows:(20)ab=BTB-1BTyn,B=-12x11+x211-12x21+x311⋮⋮-12xn-11+xn11,yn=x20x30⋮xn0.The results of a and b are −0.0295 and 75.6988, and the grey forecast model is(21)x^t+11=x10-bae-at+ba=2641.9629e0.0295t+75.6988.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.Figure 3
Comparison of grey forecast results and actual values.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. Inputtingt=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 regression (QR), and cubic curve regression (CR) can also be applied to fit the relation between the fixed assets investment in primary industry and the annual sequence number. The regression results are obtained as follows:(22)ER:FI1=74.245e0.0301t,R2=0.9691,PR:FI1=69.661t0.1636,R2=0.8773,QR:FI1=0.0274t2+3.2653t+71.391,R2=0.9736,CR:FI1=-0.0191t3+0.4308t2+0.2372t+76.065,R2=0.9827.All the correlation coefficient square values are greater than 0.8, so the above regression equations are significant. Inputtingt=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 x1 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.Table 2
Forecast results of decision variables in 2015.
Forecast method Fixed assets investment (108 Yuan) Total energy consumption (104 tce) Investment in environmental pollution control (108 Yuan) Primary industry Secondary industry Tertiary industry LR 117.9 1797.7 1969.1 43613.5 851.2 GF 119.6 1824.8 1998.7 45597.1 1112.7 ER 120.2 1833.1 2008.0 51777.5 1211.8 PR 109.6 1672.5 1831.9 40551.8 735.2 QR 130.7 1778.8 1948.3 38334.4 984.6 CR 111.9 1707.3 1870.3 33948.2 970.3 Lower bound 109.6 1672.5 1831.9 33948.2 735.2 Upper bound 130.7 1833.1 2008.0 51777.5 1211.8The added values per unit investment of three industries from 2000 to 2014 are shown in Figure4, and all three curves show rising trends. By using regression and grey forecast methods, the forecast values in 2015 are shown in Table 3.Table 3
Forecast results of the added values per unit investment of three industries in 2015 (unit: Yuan/Yuan).
Forecast method The added values per unit investment Primary industry Secondary industry Tertiary industry LR 20.57 14.77 7.48 GF 20.64 16.68 8.31 ER 20.66 17.31 8.37 PR 19.74 12.72 6.38 QR 20.80 15.83 8.14 CR 21.21 15.62 8.01Figure 4
The added values per unit investment of three industries from 2000 to 2014.The production coefficients of four pollutants and CO2 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, CO2 emission coefficients are forecasted based on the samples of 2005–2014.Table 4
The production coefficients of pollutants and CO2 emission coefficients from 2000 to 2014.
Year The production coefficients of pollutants (t/108 Yuan) CO2 emission coefficients (t/104 Yuan) COD Ammonia nitrogen SO2 N O X 2000 531.08 — 427.73 — 1.5588 2001 473.72 13.63 383.51 — 1.5956 2002 385.00 12.52 340.54 — 1.2807 2003 338.53 10.41 339.03 — 1.5982 2004 280.36 9.58 293.57 — 1.6750 2005 247.16 9.03 284.08 — 2.2665 2006 199.37 13.91 311.22 121.74 2.0849 2007 180.69 7.29 251.87 89.25 1.9663 2008 182.83 8.00 280.38 77.03 1.9407 2009 165.39 9.27 273.67 73.96 1.7856 2010 154.63 11.32 268.44 69.16 1.6662 2011 150.94 7.00 314.13 68.12 1.5793 2012 134.31 5.69 239.98 62.21 1.5230 2013 98.60 4.26 247.97 59.57 1.2420 2014 86.96 3.42 268.11 61.15 1.2194 Note: — indicates the lack of data.Table 5
Forecast results of the production coefficients of four pollutants and CO2 emissions coefficients in 2015.
Forecast method The production coefficients of pollutants (t/108 Yuan) CO2 emissions coefficient (t/t) COD Ammonia nitrogen SO2 N O X LR 14.33∗ 4.24 225.62 45.22 1.1026 GF 74.65 5.32 237.77 53.96 1.1842 ER 80.78 4.51 234.10 50.72 1.1640 PR 114.70 5.84 250.08 56.27 1.3464 QR 127.19 2.91 282.99 70.03 1.0765 CR 57.94 −0.27∗ 239.55 48.1 0.9863 Note:∗ indicates that the value does not conform to the reality and should be discarded.The removal coefficients of COD, ammonia nitrogen, SO2, and NOX from 2000 to 2014 are shown in Table 6, and the forecast values of them in 2015 are shown in Table 7.Table 6
The removal coefficients of four pollutants from 2000 to 2014 (unit: t/108 Yuan).
Year The removal coefficients of pollutants COD Ammonia nitrogen SO2 N O X 2000 22312.87 — 4204.09 — 2001 19102.87 386.51 4054.03 — 2002 11278.36 248.46 2884.60 — 2003 10841.76 247.79 3612.67 — 2004 9002.71 224.17 3286.14 — 2005 7506.24 195.67 3160.82 — 2006 6529.95 449.66 5691.74 13.98 2007 5667.44 204.30 4286.80 11.91 2008 4968.03 206.52 5140.90 18.51 2009 4824.66 272.07 5948.39 50.10 2010 4789.07 363.43 6514.73 27.20 2011 4409.67 195.86 7002.29 29.71 2012 3601.40 145.61 4683.00 60.08 2013 2523.13 103.74 5034.74 77.42 2014 2498.19 92.75 6540.74 84.55 Note: — indicates the lack of data.Table 7
Forecast results of the removal coefficients of pollutants in 2015 (unit: t/108 Yuan).
Forecast method The removal coefficients of pollutants COD Ammonia nitrogen SO2 N O X LR −1206.00∗ 135.80 6503.00 86.63 GF 1741.53 × × 110.45 ER 2085.52 127.80 6655.57 112.52 PR 3099.84 161.19 5659.68 72.30 QR 4924.04 67.34 6539.95 103.59 CR 72.02∗ −66.19∗ 4731.15 113.16 Note:∗ indicates that the value does not conform to the reality and × indicates the corresponding method does not pass significant test.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 variableT (u1, u2, and u3) is fully determined by the triplet (u1, u2, and u3) of crisp numbers with u1<u2<u3, whose membership function is given by(23)
μ
x
=
x
-
u
1
u
2
-
u
1
,
if
u
1
≤
x
≤
u
2
,
x
-
u
3
u
2
-
u
3
,
if
u
2
≤
x
≤
u
3
,
0
,
otherwise,where u1, u2, and u3 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 Table 8.Table 8
The values of the parameters in triangular fuzzy variables of FCCGP model.
Variable u 1 u 2 u 3 I 1 19.74 20.60 21.21 I 2 12.72 15.49 17.31 I 3 6.38 7.78 8.37 P 1 57.94 91.05 127.19 P 2 2.91 4.56 5.84 P 3 225.62 245.02 282.99 P 4 45.22 54.05 70.03 c 0.9863 1.1433 1.3464 E 1 1741.53 2962.73 4924.04 E 2 67.34 123.03 161.19 E 3 4731.15 6017.87 6655.57 E 4 72.30 99.78 113.16


### 4.2. FCCGP Model

According to the 12th five-year plan of Shandong province, the annual average growth rate of GDP is 9%, the industrial structure in 2015 is 7: 48: 45, the annual average growth rate of fixed assets investment is set about 15%, energy consumption intensity and carbon emission intensity should drop by 17% and 18% from 2010, the emissions of COD, ammonia nitrogen, SO2, and NOX should drop by 12%, 13.3%, 14.9%, and 16.1%, respectively. GDP, energy consumption intensity, carbon emission intensity, and the emissions of COD, ammonia nitrogen, SO2, and NOX in 2010 are 2855.321 billion Yuan, 56.85 tce per million Yuan, 166.62 tons per million Yuan, 295127.6 t, 15441.3 t, 1382874 t, and 1167190 t, respectively, so the targets of above indicators in 2015 are 4393.265 billion Yuan, 46.62 tce per million Yuan, 136.62 tons per million Yuan, 259712.3 t, 13387.6 t, 1176826 t, and 979272 t, respectively. Assuming the annual average growth rate of fixed assets investment fluctuates from 12% to 18%, the lower and upper bounds of the fixed assets investment of three industries are 380.13 and 400.49 billion Yuan.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 CO2 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:(24)min⁡P1d1++d2++d3-+d4-,P2d5+,P3d6++d7++d8++d9++d10+,G=∑i=13Ii×xi,CrI1×x1G-0.07≤d1+≥α1,CrI2×x2G-0.48≤d2+≥α2,Cr0.45-I3×x3G≤d3-≥α3,Cr43932.65-G×10-5≤d4-≥α4,Crx4G-0.4662≤d5+≥β,CrP1×I2×x2-E1×x5-295127.6×10-6≤d6+≥γ1,CrP2×I2×x2-E2×x5-15441.3×10-5≤d7+≥γ2,CrP3×I2×x2-E3×x5-1382874×10-7≤d8+≥γ3,CrP4×I2×x2-E4×x5-1167190×10-7≤d9+≥γ4,Crc×x4G-1.3662≤d10+≥γ5,xi-≤xi≤xi+,i=1,2,…,5,3801.3≤∑i=13xi≤4004.9d1+,d2+,d3-,d4-,d5+,d6+,d7+,d8+,d9+,d10+≥0,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.


### 4.3. Solution and Discussions

According to the above hybrid intelligent algorithm, the population sizem=50, evaluation function parameter a=0.05, crossover probability pj=0.3, mutation probability pb=0.2, and generation cycles T=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.49 tce per million Yuan, the emissions of SO2 and NOX are 1566900 t and 979302 t, and CO2 emission intensity is 139.22 tons per million Yuan.Table 9
Solution results of FCCGP model of Shandong province in economic priority scenario.
Variable Solution x 1 ∗ 130.7 x 2 ∗ 1784.5 x 3 ∗ 2008 x 4 ∗ 37470.8 x 5 ∗ 879.9 d 1 + 0 d 2 + 0 d 3 - 0 d 4 - 0 d 5 + 0.0787 d 6 + 0 d 7 + 0 d 8 + 390074.21 d 9 + 29.57 d 10 + 0.0262In 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.Table 10
Solution results of FCCGP model of Shandong province in environment priority scenario.
Variable Solution x 1 ∗ 130.7 x 2 ∗ 1716.2 x 3 ∗ 2008 x 4 ∗ 37080.1 x 5 ∗ 1307.6 d 1 + 0.0101 d 2 + 0.0204 d 3 - 0.0103 d 4 - 610.71 d 5 + 0 d 6 + 0 d 7 + 0 d 8 + 0 d 9 + 0 d 10 + 0A scenario endowing the equivalent status for three subsystems is put forward to compare with the first two scenarios. The optimal solutions are shown in Table11. 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 SO2 are 1340146 t.Table 11
Solution results of FCCGP model of Shandong province in the scenario of equivalent status.
Variable Solution x 1 ∗ 130.7 x 2 ∗ 1749.6 x 3 ∗ 2008 x 4 ∗ 44278.8 x 5 ∗ 1103.5 d 1 + 0.0007 d 2 + 0.0186 d 3 - 0.0293 d 4 - 549.72 d 5 + 0.0949 d 6 + 0 d 7 + 0 d 8 + 163320 d 9 + 0 d 10 + 0It 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 of 3Es system 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.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.


## 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 Table1. All the indicators show the upward trend; therefore, regression, grey forecasting, and other forecasting methods can be used to forecast the future values theoretically.Table 1
Data of the fixed assets investments, energy consumptions, and investment in environmental pollution control of Shandong province from 2000 to 2014.
Year Fixed assets investment (108 Yuan) Total energy consumption (104 tce) Investment in environmental pollution control (108 Yuan) Primary industry Secondary industry Tertiary industry 2000 77.1 1176.7 1288.8 11361.9 76.2 2001 78.2 1193.1 1306.9 9955.0 89.8 2002 79.1 1206.3 1321.3 14599.0 144.6 2003 81.4 1241.2 1359.6 16624.8 156.4 2004 87.4 1333.1 1460.2 19623.7 191.9 2005 89.9 1371.1 1501.8 24162.0 238.8 2006 91.5 1395.8 1528.8 26759.3 258.1 2007 95.2 1452.2 1590.6 29176.6 320.8 2008 102.5 1564.0 1713.1 30570.0 432.2 2009 99.4 1515.5 1660.0 32420.2 459.5 2010 103.0 1570.7 1720.4 34807.8 483.9 2011 110.0 1678.1 1838.1 37132.0 614.1 2012 110.9 1692.2 1853.5 38899.3 739.1 2013 111.4 1698.8 1860.7 35358.0 848.0 2014 111.7 1703.9 1866.3 36511.0 823.8 Sources:Shandong Statistics Yearbook 2001–2015 and China Energy Statistical Yearbook 2001–2015. Note: — indicates the lack of data.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:(18)FI1=2.8271t+72.629560.781821.5116R2=0.9727,F=426.7498,sig.=1.51e-11.Note: data in the parentheses denotet-test values.t-test andF value show that the regression equation is significant and can be used for short-term forecasting. Inputting t=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 vectorX1=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:(19)dx1dt+ax1=b,where a and b can be derived by using least squares estimation, and the calculation formulas are as follows:(20)ab=BTB-1BTyn,B=-12x11+x211-12x21+x311⋮⋮-12xn-11+xn11,yn=x20x30⋮xn0.The results of a and b are −0.0295 and 75.6988, and the grey forecast model is(21)x^t+11=x10-bae-at+ba=2641.9629e0.0295t+75.6988.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.Figure 3
Comparison of grey forecast results and actual values.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. Inputtingt=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 regression (QR), and cubic curve regression (CR) can also be applied to fit the relation between the fixed assets investment in primary industry and the annual sequence number. The regression results are obtained as follows:(22)ER:FI1=74.245e0.0301t,R2=0.9691,PR:FI1=69.661t0.1636,R2=0.8773,QR:FI1=0.0274t2+3.2653t+71.391,R2=0.9736,CR:FI1=-0.0191t3+0.4308t2+0.2372t+76.065,R2=0.9827.All the correlation coefficient square values are greater than 0.8, so the above regression equations are significant. Inputtingt=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 x1 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.Table 2
Forecast results of decision variables in 2015.
Forecast method Fixed assets investment (108 Yuan) Total energy consumption (104 tce) Investment in environmental pollution control (108 Yuan) Primary industry Secondary industry Tertiary industry LR 117.9 1797.7 1969.1 43613.5 851.2 GF 119.6 1824.8 1998.7 45597.1 1112.7 ER 120.2 1833.1 2008.0 51777.5 1211.8 PR 109.6 1672.5 1831.9 40551.8 735.2 QR 130.7 1778.8 1948.3 38334.4 984.6 CR 111.9 1707.3 1870.3 33948.2 970.3 Lower bound 109.6 1672.5 1831.9 33948.2 735.2 Upper bound 130.7 1833.1 2008.0 51777.5 1211.8The added values per unit investment of three industries from 2000 to 2014 are shown in Figure4, and all three curves show rising trends. By using regression and grey forecast methods, the forecast values in 2015 are shown in Table 3.Table 3
Forecast results of the added values per unit investment of three industries in 2015 (unit: Yuan/Yuan).
Forecast method The added values per unit investment Primary industry Secondary industry Tertiary industry LR 20.57 14.77 7.48 GF 20.64 16.68 8.31 ER 20.66 17.31 8.37 PR 19.74 12.72 6.38 QR 20.80 15.83 8.14 CR 21.21 15.62 8.01Figure 4
The added values per unit investment of three industries from 2000 to 2014.The production coefficients of four pollutants and CO2 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, CO2 emission coefficients are forecasted based on the samples of 2005–2014.Table 4
The production coefficients of pollutants and CO2 emission coefficients from 2000 to 2014.
Year The production coefficients of pollutants (t/108 Yuan) CO2 emission coefficients (t/104 Yuan) COD Ammonia nitrogen SO2 N O X 2000 531.08 — 427.73 — 1.5588 2001 473.72 13.63 383.51 — 1.5956 2002 385.00 12.52 340.54 — 1.2807 2003 338.53 10.41 339.03 — 1.5982 2004 280.36 9.58 293.57 — 1.6750 2005 247.16 9.03 284.08 — 2.2665 2006 199.37 13.91 311.22 121.74 2.0849 2007 180.69 7.29 251.87 89.25 1.9663 2008 182.83 8.00 280.38 77.03 1.9407 2009 165.39 9.27 273.67 73.96 1.7856 2010 154.63 11.32 268.44 69.16 1.6662 2011 150.94 7.00 314.13 68.12 1.5793 2012 134.31 5.69 239.98 62.21 1.5230 2013 98.60 4.26 247.97 59.57 1.2420 2014 86.96 3.42 268.11 61.15 1.2194 Note: — indicates the lack of data.Table 5
Forecast results of the production coefficients of four pollutants and CO2 emissions coefficients in 2015.
Forecast method The production coefficients of pollutants (t/108 Yuan) CO2 emissions coefficient (t/t) COD Ammonia nitrogen SO2 N O X LR 14.33∗ 4.24 225.62 45.22 1.1026 GF 74.65 5.32 237.77 53.96 1.1842 ER 80.78 4.51 234.10 50.72 1.1640 PR 114.70 5.84 250.08 56.27 1.3464 QR 127.19 2.91 282.99 70.03 1.0765 CR 57.94 −0.27∗ 239.55 48.1 0.9863 Note:∗ indicates that the value does not conform to the reality and should be discarded.The removal coefficients of COD, ammonia nitrogen, SO2, and NOX from 2000 to 2014 are shown in Table 6, and the forecast values of them in 2015 are shown in Table 7.Table 6
The removal coefficients of four pollutants from 2000 to 2014 (unit: t/108 Yuan).
Year The removal coefficients of pollutants COD Ammonia nitrogen SO2 N O X 2000 22312.87 — 4204.09 — 2001 19102.87 386.51 4054.03 — 2002 11278.36 248.46 2884.60 — 2003 10841.76 247.79 3612.67 — 2004 9002.71 224.17 3286.14 — 2005 7506.24 195.67 3160.82 — 2006 6529.95 449.66 5691.74 13.98 2007 5667.44 204.30 4286.80 11.91 2008 4968.03 206.52 5140.90 18.51 2009 4824.66 272.07 5948.39 50.10 2010 4789.07 363.43 6514.73 27.20 2011 4409.67 195.86 7002.29 29.71 2012 3601.40 145.61 4683.00 60.08 2013 2523.13 103.74 5034.74 77.42 2014 2498.19 92.75 6540.74 84.55 Note: — indicates the lack of data.Table 7
Forecast results of the removal coefficients of pollutants in 2015 (unit: t/108 Yuan).
Forecast method The removal coefficients of pollutants COD Ammonia nitrogen SO2 N O X LR −1206.00∗ 135.80 6503.00 86.63 GF 1741.53 × × 110.45 ER 2085.52 127.80 6655.57 112.52 PR 3099.84 161.19 5659.68 72.30 QR 4924.04 67.34 6539.95 103.59 CR 72.02∗ −66.19∗ 4731.15 113.16 Note:∗ indicates that the value does not conform to the reality and × indicates the corresponding method does not pass significant test.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 variableT (u1, u2, and u3) is fully determined by the triplet (u1, u2, and u3) of crisp numbers with u1<u2<u3, whose membership function is given by(23)
μ
x
=
x
-
u
1
u
2
-
u
1
,
if
u
1
≤
x
≤
u
2
,
x
-
u
3
u
2
-
u
3
,
if
u
2
≤
x
≤
u
3
,
0
,
otherwise,where u1, u2, and u3 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 Table 8.Table 8
The values of the parameters in triangular fuzzy variables of FCCGP model.
Variable u 1 u 2 u 3 I 1 19.74 20.60 21.21 I 2 12.72 15.49 17.31 I 3 6.38 7.78 8.37 P 1 57.94 91.05 127.19 P 2 2.91 4.56 5.84 P 3 225.62 245.02 282.99 P 4 45.22 54.05 70.03 c 0.9863 1.1433 1.3464 E 1 1741.53 2962.73 4924.04 E 2 67.34 123.03 161.19 E 3 4731.15 6017.87 6655.57 E 4 72.30 99.78 113.16


## 4.2. FCCGP Model

According to the 12th five-year plan of Shandong province, the annual average growth rate of GDP is 9%, the industrial structure in 2015 is 7: 48: 45, the annual average growth rate of fixed assets investment is set about 15%, energy consumption intensity and carbon emission intensity should drop by 17% and 18% from 2010, the emissions of COD, ammonia nitrogen, SO2, and NOX should drop by 12%, 13.3%, 14.9%, and 16.1%, respectively. GDP, energy consumption intensity, carbon emission intensity, and the emissions of COD, ammonia nitrogen, SO2, and NOX in 2010 are 2855.321 billion Yuan, 56.85 tce per million Yuan, 166.62 tons per million Yuan, 295127.6 t, 15441.3 t, 1382874 t, and 1167190 t, respectively, so the targets of above indicators in 2015 are 4393.265 billion Yuan, 46.62 tce per million Yuan, 136.62 tons per million Yuan, 259712.3 t, 13387.6 t, 1176826 t, and 979272 t, respectively. Assuming the annual average growth rate of fixed assets investment fluctuates from 12% to 18%, the lower and upper bounds of the fixed assets investment of three industries are 380.13 and 400.49 billion Yuan.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 CO2 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:(24)min⁡P1d1++d2++d3-+d4-,P2d5+,P3d6++d7++d8++d9++d10+,G=∑i=13Ii×xi,CrI1×x1G-0.07≤d1+≥α1,CrI2×x2G-0.48≤d2+≥α2,Cr0.45-I3×x3G≤d3-≥α3,Cr43932.65-G×10-5≤d4-≥α4,Crx4G-0.4662≤d5+≥β,CrP1×I2×x2-E1×x5-295127.6×10-6≤d6+≥γ1,CrP2×I2×x2-E2×x5-15441.3×10-5≤d7+≥γ2,CrP3×I2×x2-E3×x5-1382874×10-7≤d8+≥γ3,CrP4×I2×x2-E4×x5-1167190×10-7≤d9+≥γ4,Crc×x4G-1.3662≤d10+≥γ5,xi-≤xi≤xi+,i=1,2,…,5,3801.3≤∑i=13xi≤4004.9d1+,d2+,d3-,d4-,d5+,d6+,d7+,d8+,d9+,d10+≥0,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.


## 4.3. Solution and Discussions

According to the above hybrid intelligent algorithm, the population sizem=50, evaluation function parameter a=0.05, crossover probability pj=0.3, mutation probability pb=0.2, and generation cycles T=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.49 tce per million Yuan, the emissions of SO2 and NOX are 1566900 t and 979302 t, and CO2 emission intensity is 139.22 tons per million Yuan.Table 9
Solution results of FCCGP model of Shandong province in economic priority scenario.
Variable Solution x 1 ∗ 130.7 x 2 ∗ 1784.5 x 3 ∗ 2008 x 4 ∗ 37470.8 x 5 ∗ 879.9 d 1 + 0 d 2 + 0 d 3 - 0 d 4 - 0 d 5 + 0.0787 d 6 + 0 d 7 + 0 d 8 + 390074.21 d 9 + 29.57 d 10 + 0.0262In 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.Table 10
Solution results of FCCGP model of Shandong province in environment priority scenario.
Variable Solution x 1 ∗ 130.7 x 2 ∗ 1716.2 x 3 ∗ 2008 x 4 ∗ 37080.1 x 5 ∗ 1307.6 d 1 + 0.0101 d 2 + 0.0204 d 3 - 0.0103 d 4 - 610.71 d 5 + 0 d 6 + 0 d 7 + 0 d 8 + 0 d 9 + 0 d 10 + 0A scenario endowing the equivalent status for three subsystems is put forward to compare with the first two scenarios. The optimal solutions are shown in Table11. 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 SO2 are 1340146 t.Table 11
Solution results of FCCGP model of Shandong province in the scenario of equivalent status.
Variable Solution x 1 ∗ 130.7 x 2 ∗ 1749.6 x 3 ∗ 2008 x 4 ∗ 44278.8 x 5 ∗ 1103.5 d 1 + 0.0007 d 2 + 0.0186 d 3 - 0.0293 d 4 - 549.72 d 5 + 0.0949 d 6 + 0 d 7 + 0 d 8 + 163320 d 9 + 0 d 10 + 0It 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 of 3Es system 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.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.


## 5. 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.



---

*Source: 2675759-2016-08-09.xml*