Fuzzy Complex Proportional Assessment Applied in Location Selection for Installation of Photovoltaic Plants

: This article proposes the application of a multi-criteria decision-making method based on fuzzy complex proportional assessment (COPRAS-F), applying business sustainability criteria, to rank alternative sites for the installation of photovoltaic plants, as selected from a geographic suitability index map of vast and heterogeneous regions in Southern Brazil, obtained from a geographic information system. As a key feature, the methodology serves the interests of stakeholders during decision-making meetings in a personalized manner, eliminating the need for consensus to overcome conﬂicts and antagonistic positions, allowing them to maintain a strictly professional view and to communicate in business language. The results indicate that the best location from a geographical viewpoint may not be the best location from a business viewpoint.


Introduction
The increasing demand for electric power in a global economy, characterized by an alarming recession of global fossil fuel reserves, has imperatively induced a collapse in energy generation balance [1]. Therefore, generation expansion planning for renewable and unrenewable energy has been used to define the least costly capacity expansion plan to meet forecasted demand via a pre-defined reliability criterion and emission constraint through a long-term plan [2].
Photovoltaic energy stands out among renewable energy sources for being environmentally and socioeconomically beneficial, if its implementation takes place in appropriate locations through efficient planning and management [3]. It has generally been met with public acceptance.
Currently, making photovoltaic energy system projects more widely and swiftly feasible depends on adequate incentive policies or innovative business models that are capable of mitigating high initial costs, as well as lowering uncertainty regarding long-term financial return [3].
These factors add complexity to the high-risk decision-making processes of locating possible installation sites for photovoltaic plants, and of implementing models of geographic information systems (GIS), while applying geographic criteria and restrictions. In the face of such complexity, policymakers, specialists, and researchers in the field have proposed the use of models based on multi-criteria decision-making (MCDM) methods that include business-oriented criteria.

Literature Review
From the publications of "State of the Art", it is possible to notice that several articles, which focus on energy planning-applied GIS and SIG-MCDM methods with geographic subcriteria in the selection and rank of alternative sites for the location of photovoltaic power plants, feature weaknesses when not jointly applying business criteria.
By further analyzing these methods, attention is drawn to the conclusions of the following researches, which assessed photovoltaic energy projects: • Tomosk et al. [4] demonstrated that only the use of geographical criteria for the location of photovoltaic plants is insufficient. The authors concluded that financial return rate maps differ substantially from solar irradiation maps in the United States. Thus, the prediction that climate is always directly related to the rate of return of a solar project is mistaken. The reason for that is because, in this case, the rate of economic return is directly related to the tax model and time of use; • Vafaeipour et al. [5] stated, "that in contrast to common assumptions, the identification of suitable geographical locations for the implementation of large photovoltaic plants is not only related to the amount of solar radiation received, but there are also many other economic, environmental, social, and risk factors that must be taken into account".
Changes are necessary to overcome these deficiencies in the GIS and SIG-MCDM decision-making methods, which only take into account the geographic subcriteria. The proposed solution is the application of MCDM methods with business subcriteria to rank the alternatives of sites as selected by the geographic suitability index, allowing the use of the investors' everyday business language in choosing between more geographically compatible locations for photovoltaic plant installation.
There has been an evolutionary process in the lines of research involving decision making in renewable energy projects. This process began with the application of overlay operations in geographic information systems (GIS) of geographic subcriteria fuzzified maps (climate, geomorphology, location, and land use), with the elimination of geographical restrictions, so that a geographic suitability index map of alternative areas for the installation of renewable energy plants ( Figure 1) could be generated. The works listed below (Table 1) stand among those focusing on photovoltaic plants: Table 1. GIS.

Authors Results
Brewer et al. [6] Generated maps of the solar energy potential of southwestern USA, including social acceptance data.
Janke [7] Identified 191 km 2 of highly suitable areas for solar farms in Northwest Colorado, in the USA.
Sabo et al. [8] Identified 10,092.08 km 2 for the installation of large photovoltaic plants in Malaysia.
Wang et al. [9] Applied GIS to select 4005 sites suitable for photovoltaic energy projects, with the majority in the Shigatse and Ngari regions in Tibet, China.
Maleki et al. [10] Applied GIS to assess the size and location of rural areas for photovoltaic projects, considering energy quality and cost for the Birjand region in Iran.
Another concomitant research has introduced the GIS-MCDM method by applying weighted linear combination (WLC) or ordered weighted average (OWA) when adding fuzzy maps of geographic subcriteria during the GIS overlap process ( Figure 2) to identify alternative photovoltaic plant locations. The works listed below (Table 2) stand among those focusing on photovoltaic plants: Table 2. Typical GIS-MCDM.

Charabi and
Gastli [11] Applied the MCDM ordered weighted average fuzzy (FLOWA) to define weights used in GIS when generating a map of Oman that showed high suitability for the installation of large photovoltaic plants in 0.5% of the country's area Al Garni and Awasthi [12] Applied the MCDM analytic hierarchy process (AHP) to define weights used in GIS-MCDM (WLC) when generating a suitability map of Saudi Arabia that showed that 16% of the country is suitable for installing small-scale photovoltaic plants Aydin et al. [13] Applied GIS-MCDM (OWA) to generate maps of global performance indexes for environmental objectives and economic viability, indicating the provinces of Aydin, Usak, Burdur, Denizli, and Mugla in Turkey as priorities for photovoltaic energy sites.

Wanderer and
Herle [14] Applied MCDM (AHP) to define the weights used in GIS-MCDM (WLC and OWA) methods, generating a map of suitability for Andalusia, in Spain, regarding solar energy, presenting more sites that are suitable-located in the Guadalquivir Valley.
Watson et al. [15] Applied GIS-MCDM (AHP) to generate a map of the south of England, which showed an added 345 km 2 of suitable areas.
Merrouni et al. [16] Applied MCDM (AHP) to define the weights used to aggregate the suitability index of categorized subcriteria maps, when determining that eastern Morocco is highly suitable, with 19% of the land proven as appropriate for the installation of large photovoltaic plants.
The same period also showed the increased application of MCDM methods with technical-economic and socio-environmental business criteria (Figure 3), more specifically in energy planning. The following (Table 3) are highlights among these works: Table 3. MCDM.

Mourmouris
and Potolias [17] Applied the REGIME method with business criteria to classify renewable energy systems based on energy efficiency for installation on the island of Thassos, in Greece Vafaeipour [5] Applied step-wise weight assessment ratio analysis-weighted aggregates sum product assessment (SWARA-WASPAS) with business criteria in prioritizing city alternatives for photovoltaic plant location in Iran Wang et al. [18] Applied MCDM methods fuzzy analytic hierarchy process-technique for order of preference by similarity to ideal solution (FAHP-TOPSIS) to rank locations for the construction of a solar power plant in Vietnam.
Barin et al. [19] Applied a novel fuzzy-based expert system with business criteria for selecting renewable energy technologies, along with technical-economic and environmental criteria in Brazil.
Karasan et al. [20] Applied the combinative distance based assessment (CODAS) method, crisp and ordinary fuzzy CODAS methods, with business criteria to select sites for the construction of a wind farm in Turkey.

Samanlioglu
and Ayag [21] Applied MCDM Fuzzy AHP (F-AHP) and fuzzy preference ranking organization method enrichment evaluation (F-PROMETHEE) methods to evaluate alternative solar plant locations in Turkey.
It has been observed that the GIS and GIS-MCDM methods that apply geographic subcriteria, and MCDM that apply business subcriteria, are complementary, because: • The GIS and GIS-MCDM methods present strengths when selecting and classifying alternative areas from a geographic point of view, but they present weaknesses when technical-economic and socio-environmental characteristics are not considered; • The MCDM methods present strengths when selecting and classifying alternative areas from a business point of view, but they present weaknesses when climatic, geomorphological, location and land use characteristics are not considered.
In order to eliminate the weaknesses of the previous methods, this article suggests the following model ( Figure 4) for the selection and classification of alternative sites for the installation of photovoltaic plants in vast areas of land, generally characterized by their great heterogeneity of the climate, geomorphology, location, land use, technical-economic, and socio-environmental conditions: • Selection of alternative sites from a geographic suitability index map developed in current GIS-MCDM from geographic subcriteria; • Application of the MCDM method based in business criteria to rank specific alternatives of sites for the installation of photovoltaic plants.

The Model
The model proposed in this article was developed from the identification and description of the following elements:

•
The problem: consists in the selection and classification of alternatives of suitable areas for the installation of photovoltaic plants, localized in different regions of a large territory, presenting climatic, geomorphological, location, land use, economic, social, and environmental heterogeneities; • The alternatives: are selected from the geographic suitability index map of the area under evaluation, obtained by overlaying fuzzified geographical subcriteria maps; • The stakeholders: the choice of stakeholders must take into account the complexity of the decision-making process and must come from different groups of stakeholders with different scenarios of interests, such as public and private investors, different spheres of public administration, environmentalists, and communities; • The business criteria and subcriteria: stakeholders are responsible for evaluating, based on the literature, which economic, technical and environmental criteria of the business with their respective subcriteria are in the interest of business, considering the availability of public or private databases, for research on the values of the business subcriteria; • The MCDM method: the application of MCDM methods using business criteria is widely used in energy planning [22]. Among other advantages, it allows stakeholders to reach an agreement based on the evaluation of numerous conflicting criteria from multiple dimensions. In the decision-making process for locating decentralized energy projects, typical for several renewable energies, some factors significantly influence the choice of the multi-criteria method: Energies 2020, 13, x; doi: FOR PEER REVIEW 5 of 19 • The problem: consists in the selection and classification of alternatives of suitable areas for the installation of photovoltaic plants, localized in different regions of a large territory, presenting climatic, geomorphological, location, land use, economic, social, and environmental heterogeneities; • The alternatives: are selected from the geographic suitability index map of the area under evaluation, obtained by overlaying fuzzified geographical subcriteria maps; • The stakeholders: the choice of stakeholders must take into account the complexity of the decisionmaking process and must come from different groups of stakeholders with different scenarios of interests, such as public and private investors, different spheres of public administration, environmentalists, and communities; • The business criteria and subcriteria: stakeholders are responsible for evaluating, based on the literature, which economic, technical and environmental criteria of the business with their respective subcriteria are in the interest of business, considering the availability of public or private databases, for research on the values of the business subcriteria; • The MCDM method: the application of MCDM methods using business criteria is widely used in energy planning [22]. Among other advantages, it allows stakeholders to reach an agreement based on the evaluation of numerous conflicting criteria from multiple dimensions. In the decision-making process for locating decentralized energy projects, typical for several renewable energies, some factors significantly influence the choice of the multi-criteria method:  The context of the decision-making process of companies investing in the energy sector, in which there are usually multiple decision makers or investors, thus hindering the reaching of a consensus and leading decision makers to act more individually due to economic issues;  Uncertainty plays an important role in the project's decision-making problems that involve incomplete and imperfect information of real-world situations [23];  Conflicts among decision makers and investors regarding the real importance of interests from multiple dimensions, represented in practice by those criteria with accompanying business-related subcriteria;  The uncertainties surrounding the real values of the multiple interest dimensions regarding alternative sites for photovoltaic plants.
The choice of the MCDM method for the installation of photovoltaic plants in vast and large regions with geo-economic and socio-environmental heterogeneities must take into account all the aspects that were mentioned above.

Identification of the Problem
The area under evaluation for selecting sites alternatives for installing photovoltaic plants spreads throughout 199,315 km 2 in Southern Brazil ( Figure 5). It consists of 10 mesoregions with a population (estimated by the Brazilian Institute of Geography and Statistics-IBGE) of 11,348,937 inhabitants distributed across 399 municipalities, 335 of which with populations under 30,000. The area comprises a large-scale energy distribution network, mostly with 34.5/13.8 kV substations, representing an ideal framework for small-scale distribution of photovoltaic power generation.

Identification of Alternative Sites Based on a Geographic Suitability Index Map
The geographic suitability index map of the area under evaluation ( Figure 5) was prepared in Quantum GIS (QGIS) by the Department of Cartographic Engineering of the Federal University of Paraná (UFPR). This was prepared through weighted aggregation, applying the MCDM weighted linear combination (WLC) in the overlay of geographic subcriteria suitability index maps, which were obtained from the fuzzification of the irradiation, temperature, land slope, land aspect, and distance The context of the decision-making process of companies investing in the energy sector, in which there are usually multiple decision makers or investors, thus hindering the reaching of a consensus and leading decision makers to act more individually due to economic issues; Energies 2020, 13, x; doi: FOR PEER REVIEW 5 of 19 • The problem: consists in the selection and classification of alternatives of suitable areas for the installation of photovoltaic plants, localized in different regions of a large territory, presenting climatic, geomorphological, location, land use, economic, social, and environmental heterogeneities; • The alternatives: are selected from the geographic suitability index map of the area under evaluation, obtained by overlaying fuzzified geographical subcriteria maps; • The stakeholders: the choice of stakeholders must take into account the complexity of the decisionmaking process and must come from different groups of stakeholders with different scenarios of interests, such as public and private investors, different spheres of public administration, environmentalists, and communities; • The business criteria and subcriteria: stakeholders are responsible for evaluating, based on the literature, which economic, technical and environmental criteria of the business with their respective subcriteria are in the interest of business, considering the availability of public or private databases, for research on the values of the business subcriteria; • The MCDM method: the application of MCDM methods using business criteria is widely used in energy planning [22]. Among other advantages, it allows stakeholders to reach an agreement based on the evaluation of numerous conflicting criteria from multiple dimensions. In the decision-making process for locating decentralized energy projects, typical for several renewable energies, some factors significantly influence the choice of the multi-criteria method:  The context of the decision-making process of companies investing in the energy sector, in which there are usually multiple decision makers or investors, thus hindering the reaching of a consensus and leading decision makers to act more individually due to economic issues;  Uncertainty plays an important role in the project's decision-making problems that involve incomplete and imperfect information of real-world situations [23];  Conflicts among decision makers and investors regarding the real importance of interests from multiple dimensions, represented in practice by those criteria with accompanying business-related subcriteria;  The uncertainties surrounding the real values of the multiple interest dimensions regarding alternative sites for photovoltaic plants.
The choice of the MCDM method for the installation of photovoltaic plants in vast and large regions with geo-economic and socio-environmental heterogeneities must take into account all the aspects that were mentioned above.

Identification of the Problem
The area under evaluation for selecting sites alternatives for installing photovoltaic plants spreads throughout 199,315 km 2 in Southern Brazil ( Figure 5). It consists of 10 mesoregions with a population (estimated by the Brazilian Institute of Geography and Statistics-IBGE) of 11,348,937 inhabitants distributed across 399 municipalities, 335 of which with populations under 30,000. The area comprises a large-scale energy distribution network, mostly with 34.5/13.8 kV substations, representing an ideal framework for small-scale distribution of photovoltaic power generation.

Identification of Alternative Sites Based on a Geographic Suitability Index Map
The geographic suitability index map of the area under evaluation ( Figure 5) was prepared in Quantum GIS (QGIS) by the Department of Cartographic Engineering of the Federal University of Paraná (UFPR). This was prepared through weighted aggregation, applying the MCDM weighted linear combination (WLC) in the overlay of geographic subcriteria suitability index maps, which were Uncertainty plays an important role in the project's decision-making problems that involve incomplete and imperfect information of real-world situations [23]; The MCDM method: the application of MCDM methods using business criteria is widely used in energy planning [22]. Among other advantages, it allows stakeholders to reach an agreement based on the evaluation of numerous conflicting criteria from multiple dimensions. In the decision-making process for locating decentralized energy projects, typical for several renewable energies, some factors significantly influence the choice of the multi-criteria method:  The context of the decision-making process of companies investing in the energy sector, in which there are usually multiple decision makers or investors, thus hindering the reaching of a consensus and leading decision makers to act more individually due to economic issues;  Uncertainty plays an important role in the project's decision-making problems that involve incomplete and imperfect information of real-world situations [23];  Conflicts among decision makers and investors regarding the real importance of interests from multiple dimensions, represented in practice by those criteria with accompanying business-related subcriteria;  The uncertainties surrounding the real values of the multiple interest dimensions regarding alternative sites for photovoltaic plants.
The choice of the MCDM method for the installation of photovoltaic plants in vast and large regions with geo-economic and socio-environmental heterogeneities must take into account all the aspects that were mentioned above.

Identification of the Problem
The area under evaluation for selecting sites alternatives for installing photovoltaic plants spreads throughout 199,315 km 2 in Southern Brazil ( Figure 5). It consists of 10 mesoregions with a population (estimated by the Brazilian Institute of Geography and Statistics-IBGE) of 11,348,937 inhabitants distributed across 399 municipalities, 335 of which with populations under 30,000. The area comprises a large-scale energy distribution network, mostly with 34.5/13.8 kV substations, representing an ideal framework for small-scale distribution of photovoltaic power generation.

Identification of Alternative Sites Based on a Geographic Suitability Index Map
The geographic suitability index map of the area under evaluation ( Figure 5) was prepared in The MCDM method: the application of MCDM methods using business criteria is widely used in energy planning [22]. Among other advantages, it allows stakeholders to reach an agreement based on the evaluation of numerous conflicting criteria from multiple dimensions. In the decision-making process for locating decentralized energy projects, typical for several renewable energies, some factors significantly influence the choice of the multi-criteria method:  The context of the decision-making process of companies investing in the energy sector, in which there are usually multiple decision makers or investors, thus hindering the reaching of a consensus and leading decision makers to act more individually due to economic issues;  Uncertainty plays an important role in the project's decision-making problems that involve incomplete and imperfect information of real-world situations [23];  Conflicts among decision makers and investors regarding the real importance of interests from multiple dimensions, represented in practice by those criteria with accompanying business-related subcriteria;  The uncertainties surrounding the real values of the multiple interest dimensions regarding alternative sites for photovoltaic plants.
The choice of the MCDM method for the installation of photovoltaic plants in vast and large regions with geo-economic and socio-environmental heterogeneities must take into account all the aspects that were mentioned above.

Identification of the Problem
The area under evaluation for selecting sites alternatives for installing photovoltaic plants spreads throughout 199,315 km 2 in Southern Brazil ( Figure 5). It consists of 10 mesoregions with a population (estimated by the Brazilian Institute of Geography and Statistics-IBGE) of 11,348,937 inhabitants distributed across 399 municipalities, 335 of which with populations under 30,000. The area comprises a large-scale energy distribution network, mostly with 34.5/13.8 kV substations, representing an ideal framework for small-scale distribution of photovoltaic power generation.
The uncertainties surrounding the real values of the multiple interest dimensions regarding alternative sites for photovoltaic plants.
The choice of the MCDM method for the installation of photovoltaic plants in vast and large regions with geo-economic and socio-environmental heterogeneities must take into account all the aspects that were mentioned above.

Identification of the Problem
The area under evaluation for selecting sites alternatives for installing photovoltaic plants spreads throughout 199,315 km 2 in Southern Brazil ( Figure 5). It consists of 10 mesoregions with a population (estimated by the Brazilian Institute of Geography and Statistics-IBGE) of 11,348,937 inhabitants distributed across 399 municipalities, 335 of which with populations under 30,000. The area comprises a large-scale energy distribution network, mostly with 34.5/13.8 kV substations, representing an ideal framework for small-scale distribution of photovoltaic power generation.

Identification of Alternative Sites Based on a Geographic Suitability Index Map
The geographic suitability index map of the area under evaluation ( Figure 5) was prepared in Quantum GIS (QGIS) by the Department of Cartographic Engineering of the Federal University of Paraná (UFPR). This was prepared through weighted aggregation, applying the MCDM weighted linear combination (WLC) in the overlay of geographic subcriteria suitability index maps, which were obtained from the fuzzification of the irradiation, temperature, land slope, land aspect, and distance from roads maps, with the subsequent elimination of crisp geographic restrictions of conservation units, land use and distance from cities and towns.
The geographic suitability index map of the area consists of pixels with degrees of fuzzy pertinence between 0 and 1, classified into categories from unsuitable to most suitable (unsuitable < 0.188, 0.1889 ≤ low suitability < 0.376, 0.3777 ≤ moderate suitability < 0.584, 0.585 ≤ suitable < 0.752, most suitable ≥ 0.753). The geographic suitability index map of the area under evaluation allowed the identification of the 10 alternative sites (Table 4) with the highest suitability indexes. These sites, however, presented different energy production capacities: across the 10 mesoregions, they were located next to paved roads assisted by 13.8 or 34.5 kV electricity networks and with areas between 20,000 m 2 and 100,000 m 2 , sufficient for small-scale photovoltaic generation from 1 to 5 MW.

Identification of Stakeholders
The following group of 10 stakeholders with 20 to 35 years of experience was invited to participate in the surveys: • Two male and two female engineers from the superintendence of energy distribution of an energy utility, with technical-economic and socio-environmental experience and a point of view aligned with government and investor interests alike; • Four managerial engineers from three multinational companies that supply the energy sector, with technical-economic experience and a point of view aligned with the interests of private supplier industries; • A state social worker with a point of view aligned with the government's and society's interests alike; • A chemical industry director engineer with technical administrative experience with a point of view aligned with environmental interests.
All of them accumulate strategic-level decision-making experience in large national and multinational companies with thousands of employees. The group is also experienced in state contracts and projects, mainly in the energy sector, representing the interests of the government, investors and society.
Considering the homogeneity of experience of the stakeholder group, the information collected in the surveys was treated with the same degree of importance.

Identification of Criteria and Subcriteria
The stakeholders agreed on the use of everyday business-related subcriteria drawn from the literature in the field and with the availability of information in public or private databases. After further evaluation of the business subcriteria by the stakeholders, the following were eliminated or replaced: • The economic subcriteria of net present cost and cost-benefit were eliminated because they need extensive and complex economic feasibility studies for each site alternative, in addition to the limited availability of national public data for smaller projects involving photovoltaic plants;

•
The technical subcriterion of net generation efficiency was eliminated for presenting calculations with inaccurate future forecasts. The annual availability factor was replaced by the duration of interruption per consumer unit (DEC), for having values for each alternative site in a public database; • Technical subcriteria for reliability, maturity and safety of solar technology were eliminated due to the subcriterion Duration of interruption per consumer unit (DEC) having already encompassed failures of this type of power generation; • Political and legal risk subcriteria were eliminated because it is a region with the same political and legal risks.
The subcriteria that were maintained were the ones that stemmed from literature and were available in a public access database of private or state source, or of simple calculation, mandatory for administrative daily use in the business (Table 5).

Identification of the Multi-Criteria Decision-Making Method
Taking into account, as previously mentioned, the interests and conflicts among stakeholders, also considering the insecurity and risks of investors in the energy sector, a qualitative and diffuse MCDM method was applied in this study. The method considers the importance and values of photovoltaic plant location criteria, regarding a compensatory nature in their internal aggregation procedure.
The option for the fuzzy complex proportional assessment (COPRAS-F) method relied on its easy application as well as its linguistic and diffuse nature, which facilitate the resolution of conflicts during the decision-making process. Furthermore, this method allows for good performance in some subcriteria to counterbalance poor performance in other subcriteria, and has been applied by many researchers in the last decade [23].

Description of the Fuzzy Complex Proportional Assessment Method
Zavadskas and Antucheviciene introduced COPRAS-F in 2007. It is one of the fuzzy multicriteria decision-making techniques used to determine the importance of each criterion and classify alternatives through linguistic terms associated with fuzzy membership functions [35][36][37]. The method presents the following characteristics:

•
It allows benefit and cost criteria to be used in the decision matrix [38] and criteria values to be normalized so that different units of measure can be used and compared; • It determines the best among a set of viable alternatives and the remaining alternatives as reasons for the best one, assuming direct and proportional dependence on the degrees of significance and utility of the alternatives in a criteria system [38]; • It uses fuzzy logic that can model the existing uncertainty by applying a linguistic variable associated with fuzzy numbers, identifying the degree of relevance in a closed range from 0 to 1, which makes it a useful concept to deal with complex or not well-defined situations; A fuzzy set is a powerful mathematical tool for handling the existing uncertain in decision making, and is a general form of a crisp set. A fuzzy number belongs to the closed interval 0 and 1, in which 1 addresses full membership and 0 expresses non-membership; whereas, crisp sets only allow 0 or 1 [37].
To apply the MCDM COPRAS-F method, the model follows three steps, as shown in Figures 6,  8 and 9 respectively. In the first step, the decision fuzzy matrix ( Figure 6) must be structured based on the choice of the fuzzy membership function, followed by the structuring of the linguistic terms of criteria to apply the fuzzy searches of criteria. Different fuzzy membership functions may be applied from the fuzzy modeling of the problem. An option was made for the application of fuzzy triangular membership functions (Figure 7) [37], which presents the advantages of facilitating the surveys with stakeholders, making it easier to implement and process information in a fuzzy environment than other fuzzy pertinent functions [37]. To represent the triangular membership function (Figure 6), the fuzzy numbers are taken spread left d1, center d2 and spread right d3, being: Chatterjee and Bose [35] described the COPRAS-F mathematical concept as follows: • Choice of the linguistic labels for criteria and for alternatives (Table 6) [35]: both the importance of the criteria weights and the qualitative categorization of alternatives are collected through research along the stakeholders, using linguistic variables expressed in triangular fuzzy numbers. Table 6. Linguistic terms for criteria and alternatives.

Linguistic Terms for Criteria Fuzzy Numbers Linguistic Terms for Alternatives Fuzzy Numbers
Very High (VH) (0.9, 1.0, 1. • Development of the fuzzy decision matrix: for each decision maker, a fuzzy number decision matrix will be obtained with m rows and n columns, considering that the number of criteria is n and the number of alternatives is m: .
In such case: • DM k is the kth decision maker • A m is the mth alternative • C n is the nth criterion During the second step, after obtaining the fuzzy decision matrix, the weighted normalized defuzzied decision matrix (Figure 8) must be obtained through calculations of fuzzy (weights and alternatives) criteria average followed by the defuzzification process by the center of the area method with final normalization.

•
Subsequently, applying the fuzzy weighted averaging (FWA) operator, the responses from decision makers are aggregated (Xu and Da, 2003), considering q k the diffuse significance coefficient of the k th decision maker. If k is homogeneous, then q k = (1/k, 1/k, 1/k). Therefore, the aggregated fuzzy decision matrix using Equation (2) will be represented as follows: where: • Defuzzification: to defuzzify the fuzzy weight of each criterion into crisp values, the centre of area (COA) defuzzification method is applied. This method is simple and practical, with no need of including the preferences of any evaluator. The value of best non-fuzzy performance (BNP) for the fuzzy number can be found by using the following equation for a triangular fuzzy function: • Normalization of defuzzified decision-making matrix: the normalized values of this matrix are calculated by the equation: • Calculation of criteria weight: defuzzify the fuzzy weight of each criterion (w 1 , w 2 , w 3 ) in crisp values (w j ): • Calculation of the weight of the j th criterion is depicted as follows: • Calculation of the weighted normalized fuzzified decision matrix: the weighted normalized fuzzified values are calculated by the equation: In which: q j is the weight of the j-th criterion. During the third step, after obtaining the weighted normalized defuzzified matrix, the significance and definition of priority through calculations of the sum Pj, sum Ri and optimization criterion must be obtained (Figure 9).

•
Calculation of the sum P j : for each alternative of the fuzzified normalized and weighted values of the criteria of benefits, as they are preferable in optimization aimed at maximization In which: l is the number of criteria that must be maximized.
• Calculation of the sum R i : for each alternative of the fuzzified normalized and weighted values of the criteria of costs, as they are preferable in optimization aimed at minimization: In which: (m − l) is the number of criteria that must be minimized.
• Calculation of significance of the alternatives or relative weight of each alternative of Q i project: • Calculation of the optimization criterion K, i.e., the maximum Q i (Qmax) • Definition of priority regarding design alternative: the greater the relative weight of the alternative Q i , the higher the priority or significance of the design alternative. This way, the Q max significance of the more rational alternative will always be that of the highest degree of satisfaction.

•
Calculation of the degree of utility N j from 0% to 100% of the alternatives: calculated by comparing each alternative analyzed A j with the most efficient alternative, which has a degree of utility equal to 100%, so it is possible to assess the extent to which an alternative A j is better or worse than the others being compared: In which Q i and Q max are the significance of the project alternatives under evaluation. As the priority of the analyzed alternative increases or decreases, its degree of utility fluctuates as well.

Application of the Fuzzy Complex Proportional Assessment Method
To implement the method Complex Proportional Fuzzy Assessment (COPRAS-F) the following phases described in the items below are necessary.

Research of Linguistic Classifications for Criteria Importance
A survey for linguistic classifications of criteria importance, based on linguistic variables expressed in triangular fuzzy numbers (Table 6), was conducted with the four stakeholders from the energy utilities (D1, D2, D3 and D4), from which we reached the linguistic matrix of the importance of subcriteria weights (Appendix A).
It is important to understand that to apply the surveys, a limited number (four) decision-makers were chosen due to the requirement of extensive experience in strategic decision-making in a utility with thousands of employees, billions of dollars in revenue, and millions of consumers, occupying or having already occupied a high position in the utility.

Survey of Linguistic Classifications for Alternatives
We conducted a survey of linguistic classifications for alternatives based on linguistic variables expressed in triangular fuzzy numbers (Table 6) with the four stakeholders from the energy utilities (D1, D2, D3 and D4), from which we reached the linguistic matrix of alternative (Appendix B).

Fuzzified Linguistic Matrix of Weight Importance and Defuzzified by the Average with Calculation of Criteria Weight
From the results of the survey for linguistic classifications of criteria importance with the decision makers, a fuzzified linguistic matrix for the criteria was generated, later defuzzified by the mean, in which the center of area method (COA) was applied, obtaining the normalized weights (Table 7). Table 7. Calculation of criteria weights using the center of area (COA). From the linguistic matrix of alternatives (Appendix B), the fuzzified matrix was generated, later defuzzified, normalized and weighted decision matrix (Table 8). Finally, the calculations were performed in the defuzzified, normalized and weighted matrix to determine the ranking of the alternatives of locations for the installation of photovoltaic plants (Table 9).

Discussion of Application and Results of the Application of the MCDM COPRAS-F METHOD
The facilities in decision making were analyzed, which became clear during the research process and the format of the results obtained by applying the COPRAS-F method (Table 9), proving the practicality of the COPRAS-F method when combining structures from: • COPRAS: due to the direct application of benefit and cost subcriteria, in addition to calculating the degree of utility (significance) of the alternatives. This demonstrated, in percentage, to what extent an alternative is better or worse than other alternatives, which visually helps in decision making; • FUZZY: practicality of calculating fuzzy logic in resolving conflicts between decision makers, without the need for consensus meetings.
The COPRAS-F method has demonstrated, during the stakeholders' survey process, great applicability for companies focused on energy sector projects, in which there are multiple decision-makers with different interests, and that often, when failing to reach a consensus, prefer to take more personalized decisions during meetings.
Regarding the results of the COPRAS-F classification (Table 10), it has been observed that the alternative (A7) presented the best solution with 8.61 GWh/year and a utility degree being 4% higher than the second-best solution (A1) with 8.47 GWh/year, 7% higher than the third best solution (A6) and 38% higher in relation to the worst solution (A3). The technical priority of power generation in relation to the other business criteria is evident.

Conclusions
This section covers the decision-making model proposed here, and suggested improvements for future models.

Analysis and Discussion of the Proposed Decision-Making Model
Aiming at this advance in relation to current methods, the Fuzzy COPRAS method was applied in this article with business criteria in alternatives of sites identified from a geographic suitability index map of Southern Brazil, thus establishing a ranking among alternative areas with the best geographic characteristics of 10 mesoregions (Table 10).
The differences observed between the GIS-MCDM ranking, established by the geographic suitability index, and the ranking obtained when applying COPRAS-F (Table 11), demonstrate the need for integrating both the geographic and business points of view.
The best locations from the geographic point of view for the installation of photovoltaic plants were ranked based on the geographic suitability index, which was obtained through the GIS-MCDM, and they were not the best locations from the business point of view, as ranked by MCDM Fuzzy COPRAS method (Table 11). This occurs especially in the case of large territorial extensions that present heterogeneity in economic, social, and environmental aspects. The contributions of the proposed model for the location of areas in vast and heterogeneous extensions of land for the installation of photovoltaic plants are:

•
Eliminate the weaknesses of the GIS and GIS-MCDM models in the impossibility of mapping business subcriteria, and the weaknesses of the MCDM method in the impossibility of mapping geographical subcriteria; • Demonstrate the complementarity between geographical subcriteria, which can be excellent in the selection of geographically suitable areas, and business subcriteria, which can be excellent for the classification of site alternatives regarding the viability of an enterprise.
On the other hand, the application of the present model has not been shown to make sense in selecting and classifying alternative areas for installation of photovoltaic plants in:

•
Regions characterized by economic, social, and environmental homogeneity, in which case the direct use of the GIS or GIS-MCDM methods prevails; • Regions characterized by climatic, geomorphological, location and land use homogeneity, in which case the direct use of MCDM methods prevails.

Analysis and Discussion of Future Decision-Making Models
Here are some suggestions for the next steps in improving this MCDM model: • Develop research involving structuring criteria for localization applications for distributed renewable energies, primarily wind energy, wind-photovoltaic and hydraulic hybrids aimed at small hydroelectric plants;

•
Apply other MCDM or MCDM-FUZZY methods with new business sub-criteria determined by business literature or policies or consensus of decision makers [36,[39][40][41][42], such as: political and economic instability, changes in policy changes in policy premiums, financing policy, environmental releases.