Optimizing sustainable development problem under uncertainty: Robust vs fuzzy optimization methods

Jia, Ruru; Bai, Xuejie; Song, Fengxuan; Liu, Yankui

doi:10.3233/JIFS-182763

Optimizing sustainable development problem under uncertainty: Robust vs fuzzy optimization methods

Article type: Research Article

Authors: Jia, Ruru^a | Bai, Xuejie^{b; c; *} | Song, Fengxuan^{c; *} | Liu, Yankui^a

Affiliations: [a] Hebei Key Laboratory of Machine Learning and Computational Intelligence, College of Mathematics and Information Science, Hebei University, Baoding, Hebei, China | [b] College of Science, Hebei Agricultural University, Baoding, Hebei, China | [c] College of Management Science, Hebei University, Baoding, Hebei, China

Correspondence: [*] Corresponding author. Xuejie Bai, E-mail: [email protected]; Fengxuan Song, E-mail: [email protected].

Keywords: Sustainable development, multiple objectives, robust optimization, fuzzy optimization

DOI: 10.3233/JIFS-182763

Journal: Journal of Intelligent & Fuzzy Systems, vol. 37, no. 1, pp. 1311-1326, 2019

Published: 09 July 2019

Get PDF

Abstract

Sustainable development problem involves several conflicting criteria on economy, energy, environment and social aspects and some sources of uncertainty due to its attribute. It is necessary to mitigate risks in sustainable development schedule. However, the existing literature includes limited quantitative research on robust strategies for sustainable development under uncertainty. Under these concerns, this paper presents a robust multi-objective optimization formulation for allocating labor across economic sectors to simultaneously satisfy economy (gross domestic product, GDP), environment (greenhouse gas (GHG) emissions), energy (electricity) and society (labor), in which the per capita GDP, per capita electricity consumption, per capita GHG and per capita rate of unemployment are robust uncertainty. There are three conflicting objectives in our sustainable development problem. The first objective maximizes the sectoral GDP, the second objective minimizes the sectoral electricity consumption and the third objective minimizes the sectoral GHG emissions. We adapt ɛ-constraint method to deal with the multiple objectives for the sustainable development problem, and obtain a more flexible result via an interactive decision-making process. Since uncertainty has often led to computationally intractable models, we reformulate the proposed model using robust optimization method into the tractable robust counterpart (RC) forms under two types of uncertainty sets. Finally, to demonstrate the effectiveness and applicability for our model, we conduct a case study for meeting year 2030 sustainable development of the United Arab Emirates (UAE). The numerical results show: (a) the robust multi-objective model is effective in uncertain environment and provides a reliable decision tool than deterministic model for integrated multi-objective sustainable development problem; (b) when the possibility distributions of uncertain parameters are available, fuzzy optimization model can provide better decision-making than robust model under box uncertainty set.

1Introduction

Sustainable development over economy, energy, environment and sociology, started in the 1990s for balancing energy-economy-environmental system and resource allocation [25, 32], has always been an important direction of efforts for global nations in the recent years. The United Nations agenda for sustainable development by the year 2030 proposes 17 sustainable development goals which include access to affordable, reliable and clean energy, sustained economic growth with full productive employment, and urgent action to mitigate environmental degradation [16]. The rapid development of economy has led to a remarkable rise in energy consumption and environmental pollution level. Economic development increases energy consumption that affects all economic and industrial activities, while the increased energy consumption is the primordial themes impacting any national economy. Furthermore, energy consumption, still mainly dependent on fossil fuels, such as oil, coal, petroleum and natural gas products, contributes to the increase of GHG emissions. At the same time, economic development affects population levels, while rapid population growth has also contributed the level of GHG emissions [34]. Therefore, in order to achieve sustainable development, the relationship between economic growth, energy consumption, environmental pollution and population growth is the key point that country should be concerned about. This is principally because of the enormous implications that this relationship has for policy prescriptions relative to sustainable development. Therefore, countries need to formulate appropriate policies to jointly address the problems of improving economic development, controlling energy consumption, reducing GHG emissions and meeting workforce allocation [18].

Inevitably, the key factor to the success of sustainable development is to commendably comprehend the interplay between economy, energy, environmental and social aspects simultaneously for finding emergency measures. Jayaraman et al. [18] studied sustainable development on energy-environment-economic system by building goal programming model. This paper extends the literature by investigating the multi-objective integer programming model for allocating labor resource to achieve the sustainable development. Our work pays attention to three objectives: maximum GDP growth, minimum electricity consumption and minimum GHG emissions, and models workforce allocation as constraint. In addition, the numbers of employees in each economic sector are chosen as decision variables. Subsequently, the proposed model is applied to key economic sectors of UAE to achieve sustainability related objectives by the year 2030. The results of the model can clearly show the causal relationship between the three objectives in UAE, which can provide a quantitative basis for policy makers and economists’ policy planning and regulation. Moreover, the results are also used to study the trade-off between resource allocation among different economic sectors and the impact of long-term investments on energy sustainability.

Policy makers need to face not only uncertain per capita contributions but also the rate of unemployment in strategic planning. Jayaraman et al. [16, 17] considered just sustainable development problem under certain environment. Due to the variety of per capita contribution, there often exist some impreciseness and ambiguousness in acquiring accurate data. Accordingly, policy makers have only access to the estimated values or the nominal value, but not enough to precisely pin down the exact data. At this point, the proposed optimization model in this paper incorporates the uncertainties encountered in optimizing sustainable development, such as per capita GDP, per capita electricity consumption, per capita GHG emissions and the rate of unemployment. In the presence of uncertainty, however, the problem becomes intractable. There exist several ways to treat data uncertainty in optimization, for example, robust optimization with distribution-free property, fuzzy technique and stochastic method. Jayaraman et al. [18, 19] employed stochastic and fuzzy methods to build goal programming models for sustainable development problem, respectively. In this paper, we utilize a relatively novel robust optimization method to address the uncertainty. By combining the aforementioned aspects, our multi-objective model is provided with robust uncertainty depicted with interval. We assume that uncertain parameters are taking value in two refined uncertainty sets parameterized in an affine way. As a result, we can transform the model with uncertainty to computational tractable formulation by deriving robust counterpart. When the possibility distributions of uncertain parameters are available, we leverage fuzzy optimization method for planning the sustainable development problem.

This paper contributes to the sustainable development problem in the following four aspects. (i) We propose a new robust multi-objective programming model for planning the sustainable development under uncertainty that involves three conflicting objectives. The first objective maximizes the sectoral GDP. The second objective minimizes the sectoral electricity consumption, and the third objective minimizes the sectoral GHG emissions. (ii) This paper assumes that the rate of unemployment exists within each economic sector in the sustainable development problem, which is obviously in line with the reality. In contrast to the related literature, the rate of unemployment is the first considered in planning the sustainable development problem. (iii) We characterize the per capita GDP, electricity consumption, GHG emissions and the rate of unemployment with robust uncertainty depicted with interval, which is different from the previous related literature. To obtain tractable formulation, we reformulate the proposed model into the equivalent robust counterpart forms under two types of uncertainty sets. (iv) We carry out case study based on the UAE to demonstrate the practicability of the model as well as the effectiveness of the approach. When the possibility distributions of uncertain parameters are available, we validate the fuzzy model and find that it can offer better decision than robust model under box uncertainty set. Further analysis provides the significant justification about studying the trade-off between resource allocation and long-term investments towards sustainability.

The structure of the remaining paper is organised as follows. Section 2 briefly reviews related literature. Section 3 gives the problem statements in detail and builds a robust multi-objective optimization model. Section 4 transforms the proposed model into tractable robust counterpart form. Section 5 conducts a case study based on UAE to verify the effectiveness of the model and approach. Section 6 draws some conclusions and outlines directions for further research.

2Literature review

The issue of sustainable development has always been a hot research topic of global nations in recent years. Chapple [11] provided a comprehensive framework on how to plan sustainable cities and regions, and used the case of the San Francisco Bay area as an example, and also provided an extremely interesting attempt to develop a broad conceptualisation of how to plan more sustainable cities and regions. Asif et al. [1] investiged the relationship between economic growth, urbanization, energy consumption, and CO₂ emission in GCC countries. Aznar-Marquez and Ruiz- Tamarit [2] studied economic growth in the presence of two environmental externalities: local and gives account of the marginal damage from the emissions flow and aggregate, or global, and relates to the extreme damage and found and completely characterized the socially optimal equilibrium, which showed that sustained growth is feasible along the balanced path. Jayaraman et al. [16] introduced a weighted goal programming model involving criteria on the economic development, the electricity consumption, the greenhouse gas emissions, and the total number of employees to determine optimal labor allocation across various economic sectors and applied the model to Gulf Cooperation Council Countries. Barreto [6] presented a theoretical framework that incorporates energy within an endogenous growth model considering the interaction and substitution between fossil fuels to reach the steady state where alternative energy production fuels the entire economy. Although the above literatures review the related study of sustainable development without considering uncertainty, some of the more important problems related to sustainability involve uncertainty. Thus, real-world optimization should consider uncertainty and find the method to deal with uncertainty.

Taking uncertainty into consideration, some researchers studied the topic on sustainability. Nomani et al. [27] utilized fuzzy goal programming to analyse environmental, energy and sustainability goals of India by the year 2030 with reference to the key economic sectors of India. Jayaraman et al. [19] studied the sustainable development of UAE by the year 2030 considering economy, population, energy and environment by fuzzy goal programming. Gupta et al. [15] analysed the economic development and sustainable goals on GDP, the electricity consumption and the GHG emissions of India and provides an optimal allocation of resources by building fuzzy goal programming model. Jayaraman et al. [18] utilized stochastic goal programming to consider the sustainable development planning for economic development, energy consumption, GHG emissions and population growth of UAE by the year 2030. Rodriguez-Gonzalez et al. [33] presented a stochastic integrated ecological-economical-social model which includes a crucial social aspect in the sustainability study such as inequity, and analyzed four typical kinds of societies with different levels of development and inequity. For a thorough coverage of development and recent advances in global research on sustainability and sustainable development, the interested reader can refer to Olawumi and Chan [28]. Some of the above literatures on sustainable development mainly considered the uncertainty of goals and used fuzzy or stochastic goal programming method to study the sustainable development. However, the uncertainty of goals is mainly due to the presence of uncertainty about input data. Unlike these documents, this paper assumes that input data including per capita GDP, per capita electricity consumption, per capita GHG emissions and the rate of unemployment are modeled as uncertain parameters.

When the possibility distributions of uncertain input data are available, we study sustainable development problem by fuzzy optimization method. Fuzzy set, pioneered by Zadeh [39], is a strong tool to deal with subjective uncertainty. Since then, fuzzy theory has been explored in the literature. Liu and Liu [21] presented a novel concept of credibility measure and expected value of fuzzy variable. Qiu and Zhang [31] discussed the properties of symmetric fuzzy numbers, showed an equivalent characterization of convex fuzzy sets, and presented a way to construct a symmetric convex fuzzy set with a convex fuzzy set. Qiu et al. [30] obtained some algebraic proper ties and topological properties of the quotient space of fuzzy numbers with respect to the equivalence relation and introduced a new concept of convergence under which the quotient space was complete. At the same time, fuzzy theory integrated with optimization method has been applied successfully to more and more application-oriented fields [3, 22, 23]. The inclusion of fuzzy uncertainty in the context of sustainable development is a challenging issue in terms of modeling and solution. Although several studies may have considered the related problem, the uncertain parameters considered in those documents are characterized by membership functions, while the uncertain data in our paper are measured by possibility distribution functions and a credibilistic sustainable development model is built.

In reality, the possibility distributions of uncertain model parameters are often not easily known. Under the consideration, this paper studies the sustainable development problem in the presence of uncertainty by a relatively novel method - Robust optimization. Robust optimization, first proposed by Soyster [35], is a new approach to deal with the programming problems under uncertain input data without the distribution information, which is different from the fuzzy or stochastic method in nature. Ben-Tal and Newmirovsky [7, 8] successively proved that the robust convex formulation of programming problems with uncertainty is tractable when the uncertainty set is described as a box or an ellipsoid. An extension of this method can be referred to the literature [9, 10, 20]. For a thorough coverage of developments and recent advances in robust optimization, the interested reader can refer to [13]. Due to the superiority of robust optimization, robust optimization approach has been applied in many areas, including data envelopment analysis [5], energy system planning [12], sustainable energy supply systems [24], portfolio of R&D projects [36], single-period inventory problem [14], p-hub median problem [37] and so on. This implies that robust optimization method in uncertain programming problems can address uncertainty and the computational intractable issue in the modeling process, which motivates us to consider the sustainable development problem over economic development, energy consumption, environment protection and population growth by robust optimization method. In this paper, we take advantage of the robust optimization method to resist the uncertain per capita contributions on GDP, electricity consumption, GHG emissions and the rate of unemployment when these uncertain data take their values in box+ellipsoidal and box+generalized budgeted uncertainty sets.

3Robust multi-objective sustainable development model

3.1Problem description

We consider a problem of planning sustainable development. Our problem stems from the increased concern about the relation between energy consumption, global population growth and GHG emissions and the effects on long-run sustainability. The rapid growth of economy requires the massive use of labor resources and energy consumption, in turn leading to high-level pollution. The long-run dynamics of the sustainability criteria offers unique challenges and opportunities to better plan resource allocation. Therefore, we hope to allocate resources (labor) to meet economic growth, reduce GHG emissions and energy consumption in order to achieve sustainable development. Here we take economic sectors as units to allocate labors to simultaneously satisfy GDP growth, electricity consumption, GHG emissions and the number of employees in each economic sector goals necessary to meet an economic sustainability. Therefore, we consider GDP, electricity consumption and GHG emissions as objective functions optimized simultaneously, and the workforce allocation as a constraint satisfying a specified level. Here we consider several economic sectors and following assumptions:

• Objective functions are assumed to be linear: Each objective is supposed to be linearly dependent on each of decision variables; The same as constraints.
• Rate of unemployment in each sector is considered, which is uncertain prior to decision making.
• Set a specific time period for one year.

3.2Notations on the sustainable development problem

To formally characterize the sustainable development problem, we first list all symbols and parameters used in the problem formulation, as shown in Table 1.

Table 1

Notation

Notation	Description
Sets
[I] = {i = a, b, c, d}	a set of uncertain parameters
[J] = {j = 1, . . . , n}	a set of economic sectors
Ui	uncertainty set i ∈ [I]
Parameters
a_j	the per capita GDP contribution of each employee in the economic sector j
b_j	the per capita electricity consumption of each employee in the economic sector j
c_j	the per capita GHG emissions contribution of each employee in the economic sector j
d_j	the per capita workforce contribution of each employee in the economic sector j
e_j	the current number of employees each economic sector j
t	the predetermined upper bound on the total number of employees
z_i	an L-dimensional perturbation vector i ∈ [I]
Decision variables
x_j	the number of employees in the economic sector j

3.3Mathematical model

First we give a detailed account of three objectives in our problem.

• The formulation of the total sectorial GDP under resulting optimal resources allocation across whole country in a year is as follows:

(1)

F1(x,a(za))=a(za)Tx,

where x = (x₁ ; x₂ ; . . . ; x_n), a (z_a) = (a₁ (z_a); a₂ (z_a); … ; a_n (z_a)).

• The cumulative quantity F₂ (x, b (z_b)) of per capita electricity consumption (in Giga watt hours (Gwh)) across all economic sectors is as follows:

(2)

F2(x,b(zb))=b(zb)Tx,

where b (z_b) = (b₁ (z_b) ; b₂ (z_b) ; … ; b_n (z_b)).

• F₃ (x, c (z_c)), the total level of sectorial GHG emissions (in Giga gram equivalent of CO₂), is expressed as follows:

(3)

F3(x,c(zc))=c(zc)Tx,

where c (z_c) = (c₁ (z_c) ; c₂ (z_c) ; … ; c_n (z_c)).

The per capita GDP contribution a, per capita electricity consumption b and per capita GHG emissions c are provided with a certain degree of uncertainty. We may not be able to evaluate them accurately, which lead to the uncertainty of the total sectorial GDP F₁ (x, a (z_a)), the total sectorial electricity consumption F₂ (x, b (z_b)) and the total sectorial GHG emissions F₃ (x, c (z_c)).

By introducing auxiliary decision variables F10 , F20 and F30 , the uncertain objectives (1)-(3) can be converted to certain objectives, which are defined in expressions (4)-(9).

Next we introduce the constraints of our optimization problem. Constraint (10) guarantees that for all economic sectors the total employed work force does not exceed t inducing thus an upper bound on the total number of employees, in which d (z_d) = (d₁ (z_d) ; d₂ (z_d) ; … ; d_n (z_d)) is also uncertain. We impose the constraint (11) to ensure that the worker number x_j of each economic sector j ∈ [J] is integer and take into account that for each economic sector j ∈ [J] the employed worker number satisfies the current number of employees per sector, which guarantees the sustainability.

The robust multi-objective sustainable development problem can be modeled as follows:

(4)

maxF10

(5)

minF20

(6)

minF30

(7)

s.t.a(za)Tx≥F10,∀a∈Ua

(8)

b(zb)Tx≤F20,∀b∈Ub

(9)

c(zc)Tx≤F30,∀c∈Uc

(10)

d(zd)Tx≤t,∀d∈Ud

(11)

xj≥ej,xj∈Z+,j∈[J],

which is a multi-objective semi-infinite programming problem.

3.4Model transformation into single objective

In this subsection, we consider transforming this multi-objective programming into a single objective programming. Here we use ɛ-constraint method (see [38]). According to the preferences of policy makers, GDP is considered as the objective of the single-objective programming model. As a result, the multi-objective programming problem can be transformed to the following single-objective programming:

(12)

maxF10s.t.F20≤(1-ɛE)EmaxF30≤(1-ɛG)Gmaxconstraints (7)-(11).

where ɛ_E represents target decrease ratio of total electricity consumption, and ɛ_G represents target decrease ratio of the total GHG emissions. E_max and G_max denote the maximum total sectorial electricity consumption and maximum total sectorial GHG emissions, respectively, when the total sectorial GDP is maximum and the energy and environmental effects are not considered.

The model (12) can well resist the uncertainty. However, the model faces the computationally intractable challenge in virtue of the semi-infinite constraints of the model. In order to circumvent this difficulty of the model, our goal now is to built a representation capable of expressing equivalently the semi-infinite linear constraints of the model as a finite system of explicit convex constraints. The outlined strategy allows us to focus on robust optimization method. Robust optimization method can provide the RC of model (12) under uncertainty set, which is a finite computationally tractable system of convex constraints.

4Reformulation

In the section, we will reformulate the computationally tractable forms of our robust model (12) under two types of uncertainty sets. The uncertain per capita contributions a, b, c and d in problem (12) are varying in the given uncertainty set parameterized in an affine way by perturbation vector zi=(zi1,zi2,...,ziL)T (i ∈ [I]) varying in given perturbation sets Zi, i∈[I] . For instance, constraint (7) can be written as follows:

a(za)Tx≥F10,Ua={a(za)=[a0]T+∑l∈[L]alzal:∀ za ∈ Za}.

Similarly, constraints (8)-(10) can also be equivalently written as this kind of form. The uncertain parameters [a ; b ; c ; d] represent basic perturbations from the nominal data [a⁰ ; b⁰ ; c⁰ ; d⁰].

4.1RC model under box+ellipsoidal uncertainty

In this subsection we study the RC model of problem (12) with uncertain per capital contributions a, b, c, d varying in box+ellipsoidal uncertainty set parameterized in an affine way by the following box+ellipsoidal perturbation set:

(13)

Zi={zi∈ℝL, -1≤zil≤1,∑l∈[L](zilσl)2≤Ωi, l∈[L]},

where σ_l > 0 is the given parameter, and Ω_i > 0 (i ∈ [I]) is the adjustable safe parameter controlling the size of the uncertainty set.

Theorem 1. For sustainable development problem (12), suppose that the uncertain per capital contributions a, b, c, d are parameterized by perturbation vector z_i belonging to the box+ellipsoidal perturbation set (13). Then the RC model of sustainable development problem (12) can be represented by the following system:

(14)

max F10s.t. F10≤(1-ɛE)EmaxF30≤(1-ɛG)Gmax-∑l∈[L]|rl|-Ωa∑l∈Lσl2sl2+[a0]Tx≥F10rl+sl=[al]Tx,l∈[L]

∑l∈[L]|ql|+Ωb∑l∈[L]σl2ωl2+[b0]Tx≤F20ql+ωl=-[bl]Tx,l∈[L]∑l∈[L]|gl|+Ωc∑l∈[L]σl2hl2+[c0]Tx≤F30gl+hl=-[cl]Tx,l∈[L]∑l∈[L]|fl|+Ωd∑l∈[L]σl2kl2+[d0]Tx≤tfl+kl=-[dl]Tx,l∈[L]constraint (11).

Proof. First we build the RC form of constraint (7) on GDP. The perturbation vector z_a affecting constraint (7) on GDP belongs to following perturbation set

Za={za∈ℝL, -1≤zal≤1,∑l∈[L](zalσl)2≤Ωa, l∈[L]}.

The conic representation of the perturbation set becomes:

Za={za∈ℝL:A1za+a1∈K1,A2za+a2∈K2},

where

A₁z_a ≡ [z_a ; 0], a₁ = [0_L×1 ; 1] and K1={(y,t)∈ℝL×ℝ:t≥∥y∥∞} , whence K*1={(y,t)∈ℝL×ℝ:t≥∥y∥1} , which is referred to as the cone dual to K¹.

A₂z_a = [∑^-1z_a ; 0] with ∑ = Diag {σ₁, . . . , σ_L}, a₂ = [0_L×1 ; Ω_a] and K2={(y,t)∈ℝL×ℝ:t≥∥y∥2} , whence K*2=K2 .

Setting φ¹ = [τ₁ ; π₁], φ² = [τ₂ ; π₂] with one-dimensional π₁, π₂ and L-dimensional τ₁, τ₂, the inequality (7) becomes the following system of constraints in variables π, τ, x:

π1+Ωaπ2-[a0]Tx≤-F10,(τ1+Σ-1τ2)l=[al]Tx,l∈[L],∥τ1∥1≤π1[⇔[τ1;π1]∈K*1],∥τ2∥∞≤π2[⇔[τ2;π2]∈K*2].

For every feasible solution to the system, one have π1≥π¯1≡∥τ1∥1 , π2≥π¯2≡∥τ2∥2 . We can eliminate the variables π₁, π₂ from this system, and the solution obtained still is feasible when replacing π₁, π₂ with π¯1 , π¯2 . Based on the above analysis, the reduced system in variables x, r = τ₁, s = Σ^-1τ₂ reads:

(15)

-∑l∈[L]|rl|-Ωa∑l∈Lσl2sl2+[a0]Tx≥F10rl+sl=[al]Tx,l∈[L]

in variables x, r, s, and every feasible solution to this system is feasible for Eq. (7).

The above proves that the system (15) is the RC of Eq. (7). Similar to the above derivation process, constraints (8)-(10) with box+ellipsoidal uncertainty set can also be turned into the kind of RC forms.

As a result, the uncertain sustainable development problem (12) can be equivalently transformed into the RC model (14).□

Remark 2. The RC model (14) with finite constraints is computationally tractable. Moreover, any x satisfying the RC model (14) is a solution of problem (12).

4.2RC model under box+generalized budgeted uncertainty

In this subsection we study the RC model of problem (12) with uncertain per capital contributions a, b, c, d varying in box+generalized budgeted uncertainty set parameterized in an affine way by the following box+generalized budgeted perturbation set:

(16)

Zi={zi∈ℝL:-1≤zil≤1,∑l∈[L]|zilσl|≤γi,l∈[L]},

where σ_l is the given parameter, and γ_i > 0 (i ∈ [I]) is the adjustable safe parameter controlling the size of the uncertainty set.

Theorem 2. For the sustainable development problem (12), suppose that the uncertain per capital contributions a, b, c, d are parameterized by perturbation vector z_i belonging to the box+generalized budgeted perturbation set (16). Then the RC model of sustainable development problem (12) can be represented by the following system:

(17)

Proof. This proof is similar to Theorem 1.□

Remark 2. The RC model (17) is also computationally tractable. Every feasible solution to the RC model (17) is feasible for uncertain problem (12).

5Case study

In the section, we present a case study about year 2030 sustainability goals of UAE from two perspectives: (a) uncertain per capita contributions are distribution-free; (b) the possibility distribution of uncertain per capita contributions are available. All mathematical models are solved by CPLEX studio1263.win-x86-64 on personal computer (Intel(R) Core(TM) i5-4200M 2.50GHz CPU and RAM 4.00GB).

5.1Data source and analysis

The UAE vision 2021 highlights an ambition to attain a ‘competitive and resilient economy’ and a ‘sustainable environment’ [19]. Mokri et al. [26] pointed out that UAE is the world’s 10th largest per capita electricity consumer with an average annual increase rate of about 8.8% during the last decade. During the years 2006–2011, the annual increase in the electricity demand in UAE was as high as 10.8%, which was closely followed the trend in annual population growth of 11% during the same period [26]. In addition, 97.5% of electricity generation in the UAE relies on natural gas-powered plants [29], which naturally leads to the environmental pollution including the production of GHG emissions and other particulate matter. Considering the currently high, yet increasing energy consumption, GHG emissions and population in the UAE, an expected harmonious sustainable development towards economy, environment, energy and population, it is imperative to explore sustainability of the UAE. Accordingly, we find the data on the per capita contributions in several key sectors of the UAE presented in Jayaraman et al. [18]. There are 8 economic sectors identified by Vellinga are considered: (1) agriculture, (2) crude oil, natural gas and quarrying, (3) manufacturing and electricity, (4) construction and real estate, (5) trade and transport, (6) trade and transport, (7) banking and financial corporations and (8) government, social and personal services. According to Jayaraman et al. [18], the population size t are estimated as 9452000 for 2030 year.

5.2Robust method for sustainable development experiments

In the problem the uncertain per capita contributions a, b, c and d with distribution-free property, basic perturbations around nominal value a⁰, b⁰, c⁰ and d⁰, are varying in box+ellipsoidal and box+generalized budgeted uncertainty sets. We assume the nominal data [a⁰ ; b⁰ ; c⁰] of the uncertain parameters [a ; b ; c] are provided in Jayaraman et al. [18], and the uncertain parameter d_j takes its value in the segment [0.98, 1.0] with dj0=0.99 (j = 1, . . . , 8) as nominal value. The problem instances are generated by setting Ω to 0.5, 1.0, 1.5, 2.0 and γ to 1.0, 2.0, 3.0 and 4.0. For each assigned safe parameter values Ω and γ, for the sake of presentation, we assume the parameter σ_l are equal for each l and set it to 0.1, 0.3, 0.5 and 0.7, at the same time, the data variability of uncertain parameters [a ; b ; c] are set to 10% relative to their nominal values [a⁰ ; b⁰ ; c⁰] and (ɛ_E, ɛ_G) = (0.25, 0.60). Table 2 summaries the perturbation of data.

Table 2

Perturbation values of uncertain parameters

Sectors	Perturbation values
	ajl	bjl	cjl	djl
Agriculture	0.003521739	0.000478696	0.001728696	0.0100
Crude oil, natural gas and quarrying	0.46969697	0.005912121	0.171707576	0.0100
Manufacturing and electricity	0.018134206	0.002502291	0.006629133	0.0100
Construction and real estate	0.00838565	0.001873543	0.000267227	0.0100
Trade and transport	0.017690457	0.001614274	0.000627506	0.0100
Restaurant and hotels	0.008095238	0.000738571	0.000258095	0.0100
Banking and financial corporations	0.105138889	0.014509722	0.003349306	0.0100
Government, social and personal services	0.009569444	0.000872083	0.000305000	0.0100

Table 3

Computational results with box+ellipsoidal uncertainty set

Ω	σ_l	x₁	x₂	x₃	x₄	x₅	x₆	x₇	x₈
0.5	0.1	230000	1999329	611000	1338000	4072813	210000	363884	720000
	0.3	230000	1977456	611000	1338000	4095152	210000	358507	720000
	0.5	230000	1956015	611000	1338000	4117006	210000	353165	720000
	0.7	230000	1934993	611000	1338000	4138389	210000	347855	720000
1.0	0.1	230000	1988338	611000	1338000	4084044	210000	361191	720000
	0.3	230000	1945453	611000	1338000	4127755	210000	350506	720000
	0.5	230000	1904220	611000	1338000	4169611	210000	339949	720000
	0.7	230000	1864549	611000	1338000	4209619	210000	329514	720091
1.5	0.1	230000	1977456	611000	1338000	4095152	210000	358507	720000
	0.3	230000	1914379	611000	1338000	4159315	210000	342577	720000
	0.5	230000	1854864	611000	1338000	4219481	210000	326915	720000
	0.7	230000	1806035	611000	1338000	4271191	210000	308867	720013
2.0	0.1	230000	1966682	611000	1338000	4106139	210000	355832	720000
	0.3	230000	1884195	611000	1338000	4189876	210000	334715	720000
	0.5	230000	1807780	611000	1338000	4266809	210000	314040	720000
	0.7	230000	1799913	611000	1338000	4297493	210000	273229	720000

In Table 2, al=ajl(0j-1,1;1;08-j,1) , bl=bjl (0_j-1,1 ; 1 ;0_8-j,1), cl=cjl(0j-1,1;1;08-j,1) and dl=djl(0j-1,1 ; 1; 0_8-j,1) 1 ≤ j, l ≤ 8, and we assume l = j.

5.2.1Computational results

In order to assess the applicability of the proposed model and the approach, 32 test problems in terms of the different sizes of box+ellipsoidal and box+generalized budgeted uncertainty sets and parameter σ_l are generated when data variability of uncertain parameters are set to 10% relative to their nominal values and (ɛ_E, ɛ_G) = (0.25, 0.60). The summary of the results on the test problems for box+ellipsoidal and box+generalized budgeted uncertainty sets are listed in Tables 3 and 4, Figs. 1–3.

Fig.1

Impacts of parameter σ_l and uncertainty size on the GDP F10 .

Fig.2

Impacts of parameter σ_l and uncertainty size on the electricity consumption F20 .

Fig.3

Impacts of parameter σ_l and uncertainty size on the GHG emissions F30 .

Table 4

Computational results with box+generalized budgeted uncertainty set

γ	σ_l	x₁	x₂	x₃	x₄	x₅	x₆	x₇	x₈
1.0	0.1	230000	1988464	611000	1338000	4082356	210000	363531	720000
	0.3	230000	1945823	611000	1338000	4122557	210000	357602	720000
	0.5	230000	1904824	611000	1338000	4160733	210000	351903	720000
	0.7	230000	1865376	611000	1338000	4197000	210000	346422	720000
2.0	0.1	230000	1966932	611000	1338000	4102717	210000	360537	720000
	0.3	230000	1884911	611000	1338000	4179099	210000	349136	720000
	0.5	230000	1808926	611000	1338000	4248058	210000	338580	720000
	0.7	230000	1806446	611000	1338000	4264391	210000	317262	720002
3.0	0.1	230000	1945823	611000	1338000	4122557	210000	357602	720000
	0.3	230000	1827393	611000	1338000	4231468	210000	341145	720000
	0.5	230000	1805826	611000	1338000	4268506	210000	311906	720000
	0.7	230000	1802944	611000	1338000	3180032	210000	353604	1769774
4.0	0.1	230000	1925124	611000	1338000	4141890	210000	354725	720000
	0.3	230000	1807686	611000	1338000	4256106	210000	327947	720092
	0.5	230000	1803425	611000	1338000	3200765	210000	356099	1747638
	0.7	230000	1798327	611000	1338000	4306699	210000	260879	720090

As can be seen in Tables 3 and 4, adjusting the value of σ_l impacts the numbers of labor allocation x_j for each economic sector j. On the one hand, it can be observed the distribution of labor force across 8 economic sectors under box+ellipsoidal uncertainty set. For example, under Ω = 0.5, by altering the parameter value of σ_l for 0.3 and 0.5, x₂ changes from 1977456 to 1956015, x₅ changes from 4095152 to 4117006, and x₇ changes from 358507 to 353165, while x₁, x₃, x₄, x₆ and x₈ keep unchanged. This event similarly happens to other cases (e.g., σ_l=0.1, 0.7). On the other hand, we can see that the allocation of labor force in each economic sector under box+generalized budgeted uncertainty set. For example, under γ = 4, by manipulating the parameter value σ_l for 0.1 and 0.3, x₂ and x₇ are down, x₅ and x₈ go up, while x₁, x₃, x₄ and x₆ remain immovability. As a consequence, for each Ω and γ, the x₂ and x₇ decrease and x₅ increases when σ_l is decreasing. The fact demonstrates that policy makers can easily adjust the parameter value σ_l in order to obtain better allocation of labor force across 8 economic sectors to achieve sustainable development.

According to Tables 3 and 4, it can be also observed that by adjusting the safe parameters Ω and γ, the number of labor allocation for each economic sector are different. For each σ_l, if the safe parameter Ω takes its value as 0.5, 1.0, 1.5 and 2.0, respectively, the x₂ and x₇ decrease evidently and x₅ is gradually increasing, while the numbers of other sectors are largely remain unchange. For example, under σ_l = 0.1, by altering the safe parameter Ω for 0.5 and 0.1, x₂ decreases from 1999329 to 1988338, x₇ decreases from 363884 to 361191, and x₅ increases from 4082356 to 4197000. In the same case, when the safe parameter γ takes its value as 1.0, 2.0, 3.0 and 4.0, x₂ and x₇ decrease evidently, and x₅ gradually increases, while the numbers of employees in other sectors are largely remain immovability. These imply that policy makers can adjust Ω or γ to better allocate the labor force to achieve sustainable development. By comparing Tables 3 and 4, it is found that the numbers of labor allocation are different under different uncertainty set. For example, under σ_l= 0.5, x₂ are 1807780 and 1808926, respectively, when Ω= 2 and γ= 2. This event similarly happens to σ_l= 0.1, 0.3 and 0.7. Based on the analysis, policy makers can choose appropriate uncertainty set on the basis of the situation.

Next the alteration of parameter σ_l and the different sizes of box+ellipsoidal and box+generalized budgeted uncertainty sets with respect to the three objective values, respectively, are plotted in Figs. 1–3, where the left horizontal axis corresponds to parameter σ_l, the right horizontal axis corresponds to different levels of uncertainty set and the vertical axis corresponds to the objective value. As expect, the values of GDP, electricity consumption and GHG emissions decrease monotonically by enlarging the size of the uncertainty set when parameter σ_l is enforced. Moreover, the three objective values are also less at higher parameter σ_l under any choice of the level of uncertainty sets. Furthermore, comparing the three objective values under box+ellipsoidal and box+generalized budgeted uncertainty sets when Ω = γ = 2, σ_l = 0.5 reveals that the three objective values are less under box+ellipsoidal uncertainty than the values under box+generalized budgeted uncertainty, respectively. These findings appear the impacts of parameter σ_l, the different uncertainty sets and its size on GDP, electricity consumption and GHG emissions, which helps policy makers make better decisions.

5.2.2Model comparison

To evaluate the advantages of the proposed model, some comparison results are obtained in two different environments: deterministic model and robust model. In the case that a, b, c and d take its nominal values a⁰, b⁰, c⁰ and d⁰, respectively, the uncertain parameters are reduced to the deterministic ones. Here (ɛ_E, ɛ_G)= (0.25, 0.60). Table 5 and Fig. 4 summarize the comparison of solution results.

Table 5

Model comparison results

	Nominal value	Ω = 2.5 (σ_l = 0.7)	Ω = 3 (σ_l = 0.5)	γ = 5 (σ_l = 0.7)	γ = 6 (σ_l = 0.5)
x₁	230000	230000	230000	230000	230000
x₂	2010432	1795730	1798637	1795824	1797216
x₃	611003	611000	611000	611000	611000
x₄	1338000	1338000	1338000	1338000	1338000
x₅	4061453	4316998	4303804	4317165	4310429
x₆	210000	210000	210000	210000	210000
x₇	366586	246455	264455	246443	255478
x₈	720000	720001	720249	720143	720143
F20	286900	283300	283600	283330	283480
F30	3540600	3484900	3489500	3485200	3487700

Fig.4

Comparison results under different uncertainty set, σ_l = 0.5.

By Table 5, we can clearly see that the solutions for the deterministic model are a little different from the proposed model in this paper. For x₂, under Ω = 2.5, σ_l = 0.7, the nominal solution 2010432 is higher than robust solution 1795730. The model which includes robustness has substantial optimal decision than the one without uncertainty. This is principally because the uncertainty of parameters profoundly reflects the optimal solution. In other words, the uncertainty of data can not brook being ignored.

In addition, Fig. 4 reports the comparison results of objective GDP under box+ellipsoidal and box+ generalized budgeted uncertainty sets, respectively. It is easily found that robust GDP is smaller than nominal GDP. Not surprisingly, although the robust model incurs lower GDP, the electricity consumption and GHG emissions shown in Table 5 are a little higher than our robust values. That is, our robust model achieves lower electricity consumption and GHG emissions. Therefore, policy makers should be able to formulate a good strategic planning for the sustainable development by 2030 year for UAE.

5.2.3Sensitivity analysis

In order to recognize the most significant parameters of robust model, several sensitivity analysis are carried out and the impact of parameter alteration on objective GDP is investigated. We perform sensitivity analysis on parameters (ɛ_E, ɛ_G) and the different data variability of uncertain parameters a, b and c, while d still takes its value in segment [0.98, 1.0] with nominal value 0.99. Here Ω = γ = 2, σ_l=0.5. The analysis results are provided in Fig. 5 and Table 6.

Fig.5

Impacts of data variability on the GDP F10 .

Table 6

Sensitivity of GDP, σ_l = 0.5

ɛ_E	F10	F10	ɛ_G	F10	F10
(ɛ_G = 0.6)	Ω = 2	γ = 2	(ɛ_E = 0.25)	Ω = 2	γ = 2
0.05	9493200	9538300	0.45	11657000	11683000
0.10	9381900	9429800	0.50	10789000	10812000
0.15	9269500	9312900	0.55	9915100	9941800
0.20	9155500	9192100	0.60	9039800	9071300
0.25	9039800	9071300	0.65	8162800	8200800
0.30	8922000	8950500	0.70	7283700	7330400
0.35	8785400	8815700	0.75	6402200	6447100

Figure 5 depicts the sensitivity of objective GDP to the alteration of the data variability of the uncertain parameters under box+ellipsoidal and box+generalized budgeted uncertainty sets, respectively. The horizontal axis corresponds to the data variability (%) and the vertical axis corresponds to the GDP. It can be inferred that by an increase of data variability, the value of GDP gets significantly decrease for any choice of box+ellipsoidal and box+generalized budgeted uncertainty sets.

The effect of parameters ɛ_E and ɛ_G on the objective GDP is shown in Table 6. It is intuitive that as the ɛ_E and ɛ_G become higher, respectively, the GDP is lower. By Table 6, if the ɛ_E takes its value 0.1, which corresponds to the ellipsoidal uncertainty set, the optimal objective GDP = 938190 is higher than the value for ɛ_E = 0.15, 0.20, 0.25, 0.30, 0.35. Since the electricity consumption and GHG emissions decrease as the ɛ_E and ɛ_G increase, according to the relation on the three objectives, we find that the results are stand to reason.

5.3Fuzzy method for sustainable development experiments

In this subsection, we study a sustainable development problem by fuzzy optimization method [21, 23] instead of robust optimization, once the possibility distributions of uncertain per capita contributions a, b, c and d are available. In this case, we build the following fuzzy sustainable development model in which the uncertain per capita contributions a, b, c and d are assumed as fuzzy variables.

(18)

maxF10s.t.F20≤(1-ɛE)EmaxF30≤(1-ɛG)GmaxCr{aTx≥F10}≥1-α

Cr{bTx≤F20}≥1-βCr{cTx≤F30}≥1-γCr{dTx≤t}≥1-ηconstraint (11),

where α, β, γ and η ∈ (0, 0.5) are predetermined violation levels. In order to obtain a meaningful solution to the fuzzy sustainable development problem (18), we assume that the per capita contributions a_j, b_j, c_j and d_j are characterized by triangular fuzzy variables, i.e., aj=Tri[ajl0;aj0;ajr0] , bj=Tri[bjl0;bj0;bjr0] , cj=Tri[cjl0;cj0;cjr0] and dj=Tri[djl0;dj0;djr0] , which drives naturally a^Tx, b^Tx, c^Tx and d^Tx to be fuzzy variables. Based on the above analysis, model (18) can be transformed as follows:

(19)

max F10s.t. F20≤(1-ɛE)EmaxF30≤(1-ɛG)Gmax

F10-∑j=18ajl0xj≤2α[∑j=18(aj0-ajl0)xj]∑j=18bjr0xj-F20≤2β[∑j=18(bjr0-bj0)xj]∑j=18cjr0xj-F30≤2γ[∑j=18(cjr0-cj0)xj]∑j=18djr0xj-t≤2η[∑j=18(djr0-dj0)xj]constraint (11),

where ajl0=(1-Δjla)aj0 , bjr0=(1+Δjrb)bj0 , cjr0=(1+Δjrc)cj0 and djr0=dj0+Δjrd .

5.3.1Computational results

Let us set (ɛ_E, ɛ_G) = (0.25, 0.60) and consider the values of aj0, bj0 , cj0 and t the same as mentioned before. Six test problems in terms of six series of the tolerance parameters [Δla, Δrb, Δrc, Δrd] are generated when violation levels α = β = γ = η = 0.05: (I) [Δla, Δrb, Δrc, Δrd]=[0.05, 0.05, 0.05, 0.02] ; (II) [Δla, Δrb, Δrc, Δrd]=[0.05, 0.1, 0.150.02] ; (III) [Δla, Δrb, Δrc, Δrd]=[0.1, 0.1, 0.1, 0.01] ; (IV) [Δla, Δrb, Δrc, Δrd]=[0.15, 0.1,0.05, 0.01] ; (V) [Δla, Δrb, Δrc, Δrd]=[0.1, 0.1, 0.15, 0.01] ; (VI) [Δla, Δrb, Δrc, Δrd]=[0.15, 0.15, 0.15, 0.02] . The summary of the results on the test problems are listed in Table 7. In addition, 6 test problems in terms of six series of the violation levels [α, β, γ, η] are generated: (i) α = β = 0.5, γ = η = 0.1; (ii) α = β = γ = η = 0.1; (iii) α = β = 0.1, γ = η = 0.15; (iv) α = β = γ = η = 0.15; (v) α = β = 0.15, γ = η = 0.20; (vi) α = β = γ = η = 0.20. The computational results are provided in Table 8.

Table 7

Computational results with different tolerance parameters with violation levels 0.05

	(I)	(II)	(III)	(IV)	(V)	(VI)
x₁	230000	230000	230000	230000	230000	230000
x₂	1900587	1739197	1810307	1891133	1735890	1740476
x₃	611002	611003	611004	611004	611004	611004
x₄	1338000	1338000	1338000	1338000	1338000	1338000
x₅	4163307	4364020	4308933	4255045	4358548	4443135
x₆	210000	210000	210000	210000	210000	210000
x₇	298045	258721	233217	206279	258019	178326
x₈	720000	720000	720000	720000	720000	720000
F10	9830600	9101100	8943000	8796400	8656600	8187600
F20	284350	283510	283090	283090	283090	283510
F30	3501100	3488100	3481500	3481500	3481500	3488100

Table 8

Computational results with different violation levels with tolerance parameters (III)

	(i)	(ii)	(iii)	(iv)	(v)	(vi)
x₁	230000	230000	230000	230000	230000	230000
x₂	1831165	1830866	1852123	1851818	1873483	1873172
x₃	611004	611003	611004	611004	611003	611003
x₄	1338000	1338000	1338000	1338000	1338000	1338000
x₅	4302703	4284227	4277719	4258869	4252076	4232842
x₆	210000	210000	210000	210000	210000	210000
x₇	228069	246845	241595	260750	255397	274942
x₈	720000	720000	720000	720000	720000	720000
F10	9026200	9139200	9225000	9339500	9427800	9543900
F20	283510	283510	283930	283930	284350	284350
F30	3488100	3488100	3494600	3494600	3501100	3501100

From the computational results in Tables 7 and 8, we find that with the different tolerance parameters and violation levels, the solutions and objective values are different. On the one hand, the labor allocation decision is relevant to the tolerance parameters and violation levels. Specifically, the second sector, the fifth sector and the seventh sector have noticeable changes in the number of workers allocated as the tolerance parameters and violation levels change. For example, the employee number for the second sector changes from 1900587 to 1739197 by altering the tolerance parameters from (I) to (II) for the test problems in Table 7, and the employee number for the second sector changes from 1831165 to 1830866 by altering violation levels from (i) to (ii) for the test problems in Table 8. On the other hand, the objective values are also relevant to the tolerance parameters and violation levels. For example, the GDP decreases from 9101100 to 8943000, the electricity consumption decreases from 283510 to 283090, the GHG emissions decreases from 3488100 to 3481500 by changing the tolerance parameters from (II) to (III) in Table 7, and when violation levels change from (ii) to (iii) in Table 8, the GDP increases from 9139200 to 9225000, the electricity consumption increases from 283510 to 283930, the GHG emissions increases from 3488100 to 3494600.

5.3.2Comparison study

Fuzzy model with fixed distributions is less conservative than the model with only supporting information. To see this, we conduct the comparative experiments by compare fuzzy model (19) with triangular probability distribution and the robust model under box uncertainty. For the robust model, we consider the supporting information of a_j, b_j, c_j and d_j, respectively, are [(1-Δjla)aj0, (1+Δjla)aj0] , [(1-Δjrb)bj0, (1+Δjrb)bj0] , [(1-Δjrc)cj0, (1+Δjrc)cj0] and [dj0-Δjrd, dj0+Δjrd] . Table 9 summarizes the comparison results under case (I), (III) and (VI).

Table 9

Comparison results with violation levels 0.05

	(I)		(III)		(VI)
	fuzzy	robust	fuzzy	robust	fuzzy	robust
x₁	230000	230000	230000	230000	230000	230000
x₂	1900587	1877012	1810307	1790128	1740476	1710800
x₃	611002	611003	611004	611004	611004	611003
x₄	1338000	1338000	1338000	1338000	1338000	1338000
x₅	4163307	4180190	4308933	4333008	4443135	4472538
x₆	210000	210000	210000	210000	210000	210000
x₇	298045	285795	233217	219860	178326	159659
x₈	720000	720000	720000	720000	720000	720000
F10	9830600	9664600	8943000	8750600	8187600	7914900
F20	284350	282670	283090	282670	283510	282670
F30	3501100	3475000	3481500	3475000	3501100	3475000

From the computational results, we can clearly observe that the obtained optimal decision of our fuzzy model and the one of robust model under box uncertainty are completely different. For the objective value F10 , we can see that fuzzy objective value is always large than robust objective value, which illustrates that the robust model under box uncertainty is more conservative than the fuzzy model. On the other hand, we choose the optimal solution of robust model under box uncertainty corresponding to the objectives F10 , F20 and F30 , and submit them to the fuzzy optimization model. We find that this solution is always feasible for our fuzzy optimization model, i.e., the solution can ensure the validness of credibility constraints in the fuzzy optimization model. From this observation, we conclude that the obtained optimal solutions to the robust model under box uncertainty set are also the solutions but not the optimal solution for fuzzy model with fixed distributions. Accordingly, if the possibility distributions are available, fuzzy model can be regarded as an improvement of robust model under box uncertainty set, Obviously, fuzzy model can provide more substantial and better decision-making.

6Conclusions

In this paper we introduced a robust multi-objective programming model to study the sustainable development over economic development, electricity consumption, GHG emissions and workforce allocation under uncertain environment. Specifically, we allocated labor across economic sectors to simultaneously satisfy GDP growth, electricity consumption, GHG emissions optimization objectives. After that, we utilized ɛ-constraint method to transform the multi-objective model into single-objective model. For uncertainty involved in the problem, we assumed that the per capita contribution on sectorial GDP, electricity consumption, GHG emissions and the rate of unemployment were characterized with robust uncertainty. That is, these uncertain parameters were known to belong to the given uncertainty sets. However, it would give rise to an overly conservative or computationally intractable robust model along with poorly uncertainty sets. To address the uncertainty and obtain the tractable model, we derived the proposed model to its robust counterpart forms under box+ellipsoidal and box+generalized budgeted uncertainty sets.

We implemented case study based on the UAE to demonstrate the applicability of the proposed model. The computational results showed a mathematical validation to satisfy conflicting objectives on GDP growth, electricity consumption and GHG emissions to achieve UAE’ sustainable development for year 2030. The results indicated that the method yielded good practical efficacy, which were different from previous studies on this topic. Decision makers can commendably make strategic planning and investment allocations towards sustainable development, and can also adjust their parameters according to their preferences through leveraging this method. In addition, the results provided a quantitative decision tool and advice for policy makers on investment and allocation in economic sectors to achieve the sustainability of year 2030. Furthermore, we formulated a fuzzy model for UAE when the possibility distribution of uncertain input data was available.

Future research of this model presented can incorporate additional sustainability criteria on: (i) water resources management and its impact on energy-environment system, (ii) energy saving, recycling and green energy. The model developed in the case study only analysed 8 sectors across different industries. Hence, it would be another meaningful attempt to extend the model to be able to simulate the energy consumption and GHG emissions for more sectors. For uncertainty, in this paper we characterized uncertain per capita GDP, per capita electricity consumption, per capita GHG emission and the unemployment rate by finite closed interval and valivated a robust model for planning the sustainability of UAE in case study. If the possibility distribution of uncertain input data is available, a credibilistic optimization is a better choose for formulating and solving. For the further work we will suggest to use other fuzzy techniques [3, 4, 22] or stochastic methods to address uncertainty.

Acknowledgments

This work are supported partially by the National Natural Science Foundation of China (Nos. 61773150, 61374184) and the Youth Natural Science Foundation of Hebei Province (No. A2016204057). Ms. Ruru Jia is supported by the Post-graduate’s Innovation Fund Project of Hebei University (No. hbu2019ss029).

References

[1]	Asif M. , Sharma R.B. and Adow A.H.E. , An empirical investigation of the relationship between economic growth, urbanization, energy consumption, and CO₂ emission in GCC countries: A panel data analysis, Asian Social Science 11: (21) ((2015) ), 270–284.
[2]	Aznar-Marquez J. and Ruiz-Tamarit J.R. , Environmental pollution, sustained growth, and suffcient conditions for sustainable development, Economic Modelling 54: ((2016) ), 439–449.
[3]	Bai X. , Optimal decisions for prepositioning emergency supplies problem with type-2 fuzzy variables, Discrete Dynamics in Nature and Society, Article ID: 9275192 ((2016) ), 17 pages.
[4]	Bai X. and Liu Y.K. , CVAR reduced fuzzy variables and their second order moments, Iranian Journal of Fuzzy Systems 12: (5) ((2015) ), 45–75.
[5]	Bai X. , Zhang F. and Liu Y.K. , Modeling fuzzy data envelopment analysis under robust input and output data, RAIRO-Operations Research 52: (2) ((2018) ), 619–643.
[6]	Barreto R.A. , Fossil fuels, alternative energy and economic growth, Economic Modelling 75: ((2018) ), 196–220.
[7]	Ben-Tal A. and Newmirovsky A. , Robust convex optimization, Mathematics of Operations Research 23: (4) ((1998) ), 769–805.
[8]	Ben-Tal A. and Newmirovsky A. , Robust solutions of uncertain linear programs, Operations Research Letters 25: (1) ((1999) ), 1–13.
[9]	Bertsimas D. , Pachamanova D. and Sim M. , Robust linear optimization under general norms, Operations Research Letters 32: (6) ((2004) ), 510–516.
[10]	Bertsimas D. and Sim M. , Tractable approximations to robust conic optimization problems, Mathematical Programming 107: (1–2) ((2006) ), 5–36.
[11]	Chapple K. , Planning Sustainable Cities and Regions: Towards More Equitable Development, Routledge, (2014) .
[12]	Dong C. , Huang G.H. , Cai Y.P. , Cheng G.H. and Tan Q. , Bayesian interval robust optimization for sustainable energy system planning in Qiqihar City, China, Energy Economics 60: ((2016) ), 357–376.
[13]	Gabrel V. , Murat C. and Thiele A. , Recent advances in robust optimization: An overview, European Journal of Operational Research 235: (3) ((2014) ), 471–483.
[14]	Guo Z.Z. and Liu Y.K. , Modelling single-period inventory problem by distributionally robust fuzzy optimization method, Journal of Intelligent and Fuzzy Systems 35: (1) ((2018) ), 1007–1019.
[15]	Gupta S. , Fugenschuh A. and Ali I. , A multi-criteria goal programming model to analyze the sustainable goals of India, Sustainability 10: (3) ((2018) ), number.778.
[16]	Jayaraman R. , Colapinto C. , Torre D.L. and Malik T. , A weighted goal programming model for planning sustainable development applied to gulf cooperation council countries, Applied Energy 185: ((2017) ), 1931–1939.
[17]	Jayaraman R. , Colapinto C. , Torre D.L. and Malik T. , Multi-criteria model for sustainable development using goal programming applied to the United Arab Emirates, Energy Policy 87: ((2015) ), 447–454.
[18]	Jayaraman R. , Colapinto C. , Liuzzi D. and Torre D.L. , Planning sustainable development through a scenario-based stochastic goal programming model, Operational Research 17: ((2017) ), 789–805.
[19]	Jayaraman R. , Liuzzi D. , Colapinto C. and Malik T. , A fuzzy goal programming model to analyze energy, environmental and sustainability goals of the United Arab Emirates, Annals of Operations Research 251: ((2017) ), 255–270.
[20]	Jia R.R. and Bai X.J. , Robust optimization approximation of amiguous P-model ang its application, Mathematical Problems in Engineering Article ID: 5203127 ((2018) ), 14 pages.
[21]	Liu B.D. and Liu Y.K. , Expected value of fuzzy variable and fuzzy expected value models, IEEE Transactions on Fuzzy Systems 10: (4) ((2002) ), 445–450.
[22]	Liu N.Q. , Chen Y.J. and Liu Y.K. , Optimization portfolio selection problems under credibilistic CVaR criterion, Journal of Intelligent and Fuzzy Systems 34: ((2018) ), 335–347.
[23]	Liu Y.K. , Bai X.J. , Yang K. , Parametric Credibilistic Optimization Methods, Science Press, Bejing, (2017) .
[24]	Majewski D.E. , Wirtz M. , Lampe M. and Bardow A. , Robust multi-objective optimization for sustainable design of distributed energy supply systems, Computers and Chemical Engineering 102: ((2017) ), 26–39.
[25]	Mezher T. , Chedid R. and Zahabi W. , Energy resource allocation using multi-objective goal programming: The case of Lebanon, Applied Energy 61: (4) ((1998) ), 175–192.
[26]	Mokri A. , Ali M.A. and Emziane M. , Solar energy in the United Arab Emirates: A review, Renewable and Sustainable Energy Reviews 28: ((2013) ), 340–375.
[27]	Nomani M.A. , Ali I. , Fugenschuh A. and Ahmed A. , A fuzzy goal programming approach to analyse sustainable development goals of India, Applied Economics Letters 24: (7) ((2017) ), 443–447.
[28]	Olawumi T.O. and Chan D.W.M. , A scientometric review of global research on sustainability and sustainable development, Journal of Cleaner Production 183: ((2018) ), 231–250.
[29]	Omri A. , CO₂ emissions, energy consumption and economic growth nexus in MENA countries: Evidencec from similtaneous equations models, Energy Economics 40: ((2013) ), 657–664.
[30]	Qiu D. , Lu C. , Zhang W. and Lan Y. , Algebraic properties and topological properties of the quotient space of fuzzy numbers based on Mares equivalence relation, Fuzzy Sets and Systems 245: ((2014) ), 63–82.
[31]	Qiu D. and Zhang W.Q. , Symmetric fuzzy numbers and additive equivalence of fuzzy numbers, Soft Computing 17: ((2013) ), 1471–1477.
[32]	Ramathan R. and Ganesh L.S. , Energy resource allocation incorporating qualitative and quantitative criteria: An integrated model using goal programming and AHP, Socio-Economic Planning 29: (3) ((1995) ), 197–218.
[33]	Rodriguez-Gonzalez P.T. , Rico-Martinez R. and Rico-Ramirez V. , An integrated stochastic economic-ecological-social model with stratified-population, Ecological Modelling 368: ((2018) ), 15–26.
[34]	Shahbaz M. and Lean H.H. , Does financial development increase energy consumption? The role of industrialization in Tunisia, Energy Policy 40: ((2012) ), 473–479.
[35]	Soyster A.L. , Convex programming with set-inclusive constraints and applications to inexact linear programming, Operations Research 21: (5) ((1973) ), 1154–1157.
[36]	Sun B. , Liu Y. and Yang G.Q. , A robust pharmaceutical R&D project portfolio optimization problem under cost and resource uncertainty, Journal of Uncertain Systems 11: (3) ((2017) ), 205–220.
[37]	Yin F.H. and Chen Y.J. , The robust optimization for p-hub Median problem under carbon emissions uncertainty, Journal of Uncertain Systems 12: (1) ((2018) ), 22–35.
[38]	Yoshimura M. , System Design Optimization for Product Manufacturing, Springer Publishing Company, Incorporated, (2010) .
[39]	Zadeh L.A. , Fuzzy set as a basis for a theory of possibility, Fuzzy Sets and Systems 1: ((1978) ), 3–28.