# Electric Vehicle Charging Station Location Decision Analysis for a Two-Stage Optimization Model Based on Shapley Function

Authors: Lifeng Yang; Zhongwei Cheng; Baojie Zhang; Fengyun Ma
Journal: Journal of Mathematics
 (2021)

---

## Abstract

The promotion of electric vehicles and their charging facilities to achieve carbon emission reduction is a research hotspot in the field of transportation. Aiming at the comprehensive decision of electric vehicle charging station (EVCS) location, this paper constructs an EVCS location evaluation index system that includes five indexes of grid load, traffic facilities, user preference, construction cost, and service radius. Firstly, we convert the exact number into interval judgment matrix, introduce Shapley fuzzy measure to calculate the weight of factors, and use the two-stage optimization model to further optimize the weight. Then, we combine the multiple criteria decision-making (MCDM) method in the Pythagorean fuzzy environment with partitioned normalized weighted Bonferroni mean (PFPNWBM) operator, and calculate the optimal ranking of alternatives according to the performance function and the accuracy function. Finally, a numerical example is used to analyze the difference between first-order linear optimization and two-stage optimization in alternative scheme evaluation, and the practical value of using model to evaluate EVCS location is verified.


---

## Body

## 1. Introduction

Under the pressure of carbon emission reduction and carbon neutrality, the promotion and use of new energy vehicles represented by electric vehicles are becoming increasingly popular. The location and construction of EVCS have attracted extensive attention in the academic community [1]. EVCS location is affected by many objective and subjective factors such as service radius [2], electric vehicle (EV) mileage [3], and user preference [4]. There is great uncertainty in the decision-making process [5], and it becomes increasingly impossible to assign accurate values in the evaluation scheme [6]. Fuzzy set theory [7] and MCDM model [8], which rely on subjective judgment data selection, are often used in the scheme selection of such problems [9].The evaluation form of ordinary fuzzy sets cannot adequately describe the actual situation related to a large number of fuzziness and uncertainties [10]. To cope with this challenge, Zeng et al. [11] proposed the concepts of weighted double hesitation fuzzy set (WDHFS) and weighted double hesitation fuzzy element (WDHFE) and applied them to group decision-making. Fan and Xiao [12] proposed a two-dimensional intuitionistic fuzzy set (TDIFS) model to integrate the uncertainty and reliability expression of intuitionistic fuzzy set (IFS), which can carry more information than simple indexes. In order to accurately reflect the relative relationship between alternatives, Tao et al. [13] extended the classical alternative queuing method (AQM) by using fuzzy preference relation (FPR). However, since IFS can only describe the case where the sum of membership and nonmembership does not exceed 1 [14], Yager and Abbasov [15] successfully extended IFS to Pythagorean fuzzy set (PFS), which allows a larger space to describe certain, uncertain, and hesitant information. PFS method is used in the selection of electric bicycle sharing recycling suppliers [16], on-site evaluation of household waste processing plants [17], and other optimization scheme selection decisions. Due to the uncertainty of the reality in the Pythagorean fuzzy environment, the choice of weight is often restricted by objective constraints. Based on this, Nie et al. [18] defined the constraints of incomplete known weights and built a linear programming model with the maximum performance value of the alternatives schemes as the optimal goal. The weight acquisition process of Shapley fuzzy metric was optimized.In addition, Wang et al. [19] developed a rapid charging facility planning model by considering battery deterioration and vehicle heterogeneity within the driving range. The corresponding MCDM conceptual framework is also widely used in location decisions. Ni et al. [20], in the framework of MCDM, proposed a comprehensive model to characterize vaccination decision-making process influenced by society. Gao et al. [21] suggested a MCDM method based on probabilistic language term set and regret theory. For the significance of Shapley fuzzy measure method in measuring the importance of elements and the interaction characteristics of combinations between elements, Tan et al. [22] used Shapley comprehensive authorization benefit distribution method to analyze the energy storage and charging pile project. And this method is often used in conjunction with MCDM to study program decision-making problems. Jing et al. [23] proposed a decision-making method of integrated product conceptual scheme based on Shapley value method and fuzzy logic model.Yang and Ma [24] believed that, with the increase of the order of the judgment matrix, the granularity of the scale is not meticulous enough, which limits the consistency of the judgment matrix. There is a trade-off between simplicity and precision when choosing a digital scale. Pasman and Rogers [25] proposed that expert opinions are enveloped in uncertainty, which includes accidental uncertainty (lack of accuracy and precision observation means) and cognitive uncertainty (quantity and quality of data). Xu and Hua [26] found that, in view of the differences in professional knowledge, problem familiarity, social background, and other aspects of the evaluation team members, they are more inclined to represent the evaluation situation by their own preferred information forms, so that the evaluation results often present a complex situation where multiple types of information coexisted. Due to the difficulty and contradiction of subjective judgment, as well as the time pressure, incomplete relative information of factors in the decision-making process, limitations of experts’ expertise, and other reasons [27], Teng and Liu [28] proposed a new large group decision method in probabilistic language to analyze the interaction between interest subgroups and the relationship between criteria. Li [29] believed that because different quantitative attributes have different physical dimensions, normalization processing is needed to transform the exact number into an interval fuzzy number with membership degree and nonmembership degree.In Section2, we briefly review the concepts of PFS and Shapley fuzzy measure on the basis of existing research results. In Section 3, based on the review of relevant research literature, a systematic and complete evaluation index system of EVCS is constructed, and then an empirical study is conducted on 5 regions with large charging demand. In Section 4, we establish a two-stage optimization model to identify the optimal Shapley fuzzy measure under the background that incomplete weight information is known, which is the basis for determining the standard weight. Then, based on the PFPNWBM operator, performance, and accuracy function, the comprehensive performance value of each alternative scheme is calculated, so as to determine the optimal scheme ranking. Section 5 provides detailed resolution steps. Section 6 illustrates the availability of the proposed method through an example. Finally, Section 7 compares the differences between the optimal Shapley fuzzy measure identified by the two-stage dynamic optimization model and the first-order linear optimization model through a numerical example. The research proves that the evaluation index system and the evaluation method have good practicability, accuracy, and flexibility.In summary, from the perspective of EVCS location decision analysis, this paper constructs a comprehensive and systematic EVCS location evaluation index system that includes five indexes of grid load, traffic facilities, user preference, construction cost, and service radius. And the heterogeneity of experts and the incomplete information of the weights in the evaluation are well dealt with by establishing a two-stage optimization model to identify the optimal Shapley fuzzy measure, which provides some theoretical support and practical methods for EVCS location construction.


## 2. Preliminaries

In this section, some basic concepts and theoretical methods are introduced, which conclude IFS, PFS, and Shapley function.


### 2.1. Intuitionistic Fuzzy Set

Definition 1.
(see [30]). Let X be a finite universe of discourse. A=x,μAx,vAx|x∈X is defined as an IFS with the condition:(1)μAx,vAx:X⟶0,1,0≤μAx+vAx≤1.
TheμAx, vAx is the membership degree of x∈X indicating the degree of support and unsupported.


### 2.2. Pythagorean Fuzzy Set

Definition 2.
(see [31]). Let X be the domain of argument. A=x,AμAx,vAx,|x∈X is a PFS on X, where μAx,vAx:X⟶0,1 represents the membership and nonmembership of element X to A and 0≤μAx2+vAx2≤1, respectively. The degree of indeterminacy πAx is expressed by πAx=1−μAx2−vAx2.Definition 3.
(see [32]). Let a=Aμa,va be a PFN, and its score function and accuracy function are defined as follows:(2)Sa=μa2−va2,Ha=μa2+va2,whereSa,Ha∈0,1.Definition 4.
(see [33]). Let a1=Aμa1,va1 and a2=Aμa2,va2 be two PFNs. Then, the distance between a1 and a2is expressed as(3)da1,a2=12μa12−μa222+va12−va222+πa12−πa222.
As shown in Figure1, the applicable scope of PFS and IFS is 1/4 circle and right triangle in the first quadrant of the plane Cartesian coordinate system, respectively. Obviously, the use range of PFS is larger than that of IFS, so PFS has a wider computing space.Figure 1
IFS and PFS applicable scope comparison.


### 2.3. Shapley Function

Definition 5.
(see [34]). Let X be a finite set, and Px is the power set of X. The fuzzy measure μ defined on X:Px⟶0,1 satisfies the following conditions:(1)
μΦ=0,μX=1(2)
IfA,B∈Px and A⊆B, then μA⊆μB
Shapley function is used to measure the income distribution of each member alliance in the game. With respect to the fuzzy measureμA on the finite set defined as above, it can be expressed as(4)φiμ,N=μi=∑S⊆N/in−s−1!s!n!μS∪i−μi,∀i∈N.


## 2.1. Intuitionistic Fuzzy Set

Definition 1.
(see [30]). Let X be a finite universe of discourse. A=x,μAx,vAx|x∈X is defined as an IFS with the condition:(1)μAx,vAx:X⟶0,1,0≤μAx+vAx≤1.
TheμAx, vAx is the membership degree of x∈X indicating the degree of support and unsupported.


## 2.2. Pythagorean Fuzzy Set

Definition 2.
(see [31]). Let X be the domain of argument. A=x,AμAx,vAx,|x∈X is a PFS on X, where μAx,vAx:X⟶0,1 represents the membership and nonmembership of element X to A and 0≤μAx2+vAx2≤1, respectively. The degree of indeterminacy πAx is expressed by πAx=1−μAx2−vAx2.Definition 3.
(see [32]). Let a=Aμa,va be a PFN, and its score function and accuracy function are defined as follows:(2)Sa=μa2−va2,Ha=μa2+va2,whereSa,Ha∈0,1.Definition 4.
(see [33]). Let a1=Aμa1,va1 and a2=Aμa2,va2 be two PFNs. Then, the distance between a1 and a2is expressed as(3)da1,a2=12μa12−μa222+va12−va222+πa12−πa222.
As shown in Figure1, the applicable scope of PFS and IFS is 1/4 circle and right triangle in the first quadrant of the plane Cartesian coordinate system, respectively. Obviously, the use range of PFS is larger than that of IFS, so PFS has a wider computing space.Figure 1
IFS and PFS applicable scope comparison.


## 2.3. Shapley Function

Definition 5.
(see [34]). Let X be a finite set, and Px is the power set of X. The fuzzy measure μ defined on X:Px⟶0,1 satisfies the following conditions:(1)
μΦ=0,μX=1(2)
IfA,B∈Px and A⊆B, then μA⊆μB
Shapley function is used to measure the income distribution of each member alliance in the game. With respect to the fuzzy measureμA on the finite set defined as above, it can be expressed as(4)φiμ,N=μi=∑S⊆N/in−s−1!s!n!μS∪i−μi,∀i∈N.


## 3. The Evaluation Index System of EVCS Location

With the rapid development of EV technology, EV has become a major participant in the future development of new energy vehicles. However, due to the limited number of fast chargers and the increase of residents’ charging demands, EV charging operation may encounter the problems of increased queuing delay and increased charging operation costs [35]. To this end, governments around the world are promoting the use of EV in possible ways and issuing various preferential policies to encourage EV enterprises to cooperate with energy companies to build EVCS [4]. Therefore, establishing a comprehensive and complete EVCS location evaluation index system is the most effective way to promote the popularization of EV and solve the contradiction between supply and demand. Based on relevant literature review and analysis, we select the following five indicators from 10 indicators, such as grid load, traffic facilities, user preference, construction cost, service radius, EV ownership, driving mileage, policy support, regional service level, and elastic demand [2–4, 35–37].


### 3.1. Grid Load

Strong grid load capacity is the basic requirement of EVCS construction. Grid load refers to the ability of regional grid to withstand voltage load. The simultaneous charging of many vehicles in a fixed station by quick charging technology may produce huge additional demands on the power grid and may lead to power loss, which is bound to have certain requirements on the bearing capacity and supply level of grid.


### 3.2. Traffic Facilities

The location of EVCS directly affects the cost of driving on the road for vehicle owners (time and energy consumption on the way to the charging station). Electric vehicle users will give priority to EVCS with short travel distance and convenient road passage, so as to reduce driving costs. Therefore, the perfection of regional traffic facilities, namely, the minimization of customer travel costs is an important factor affecting EVCS location.


### 3.3. User Preference

User preference is mainly for users themselves, that is, drivers may form the habit of charging battery electric vehicle (BEV) after daily travel, including the battery capacity of the purchased electric vehicles, the acceptable charging time, the familiarity with public charging stations and travel mode (total travel duration of the next travel day, day interval between the start of the next travel and the start of the current charging event), etc.


### 3.4. Construction Cost

EVCS location should minimize the total cost and network loss cost of the charging station during the planning period as the goal of the optimal planning of the charging station [38]. The construction cost of charging station location includes investment cost, operating cost, and maintenance cost. Among them, investment cost refers to the investment cost per unit capacity of transformers and charging equipment, investment cost per unit capacity of other equipment except transformers and charging equipment, and land use cost. Operating cost includes charging cost, electric energy consumption cost of electrical equipment, filtering compensation cost of charging station, and labor cost. Maintenance cost refers to the maintenance cost per unit capacity of transformers and charging equipment and the maintenance cost per unit capacity of other equipment except transformers and charging equipment.


### 3.5. Service Radius

When the number of charging station is large, the smaller the Voronoi diagram area of the service area of the corresponding charging stations is, the lower the driving cost of the vehicle owner will be. In order to maximize EV battery cycle life and EVCS resources, on the premise of satisfying users’ charging needs, the service radius of charging station should be less than the actual distance between two adjacent charging stations under the premise of meeting the charging needs of the user and at most equal to the reasonable continuous driving distance of electric vehicle at a constant speed.


## 3.1. Grid Load

Strong grid load capacity is the basic requirement of EVCS construction. Grid load refers to the ability of regional grid to withstand voltage load. The simultaneous charging of many vehicles in a fixed station by quick charging technology may produce huge additional demands on the power grid and may lead to power loss, which is bound to have certain requirements on the bearing capacity and supply level of grid.


## 3.2. Traffic Facilities

The location of EVCS directly affects the cost of driving on the road for vehicle owners (time and energy consumption on the way to the charging station). Electric vehicle users will give priority to EVCS with short travel distance and convenient road passage, so as to reduce driving costs. Therefore, the perfection of regional traffic facilities, namely, the minimization of customer travel costs is an important factor affecting EVCS location.


## 3.3. User Preference

User preference is mainly for users themselves, that is, drivers may form the habit of charging battery electric vehicle (BEV) after daily travel, including the battery capacity of the purchased electric vehicles, the acceptable charging time, the familiarity with public charging stations and travel mode (total travel duration of the next travel day, day interval between the start of the next travel and the start of the current charging event), etc.


## 3.4. Construction Cost

EVCS location should minimize the total cost and network loss cost of the charging station during the planning period as the goal of the optimal planning of the charging station [38]. The construction cost of charging station location includes investment cost, operating cost, and maintenance cost. Among them, investment cost refers to the investment cost per unit capacity of transformers and charging equipment, investment cost per unit capacity of other equipment except transformers and charging equipment, and land use cost. Operating cost includes charging cost, electric energy consumption cost of electrical equipment, filtering compensation cost of charging station, and labor cost. Maintenance cost refers to the maintenance cost per unit capacity of transformers and charging equipment and the maintenance cost per unit capacity of other equipment except transformers and charging equipment.


## 3.5. Service Radius

When the number of charging station is large, the smaller the Voronoi diagram area of the service area of the corresponding charging stations is, the lower the driving cost of the vehicle owner will be. In order to maximize EV battery cycle life and EVCS resources, on the premise of satisfying users’ charging needs, the service radius of charging station should be less than the actual distance between two adjacent charging stations under the premise of meeting the charging needs of the user and at most equal to the reasonable continuous driving distance of electric vehicle at a constant speed.


## 4. EVCS Location Evaluation Method

### 4.1. Conversion of Exact Numbers

Comparison scale refers to the pairwise comparison of the relative importance of the indicators and gives the exact number judgment result according to one’s own preference. We use the 1–9 scale: the range ofaij is going to be in the range 1,2,…,9. When making qualitative pairwise comparisons, people usually have five distinct levels in mind. Assume that m schemes comprise scheme set X=x1,x2,…,xm and n evaluation indicators are divided into index set C=c1,c2,…,cn, and then the exact scoring value of the j evaluation object given by experts on the i indicator performance is aij. It can be transformed into an intuitionistic fuzzy set by the following formula:(5)μij=αiaijaimax,i∈Ωb,δiaiminaij,i∈Ωc,aimin≠0,δi1−aijaimaxi∈Ωc,aimin=0,vij=βiaijaimax,i∈Ωb,γiaiminaij,i∈Ωc,aimin≠0,γi1−aijaimax,i∈Ωc,aimin=0,where Ωb and Ωc represent the set of cost-type and benefit-type indicators, respectively. The determination method of aimin and aimax is as follows: (6)aimin=min1≤j≤naij,aimax=min1≤j≤naij.αi, βi, γi, and δi are four adjustment parameters, which are given by the actual decision-making environment and decision-maker’s experience [29]. Their values are all between 0 and 1, and they meet 0≤αi+βi≤1 and 0≤αi+βi≤1. Accordingly, the exact number aij can be transformed into an interval fuzzy number Aij=μij,vij.


### 4.2. Identifying the Optimal Shapley Fuzzy Measure

In practical MCDM issues, the weights of the attributes and decision makers (DM) are always partially known or completely unknown. To overcome this limitation, Wang et al. [39] developed the maximized deviation method (MDM) to handle more complex decision applications. It is well known that it is easier to choose better alternatives by ranking the fusion value, so as to more clearly reflect the differences between all alternatives. The main principle of MDM is that when an attribute has a significant impact on the ranking priority results, it is obvious that the attribute should be given greater weight in influencing the decision results. Conversely, if an attribute causes only a small change in the order of alternatives, it means that the attribute has less impact on the decision process and should be given less weight.Set the scheme setX=x1,x2,…,xm, n evaluation index set C=c1,c2,…,cn, and weighted vector set w=w1,w2,…,wn. Under the same criteria, the set of deviation values of the pairwise alternatives is daij,akj. Then we can calculate the degree of difference between each alternative through MDM. The calculation formula of deviation is as follows:(7)Djφ=∑i=1mDijφ=∑i=1m∑k=1mdaij,akjφjμ,Ph,j=1,2,…,n.In order to determine the optimal Shapley fuzzy measure, this paper establishes a two-stage dynamic programming model, which aims to achieve the maximum deviation of all standards.(8)maxDφ=∑j=1n∑i=1m∑k=1mdaij,akjφjμ,Phs.t.μPh=1,μS≤μT,∀S,T⊆Ph,S⊆T,μj∈Uj,∀j∈Ph.


### 4.3. Calculating the Comprehensive Performance Value

Letr,q≥0,ai=Aμai,vaii=1,2,…,m be a collection of PFNs, the aggregated results derived by PFPNWBM [18] operator are still a PFN, and it is expressed as(9)PFPNWBMa1,a2,…,am=1−∏h=1d1−1−1−μai2rμaj2qφi,hμ,Phφj,hμ,Ph/1−φi,hμ,Ph1/r+q1/d,∏h=1d1−1−∏i,j∈Phi≠j1−1−vai2r1−vaj2φi,hμ,Phφj,hμ,Ph/1−φi,hμ,Ph1/r+q1/d.


## 4.1. Conversion of Exact Numbers

Comparison scale refers to the pairwise comparison of the relative importance of the indicators and gives the exact number judgment result according to one’s own preference. We use the 1–9 scale: the range ofaij is going to be in the range 1,2,…,9. When making qualitative pairwise comparisons, people usually have five distinct levels in mind. Assume that m schemes comprise scheme set X=x1,x2,…,xm and n evaluation indicators are divided into index set C=c1,c2,…,cn, and then the exact scoring value of the j evaluation object given by experts on the i indicator performance is aij. It can be transformed into an intuitionistic fuzzy set by the following formula:(5)μij=αiaijaimax,i∈Ωb,δiaiminaij,i∈Ωc,aimin≠0,δi1−aijaimaxi∈Ωc,aimin=0,vij=βiaijaimax,i∈Ωb,γiaiminaij,i∈Ωc,aimin≠0,γi1−aijaimax,i∈Ωc,aimin=0,where Ωb and Ωc represent the set of cost-type and benefit-type indicators, respectively. The determination method of aimin and aimax is as follows: (6)aimin=min1≤j≤naij,aimax=min1≤j≤naij.αi, βi, γi, and δi are four adjustment parameters, which are given by the actual decision-making environment and decision-maker’s experience [29]. Their values are all between 0 and 1, and they meet 0≤αi+βi≤1 and 0≤αi+βi≤1. Accordingly, the exact number aij can be transformed into an interval fuzzy number Aij=μij,vij.


## 4.2. Identifying the Optimal Shapley Fuzzy Measure

In practical MCDM issues, the weights of the attributes and decision makers (DM) are always partially known or completely unknown. To overcome this limitation, Wang et al. [39] developed the maximized deviation method (MDM) to handle more complex decision applications. It is well known that it is easier to choose better alternatives by ranking the fusion value, so as to more clearly reflect the differences between all alternatives. The main principle of MDM is that when an attribute has a significant impact on the ranking priority results, it is obvious that the attribute should be given greater weight in influencing the decision results. Conversely, if an attribute causes only a small change in the order of alternatives, it means that the attribute has less impact on the decision process and should be given less weight.Set the scheme setX=x1,x2,…,xm, n evaluation index set C=c1,c2,…,cn, and weighted vector set w=w1,w2,…,wn. Under the same criteria, the set of deviation values of the pairwise alternatives is daij,akj. Then we can calculate the degree of difference between each alternative through MDM. The calculation formula of deviation is as follows:(7)Djφ=∑i=1mDijφ=∑i=1m∑k=1mdaij,akjφjμ,Ph,j=1,2,…,n.In order to determine the optimal Shapley fuzzy measure, this paper establishes a two-stage dynamic programming model, which aims to achieve the maximum deviation of all standards.(8)maxDφ=∑j=1n∑i=1m∑k=1mdaij,akjφjμ,Phs.t.μPh=1,μS≤μT,∀S,T⊆Ph,S⊆T,μj∈Uj,∀j∈Ph.


## 4.3. Calculating the Comprehensive Performance Value

Letr,q≥0,ai=Aμai,vaii=1,2,…,m be a collection of PFNs, the aggregated results derived by PFPNWBM [18] operator are still a PFN, and it is expressed as(9)PFPNWBMa1,a2,…,am=1−∏h=1d1−1−1−μai2rμaj2qφi,hμ,Phφj,hμ,Ph/1−φi,hμ,Ph1/r+q1/d,∏h=1d1−1−∏i,j∈Phi≠j1−1−vai2r1−vaj2φi,hμ,Phφj,hμ,Ph/1−φi,hμ,Ph1/r+q1/d.


## 5. The Decision Steps and the Characteristics

In this section, the detailed steps of the evaluation index system of EVCS location are provided:Step 1: establish a normalized judgment matrix.We adopt a 1–9 comparative scale and invited experts to give an exact score ofaij for the performance of the Jth evaluation object under the Ith indicator according to their preferences. If there are cost criteria, the decision matrix A=aijm×n need to be normalized by the following equations [18]:(10)aij∼=aij,for benefit criteria,aijc=va,ua,for cost criteria.Step 2: calculate the distance of each alternative.The distance between the alternative schemeI and k under the standard j can be calculated by the formula.Step 3: determine the optimal Shapley fuzzy measure.Using the above distance value and formula (7), the linear programming model is solved and the optimal Shapley fuzzy measure is obtained.Step 4: use formula (9) to calculate the comprehensive performance value of each alternative.Step 5: use formula (2) to calculate the value of score and accuracy function.The decision frame is shown in Figure2.Figure 2
EVCS site selection decision diagram.


## 6. An Illustrative Example

To test the practicality of the EVCS location evaluation index system, we select five optional areasX=x1,x2,x3,x4,x5: electric refueling area x1, high-speed service area x2, user residential area x3, commercial shopping area x4, and civic work area x5. The index set consists of five evaluation indicators: grid load C1, traffic facilities C2, user preference C3, construction cost C4, and service radius C5. Based on the correlation pattern, it is assumed that the criteria are divided into two categories P1=c1,c2,c3 and P2=c4,c5. The above potential alternatives are evaluated by experts according to these criteria using Pythagorean fuzzy information. This assessment provides incomplete information on importance weights:(11)U1=0.25,0.4,U2=0.4,0.5,U3=0.35,0.4,U4=0.45,0.6,U5=0.3,0.5.Now the method in the paper can be used to solve the problem. The detailed decision steps are shown as follows :Step 1: establish the normalized judgment matrix.The exact quantitative evaluation matrix given by experts is shown in Table1, which is converted into Table 2 by formula (10). As the cost types, they should be normalized into Table 3 by formula (5).Step 2: calculate the distance of the alternative.Under the same criteria, the distance values of pairwise schemes obtained by using the formula (3) are shown in Table 4.Step 3: determine the optimal Shapley fuzzy measure.The indexes are divided intop1 and p2, and their Shapley fuzzy measure is calculated, respectively. Taking the criteria in the first category p1 as an example, Shapley fuzzy measure of each criterion can be obtained from equation (4):(12)φ1μ,Ph=13μc1+16μc1,c2−μc2+μc1,c3−μc3+131−μc2,c3,φ2μ,Ph=13μc2+16μc1,c2−μc1+μc2,c3−μc3+131−μc1,c3,φ3μ,Ph=13μc3+16μc1,c3−μc1+μc2,c3−μc2+131−μc1,c2.Based on the above derivation, a linear programming model can be established as follows:(13)max−0.4165μc1−0.0775μc2+0.4939μc3−0.4939μc1,c2+0.0775μc1,c3+0.4165μc2,c3+1.4758s.t.μ1,2,3=1,μ1≤μ1,2,μ1≤μ1,2,μ1≤μ1,3,μ2≤μ1,3,μ2≤μ2,3,μ3≤μ2,3,μ1∈0.25,0.4,μ2∈0.4,0.5,μ3∈0.35,0.4.By solving the model, we can obtain(14)μ1=0.25,μ2=0.4,μ3=0.35,μ1,2=0.5,μ1,3=μ2,3=1.After a round of anonymous ratings by experts, the proceeds of the collaboration are returned to all parties.μ1+μ2≥μ1,2 is found, which means that the profit value after the two cooperation is lower than the value created separately. So they will choose the final cooperation. In order to guarantee the firmness and reliability of alliance cooperation, the benefits of all parties after cooperation should be at least greater than or equal to the sum of the benefits of all parties without cooperation. Therefore, on the basis of the original linear programming, we add constraint conditions:(15)μSm+μSn≤μT,∀Sm,Sn,T⊆Ph,Sm,Sn⊆T.We can obtainμ1=0.25,μ2=0.4,μ3=0.35,μ1,2=1,μ1,3=0.6,andμ2,3=1.According to formula (4), the Shapley fuzzy measure is calculated as φ1μ,Ph=0.2083,φ2μ,Ph=0.2833,andφ3μ,Ph=0.5084.Similar to the processing procedure of the criteria inp1, for Shapley’s fuzzy measure of C3 and C4 in p2, we can get(16)φ4μ,Ph=0.2,φ5μ,Ph=0.3.Step 4: use formula (9) to calculate the comprehensive performance value of each alternative.Whenr=q=1, the comprehensive performance value of each alternative can be calculated as(17)a1=0.1609,0.4738,a2=0.1065,0.7398,a3=0.0354,0.6870,a4=0.1628,0.5836,a5=0.0941,0.7273.Step 5: the score and accuracy function values are calculated using equation (2).Table 1
Accurate number rating.
c1c2c3c4c5x184367x245627x334265x477899x521435Table 2
Interval fuzzy number rating.
c1c2c3c4c5x1(0.26, 0.17)(0.20, 0.13)(0.12, 0.80)(0.43, 0.29)(0.60, 0.40)x2(0.33, 0.22)(0.33, 0.22)(0.60, 0.40)(0.33, 0.22)(0.43, 0.29)x3(0.20, 0.13)(0.20, 0.13)(0.30, 0.20)(0.60, 0.40)(0.60, 0.40)x4(0.12, 0.08)(0.12, 0.08)(0.12, 0.08)(0.60, 0.40)(0.30, 0.20)x5(0.12, 0.08)(0.12, 0.08)(0.60, 0.40)(0.60, 0.40)(0.60, 0.40)Table 3
Normalized rating.
c1c2c3c4c5x1(0.17, 0.26)(0.40, 0.60)(0.12, 0.80)(0.29, 0.43)(0.60, 0.40)x2(0.22, 0.33)(0.40, 0.60)(0.60, 0.40)(0.22, 0.33)(0.43, 0.29)x3(0.13, 0.20)(0.27, 0.41)(0.30, 0.20)(0.40, 0.60)(0.60, 0.40)x4(0.08, 0.12)(0.32, 0.48)(0.12, 0.08)(0.40, 0.60)(0.30, 0.20)x5(0.08, 0.12)(0.40, 0.60)(0.60, 0.40)(0.40, 0.60)(0.60, 0.40)Table 4
The same standard pairwise alternative distance.
c1c2c3c4c5da1j,a2j0.07230.00000.44290.11810.2230da1j,a3j0.03520.24720.09690.22300.0000da1j,a4j0.06730.16610.44290.22300.3460da1j,a5j0.06730.00000.44290.22300.0000da2j,a3j0.09650.24720.34600.32180.2230da2j,a4j0.12310.16610.00000.32180.1230da2j,a5j0.12310.00000.00000.32180.2230da3j,a4j0.03220.08110.34600.00000.3460da3j,a5j0.03220.24720.34600.00000.0000da4j,a5j0.00000.16610.00000.00000.3460Performance function and accuracy function are shown in Table5.Table 5
Score and accuracy ranking.
SaiHaiRankingx1−0.19860.25042x2−0.05660.21814x3−0.62270.62731x4−0.31410.36715x5−0.44250.47503


## 7. Comparison and Analysis

The application case shows that the constructed EVCS location evaluation system and evaluation method can prioritize the final alternatives and can better deal with problems such as the heterogeneity of experts and the incomplete information of the weights in the evaluation.From the point of view of indicator weight, the importance of the five indicators of “grid loadC1, traffic facilities C2, user preference C3, construction cost C4, and service radius C5” in EVCS location evaluation system is significantly different. The importance is ranked as follows: C4≺C1≺C2≺C5≺C3. It can be seen that, for government EVCS location, understanding the charging preference and travel modes of regional users is the key to EVCS obtaining maximum profits. The service radius ranks second in importance, showing that regional charging demand is the basic requirement for EVCS to perform its service functions.Judging from the evaluation results, the location of the 5 EVCS areas obtained by the two-stage optimization planning model is relatively good, especially the expressway, which has perfect infrastructure, convenient transportation, and huge passenger flow. Only through the first stage of decision-making, we get the comprehensive performance value of the five alternatives:(18)a1=0.2705,0.8107,a2=0.1065,0.7398,a3=0.0354,0.6870,a4=0.1326,0.6541,a5=0.0941,0.7273.The optimal ranking isx3≺x5≺x4≺x1≺x2.The results show that, due to the special complexity of the decision state, the two-stage optimization modifies the Shapley fuzzy measure, changes the original weight value, and then affects the optimal ranking of the final alternative. This supports the EVCS location evaluation method constructed by this research.


## 8. Conclusions

Using scientific methods to evaluate EVCS location is of great significance for improving the design, development, promotion, and application of EV. This paper builds a systematic EVCS location evaluation index system on the basis of comprehensive domestic and foreign research results on EV evaluation and various charging station location evaluation.First of all, in order to better deal with the heterogeneity, fuzziness, and uncertainty of expert evaluation information, this study introduces 1–9 evaluation scale for qualitative paired comparison and constructs interval fuzzy quantitative evaluation matrix by converting the exact number into the intuitionistic fuzzy set and other relevant theoretical results. Moreover, a two-stage optimization model to identify the optimal Shapley fuzzy measure is established by fully considering the incomplete knowledge of weights in the actual decision-making environment. And five representative EVCS locations are taken as the research objects for empirical application. Finally, the practicability, effectiveness, and flexibility of the evaluation index system and evaluation method are verified by analyzing the difference of the evaluation results between the first-order linear optimization and the two-stage optimization.Admittedly, this study also has some limitations, such as relatively small selection of research samples, insufficient geographical coverage, relatively small number of evaluators, and sample data collected. Although the heterogeneity of expert groups in the evaluation process has been considered, the diversity of evaluation criteria and the multigranularity of decision-making information have not been considered. In follow-up research, we will further expand research objects, enrich sample data, optimize research methods, and consider researching EVCS construction issues from the perspective of consumer experience, so as to enhance consumer experience and make it better to use charging services.



---

*Source: 5098378-2021-08-02.xml*