Location Selection for Wind Farms Using GIS Multi-Criteria Hybrid Model: An Approach Based on Fuzzy and Rough Numbers

: This paper presents spatial mathematical model in order to identify sites for the wind farms installment which can have significant support for the planners in the area of strategy and management of wind power use. The suggested model is based on combined use of Geographical Information Systems (GIS) with multi-criteria techniques of Best-Worst method (BWM) and MultiAttributive Ideal-Real Comparative Analysis (MAIRCA). Rough numbers and fuzzy logic are used to exploit uncertainty during data analysis in spatial mathematical model. The model is applied on the case study. Rough BWM model is used to determine weight coefficients of the criteria and rough MAIRCA method is used to rank separated sustainable locations. The implementation of MAIRCA method has shown that the location L3 is the most suitable for the wind farm in the area covered in the case study. Therefore, the suggested spatial mathematical model can be successfully used to identify the potential suitable sites for the wind farms in other areas with similar geographic conditions.


Introduction
The development of renewable energy sources (RES), together with measures aimed at more efficient use of energy, are priorities at a national and European level, as both the basis of environmental protection and of energy policy.According to Directive 2009/28/EZ [1] and the Decision of the Ministerial Council of the Energy Community [2] a binding target was determined that by 2020, renewable energy would make up 27% of the gross share of energy consumption in the Republic of Serbia.In order to achieve this level and to encourage use of RES, the Republic of Serbia has adopted a number of bylaws in the field of renewable energy which define the conditions for using renewable energy sources [3], and it has developed a National action plan for the use of sources of renewable energy [4].
States of the European Union (EU) have determined strategy that until 2020 increase the share of RES to 20% of total consumed energy aiming greater energetic independency [5].Apart of being inexhaustible these energy sources are the most acceptable according to the environmental protection standards and global climate changes.The use of wind energy has a low impact on the environment compared to many serious effects of conventional electricity production that contribute to climate changes and disrupt the natural balance [6][7][8][9][10][11]. Wind power is renewable and clean, does not pollute the air or emit CO2 and does not cause acid rains.In addition, it does not radiate and does not destroy the ozone layer [8][9][10].
The choice of suitable locations for wind farms is a complex issue that requires the careful and combined analysis of numerous criteria [12][13][14][15][16][17].On a global level, a large number of wind farms follow mainly economic criteria, ignoring any negative effects on the environment [18,19].The optimal selection of locations for wind farms requires the resolution of mutually conflicting factors of an economic and technological nature, with ecological and social limitations, while also respecting public opinion [20][21][22][23][24][25][26].In other words, spatial planners are faced with a double challenge, since they must answer for the wind farm projects they work on, which must be able to generate economic benefit on the one hand, while minimizing risks to the environment and the interests of stakeholders on the other [27][28][29][30][31][32][33][34][35][36][37][38][39].In order to achieve this, it is essential to establish certain rules that will evaluate different locations on the basis of a range of environmental, economic and social criteria and constraints [18,40].
This paper presents a new model for identifying the best locations for wind farms.The model considers 7 evaluation criteria.In the hybrid GIS-MCDA model, the rough Best-Worst method [41] is implemented to determine criteria weights.Weight coefficients of the criteria are obtained as the output values from the rough Best-Worst method (BWM), which are then used in the GIS for obtaining the final map of most suitable locations for wind farms.After determining the most suitable locations, the rough MAIRCA method is used to select the most suitable location from them.Application of the model is presented as a case study of the region of Vojvodina in Serbia.
The paper is organized in three sections.The introductory section presents the importance of renewable sources of energy.The second section of the paper presents and describes the phases of the proposed rough GIS-MCDA model.In the third and final section the testing of the rough GIS-MCDA model through empirical study in the region of Vojvodina is presented along with the suggestions for future research in this area.

GIS-Multi-Criteria Model Based on Rough Numbers
The location selection for the construction of wind farms represent a spatial problem that requires manipulating a large number of geophysical data, data from the environmental and socioeconomic data.Methodological hierarchical model in this paper is based on GIS-Multicriteria decision analysis structure (MCDA).Integrating GIS with the techniques for decision making creates a powerful tool for solving the problem of selecting optimal wind farm locations [42][43][44].
The rough approach was tested by means of a multi-criteria decision making (MCDM) model implemented in two phases: (1) determining the weight coefficients of the evaluation criteria; and (2) evaluation of the alternatives.To determine the weight coefficients of the criteria the BWM was modified using the rough approach.The BWM [41][42][43][44][45] is among the more recent methods.The primary advantages of the BW method suggested by the authors are as follows: (1) Compared with the Analytic Hierarchy Process (AHP) method, which in the literature is most commonly used to determine the weight coefficients [46,47], it requires significantly less pair wise comparison (the AHP method requires n(n − 1)/2 comparison, BWM 2n − 3 comparison); (2) the values of the weight coefficients obtained by the BWM are more reliable because comparison in the BWM is carried out with a higher consistency ratio compared with the AHP method; (3) while for the majority of MCDM models (e.g., AHP) the consistency ratio is a test of whether the comparison of criteria is consistent or not, in the BWM the consistency ratio is used to determine the level of confidence since the outputs from the BWM are always consistent; (4) the BWM for pair wise comparison of criteria uses only integers as opposed to other MCDM methods (e.g., AHP) which also require the use of fractional numbers.
After applying the BWM, in the second phase of the MCDM model for evaluating the alternatives is an original modification of the MAIRCA method [48][49][50] based on a rough approach.The authors chose the MAIRCA method because of its many advantages: (1) the mathematical framework of the method remains the same regardless of the number of alternatives and criteria; (2) it can be applied to greater numbers of alternatives and criteria; (3) it has a clearly defined rank of alternatives given in numerical values, enabling easier understanding of the results; (4) it is applicable for qualitative and quantitative criteria and (5) it gives stable solutions regardless of a change in the measuring scale for qualitative criteria or a change in the method of formulating quantitative criteria [48].
This approach uses capabilities of GIS in the management of geospatial data and flexibility of MCDA to combine factual information (e.g., land use, slope, communications, etc.) with a valuebased information (e.g., expert opinion, standards, surveys, etc.) [51].From a methodological point of view, the proposed GIS-MCDA model for selecting optimal locations for wind farms is comprising the following steps (Figure 1): After defining the problem and establishing the model, the constraints and evaluation criteria are identified.The constraints are based on criteria which limit (dismiss) possible alternatives and which are based on the Boolean relation (true/false), while the evaluation criteria can be quantified according to their degree of suitability for all feasible alternatives [50,51].
After this, individual assessment of the evaluation criteria and standardization of the criteria takes place.The criteria on the maps are presented in the form of GIS layers in different ways and in different forms.Use of the method of weighted linear combination (WLC) requires that all data sets be standardized [51] or transformed into units that can be compared.Use of the WLC method requires normalization of the weights.When determining the weight coefficients of the criteria, the rough BWM is used for calculating the normalized weight criteria and for the final application of the WLC method.The next section gives a detailed algorithm for the application of the GIS-MCDM model.

Rough Numbers
In group decision making problems, the priorities are defined on multi-expert's aggregated decision and process subjective evaluation of expert's decisions.Rough numbers consisting of upper, lower and boundary interval respectively, determine intervals of their evaluations without requiring additional information by relying only on original data [12,[51][52][53][54].Hence, obtained expert decision makers (DMs) perceptions objectively present and improve their decision making process.According to [53], the definition of rough number is shown below.
Let's U be a universe containing all objects and X be a random object from U .Then we assume that there exists set build with k classes representing DMs preferences, Apr J , upper approximation ( ) q Apr J and boundary interval ( ) q Bnd J are determined, respectively, as follows [55]: { } ( ) / ( ) The object can be presented with rough number (RN) defined with lower limit ( ) q Lim J and upper limit ( ) q Lim J , respectively: where L M and U M represent the sum of objects contained in the lower and upper object approximation of q J , respectively.For object q J , rough boundary interval ( ) IRBnd J presents interval between lower and upper limit as: Rough boundary interval presents measure of uncertainty.The bigger ( ) q IRBnd J value shows that variations in experts' preferences exist, while smaller values show that experts had harmonized opinions without major deviations.In ( ) are comprised all objects between lower limit ( ) q Lim J and upper limit ( ) q Lim J of rough.

Number
( ) q RN J .That means that ( ) q RN J can be presented using ( ) q Lim J and ( ) q Lim J .
( ) ( ), ( ) Since rough numbers belong to the group of interval numbers, arithmetic operations applied in interval numbers is also appropriate for rough numbers [55].
Since rough numbers belong to the group of interval numbers, arithmetic operations applied in interval numbers is also appropriate for rough numbers [15].If A and B presents two rough numbers  , k denotes constant, 0 k > , then the arithmetic operations with ( ) RN A , ( ) RN B and k are as follows: (1) Addition of rough numbers "+" (2) Subtraction of rough numbers "−" (3) Multiplication of rough numbers "×" (4) Dividing of rough numbers "/" (5) Scalar multiplication of rough numbers, where Ranking rule of rough numbers.( ) Lim A , ( ) Lim B represent their lower and upper limits, respectively, are ranked by the use of the following rules [52,53]: If the rough boundary interval of a rough number is not strictly bound by another, then the ranking order is easily determined, i.e., (a) If If the rough boundary interval of a rough number is strictly bound by another, then ranking becomes awkward and medians

Rough Based Best-Worst Method (R-BWM)
In order to take into account the subjectivity that appears in group decision making more comprehensively, in this study a modification of the Best-Worst method (BWM) was carried out using rough numbers (RN).The application of RN eliminates the necessity for additional information when determining uncertain intervals of numbers.In this way the quality of the existing data is retained in group decision making and the perception of experts is expressed in an objective way in aggregated Best-to-Others (BO) and Others-to-Worst (OW) matrices.Since the method is very recent, the literature so far only has the traditional (crisp) BWM [44,45,56] and modification of the BWM carried out using fuzzy numbers [57,58].The approach in this section introduces RN which enable more objective expert evaluation of criteria in a subjective environment.The proposed modification of the BWM using RN (R-BWM) makes it possible to take into account doubts that occur during the expert evaluation of criteria.R-BWM makes it possible to bridge the existing gap in the BWM methodology with the application of a novel approach in the treatment of uncertainty based on RN.The following section presents the algorithm for the R-BWM that includes the following steps: Step 1. Determining the set of evaluation criteria.This starts from the assumption that the process of decision making involves m experts.In this step, experts consider the set of evaluation criteria and select the final set of criteria , where n represents the total number of criteria.
Step 2. Determining the most significant (most influential) and worst (least significant) criteria.
The experts decide on the best and the worst criteria from the set of criteria . If the experts decide on two or more criteria as the best, or worst, the best and worst criteria are selected arbitrarily.
Step 3. Determining the preferences of the most significant (most influential) criteria (B) from set C over the remaining criteria from the defined set.Under the assumption that there are m experts and n criteria under consideration, each expert should determine the degree of influence of the best criterion B on the criteria j ( ).This is how we obtain a comparison between the best criterion and the other criteria.The preference of criterion B compared to the j-th criterion defined by the e-th expert is denoted with .As a result a Best-to-Others (BO) vector is obtained: .As a result an Others-to-Worst (OW) vector is obtained: where e jW a represents the influence (preference) of criterion j in relation to criterion W, whereby 1 e WW a = .This is how we obtain OW matrices Step 5. Determining the rough BO matrix for the average answers of the experts.Based on the BO matrices of the experts' answers ; where { } where , ,..., Step 6. Determining the rough OW matrix of average expert responses.Based on the WO matrices of the expert responses  , as with the rough BO matrices, for each element e jW a we form matrices of the aggregated sequences of the experts *e where { }  , ,..., Step 7. Calculation of the optimal rough values of the weight coefficients of the criteria . On this basis we can conclude that in the case of the rough of the weight coefficients of the criteria the condition is met that The previously defined limits will be presented in the following min-max model: where is the rough weight coefficient of a criterion.
Model (22) is equivalent to the following model:  ) , respectively represent the values from the average rough OW and rough BO matrices (see Equations ( 17) and ( 20)).By solving model (23) we obtain the optimal values of the weight coefficients of the evaluation criteria

The Consistency Ratio of the R-BWM
The consistency ratio is a very important indicator by means of which we check the consistency of the pair wise comparison of the criteria in the rough BO and rough OW matrices.respectively represent the preference of the best criterion over criterion j, the preference of criterion j over the worst criterion, and the preference of the best criterion over the worst criterion.
However, when comparing the criteria it can happen that some pairs of criteria j are not completely consistent.Therefore the next section defines the consistency ratio (CR), which gives us information on the consistency of the comparison between the rough BO and rough OW matrices.In order to show how CR is determined we start from a calculation of the minimum consistency when comparing the criteria, which is explained in the following section.
As previously indicated, pair wise comparison of the criteria is carried out based on a predefined scale in which the highest value is 9 or any other maximum from a scale defined by the decision maker.The consistency of the comparison decreases when , which continues to affect the value of ξ .Based on these relationships we can conclude that As the largest inequality occurs when .Thus, we obtain Equation (25) [ ] Since for the minimum consistency Since we are using rough numbers, and if there is no consensus between the DM on their preferences of the best criterion over the worst criterion, then  , Table 1.1,2,...,9 which is defined in [45].On the basis of CI (Table 1) we obtain the consistency ratio (CR) The CR takes values from interval [ ] 0,1 , where values closer to zero show high consistency, while the values of CR closer to one show low consistency.

Rough MAIRCA Method
The basic assumption of MAIRCA method is to determine the gap between ideal and empirical weights.The summation of gaps for each criterion gives the total gap for every observed alternative.Finally, alternatives will be ranked, and the best ranked alternative is the one with the smallest value of the total gap.MAIRCA method shall be carried out in 6 steps [12,48,49]: Step 1.Formation of the initial decision matrix ( Y ).The first step includes evaluation of l alternatives per n criteria.Based on response matrices Yk = [y k ij]l×n by all m experts we obtain matrix * Y of aggregated sequences of experts where { } where e denotes e-th expert ( In such a way, rough vectors where l denotes number of alternatives, and n denotes total sum of criteria. Step 2. Define preferences according to selection of alternatives i A P .When selecting alternative, a decision maker (DM) is neutral, i.e., does not have preferences to any of the proposed alternatives.Since any alternative can be chosen with equal probability, preference per selection of one of l possible alternatives is as follows where l denotes number of alternatives. Step where ndenotes number of criteria, where ( ) pij RN t denotes elements of theoretical assessment matrix, and . Normalization of the mean initial decision matrix ( 31) is done by applying Equations ( 37) and ( 38) (a) For the "benefit" type criteria (larger criterion value is preferable)

y y y y RN y Lim y Lim y y y y y y y
(b) For the "cost" type criteria (lower criterion value is preferable) Lim y Lim y y y y y y y where i y − and i y + denote minimum and maximum values of the marked criterion by its alternatives, respectively: Step 5. Calculation of the total gap matrix ( G ). Elements of G matrix are obtained as difference (gap) between theoretical ( pij t ) and real evaluations ( rij t ), or by actually subtracting the elements of the theoretical evaluation matrix ( p T ) with the elements of the real evaluation matrix ( r T ) where n denotes number of criteria, l denotes number of alternatives, and gij represents the obtained gap of alternative i as per criterion j.Gap gij takes values from the interval rough number according to Equation ( 42) It is preferable that

RN g RN t ≈
. This means that alternative i A for criterion i C is the worst (anti-ideal) alternative.
Step 6. Calculation of final values of the criteria functions ( i Q ) per alternatives.Values of criteria functions are obtained by summing the gaps from matrix (41) for each alternative as per evaluation criteria, i.e., by summing matrix elements ( G ) per columns as shown in Equation ( 43) where n denotes number of criteria, m denotes number of the chosen alternatives.
Ranking of alternatives can be done by applying rules governing ranking of rough numbers described in Section 2.1.

Aggregation of the Criteria
After determining the weight coefficients of the clusters/criteria, it is necessary to carry out aggregation of the criteria.WLC is used in the process of aggregating the criteria.The WLC aggregation method multiplies each standardized factor map (i.e., each raster cell within each map) by its factor weight and then sums the results.The following mathematical expression was used to combine the evaluation criteria (factors) according to the WLC method: where S is suitability, wi is the normalized value of the weight of factor i, and xi is the criterion score of factor i.
All GIS software systems provide the basic tools for evaluation of such a model [51].WLC is relatively easy to understand and can be applied in many different situations.In addition, it is compensatory, meaning that low scores in one criterion can be compensated for by high scores in another one, which is desired for this particular decision problem.For these reasons, WLC was selected as the method of aggregation.As the final result of applying WLC and the aggregation of criteria maps, a suitability map is obtained.
In the next step it is necessary to separate the suitable locations, by separating the cells with the highest values from the suitability map, which represent alternatives for the location of wind farms.The cells are filtered by combining the application of the GIS arithmetic operations and questions, by means of which feasible locations for installing wind farms are defined.

Selected Area Case Study
Based on the research of energy resources of wind and solar radiation, conducted by the Institute of Multidisciplinary Research (IMSI), which were published in the "Atlas of the energy potential of the sun and wind in Serbia" [59], it is recognized that in Serbia exist potentially suitable locations for the wind energy utilization.The results of study "Wind atlas of Vojvodina" [60] have shown that micro region of the southern Banat has the greatest wind energy potential in Serbia, with average annual rates of more than 6 m/s, is considered favorable for the realization of a wind farm project from the economic aspect.
South Banat is a micro region that is located in the southeast of the autonomous province of Vojvodina.It spreads between 44°39′ and 45°0′ N, 20°25′ and 21°34′ E, with area of 4252 km 2 (Figure 2).It is an integral part of the Pannonian part of Serbia, with no relief barriers.Located at the entrance to the Djerdap gorge from which the east wind blows strongest in Serbia, so called košava.For the full utilization of the energy potential of Košava wind crucial question is the choice of location of wind turbines.An additional benefit of installing and maintaining wind farms is developed network of roads and railways in the region.This study was designed in order to submit a reliable spatial multicriteria model for the selection of suitable sites for the installation of wind farms in the area of micro region of South Banat

Identification of Criteria that Influence the Wind Farms Location Selection
Identification of the criteria for the wind farms location selection is an important step in the proposed model [61].Based on the geographical features of the Southern Banat, analysis of legislative acts in the field of environmental protection, protection of forests, water protection, life safety standards, review the recent literature and considering the opinion of decision-makers, who dealt with similar problems, the following criteria (Table 1) were adopted [18]:

C1
Wind speed average.The most important economic factor to consider when building a wind power plant on the wind speed at the monitored site.Locations with an average wind speed of more than 3 m/s are evaluated continuously, with the recommendation of the selection of sites to the average wind speed of 20 m/s.At speeds over 20 m/s it is possible that wind will damage the equipment.

C2
Land cover/use.Land is one of the most important factors for energy investments.
Wind farms should be installed in the area in which they negligibly interfere with existing land use outside protected areas, artificial surfaces, wetlands, aquatic and forest areas.The provided information on the land use is completed using CORINE project that provides sufficient comprehensiveness and consistency in European states.

C3
Distance from main communication.The specificity of the installation equipment transport and subsequent maintenance of the wind farm is reflected in the need for secure road network to the location of wind turbines, with maximum use of the existing road network.In order to reduce the cost of transport is desirable that the distance from the roads outside the safety zone of 200 m, is as small as possible.

C4
Slope.In a physical sense the wind farm should be located at suitable locations.Areas with steep slopes are generally considered less suitable, because of the extra costs of construction and maintenance.

C5
Orientation Aspect.Orientation of slopes relative to the direction of the wind is an important criterion from the point of full utilization of the wind potential.It is desirable that the orientation of the slope on the location of the wind farm is toward the dominant wind blow.

C6
Distance from urban places.Wind farms near urban areas can cause negative impacts on the local population.In order to minimize negative impacts on the environment it is recommended that wind farms are at the greatest possible distance from urban areas, at least 500 m.

C7
Distance from power lines.One of the most important issues when it comes to reducing the cost of building wind farms is reducing costs for the construction of energy transmission infrastructure.When choosing a location, it is necessary to consider the availability of and access to existing electrical power lines.Due to the risk of potential accidents on power lines, the minimum distance of 100 m should be provided with adequate space between the columns and cable of power lines.

Data Collection and Entry into GIS
Entering GIS, each of the criteria is presented in the form of spatially defined maps, which represent the attributes of cells which have different ranges or scale value.The data used in this study have been compiled from various sources.All the processes of transformation and GIS data modeling were derived using integrated tools ESRI ArcGIS 10.2 software.Map of land use is obtained on the basis of CORINE Land Cover 2006 (CLC2006) database gathered through the CORINE (Coordination of information on the environment program.Data on wind speed are shown using georeferencing method based on a study "Wind atlas of Serbia".Further information is downloaded in digital format from the Digital Terrain Model DMT (25 × 25 m.) and Digital topographic map 1:300,000made in Military Geographical Institute in Belgrade.

GIS-MCDA
Given that the data on the criteria were collected in different ways and have different formats, it is necessary that data sets are standardized and converted into units that can be compared with each other.There are a number of approaches that can be used for attribute layers criteria to be comparable.In this case, with the criteria set whose elements are categorical values (Land use), the discrete classification is used, in which experts gave out value of the elements of fuzzy sets directly.In all other categories, which are the values of the gradual change from one location to another, elements of the set are standardized using fuzzy concept based on linear or sigmoidal membership function.For fuzzification scale ranges from 0 to 1 byte was used, where closer to 0 is the least, and closer to 1 the most appropriate value set in relation to the suitability of the location of the wind farm.Values with 0 are excluded from the analysis.Standardized criteria for evaluating the fuzzy membership functions and forms are shown in Table 2 and Figure 3.
It should be noted that the adopted criteria are not final and can be extended to a larger number of objective circumstances, taking into account application areas.
In this case study, when determining the weights of criteria, ten experts were involved with rich experience in the design of renewable energy sources and capital projects.Experts who have an experience of minimum five years in providing consulting services in the design of renewable energy sources were identified.After interviewing experts, the collected data are processed and aggregation of their opinions was executed.
In the first step of applying BWM experts used the following scale for the evaluation criteria: 1-No impact; 2-Very little impact; ...; 9-Very large impact.Seven experts took part in the research.The R-BWM was used to determine the weight coefficients of the criteria.After defining the criteria for evaluation, the experts also determined the best (B) and worst (W) criterion.On this basis, the experts determined the BO and OW matrices in which the preferences of the B and W over the criteria were considered for the remaining criteria from the defined set.Evaluation of the criteria was carried out using a scale [ ] 1,9 e ij a ∈ [44].The BO and OW matrices are presented in Table 3.

8.92
U BW a = on the basis of the aggregated decisions of the experts, and they affect the interval of the RN, it is not possible to predefine the values of the consistency index ξ .Using Equation ( 27), the values of the consistency index (

The Summing of Pondered Maps into Suitability Map
WLC method integrated into ArcGIS Advanced 10.2 software and ArcGIS Spatial Analyst extension were used in this step to sum normalized weighted criteria values.WLC aggregation of map criteria into the final suitability map of suitable location for wind farms is presented in the same value range from 0 to 1 was completed on the basis of the adopted criteria and particular weights determined by experts (Figure 4a).The greater cell values represent the area more suitable for the wind farm location.
The final step is to extract the cells with the highest values in the final suitability map.These are the areas which pixels have a value of 0.750 to 1 representing the most favorable locations for wind farms.There is 89.7 km 2 of very suitable area for the construction of wind farms, obtained from the analysis in the Southern Banat, differentiated in eight locations (Figure 4b).These sites can be used as the main locations for the construction of wind farms in the Southern Banat.

The Location Ranking
Taking into account that the main purpose of this analysis is identification of the most suitable wind farm locations, this step involves ranking of the locations based on criteria parameters calculated using Zonal statistical algorithm ArcGis 10.2 software for each location (Table 8).After identification of suitable locations, their ranking and selection of the best location using MAIRCA method is conducted.Locations evaluation and ranking were obtained using Expressions ( 29)-( 43), Table 9. Ranking of alternatives using MAIRCA method was carried out on the basis of the total value of the gap, whereby as the best alternative is taken one that has the lowest value of the gap, or the one that is closest to the ideal alternative.An alternative that is the worst is the one that has the highest value of the gap, or the one that is farthest from the ideal alternative.It is desirable that alternative has the lowest value of the total gap.Based on the conditions for ranking locations by using MAIRCA method it can be concluded that the best ranked alternative is L3.

Conclusions
This paper proposed the combination of rough sets in combination with multi-criteria methods within GIS software environment for the determination of suitable locations for the construction of wind farms.Developed rough BWM model that was used in this study is used to determine the weight or importance of each criterion.The WLC is used to sum up the weights and to identify the final suitable locations for wind farms.Application of GIS multi-criteria model in choosing locations for wind farms is carried out on a real case of micro region of the Southern Banat, Vojvodina, Serbia.The final suitability map is obtained defining eight locations that are most suitable for the construction of wind farms.In the following procedure, using rough MAIRCA model, the locations were ranked and the best one is proved to be location L3.
It should be pointed out that this procedure allows the inclusion of other criteria that were not included in this study.This model extends the theoretical framework of knowledge in location selection.The current problem is discussed using the new methodology, which creates the basis for further theoretical and practical upgrade.

Figure 1 .
Figure 1.Model for wind farm location selection.

a 1 B A , 2 BStep 4 .
represents the influence (preference) of the best criterion B over criterion j, whereby 1 e BB a = .This is how we obtain BO matrices Determining the preferences of the criteria from set C over the worst criterion (W) from the defined set.Each expert should determine the degree of influence of criterion j ( 1, 2,..., j n = ) in relation to criterion W. The preference of criterion j in relation to criterion W defined by the e-th expert is denoted as e represents the e-th expert ( sequences.We thus obtain the averaged rough BO matrix of average responses B A where e represents the e-th expert ( the rough sequences.Thus, we obtain the averaged rough OW matrix of average responses W A

.
In this way the condition is met that the weight coefficients are found at interval [0,1], ( 1, 2,..., ) represents the optimum values of the weight coefficients,( ) [ , ]represent the weight coefficients of the best and worst criterion respectively, while ( values, then we need to subtract the value ξ from ( )

a
the maximum possible values of ξ , which is the CI for the R-BW method.Since we obtain the values of ( ) on the basis of the aggregated decisions of the DM, and these change the IVFRN interval, it is not possible to predefine the values of ξ .The values of ξ depend on uncertainties in the decisions, since uncertainties change the RN interval.As explained in the algorithm for the R-BW method, interval , on uncertainties in evaluating the criteria.If the DM agree on their preference for the best criterion over the worst then BW a represents the crisp value of BW a from the defined scale and then the maximum values of ξ apply for different values of { } 1, 2,...,9 BW a ∈

Step 4 .
Determination of real evaluation ( r T ).Calculation of the real evaluation matrix elements ( r T ) is done by multiplying the real evaluation matrix elements ( p T ) and elements of the initial decision making matrix ( X ) according to the following equation:

Figure 2 .
Figure 2. The geographical position of study area.

Figure 3 .
Figure 3. Standardized criteria maps with fuzzy membership sets.

Figure 4 .
Figure 4. Final suitability map (a) and most suitable locations (b) for the wind farm in the Southern Banat, Vojvodina, Serbia.

Table 1 .
Values of the consistency index (CI).
CI ( max ξ ) 0.00 0.44 1.00 1.63 2.30 3.00 3.73 4.47 5.23In Table1the values BW a are taken from the scale { } , where m denotes number of experts.Therefore, for the group of rough matrices Y 1 , Y 2 , …,Y m we obtain rough sequences (7)ote sequences for describing of relative importance of criterion i in relation to alternative j.By applying Equations (1) through(7), sequences ij m RN y .Consequently, rough matrices Y 1L , Y 2L , …, Y mL will be obtained for the

Table 2 .
Fuzzy standardization of criteria.

Table 3 .
The BO and OW expert evaluation matrices.

Table 6 .
Aggregated rough BO and rough OW matrices.

Table 7 .
Optimal values of the criteria.

Table 8 .
Initial matrix for location evaluation.