Next Article in Journal
A New Hybrid Approach Using the Simultaneous Perturbation Stochastic Approximation Method for the Optimal Allocation of Electrical Energy Storage Systems
Next Article in Special Issue
Characterization of a Thermoelectric Generator (TEG) System for Waste Heat Recovery
Previous Article in Journal
An Optimization Method for Local Consumption of Photovoltaic Power in a Facility Agriculture Micro Energy Network
Previous Article in Special Issue
Exploiting Game Theoretic Based Coordination Among Appliances in Smart Homes for Efficient Energy Utilization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Criteria Decision Making (MCDM) Approaches for Solar Power Plant Location Selection in Viet Nam

by
Chia-Nan Wang
1,2,*,
Van Thanh Nguyen
1,3,*,
Hoang Tuyet Nhi Thai
3 and
Duy Hung Duong
3
1
Department of Industrial Engineering and Management, National Kaohsiung University of Science and Technology, Kaohsiung 80778, Taiwan
2
Department of Industrial Engineering and Management, Fortune Institute of Technology, Kaohsiung 83160, Taiwan
3
Department of Industrial Systems Engineering, CanTho University of Technology, Can Tho 900000, Viet Nam
*
Authors to whom correspondence should be addressed.
Energies 2018, 11(6), 1504; https://doi.org/10.3390/en11061504
Submission received: 21 May 2018 / Revised: 5 June 2018 / Accepted: 6 June 2018 / Published: 8 June 2018
(This article belongs to the Special Issue Energy Economy, Sustainable Energy and Energy Saving)

Abstract

:
The ongoing industrialization and modernization period has increased the demand for energy in Viet Nam. This has led to over-exploitation and exhausts fossil fuel sources. Nowadays, Viet Nam’s energy mix is primarily based on thermal and hydro power. The Vietnamese government is trying to increase the proportion of renewable energy. The plan will raise the total solar power capacity from nearly 0 to 12,000 MW, equivalent to about 12 nuclear reactors, by 2030. Therefore, the construction of solar power plants is needed in Viet Nam. In this study, the authors present a multi-criteria decision making (MCDM) model by combining three methodologies, including fuzzy analytical hierarchy process (FAHP), data envelopment analysis (DEA), and the technique for order of preference by similarity to ideal solution (TOPSIS) to find the best location for building a solar power plant based on both quantitative and qualitative criteria. Initially, the potential locations from 46 sites in Viet Nam were selected by several DEA models. Then, AHP with fuzzy logic is employed to determine the weight of the factors. The TOPSIS approach is then applied to rank the locations in the final step. The results show that Binh Thuan is the optimal location to build a solar power plant because it has the highest ranking score in the final phase of this study. The contribution of this study is the proposal of a MCDM model for solar plant location selection in Viet Nam under fuzzy environment conditions. This paper also is part of the evolution of a new approach that is flexible and practical for decision makers. Furthermore, this research provides useful guidelines for solar power plant location selection in many countries as well as a guideline for location selection of other industries.

1. Introduction

The Earth is facing global warming and climate change challenges. Many studies have been done to find solutions by exploiting renewable energy. This is the best way of reducing fossil fuel use, reducing greenhouse gas emissions and maintaining the Earth’s temperature increase under 2 °C [1].
Solar energies, including solar energy, photovoltaic, and solar thermal, have many positive effects on the environment, contributing to the sustainable development of society and improving the quality of human life [2]. Nowadays, building solar power plants is becoming easier because the price of solar panels is decreasing [3,4]. The advantages of solar energies are increased CO2 mitigation, they do not make noise, they minimize toxic wastes and they do not require environmental remediation treatments [5]. In this study, we introduced a MCDM approach including DEA, Fuzzy AHP and TOPSIS to select the best location for building a solar power plant.
Many studies have applied the MCDM approach to various fields of science and engineering and their number has been increasing over the past years. One of the fields where the MCDM model has been employed is the location selection problem. Location selection is an important use of MCDM models. Huang et al. [6], Loken [7] used MCDM models for selecting construction locations in the energy sector. The general procedure of MCDM is shown in Figure 1 [8].
The AHP model was proposed by Saaty [9] in the 1980s.There are six steps in the AHP procedure as follows:
(1)
Specifying the problem;
(2)
Constructing the AHP hierarchy;
(3)
Building a pairwise comparison matrix
(4)
Defining the weight of factors.
(5)
Checking Consistency Index.
(6)
Obtaining the overall rating and making decision
The AHP model has many advantages, however, the AHP model cannot accommodate uncertainty and inaccuracies between the perceptions and judgment of the decision makers. Thus, the AHP model with fuzzy logic is proposed to address this problem. In the FAHP model, decision makers can approximate input data by using fuzzy numbers. As with capacity planning, decision makers need to follow a four step procedure when making location selection. These steps are as shown in the following Figure 2 [10,11]:
Stage 1. Identifying location criteria. In this stage, all criteria that affect a business will be defined.
Stage 2. Developing site alternatives. Once decision makers know what criteria affect a business, they can identify location options that satisfy the selected criteria.
Stage 3. Evaluating location options. After a set of location options are defined, decision makers will evaluate and rank options by quantitative or qualitative methods.
Stage 4: Making final decisions. The best location with the highest ranking score will be selected.
Azadeh et al. [12] proposed a hybrid MCDM model including DEA, PCA and NT for selecting solar power plant sites. A. Azadeh et al. [13] also presented a hybrid ANN and fuzzy DEA approach to select solar power plant locations. Lee et al. [14] introduced a MCDM model to select PV solar plant locations. Ali et al. combined GIS and MCDM approaches to determine the best place for wind farm location [15].
Gao et al. [16] determined the best enterprise location by using AHP and DEA models. Yang and Kuo [17] proposed an analytic hierarchy process and data envelopment analysis model for location selection. Kabir and Hasin [18] proposed a hybrid fuzzy analytic hierarchy process (FAHP) and PROMEETHE for locating power substation sites. Lee, Kang and Liou [19] proposed a hybrid model including ISM, FANP and VIKOR to select the most suitable PV solar plant site. Suh and Brownson applied GISFS and AHP approaches to select PV solar plant locations [20].
Noorollahi et al. [21] used GIS and FAHP for land analyses in solar farm location. Gan et al. [22] analyzed economic feasibility for renewable energy projects by using integrated TFN-AHP-DEA approaches. Liu et al. proposed a hybrid MCDM model for evaluating the total factor energy efficiency by combining DEA and the Malmquist index in the thermal power industry [23].
Samanlioglu and Ayag˘ [24] used the FAHP and F-PROMETHEE II for solar power plant location selection. Nazari, Aslani and Ghasempour [25] proposed TOPSIS approaches for analysising of solar farm site selection options. Al Garni and Awasthi [26] selected the best location for utility scale solar PV projects by using GIS and an MCDM approach. Merrouni et al. [27] used GIS and the AHP to assess the capacity of Eastern Morocco to host large-scale PV farms. Lozano, García-Cascales and Lamata [28] proposed a comparative TOPSIS-ELECTRE TRI method for photovoltaic solar farm site selection. Beltran et al. [29] used ANP for selection of photovoltaic solar power project sites.
The remainder of the paper provides background materials to assist in developing the MCDM model. Then, a hybrid DEA-FAHP-TOPSIS approach is proposed to select the best location for construction of a solar power plant from among 46 potential locations in Viet Nam. Discussions and the main contributions of this research are presented at the end of this article.

2. Material and Methodology

2.1. Research Development

In this study, the authors present a MCDM model including DEA, Fuzzy AHP and TOPSIS to select the best location for building a solar power plant in Viet Nam. There are four steps in our research, as shown in Figure 2:
Step 1: Determining evaluate criteria. In this step, the criteria for selecting the best location will be defined. The key criteria and sub-criteria have built through expert interviews and the results from others’ research. All of the criteria are shown in Figure 3.
Step 2: Employing the DEA model. There are 46 location options that can be highly effective for a solar power plant construction. In this step, several DEA methods including the CCR model, BCC model, and SBM model are applied to rank all options. The options that reach EFF = 1 in all models are potential locations and will be considered in the next step.
Step 3: Applying FAHP model. The FAHP model is the most effective tool for addressing complex problems of decision making with a connection to various qualitative criteria. The weight of criteria will be defined in this step.
Step 4: Implementing the TOPSIS model. The TOPSIS model is employed to rank potential locations. The optimal options have the shortest geometric distance from the positive ideal solution (PIS) and the longest geometric distance from the negative ideal solution (NIS). The best potential site will be presented in this stage.

2.2. Methodology

2.2.1. Data Envelopment Analysis Model

(1) Charnes-Cooper-Rhodes model (CCR model)
The basis DEA model is CCR model [30], CCR model is defined as follows:
max c . a ξ = a V y 0 c V x 0 S . t : a V y e   c V x e   0 ,   e = 1 ,   2 , , n a 0 c 0
Constraints mean that the ratio of virtual output to virtual input cannot exceed 1 per DMU. The goal is obtain a rate of weighted output for weighted inputs. Due to constraints, the optimal goal value ξ* is at most 1.
DMU0 is CCR efficient if ξ * = 1 and the result must have at least 1 optima a* > 0 and c* > 0. In addition, the fractional program can be described as a linear program (LP) as follows [31]:
max c . a ξ = a v y 0 S . t : c v x 0 1 = 0 a v y e c v x e 0 ,   e = 1 ,   2 , . . , n c 0 a 0
The fractional program (1) is equal to the linear program (2) [32]. The Farrell model of linear program (2) with variable ξ and a nonnegative vector α = α 1 , α 2 , α 3 , , α f as [31]:
m a x b = 1 m s i + r = 1 q s r + S . t : e = 1 n x b e α b + s b = ξ x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q α e 0 ,   e = 1 ,   2 , , n s b 0 , b = 1 , 2 , , p s r + 0 , r = 1 , 2 , , q
The model (3) has a feasible solution, ξ = 1 , α 0 * = 1 , α j * = 0 , ( j 0 ) ,   which   affects optimal value when ξ * is not greater than 1. The optimal solution, ξ * , provides an effective point for a specific DMU. The process will be repeated for each DMUe, e = 1, 2,…, n. DMUs are inefficient when ξ * < 1, while DMUs are boundary points if ξ * = 1 . We can avoid weakly efficient frontier points by invoking a linear program as follows [31]:
m a x b = 1 m s b + r = 1 s s r + S . t : e = 1 n x b e α e + s b = ξ x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q α e 0 ,   e = 1 ,   2 , , n s b 0 , b = 1 , 2 , , p s r + 0 , r = 1 , 2 , , q
In this case, we note that the choices of s b and s r + do not affect the optimal ξ * .
DMU0 achieves 100% efficiency if and only if both (1) ξ = 1 and (2) s b * = s r + = 0 . The performance of DMU0 is weakly efficient if and only if both (1) ξ * = 1 and (2) s b * 0 and s r + 0 for i or r in optimal options. Thus, the preceding development amounts to solving the problem as follows [31]:
m i n θ μ ( b = 1 m s b + r = 1 s s r + ) S . t : b = 1 n x b e α e + s b = ξ x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q α e 0 ,   e = 1 ,   2 , , n s b 0 , b = 1 , 2 , , p s r + 0 , r = 1 , 2 , , q
In this case, s b and s r + variables will be used to convert the inequalities into equivalent equations. This is similar to solving (3) by minimizing ξ in first stage and then fixing 𝜉 = ξ * as in (4), where the slacks variables achieve a maximum value but do not affect to previously determined value of ξ = ξ * . The objective will be converted from max to min, as in (1), to obtain [31]:
max c . a ξ = c V x 0 a V y e S . t : a V x 0 c V y e ,   e = 1 ,   2 , , n c ε > 0 a ε > 0
If the ε > 0 and the non-Archimedean element is defined, the input models are similar to model (2) and (5) as follows [31]:
max c . a ξ = c V x 0 S . t : a V y 0 = 1 c V x o a V y e 0 ,   e = 1 ,   2 , , n c ε > 0 a ε > 0
and:
m a x ϕ ε ( b = 1 m s i + r = 1 s s r + ) S . t : e = 1 n x b e α e + s b = x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q α e 0 ,   e = 1 ,   2 , , n s b 0 , b = 1 , 2 , , p s r + 0 , r = 1 , 2 , , q
The CCR input-oriented (CCR-I) has the dual multiplier model of is expressed as [31]:
max z = r = 1 q r y r 0 S . t : r = 1 q r y r e r = 1 q a r y r e 0 b = 1 p a b x b 0 = 1 c r , a b ε > 0
The CCR output-oriented (CCR-O) has the dual multiplier model of is expressed as [31]:
min q = b = 1 p a b x b 0 S . t : b = 1 p a b x b e r = 1 q r y r e 0 r = 1 q r y r 0 = 1 c r , a b ε > 0
(2) Banker Charnes Cooper model (BCC Model)
Banker et al. introduced input-oriented BBC model (BCC-I) [30], which is able to assess the efficiency of DMU0 by solving the following linear program (11) [31]:
ξ B = m i n ξ S . t : e = 1 n x b e α e + s b = ξ x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q k = 1 n α k = 1 α k 0 , k = 1 , 2 , , n
We can avoid the weakly efficient frontier points by invoking a linear program as follows [31]:
m a x b = 1 m s b + r = 1 s s r + S . t : e = 1 n x b e α e + s b = ξ x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , , n s b 0 , b = 1 , 2 , , p s r + 0 , r = 1 , 2 , , q
Therefore, this is the first multiplier form to the solve problem as follows [31]:
m i n ξ ε ( b = 1 m s b + r = 1 s s r + ) S . t : e = 1 n x b e α e + s b = ξ x i 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = y r 0 ,   r = 1 ,   2 , , q k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , , n d b 0 , b = 1 , 2 , , p d r + 0 , r = 1 , 2 , , q
The linear program (12) gives us the second multiplier form, which is expressed as [31]:
max c . a , a 0 ξ B = a V y 0 a 0 S . t : c V x 0 = 1 a V y e c V x e a 0 0 ,   e = 1 ,   2 , , n c 0 a 0
In this case v and u, which are mentioned in the Formula (14), are vectors, and the scalar v 0 may be positive or negative or zero. Thus, the equivalent BCC fractional program is got from the dual program (14) as [31]:
max c . a ξ = a V y 0 a 0 c V x 0 S . t : a V y e a 0 c V x e 1 ,   e = 1 ,   2 , , n c 0 a 0
The DMU0 can be called BCC efficient if an optimal solution ( ξ B * ,   s * , s + * )   as claimed in this two phase processes for model (9) satisfies ξ B * = 1 and has no slack s * =   s + * = 0. On the other hand, it is BCC inefficient.
The improved activity ( ξ * x   s * , y + s + * ) also can be illustrated as BCC efficient [31]. A DMU, which has a minimum input value for any input item, or a maximum output value for any output item, is BCC-efficient.
The output-oriented BCC model (BCC-O) is:
max η S . t : e = 1 n x b e α e + s b = ξ x b 0 ,   b = 1 ,   2 , , p e = 1 n y r e α e s r + = η y r 0 ,   r = 1 ,   2 , , q k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , , g
From the linear program (16), we have the associate multiplier form, which is expressed as [31]:
min c . a , c 0 a V y 0 a 0 S . t : a V y 0 = 1 c V x e a V y e a 0 0 ,   e = 1 ,   2 , , n c 0 a 0
In the envelopment model, the v0 is the scalar combined with k = 1 n α k = 1 . In conclusion, the authors achieve the equivalent (BCC) fractional programming formulation for model (17) [31]:
min c . a , c 0 c V x 0 c 0 a V y 0 S . t : c V x e a 0 a V y e 1 ,   e = 1 ,   2 , , n c 0 a 0
(3) Slacks Based Measure model (SBM Model)
The SBM model is developed by Tone [33,34]. It has three elements, i.e., input-oriented, output-oriented, and non-oriented.

Input-Oriented SBM (SBM-I-C)

Input-oriented SBM under constant-returns-to-scale-assumption [31]:
ρ I * =   min α ,     s , s + 1   1 m b = 1 m s b x b h S . t : x b c = e = 1 m x b c α b +   s b   ,   b = 1 ,   2 ,   p y r c = e = 1 m y r c α e   s r +   ,   r = 1 ,   2 ,   q α e 0 ,   k   ( j ) ,   s b 0   ( e ) , s r + 0   ( e )
is called the SBM input efficiency.

Output-Oriented SBM (SBM-O-C)

The output-oriented SBM efficiency ρ O * of DMUz = (xz, yz) is defined by [SBM-O-C] [32]:
1 ρ O * = max α , s , s + 1 + 1 s r = 1 q s r + y r z S . t : x b z = e = 1 n x b e α e +   s e ( b = 1 , . . p ) y b z = e = 1 n y b e α e + s b + ( b = 1 , p ) α e 0 ( e ) ,   s e 0 ( b ) ,   s e + 0   ( r )

2.2.2. Fuzzy Analytic Hierarchy Process (FAHP)

(1) Fuzzy sets and fuzzy number
Fuzzy sets were proposed by Zadeh in 1965 [35] for solving problems existing in uncertain environments. A fuzzy set is a function that shows a degree of dependence of one fuzzy number on a set number, where each value of the membership function is between [0, 1] [36,37]. The triangular fuzzy number (TFN) can be defined as (l, m, u). The parameters l, m and u (lmu), indicate the smallest, the promising and the largest value. A triangular fuzzy number are shown in Figure 4.
TFN can be defined as:
μ ( x M ˜ ) = { 0 , x l m l u x u m 0 ,     x < l , l x m , m x u , x > u ,
A fuzzy number is given by the representatives of each level of membership are the following:
M ˜ = ( M l ( y ) , M r ( y ) ) = [ l + ( m l ) y ,   u + ( m u ) y ] , y   [ 0 , 1 ]
l(y), r(y) indicates both the left side and the right side of a NF. Two positive triangular fuzzy numbers (l1, m1, u1) and (l1, m2, u2) are introduced as below:
( l 1 ,   m 1 ,   u 1 ) + ( l 2 ,   m 2 ,   u 2 ) = ( l 1 + l 2 ,   m 1 + m 2 , u 1 + u 2 ) ( l 1 ,   m 1 ,   u 1 ) ( l 2 ,   m 2 ,   u 2 ) = ( l 1 l 2 ,   m 1 m 2 , u 1 u 2 ) ( l 1 ,   m 1 ,   u 1 ) × ( l 2 ,   m 2 ,   u 2 ) = ( l 1 × l 2 ,   m 1 × m 2 , u 1 × u 2 ) ( l 1 ,   m 1 ,   u 1 ) ( l 2 ,   m 2 ,   u 2 ) = ( l 1 / u 2 ,   m 1 / m 2 , u 1 / l 2 )
(2) Analytical hierarchy process (AHP)
The AHP method uses a pairwise comparisons maxtrix for determing priorities on each level of the hierarchy that are quantified using a 1–9 scale are shown in Table 1.
(3) Fuzzy AHP
In this section, the weight of criteria are dedined by fuzzy AHP. There are eight steps in this process, as follows:
Step 1: Calculation of TFNs
A pairwise comparison of the criteria will be performed. Instead of a numerical value, the fuzzy analytical hierarchy process is a range of values that are combined for evaluating criteria in this step [38]. This scale is applied in Parkash’s [39] fuzzy prioritization method. The fuzzy conversion scales are shown in Table 2.
Step 2: Calculation of P ˜ 1
A pairwise comparison and relative scores as (24):
P a ˜ = ( l a , m a , u a )
l a = ( l a 1 l a 2 l a i ) 1 i ,   a = 1 ,   2 ,   i
m a = ( m a 1 m a 2 m a i ) 1 i ,   a = 1 ,   2 ,   i
u a = ( u a 1 u a 2 u a i ) 1 i ,   i = 1 ,   2 ,   i
Step 3: Calculation of P ˜ Y
The geometric fuzzy mean was established by (28):
P ˜ Y =   ( a = 1 i l a , a = 1 i m a , a = 1 i u a )
Step 4: Calculation of R ˜
The fuzzy geometric mean was determined as:
R ˜ = P ˜ a P ˜ Y = ( l a , m a , u a ) a = 1 i l a , a = 1 i m a , a = 1 i u i = [ l a a = 1 i u a , m a a = 1 i m a , u a a = 1 i l a ]
Step 5: Calculation of W a β l
The criteria depending on β cut values are defined for the calculated β. The fuzzy priorities will apply for lower and upper bounds for each β value:
W a β l = ( W a l β l , W a u β l ) ; a = 1 ,   2 ,   i ; l = 1 ,   2 ,   L
Step 6: Calculation of Wal, Wau
Values of Wal, Wau are calculated by combining the lower and the upper values, and dividing them by the total β values:
W a l = a = 1 i β ( W a l ) l l = 1 L β l ; a = 1 ,   2 ,   i ; l = 1 ,   2 ,   L
W a u = a = 1 i β ( W a u ) l l = 1 L β l ; a = 1 ,   2 ,   i ; l = 1 ,   2 ,   L
Step 7: Calculation of Xad
Combining the upper and the lower bounds values by using the optimism index ( γ ) to order to defuzzify:
W a d = γ × W a u + ( 1 γ ) × W a l ; γ   [ 0 , 1 ] ; a = 1 ,   2 , i
Step 8: Calculation of Waz
The defuzzification values priorities are normalization by:
W a z = W a d a = 1 i W a d ; a = 1 , 2 , . . i

2.2.3. Technique for Order Preference by Similarity to Ideal Solution (TOPSIS)

Hwang and Yoon [40] is presented the TOPSIS approach in 1981. The main concept of TOPSIS is that the best options should have the shortest geometric distance from the positive ideal solution (PIS) and the longest geometric distance from the negative ideal solution (NIS) [41]. There are m alternatives and n criteria and the result of TOPSIS model shows the score of each option [42]. The method is illustrated below:
Step 1: Determine the normalized decision matrix, raw values (xij) are converted to normalized values (nij) by:
h a b = y a b a g y a b 2 ,   a = 1 , g ; b = 1 , . . h .
Step 2: Calculate the weight normalized value (vij), by:
l a b = P a b h a b ,   a = 1 , . , g ; b = 1 , , h .
where Pj is the weight of the a t h criterion and b = 1 h p p = 1 .
Step 3: Calculate the PIS ( B + ) and PIS ( B ), where l a + indicate the maximum values of l a b and l a indicates the minimum value l a b :
B + = { l 1 + , , l h + } = { ( max b l a b | a A ) , ( min b l a b | a A ) } ,
B = { l 1 , , l n } = { ( min b l a b | a A ) , ( max b l a b | b B ) } ,
where A is related with profit criteria, and B is related with cost criteria.
Step 4: Determine a distance of the PIS ( S a + ) separately by:
S a + = { b = 1 h ( l a b l b + ) 2 } 1 2 ,   a = 1 ,   . ,   g
Similarly, the separation from the NIS ( S i ) is given as:
S a = { b = 1 h ( l a b l b ) 2 } 1 2 ,   a = 1 ,   . ,   g
Step 5: Determine the relationship proximal to the problem solving approaches, proximal relationship from option B a to option B + :
C a = S a S a + + S a ,   a = 1 , , g .
Step 6: Rank alternatives to determine the best option with the maximum value of C a

3. Case Study

According to the research results of many scientists, Viet Nam is the best place with natural conditions and favorable terrain to develop renewable energy. Viet Nam is located in the tropical monsoon region, the average number of sunshine hours in the year ranges from 2500 to 3000 h and the average temperature of over 21 °C. In addition, Viet Nam has abundant solar radiation sources. Viet Nam’s solar map is shown in Figure 5.
The authors collected data from 46 potential sites, which are able to invest in solar power plants as shown in Table 3.
DEA is a mathematical programming technique that determines the relative effectiveness of multiple input and output decision makers (DMUs) [43]. There are two inputs, two outputs in including Temperature, Wind Speed, Sunshine hours, Elevation [14]. Inputs and Outputs of DMUs are shown in Figure 6.
Some additional data about the 46 locations are shown in Table 4.
For selecting the best potential location, several DEA models including CCR-I; CCR-O; BCC-I; BCC-O and SBM-I-C are applied in this step. The results of the DEA models are shown in Table 5.
As the results in Table 5 show, there are seven DMUs that are potential locations, including DMU 2, DMU 11, DMU 21, DMU 22, DMU 31, DMU 36 and DMU 38. These DMUs will be evaluated in next step of this research. The FAHP model is applied in this stage. The weight of criteria are defined by the comparison matrix. Criteria structures are built based on qualitative and quantitative factors. The Hierarchical structures for the FAHP approach are shown in Figure 7.
To convert the fuzzy numbers to real numbers we proceed to solve the fuzzy clusters using the triangular fuzzy method. During the defuzzification we obtain the coefficients α = 0.5 and β = 0.5 (Tang and Beynon) [44]. In it, α represents the uncertain environment, β represents the attitude of the evaluator is fair:
g 0.5 , 0.5 ( a E N , E C ¯ ) = [ ( 0.5 × 2.5 ) + ( 1 0.5 ) × 3.5 ] = 3 f 0.5 ( L E N , E C ) = ( 3 2 ) × 0.5 + 2 = 2.5 f 0.5 ( U E N , E C ) = 4 ( 4 3 ) × 0.5 = 3.5 g 0.5 , 0.5 ( a E C , E N ¯ ) = 1 / 3
The remaining calculations are similar to the above, as well as the fuzzy number priority points. The real number priorities when comparing the main criteria pairs are presented in Table 7.
To calculate the maximum individual values as follows:
GM1 = (1 × 1/3 × 2 × 2 × 1/2)1/5 = 0.92
GM2 = (3 × 1 × 3 × 2 × 1)1/5 = 1.78
GM3 = (1/2 × 1/3 × 1 × 3 × 1)1/5 = 0.87
GM4 = (1/2 × 1/2 × 1/3 × 1 × 1/2)1/5 = 0.53
GM5 = (2 × 1 × 1 × 2 × 1)1/5 = 1.32
GM = GM1 + GM2 + GM3 + GM4 + GM5 = 5.42
ω 1 = 0.92 5.42 = 0.17 ω 2 = 1.78 5.42 = 0.33 ω 3 = 0.87 5.42 = 0.16 ω 4 = 0.53 5.42 = 0.1 ω 4 = 1.32 5.42 = 0.24
[ 1 1 / 3 2 2         1 / 2 3 1 3 2           1 1 / 2 1 / 3 1 3       1 1 / 2 1 / 2 1 / 2 1       1 / 2 2 1 1 2 1 ] × [ 0.17 0.33 0.16 0.1 0.24 ] = [ 0.92 1.76 0.9 0.55 1.27 ] [ 0.92 1.76 0.9 0.55 1.27 ] / [ 0.17 0.33 0.16 0.1 0.24 ] = [ 5.41 5.33 5.63 5.5 5.29 ]
With the number of criteria is 5, we get n = 5, λmax and CI are calculated as follows:
λ m a x = 5.41 + 5.33 + 5.63 + 5.5 + 5.29 5 = 5.432
C I = 5.43 5 5 1 = 0.1075
For CR, with n = 5 we get RI = 1.12:
C R = 0.1075 1.12 = 0.09598
We have CR = 0.09598 ≤ 0.1, so the pairwise comparison data is consistent and does not need to be re-evaluated. The results of the pair comparison between the main criteria are presented in Table 8.
The comparison matrix of sub-criteria based on the alternatives is shown below:
The final weight of the criteria are shown in Table 26.
The weights of alternative locations with respect to all sub criteria and the Normalized Decision Matrix of sub criteria are shown in Table 27, Table 28 and Table 29.
In the final stage, all the potential locations will be ranked by the TOPSIS model. The weight of sub-criteria can be used from the result of the fuzzy AHP approach. The normalized weight matrix values are shown in Table 30.
In Figure 8, DMU 2 has shortest geometric distance from the PIS and the longest geometric distance from the NIS.
The results of TOPSIS model are shown in Figure 9, based on the final performance score C a , the final location ranking list are DMU 2, DMU 38, DMU 21, DMU 22, DMU 11, DMU 31 and DMU 36. The results show that DMU 2 (Binh Thuan) is the most optimal location for building a solar plant in Viet Nam.

4. Conclusions

Studies applying the MCDM approach to various fields of science and engineering have been increasing in number over the past years. One of the fields where the MCDM model has been employed ii in location selection problems. Especially in the renewable energy sector, decision makers have to evaluate both qualitative and quantitative factors. Although some studies have reviewed applications of MCDM approaches in solar power plant location selection, very few works has focused on this problem in a fuzzy environment. This is a reason why in this work MCDM model including AHP with fuzzy logic, DEA and TOPSIS is proposed for solar power plant location selection in Viet Nam. The goal of this study was to design a MCDM approach for building solar power plants based on natural and social factors. In the first step of the research, proper areas are defined by using several DEA models, then a fuzzy analytical hierarchy process is proposed for evaluating the weight of criteria. The FAHP can be applied for ranking alternatives but the number of sites selected is practically limited because of the number of pairwise comparisons that need to be made, and a disadvantage of the FAHP approach is that the input data, expressed in linguistic terms, depend on the experience of decision makers and thus involves subjectivity. This is a reason why we proposed the TOPSIS model for ranking alternatives in the final stage. Also, TOPSIS is presented to reaffirm it as a systematic method and solve the disadvantages of the FAHP model as mentioned above. As a results, the site with the best potential DMU 2 (Binh Thuan) because it has the highest ranking score in the final stage.
The contribution of this study is the presentation of a multi-criteria decision making model (MCDM) for solar plant site selection in Viet Nam under fuzzy environment conditions. This paper also represents the evolution of a new approach that is flexible and practical for the decision maker. This research also provides a useful guideline for solar power plant location selection in other countries as well as a guideline for location selection in other industries.
In future research, this MCDM model also can be applied to many different countries. In addition, different methods, such as FANP or PROMETHEE, etc., could also be combined for different scenarios.

Author Contributions

In this research, C.N.-W. built the research ideas, designed the theoretical verifications, and reviewed manuscript. V.T.N. contributed the research ideas, designed the frameworks, collected data, analyzed the data, summarized and wrote the manuscript. H.T.N.T. collected data, write a manuscript, D.H.D. collected data, wrote and format manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. United Nations. Adoption of the Paris Agreement. In Proceedings of the Conference of the Parties, Paris, France, 30 November–11 December 2015. [Google Scholar]
  2. Tsoutsos, T.; Frantzeskaki, N.; Gekas, V. Environmental impacts from the solar energy technologies. Energy Policy 2005, 33, 289–296. [Google Scholar] [CrossRef] [Green Version]
  3. Candelise, C.; Winskel, M.; Gross, R.J. The dynamics of solar PV costs and prices as a challenge. Renew. Sustain. Energy Rev. 2013, 26, 96–107. [Google Scholar] [CrossRef]
  4. Honrubia-Escribano, A.; Javier Ramirez, F.; Gómez-Lázaro, E.; Garcia-Villaverde, P.M.; Ruiz-Ortega, M.J.; Parra-Requena, G. Influence of solar technology in the economic performance of PV power. Renew. Sustain. Energy Rev. 2018, 82, 488–501. [Google Scholar] [CrossRef]
  5. Alsema, E.A. Energy requirements and CO2 mitigation potential of PV systems. In Proceedings of the BNL/NREL Workshop, Colorado, CO, USA, 23–24 July 1998. [Google Scholar]
  6. Huang, J.P.; Poh, K.L.; Ang, B.W. Decision analysis in energy and environmental modeling. Energy 1995, 20, 843–855. [Google Scholar] [CrossRef]
  7. Løken, E. Use of multicriteria decision analysis methods for energy planning problems. Renew. Sustain. Energy Rev. 2007, 11, 1584–1595. [Google Scholar] [CrossRef]
  8. Belton, V.; Stewart, T. Multiple Criteria Decision Analysis: An Integrated Approach; Kluwer Acedemic: Norwell, MA, USA, 2002. [Google Scholar]
  9. Saaty, T.L. The Analytic Hierarchy Process: Planning, Priority Setting, Resources Allocation; McGraw-Hill: New York, NY, USA, 1980. [Google Scholar]
  10. Stevenson, W.J. Operations Management; McGraw-Hill: New York, NY, USA, 2002. [Google Scholar]
  11. Sanders, R.D.R.a.N.R. Operations Management, 4th ed.; Wiley: Hoboken, NJ, USA, 2009. [Google Scholar]
  12. Azadeh, A.; Ghaderi, S.; Maghsoudi, A. Location optimization of solar plants by an integrated hierarchical. Energy Policy 2008, 36, 3993–4004. [Google Scholar] [CrossRef]
  13. Azadeh, A.; Sheikhalishahi, M.; Asadzadeh, S. A flexible neural network-fuzzy data envelopment analysis approach for location optimization of solar plants with uncertainty and complexity. In Proceedings of the 2015 International Conference on Industrial Engineering and Operations Management, Dubai, UAE, 3–5 March 2015. [Google Scholar]
  14. Lee, A.H.I.; Kang, H.-Y.; Lin, C.-Y.; Shen, K.-C. An Integrated Decision-Making Model for the Location of a PV. Sustainability 2015, 7, 13522–13541. [Google Scholar] [CrossRef]
  15. Ali, S.; Lee, S.-M.; Jang, C.-M. Determination of the Most Optimal On-Shore Wind Farm Site Location Using a GIS-MCDM Methodology: Evaluating the Case of South Korea. Energies 2017, 10, 2072. [Google Scholar] [CrossRef]
  16. Gao, Z.; Yoshimoto, K.; Ohmori, S. Application of AHP/DEA to Facility Layout Selection. In Proceedings of the Third International Joint Conference on Computational Science and Optimization, Huangshan, China, 28–31 May 2010. [Google Scholar]
  17. Yang, T.; Kuo, C. A hierarchical AHP/DEA methodology for the facilities layout design problem. Eur. J. Oper. Res. 2003, 147, 128–136. [Google Scholar] [CrossRef]
  18. Kabir, G.; Hasin, M.A.A. Integrating fuzzy AHP with TOPSIS method for optimal power substation location selection. Int. J. Logist. Econ. Glob. 2013, 5, 312–331. [Google Scholar] [CrossRef]
  19. Lee, A.H.I.; Kang, H.-Y.; Liou, Y.-J. A Hybrid Multiple-Criteria Decision-Making Approach for Photovoltaic Solar Plant Location Selection. Sustainability 2017, 9, 184. [Google Scholar] [CrossRef]
  20. Suh, J.; Brownson, J. Solar Farm Suitability Using Geographic Information System Fuzzy Sets and Analytic Hierarchy Processes: Case Study of Ulleung Island, Korea. Energies 2016, 9, 648. [Google Scholar] [CrossRef]
  21. Noorollahi, E.; Fadai, D.; Shirazi, M.A.; Ghodsipour, S.H. Land Suitability Analysis for Solar Farms Exploitation Using GIS and Fuzzy Analytic Hierarchy Process (FAHP)—A Case Study of Iran. Energies 2016, 9, 643. [Google Scholar] [CrossRef]
  22. Gan, L.; Xu, D.; Hu, L.; Wang, L. Economic Feasibility Analysis for Renewable Energy Project Using an Integrated TFN–AHP–DEA Approach on the Basis of Consumer Utility. Energies 2017, 10, 2089. [Google Scholar] [CrossRef]
  23. Liu, J.-P.; Yang, Q.-R.; He, L. Total-Factor Energy Efficiency (TFEE) Evaluation on Thermal Power Industry with DEA, Malmquist and Multiple Regression Techniques. Energies 2017, 10, 1039. [Google Scholar] [CrossRef]
  24. Samanlioglu, F.; Ayag, Z. A fuzzy AHP-PROMETHEE II approach for evaluation of solar power plant location alternatives in Turkey. J. Intell. Fuzzy Syst. 2017, 33, 859–871. [Google Scholar] [CrossRef]
  25. Nazari, M.A.; Aslani, A.; Ghasempour, R. Analysis of Solar Farm Site Selection Based on TOPSIS Approach. Int. J. Soc. Ecol. Sustain. Dev. 2018, 9, doi–10. [Google Scholar] [CrossRef]
  26. Garni, H.Z.A.; Awasthi, A. Solar PV power plant site selection using a GIS-AHP based approach with application in Saudi Arabia. Appl. Energy 2017, 206, 1225–1240. [Google Scholar] [CrossRef]
  27. Merrouni, A.A.; Elalaoui, F.E.; Mezrhab, A.; Mezrhab, A.; Ghennioui, A. Large scale PV sites selection by combining GIS and Analytical Hierarchy Process. Case study: Eastern Morocco. Renew. Energy 2018, 119, 863–873. [Google Scholar] [CrossRef]
  28. Sanchez-Lozano, J.; García-Cascales, M.; Lamata, M. Comparative TOPSIS-ELECTRE TRI methods for optimal sites for photovoltaic solar farms. Case study in Spain. J. Clean. Prod. 2016, 127, 387–398. [Google Scholar] [CrossRef]
  29. Beltran, P.A.; Gonzalez, F.C.; Ferrando, J.P.; Pozo, F.R. An ANP-based approach for the selection of photovoltaic solar power plant investment projects. Renew. Sustain. Energy Rev. 2010, 14, 249–264. [Google Scholar] [CrossRef]
  30. Charnes, A.; Cooper, W.; Rhodes, E. Measuring the efficiency of decision making units. Eur. J. Oper. Res. 1978, 2, 429–444. [Google Scholar] [CrossRef]
  31. Wen, M. Uncertain Data Envelopment Analysis. In Uncertainty and Operations Research; Springer: Heidelberg, Berlin, 2015. [Google Scholar]
  32. Farrell, M.J. The Measurement of Productive Efficiency. J. R. Stat. Soc. 1957, 120, 253–281. [Google Scholar] [CrossRef]
  33. Tone, K. A slacks-based measure of efficiency in data envelopment analysis. Eur. J. Oper. Res. 2001, 130, 498–509. [Google Scholar] [CrossRef] [Green Version]
  34. Pastor, J.T.; Ruiz, J.L.; Sirvent, I. An enhanced DEA Russell graph efficiency measure. Eur. J. Oper. Res. 1999, 115, 596–607. [Google Scholar] [CrossRef]
  35. Zadeh, L. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef] [Green Version]
  36. Shu, M.S.; Cheng, C.H.; Chang, J.R. Using intuitionistic fuzzy set for fault-tree analysis on printed circuit board assembly. Microelectron. Reliab. 2006, 46, 2139–2148. [Google Scholar] [CrossRef]
  37. Kahraman, Ç.; Ruan, D.; Ethem, T. Capital budgeting techniques using discounted fuzzy versus probabilistic cash. Inf. Sci. 2002, 42, 57–76. [Google Scholar] [CrossRef]
  38. Kuswandari, R. Assessment of Different Methods for Measuring the Sustainability of Forest Management Retno Kuswandari. Master’s Thesis, University of Twente, Enschede, The Netherlands, December 2004. [Google Scholar]
  39. Prakash, T.N. Land Suitability Analysis for Agricultural Crops: A Fuzzy Multi Criteria Decision Making Approach. Master’s Thesis, University of Twente, Enschede, The Netherlands, December 2003. [Google Scholar]
  40. Hwang, C.L.; Yoon, K.; Hwang, C.L.; Yoon, K. Multiple Attribute Decision Making: Method and Applications; Springer: New York, NY, USA, 1981. [Google Scholar]
  41. Assari, A.; Maheshand, T.; Assari, E. Role of public participation in sustainability of historical city: Usage of TOPSIS method. Indian J. Sci. Technol. 2012, 5, 2289–2294. [Google Scholar]
  42. Jahanshahloo, G.R.; Lotfi, F.H.; Izadikhah, M. Extension of the TOPSIS Method for Decision-Making Problems with Fuzzy Data. Appl. Math. Comput. 2006, 181, 1544–1551. [Google Scholar] [CrossRef]
  43. Liu, J.; Ding, F.-Y.; Lall, V. Using data envelopment analysis to compare suppliers for supplier selection and performance improvement. Supply Chain Manag. Int. J. 2000, 5, 143–150. [Google Scholar] [CrossRef]
  44. Tang, Y.-C.; Beynon, M.J. Application and Development of a Fuzzy Analytic Hierarchy Process within a Capital Investment Study. J. Econ. Manag. 2005, 1, 207–230. [Google Scholar]
Figure 1. General methodology of MCDM procedure.
Figure 1. General methodology of MCDM procedure.
Energies 11 01504 g001
Figure 2. General methodology of location selection problem.
Figure 2. General methodology of location selection problem.
Energies 11 01504 g002
Figure 3. Research methodologies.
Figure 3. Research methodologies.
Energies 11 01504 g003
Figure 4. Triangular Fuzzy Number.
Figure 4. Triangular Fuzzy Number.
Energies 11 01504 g004
Figure 5. Viet Nam’s solar resources map. (Source: The World Bank).
Figure 5. Viet Nam’s solar resources map. (Source: The World Bank).
Energies 11 01504 g005
Figure 6. Inputs and Output of DEA Models.
Figure 6. Inputs and Output of DEA Models.
Energies 11 01504 g006
Figure 7. The Hierarchical structures for the FAHP approach.
Figure 7. The Hierarchical structures for the FAHP approach.
Energies 11 01504 g007
Figure 8. Geometric distance from PIS and NIS.
Figure 8. Geometric distance from PIS and NIS.
Energies 11 01504 g008
Figure 9. The result from TOPSIS model.
Figure 9. The result from TOPSIS model.
Energies 11 01504 g009
Table 1. 1–9 Saaty Scale.
Table 1. 1–9 Saaty Scale.
Importance IntensityDefinition
1Equally importance
3Moderate importance
5Strongly more importance
7Very strong more importance
9Extremely importance
2,4,6,8Intermediate values
Table 2. The fuzzy conversion scale.
Table 2. The fuzzy conversion scale.
Importance IntensityTriangular Fuzzy ScaleImportance IntensityTriangular Fuzzy Scale
1(1, 1, 1)1/1(1, 1, 1)
2(1, 2, 3)1/2(1/3, 1/2, 1/1)
3(2, 3, 4)1/3(1/4, 1/3, 1/2)
4(3, 4, 5)1/4(1/5, 1/4, 1/3)
5(4, 5, 6)1/5(1/6, 1/5, 1/4)
6(5, 6, 7)1/6(1/7, 1/6, 1/5)
7(6, 7, 8)1/7(1/8, 1/7, 1/6)
8(7, 8, 9)1/8(1/9, 1/8, 1/7)
9(9, 9, 9)1/9(1/9, 1/9, 1/9)
Table 3. List of the 46 locations identified in Viet Nam.
Table 3. List of the 46 locations identified in Viet Nam.
No.SitesDMUsNo.SitesDMUs
1Bac GiangDMU 124My ThoDMU 24
2Binh ThuanDMU 225Nam DinhDMU 25
3Buon Ma ThuocDMU 326Nha TrangDMU 26
4Ca MauDMU 427Ninh BinhDMU 27
5Cam RanhDMU 528Phu LienDMU 28
6Can ThoDMU 629Phu QuocDMU 29
7Cang LongDMU 730Phuoc LongDMU 30
8Chau DocDMU 831PleikuDMU 31
9Con SonDMU 932Quang NgaiDMU 32
10Da NangDMU 1033Quy NhonDMU 33
11Dien BienDMU 1134Rach GiaDMU 34
12Dong HaDMU 1235Soc TrangDMU 35
13Dong HoiDMU 1336Son LaDMU 36
14Ha TinhDMU 1437Tan Son NhatDMU 37
15Hai DuongDMU 1538Tay NinhDMU 38
16Hoa BinhDMU 1639Thai BinhDMU 39
17Hoang SaDMU 1740Thanh HoaDMU 40
18Hong GaiDMU 1841Tuy HoaDMU 41
19HueDMU 1942Uong BiDMU 42
20Hung YenDMU 2043Viet TriDMU 43
21Kon TumDMU 2144VinhDMU 44
22Lai ChauDMU 2245Vinh YenDMU 45
23Moc HoaDMU 2346Vung TauDMU 46
Table 4. Data set of the 46 DMUs.
Table 4. Data set of the 46 DMUs.
DMUsTemperature Wind Speed Sunshine Hours Elevation
DMU 123.401.801695.0029.00
DMU 226.802.202878.0010.00
DMU 323.602.802460.00467.00
DMU 426.801.302300.001.00
DMU 526.902.802672.0020.00
DMU 626.601.502561.001.00
DMU 726.801.802621.001.00
DMU 827.201.702589.002.00
DMU 927.002.602351.00120.00
DMU 1025.801.502182.005.00
DMU 1122.000.902034.00490.00
DMU 1225.102.601910.0010.00
DMU 1324.502.501857.0013.00
DMU 1423.901.501664.007.00
DMU 1523.402.401658.001.00
DMU 1623.401.001641.0023.00
DMU 1726.804.802788.0038.00
DMU 1823.102.701690.003.00
DMU 1925.201.501970.003.00
DMU 2023.301.701625.004.00
DMU 2123.501.402374.00530.00
DMU 2223.000.801824.00213.21
DMU 2327.302.002686.001.00
DMU 2427.001.702645.001.00
DMU 2523.502.201619.001.00
DMU 2626.602.402540.003.00
DMU 2723.601.901611.0012.00
DMU 2823.103.001693.008.00
DMU 2927.103.002364.0053.00
DMU 3025.501.602521.00192.00
DMU 3121.702.702412.00756.00
DMU 3225.701.302248.0014.00
DMU 3326.901.902470.008.00
DMU 3427.402.802470.001.00
DMU 3526.801.702423.001.00
DMU 3621.101.102000.00673.00
DMU 3727.402.802489.004.00
DMU 3826.901.502672.0020.00
DMU 3923.302.101639.003.00
DMU 4023.601.701690.0018.00
DMU 4126.502.202467.002.00
DMU 4223.501.901920.0037.00
DMU 4323.501.501601.0020.00
DMU 4423.901.801677.0010.00
DMU 4523.801.601670.0018.00
DMU 4626.703.002728.001.00
Table 5. The results of the DEA models.
Table 5. The results of the DEA models.
DMUsDEA MODEL
CCR-ICCR-OBCC-IBCC-OSBM-I-C
DMU 10.683830.683830.901710.696050.60007
DMU 211111
DMU 30.942150.942150.942430.954050.82069
DMU 40.894230.894230.894270.935210.85555
DMU 50.908420.908420.912470.930740.76786
DMU 60.966550.966550.966630.9670.93596
DMU 70.947960.947960.949290.951090.89048
DMU 80.934560.934560.936930.948050.87444
DMU 90.801410.801410.820140.841110.67095
DMU 100.842490.842490.860320.843950.77921
DMU 1111111
DMU 120.696210.696210.840640.701520.57406
DMU 130.694340.694340.861220.696070.57425
DMU 140.68310.68310.882850.683190.62196
DMU 150.64880.64880.901710.651910.53603
DMU 160.751640.751640.933690.76670.74231
DMU 170.935920.935920.963250.974860.69601
DMU 180.662190.662190.913420.665370.53413
DMU 190.774020.774020.83730.776120.71333
DMU 200.665440.665440.905580.675320.58865
DMU 2111111
DMU 2211111
DMU 230.9380.9380.938290.952770.86668
DMU 240.960610.960610.963690.968560.91935
DMU 250.635470.635470.897870.642930.53539
DMU 260.883240.883240.894320.88820.76302
DMU 270.638320.638320.894070.650910.55676
DMU 280.659370.659370.913420.666560.52121
DMU 290.793510.793510.804160.83110.64609
DMU 300.970020.970020.971370.972450.93655
DMU 3111111
DMU 320.902380.902380.902660.914070.85562
DMU 330.881640.881640.888620.88540.78419
DMU 340.82570.82570.833480.858230.68268
DMU 350.88540.88540.886980.88910.80428
DMU 3611111
DMU 370.832050.832050.838920.864840.68793
DMU 3811111
DMU 390.650710.650710.905580.661330.55309
DMU 400.684970.684970.894070.692920.60721
DMU 410.866230.866230.884880.867050.75471
DMU 420.763660.763660.897870.779190.66542
DMU 430.666450.666450.897870.669310.60457
DMU 440.665270.665270.882850.673580.58559
DMU 450.680020.680020.886550.683650.61039
DMU 460.928170.928170.942260.95090.78322
Table 6. Fuzzy comparison matrix for criteria.
Table 6. Fuzzy comparison matrix for criteria.
CriteriaECENSCSOTE
EC(1, 1, 1)(1/4, 1/3, 1/2)(1, 2, 3)(1, 2, 3)(1/3, 1/2, 1)
EN(2, 3, 4)(1, 1, 1)(2, 3, 4)(1, 2, 3)(1, 1, 1)
SC(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1, 1, 1)(2, 3, 4)(1, 1, 1)
SO(1/3, 1/2, 1)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1, 1, 1)(1/3, 1/2, 1)
TE(1, 2, 3)(1, 1, 1)(1, 1, 1)(1, 2, 3)(1, 1, 1)
Table 7. Real number priority.
Table 7. Real number priority.
CriteriaECENSCSOTE
EC11/3221/2
EN31321
SC1/21/3131
SO1/21/21/311/2
TE21121
Table 8. Fuzzy comparison matrices for criteria.
Table 8. Fuzzy comparison matrices for criteria.
CriteriaECENSCSOTEWeight
EC(1, 1, 1)(1/4, 1/3, 1/2)(1, 2, 3)(1, 2, 3)(1/3, 1/2, 1)0.17201
EN(2, 3, 4)(1, 1, 1)(2, 3, 4)(1, 2, 3)(1, 1, 1)0.32965
SC(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1, 1, 1)(2, 3, 4)(1, 1, 1)0.16526
SO(1/3, 1/2, 1)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1, 1, 1)(1/3, 1/2, 1)0.09694
TE(1, 2, 3)(1, 1, 1)(1, 1, 1)(1, 2, 3)(1, 1, 1)0.23614
Total1
CR = 0.09598
Table 9. Comparison matrix for environmental criteria.
Table 9. Comparison matrix for environmental criteria.
CriteriaEN 3EN 2EN 1Weight
EN3(1, 1, 1)(1/3, 1/2, 1)(1/3, 1/2, 1)0.19580
EN2(1, 2, 3)(1, 1, 1)(1, 2, 3)0.49339
EN1(1, 2, 3)(1/3, 1/2, 1)(1, 1, 1)0.31081
Total1
CR = 0.05156
Table 10. Comparison matrix for site characteristics criteria.
Table 10. Comparison matrix for site characteristics criteria.
CriteriaSC 3SC 2SC 1Weight
SC 3(1, 1, 1)(1, 2, 3)(2, 3, 4)0.52784
SC 2(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)0.33251
SC 1(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1, 1, 1)0.13965
Total1
CR = 0.05156
Table 11. Comparison matrix for social criteria.
Table 11. Comparison matrix for social criteria.
CriteriaSO 3SO 2SO 1Weight
SO 3(1, 1, 1)(1/3, 1/2, 1)(1/4, 1/3, 1/2)0.15706
SO 2(1, 2, 3)(1, 1, 1)(1/4, 1/3, 1/2)0.24931
SO 1(2, 3, 4)(2, 3, 4)(1, 1, 1)0.59363
Total1
CR = 0.05156
Table 12. Comparison matrix for economic criteria.
Table 12. Comparison matrix for economic criteria.
CriteriaEC 1EC 3EC 2Weight
EC 1(1, 1, 1)(1/3, 1/2, 1)(1, 2, 3)0.31081
EC 3(1, 2, 3)(1, 1, 1)(1, 2, 3)0.49339
EC 2(1/3, 1/2, 1)(1/3, 1/2, 1)(1, 1, 1)0.19580
Total1
CR = 0.05156
Table 13. Comparison matrix for technological criteria.
Table 13. Comparison matrix for technological criteria.
CriteriaTE 1TE 2TE 3Weight
TE 1(1, 1, 1)(1/3, 1/2, 1)(1/3, 1/2, 1)0.19580
TE 2(1, 2, 3)(1, 1, 1)(1/3, 1/2, 1)0.31081
TE 3(1, 2, 3)(1, 2, 3)(1, 1, 1)0.49339
Total1
CR = 0.05156
Table 14. Comparison matrix for TE1 based on the alternatives.
Table 14. Comparison matrix for TE1 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(2, 3, 4)(4, 5, 6)(2, 3, 4)(6, 7, 8)(2, 3, 4)(1, 2, 3)0.354451
DMU 11(1/4, 1/3, 1/2)(1, 1, 1)(1, 2, 3)(4, 5, 6)(2, 3, 4)(1, 2, 3)(3, 4, 5)0.208178
DMU 21(1/6, 1/5, 1/4)(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)(1, 2, 3)(1, 1, 1)(2, 3, 4)0.125652
DMU 22(1/4, 1/3, 1/2)(1/6, 1/5, 1/4)(1/4, 1/3, 1/2)(1, 1, 1)(1/5, 1/4, 1/3)(1/3, 1/2, 1)(1/4, 1/3, 1/2)0.045139
DMU 31(1/8, 1/7, 1/6)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(3, 4, 5)(1, 1, 1)(1/3, 1/2, 1)(1/3, 1/2, 1)0.069171
DMU 36(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)(1, 2, 3)(1, 2, 3)(1, 1, 1)(1, 1, 1)0.100762
DMU 38(1/3, 1/2, 1)(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(2, 3, 4)(1, 2, 3)(1, 1, 1)(1, 1, 1)0.096648
Total1
CR = 0.08404
Table 15. Comparison matrix for TE2 based on the alternatives.
Table 15. Comparison matrix for TE2 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(1, 2, 3)(3, 4, 5)(2, 3, 4)(5, 6, 7)(4, 5, 6)(2, 3, 4)0.333294
DMU 11(1/3, 1/2, 1)(1, 1, 1)(3, 4, 5)(1, 2, 3)(2, 3, 4)(2, 3, 4)(5, 6, 7)0.233267
DMU 21(1/5, 1/4, 1/3)(1/5, 1/4, 1/3)(1, 1, 1)(1, 2, 3)(1, 1, 1)(1, 2, 3)(2, 3, 4)0.109613
DMU 22(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/3, 1/2, 1)(1, 1, 1)(3, 4, 5)(2, 3, 4)(1, 2, 3)0.13301
DMU 31(1/7, 1/6, 1/5)(1/4, 1/3, 1/2)(1, 1, 1)(1/5, 1/4, 1/3)(1, 1, 1)(3, 4, 5)(2, 3, 4)0.092202
DMU 36(1/6, 1/5, 1/4)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)(1, 1, 1)(1, 2, 3)0.052846
DMU 38(1/4, 1/3, 1/2)(1/7, 1/6, 1/5)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)0.045768
Total1
CR = 0.09594
Table 16. Comparison matrix for TE3 based on the alternatives.
Table 16. Comparison matrix for TE3 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(1, 2, 3)(3, 4, 5)(1/3, 1/2, 1)(5, 6, 7)(3, 4, 5)(1, 2, 3)0.293807
DMU 11(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)(1/5, 1/4, 1/3)(1, 2, 3)(2, 3, 4)(1, 2, 3)0.214438
DMU 21(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1, 1, 1)(1/3, 1/2, 1)(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1/3, 1/2, 1)0.049645
DMU 22(1/3, 1/2, 1)(1/5, 1/4, 1/3)(1, 2, 3)(1, 1, 1)(1, 2, 3)(1, 1, 1)(1/5, 1/4, 1/3)0.086324
DMU 31(1/7, 1/6, 1/5)(1/3, 1/2, 1)(3, 4, 5)(1/3, 1/2, 1)(1, 1, 1)(1/3, 1/2, 1)(1/6, 1/5, 1/4)0.071619
DMU 36(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(2, 3, 4)(1, 1, 1)(1, 2, 3)(1, 1, 1)(1/4, 1/3, 1/2)0.087677
DMU 38(1/3, 1/2, 1)(1/3, 1/2, 1)(1, 2, 3)(3, 4, 5)(4, 5, 6)(2, 3, 4)(1, 1, 1)0.196491
Total1
CR = 0.08852
Table 17. Comparison matrix for SC1 based on the alternatives.
Table 17. Comparison matrix for SC1 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(1/3, 1/2, 1)(2, 3, 4)(1, 1, 1)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(2, 3, 4)0.13004
DMU 11(1, 2, 3)(1, 1, 1)(1, 2, 3)(2, 3, 4)(1/3, 1/2, 1)(1, 2, 3)(2, 3, 4)0.197223
DMU 21(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)(1, 2, 3)(1/6, 1/5, 1/4)(1, 2, 3)(1, 2, 3)0.113704
DMU 22(1, 1, 1)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)(1/3, 1/2, 1)(1/3, 1/2, 1)(1, 1, 1)0.077497
DMU 31(1, 2, 3)(1, 2, 3)(4, 5, 6)(1, 2, 3)(1, 1, 1)(1, 2, 3)(2, 3, 4)0.274544
DMU 36(2, 3, 4)(1/3, 1/2, 1)(1/3, 1/2, 1)(1, 2, 3)(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)0.149847
DMU 38(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1, 1, 1)0.057146
Total1
CR = 0.09079
Table 18. Comparison matrix for SC2 based on the alternatives.
Table 18. Comparison matrix for SC2 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(2, 3, 4)(4, 5, 6)(2, 3, 4)(3, 4, 5)(1, 2, 3)(2, 3, 4)0.3033
DMU 11(1/4, 1/3, 1/2)(1, 1, 1)(4, 5, 6)(3, 4, 5)(1, 2, 3)(1/3, 1/2, 1)(2, 3, 4)0.189088
DMU 21(1/6, 1/5, 1/4)(1/6, 1/5, 1/4)(1, 1, 1)(1/5, 1/4, 1/3)(1, 1, 1)(1/4, 1/3, 1/2)(1/3, 1/2, 1)0.044302
DMU 22(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)(3, 4, 5)(1, 1, 1)(3, 4, 5)(1/3, 1/2, 1)(2, 3, 4)0.132378
DMU 31(1/5, 1/4, 1/3)(1/3, 1/2, 1)(1, 1, 1)(1/5, 1/4, 1/3)(1, 1, 1)(1/5, 1/4, 1/3)(2, 3, 4)0.071387
DMU 36(1/3, 1/2, 1)(1, 2, 3)(2, 3, 4)(1, 2, 3)(3, 4, 5)(1, 1, 1)(1, 2, 3)0.197566
DMU 38(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1, 2, 3)(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)0.06198
Total1
CR = 0.09473
Table 19. Comparison matrix for SC3 based on the alternatives.
Table 19. Comparison matrix for SC3 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(3, 4, 5)(2, 3, 4)(1, 2, 3)(3, 4, 5)(6, 7, 8)(1, 2, 3)0.314672
DMU 11(1/5, 1/4, 1/3)(1, 1, 1)(1, 1, 1)(2, 3, 4)(3, 4,5 )(2, 3, 4)(1/3, 1/2, 1)0.144856
DMU 21(1/4, 1/3, 1/2)(1, 1, 1)(1, 1, 1)(1, 2, 3)(3, 4, 5)(4, 5, 6)(1, 2, 3)0.176472
DMU 22(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)(5, 6, 7)(1/4, 1/3, 1/2)0.107392
DMU 31(1/5, 1/4, 1/3)(1/5, 1/4, 1/3)(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1, 1, 1)(2, 3, 4)(1/5, 1/4, 1/3)0.051436
DMU 36(1/8, 1/7, 1/6)(1/4, 1/3, 1/2)(1/6, 1/5, 1/4)(1/7, 1/6, 1/5)(1/4, 1/3, 1/2)(1, 1, 1)(1/3, 1/2, 1)0.036723
DMU 38(1/3, 1/2, 1)(1, 2, 3)(1/3, 1/2, 1)(2, 3, 4)(3, 4, 5)(1, 2, 3)(1, 1, 1)0.168449
Total1
CR = 0.09232
Table 20. Comparison matrix for SO1 based on the alternatives.
Table 20. Comparison matrix for SO1 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(4, 5, 6)(3, 4, 5)(1, 2, 3)(2, 3, 4)(4, 5, 6)(3, 4, 5)0.369784
DMU 11(1/6, 1/5, 1/4)(1, 1, 1)(2, 3, 4)(1, 2, 3)(3, 4, 5)(2, 3, 4)(1, 2, 3)0.183518
DMU 21(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1, 1, 1)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/3, 1/2, 1)(1/5, 1/4, 1/3)0.044326
DMU 22(1/3, 1/2, 1)(1/3, 1/2, 1)(2, 3, 4)(1, 1, 1)(1/3, 1/2, 1)(1, 2, 3)(1/4, 1/3, 1/2)0.095751
DMU 31(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)(1, 2, 3)(1, 2, 3)(1, 1, 1)(1/3, 1/2, 1)(1/4, 1/3, 1/2)0.079245
DMU 36(1/6, 1/5, 1/4)(1/4, 1/3, 1/2)(1, 2, 3)(1/3, 1/2, 1)(1, 2, 3)(1, 1, 1)(1/3, 1/2, 1)0.07467
DMU 38(1/5, 1/4, 1/3)(1/3, 1/2, 1)(3, 4, 5)(2, 3, 4)(1, 2, 3)(1, 2, 3)(1, 1, 1)0.152707
Total1
CR = 0.09669
Table 21. Comparison matrix for SO2 based on the alternatives.
Table 21. Comparison matrix for SO2 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(1, 2, 3)(5, 6, 7)(1, 2, 3)(3, 4, 5)(4, 5, 6)(2, 3, 4)0.305459
DMU 11(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)(1, 2, 3)(3, 4, 5)(5, 6, 7)(2, 3, 4)0.243481
DMU 21(1/7, 1/6, 1/5)(1/4, 1/3, 1/2)(1, 1, 1)(1/6, 1/5, 1/4)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/4, 1/3, 1/2)0.03918
DMU 22(1/3, 1/2, 1)(1/3, 1/2, 1)(4, 5, 6)(1, 1, 1)(2, 3, 4)(4, 5, 6)(2, 3, 4)0.191699
DMU 31(1/5, 1/4, 1/3)(1/5, 1/4, 1/3)(2, 3, 4)(1/4, 1/3, 1/2)(1, 1, 1)(1, 2, 3)(2, 3, 4)0.098135
DMU 36(1/6, 1/5, 1/4)(1/7, 1/6, 1/5)(1, 2, 3)(1/6, 1/5, 1/4)(1/3, 1/2, 1)(1, 1, 1)(1, 1, 1)0.052314
DMU 38(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(2, 3, 4)(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1, 1, 1)(1, 1, 1)0.069733
Total1
CR = 0.05496
Table 22. Comparison matrix for SO3 based on the alternatives.
Table 22. Comparison matrix for SO3 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(1, 2, 3)(2, 3, 4)(4, 5, 6)(3, 4, 5)(1, 2, 3)(2, 3, 4)0.307084
DMU 11(1/3, 1/2, 1)(1, 1, 1)(1, 1, 1)(4, 5, 6)(3, 4, 5)(2, 3, 4)(1, 2, 3)0.218813
DMU 21(1/4, 1/3, 1/2)(1, 1, 1)(1, 1, 1)(1, 2, 3)(2, 3, 4)(1, 2, 3)(3, 4, 5)0.165351
DMU 22(1/6, 1/5, 1/4)(1/6, 1/5, 1/4)(1/3, 1/2, 1)(1, 1, 1)(2, 3, 4)(1, 2, 3)(3, 4, 5)0.112322
DMU 31(1/5, 1/4, 1/3)(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1, 1, 1)(1, 2, 3)(1, 1, 1)0.064014
DMU 36(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/3, 1/2, 1)(1/3, 1/2, 1)(1, 1, 1)(1, 2, 3)0.075672
DMU 38(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1/5, 1/4, 1/3)(1/5, 1/4, 1/3)(1, 1, 1)(1/3, 1/2, 1)(1, 1, 1)0.056743
Total1
CR = 0.09264
Table 23. Comparison matrix for EN1 based on the alternatives.
Table 23. Comparison matrix for EN1 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(2, 3, 4)(4, 5, 6)(6, 7, 8)(4, 5, 6)(3, 4, 5)(2, 3, 4)0.379402
DMU 11(1/4, 1/3, 1/2)(1, 1, 1)(4, 5, 6)(1, 2, 3)(1, 2, 3)(2, 3, 4)(3, 4, 5)0.216035
DMU 21(1/6, 1/5, 1/4)(1/6, 1/5, 1/4)(1, 1, 1)(2, 3, 4)(1, 1, 1)(2, 3, 4)(1, 2, 3)0.107386
DMU 22(1/8, 1/7, 1/6)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1, 1, 1)(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1/3, 1/2, 1)0.042011
DMU 31(1/6, 1/5, 1/4)(1/3, 1/2, 1)(1, 1, 1)(3, 4, 5)(1, 1, 1)(1, 1, 1)(3, 4, 5)0.114659
DMU 36(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(2, 3, 4)(1, 1, 1)(1, 1, 1)(1, 2, 3)0.082888
DMU 38(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)(1/3, 1/2, 1)(1, 2, 3)(1/5, 1/4, 1/3)(1/3, 1/2, 1)(1, 1, 1)0.057619
Total1
CR = 0.09124
Table 24. Comparison matrix for EN2 based on the alternatives.
Table 24. Comparison matrix for EN2 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(2, 3, 4)(3, 4, 5)(1, 2, 3)(4, 5, 6)(2, 3, 4)(3, 4, 5)0.330479
DMU 11(1/4, 1/3, 1/2)(1, 1, 1)(5, 6, 7)(3, 4, 5)(3, 4, 5)(1, 2, 3)(2, 3, 4)0.233849
DMU 21(1/5, 1/4, 1/3)(1/8, 1/7, 1/6)(1, 1, 1)(1, 1, 1)(2, 3, 4)(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)0.058942
DMU 22(1/3, 1/2, 1)(1/5, 1/4, 1/3)(1, 1, 1)(1, 1, 1)(1, 2, 3)(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)0.068136
DMU 31(1/6, 1/5, 1/4)(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(1/3, 1/2, 1)(1, 1, 1)(1/3, 1/2, 1)(1/5, 1/4, 1/3)0.041927
DMU 36(1/4, 1/3, 1/2)(1/3, 1/2, 1)(2, 3, 4)(2, 3, 4)(1, 2, 3)(1, 1, 1)(1, 2, 3)0.14113
DMU 38(1/5, 1/4, 1/3)(1/4, 1/3, 1/2)(3, 4, 5)(2, 3, 4)(3, 4, 5)(1/3, 1/2, 1)(1, 1, 1)0.125536
Total1
CR = 0.08438
Table 25. Comparison matrix for EN3 based on the alternatives.
Table 25. Comparison matrix for EN3 based on the alternatives.
DMUsDMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38Weight
DMU 2(1, 1, 1)(4, 5, 6)(2, 3, 4)(3, 4, 5)(1, 2, 3)(5, 6, 7)(4, 5, 6)0.3751
DMU 11(1/6, 1/5, 1/4)(1, 1, 1)(2, 3, 4)(1, 1, 1)(2, 3, 4)(1, 2, 3)(1, 2, 3)0.15814
DMU 21(1/4, 1/3, 1/2)(1/4, 1/3, 1/2)(1, 1, 1)(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)(1, 2, 3)(1/4, 1/3, 1/2)0.055333
DMU 22(1/5, 1/4, 1/3)(1, 1, 1)(2, 3, 4)(1, 1, 1)(1, 2, 3)(2, 3, 4)(1, 1, 1)0.135332
DMU 31(1/3, 1/2, 1)(1/4, 1/3, 1/2)(3, 4, 5)(1/3, 1/2, 1)(1, 1, 1)(3, 4, 5)(1, 2, 3)0.133423
DMU 36(1/7, 1/6, 1/5)(1/3, 1/2, 1)(1/3, 1/2, 1)(1/4, 1/3, 1/2)(1/5, 1/4, 1/3)(1, 1, 1)(1/5, 1/4, 1/3)0.040685
DMU 38(1/6, 1/5, 1/4)(1/3, 1/2, 1)(2, 3, 4)(1, 1, 1)(1/3, 1/2, 1)(3, 4, 5)(1, 1, 1)0.101987
Total1
CR = 0.08831
Table 26. The weight of criteria.
Table 26. The weight of criteria.
SymbolCriteriaWeight
SO 1Support mechanisms (SO 1)0.05755
SO 2Protection law (SO 2)0.02417
SO 3Legal and Regulatory compliance (SO 3)0.01522
EN 1Temperature (EN 1)0.10246
EN 2Sunshine hours (EN 2)0.16264
EN 3Humidity (EN 3)0.06454
TE 1Distance from major road (TE 1)0.04624
TE 2Distance from power network (TE 2)0.07339
TE 3Potential demand (TE 3)0.1165
EC 1Constructions cost (EC 1)0.05346
EC 2Operations and Maintenances Cost (EC 2)0.03368
EC 3New feeder cost (EC 3)0.08487
SC 1Ecology (SC 1)0.05495
SC 2Elevation (SC 2)0.02308
SC 3Approachability (SC 3)0.08723
Table 27. The weights of alternative locations with respect to sub criteria.
Table 27. The weights of alternative locations with respect to sub criteria.
Sub-CriteriaDMUs
DMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38
EN 10.319780.212270.050780.124540.178700.069680.04426
EN 20.330480.233850.058940.068140.041930.141130.12554
EN 30.375100.158140.055330.135330.133420.040690.10199
SC 10.130040.197220.113700.077500.274540.149850.05715
SC 20.303100.189090.0443010.132380.071390.197570.06198
SC 30.314670.144860.176470.107390.051430.036720.16845
SO 10.369780.183520.044320.095750.079240.074670.15271
SO 20.305460.2434820.039180.191700.098130.052310.06973
SO 30.307080.2188120.165350.112320.064010.075670.05674
EC 10.319780.212270.050780.124540.178700.069680.04426
EC 20.328010.102820.046740.064510.194370.160460.10306
EC 30.376790.215220.065110.127520.035460.058190.12170
TE 10.354450.208180.125650.045140.069170.100760.09665
TE 20.333290.233270.109610.133010.092200.052850.04577
TE 30.297820.212700.049600.077310.077130.088540.19689
Table 28. Normalized Decision Matrix.
Table 28. Normalized Decision Matrix.
Sub-CriteriaDMUs
DMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38
EN 10.319780.212270.050780.124540.178700.069680.04426
EN 20.330480.233850.058940.068140.041930.141130.12554
EN 30.3750100.158140.055330.135330.133420.040690.10199
SC 10.130040.197220.113700.077100.274540.149850.05715
SC 20.303100.189090.044300.132380.071390.197570.06198
SC 30.314670.144860.176470.107390.051440.036720.16845
SO 10.369780.183520.044320.095750.079240.074670.15271
SO 20.305460.243480.039180.191700.098130.052310.06973
SO 30.307080.218810.165350.112320.064010.075670.05674
EC 10.319780.212270.050780.124540.178700.069660.04426
EC 20.328010.102820.046740.064510.194370.160480.10306
EC 30.376790.215220.065110.127520.035450.058190.12170
TE 10.354450.208180.125650.045140.069170.100760.09665
TE 20.333290.233270.109610.133010.092200.052850.04577
TE 30.297830.212700.049100.077310.077130.088540.19689
Table 29. The weighted Normalized Decision Matrix.
Table 29. The weighted Normalized Decision Matrix.
DMUsMain-Criteria
EN (0.0992)SC (0.0536)SO (0.0916)EC (0.0929)TE (0.0907)
DMU 20.319990.235240.370840.125910.26166
DMU 110.076510.115560.068480.215900.07157
DMU 210.028160.034640.100100.065390.04642
DMU 220.091990.143940.090650.028590.07826
DMU 310.039710.040450.056070.033320.04939
DMU 360.140220.371210.047450.097610.10414
DMU 380.303430.058960.266410.433270.38857
Table 30. Normalized Weight Matrix.
Table 30. Normalized Weight Matrix.
CriteriaDMUs
DMU 2DMU 11DMU 21DMU 22DMU 31DMU 36DMU 38
SO 10.022100.015800.015800.015800.022100.022100.02840
SO 20.010700.008400.006000.006000.008400.009500.01070
SO 30.040400.030300.035300.025200.025200.040400.03030
EN 10.045500.035400.040400.040400.030300.030300.04550
EN 20.079400.052900.061700.044100.061700.052900.07060
EN 30.013500.033600.013500.026900.033600.020200.02020
TE 10.013500.013500.020200.006700.020200.027000.01350
TE 20.026100.026100.034800.017400.017400.043500.01740
TE 30.038700.045100.051500.038700.058000.032200.03870
EC 10.015400.023100.023100.023100.015400.015400.02310
EC 20.009300.009300.013900.013900.013900.013900.01390
EC 30.024000.036000.024000.036000.048000.024000.02400
SC 10.024400.021700.024400.019000.021700.016300.01630
SC 20.011900.009500.009500.002400.004800.009500.00950
SC 30.033900.043600.024200.043600.019400.029100.02910

Share and Cite

MDPI and ACS Style

Wang, C.-N.; Nguyen, V.T.; Thai, H.T.N.; Duong, D.H. Multi-Criteria Decision Making (MCDM) Approaches for Solar Power Plant Location Selection in Viet Nam. Energies 2018, 11, 1504. https://doi.org/10.3390/en11061504

AMA Style

Wang C-N, Nguyen VT, Thai HTN, Duong DH. Multi-Criteria Decision Making (MCDM) Approaches for Solar Power Plant Location Selection in Viet Nam. Energies. 2018; 11(6):1504. https://doi.org/10.3390/en11061504

Chicago/Turabian Style

Wang, Chia-Nan, Van Thanh Nguyen, Hoang Tuyet Nhi Thai, and Duy Hung Duong. 2018. "Multi-Criteria Decision Making (MCDM) Approaches for Solar Power Plant Location Selection in Viet Nam" Energies 11, no. 6: 1504. https://doi.org/10.3390/en11061504

APA Style

Wang, C. -N., Nguyen, V. T., Thai, H. T. N., & Duong, D. H. (2018). Multi-Criteria Decision Making (MCDM) Approaches for Solar Power Plant Location Selection in Viet Nam. Energies, 11(6), 1504. https://doi.org/10.3390/en11061504

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop