Abstract
The extensions of ordinary fuzzy sets such as intuitionistic fuzzy sets (IFS), Pythagorean fuzzy sets (PFS), and neutrosophic sets (NS), whose membership functions are based on three dimensions, aim at collecting experts’ judgments more informatively and explicitly. In the literature, generalized three-dimensional spherical fuzzy sets have been developed by Kutlu Gündoğdu and Kahraman (2019), including their arithmetic operations, aggregation operators, and defuzzification operations. Spherical Fuzzy Sets (SFS) are a new extension of Intuitionistic, Pythagorean and Neutrosophic Fuzzy sets, a SFS is characterized by a membership degree, a nonmembership degree, and a hesitancy degree satisfying the condition that their squared sum is equal to or less than one. These sets provide a larger preference domain in 3D space for decision makers (DMs). In this paper, our aim is to extend classical VlseKriterijumska Optimizacija I Kompromisno Resenje (VIKOR) method to spherical fuzzy VIKOR (SF-VIKOR) method and to show its applicability and validity through an illustrative example and to present a comparative analysis between spherical fuzzy TOPSIS (SF-TOPSIS) and SF-VIKOR. We handle a warehouse location selection problem with four alternatives and four criteria in order to demonstrate the performance of the proposed SF-VIKOR method.
Introduction
After the presentation of ordinary fuzzy sets by Zadeh (1965), they have been very popular in almost all branches of science. Various researchers [6, 57] have developed several extensions of ordinary fuzzy sets as given in Fig. 1 with a historical order. Yager (2013) has renamed Atanassov’s intuitionistic fuzzy sets of second type (IFS2) as Pythagorean fuzzy sets (PFS). Hence PFS and IFS2 mean the same fuzzy sets thereafter [44].

Extensions of fuzzy sets.
In recent years, several researchers have utilized these extensions in the solution of multi-criteria decision making problems. A classification of some recent publications with respect to the types of fuzzy extensions is as follows:
Type-2 fuzzy sets (T2FS): The concept of a type-2 fuzzy set was introduced by Zadeh (1975) as an extension of the concept of an ordinary fuzzy set called a type-1 fuzzy set [28]. Such sets are fuzzy sets whose membership grades themselves are type-1 fuzzy sets; they are very useful in circumstances where it is difficult to determine an exact membership function for a fuzzy set [20, 45].
Intuitionistic fuzzy sets (IFS): Intuitionistic fuzzy sets introduced by Atanassov (1986) enable defining both the membership and non-membership degrees of an element in a fuzzy set [8, 58].
Hesitant fuzzy sets (HFS): Hesitant fuzzy sets can be used as a functional tool allowing many potential degrees of membership of an element to a set. These fuzzy sets force the membership degree of an element to be possible values between zero and one [1, 62].
Pythagorean fuzzy sets (PFS): Atanassov’s intuitionistic fuzzy sets of second type (IFS2) or Yager’s Pythagorean fuzzy sets are characterized by a membership degree and a nonmembership degree satisfying the condition that the square sum of its membership degree and nonmembership degree is equal to or less than one, which is a generalization of Intuitionistic Fuzzy Sets (IFS) [15–17, 65].
Neutrosophic sets (NS): Smarandache (1999) developed neutrosophic logic and neutrosophic sets (NSs) as an extension of intuitionistic fuzzy sets. The neutrosophic set is defined as the set where each element of the universe has a degree of truthiness, indeterminacy and falsity [6, 41].
The spherical fuzzy sets (SFS) have been recently introduced by Kutlu Gündoğdu and Kahraman (2019) [9]. These sets are based on the fact that the hesitancy of a decision maker can be defined independently from membership and nonmembership degrees, satisfying the following condition [9]:
On the surface of the sphere, Equation (1) becomes
The idea behind SFS is to let decision makers to generalize other extensions of fuzzy sets by defining a membership function on a spherical surface and independently assign the parameters of that membership function with a larger domain. SFS are a synthesis of PFS and NS.
The VIKOR method was developed by S. Opricovic (1998) as an MCDM method to solve a discrete multi-criteria problem with non-commensurable and conflicting criteria [49]. The method aims at determining a compromise solution for ranking the alternatives by considering those criteria. A compromise solution is a feasible solution nearest to the ideal solution [50]. The classical VIKOR method has been extended to its fuzzy versions by using various types of fuzzy sets such as type-2 fuzzy VIKOR [36], intuitionistic fuzzy VIKOR [22, 53], and hesitant fuzzy VIKOR [13], Pythagorean fuzzy VIKOR [56].
To the best of our knowledge, spherical fuzzy VIKOR method has not yet been developed in the literature. The motivation of this paper is to extend the VIKOR method under spherical fuzzy environment, assuming that membership, non-membership and hesitancy degrees are independently assigned and their squared sum is equal to at most 1. In this study, the proposed SF-VIKOR decision making model is applied to a warehouse site selection problem since the experts’ evaluations of location alternatives involve impreciseness and vagueness. An illustrative case study containing four criteriaand four alternatives is presented in the application section.
The originality of the paper comes from the presentation of a novel SF-VIKOR and the application of the proposed method to warehouse site selection problem. The SF-VIKOR enables decision makers to independently reflect their hesitancies in the decision process by using a linguistic evaluation scale based on spherical fuzzy sets.
The rest of this paper is organized as follows. Section 2 includes the introductory definitions and the preliminaries on SFS. Section 3 summarizes a literature review on VIKOR method. Section 4 includes our proposed MCDM method: Spherical Fuzzy VIKOR (SF-VIKOR). Section 5 applies SF-VIKOR method to a warehouse site selection problem. Section 6 includes a comparative analysis between SF-VIKOR and SF-TOPSIS and a sensitivity analysis for SF-VIKOR. Finally, the study is concluded in the last section.
Intuitionistic and Pythagorean fuzzy membership functions are composed of membership, non-membership and hesitancy parameters, which can be calculated by

Geometric representations of IFS, PFS, NS and SFS.
In this section, we give the definition of SFS and summarize spherical distance measurement, arithmetic operations, aggregation operators and defuzzification operations.
A spherical fuzzy set
For each u, the numbers
On the basis of relationship between SFS and PFS, Kutlu Gündoğdu & Kahraman (2019) further define some novel operations for SFS as below [9]:
λ .
One of the most popular MCDM methods, VIKOR focuses on ranking and sorting a set of alternatives against various, or possibly conflicting and non-commensurable decision criteria assuming that compromising is acceptable to resolve conflicts. Similar to some other MCDM methods, VIKOR relies on an aggregating function that represents closeness to the ideal solution. It introduces the ranking index based on the particular measure of closeness to the ideal solution and uses linear normalization to eliminate units of criterion functions [50].
As seen in Table 1, VIKOR has been integrated with different fuzzy extensions in a few fields of application in recent years. The publications on VIKOR method are summarized in Table 1.
A literature review on fuzzy VIKOR
A literature review on fuzzy VIKOR
A literature review on fuzzy VIKOR method using SCOPUS database gives 4,595 published papers in all fields. Among these, 514 papers mention fuzzy VIKOR in their titles, abstracts, or keywords. Yearly distribution of papers on fuzzy VIKOR is given in Fig. 3. As it is clearly seen, about 150 papers more per year are added to the previous year VIKOR publications.

Fuzzy VIKOR studies based on years 2010–2018.
Figure 4 illustrates the article numbers on fuzzy VIKOR with respect to their authors. Authors E. K. Zavadskas (with 121 publications) from Vilnius Gediminas Technical University, G.H. Tzeng (with 89 publications) from Harbin Institute of Technology, and Z. Xu (with 54 publications) from Nanjing University of Information Science and Technology are the most productive three researchers in this field.

Researchers studying on fuzzy VIKOR.
Fuzzy VIKOR method has been used in many different areas. These areas can be categorized as follows: engineering, computer science, business management, mathematics, environmental science, decision sciences, social sciences, energy, earth and planetary sciences, material science and other areas as represented in Fig. 5. Especially, in the computer science and engineering areas, the method has been extensively used.

Fuzzy VIKOR studies with respect to their research areas.
The proposed SF-VIKOR method integrates the superiorities of Pythagorean fuzzy sets and neutrosophic sets as follows. It presents a larger definition space for the parameters of a SFS as in PFS. Hesitancy degree in a SFS can be assigned independently as in NS. The proposed sets eliminate the argued condition that the sum of membership, non-membership and hesitancy degrees is at most 3 in NS. Besides, the proposed method involves a new spherical fuzzy distance measurement, a new SF aggregation operator and a new SF defuzzification formula.
The proposed spherical fuzzy VIKOR method is composed of several steps as given in the following. Before giving these steps, we present the flow chart of the SF-VIKOR method in Fig. 6 in order to make it easily understandable.

SF-VIKOR proposed methodology.
A MCDM problem can be expressed as a decision matrix whose elements indicate the evaluation values of all alternatives with respect to each criterion under Spherical fuzzy environment. Let X ={ x1, x2, . . . . . . x
m
} (m ⩾ 2) be a discrete set of m feasible alternatives and C ={ C1, C2, . . . . . . C
n
} be a finite set of criteria, and w ={ w1, w2, . . . . . . w
n
} be the weight vector of all criteria which satisfies 0 ⩽ w
j
⩽ 1 and
Linguistic terms and their corresponding spherical fuzzy numbers
Construct aggregated spherical fuzzy decision matrix based on the opinions of decision makers.
Denote the evaluation values of Alternative X
i
(i = 1, 2 . . . . . m) with respect to criterion C
j
(j = 1, 2 . . . . . n) by
At this point, there are two possible ways to follow. The first possible way is to follow a partially fuzzy approach as follows:
Defuzzify the aggregated criteria weights by using the score function given in Equation (21). Normalize the aggregated criteria weights by using Equation (22).
The second way is to continue without defuzzifying the criteria weights. This approach is called full fuzzy approach.
For the SF-PIS, Equation (23) is used to find the maximum scores in the decision matrix. Based on the crisp maximum scores, the corresponding SF numbers are determined as in Equation(24).
For the SF-NIS, Equation (25) is used to find the minimum scores in the decision matrix. Based on the crisp minimum scores, the corresponding SF numbers are determined as in Equation (26).
Alternatively, spherical fuzzy weights can be used to continue with the full fuzzy approach as in Equation (29).
At this point, there are three possible distance formulas to follow. Euclidean distance (Equations (32) and (33)), Xu and Zhang’s distance (Equations (34) and (35)) and spherical distance (Equations (36) and (37)) can be used in this step[5, 64].
The indices min S i and min R i are related to a maximum majority rule, and a minimum individual regret of an opponent strategy, respectively. As well, v is introduced as the weight of the strategy of the maximum group utility. v is usually assumed to be 0.5.
Propose as a compromise solution the alternative (a′) which is ranked the bet by the measure Q i (minimum) if the following two conditions are satisfied:
C1:“Acceptable advantage”:
C2: “Acceptable stability in decision making”: Alternative a′ must also be the best ranked by S or/and R.
Facility layout and location has been a well-established research area within operations research and management science. However, the originality of the paper comes from the proposed novel SF-VIKOR method in the warehouse site selection. Our proposed methodology is applied to the selection of a facility warehouse location. For this goal, proposed four subprovinces (X: Cihanbeyli, A2: Aksehir, A3: Eregli, A4: Meram) in Konya as shown in Fig. 7 are evaluated.

The map of Konya.
After a comprehensive literature review, four criteria have been determined. Criteria are infrastructure (C1), markets (C2), costs (C3), and labor characteristics (C4). In this structure, while “costs” is a non-beneficial criterion, the rest of them are beneficial. First of all, the assessments for the criteria are collected from decision maker with respect to the goal, using the linguistic terms given in Table 2. In the evaluation process, three decision makers (DM1, DM2, and DM3) are included who are three experienced engineers in supply chain and logistics management. The weights of these decision makers who have different experience levels are 0.3, 0.2 and 0.5, respectively. All assessments are given in Tables 3–5.
Assessments of DM1
Assessments of DM2
Assessments of DM3
These judgments are aggregated using SWAM operator by considering the importance levels of decision makers. Aggregated decision matrix is obtained as in Table 6.
Aggregated decision matrix
The linguistic importance weights of the criteria assigned by DMs are shown in Table 7.
Importance weights of the criteria
The weight of each criterion obtained by using SWAM operator is presented in Table 8.
Aggregation criteria weights
After the weights of the criteria have been determined, the defuzified and normalized criteria weights are calculated by utilizing Equations (21) and (22) as given in Table 9.
Defuzzified and normalized criteria weights
Based on Table 6 and Equation (21), score function values are obtained as in Table 10. The highest values represent positive ideal solution while the lowest values represent negative ideal solution.
Score function values based on SWAM operator
According to the best and worst scores, the corresponding Spherical Fuzzy Positive Ideal Solution (SF-PIS) and Spherical Fuzzy Negative Ideal Solution (SF-NIS) are given in Table 11.
SF-PIS and SF-NIS
In the next step, based on Equations (38–40), we can calculate the values of S i , R i , and Q i for each alternative based on Euclidean distance. They are given in Tables 12 and 13.
The values of S i , R i , and Q i for each alternative based on Euclidean distance
The ranking of alternatives in ascending order by S i , R i and Q i based on Euclidean distance
In the next step, based on Equations (38–40), we can calculate the values of S i , R i , and Q i for each alternative based on spherical distance. They are given in Tables 14 and 15.
The values of S i , R i , and Q i for each alternative based on spherical distance
The ranking of alternatives in ascending order by S i , R i and Q i based on spherical distance
In the next step, based on Equations (38–40), we can calculate the values of S i ,R i , and Q i for each alternative based on Xu and Zhang’s distance. They are given in Tables 16 and 17.
The values of S i ,R i , and Q i for each alternative based on Xu and Zhang’s distance
The ranking of alternatives in ascending order by S i , R i and Q i based on Xu and Zhang’s distance
According to the full fuzzy spherical Fuzzy VIKOR method, the values of S i , R i , and Q i for each alternative and the rankings based on spherical distance, Zhang & Xu’s distance and Euclidean distance are given in Tables 18–20.
The values of S i , R i , and Q i for each alternative based on spherical distance (full fuzzy approach)
The values of S i , R i , and Q i for each alternative based on Xu and Zhang’s distance (full fuzzy approach)
The values of S i , R i , and Q i for each alternative based on Euclidean distance (full fuzzy approach)
All the results show that the ranking of the alternatives from the best to the worst is X3 > X1 > X4 > X2.
Kutlu Gündoğdu and Kahraman (2019) extended TOPSIS method to Spherical fuzzy TOPSIS (SF-TOPSIS) using spherical fuzzy sets [9]. SF-TOPSIS has been proposed in their paper and applied to the performance comparison of air condition suppliers successfully. In this study, our proposed methodology is compared with SF-TOPSIS for the site selection of warehouse.
In this comparison, the same aggregated decision matrix is used as given in Table 7. Xu and Zhang’s distances to positive and negative ideal solutions to SF-PIS and SF-NIS can be used as given Table 21.
Xu and Zhang’s distances to positive and negative ideal solutions
Xu and Zhang’s distances to positive and negative ideal solutions
Based on the classical closeness ratio formula, the ratios are calculated and presented in Table 22.
Closeness ratio of each alternative
The closeness ratios based on SF-TOPSIS method indicate that the best alternative is X3 and overall ranking is X3 > X1 > X4 > X2 similar with SF-VIKOR method.
We applied a sensitivity analysis by changing v, which is introduced as the weight of the strategy of the maximum group utility and the individual regret, and observed the robustness of the given decisions. Sensitivity analysis showed that very robust decisions have been obtained from SF-VIKOR as given in Fig. 8. Although the compromise solutions (Q) changed, the ranking of alternatives remained the same.

Sensitivity analysis.
VIKOR method is a MCDM method developed for the solution of complex systems. The method determines compromise solutions for the problems including conflicting criteria. We presented a spherical fuzzy extension of VIKOR to determine fuzzy compromise solution, where both criteria weights and alternative evaluations could be spherical fuzzy sets. Spherical fuzzy VIKOR is based on the spherical fuzzy aggregation functions and different distance measurement methods. We have also presented both a partially fuzzy approach and a full fuzzy approach to spherical fuzzy VIKOR. Comparative and sensitivity analyses showed that the ranking obtained from the SF-VIKOR method is quite robust and similar to SF-TOPSIS method. Spherical fuzzy sets let the decision makers assign membership, non-membership, and hesitancy degrees independently under the constraint that their squared sum is at most 1. Spherical fuzzy sets are a kind of combination of Pythagorean (IFS2) fuzzy sets and neutrosophic fuzzy sets.
For further research, we suggest interval-valued SF-VIKOR (IVSF-VIKOR) method to be developed and the obtained results to be compared with the SF-VIKOR in this paper. Aggregation operators and distance measurement equations for IVSF sets need to be developed prior to IVSF-VIKOR method.
