Intelligent Multi-Objective Public Charging Station Location with Sustainable Objectives

: This paper investigates a multi-objective charging station location model with the consideration of the triple bottom line principle for green and sustainable development from economic, environmental and social perspectives. An intelligent multi-objective optimization approach is developed to handle this problem by integrating an improved multi-objective particle swarm optimization (MOPSO) process and an entropy weight method-based evaluation process. The MOPSO process is utilized to obtain a set of Pareto optimal solutions, and the entropy weight method-based evaluation process is utilized to select the ﬁnal solution from Pareto optimal solutions. Numerical experiments are conducted based on large-scale GPS data. Experimental results demonstrate that the proposed approach can effectively solve the problem investigated. Moreover, the comparison of single-objective and multi-objective models validates the efﬁciency and necessity of the proposed multi-objective model in public charging station location problems.


Introduction
In the last two decades, nearly three-quarters of the carbon emissions from human behavior have come from burning fossil fuels. Particularly, the continued growth of motor vehicles worldwide has inevitably affected global crude oil demand and carbon emissions [1]. The transportation sector accounts for at least a quarter of total greenhouse gas (GHG) emissions in worldwide [2]. And changes in GHG emissions are closely related to the energy consumption [3]. A clean energy revolution is thus taking place around the world. Concerns about environmental issues [4] and resource shortages have led to a focus on renewable energy and sustainable development of transportation sector [5]. The electrification of road transport by electric vehicles (EVs) such as plug-in electric vehicles (PHEVs) and battery electric vehicles (BEVs) is widely regarded as one of the most effective ways to reduce the CO 2 emissions and the petroleum dependence of transportation sector [6]. Hybrid and plug-in electric vehicles can have significant emissions advantages over conventional vehicles. The life cycle emissions of BEVs or PHEVs depend on the source of electricity used for charging [7]. EVs generate lower GHG emissions than the most efficient gasoline vehicles, particularly if electricity generation is clean enough [8]. The carbon footprint of EVs is approximately only about 40% of traditional vehicles [9]. Thus, using rechargeable EVs to replace traditional vehicles is an inevitable trend because it not only decreases petroleum dependence but also reduces GHG emissions.
However, one of the most important difficulties that the large-scale deployment of EVs faces is consumer acceptance, which seriously affects the success of EVs market. Consumer acceptance of EVs depends on battery range as well as the convenience of charging, which are their main disadvantages. Therefore, in order to improve the market acceptance of EVs, drivers' range anxiety to reach their Sustainability 2018, 10 destinations and the convenience of charging need to be settled by effectively constructing public charging stations (CSs), which is the focus of EV public charging station location problem. This problem is concerned with the optimal placement of public charging stations in a given area so that one or more objectives are optimized without violating certain constraints. The construction of CSs involves many aspects of considerations and needs to balance the interests of multiple parties. In the trend of low-carbon and sustainable development, the triple bottom line principle is widely adopted to evaluate the performance of green and sustainable development from economic, environmental and social perspectives [10,11]. In economic perspective, building a CS is cost-expensive (e.g., more than one million yuan per station in China). In environmental perspective, the environmental protection and the GHG emissions reduction are important concerns. In social perspective, drivers want more CSs and thereby have less recharge waiting time to ease their mileage anxiety. Due to these conflicting concerns, it is thus necessary to consider these objectives simultaneously in the public charging station location problem for EVs.
With the development of information technology, Global Positioning System (GPS)-enabled devices are widely used in floating vehicles around the world, which can collect massive amounts of vehicle travel trajectory data. On the basis of these data, we can obtain the accurate travel trips of floating vehicles and their charging behaviors under the constraints of the proposed model, which is helpful in making more reliable decisions in charging station location problems.
Thus, this paper investigates a multi-objective charging station location problem for EVs based on vehicles' travel trips derived from massive GPS-enabled travel trajectories of floating vehicles. Sustainable objectives derived from the triple bottom line principle will be considered. Considering economic and environmental issues, the first two objectives of this research are to minimize both the number of CSs and the CO 2 emissions generated by all trips, which are helpful to control the construction costs and reduce carbon emissions. Considering the social issue, the drivers' satisfaction is crucial and the third objective is thus to minimize the average waiting time of EVs queuing at CSs, which is beneficial for promoting EVs by improving the convenience of charging.
The main contributions of this paper are listed as follows. First, a multi-objective intelligent location model is developed to handle the multi-objective public charging station location problem with realistic travel trips. The model considers three sustainable goals simultaneously base on the triple bottom line principle. Second, we investigate the differences between results generated by single-objective and multi-objective models, and validate the necessity of considering three objectives simultaneously proposed in this study.
The remainder of this paper is organized as follows. Section 2 reviews the related literature. Section 3 describes and formulates the investigated problem. In Section 4 we detail the proposed multi-objective intelligent location approach. Section 5 introduces the numerical experiments and discussions. Concluding remarks are made in Section 6.

Literature Review
With the development of the EV industry and increasing public awareness of environmental protection, the public charging station location problems for EVs have attracted more and more researchers' attention, belonging to one type of facility location problem that has been widely investigated since the 1990s [12,13]. For more detail, the interested readers can refer to review papers [14].
Among the location problems for the charging station, most studies addressed a single-objective location problem. Some studies focused on maximizing the electrification of itineraries. Dong et al. [15] adopted travel data of gasoline vehicles to simulate the travel patterns of BEVs, and minimized the total number of missed trips which cannot be completed by electricity. Shahraki et al. [16] aimed at maximizing the amount of vehicle-miles-traveled being electrified and developed an optimization model to determine the public CS locations for PHEVs. Some studies considered cost-related objectives. He et al. [17] aimed to minimize the cost incurred by missed trips and the total driving and recharging Sustainability 2018, 10, 3760 3 of 18 time cost. Li et al. [18] proposed a multi-period multi-path location model aiming to minimize the total construction cost of CSs. Yang et al. [19] allocated chargers for BEV taxis with the objective of minimizing the infrastructure investment. Charging demand-related objectives have also been considered in some studies. Tu et al. [20] proposed a demand coverage location model with the objective of maximizing the service capacities of EVs and CSs. He et al. [21] regarded the maximum flows that charge in route as the model objective with the consideration of driving range limitation and charging time required in stations.
Due to the realistic requirements in charging station location, some researchers considered multiple realistic objectives in EVs charging stations location. Wang and Wang [22] formulated a multi-objective model to maximize the coverage and minimize the cost. Yeo et al. [23] aimed to minimize the overall annual cost of investment and energy losses as well as maximize the annual traffic flow. Shinde and Swarup [24] focused on minimizing the charging cost and maximizing the amount of charging for optimal charging station allocation to different types of EVs. Zhang et al. [25] optimized the objectives of minimizing the total charging cost and load variance based on dynamic time-of-use price. Lou et al. [26] addressed the urban CS location problem to minimize the costs of service and construction while satisfying all the charging demands. Wang [27] put forward a model to maximize the EVs charging demands and minimize the total power loss and voltage deviation. Zhang et al. [28] proposed a multi-period capacitated flow refueling location model whose goal is to maximize EV demand and flow coverage. Spieker et al. [29] addressed the charging station placement problem with two objectives of maximizing the demand location coverage and the coverage of traffic flow. In previous multi-objective charging station location studies, the objectives related to construction costs and demand coverage have been considered most frequently, followed by energy-related objectives. However, the sustainable objectives of optimizing environmental, economic, and social goals simultaneously have not been handled in the charging station location literature so far.
The GPS trajectory data of floating vehicles has been adopted in several previous studies to handle the charging station location problem. Dong et al. [15] searched for the optimal CSs location solution based on travel patterns and charging behaviors collected from 18-month GPS data of over four hundred private gasoline vehicles in Seattle. Xu et al. [30] utilized nearly 500 BEV users' GPS data and charging information over an entire year in Japan to investigate the factors influence the charging mode and location. Luo et al. [31] explored the travel patterns and related spatial-temporal features of emissions through about 13,600 taxis' real time trajectories for one month in Shanghai for CSs deployment. Yang et al. [19] collected travel activities of over 7910 taxis for one week in Changsha, China to allocate charging stations for BEV taxis.
In summary, it is desirable to investigate the public charging station location problem with multiple objectives considering cost and environmental issues as well as drivers' satisfaction, especially based on massive GPS trajectory data of floating vehicles.

Problem Description
A city plans to develop public charging stations for EVs. There are I candidate charging stations (indexed by i) dispersedly distributed in the city, which can be determined according to the requirement of EV development or the local authority's development plan. The investigated problem needs to determine the optimal sub-set of these candidate stations as the final solution to the public charging station location problem, based on the public charging demands obtained by floating vehicles' travel trips. Three sustainable objectives, involving economic, environmental and social objectives, need to be optimized. There are J floating vehicles (indexed by j) in total in the city and each of them produces K j separate trips (indexed by k) in an examined period, which can be obtained based on GPS trajectory data of floating vehicles. The driver wishes to charge the vehicle after a trip is completed, rather than charge during a trip. Thus, at the end of each trip, the driver needs to decide whether the vehicle needs to, and can, be recharged. Trips that require charging form the charging demands, which can be obtained based on travel trips extracted from GPS trajectories of floating vehicles.
The type of vehicles in the investigated problem is PHEVs, which is one type of EVs that prefers to use the battery as driving force and use gasoline to complete trips only when the electricity is exhausted. Hybrid power effectively increases the EV's driving range and prevents trips obtained in reality from being not completed. Each PHEV is fully charged at the beginning of its first trip each day. The energy consumption in each PHEV is distance-dependent only and the electricity is only generated by thermal power. Each charging station has the same number of chargers and it is fully charged once the vehicle is recharged.
In the mathematical model, the notations in Table 1 are used.
a number much larger than E j β the charging rate of PHEVs (unit: kilowatt-hour per hour) λ CO 2 emissions rate by using electricity (unit: kilograms per kilowatt-hour) γ CO 2 emissions rate by using gasoline (unit: kilograms per liter) γ CO 2 emissions rate by using gasoline (unit: kilograms per kilometer) s the number of chargers in each CS t jk the time of the jth vehicle arrives at the CS after the kth trip Intermediate variables: P jk the remaining electricity of the jth vehicle after the kth trip, if it is less than zero, the electricity is insufficient for completing the kth trip (unit: kilowatt-hour) P jk the real remaining electricity of the jth vehicle after the kth trip,P jk ≥ 0 (unit: kilowatt-hour)

Mathematical Model
The mathematical model of the presented charging station location problems is formulated as follows: Minimize Subject to Z ijk = y i D ijk r jk , ∀i, j, k Sustainability 2018, 10, 3760 6 of 18 The model achieves a balanced trade-off among three sustainable objectives ( f 1 , f 2 , f 3 ) formulated in Equations (1)-(3). f 1 is the number of CSs considering the economy factor. f 2 represents the total CO 2 emissions (unit: kilogram) generated by PHEVs in a city from an environmental perspective and generated indirectly by electricity that converted by thermal power and directly by gasoline. f 3 is the average waiting time (unit: hour) of all trips with recharged at CSs. Constraint (8) ensures that the jth vehicle is recharged at the ith CS at the end of the kth trip only if constraints (9) and (10) are satisfied simultaneously. Constraint (9) stipulates that D ijk is 1 if the distance from the end location of the kth trip of the jth vehicle to the ith CS is less than the service radius of the CS; otherwise it is 0. Constraint (10) prescribes that a PHEV needs to be recharged if its remaining electricity falls below cE j . Constraint (11) ensures that a vehicle can only go to one CS for charging after a trip. Constraint (12) calculates the CO 2 emissions produced on the round trips to CSs. If the distance of the jth vehicle to the ith CS at the end location of the kth trip is not more than the distance could afforded byP jk , the CO 2 emissions in the round-trip are all generated indirectly by electricity; otherwise, C ijk is generated by both electricity and gasoline. Constraint (13) specifies the queue condition. If the j th vehicle arrives at the ith CS after its k th trip earlier than the jth vehicle's the kth trip and has not left when the jth vehicle arrives, the k th trip of the j th vehicle needs to be recharged before the jth vehicle's the kth trip, at ijk ij k = 1; otherwise at ijk ij k = 0. Constraint (14) counts the total number of vehicles in the ith CS when the jth vehicle arrives at the ith CS after its kth trip. T ijk is calculated in constraint (15). Constraint (16) indicates binary restrictions on the 2 decision variables and 3 intermediate variables.

Intelligent Multi-Objective Location Approach
Since the three objectives to be optimized are conflicting, a single optimal solution cannot be found and a multi-objective optimization is needed to find Pareto optimal solutions. Thus, this research proposes an intelligent multi-objective location approach to solve the investigated multi-objective CS location problem based on massive GPS data. This intelligent approach can be divided into two parts: an improved multi-objective particle swarm optimization (MOPSO) process and an entropy weight method-based evaluation process. In the first part, we construct an improved MOPSO process by integrating the discrete particle swarm optimization (PSO) algorithm and the non-dominated sorting technique to produce the Pareto optimal solutions to the investigated problem. This process provides a very efficient procedure in keeping the elitism optimization process as well as preserving the diversity, which assures a good convergence towards the Pareto-optimal front without losing the solution diversity. After obtaining the Pareto optimal solutions, in the second part, we evaluate these solutions and select the one with the highest comprehensive scores based on the entropy weight method [32], which determines the objective weight according to the variability of objective indicators obtained by Pareto optimal solutions' objective function values and avoids the subjectivity of determining objectives' weights. The flowchart of this approach is shown in Figure 1.
Sustainability 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sustainability as preserving the diversity, which assures a good convergence towards the Pareto-optimal front without losing the solution diversity. After obtaining the Pareto optimal solutions, in the second part, we evaluate these solutions and select the one with the highest comprehensive scores based on the entropy weight method [32], which determines the objective weight according to the variability of objective indicators obtained by Pareto optimal solutions' objective function values and avoids the subjectivity of determining objectives' weights. The flowchart of this approach is shown in Figure 1.

Improved MOPSO Process
The improved MOPSO process is employed to generate the Pareto optimal solutions to the investigated problem. The main procedures involved of Part 1 in Figure 1 are described as follows.

Improved MOPSO Process
The improved MOPSO process is employed to generate the Pareto optimal solutions to the investigated problem. The main procedures involved of Part 1 in Figure 1 are described as follows.

Initialization of Parameters and Individuals
Initialization of parameters and individuals are implemented first. The parameters utilized in the optimization process that are initialized include: the population size N = 150, the maximal number of iterations g max = 80 and the current iteration count g = 1. Each individual has a position string and a velocity string. The position string is represented as [y 1 , y 2 , . . . , y I ] to reveals the location plan whereas the velocity string is expressed as [v 1 , v 2 , . . . , v I ] to reflect the moving velocity of the individual, where I is equal to the number of candidate stations. v i is limited in [−10, 10]. The position string and velocity string of each individual are initialized randomly in step 1.

Generation of New Individuals
For each iteration g, the nth individual's best position Pbest n in its previous optimization process and the best position Gbest of the population P g−1 should be found for updating individuals' speeds and positions in P g . If the nth individual's positions obtained during all previous iterations in its optimization process are non-dominating, then we randomly select one as Pbest n . The best position Gbest of the population P g−1 is acquired through non-dominated sorting technique introduced in c 1 = c 2 = 2, the random numbers r 1 = 0.6 and r 2 = 0.3. rand( ) represents a random number uniformly distributed between [0,1].
Through the process above, new individuals produced and they form the offspring population Q g .

Objective Function Evaluation
Given candidate individual y and trips data extracted from massive GPS trajectory data as inputs, then individual's objective function values are calculated according to the method below. The value of objective function (1) is the summation of values of elements in the individual. The value of objective function (2) is calculated by formula (2) described in Section 3.2. For any trip, the driver decides whether the PHEV needs to be recharged. If the conditions in constraint (8) are satisfied, the vehicle goes to the CS and is fully recharged; otherwise, the driver starts the next trip without recharging. If there are multiple CSs that meet the conditions, the driver selects the nearest one. The remaining electricity of a trip is related to its previous trips' electricity consumption and calculated by formulas (4)- (6). During this process, information about charging behaviors such as t jk , the index i of this CS, E jk and C ijk are recorded if the jth vehicle goes to charge after the kth trip. The value of objective function (2) is calculated when completing all trips of all vehicles successively. Then, based on the charging information obtained when calculating function (2), according to the arrival time of a vehicle and the departure time of the preceding vehicles calculated by formula (7), the total number of vehicles q i t jk in the ith CS at that time can be obtained by formulas (12) and (13). In each CS, vehicles queue to recharge by their arrival time, the congestion situation and the waiting time are computed by formula (14). After finishing all trips of all vehicles, objective function (3) is thus solved. The fitness of the nth individual is represented as f n 1 , f n 2 , f n 3 because each individual has three objective values. The lower the three values, the better the fitness.

Individual Sorting by Non-Dominated Sorting Technique
The non-dominated sorting technique proposed in the non-dominated sorting genetic algorithm (NSGA-II) is adopted in sorting candidate individuals, which contains two procedures: fast-non-dominated sorting and calculation of crowding distance. Suppose the parent population is P g , the offspring population Q g is updated from P g adopting the method described in Section 4.1.2, and the combined population is R g = P g ∪ Q g (size = 2N). R g is divided into different non-domination levels (F1, F2, etc.) by performing fast-non-dominated sorting procedure. Individuals in the same level do not dominate one another. The higher the non-domination level of the individual, the better the fitness.
Crowding distance calculation procedure is put forward to compare fitness among individuals at the same nondomination level to ensure the results are more evenly distributed in the objective functions space and maintain population diversity. The crowding distance of the nth individual is the distance between two individuals adjacent to n in the objective functions space. Detailed calculation process can refer to [33].

Generation of Populations
For the initial population P 1 , the individuals in it are generated randomly based on the representation described in Section 4.1.1. Pbest n is initialized to the nth individual itself and Gbest is initialized to a random individual. For the population P g+1 , g ∈ [1, g max ], the candidate population S g is constructed by selecting individuals of different nondomination levels layer by layer starting from F 1 until S g ≥ N.
Assume that the last level included is the F l . Individuals in F 1 , F 2 , . . . , F l−1 have been selected for P g+1 , then individuals in F l are picked by prioritizing them with large crowding distances until P g+1 = N. The new population P g+1 are produced. Gbest of population P g+1 is the individual with the largest crowding distance in F 1 accordingly.
Repeat the above process and increment the iteration count until g = g max . The Pareto optimal solutions are gained eventually.

Entropy Weight Method-Based Evaluation Process
The entropy weight method-based evaluation process is applied to select the optimal one from the Pareto optimal solutions acquired in Part 1 of the intelligent approach, which is illustrated in Part 2 of Figure 1. Detailed steps are as follows.

Standardization
Suppose there are n opt Pareto optimal location solutions. Indicator f nm is the value of the mth objective of the nth Pareto optimal solution. First, f nm needs to be normalized as f nm . If the objective function is a positive indicator, its standardization process is in accordance with formula (20); otherwise, it uses formula (21). Due to the three objectives in this research are negative indicators, their standardization applies formula (21).
Then, the proportion of f nm in each objective is calculated by formula (22).

Calculation of Objectives' Weights and Solutions' Composite Scores
Calculation of each objective's weight depends on the entropy e m of each objective and the information entropy redundancy d m that formulated in Equations (23) and (24). e m represents the information entropy provided by the mth objective. The smaller the value of e m , the more information the mth objective provides in the evaluation process.
Then, the weight w m of each objective is calculated according to formula (25). The smaller the value of e m , the greater the value of w m . Finally, each location solution among the Pareto optimal solutions has a composite score calculated by formula (26) and the location solution with the highest score is the optimal solution to the investigated problem.

Data Collection and Processing
To obtain the travel trips of floating vehicles, this research collects and uses the GPS trajectory data of taxis in Chengdu, China. The GPS-enabled device installed in each taxi collects the vehicle's positions in latitude and longitude, and travel speed in a real-time manner (roughly once per several seconds). One-month data, from 0:00 on 1 June 2015, to 23:59 on 30 June 2015, are collected. These data consist of more than 2 billion original data records, and can reflect the real-time operation statuses of more than 10,000 taxis per day.
The raw trajectory data should be processed to obtain the trips data of each taxi to simulate travel and recharge patterns of PHEVs. In the process procedure, if the parking time of the taxi is not less than 30 min, we think this parking event is the boundary between trips. Parking events with less than 30 min of parking time between any two consecutive trips are neglected and these trips are consolidated to a large trip. Thus, trips data with parking time of not less than 30 min are generated directly from the raw taxi GPS trajectory data. Then, trips data were cleaned by deleting invalid trips caused by data recording or transmission errors and irrational trips with too long (>350 km) or too short (<1 km) distances. Finally, we get the dataset of over 700,000 valid trips in total that stored by vehicle ID and each vehicle's trips are numbered sequentially by time. The number of vehicles is 10,500, and the time span is 30 days. The data fields of trip records are described in Table 2. We suppose the end location of the current trip is the start location of the next trip. If there are CSs within a given service radius around the parking position and the remaining capacity is less than cE j , drivers will go to charge. Thus, charging demands and charging behaviors are obtained based on travel trips and the deployment of CSs. As for CS candidates, due to the construction of CSs will occupy valuable land resources, especially in the first and second-tier cities, there is no sufficient land to construct them. So, this study selected some existing gasoline stations as candidate sites for CSs. There are more than 300 gasoline stations in Chengdu currently. The top 100 sites from existing gasoline stations were selected, based on the number of parking events occurring around each station because they cover 99% of parking points. Figure 2 shows the distribution of CSs candidates.
Trip-related parameters required for the mathematical model are d jk , d ijk and t jk . d jk denotes the travel distance described in Table 2. d ijk is the distance from the parking position of the jth vehicle's the kth trip to the ith CS candidate's location. t jk represents the time of the jth vehicle arrives at CS after the kth trip, but in numerical experiments, t jk is the end time described in Table 2 for simplifying calculations.
Given the location solution y and travel trips of floating vehicles, the values of objective functions can be calculated by the method described in Section 4.1.3. Sustainability 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sustainability on travel trips and the deployment of CSs. As for CS candidates, due to the construction of CSs will occupy valuable land resources, especially in the first and second-tier cities, there is no sufficient land to construct them. So, this study selected some existing gasoline stations as candidate sites for CSs. There are more than 300 gasoline stations in Chengdu currently. The top 100 sites from existing gasoline stations were selected, based on the number of parking events occurring around each station because they cover 99% of parking points. Figure 2 shows the distribution of CSs candidates.   Table 2 for simplifying calculations. Given the location solution y and travel trips of floating vehicles, the values of objective functions can be calculated by the method described in Section 4.1.3.

Results and Discussions
This section presents the experimental results and discusses the performances of different solutions. The solutions' performances are evaluated by three indicators corresponding to the three objective functions presented in Section 3.2, including the number of CSs, daily average CO2 emissions reduction (daily ACER, unit: tons per day) and average waiting time (AWT, unit: hours per trip). The number of CSs and AWT are the values of objective functions (1) and (3). Daily ACER is the CO2 emissions reduction generated by PHEVs in the case area compared to using gasoline in all trips, which is calculated by formula (27).

Results and Discussions
This section presents the experimental results and discusses the performances of different solutions. The solutions' performances are evaluated by three indicators corresponding to the three objective functions presented in Section 3.2, including the number of CSs, daily average CO 2 emissions reduction (daily ACER, unit: tons per day) and average waiting time (AWT, unit: hours per trip). The number of CSs and AWT are the values of objective functions (1) and (3). Daily ACER is the CO 2 emissions reduction generated by PHEVs in the case area compared to using gasoline in all trips, which is calculated by formula (27).

Parameter Setting and Charging Station Location Solutions
Parameters listed in Section 3.1 other than trip-related parameters (d jk , d ijk , t jk ) need to be set in the proposal approach and experiments. Some of them are set according to the literature, including W a = 0.22 kWh/km [16], λ = 0.47 kg/kwh [34], γ = 0.22 kg/km [35]. Some of them are set based on the information provided by Alternative Fuels Data Center (AFDC) of U.S. Department of Energy [36], including W b = 0.17 L/km, t 0 = 15 min, β = 14 kWh/h, the PHEVs' fuel consumption per 100 km is 10 L and γ = 2.24 kg/L, a PHEV can travel nearly 100 miles on electricity alone and E j = 35 kWh is converted based on this information and W a . We set c = 30%, r = 2 km, s = 20 and I = 100 according to the realistic requirements of the case city. M is a sufficiently large number and we set M = 100,000.
A total of 137 Pareto optimal solutions are obtained by the improved MOPSO process, which is shown in Figure 3.
The descriptive statistics of the objective values generated by these 137 Pareto optimal solutions are described in Table 3. The number of CSs varies from 28 to 67. Daily ACER varies from 0.22 to 30.36 tons per day and the average value is 16.14. AWT varies from 0.21 to 0.58 h and it has the most concentrated distribution. A total of 137 Pareto optimal solutions are obtained by the improved MOPSO process, which is shown in Figure 3. The descriptive statistics of the objective values generated by these 137 Pareto optimal solutions are described in Table 3. The number of CSs varies from 28 to 67. Daily ACER varies from 0.22 to 30.36 tons per day and the average value is 16.14. AWT varies from 0.21 to 0.58 h and it has the most concentrated distribution. Afterwards, the final optimal solution is selected by the entropy weight method-based evaluation process introduced in Subsection 4.2. Figure 4 shows the deployment of this solution, in which red points represent locations selected to install CSs and the black points represent the candidate sites not being selected.  Afterwards, the final optimal solution is selected by the entropy weight method-based evaluation process introduced in Section 4.2. Figure 4 shows the deployment of this solution, in which red points represent locations selected to install CSs and the black points represent the candidate sites not being selected.    To verify the reliability of the results obtained by our model, the effects of different parameters   To verify the reliability of the results obtained by our model, the effects of different parameters on location results and the performance comparison of single-objective and multi-objective models are presented in Sections 5.2.2 and 5.2.3. Table 5 shows the location performances generated by different settings of parameters c, s and r. As seen in this table, the charging threshold c has relatively few effects on the number of CSs, but affect significantly the other two indicators. Drivers' charging behaviors tend to be conservative with the increases of c and they are even more reluctant to withstand the inconvenience caused by the exhaustion of electricity. The number of trips that need to be recharged also increases. Therefore, with the increases of c, the possibility of completing the trip by gasoline decreases and the values of daily ACER and AWT increase. Thus, questionnaires can be utilized to investigate the charging attitudes of drivers in Chengdu and c can be set according to their charging behavioral preferences. The CS placements at different c are illustrated in Figure 5. The location solutions at c = 10%, c = 30% and c = 50% are represented by red, green and yellow points respectively. Unselected candidate stations are represented by black points. When drivers' charging preferences cannot be judged by decision makers, c can be set as 30%.  Same as c , the number of chargers s at each CS has few effects on the number of CSs but has significant effects on the other indicators. With the increase of s under the case of a similar number of CSs, the number of vehicles that are allowed to charge at the same time evidently increases and the city's overall charging service capacity improves. Thus, daily ACER increases and AWT decreases. Under the same conditions, the performance of the solution becomes better as s increases Same as c, the number of chargers s at each CS has few effects on the number of CSs but has significant effects on the other indicators. With the increase of s under the case of a similar number of CSs, the number of vehicles that are allowed to charge at the same time evidently increases and the city's overall charging service capacity improves. Thus, daily ACER increases and AWT decreases. Under the same conditions, the performance of the solution becomes better as s increases but it is limited by land resources in reality. Figure 6 shows the CSs locations of different s. The location at s = 10, 20 and 30 are consisted by red, green and yellow points, black points represent the remaining unselected candidate stations. Same as c , the number of chargers s at each CS has few effects on the number of CSs but has significant effects on the other indicators. With the increase of s under the case of a similar number of CSs, the number of vehicles that are allowed to charge at the same time evidently increases and the city's overall charging service capacity improves. Thus, daily ACER increases and AWT decreases. Under the same conditions, the performance of the solution becomes better as s increases but it is limited by land resources in reality. Figure 6 shows the CSs locations of different s . The location at 10 20 and 30 s  ， are consisted by red, green and yellow points, black points represent the remaining unselected candidate stations. Service radius r of CSs has significant influence on all indicators. With the increases of r , the number of CSs decreases, daily ACER and AWT increase. As r increases, the service range of a CS Service radius r of CSs has significant influence on all indicators. With the increases of r, the number of CSs decreases, daily ACER and AWT increase. As r increases, the service range of a CS expands, thus the charging demands that cannot be satisfied because of the service radius originally are satisfied and trips required to complete using gasoline are reduced, but the number of vehicles to be recharged in the CS increases and the queue time increases when s unchanged. An increase in r will result in a sparse distribution of CSs, which affects drivers' charging behaviors negatively. So, r preferably does not exceed 2 km in reality. The location solutions at different r are illustrated in Figure 7, in which the location at r = 1, 2 and 3 are represented by red, green and yellow points. Unselected candidate stations are represented by black points. expands, thus the charging demands that cannot be satisfied because of the service radius originally are satisfied and trips required to complete using gasoline are reduced, but the number of vehicles to be recharged in the CS increases and the queue time increases when s unchanged. An increase in r will result in a sparse distribution of CSs, which affects drivers' charging behaviors negatively. So, r preferably does not exceed 2 km in reality. The location solutions at different r are illustrated in Figure 7, in which the location at =1 2 and 3 r ， are represented by red, green and yellow points.

Effects of Different Parameters on Location Results
Unselected candidate stations are represented by black points. This section discusses the performance differences generated by single-objective and multiobjective CS location models, to address the necessity of considering multiple objectives in CS location problems. The comparison results are shown in Table 6, in which the first row shows the

Performance Comparison of Single-Objective and Multi-Objective Models
This section discusses the performance differences generated by single-objective and multi-objective CS location models, to address the necessity of considering multiple objectives in CS location problems. The comparison results are shown in Table 6, in which the first row shows the results generated by the model with the consideration of three objectives, the second to the fourth rows show the results generated by the models with objectiv1 to objective 3 respectively.  Table 6 shows that each single-objective model can generate the optimal value of its corresponding objective function, but the performances of other two indicators are not so good. The locations of objective 1 and objective 2 have the maximal and minimal number of CSs respectively. But their values of AWT are almost equal. The solution caused by objective 3 has the best performance in terms of AWT but the performances of other indicators are not obviously superior. In general, the performance of the three objectives model is more reliable and practical, the values of its indicators are slightly lower than the value of the corresponding single-objective model and higher than the values of other two models. Consequently, applying three-objectives model is the most reliable choice when planner is unable to determine the weights of different objectives. The locations of CSs generated by different models are illustrated in Figure 8, in which the location of multi-objective model is represented by red points, and the locations of objectives 1, 2 and 3 are represented by black, green and yellow points.

Conclusions
This paper presented a multi-objective model for public charging station location problem considering the triple bottom line principle under the trend of sustainable development. The mathematical model has been established with three sustainable objectives, including minimizing the number of CSs, minimizing the total CO2 emissions, and minimizing the average waiting time of trips that need to be recharged, which are especially beneficial for local governments to achieve sustainable development goals.
An intelligent multi-objective location approach is developed to tackle the investigated problem. In this approach, an improved MOPSO process, combining a discrete PSO algorithm with a nondominated sorting technique, is proposed to seek a set of Pareto optimal solutions, and an entropy

Conclusions
This paper presented a multi-objective model for public charging station location problem considering the triple bottom line principle under the trend of sustainable development. The mathematical model has been established with three sustainable objectives, including minimizing the number of CSs, minimizing the total CO 2 emissions, and minimizing the average waiting time of trips that need to be recharged, which are especially beneficial for local governments to achieve sustainable development goals.
An intelligent multi-objective location approach is developed to tackle the investigated problem. In this approach, an improved MOPSO process, combining a discrete PSO algorithm with a non-dominated sorting technique, is proposed to seek a set of Pareto optimal solutions, and an entropy weight method-based evaluation process is put forward to pick out the final preferred solution from the set of Pareto optimal solutions. The efficiency of the proposed multi-objective model has been validated by employing a large number of real trips data of floating vehicles. Experimental comparisons were conducted to examine the impacts of three parameters (charging threshold, the number of chargers and service radius of CSs) on location results and the necessity of the proposed multi-objective model.
This research assumes that all CSs have the same capacity. Our further research will consider the location problem of CSs with different capacities. Comparing the performance of the improved MOPSO process proposed in this paper with other existing multi-objective evolutionary optimization algorithms is another future direction worthy to do.