Characterization of Renewable Energy Utilization Mode for Air-Environmental Quality Improvement through an Inexact Factorial Optimization Approach

Energy-related environmental problems have been hot spot issues in regional energy system sustainable development. Thus, comprehensive planning of energy systems management is important for social and economic development, as well as environmental sustainability. In addition, uncertainties and complexities, as well as their potential interactions pose a great challenge for effective management in energy and environmental system. This study proposes a stochastic factorial energy systems management model to conduct uncertainties and risks in the energy systems, as well as handle their interaction effects among different environmental policies. The developed method can not only tackle uncertainties expressed as probability distributions and even interval values, but also be applied to determine decision alternatives associated with multiple economic penalties if the formulated environmental policy targets are violated. Meanwhile, by introducing the factorial technology, it can analyze a parameter’s impact on the system and their coordination effect. To verify the feasibility and effectiveness of the proposed method, the developed model was applied to a hypothetical case study for energy structure optimization under considering energy supply, SO2 emissions reduction, and environmental quality requirements. Multiple facilities, related environmental pollutants, and energy demand levels were taken into account. Moreover, the key factors of the system and their interaction effect were discovered. The results indicated that the developed method can resolve meritorious uncertainties in decision-making and analysis, generate effective management programming under multi-levels of the proposed energy and environmental systems. The method can be used for supporting the adjustment for allocating fossil fuels and renewable energy resources, analyzing the tradeoff between conflicting economic and environmental objectives and formulating the local policies.


Introduction
Along with the rapid development of economic, the energy demands are also increasingly growing.According to the International Energy Agency, the world energy consumption has more than doubled in the past thirty years, and it is estimated that the consumption will continue to rise by 50% over the next thirty years.The massive fossil fuel consumption produces adverse effects, including climate change, regional air-quality deterioration, and energy resource depletion.Thus, effective planning of energy systems management is important for accelerating social and economic development, as well as environmental sustainability.Comprehensive management is vital for reducing system costs, carbon emissions and risks.It is essential to keep economies growing while preventing catastrophic effects.However, the energy management processes and the related factors contain multiple uncertainties, such as distributing energy demand, planning power generation, or dealing emission reduction [1][2][3][4].In addition, the energy system contains numerous factors.It is also important to determine the key factors affecting the system and analyzing their potential interactions.Therefore, to tackle the existing challenge, it is required to develop effective decision-supporting tools for adjusting energy systems management.
Previously, there were research efforts focusing on a variety of complexities and uncertainties in energy and environmental systems management [5][6][7][8][9][10][11].Among them, stochastic programming methods are widely devised and analyzed to provide reliable assistance for decision making.In general, there are two kinds of stochastic programming approaches, which are two/multi-stage programming (T/MSP) and chance constrained programming (CCP) approaches.The T/MSP is effective for policy related problems, where coefficients in variables are provided with given probability distributions or it could be easily estimated [12][13][14][15].The CCP is most used for risk based random issues.In this case, system constraints are not required to be totally satisfied, but it should meet the requirements not less than the provided confidence level [16,17].Thus, the method combination has an advantage in reflecting complexities of system uncertainties, and it is well suited for policy analysis with pre-setting targets which might be violated during the actual operation.For example, Fleten and Kristofferson [18] developed short-term planning to demonstrate a hydropower plant operating system, where a stochastic programming method was introduced to address the uncertain beyond the operation day to maintain a balance between current profits and expected future benefits.Rong and Lahdelma [19] established a carbon emission trading optimization model with stochastic methods.The scheme shows satisfactory transaction efficiency in terms of profit turnover rate.Wang et al. [20] optimal CO 2 trading model with stochastic programming under uncertainty.Desired policies could be identified by analyzing decision alternatives.Cai et al. [21] proposed an inexact method based on community-scale energy management systems to balance the tradeoffs between system costs and violation risks in the constraint.Cheng et al. [22] discussed how energy structure optimization process could reduce air pollution in China.
However, in the energy management system, it always contains multiple parameters and components with multiple periods, multiple objectives and dynamic constraints.Many system parameters and components may be associated with different formats of uncertainties.Particularly, these uncertainties may be correlated with each other, leading to significant impacts on the resulting energy management strategies and the associated risks.The stochastic optimization methods provide a powerful scheme to address the uncertain random information.But it fails to distinguish the different degrees of importance and their relationships.It is thus advisable to develop modified tools for identifying the main factors and their potential interactions, and also reveal their impact on the system objectives.
Therefore, a stochastic factorial energy systems management (SFESM) model will be proposed for optimal energy systems management under considering environmental pollutants emission reduction, energy resources consumption control, and multiple uncertainties.The SFESM model will be able to reflect interactions among multiple uncertainties in energy and environmental management systems.In detail, two-stage stochastic programming will be employed to reflect the conflicts between developing an economy and protecting the environment; method of inexact chance-constrained programming will be used for tackling uncertainties in the system's capacity limitations.Moreover, a factorial analysis will be undertaken to locate the main effects of crucial factors and their interactions to the system objective, since factorial technique has been recommended as a useful tool for sensitive analysis [23][24][25][26][27][28][29][30][31].Through the introduction of factorial analysis to energy environmental systems research, the reflection of high-dimension complexities that are beyond the coverage of regular energy management models with the related environmental problems would become possible.The proposed SFESM model will then be applied to a case study of energy systems management for improving air-environmental quality to demonstrate its applicability.The following optimization scheme will be obtained: (I) energy purchasing plan, (II) electricity generation scheduling and (III) contamination control strategy.The results can be helpful for identifying energy allocation patterns, addressing conflicts between economic objectives and environment, as well as examining the linkage between existing policies and economic penalties.

Methodology
Optimization techniques are considered as effective tools to identify suitable strategies for management with various complexities.On the basis of traditional deterministic models, uncertain optimization programming is widely used for supporting sustainable energy system planning.

Two-Stage Programming
A general linear programming model can be formulated as follows: subject to: x ≥ 0, where f is the objective function, x is the decision variable, a, b and c are the parameters in the objective and constraints.
In the real-world applications, the biggest weakness for the linear programming model is that most variables and parameters are uncertain and cannot be expressed as a definite number [32,33].That means in instance study, a, b and c are random variables all along, decisions should be made at the proper discrete points or probability levels.TSP method is efficacious to handle uncertain dada, or an analysis of multiple scenarios are desired.In TSP, the first-stage planning is made with historical experience.After a random event, the remedial measure could be taken at the second-stage to deal with the "penalties" [34].In general, a typical TSP model can be formulated as follows: subject to: where p m is the probability of occurrence for scenario m, with p m ≤ 1 and M m=1 p m = 1.In this model, decision variables X t must be determined at the first stage.Then, the correction variable W m,t can be introduced after randomness is revealed [35].

Chance-Constrained Programming
In the management of energy systems, trade-offs exist if the manager intends to pursue the maximization profit.To meet the desired energy supply, it would require the reduction of industrial development scale.But also, taking the restrictive conservation measures over production would retard their participation in energy programs, which may lead to a decrease in economic benefits.Therefore, in the economic, environmental or policy aspect, changes provided as tolerance levels would be useful for analyzing risk-based information.The CCP method is effective for processing risk analysis.For the objective function, it could provide trade-off analysis; for the constraint, it is helpful to identify the tolerance values.The method is mostly applied to cases where the prescribed level of probability could be formulated.
A general stochastic linear programming problem can be formulated as follow: subject to: where X is a vector of decision variables, and A(t), B(t) and C(t) are sets with random elements defined on a probability space T [36,37].To solve this model, it should be converted to a deterministic vision by providing satisfaction degree.By introducing probability p i ∈ [0, 1], for each constraint i, the relevant satisfaction degree determined by decision makers should not less than 1 − p i .It can be expressed as follows: which are generally nonlinear.For feasible constraints, the distributions should be convex only with certain levels of p i .For example, (I) a ij are deterministic and b i are random (for all p i values), (II) a ij and b i are discrete random coefficients, with p i ≥ max r=1,2,...,R (1 − q r ), where q r is the probability associated with realization r, or (III) a ij and b i have Gaussian distributions, with p i ≥ 0.5 [38].When a ij are deterministic and b i are random, constraint (12) becomes linear: where b i (t) (p i ) = F −1 i (p i ), given the cumulative distribution function of b i , and the probability of violating constraint i.The problem can only reflect the case when A is deterministic [39][40][41][42].Thus, a typical CCP model can then be reformulated as follows: subject to:

Inexact Chance Constrained Two-Stage Stochastic Programming
In energy systems management model, the TSP and CCP methods are valid for tackling right-hand-side uncertainties such as energy resources availabilities that are presented as probability distributions.However, it is occasionally difficult to acquire professional information that meets the demand.Interval parameter programming (IPP) is valid for processing complexities that cannot be expressed with probability distributions.And they can be both in the objective functions and constraints.
Normally, an interval number x ± can be expressed as [x − , x + ], where x + represents the maximum value while x − is for the minimum value.When x − = x + , x ± becomes a deterministic number.IPP can be proposed to the optimization framework to handle uncertainties which cannot be expressed by distribution functions, but with given lower and upper bounds.A general IPP model can be defined as follows: max subject to: x ± ≥ 0 (19) where , R ± is a matrix composed by interval number.Integrated with TSP, CCP and IPP methods into a general framework, an inexact chance-constrained two-stage stochastic programming (ICCTSP) model could be formulated as follows: subject to: Over all, the ICCTSP method can deal with optimization management associated with multiple uncertainties.It can also handle the complex tradeoff between the conflicting objects and variables within a two-stage context.
To solve this model, according to Huang et al. [43], the uncertain ICCTSP model should be partitioned into two deterministic sub-models f + and f − , with each corresponding to the upper and lower of system objective, respectively.Through solving the corresponding sub-models, the optimal solution for the system can be expressed as: Main haze for the computational process is that it is hard to know whether sub-model f + or f − corresponds to the minimum object value [44].An intermediate variable y t is thus designed by letting X ± t = X − t + ∆X t y t , where ∆X t = X + t − X − t and 0 ≤ y t ≤ 1.Therefore, the lower bond of the objective function can firstly be formulated as follows: Pr where W − m,t and y i are decision variables, and their solutions of W − m,topt , y topt and f − opt can be obtained.Based on the providing method, the upper bound of the objective function value can be formulated as follows: Pr where W + m,t are decision variables and carrying out corresponding solutions of W + m,topt and f + opt .The optimized contaminants generation target can be determined by calculating X ± i = X − t + ∆X t y t .

Factorial Analysis
The proposed optimization method is effective in addressing uncertainties that exist in system parameters.However, in energy systems management, it is also important to reveal the contributions of individual uncertain inputs that are related to various economic and environment conditions.Therefore, to obtain more favorable decision support, it is important to address sensitivity analysis for investigating the uncertain inputs in the in the optimization process.Compared to traditional techniques, a sensitivity analysis using factorial methods are efficient for identifying the effects of two or more parameters.Based on a factorial analysis, for each complete trial of the repetition process, all possible combinations for two or multiple levels factors could be examined.A factor effect can be estimated at different levels of other factors, and the valid conclusions over a series of experimental conditions could be provided.In addition, to avoid misleading conclusions, a factorial analysis is necessary when factor interactions may exist [30,31].
The most common factorial analysis is based on the 2 k factorial design.A typical 2 k factorial design would contain k design factors, with each being expressed as interval numbers.Besides, interaction effects exist between factors, since for each factor, the effect is depending on the level chosen for other factors.In a full factorial design, all possible combinations of the factors with different levels should be investigated [45].The main effect on individual factor represents the difference between the average response at the low and high levels.The interaction effect could be indistinctive in some experiments.In some cases, there would be massive factors to be investigated, parameter selection should be processed to identify the key factors.For less significant factors, they could be eliminated to reduce unnecessary computation.
Therefore, combined the factorial analysis and the proposed uncertain method within the optimization framework, a stochastic factorial method will be developed.It can track uncertainties particularly described as probability distributions or interval values; it is valid for analysis where policy scenarios are desired.Moreover, it could reveal the contributions of uncertain variables to the system objective and their potential interactions.

Overview of the Study System
The renewable energy management system is intricate, and it consists of multiple factors such as economic, environmental and mechanical manners.The system is the combination of resource supply, power generation, capacity expansion and energy utilization process.And more importantly, those internal factors have definite connections with each other through transmission lines.Generally, primary energy supply resources include fossil fuels and renewable energy resources.Correspondingly, a variety of power generation technologies should be adopted in the energy management systems.Energy power generated should be allocated to meet the need of end users.Besides, according to the requirements of social and economic development, it is also very important to ensure all the facilities meet the environmental standard.
In this study, a hypothetical, but typical renewable energy management model is proposed for demonstrating the capability of the proposed optimization approach.The general structure of the proposed system is presented in Figure 1.The data are summarized through government reports and related literature, typical cost, and representative technical data are provided.The system contains multiple energy sources and technologies.Particularly, coal, diesel and natural gas are selected as typical fossil energy, while wind power and solar power are chosen as representative renewable energy.Decision makers should develop an effective plan to satisfy the industrial, agricultural and municipal energy demands.Regarding the environmental aspect, the power generation process will cause an inevitable pollution emission.SO 2 is chosen as the representative emission.Based on historical data, a predefined pollution control target is stipulated to each power generation facility.If the promised target is achieved, the unit could obtain certain profits as the net benefit.However, if an excess emission occurs, an economic penalty will be charged.
In this study, a hypothetical, but typical renewable energy management model is proposed for demonstrating the capability of the proposed optimization approach.The general structure of the proposed system is presented in Figure 1.The data are summarized through government reports and related literature, typical cost, and representative technical data are provided.The system contains multiple energy sources and technologies.Particularly, coal, diesel and natural gas are selected as typical fossil energy, while wind power and solar power are chosen as representative renewable energy.Decision makers should develop an effective plan to satisfy the industrial, agricultural and municipal energy demands.Regarding the environmental aspect, the power generation process will cause an inevitable pollution emission.SO2 is chosen as the representative emission.Based on historical data, a predefined pollution control target is stipulated to each power generation facility.If the promised target is achieved, the unit could obtain certain profits as the net benefit.However, if an excess emission occurs, an economic penalty will be charged.
This study aims to develop an optimization integrate sensitive analysis model for supporting the scientific management of energy system by effectively dealing with the following problems: (1) Various types of uncertainties existed in the provided system, it should be reflected accurately and resolved effectively; (2) An optimal solution of energy activities should be obtained to balance the tradeoff between economic and environmental benefit, as well as the system stability.The results should include energy generation plan, technology selection result, and capacity expansion scheme; (3) The environmental control process should be represented.The related environmental factors would restrict the application of energy resources and technologies.For example, with a stricter policy, energy with lower emission rate or higher generation efficiency will be recommended.Correspondingly, it will provoke increased system cost; (4) The key impact factors and their potential interaction should be revealed.It would provide more detailed information for future decision making.

Stochastic Factorial Energy Systems Management Model
Therefore, the objectives of the study system are: (1) Arrange the five power conversion technologies effectively to meet the energy demands while maximizing the system benefits under uncertainty; (2) determine the minimum energy supply amount; (3) coordinate energy activities and environmental standard.Through the proposed method, a stochastic factorial energy systmes management model can be formulated.The objective function of the SFESM model can be expressed as follows: This study aims to develop an optimization integrate sensitive analysis model for supporting the scientific management of energy system by effectively dealing with the following problems: (1) Various types of uncertainties existed in the provided system, it should be reflected accurately and resolved effectively; (2) An optimal solution of energy activities should be obtained to balance the tradeoff between economic and environmental benefit, as well as the system stability.The results should include energy generation plan, technology selection result, and capacity expansion scheme; (3) The environmental control process should be represented.The related environmental factors would restrict the application of energy resources and technologies.For example, with a stricter policy, energy with lower emission rate or higher generation efficiency will be recommended.Correspondingly, it will provoke increased system cost; (4) The key impact factors and their potential interaction should be revealed.It would provide more detailed information for future decision making.

Stochastic Factorial Energy Systems Management Model
Therefore, the objectives of the study system are: (1) Arrange the five power conversion technologies effectively to meet the energy demands while maximizing the system benefits under uncertainty; (2) determine the minimum energy supply amount; (3) coordinate energy activities and environmental standard.Through the proposed method, a stochastic factorial energy systmes management model can be formulated.The objective function of the SFESM model can be expressed as follows: max (1) System benefit related to environmental policy: (2) Cost for power generation (3) Purchasing cost for primary energy supply (4) Cost for capacity expansion The net benefit is to be maximized by a series of constraints.The impact factors and their interactions are also provided.The detailed constraints are expressed as follows: (1) Mass balance (energy supply should be more than resource consumption) (2) Availability of energy resources (energy usage should be less than its availability) (3) Electricity constraints (electricity generation should be more than energy demands) (4) Capacity limit (installed capacity should be more than electricity generation) (5) Constrains for controlling contamination (SO 2 emission should be less than environmental standard) (Actual SO 2 emission should be less than emission control target) (SO 2 generation amounts should be more than the excess amount, and less than given maximum) 1, if expansion with option n for generating technology is undertakan 0, otherwise , The detailed nomenclatures for the variables and parameters are provided in the Appendix A. The research target is to maximize the system benefit under uncertainty.The developed energy model can be solved through the above method by decomposed into two deterministic sub-models.Specifically, the two-stage problem can be solved by letting , where ∆WS i,t = WS + i,t − WS − i,t .According to the Formulas ( 14)-( 16) in Section 2.2, probability bounds of constraint violation under consideration are 1% and 10%.Transform the developed model into two sub-models, formulating the first sub-model which corresponds to f − opt .Then Y − opt , X − opt , ZY − opt can be obtained through the solution.Similarly, formulating the second sub-model corresponding to f + opt and substituted calculated value above.Combine the two sub-models' solutions to obtain the optimal solution of the proposed model.The data of energy activities are provided in Table 1.Pollution treatment targets and the related economic data are given in Table 2.
Furthermore, the factorial technique is introduced to a sensitive analysis of factor impact for system objective.It is predetermined that if energy activities meet the environmental standards, the specific regulation would affect net benefit directly.Therefore, factors related to pollution emission would have a more obvious influence on the evaluation of the system.The factorial design would then be simplified.It is not necessary to examine all system variables, only the factors relevant to the contamination control process should be taken into account.It would reduce the complexity of sensitive analysis and retrench the operation time.By factorial analysis technique, the key impact factors could be revealed.Meanwhile, it could demonstrate their potential interaction effect, which may have a greater impact on the system than individual variable.It would provide more detailed information for decision maker than traditional energy management methods.Future planning should be paid more attention to the key factors such as formulating stricter environmental standard or greater penalties.

Result Analysis
In this study, fifteen planning periods are considered, and each representative for one year.Through solving the developed model, the optimal primary energy supply scheme, electricity generation plan and conceivable capacity expansion options were generated.Solutions of contamination control and excess emission amount were calculated under different risk levels.It can demonstrate a basic tendency of energy activities and provide potential suggestions for policy development.
Figure 2 presents the purchasing plan for coal at each planning period.The result shows that the allocation amount had a visible downtrend over time.The allocated amount would decrease from [585.2, 631.4] × 10 3 tonnes to [559.3, 599.7] × 10 3 tonnes, fell about 4~5% during the fifteen years.It is in conformity with the basic requirement of a popular total amount control policy for coal.On the other hand, it also reveals that coal would still be the major energy source in the near future for its splendid reserve, extensive distribution, and shallow embedding.It means the coal usage should be on a declining trend yet would be less obvious.The non-renewable energy resources are a major component of energy supply in most cases.In this study, diesel and natural gas would also be selected to provide power to the end users.Figure 3 provides the energy supply pattern of diesel and natural gas.The diesel supply would increase from 225.3 to 238.4 × 10 3 tonnes during the planning period, while the deliverability of natural gas would increase from 194.8 to 225.4 × 10 3 tonnes.The growth rate would be 5.8% and 15.7% respectively.The total purchasing amount of diesel and natural gas would be nearly half of coal.The allocated amount for diesel would be relatively stable.However, there's a significant increase in natural gas.The difference between these two energy sources would be diminished over the planning period.If the tendency persists, natural gas may take coal's place as the most favorite fuel in just a few decades.In a practical environment, according to BP statistical review of world energy, the world's energy structure begins to diversify.It trends to increase thermal power generation by natural gas and decrease oil usage.That indicated the proposed case study could reflect the general trend thus provide useful information for policy making.
while the deliverability of natural gas would increase from 194.8 to 225.4 × 10 tonnes.The growth rate would be 5.8% and 15.7% respectively.The total purchasing amount of diesel and natural gas would be nearly half of coal.The allocated amount for diesel would be relatively stable.However, there's a significant increase in natural gas.The difference between these two energy sources would be diminished over the planning period.If the tendency persists, natural gas may take coal's place as the most favorite fuel in just a few decades.In a practical environment, according to BP statistical review of world energy, the world's energy structure begins to diversify.It trends to increase thermal power generation by natural gas and decrease oil usage.That indicated the proposed case study could reflect the general trend thus provide useful information for policy making.With the economy development, electrical power demand would increase every year.On the one hand, coal-fired power is the key factor in electricity generation activities.The efficiency for oilbased generators is close to coal, while the efficiency for natural gas-based generators would be almost twice higher.On the other hand, when emission control is considered, gas may be at least 25% cleaner than coal.From the environmental aspect, natural gas would be the optimal sources for energy supply.However, the operating cost for gas generators is much higher.Thus, it's important to find a balance between them.Figure 4 presents the optimized solution of electricity generation plans of fossil energy.In detail, the energy generation by coal would be [1350.0,1637.2]GWh over the planning horizon while the energy consumption would slightly decrease.Diesel and gas-fired power would be markedly increased.The electricity generated by diesel would range from 780.4 to 842.5 GWh.Gas based electricity would increase from 847.0 GWh to 979.9 GWh, and the median would be 900.7 GWh.The growth rate for gas-fired power would be 15.7%, while the growth for diesel was modest.In this study, renewable energy sources were also employed.In this study, renewable energy sources were also considered.The installed capacities for wind/solar power were set under 50 MW, since the large-scale applications of renewable energy are still not mature and the maintenance costs are much higher (Table 1).The initially installed capacities for fossil energy were set to 300 MW.It turns out that contribution rate for renewable energy would be 1%.Among them, wind power would occupy 58%, while solar energy would be 42%.The results of power generation pattern are provided in Figure 5 The result also shows that they could meet the demands of endusers; hence capacity expansion would not be needed.With the economy development, electrical power demand would increase every year.On the one hand, coal-fired power is the key factor in electricity generation activities.The efficiency for oil-based generators is close to coal, while the efficiency for natural gas-based generators would be almost twice higher.On the other hand, when emission control is considered, gas may be at least 25% cleaner than coal.From the environmental aspect, natural gas would be the optimal sources for energy supply.However, the operating cost for gas generators is much higher.Thus, it's important to find a balance between them.Figure 4 presents the optimized solution of electricity generation plans of fossil energy.In detail, the energy generation by coal would be [1350.0,1637.2]GWh over the planning horizon while the energy consumption would slightly decrease.Diesel and gas-fired power would be markedly increased.The electricity generated by diesel would range from 780.4 to 842.5 GWh.Gas based electricity would increase from 847.0 GWh to 979.9 GWh, and the median would be 900.7 GWh.The growth rate for gas-fired power would be 15.7%, while the growth for diesel was modest.In this study, renewable energy sources were also employed.In this study, renewable energy sources were also considered.The installed capacities for wind/solar power were set under 50 MW, since the large-scale applications of renewable energy are still not mature and the maintenance costs are much higher (Table 1).The initially installed capacities for fossil energy were set to 300 MW.It turns out that contribution rate for renewable energy would be 1%.Among them, wind power would occupy 58%, while solar energy would be 42%.The results of power generation pattern are provided in Figure 5 The result also shows that they could meet the demands of end-users; hence capacity expansion would not be needed.would be 900.7 GWh.The growth rate for gas-fired power would be 15.7%, while the growth for diesel was modest.In this study, renewable energy sources were also employed.In this study, renewable energy sources were also considered.The installed capacities for wind/solar power were set under 50 MW, since the large-scale applications of renewable energy are still not mature and the maintenance costs are much higher (Table 1).The initially installed capacities for fossil energy were set to 300 MW.It turns out that contribution rate for renewable energy would be 1%.Among them, wind power would occupy 58%, while solar energy would be 42%.The results of power generation pattern are provided in Figure 5 The result also shows that they could meet the demands of endusers; hence capacity expansion would not be needed.Normally, a predefined SO2 treatment target is promised to meet the environmental standards.If the commitment amounts are not delivered, due to the insufficient pollution control availabilities, excess emission will occur thus cause penalty for the system benefit.Under normal circumstances, the pre-regulated amounts will meet the demand.The higher probability is formulated for this case.However, some artificial subjective reasons like budgeting control or some uncontrollable factors, pollution treatment facilities are not effectively operated.It could affect pollution treatment availabilities.Therefore, a relatively lower probability is provided.In optimization studies, constraints could be partially satisfied.In order to obtain a more acceptable result, constraint violation may be allowed at certain risk levels.In this study, two risk levels were set to 0.01 and 0.1, respectively.Correspondingly, each represents the satisfaction degree for the constraints should be at least 99% or 90%.The detailed information is provided in Table 3. Solutions were generated for pollution control recommendations.As shown in Table 4, optimized target treatment amounts were firstly obtained for three sources.The solutions of WSiopt revealed that the optimal treatment targets would be 5.0 × 10 3 tonnes for SO2 generated by the burning of coal, 2.8 × 10 3 tonnes for diesel-fired SO2, and 180 tonnes for gas-fired SO2, respectively.The existing Normally, a predefined SO 2 treatment target is promised to meet the environmental standards.If the commitment amounts are not delivered, due to the insufficient pollution control availabilities, excess emission will occur thus cause penalty for the system benefit.Under normal circumstances, the pre-regulated amounts will meet the demand.The higher probability is formulated for this case.However, some artificial subjective reasons like budgeting control or some uncontrollable factors, pollution treatment facilities are not effectively operated.It could affect pollution treatment availabilities.Therefore, a relatively lower probability is provided.In optimization studies, constraints could be partially satisfied.In order to obtain a more acceptable result, constraint violation may be allowed at certain risk levels.In this study, two risk levels were set to 0.01 and 0.1, respectively.Correspondingly, each represents the satisfaction degree for the constraints should be at least 99% or 90%.The detailed information is provided in Table 3. Solutions were generated for pollution control recommendations.As shown in Table 4, optimized target treatment amounts were firstly obtained for three sources.The solutions of WS iopt revealed that the optimal treatment targets would be 5.0 × 10 3 tonnes for SO 2 generated by the burning of coal, 2.8 × 10 3 tonnes for diesel-fired SO 2 , and 180 tonnes for gas-fired SO 2, respectively.The existing SO 2 processing capacity could not meet requirements, excess emission occurred at both risk levels.For risk level at 0.01, when the availability is relatively low with a probability of 40%, excess emission for the three sources would be 5.0, 0.3 and 0.2 × 10 3 tonnes, respectively; when the availability is higher with a probability of 60%, excess SO 2 emission for coal and natural gas would be 1.9 and 0.2 × 10 3 tonnes, while no excess emission would occur for diesel-based power plant.For risk level at 0.1, the availability increases under this circumstance.Generally, with stricter environmental policies, excess emission would decrease.When at the lower probability, excess emission for diesel-fired plants would drop to 0.1 × 10 3 tonnes; while at the higher probability, excess emission for coal-based plants would drop to 1.7 × 10 3 tonnes.Other data would be unchanged at two risk levels.The results also indicate that when pollution treatment availability is insufficient, the SO 2 produced by diesel plants would be firstly treated, and secondly for coal-fired plants.This is because the penalty for diesel-fired pollution is the highest and the relevant units could also bring substantial benefit.For gas-fired plants, SO 2 discharge would not be controlled, since the total amount is quite small.The actual treated amount for each type of power plant at different risk levels can be calculated from treatment targets and excess emissions, all the results are also expressed in Table 4.

Probability (q m ) Source
Coal-Fired Diesel-Fired Gas-Fired To address uncertainties in a thorough manner and provide a comprehensive outcome, a further sensitive analysis is necessary.The object of the proposed system is to address the maximum net benefits while meeting the constraints.It is important to reveal the main factors related to the study objectives.Factorial analysis is therefore introduced to examine these uncertain parameters and their interactions.In this study, the system benefit is calculated based on emission production amount, while the penalty will be taken when violated the environmental policy.Thus, seven factors related to the pollution control process are selected.These factors are donated as A, B, C, D, E, F and G.All the factors are presented at two levels.For factors A: F, they refer to the lower and upper bound of the chosen variable, since they are all interval numbers.For factor G, they represent the two risk levels.Table 5 shows all uncertain parameters.The proposed two-level factorial design with seven factors would require 128 experimental runs.The major purpose of this design would be verifying and evaluating the factors influencing on total system benefits.Figure 6 presents the half-normal plot of the effects of all the factors.The plot is an effective graphical technique to help identify the key factors.The larger the distance represents the more significant impact of the factor.The yellow outcome represents positive effects, which means increasing the net benefits.Correspondingly, the blue square representing penalty to the system.The results indicate that factors A, D, B and G would be the key effect factors.Besides, for some interaction factors, the effect is significant.For example, the effects for interaction DE, DF, DG, CF and CG are higher than the main effect of factor E. It implies that an important interrelationship may exist between the factors and their interaction plot should therefore be analyzed.
Figure 7 provides the interaction plot of factors D and E. It implies that the factor D would have a negative effect on the net benefit.When the value for factor D approach to its upper bound, the total revenue of the system will decrease gradually.This trend would become more prominent when factor E is at its higher bound.Figure 8 presents the interactions plot matrix for factor C/D with F. It indicates that factor C would have a positive effect on the system, but the impact is unapparent.The impact of factor D would become more significant with F at the lower bound, while the impact of C would be more obvious when F at the upper bound.Similarly, Figure 9 shows the interaction plot of factors D and G.It states that the impact of factor D would become more significant when G is at the lower level.The results show that system benefits would be improved with a higher violation risk.When the profits surpass the punishment, it would contribute significantly to excess emission.That indicates that policies regarding the environment will not be loosened in order to allow more profits to pass environmental examinations.It reveals the tradeoffs between conflicting economic development and environment protection.It is a very important responsibility for decision makers to protect the environment during economics developing.
To demonstrate the superiority of the proposed method, a compared result is obtained with a traditional deterministic method.The proposed energy management system contains multiple planning affairs and multiple risk levels, due to the space limitation, only electricity generation plans with risk level α = 0.1 is emphasized in Table 6.The programming approach was solved by replacing the uncertain variables with their mid-point values.Differing from the proposed method, with various uncertain inputs, the application of deterministic programming can only provide a single response.In realistic energy system planning, an accurate value could hardly provide the guidance for decision makers.Similarly, by solving the maximum and minimum values of the uncertain parameters, solutions under best/worst scenarios can be obtained.It can be used for judging the capability of the system, but hardly establish a stable interval for policy or strategy.Besides, the further sensitive analysis could be employed, but solutions of the deterministic method cannot reveal the interaction effects among different factors.Thus, it hardly provides useful analysis for decision variables.Figure 7 provides plot of factors D and E. It implies that the factor D would have a negative effect on the net benefit.When the value for factor D approach to its upper bound, the total revenue of the system will decrease gradually.This trend would become more prominent when factor E is at its higher bound.Figure 8 presents the interactions plot matrix for factor C/D with F. It indicates that factor C would have a positive effect on the system, but the impact is unapparent.The impact of factor D would become more significant with F at the lower bound, while the impact of C would be more obvious when F at the upper bound.Similarly, Figure 9 shows the interaction plot of factors D and G.It states that the impact of factor D would become more significant when G is at the lower level.The results show that system benefits would be improved with a higher violation risk.When the profits surpass the punishment, it would contribute significantly to excess emission.That indicates that policies regarding the environment will not be loosened in order to allow more profits to pass environmental examinations.It reveals the tradeoffs between conflicting economic development and environment protection.It is a very important responsibility for decision makers to protect the environment during economics developing.Figure 7 provides the interaction plot of factors D and E. It implies that the factor D would have a negative effect on the net benefit.When the value for factor D approach to its upper bound, the total revenue of the system will decrease gradually.This trend would become more prominent when factor E is at its higher bound.Figure 8 presents the interactions plot matrix for factor C/D with F. It indicates that factor C would have a positive effect on the system, but the impact is unapparent.The impact of factor D would become more significant with F at the lower bound, while the impact of C would be more obvious when F at the upper bound.Similarly, Figure 9 shows the interaction plot of factors D and G.It states that the impact of factor D would become more significant when G is at the lower level.The results show that system benefits would be improved with a higher violation risk.When the profits surpass the punishment, it would contribute significantly to excess emission.That indicates that policies regarding the environment will not be loosened in order to allow more profits to pass environmental examinations.It reveals the tradeoffs between conflicting economic development and environment protection.It is a very important responsibility for decision makers to protect the environment during economics developing.To demonstrate the superiority of the proposed method, a compared result is obtained with a traditional deterministic method.The proposed energy management system contains multiple planning affairs and multiple risk levels, due to the space limitation, only electricity generation plans with risk level α = 0.1 is emphasized in Table 6.The programming approach was solved by replacing the uncertain variables with their mid-point values.Differing from the proposed method, with various uncertain inputs, the application of deterministic programming can only provide a single response.In realistic energy system planning, an accurate value could hardly provide the guidance for decision makers.Similarly, by solving the maximum and minimum values of the uncertain parameters, solutions under best/worst scenarios can be obtained.It can be used for judging the capability of the system, but hardly establish a stable interval for policy or strategy.Besides, the further sensitive analysis could be employed, but solutions of the deterministic method cannot reveal the interaction effects among different factors.Thus, it hardly provides useful analysis for decision variables.Generally, the above analysis indicates that solutions of the SFESM model can provide an effective relevance with pre-regulated energy policies and the related penalties.It facilitates the settlement of multiple types of uncertainties.Optimal primary energy supply, electricity generation, capacity expansion and pollution control plans are generated.The interval results under different risk levels are operable and can help decision makers obtain diversified decision alternatives.Besides,  Generally, the above analysis indicates that solutions of the SFESM model can provide an effective relevance with pre-regulated energy policies and the related penalties.It facilitates the settlement of multiple types of uncertainties.Optimal primary energy supply, electricity generation, capacity expansion and pollution control plans are generated.The interval results under different risk levels are operable and can help decision makers obtain diversified decision alternatives.Besides, techniques of sensitive analysis can be applied for supporting the further adjustment of the optimization model and promoting the commonality to the practical situation.

Conclusions
In this study, a stochastic factorial programming is proposed for reforming regional energy structure management, and conducting uncertainties and risks, as well as handle their interaction effects among different environmental policies in the energy systems.The proposed method is able to tackle uncertainties expressed as interval values and probability distributions and can be further used for examining all possible decision options which have relevance to various levels of economic penalties if the proposed policy targets are violated.Optimal decisions of emission control schemes, primary energy supply, electricity generation, and capacity expansion can be generated.Particularly, it can help examine uncertain parameters and their interactions to analyze their impact on the system performance through factorial analysis.Compared with the conventional energy and environmental systems management, the proposed method could not only handle the uncertainties expressed as interval and random variables, but also provide more specific results of parameter effects and their potential interactions on the system performance.
The developed method has been applied to a case of environmentally-friendly oriented planning of renewable energy system management system within a multi-facility, multi-period and multi-demand-level context.This study identified significant factors (e.g., environmental control factors) and reflected their interactions in the energy management model.It proves that reducing coal usage through economic measures would be not effective, policies like total quantity control of coal are necessary to reduce SO 2 emissions.The results indicated that the proposed method would incorporate significant uncertain information into the decision-making process and capture a technically feasible solution at different levels of system reliability.The results can be used for supporting the adjustment for

Figure 1 .
Figure 1.Interactive relationships of the energy system.

Figure 1 .
Figure 1.Interactive relationships of the energy system.

Figure 2 .
Figure 2. Energy purchasing plan for coal at each planning period.

Figure 2 . 19 Figure 3 .
Figure 2. Energy purchasing plan for coal at each planning period.Sustainability 2019, 11, x FOR PEER REVIEW 11 of 19

Figure 3 .
Figure 3.Comparison result of energy purchasing plan for diesel and natural gas.

Figure 4 .
Figure 4. Electricity generation plans of fossil energy-based power plant.Figure 4. Electricity generation plans of fossil energy-based power plant.

Figure 4 . 19 Figure 5 .
Figure 4. Electricity generation plans of fossil energy-based power plant.Figure 4. Electricity generation plans of fossil energy-based power plant.Sustainability 2019, 11, x FOR PEER REVIEW 12 of 19

Figure 7 .Figure 6 .
Figure 7. Interaction plot of factors D and E.

Figure 6 .
Figure 6.Half-normal plot of the effects.

Figure 8 .
Figure 8. Interaction plot of factors C/D with F.

Figure 8 .
Figure 8. Interaction plot of factors C/D with F.

Figure 9 .
Figure 9. Interaction plot of factors D and G.

Figure 9 .
Figure 9. Interaction plot of factors D and G.

Table 1 .
Partial economic data of energy activities a .

Table 3 .
Pollution control availabilities under different risk levels (tonne).

Table 3 .
Pollution control availabilities under different risk levels (tonne).

Table 5 .
Investigated factors with two levels.