Next Article in Journal
A Study on an Out-of-Step Detection Algorithm Using the Time Variation of Complex Power: Part I, Mathematical Modeling
Previous Article in Journal
Lessons to Learn from Post-Installation Pollution Levels Assessment of Some Distribution Insulators
Previous Article in Special Issue
Vector Speed Regulation of an Asynchronous Motor Based on Improved First-Order Linear Active Disturbance Rejection Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Site Selection for a Solar Power Plant in the Mekong Delta Region of Vietnam

1
Department of Industrial Engineering and Management, National Kaohsiung University of Science and Technology, Kaohsiung 80778, Taiwan
2
Department of Logistics and Supply Chain Management, Hong Bang International University, Ho Chi Minh 723000, Vietnam
*
Authors to whom correspondence should be addressed.
Energies 2020, 13(16), 4066; https://doi.org/10.3390/en13164066
Submission received: 21 May 2020 / Revised: 23 July 2020 / Accepted: 30 July 2020 / Published: 6 August 2020

Abstract

:
Following the recent development trend in the struggle for cleaning the earth’s environment, solar is the one of most promising area that can partially be used as a replaceable energy from non-renewable fuel sources. As such, it plays a significant role in protecting the environment from global warming. As solar power does not emit harmful gases into the atmosphere, its production, distribution, setup, and operation are vital should the production remain constant. Even solar energy waste emissions are small; when compared to current energy sources, the amount of harmful gases is negligible. This paper presented an integrated approach for site of solar plants by using data envelopment analysis (DEA) and Fuzzy Analytical Network Process (FANP). Furthermore, these integrated methodologies, incorporated with the most relevant parameters of requirements for solar plants, are introduced. First, the paper considers an integrated hierarchical DEA and FANP model for the optimal geographical location of solar plants in Mekong Delta Region, Vietnam. Using the proposed model for implementation would allow the renewable energy policy makers to select and control the optimal location for allocating and constructing a solar energy power plant in Vietnam. This is the preferred strategy for location optimization problems associated with solar plant units in Vietnam and around the world.

1. Introduction

Energy has been a fundamental driver of global economies in order maintain development in almost any industry. It is compulsory that the sources of energy must be maintained and continue developing to meet the demand. With fossil fuel continuing to be proven to be detrimental to the environment, solar energy is continuing to become a modern, growing, and potential source of renewable energy. According to the International Energy Agency (IEA) in 2015, the consumption of fossil fuels coming from coal, oil, and natural gas accounts for 81.6% of the total energy consumption amongst the 38 leading countries. From the same data, wind, and solar energy only accounts for 1.5% of the energy consumption [1]. As also shown in Figure 1, solar energy has the highest development compared to other renewable energies such as hydropower, onshore and offshore wind energy, and bioenergy [2]. According to NASA’s physics and astronomy calculations, solar energy is still available for about 6.5 billion years before the sun fades away, which means that the energy source can be utilized for a longer and more durable time before it is projected to deplete [3].
Solar energy takes the direct sun rays from the sun and transforms them into usable energy consumed for many purposes. The major limitations of solar energy include geographical conditions, high dependency on technology development, large installation area with a high initial cost. Additionally, with such limitations come additional, big environmental changes where modifications could leave an impact on local ecosystems should the development a plant be initialized with the trend of sustainable development, and maintaining the balance of the triple bottom line of environmental, social, and economical elements must be considered [5].
Solar energy is still not a popular energy source amongst the economic and social development of the provinces in the Mekong Delta region, Vietnam. However, according to Vietnam Academy of Science and Technology (VAST), Vietnam is one of the potential countries for solar energy development when looking at the data for the amount of sun rays the world receives annually, shown in Figure 2 with a second highest overall evaluation of 4–5 Kwh/m2/day [6].
In fact, if more than 50% [7] of the energy supply for a national electricity system is imported from other countries, the risk of energy unavailability will be extremely high. Once that supply is depleted or affected by natural disasters, accidents, armed conflicts, sabotage, policy changes or strikes within coal mines, energy security will be directly affected. Besides, the coal transportation from a remote area to the electric production factory, which is done mainly by sea, also has many potential implications for the environment, social security, and economy. The pollution of dust, noise, and ash accumulation from the operation of coal-fired power plants can cause frustration in the local community and sow the seeds for protests, creating long-term instability in terms of social security [8]. Located in the tropical and monsoon region, with an abundant solar energy, the Mekong Delta of Vietnam is entering a period of strong solar power development. Solar energy has been selected as one of the major alternative renewable energy sources in the energy development strategy of the Mekong Delta region. Solar Radiation Maps of Vietnam is shown in Figure 3.
Currently, there is very little research for developing a solar power plant in Vietnamese Mekong Delta region. Determining an optimal location, which is sustainable and satisfies the triple bottom line concepts, is essential for the growth of the country’s power supply and alleviates the environmental damage from fossil fuel consumption by promoting a renewable alternative energy supply. Therefore, this problem is tackled as a Multi-Criteria Decision Making (MCDM) problem that requires both quantitative and qualitative information to be considered. The MCDM problem is commonly defined as determining an optimal decision based on a finite set of available decision alternatives with multiple, and potentially conflicting, criteria.
As mentioned, the purpose of this paper is now to propose a suitable MCDM methodology for localizing and selecting an optimal area in Vietnam in developing a solar power plant. In the next section, the literature will be presented to discuss the first stage of the paper. Discussions and results’ analysis are explained towards the end of the paper.

2. Literature Review

Nowadays, there are many researchers applied MCDM models to any fields of sciences. A. Kengpol et al. [10] initially conducted a study on proposing a solar power plant location to avoid flooding in Thailand using geographical information system (GIS) to determine an optimum site condition. Upon the study, a number of geographical conditions had to be considered which developed into an MCDM problem of ranking the priorities amongst the conditions using Fuzzy Analytic Hierarchy Process (FAHP) and the Technique of Order Preference by Similarity to Ideal Solution (TOPSIS).
A. Azadeh et al. [11] integrated a hierarchical model for solar plants’ location selection by data envelopment analysis (DEA), principal component analysis (PCA) and numerical taxonomy (NT). Anuchit Thongpun et al. [12] focused on a building location selection approach for locating a solar power plant in Thailand based on DEA models. A. Azadeh et al. [13] applied an artificial neural network (ANN) and fuzzy data envelopment analysis (FDEA) for location optimization of solar plants with uncertainty and complexity.
Amy H. I. Lee et al. [14] integrated MCDM model to set the assurance region (AR) of the quantitative factors, and the AR is incorporated into data envelopment analysis (DEA), additionally adopting a fuzzy analytic hierarchy process (FAHP) for the location of a PV solar plant. Adnan Sozen et al. [15] presented an approach for the location of solar plants by data envelopment analysis (DEA). Ali Azadeh et al. [16] presented an integrated fuzzy DEA model for decision making on wind plant locations. Shinya Yokota et al. [17] applied data envelopment analysis (DEA) for the optimal allocation of mega-solar.
Chabuk et al. [18] discussed two MCDM methods to propose two alternative landfill sites for each area in the studied location, which were AHP and Ratio Scale Weighting (RSW), using multiple environmental factors as criteria for evaluation. The methods successfully proposed a geographical map that addressed the suitable locations throughout the studied sector. Chakraborty et al. [19] also proposed multiple MCDM methods including AHP for assigning criteria weights; Grey Rational Analysis (GRA), Multi Objective Optimization on the Basis of Ratio Analysis (MOORA), Elimination of Choice Translating Reality (ELECTRE II), and Operational Competitiveness Rating Analysis (OCRA) for ranking the results of each location alternative; and Spearman’s rank correlation coefficient, Kendall’s rank correlation coefficient, agreement between the top three ranked alternatives were used to compare each ranking methodology with another; finally, the REGIME was used for the final evaluation of each methodology after all the applied evaluating methodologies were calculated. With a huge comparison study between the methods, the author was able to utilize which methodology was most suitable for choosing the most suitable supplier.
Amy H. I. Lee et al. [20] proposed MCDM model to decide the most suitable photovoltaic solar plant allocation by using the interpretive structural modeling (ISM), the Serbian VlseKriterijumska OptimizacijaI Kompromisno Resenje (VIKOR), meaning multi-criteria optimization with a compromise solution, fuzzy analytic network process (FANP). A. Azadeh et al. [21] integrated hierarchical Data Envelopment Analysis for the location optimization of wind plants in Iran. This model would enable the energy policy makers to select the best possible location for construction of a wind power plant with the lowest possible cost. Lei Fang et al. [22] applied the DEA model and goal programming to evaluate the relative efficiency of each potential location.
Chia-Nan Wang et al. [9] proposed a similar MCDM approach to determine a solar power plant location for the entire country of Vietnam using the methods of DEA, FAHP, and TOPSIS to evaluate qualitative and quantitative criteria. Ehsan Dehghani et al. [23] evaluated different areas for solar plants according to a set of social, geographical, and technical criteria through a data envelopment analysis (DEA) model. In this study, the DEA model considers both information of the efficient and anti-efficient frontiers to raise discrimination power in DEA analysis. Ali Mostafaeipour et al. [24] applied Data Envelopment Analysis (DEA) methodology to prioritize cities for installing the solar-hydrogen power plant so that one candidate location was selected for each city. A. Azadeh et al. [25] presented a technical and economic research for allocation of solar plants by using multivariable methods namely, Data Envelopment Analysis (DEA) and Principle Component Analysis (PCA). A hybrid model for the allocation of solar plants was presented by the utilization of most related parameters to solar plants and an integrated DEA-PCA approach. Seong Kon Lee et al. [26] proposed a hybrid model including a fuzzy Analytic Hierarchy Process (AHP)/Data Envelopment Analysis (DEA) for efficiently allocating energy R&D resources in the case of energy technologies against high oil prices.
As the literature review shows, the amount of studies that apply the MCDM approach to various fields of science and engineering has increased in number in recent years. Location selection is one of the fields where the MCDM model has been employed, especially in the renewable energy sector, where decision makers must evaluate both qualitative and quantitative factors. Although some studies have reviewed the applications of MCDM approaches in solar power plant location selection, very few works have focused on this problem in a fuzzy environment. This is a reason why we proposed a fuzzy MCDM model in this study.

3. Methodology

3.1. Basic Theory

3.1.1. Fuzzy Set

A fuzzy set is used to address the issues that exist in an uncertain environment. It was proposed by Zadeh when attempting to solve multiple criteria decision-making problems. If the set Z ˜ is a triangular fuzzy number (TFN), each value of the membership function is in the range [0, 1]:
μ Z ˜ x = x j k j j x k l x l k k x l 0
Each value of membership, which includes the left (the lowest value to evaluate) and right (the highest value to evaluate) representatives, of a TFN is shown below:
N ˜ = ( N 1 y ,   N r y ) = ( j ( k j )   y ,   l + ( k d )   y ) ,   y [ 0 ,   1 ] .
A triangular fuzzy number can be described as in Figure 4.

3.1.2. Fuzzy Analytic Network Process

Different to the AHP model, a strict hierarchical structure is not required for the ANP model. ANP models allow control elements, that can be controlled by attribute clusters or difference levels. When the elements are at the same level, they can present multiple control factors. A systematic approach in the interaction or feedback between factors is explained using the degree of interaction and interdependence between factors.
Using the rating scale of expert as the basis and data, the decision maker will pair each factor where the weight of each factor is matrixed and determined for the ANP model.
The AHP model then uses a quantitative pair comparison with a priority ratio of 1–9 to set the priority level for each level of the system, where 1 is the lowest priority and 9 is the highest priority. Meanwhile, the ANP model allows making complex relationships between criteria and their rank in the system. The 1–9 scale of the AHP model is shown in Table 1.
The display and goals in the pairwise comparison process, which were considered disadvantages of ANP, were overcome with the development of the FANP model. The FANP model uses a set of values to incorporate decision makers in an uncertain environment, while a crisp value is presented in the ANP model. The Saaty’s model is used to convert the values for the fuzzy prioritization model to easily fix the conversion values, where Oab = ( O ab x , O ab o , O ab v ) is a TFN with the core O ab o , the support [ O ab x , O ab v ], and the TFN, as shown in Figure 5.
Table 2 presents the 1–9 fuzzy conversion scale.
At the reverse level to Oab, expressing the non-preference is also shown by a TFN: (1/ O a b v , 1/ O a b o ,   1 / O a b x ). Using the fuzzy Saaty’s matrix, the weights of the criteria can be determined into four steps that are used to input the data:
  • Using Equations (2)–(4), the fuzzy synthetic extensions K a k a x , k a o , k a v calculation can be transformed into TNT.
    K a = b = 1 n O a b a = 1 n b = 1 n O a b 1
    j = 1 n O a b = b = 1 n O a b x , b = 1 n O a b o , b = 1 n O a b v
    O a b 1 = 1 / O a b v , 1 / O a b o , 1 / O a b x
    O N = ( O x . N x ,   O 0 . N 0 ,   O v . N v )
    Let a = 1, 2, …, n, in which a and b specifically are TFN (Ox, Oo, Ov) and (Nx, N0, Nv) where x is the minimum value, o is the average value, and v is the maximum value;
  • Using relationships with the fuzzy valued for addressing the weight of criteria. In order to determine certain fuzzy extensions, they are calculated by using the minimum fuzzy extension of valued relation ≤ at [5] with weights Qa calculated as follows
    Q a = m i n b k b v k a v k a o k a v k b o k b x
    where a, b = 1, 2, …., n;
  • Standardize the weights. If the decision maker expects to obtain the total weights in one matrix equal to 1, final weights qi are solved by (7)
    q a = Q a / a = 1 n Q a
    where a, b = 1, 2, …, n;
  • An evaluation of a Saaty’s matrix is used to test for its consistency. The matrix of the weights and criteria are consistent and sufficient if inequality of the Consistency Ratio (CR) from Equation (8) is defined as follows using the Consistency Index (CI) and Random Index (RI):
    C R = C I R I = λ ¯ n n 1 × R I 0.1

3.2. Data Envelopment Analysis Model

3.2.1. Charnes–Cooper–Rhodes model (CCR model)

The CCR model, as the fundamental model for the DEA model, is defined as follows:
max c . a   ξ = a Z y 0 c Z x 0 S . t :                      a Z y e c Z x e 0 ,   e = 1 ,   2 , , n a 0 c 0
The defined constraints ensure that the ratio of virtual output to virtual input cannot exceed 1 per decision making unit (DMU). The objective is to obtain a rate of weighted output for every weighted inputs. Subject to the constraints, the optimal goal value ξ* can only reach a maximum of 1.
DMU0 is CCR’s efficient if ξ * = 1 . The result must include a minimum of 1 optima a* > 0 and c* > 0. In addition, the fractional program can be defined as a linear programming problem (LP) as follows:
max c . a ξ = a Z y 0 S . t :                      c Z x 0 1 = 0 a Z y e c Z x e 0 ,   e = 1 ,   2 , , n c 0 a 0
The linear program (10) provides an equal result to the fractional program (9). The linear program from the Farrell model (10) has a variable ξ and a nonnegative vector α = α 1 , α 2 , α 3 , , α f as:
m a x d = 1 m s i + g = 1 o s g + S . t :                      e = 1 n x d e α d + s d = ξ x d 0 ,   b = 1 ,   2 , , p e = 1 n y g e α e s g + = y g 0 ,   g = 1 ,   2 , , o α e 0 ,   e = 1 ,   2 , ,   n s d 0 ,   d = 1 ,   2 , ,   p s g + 0 ,   g = 1 ,   2 , ,   o
The model (3) provides a feasible solution, ξ = 1 , α 0 * = 1 , α j * = 0 , j 0 , in which the optimal solution is affected when ξ * is not higher than 1. A specific DMU is provided when the optimal solution, ξ * , is calculated. For each DMUe, the process will repeat for every e = 1, 2, …, n. When ξ * < 1, the DMUs are inefficient. If ξ * = 1 , the DMUs then are classified as boundary units. By invoking a linear program, we can prevent weakly efficient points d as follows:
m a x d = 1 m s d + g = 1 s s g + S . t :                      e = 1 n x d e α e + s d = ξ x d 0 ,   d = 1 ,   2 , ,   p e = 1 n y g e α e s g + = y r 0 ,   g = 1 ,   2 , ,   o α e 0 ,   e = 1 ,   2 , ,   n s d 0 ,   d = 1 ,   2 , ,   p s g + 0 ,   g = 1 ,   2 , ,   o
For this situation, we clarify that the optimal solution, ξ * , is not affected by the results from s d and s g + .
For both (1) ξ = 1 and (2) s d * = s g + = 0 , DMU0 achieves 100% accuracy and efficiency. For both (1) ξ * = 1 and (2) s d * 0 and s g + 0 for d or g in optimal options, the performance of DMU0 is weakly efficient. Thus, following the development procedure to solve the problem is as follows:
m i n θ μ d = 1 m s d + g = 1 s s g + S . t :                      b = 1 n x d e α e + s d = ξ x d 0 ,   d = 1 ,   2 , ,   p e = 1 n y g e α e s g + = y g 0 ,   g = 1 ,   2 , ,   o α e 0 ,   e = 1 ,   2 , ,   n s d 0 ,   d = 1 ,   2 , ,   p s g + 0 ,   g = 1 ,   2 , ,   o
In this case, s b and s r + variables are first used to transform the inequalities into equivalent equations. For (13), it is the same in terms of methods as when we solve (3) by minimizing ξ in the first stage and then fixing 𝜉 = ξ * as in (4), where the slacks’ variables provide the highest values but the previously determined value of ξ = ξ * is not affected. The objective is converted from maximum to minimum, as in (9), to obtain the following:
max c . a ξ = c Z x 0 a Z y e S . t :                      a Z x 0 c Z y e ,   e = 1 ,   2 , , n c ε > 0 a ε > 0
If the non-Archimedean value and the 𝜀 > 0 are displayed, the input models are similar to models (10) and (13) as follows:
max c . a ξ = c Z x 0 S . t :                      a Z y 0 = 1 c Z x o a Z y e 0 ,   e = 1 ,   2 , , n c ε > 0 a ε > 0
and:
m a x ϕ ε d = 1 m s i + g = 1 s s g + S . t :                      e = 1 n x d e α e + s d = x d 0 ,   d = 1 ,   2 , , p e = 1 n y g e α e s g + = y g 0 ,   g = 1 ,   2 , ,   q α e 0 ,   e = 1 ,   2 , ,   n s d 0 ,   d = 1 ,   2 , ,   p s g + 0 ,   g = 1 ,   2 , ,   o
A dual multiplier model of the CCR input-oriented (CCR-I) is expressed as:
max z = g = 1 q g y g 0 S . t :                      g = 1 o g y r e g = 1 o a g y g e 0 d = 1 p a d x d 0 = 1 c g , a d ε > 0
A dual multiplier model of the CCR output-oriented (CCR-O) is also expressed as:
min o = d = 1 p a d x d 0 S . t :                      d = 1 p a d x d e g = 1 o g y g e 0 g = 1 o g y g 0 = 1 c g , a d ε > 0

3.2.2. Banker Charnes Cooper Model (BCC Model)

The input-oriented BBC model (BCC-I) was introduced by Banker et al., in which the efficiency of DMU0 is assessed by solving the following LP (19):
ξ B = m i n ξ S . t :                      e = 1 n x d e α e + s d = ξ x d 0 ,   d = 1 ,   2 , ,   p e = 1 n y g e α e s g + = y g 0 ,   g = 1 ,   2 , ,   o k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , ,   n
By invoking a linear program, we can prevent the weakly efficient points as in the following:
m a x d = 1 m s d + g = 1 s s g + S . t :                      e = 1 n x d e α e + s d = ξ x d 0 ,   d = 1 ,   2 , ,   p e = 1 n y g e α e s g + = y g 0 ,   g = 1 ,   2 , ,   o k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , ,   n s d 0 ,   b = 1 ,   2 , ,   p s g + 0 ,   g = 1 ,   2 , ,   o
The first multiplicative form to the solve problem is as follows:
m i n ξ ε d = 1 m s d + g = 1 s s g + S . t :                      e = 1 n x d e α e + s d = ξ x i 0 ,   d = 1 ,   2 , ,   p e = 1 n y g e α e s g + = y g 0 ,   g = 1 ,   2 , ,   o k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , ,   n s d 0 ,   b = 1 ,   2 , ,   p s g + 0 ,   r = 1 ,   2 , ,   o
A second multiplier form of the linear program is expressed as:
max c . a , a 0 ξ B = a Z y 0 a 0 S . t :                      c Z x 0 = 1 a Z y e c Z x e a 0 0 ,   e = 1 ,   2 , ,   n c 0 a 0
As mentioned in Formula (14), in this case Z and u are vectors, and the scalar v 0 may be positive or negative or zero. The dual program for the equivalent BCC fractional program (12) can be obtained as follows:
max c . a ξ = a Z y 0 a 0 c Z x 0 S . t :                      a Z y e a 0 c Z x e 1 ,   e = 1 ,   2 , ,   n c 0 a 0
A BCC-efficient solution is a solution where if an optimal solution, DMU0, ( ξ B * ,   s * , s + * ) as solved in this two phase processes for model satisfies ξ B * = 1 and has no slack values where s * = s + * = 0. Otherwise, the model would be considered as BCC-inefficient.
The improved activity, ξ * x s * , y + s + * , also can be claimed as BCC efficient. A DMU, which is a minimized input value for any input item, or a maximized output value for any output item, is BCC-efficient.
The output-oriented BCC model (BCC-O) is defined as the following:
max η S . t :                      e = 1 n x d e α e + s d = ξ x d 0 ,   b = 1 ,   2 , ,   p e = 1 n y g e α e s g + = η y g 0 ,   g = 1 ,   2 , ,   o k = 1 n α k = 1 α k 0 ,   k = 1 ,   2 , ,   f
A multiplier form of the linear program (24) can be expressed as [25]:
min c . a , c 0 a Z y 0 a 0 S . t :                      a Z y 0 = 1 c Z x e a Z y e a 0 0 ,   e = 1 ,   2 , ,   n c 0 a 0
In the envelopment model, the Z0 is the scalar combined with k = 1 n α k = 1 . Conclusively, the equivalent fractional programming formulation for the BCC model was achieved by the authors (25):
min c . a , c 0 c Z x 0 c 0 a Z y 0 S . t :                      c Z x e a 0 a Z y e 1 ,   e = 1 ,   2 , ,   n c 0 a 0

3.2.3. Slacks Based Measure Model (SBM Model)

The SBM model was developed by Tone, which has three elements, input-oriented, output-oriented. Input-Oriented SBM (SBM-I-C).
The following model can be defined as the Input-oriented SBM under constant-returns-to-scale-assumption:
ρ I * = min α ,   s , s + 1 1 m d = 1 m s d x d h S . t :                      x d c = e = 1 m x d c α d + s d ,   d = 1 ,   2 , ,   p y g c = e = 1 m y g c α e s g + ,   g = 1 ,   2 , ,   o α e 0 ,   k   j ,   s d 0   e , s g + 0   e

4. Case Study

Solar power in Vietnam belongs to the emerging energy industry group, which is engaged in the development of the world’s renewable energy sources, the import of science and technology, and meets the demand for developing power sources. When large hydropower sources are fully exploited, small hydroelectricity will not guarantee benefits compared to the environmental damage. Vietnam, on the other hand, has great potential for solar and wind power, due to its proximity to the equator and the existence of dry, sunny regions like the southern central provinces. The study results are shown in Figure 6. Therefore, solar power, together with wind power, is being encouraged to develop by the State of Vietnam, reflected in the Prime Minister’s Decision No. 2068/QD-TTg of 25 November 2015, approving the Renewable Energy Development Strategy. Vietnam’s electricity generation until 2030, with a vision of 2050, ensures the development of electricity sources when stopping nuclear power projects and reducing fossil-fired thermal power plants.
With the growth of economic sectors in general and agriculture, the demand for electricity in the Mekong Delta provinces is increasing. However, fossil fuel for electricity generation is declining, so the development of renewable energy is one of the sustainable solutions for the Mekong Delta to have electricity to ensure domestic living and economic development.
In this research, the authors present an integrated approach for the site of solar plants in Mekong Delta, Vietnam, by using data envelopment analysis (DEA) and Fuzzy Analytical Network Process (FANP). In the first stage, DEA model is applied to select some potential location, then FANP model is used for ranking these potential locations. The authors collected data from thirteen locations in Mekong Delta, Vietnam, which can invest in solar power plants, as shown in Table 3.
DEA model is a quantitative technique that determines the relative effectiveness of multiple inputs and outputs decision makers. Halasah et al. [28] employed life-cycle assessment to evaluate the energy-related impacts of PV systems at different scales of integration. The input parameters included panel efficiency, temperature coefficient, shading losses, ground cover ratio and latitude, and the input data included hourly solar radiation, wind speed and temperature. Wang and Amy [9,14] using DEA model for ranking potential location for building solar power plant. In their research, the output data included sunshine hour and elevation. Due to the information accessibility of various sites and the importance of various factors, we select two inputs and two outputs for the quantitative factors. The two inputs are temperature (I1) and wind speed (I2). The two outputs are sunshine hours (O1) and elevation (O2) [14,20]. The definition of the inputs and outputs are defined in Table 4 [20]. Raw data of inputs and outputs of DMUs are demonstrated in Figure 7.
For applied DEA model, some additional data about 13 locations in Mekong Delta, Vietnam are show in Table 5.
To select some potential sites in Mekong Delta, Vietnam, there are several DEA models, including SBM-I-C; CCR-I; BCC-O; CCR-O and BCC-I, applied in this step. The results of the DEA models are shown in Table 6.
As the results in Table 6 show that there are seven DMUs that are potential locations for building a solar power plant in Mekong Delta, Vietnam including DEL01, DEL03, DEL05, DEL09 and DEL11. These DMUs will be evaluated in the next step of this research by using the FANP model.
In the final stage, all the potential locations will be ranked by the FANP model. Some criteria affecting the location selection are shown in Figure 8.
Fuzzy comparison matrix of EC from the FANP model is shown in Table 7.
The fuzzy numbers were converted to real numbers by using the TFN. During the defuzzification, the authors obtain the coefficients α = 0.5 and β = 0.5. In this, α represents the uncertain environment conditions, and β represents whether the attitude of the evaluator is fair.
g 0.5 , 0.5 a E L , S L ¯   = [ ( 0.5   ×   4.5 ) + ( 1     5.5 )   ×   2.5 ] = 5
f0.5(LEL,SL) = (5 − 4) × 0.5 + 4 = 4.5
f0.5(USL,SL) = 6 − (6 − 5) × 0.5 = 5.5
g 0.5 , 0.5 a S L , E L ¯   = 1 / 5
The remaining calculations for other criteria are similar to the above calculation. The real number priority when comparing the main criteria pairs is shown in Table 8.
The following are used for calculating the maximum individual value:
P1 = (1 × 6 × 2 × 4)1/4 = 2.63
P2 = (1/5 × 1 × 1/4 × 3)1/4 = 0.62
P3 = (1/2 × 4 × 1 × 5)1/4 = 1.78
P4 = (1/4 × 1/3 × 1/5 × 1)1/4 = 0.36
Q = 5.39
ω 1 = 2.63 5.39 = 0.49
ω 2 = 0.62 5.39 = 0.12
ω 3 = 1.78 5.39 = 0.33
ω 4 = 0.36 5.39 = 0.07
1 6 2 4 1 / 6 1 1 / 4 3 1 / 2 4 1 5 1 / 4 1 / 3 1 / 5 1 × 0.49 0.12 0.33 0.07 = 2.15 0.49 1.41 0.30
2.15 0.49 1.41 0.30 / 0.49 0.12 0.33 0.07 = 4.39 4.08 4.27 4.29
Based on the number of main criteria, the authors get n = 5; λmax and CI are calculated as follows:
λ m a x = 4.39 + 4.08 + 4.27 + 4.29 4 = 4.23
C I = λ m a x n n 1 = 4.23 4 4 1 = 0.077
To calculate CR value, we get RI = 0.9 with n = 4.
C R = C I R I = 0.077 0.9 = 0.08556
As CR = 0.08556 ≤ 0.1, we do not need to re-evaluate.

5. Discussion

There have been many studies presenting MCDM models using fuzzy set theory to evaluate and select the location of a renewable energy plant, particularly solar energy. However, the literature review shows that most of these studies only use economic, or technological, or site characteristics, or environmental or social standards to assess potential locations. There are few studies that simultaneously use economic, technological, site characteristics, environmental and social factor standards in the process of location selection. In addition, there seems to be no research on developing an MCDM model to choose a location for a solar power plant in the Mekong Delta.
This is a reason the author proposed an MCDM model including data envelopment analysis (DEA) and fuzzy analytical network process (FANP) for solar power plant location selection in Mekong Delta, Vietnam in this work. In the first step, some potential location is selected from 13 location in Mekong Delta by using the DEA model. Then, the FANP model is applied for rank potential sites. In a renewable energy project, location decision is the most importance activity; location decision is always a multicriteria decision making factor, as the decision maker must consider both quantitative and qualitative factors. In any renewable energy project, location decision is a very important decision, decision maker must consider qualitative and quantitative factors. Thus, solar power plant location selection is a multicriteria decision making (MCDM). For selection optimal location, the proposed model considered five main factors: economic factors (construction cost, operation and management cost, new feeder cost), technological factors (distance from major road, distance from power network, potential demand), site characteristics factors (ecology, elevation, approachability), environmental factors (temperature, sunshine hours, humidity) and social factors (supports mechanisms, protection law and legal regulatory compliance) As per the results shown in Table 9, Ben Tre (DEL03) is the optimal solution for investing solar power plants in Mekong Delta, Vietnam.

6. Conclusions

Facing energy demand challenges to ensure sustainable economic growth, especially the exhaustion of fossil fuels, environmental pollution, and climate change. Incentive policies to take advantage of renewable energy supply opportunities in Vietnam. This approach to taking advantage of this new energy source not only contributes to timely supply of the energy needs of the society but also helps to save electricity and reduce environmental pollution.
However, when considering investing in solar power plants, decision makers must consider many qualitative and quantitative factors including economic, technological, site characteristics, environmental and social factors. Therefore, this problem is tackled as a Multi-Criteria Decision Making (MCDM) problem that requires that both quantitative and qualitative information are considered. This paper presented an integrated approach for site of solar plants by data envelopment analysis (DEA) and Fuzzy Analytical Network Process (FANP). Furthermore, this integrated approach, incorporating the most relevant parameters of solar plants, is introduced.
The contributions of this research include a fuzzy multicriteria decision-making (F-MCDM) approach for solar power plant site selection in Vietnam. This research also utilizes the evolution of a new approach that is flexible and practical for the decision-maker and provides useful guidelines for solar power plant location selection in many countries around the world.
For further studies regarding this topic, the study can be expanded to other MCDM approaches such as VIKOR, ELECTRE I, ELECTRE III.

Author Contributions

Conceptualization, C.-N.W., T.P.H., V.T.N. and S.T.H.; Data curation, V.T.H.V., T.P.H. and V.T.N.; Formal analysis, C.-N.W., V.T.H.V., T.P.H., V.T.N. and S.T.H.; Funding acquisition, C.-N.W. and V.T.H.V.; Investigation, C.-N.W., T.P.H., V.T.N. and S.T.H.; Methodology, C.-N.W., V.T.H.V., T.P.H. and V.T.N.; Project administration, V.T.H.V., T.P.H. and V.T.N.; Resources, V.T.N.; Visualization, T.P.H. and S.T.H.; Writing—original draft, V.T.H.V., T.P.H.,V.T.N. and S.T.H.; Writing—review and editing, C.-N.W., V.T.H.V. and V.T.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Renewable Capacity Growth between 2019 and 2024 by Technology. Available online: https://www.iea.org/data-and-statistics/charts/renewable-capacity-growth-between-2019-and-2024-by-technology (accessed on 22 March 2020).
  2. Ưu Khuyết Điểm Của Năng Lượng Mặt Trời. Available online: https://www.iea.org/data-and-statistics/charts/renewable-capacity-growth-between-2019-and-2024-by-technology (accessed on 10 March 2020).
  3. World Summit Outcome, Resolution A/60/1, adopted by the General Assembly on 15 September 2005. Available online: https://www.un.org/en/development/desa/population/migration/generalassembly/docs/globalcompact/A_RES_60_1.pdf (accessed on 5 March 2020).
  4. MRC Strategic Environmental Assessment (SEA) of Hydropower on the Mê Công Mainstream. Available online: http://www.mrcmekong.org/assets/Publications/Consultations/SEA-Hydropower/SEA-FR-summary-13oct.pdf (accessed on 5 March 2020).
  5. Đồng Bằng Sông Cửu Long Trước Nguy Cơ Lớn Từ Nhiệt điện Than. Available online: https://www.thiennhien.net/2018/06/22/dong-bang-song-cuu-long-truoc-nguy-co-lon-tu-nhiet-dien-than (accessed on 22 March 2020).
  6. Gorjian, S.; Hashjin, T.T.; Ghobadian, B. Estimation of mean monthly and hourly global solar radiation on surfaces tracking the sun: Case study: Tehran. In Proceedings of the Second Iranian Conference on Renewable Energy and Distributed Generation, Tehran, Iran, 6–8 March 2012; pp. 172–177. [Google Scholar]
  7. Điện Mặt Trời Phát Triển Mạnh Ở Đồng Bằng Sông Cửu Long. Available online: http://ceid.gov.vn/dien-mat-troi-phat-trien-manh-o-dong-bang-song-cuu-long (accessed on 22 March 2020).
  8. Gloal Solar Atlas. Available online: https://globalsolaratlas.info/map?c=-29.22889,-1.40625,2 (accessed on 22 March 2020).
  9. 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. [Google Scholar] [CrossRef] [Green Version]
  10. Kengpol, A.; Rontlaong, P.; Tuominen, M. A decision support system for selection of solar power plant locations by applying fuzzy ahp and topsis: An empirical study. J. Softw. Eng. Appl. 2013, 6, 470–481. [Google Scholar] [CrossRef] [Green Version]
  11. Azadeh, A.; Ghaderi, S.F.; Maghsoudi, A. Location optimization of solar plants by an integrated hierarchical DEA PCA approach. Energy Policy 2008, 36, 3993–4004. [Google Scholar] [CrossRef]
  12. Thongpun, A.; Nasomwart, S.; Peesiri, P.; Nananukul, N. Decision support model for solar plant site selection. In Proceedings of the 2017 IEEE International Conference on Smart Grid and Smart Cities (ICSGSC), Singapore, 23–26 July 2017; pp. 50–54. [Google Scholar]
  13. Azadeh, A.; Sheikhalishahi, M.; Asadzadeh, S.M. A flexible neural network-fuzzy data envelopment analysis approach for location optimization of solar plants with uncertainty and complexity. Renew. Energy 2011, 36, 3394–3401. [Google Scholar] [CrossRef]
  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 solar plant. Sustainability 2015, 7, 13522–13541. [Google Scholar] [CrossRef] [Green Version]
  15. Sozen, A.; Mirzapour, A.; Çakir, M.T. Selection of the best location for solar plants in Turkey. J. Energy S. Afr. 2015, 26, 52–63. [Google Scholar] [CrossRef] [Green Version]
  16. Azadeh, A.; Rahimi-Golkhandan, A.; Moghaddam, M. Location optimization of wind power generation–transmission systems under uncertainty using hierarchical fuzzy DEA: A case study. Renew. Sustain. Energy Rev. 2014, 30, 877–885. [Google Scholar] [CrossRef]
  17. Yokota, S.; Kumano, T. Mega solar optimal allocation using data envelopment analysis. Electr. Eng. Jpn. 2013, 183, 24–32. [Google Scholar] [CrossRef]
  18. Chabuk, A.; Al-Ansari, N.; Hussain, H.M. Landfill sites selection using MCDM and comparing method of change detection for Babylon Governorate. Environ. Sci. Pollut. Res. 2019, 26, 35325–35339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Chakraborty, R.; Ray, A.; Dan, P. Multi criteria decision making methods for location selection of distribution centers. Int. J. Ind. Eng. Comput. 2013, 4, 491–504. [Google Scholar] [CrossRef]
  20. Lee, A.H.; 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]
  21. Azadeh, A.; Ghaderi, S.F.; Nasrollahi, M.R. Location optimization of wind plants in Iran by an integrated hierarchical data envelopment analysis. Renew. Energy 2011, 36, 1621–1631. [Google Scholar] [CrossRef]
  22. Fang, L.; Li, H. Multi-criteria decision analysis for efficient location-allocation problem combining DEA and goal programming. RAIRO Oper. Res. 2015, 49, 753–772. [Google Scholar] [CrossRef]
  23. Dehghani, E.; Jabalameli, M.S.; Pishvaee, M.S.; Jabarzadeh, A. Integrating information of the efficient and anti-efficient frontiers in DEA analysis to assess location of solar plants: A case study in Iran. J. Ind. Syst. Eng. 2018, 11, 163–179. [Google Scholar]
  24. Mostafaeipour, A.; Sedaghat, A.; Qolipour, M.; Rezaei, M.; Arabnia, H.R.; Saidi-Mehrabad, M.; Alavi, O. Localization of solar-hydrogen power plants in the province of Kerman. Adv. Energy Res. 2017, 5, 179. [Google Scholar]
  25. Azadeh, A.; Ghaderi, S.F.; Maghsoudi, A. Location optimization of solar plants by an integrated multivariable DEA-PCA method. In Proceedings of the 2006 IEEE International Conference on Industrial Technology, Mumbai, India, 15–17 December 2006; pp. 1867–1872. [Google Scholar]
  26. Lee, S.K.; Mogi, G.; Hui, K.S. A fuzzy analytic hierarchy process (AHP)/data envelopment analysis (DEA) hybrid model for efficiently allocating energy R&D resources: In the case of energy technologies against high oil prices. Renew. Sustain. Energy Rev. 2013, 21, 347–355. [Google Scholar]
  27. Tuan, L.A. An overview of the renewable energy potentials in the Mekong river Delta, Vietnam. Can Tho Univ. J. Sci. 2016, 4, 70–79. [Google Scholar]
  28. Halasah, S.A.; Pearlmutter, D.; Feuermann, D. Field installation versus local integration of photovoltaic systems and their effect on energy evaluation metrics. Energ. Pol. 2013, 52, 462–471. [Google Scholar] [CrossRef]
  29. Radziemska, E. The effect of temperature on the power drop in crystalline silicon solar cells. Renew. Energy 2003, 28, 1–12. [Google Scholar] [CrossRef]
  30. Solar Panel Temperature Affects Output—Here’s What You Need to Know. Available online: http://www.solar-facts-and-advice.com/solar-panel-temperature.html (accessed on 5 March 2020).
  31. Guide to Meteorological Instruments and Methods of Observation. Available online: http://www.wmo.int/pages/prog/gcos/documents/gruanmanuals/CIMO/CIMO_Guide-7th_Edition2008.pdf (accessed on 5 March 2020).
  32. Panjwani, M.; Narejo, D. Effect of altitude on the efficiency of solar panel. Int. J. Eng. Res. Gen. Sci. 2014, 2, 461–463. [Google Scholar]
Figure 1. Comparison of different renewable energy sources [4].
Figure 1. Comparison of different renewable energy sources [4].
Energies 13 04066 g001
Figure 2. Average amount of potential solar energy output [6].
Figure 2. Average amount of potential solar energy output [6].
Energies 13 04066 g002
Figure 3. Solar Radiation Maps of Vietnam [9].
Figure 3. Solar Radiation Maps of Vietnam [9].
Energies 13 04066 g003
Figure 4. A triangular fuzzy number (TFN).
Figure 4. A triangular fuzzy number (TFN).
Energies 13 04066 g004
Figure 5. A triangular fuzzy number (TFN).
Figure 5. A triangular fuzzy number (TFN).
Energies 13 04066 g005
Figure 6. Irradiation Map of Mekong Delta, Vietnam [27].
Figure 6. Irradiation Map of Mekong Delta, Vietnam [27].
Energies 13 04066 g006
Figure 7. Inputs and outputs of DEA models.
Figure 7. Inputs and outputs of DEA models.
Energies 13 04066 g007
Figure 8. Main criteria and sub-criteria in Fuzzy Analytical Network Process (FANP) model.
Figure 8. Main criteria and sub-criteria in Fuzzy Analytical Network Process (FANP) model.
Energies 13 04066 g008
Table 1. The 1–9 scale [6].
Table 1. The 1–9 scale [6].
Importance IntensityDefinition
1.Equal significance
3.Moderate significance
5.Strongly higher significance
7.Very strong higher significance
9.Extremely high significance
2, 4, 6, 8.Intermediate values
Table 2. The 1–9 fuzzy conversion scale.
Table 2. The 1–9 fuzzy conversion scale.
Importance Intensity.Triangular. Fuzzy Scale.
1.(1.0, 1.0, 1.0).
2.(1.0, 1.0, 2.0).
3.(1.0, 2.0, 3.0).
4.(2.0, 3.0, 4.0).
5.(3.0, 4.0, 5.0).
6.(4.0, 5.0, 6.0).
7(5.0, 6.0, 7.0).
8(7.0, 8.0, 9.0).
9(9.0, 9.0, 9.0).
Table 3. List of the 13 locations identified in Mekong Delta, Viet Nam.
Table 3. List of the 13 locations identified in Mekong Delta, Viet Nam.
NoLocationsSymbolsNoLocationsSymbols
1An GiangDEL018Kien GiangDEL08
2Bac LieuDEL029Long AnDEL09
3Ben TreDEL0310Soc TrangDEL10
4Ca MauDEL0411Tien GiangDEL11
5Can ThoDEL0512Tra VinhDEL12
6Dong ThapDEL0613Vinh LongDEL13
7Hau GiangDEL07---
Table 4. Definition of inputs and outputs.
Table 4. Definition of inputs and outputs.
FactorsDefinition
Inputs
Temperature (I1)This is a measure of how high or low the heat radiation of the current environment is. Radziemska [29] discussed that should the temperature increase, the overall output power and the conversion efficiency of the PV module would decrease. A solar panel that receives more sunlight would lead to an increase in temperature on the panel, hence a decrease in the conversion efficiency. This decrease can vary from 10% to 25% depending on the location and material of the equipment. [30]
Wind speed (I2)Wind is the current of gas particles that flows randomly. Depending on the strength of the current, small particles from solid sediments to suspended objects could obstruct and affect the solar panels and other equipment. The wind can cause erosion and potential failures in operation for solar power plants.
Outputs
Sunshine hours (O1)An indicator that measures the average duration of sunshine in a time period (here, annually) for a given location. This period is when direct solar irradiance shows a measurement higher than 120 W/m2 [31]. A higher output power would be produced should a higher sunshine hour be recorded.
Elevation (O2)The geographical height above sea level of a given location (here, meters). If the geographical elevation is high, the distance between the ground to the sun is decreased, which means the solar radiant takes a shorter time to reach the ground. Additionally, lower distance means a lower dispersion of solar radiation, so the intensity is projected to be higher. Higher intensity yields higher solar energy output. Panjwani and Narejo discussed how elevation generated a 7–12% increase in power by testing 3 solar panels at a 27.432 m elevation [32].
Table 5. Raw data of DEA model.
Table 5. Raw data of DEA model.
DMUsTemperature (I1)Wind Speed (I2)Sunshine hours (O1)Elevation (O2)
DEL0127.201.7025892.00
DEL0227.307.2027001.50
DEL0327.001.0023412.00
DEL0426.801.3023001.00
DEL0526.902.8026722.00
DEL0627.402.5025002.00
DEL0727.307.0027000.80
DEL0827.506.0025631.20
DEL0927.702.8028001.60
DEL1026.806.4022920.90
DEL1127.001.7026451.00
DEL1227.606.8025560.80
DEL1328.001.4024001.60
Table 6. Ranking results from some DEA models.
Table 6. Ranking results from some DEA models.
LocationsDEA Approaches
SBM-I-CCCR-IBCC-OCCR-OBCC-I
DEL0111111
DEL020.65460.97840.98680.97840.9918
DEL0311111
DEL040.87280.9360.99990.9361
DEL0511111
DEL060.77670.982410.98240.9824
DEL070.64510.97840.98680.97840.9918
DEL080.61590.9220.92590.92200.9771
DEL0911111
DEL100.55210.84610.99610.84611
DEL1111111
DEL120.59350.91620.91810.91620.9735
DEL130.86040.93250.95950.93250.962
Table 7. Fuzzy comparison matrices for Economic factor (EC).
Table 7. Fuzzy comparison matrices for Economic factor (EC).
Technological (TL)Site Characteristics (ST)Environmental (EL)Social (SL)
TL(1,1,1)(5,6,7)(1,2,3)(3,4,5)
ST(1/7,1/6,1/5)(1,1,1)(1/3,1/4,1/5)(2,3,4)
EL(1/3,1/2,1)(5,4,3)(1,1,1)(4,5,6)
SL(1/5,1/4,1/3)(1/4,1/3,1/2)(1/6,1/5,1/4)(1,1,1)
Table 8. Real number priority.
Table 8. Real number priority.
TLSTELSL
TL1624
ST1/611/43
EL1/2415
SL1/41/31/51
Table 9. Final ranking from the FANP model.
Table 9. Final ranking from the FANP model.
NoLocation-SymbolsScoreRanking Order
1An Giang-DEL010.2263
2Ben Tre-DEL030.2461
3Can Tho-DEL050.2382
4Long An-DEL090.1074
5Tien Giang-DEL110.1835

Share and Cite

MDPI and ACS Style

Wang, C.-N.; Viet, V.T.H.; Ho, T.P.; Nguyen, V.T.; Husain, S.T. Optimal Site Selection for a Solar Power Plant in the Mekong Delta Region of Vietnam. Energies 2020, 13, 4066. https://doi.org/10.3390/en13164066

AMA Style

Wang C-N, Viet VTH, Ho TP, Nguyen VT, Husain ST. Optimal Site Selection for a Solar Power Plant in the Mekong Delta Region of Vietnam. Energies. 2020; 13(16):4066. https://doi.org/10.3390/en13164066

Chicago/Turabian Style

Wang, Chia-Nan, Van Tran Hoang Viet, Thanh Phong Ho, Van Thanh Nguyen, and Syed Tam Husain. 2020. "Optimal Site Selection for a Solar Power Plant in the Mekong Delta Region of Vietnam" Energies 13, no. 16: 4066. https://doi.org/10.3390/en13164066

APA Style

Wang, C. -N., Viet, V. T. H., Ho, T. P., Nguyen, V. T., & Husain, S. T. (2020). Optimal Site Selection for a Solar Power Plant in the Mekong Delta Region of Vietnam. Energies, 13(16), 4066. https://doi.org/10.3390/en13164066

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