Priority ranking for energy resources in Turkey and investment planning for renewable energy resources

Energy has been an essential factor in determining the governments’ policies. The countries had to produce their own energy to decrease dependency on external resources. That also provided to gain a great importance on investment on power plants. In this study, a multi-objective Mixed-Integer Linear Programming (MILP) model synchronously optimizing five targets determined as decreasing carbon dioxide (CO2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{2})$$\end{document} emission, increasing power consumption, increasing power plants, increasing energy generation and installed capacity was used. In described model, it was solved by considering renewable power plants in Turkey and fossil fuel-based power plants having most share in Turkey. By trying to minimize deviation values of Turkey’s 2023 targets, it aimed to determine which power plants need to be increased. To determine the priorities of these targets, Ranking Approach for fuzzy numbers by Liou and Wang (1992) was used. Besides, Fuzzy Analytic Hierarchy Process (AHP) was used to prioritize investment planning of renewable power plants in Turkey and five different kinds of power plants under the visual pollution criteria based on amount of CO2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{2}$$\end{document} emission released, environmental damage, capital costs, space requirement and provided employment were evaluated.


Introduction
Energy is one of the most essential factors for a country's social welfare and sustainable development. Recently, population growth, industry development, and income growth, especially in developing countries, cause the energy need increase.
In view of the distribution of power plants in Turkey and their effects of technologic, economic, socioeconomic and life quality, decision makers have to choose the best alternatives to accomplish their goals to decide which power plant will be built [2].
In this article, Fuzzy Analytic Hierarchy Process (AHP), which is one of the multiple criteria decision-making methods, was used. By considering environmental damage of power plants, capital costs, space requirement, their advantages and used technological development under the main criteria of social acceptance, the power plants in Turkey were ranked and an order of preference was composed. In addition, how many power plants should be built was determined with fuzzy-weighted goal programming method by considering Turkey's installed capacity, number of power plants, energy generation and targets of demands in 2023. The study consisted of an introduction, literature search, multiple criteria decision-making methods and an application. Analytic Hierarchy Process and goal programming that are multiple criteria decision-making methods and fuzzy conditions of these methods were described. Later, it continued with an application and the results were included.

Literature search
Kowalski et al. worked on determining economically best option for sustainable energy using multiple criteria decision-making methods for several gas emissions released by power plants. According to the results obtained, natural gas was determined as the best fuel [16]. Pilavachi et al. made a multiple criteria evaluation on nine types of power plants using hydrogen and natural gas as a fuel with AHP method under seven criteria. In this study, combined cycle plant with natural gas was determined as a better alternative [24]. While choosing the best alternative among renewable power plants in Turkey, Kahraman et al. used AHP method. At the end of the study, wind power plant was the best alternative for renewable energy [13]. Liu [20]. Shen et al. applied AHP method for determining economically and environmentally the energy value of investment selection for renewable energy resources within the frame of Taiwan Government's sustainable energy policies [29]. Kahraman and Kaya used fuzzy AHP method and VIKOR method to determine the choice of best alternative power plants forİstanbul [15]. Wang et al. applied AHP method for energy resources selection in China [32]. Lee et al. used fuzzy AHP in making a strategic route map of energy technologies to use alternative energy technologies against high oil prices [17]. Çebi et al. used fuzzy AHP and axiomatic design techniques for choosing the best alternative for investment and comparing renewable energy resources [14]. Jeberaj and Iniyan applied fuzzy AHP and AHP method for the choice of sustainable energy [11]. Liang et al. used the general fuzzy theory for energy planning, supply demand, assumption and renewable energy modeling. In the article, especially on the basis of the current situation, they improved a model for evaluating manufacturing projects about the usage of Analytic Hierarchy Process (AHP) and fuzzy evaluation. Later, while choosing manufacturing projects in China, a sample study for the evaluation of optimum manufacturing project showing the activity of the model was presented [18]. Anane et al. explained an innovator method to range the proper lands for watering using fuzzy AHP based on geographical information systems in the basin of TWW Nabeul-Hammamet (Tunisia) chosen as a target region [1]. Nixon et al. evaluated the main current collection of thermal solar energy technologies within the frame of AHP. These technologies were compared with technical, economic and environmental criteria [22]. Jinturkar and Deshmukh. Improved a fuzzy mixed-integer goal programming method for rural sintering and heating energy planning in the center of India [12]. Daim et al., improved a fuzzy goal programming in the state of Oregon to make a renewable portfolio giving an answer to 25% of energy demand obtained by renewable resources in 2025 [5]. San Cristobal applied a goal programming method for layout of five different renewable power plants in the region of Cantabria in Northern Spain by achieving seven targets about maintenance, management, investment, distances between power plants and energy generation [28]. Ramanathan and Ganesh used a goal programming model to ideally criticize seven energy resources used for lighting at twelve targets-driven houses representing energy, economic and environmental systems [26]. By considering financial and profit risks, Han et al. offered a multi-purpose optimization model for CO 2 minimization infrastructure design and sustainable energy generation [7]. Chang recommended a multiple choice goal programming model to deal with the problem of capacity enhancement plan at the renewable energy sector [4]. Jayaraman et al. offered a goal programming model integrating plentiful allotments of labor resources to reach sustainable, economic, energy and environmental goals of United Arab Emirates [10].

Multiple criteria decision making
Multiple criteria decision making (MCDM), with the briefest definition, is a general name given for problem-solving of multiple and conflicted criteria to be achieved. MCDM explains a top-concept including designed techniques and methods to help people facing the problems being characterized by different size of criteria, and multiple and conflicted criteria [33]. The methods of MCDM are divided into two as multiple criteria and multi-objective. Multi-objective decision making is a model defined by the alternatives as a mathematic model. Multiple criteria decision making is an evaluation process using many criteria taking generally the value of different criteria and weighted, conflicted and even qualitative values on the purpose of eliminating, prioritizing, classifying, sorting and selecting a finite number of the options.

Goal programming (GP)
Goal programming was extended by Ijiri in the middle of 1960s, but it was improved by Ignizio and Lee in 1970s [8,9]. With a different viewpoint from linear programming, goal programming, instead of minimizing or maximizing an objective function, minimizes the deviations from the targets determined within the frame of current limitations. These deviations are shown as a negative deviation and a positive deviation. Since the ratio of these variables' objective functions is 0, they cannot affect the results. Namely, the purpose of the problem in goal programming is to minimize the sum of variables showing the deviation [31]. In addition to this, goal programming combines multiple goals conflicting with decision maker's options. The targets determined in the results may not be reached, so even if there are no optimal results, the most acceptable ones can be obtained. In our article, a preemptive goal programming model which is a variety of goal programming was used. Preemptive GP formulation: are priority values being associated with negative and positive deviations. The numerical value of the target achieved by decision maker is showed with g i .The number of the targets determined by decision maker is showed with p. While being done an order of priority of a hierarchical structure by decision makers, Liou and Wang's Sorting Approach determining to reach different goals for fuzzy numbers was used here [19]. Linguistic variables are given for weight/importance of decision makers and different goals in Table 1.
As showed in Table, these linguistic variables are characterized by triangle fuzzy numbers. In the method of Liou and Wang's total integral value, a ∈ [0, 1] as an optimism index; for fuzzy numbers given as Ã = (a, b, c), total integral value is calculated in this way [19].
. These inverse functions are formulated like the following equation.
While a = 0, the total integral value represents optimistic decision maker and it is calculated in the following equation.
Total integral value of a = 0.5 represents moderate decision maker and it is calculated in the following equation.
Total integral value of a = 0.5 represents optimistic decision maker and it is calculated in the following equation.
I a T Ã = a ik ith for decision maker and j. for fuzzy goal are the performances [19].

Analytic Hierarchy Process (AHP)
AHP is a decision-making technique determining the order of importance by finding the priorities according to each other's criteria and making paired comparisons with objective and subjective criteria [30]. In these paired comparisons, it is preferred in terms of which one of them is more important than the other. By determining them, it is based on a numerical  evaluation of them. AHP enables to make an order between options as well as determining the best option for a person who is about to make a decision. For the reason that this method which considers both quantitative and qualitative factors is used widely and is applied simply, it is applied easily even in the most complicated problems. In that being widely and flexible, AHP makes it a great convenient [6].

Fuzzy Analytic Hierarchy Process (Fuzzy AHP)
An enhanced fuzzy AHP method suggested by Chang was used in many problems, which fuzzy AHP was used. In this method, the cutting levels of "a" were not necessary. Besides using the values of artificial ratings, this method comes to the forefront with simple level sequencing and integrated sequencing. The most advantageous side of this method is that calculation requirement is low and it does not need any additional process by following the steps of classical AHP. The disadvantage of it is that it only uses fuzzy triangle numbers [6]. Pairwise comparisons matrices are arranged to determine the weights of criteria and these comparisons   Step 2 M 1 = (l 1 , m 1 , u 1 ) and M 2 = (l 2 , m 2 , u 2 ) are two triangle fuzzy numbers. The degree of probability is defined below; And expressed below. The intersection of triangle fuzzy numbers is given in Fig. 1   Fig. 1 M 1 and M 2 the intersection of triangle fuzzy numbers [3] Step 3 That the degree of probability of a convex fuzzy number is bigger than the convex number of M I (I = 1, 2, . . . , k) is defined below. For k = 1, 2, . . . , n; k = i, calculated as d , the weighting vector is obtained below.
Step 4 When normalizing the weighting vector given above, The vector above is obtained. Now, this W weighting vector is not a fuzzy number [27].

The Preemptive Goal Programming application
In our study, seven power plants were analyzed; S 1 as a coal plant, S 2 as a natural gas combined cycle power plant, S 3 as a hydropower plant, S 4 as wind plant, S 5 as a geothermal plant, S 6 as a solar plant and S 7 as a biomass plant were defined. Turkey's installed capacity, energy generation, energy demand and greenhouse gas emission defined as its sectorial indicators were divided into the numbers of power plants and the obtained values are given in Table 5 [21]. It was calculated using (3.4), (3.5) and (3.6)'s priority values of goals and fuzzy evaluation matrix in Table 4 below.
According to the results, in view of the numbers of the available power plant, the power plants that are required to be built must have been the solar power plant with the unit 834 and wind power plant with the unit 572. The first goal shows that the deviation in d 12 stayed under the targeted installed capacity. d 22 , i.e., high deviation value at energy generation, is shown that energy need will be increased in the long term until 2023. Non-renewable fossil resources cannot supply energy demands in the long term. According to 2023 targets, the deviation demand in d 32 will stay under 47,92500 toe (tone of oil equivalent). The deviation in fourth goal, d 42  Tone-CO 2 /Gwh and these are needed to be reduced using renewable energy resources.

Fuzzy AHP application
After the application of project investment, according to the model results of available energy alternatives in Turkey, it was seen to make an investment in renewable energy resources. In the second stage, it was analyzed which renewable energy resources were made an investment on and it was ordered. The used power plants are below. To compare these power plants, at first, the vectors were determined. Later, the vectors were compared dynamically with a point scoring system obtained with an expert opinion working at department of energy. Here below, it was given the statement of fuzzy values of the pairwise comparisons matrix. The pairwise comparison matrix among the criteria is given in Table 7.
According to the matrix above, the synthesis value of each criterion was calculated. The Eq. (3.9) was used (Table 8).
Using the values obtained, the comparisons of fuzzy values were made with Eq. (3.14) and the values below were found. To calculate the vector in Table 9, the minimum one of priority values related to the alternatives from the obtained values is taken.
Priority vector was given below. W = (0,06; 0,08; 0,23; 0,79; 1) When the result of the calculation of priority vector, the vector below is obtained. W = (0,365; 0,463; 0,027; 0,106; 0,037) After determining the weights belonging to the criteria, taking an expert opinion working at an energy company and verbalization of numeric data of criteria, pairwise comparisons of five power plants were made. After all calculations, by making a matrix applied by multiplication the criteria weights and alternatives weights, and total superiority weights of alternatives are calculated. In Table 10, the weights of alternatives were shown according to each criterion. Main criteria weights are given in Table  11. Total superiority weights of alternatives are given in Table 12.
After determining evaluation results of alternatives under the criteria 5, by making a multiplication value related to each criterion with the obtained values, the weights were  According to the values, the wind power plant was determined as the best plant to be made an investment on. Hydropower plant follows it as rank number 2 and solar power plant follows it as rank number 3.

The conclusion
At the first stage of the application, using ranking approach for fuzzy numbers and the determined targets composing of a multi-objective mixed-integer programming model of the fossil power plants consumed most and renewable energy resources in Turkey, as a result comparisons made by decision makers, the priority values were calculated with Liou and Wang's total integral value method. In the research conducted, the target of energy demand has been the highest priority. Under the limitations composed using the data targeted in 2023, the model was solved with Lindo software. According to 2023 targets, the solar power plant is a plant which is required to be made an investment on according to the results obtained by Lindo software. The wind power plant followed it as rank number two.
In the other stage of the application, by taking an expert opinion working at the department of energy under the criteria of space requirement, provided employment, investment cost, environmental damage of power plants in Turkey, Saaty's fuzzy number approach was used. The pairwise comparisons for first criteria and for second power plants for each criterion were made. Later, according to Chang's grading analysis method, the priority values of alternatives were calculated. In the criteria, weights were obtained after calculating the priority vector. Investment cost having the highest ratio was the first criteria; as second, the environmental damage followed it. The wind power plant was determined as the most proper option. Hydropower plants and solar power plants followed it.
The study can be taken into account as an important resource for decision makers, but, in addition to these, a country's geopolitical position, inter-countries treaties, and financial status are considered for the decision of investment.