1. Introduction and Related Work
EVs are key players in the energy transition that our society must undertake to solve one of our greatest problems today: the climate crisis. Thanks to their lack of local emission, they are a primordial and short-term solution to reduce the carbon footprint of cities and to address their increasing air quality problem. In addition, EVs can contribute to an important reduction of global
emissions and also to less dependence on fossil fuels. However, one of the main drawbacks for the widespread adoption of EVs in cities is currently the insufficient number of public CSs. For example, in Spain, a public–private study considered that the number of public charging points in 2019 was around 5000, while it estimated that by 2030, there will be between
and
million EVs in circulation [
1]. Therefore, a very large number of public CSs will have to be deployed in the coming years, and optimizing their deployment is a major challenge, to which a considerable amount of public resources are being allocated. One of the main problems being determining its optimal arrangement due to its thorough interconnection with the complexity of the traffic in the city. Due to the relevance of the subject, it has received considerable research attention in the previous few years from very different perspectives.
Some works investigate the effects of several factors on the EV energy consumption. In a very recent work, Cieslik et al. [
2] have experimentally measured the energy consumption of an EV in a real use case in different scenarios (urban, rural or motorway routes) and in different driving conditions (winter and summer periods). They also investigated different scenarios for using a residential photovoltaic system as an energy source for EV charging and household appliances. Zhang et al. [
3] thoroughly analyzed the factors influencing the EV energy consumption based on real-world driving data collected from a fleet of fifty-five electric taxis in Beijing. They also employed a machine learning approach to predict energy consumption under different driving conditions. Chłopek et al. [
4] have carried out detailed field studies of the effects of different driving patterns on the EV energy consumption, comparing them with those of ICEVs. Their results showed that, under dynamic driving conditions, the EV used
less energy than the ICEV on average, even with battery charging efficiency taken into account.
Many other works study the optimal allocation of CSs, taking into account different aggregated usage data, mainly via optimization techniques. Li et al. [
5] used data mining techniques over large scale trajectory data collected from EV taxis of Shenzhen City (China). They studied how to deploy CSs and how to assign charging points to each station in order to minimize the average travel time to the nearest CS and the average waiting time for an available charging point. Their results showed that it is better to deploy more small stations when the number of charging points is sufficiently large; but super stations are a wiser choice when this number is relatively small. However, their study is based on prerecorded taxi trajectories and they do not carry out simulations to study the effect of different choices on traffic. Dong et al. [
6] employed a genetic algorithm to find the optimal locations of CSs. They handled GPS-based ICEV travel data obtained in the Seattle metropolitan area, assuming that current ICEV activity patterns will not change when switching to EVs. GPS-enabled trajectory data were also used in the paper of Qi et al. [
7] for the location of CSs, with the objective of minimizing the CO
emissions, but they did not consider the waiting time for recharging, nor the CS capacities. Chunlin Guo et al. [
8] studied the planning of CSs for urban areas with tight land supply. Specifically, they studied an area of 12 square kilometres in the middle of four main roads in Beijing. Their method was based on the sales volume of gas stations in the area, and the CS location was decided by using an energy equivalence method. However, the relationship between fast and slow CSs remains to be resolved. The paper of Hong Gao et al. [
9] studied the deployment of CSs in a hybrid network with EV and ICEV users to maximize the service efficiency of charging facilities. Min Xu and Qiang Meng [
10] also studied the optimal deployment of CSs considering path deviations and nonlinear elastic demand, but only in highway networks. They did not study the urban environment, which is subjected to more constraints such as possible traffic congestion that may occur around CSs and long charging times if only low power chargers were available. More complex models that include several factors that determine the best location of the CS, such as economic costs, power of the electrical network, traffic efficiency and EV user satisfaction have been carried out by Weiwei et al. [
11]. Pruckner et al. [
12] developed a discrete event simulation model to evaluate the success of different expansion strategies in the future development of CS networks on a municipality scale (i.e., how many new CSs must be built in each municipality during a time period). Their model is based on empirical data about the timely utilization (number and duration of charging operations) of CSs located in each municipality of the German federal state of Bavaria. They used these data to derive empirical models for the start time and duration of charging events, as well as the popularity of CSs. They did not use any traffic model, so their approach does not directly take into account the peculiarities of traffic or the location of each station in the municipality. Islam et al. [
13] employed an optimization technique (binary lightning search algorithm) to study the optimal placement and sizing of fast CSs that minimize the combined cost induced by three factors: the energy lost by EVs when travelling to CS, the buildup CS cost, and the cost due to grid power loss induced by fundamental and harmonic currents introduced by the operation of chargers. In this line, the recent work of Deb et al. [
14] studies the optimal placement and sizing of fast and slow CSs as an optimization problem, in which the cost is determined by CS placement parameters (installation, operation and travel time cost) and penalties for violating power grid constraints. They used a hybrid evolutionary algorithm combining a chicken swarm optimization and a teaching–learning-based optimization algorithm, outperforming other standard optimization algorithms such as genetic algorithms, particle swarm optimization or differential algorithms.
Although the aforementioned works have obtained great advances for the problem of optimal allocation of CSs, they have the drawback that they do not take into account the effects of CS placement on traffic. In addition, many of them are mainly based on handling a large amount of data obtained from ICEV trajectories, so they do not take into consideration the activity pattern changes due to EVs. Another strategy to study traffic problems and their relationship to the deployment of EVs is direct simulation using microscopic models such as CA. In these models, vehicles are represented as individual entities and their behavior is simulated by rules that take into account the behavior of other vehicles found in their immediate neighborhood. Traffic models with CA go back to the works of K. Nagel, M. Schreckenberg [
15] and B. Chopard [
16,
17] and more recently to B. Jia [
18], Y. Zheng [
19] and Y. Xiang et al. [
20]. Xiang et al. [
20] have employed a CA model to study the electricity demand of a CS for different traffic flows. They found that their model captures the dynamic behaviors of the interactions between the traffic flow and the charging loads. However, their model is very limited in extension: they do not study a model of a whole city or a substantial portion of it. Their model includes only one main road, one junction and one CS. Consequently, they did not investigate the effect of the CS network on traffic flows. Other microscopic models that have been used for traffic simulation are agent-based models (also referred to as multi-agent models); for example, see the works of Chen [
21], Waraich [
22] and Viswanathan [
23]. The optimal allocation of CSs has been investigated in [
23] using simulation results from an agent-based model, with application to the charging infrastructure in the state city of Singapore, in that model vehicles had fixed routes (based on daily traffic data) from an origin to a destination. This approach offered very good results for a static planning. However, it lacks flexibility and adaptability because vehicle routes were fixed. Finally, the work of Zhai et al. [
24] is maybe the most similar to the work presented here. They present a charging demand simulation method based on an agent–CA model, but unlike the present work, they do not use a microscopic model. Each cell is defined as a 1 km-long road section, which can accommodate a large number of vehicles. The dynamic fluctuations of EV traffic and charging demands are modeled by a Monte Carlo algorithm. Their case study uses a test problem with a very simple CA consisting of only 25 cells. Although this work is very interesting, due to the coarse granularity of this model, its accuracy is not sufficient to simulate the dynamic traffic processes that can take place in a real city and to capture emergent properties of the traffic as a complex system, such as traffic jams, and their effects on traffic density and charging demand. This can be reach by using a microscopic model, as shown in the present work.
In this work, we present a cellular-automata and agent-based hybrid model to analyze the layout of CS deployment in cities through microscopic traffic simulations and a software simulation tool called SIMTRAVEL that implements our model. The tool and the source code are available online at
https://github.com/amarogs/simtravel (accessed on 11 May 2021) under GPL v3 license. Researchers can download and modify the code freely to run their own particular simulations. Our model includes EVs, ICEVs and CSs. Both ICEVs and EVs coexist, since this happens in most current cities, where the percentage of ICEVs is currently much higher than that of EVs. In addition, the combination of both types of vehicles makes our tool more applicable and wide-ranging. In fact, the user has the possibility of declaring a null percentage of ICEVs if he/she wanted to simulate a city exclusively composed of EVs. Each EV is modeled as an agent that incorporates complex behaviors, like decisions about the route to destination, when to drive to a CS, or which CS to choose. We have also applied our tool to analyze the effects on the traffic flow in a synthetic city in which both traditional ICEVs and EVs circulate, depending on whether there is only a single large centralized CS or, on the contrary, CSs are scattered in smaller units throughout the city. We have considered situations of low, medium and high traffic density and different proportions between electric and combustion vehicles. Unlike [
23], in our model, each car does not have a fixed route, but can adapt its route dynamically in response to changes. These changes can consider the current traffic situation (as studied in this work) or other inputs, such as dynamic interaction with smart city traffic management systems.
The remainder of the paper is organized as follows.
Section 2 describes our model, the simulator and the experiments that have been carried out.
Section 3 presents the results and discusses their interpretation and significance. Finally,
Section 4 presents the conclusions and future research directions.
3. Results and Analysis
First, several preliminary tests were done to check the robustness of our simulator. As the street configuration is synthetic, it was proven that, without EVs that may alter the traffic around CSs, the traffic distribution was similar for different city sizes and for the same conditions.
We performed the tests using a city size sufficiently large enough to be realistic (5 km2 with six avenues), but not too large, so that the number of cells was low enough to execute a 1000-min simulation (sufficient to reach the stationary state of the city traffic) in several minutes.
This allowed us to run ten times a complete set of tests, varying the different parameters in approximately two days on a moderate size server (32 nodes summing a total of 372 cores of architecture Intel x86-64). This size was shown to be adequate to compare and find the better station arrangement for the different parameters (amount of traffic, proportion of EVs, initial EV conditions, etc.).
We have tried three different arrangements of stations, but with the same number of chargers. The three dispositions were defined as follows:
1 central station placed on an avenue with 72 chargers.
4 evenly distributed stations on avenues having 18 chargers each.
36 evenly distributed stations placed on streets (not on avenues) having 2 chargers each.
The number of vehicles in the city is taken as a proportion of the drivable cells in the city:
Besides, the number of EVs in the city is taken as a proportion of the total number of vehicles in the city:
An initial set of simulations was first carried out to find which values of can be considered as low traffic conditions, and which others lead to the first bottlenecks, mainly around the CSs. Approximately, the thresholds were found to be = and a low = . Below these thresholds, traffic conditions were almost ideal (no bottlenecks, mean speed remains constant), thus, these results are of no interest, and they are not shown here.
In order to evaluate the deployment of CSs, we have measured the following metrics for each simulation: the mean intervals that EVs spent at each state (mainly those of queuing in the station and seeking a station), the mean speed of all the vehicles (combustion and electric ones) and the cell occupancy rates that reveal bottlenecks in specific places. Graphs containing mean values are shown with deviation bars for 10 simulations with different random initial conditions.
3.1. EV State Distribution
First, we have analyzed the distribution of states of EVs during the simulation, for the three different CS layouts and different values of the parameters. This first set of tests demonstrates the correct working of our simulator and draws some first interesting results. We analyze the time spent by EVs in their two main stages that distinguish them from a combustion vehicle: (a) the mean time spent while going to the CS, that is, the mean time that a vehicle needs to reach a station once it has detected its requirement to be recharged (denoted as EV seeking time); and, (b) the elapsed time since the vehicle reached the station until a charger was available and it started to recharge (denoted as EV queuing time).
In
Figure 3, we present three graphs with the EV seeking time for different values of
. Evidently, a low proportion (10%) of EVs means that their effects on the traffic are minimized since 90% of the vehicles are ICEVs and moving around the city without the need to recharge. Obviously, the more vehicles moving in the city, the longer the search times, as a consequence of the higher traffic density. As expected, this occurs smoothly and gradually for the distributed station cases, but a remarkable effect begins to appear for one big central station: there is a threshold in the vehicle density
(or in the total number of vehicles), for which seeking times visibly rise. This critical level supposes a total number of 7000 vehicles for 10% of EVs, and decreases progressively with increasing
. This is the primary effect pointing to traffic saturation around the big central station, which does not appear when the stations are distributed. What is more, traffic saturation is mitigated somewhat by the assumption that there are infinite waiting places inside the stations, which means that many cars are in fact not contributing to traffic jams.
It is also logical that standard deviations for the set of tests are high, as the total number of vehicles approaches the threshold, which means that the traffic is problematic and has a high dependence on the chosen (random) initial conditions. This is very visible for the cases of a central station, due to the hot traffic spot found around the single station.
The next remarkable aspect is the mean time that vehicles spend queuing inside stations.
Figure 4 contains three graphs, one for each EV density, depicting this EV queuing time for each station arrangement. For
, these times are negligible (less than one minute) in all cases. However, we notice for a larger
that queues at the stations grow progressively with more vehicles. This effect is significant when multiple small-size stations are tested, due to an unbalanced distribution of EVs and stations; i.e., some stations (which differ according to the initial conditions) suffer from saturation while others are almost empty, leading to an increase in the mean queuing time. For a high EV density (
), this problem tends to be intolerable (one hour queuing), indicating that the number of chargers is clearly insufficient when there are more than 5000 vehicles (2500 EVs).
Therefore, the weak point of the distributed configuration is the saturation of some specific stations, since the more dispersed the stations are, the more likely it is to occur. It is clear that this problem should be avoided by introducing smart routing into the behavior of EV agents. Vehicles would have information about the occupancy of each station and, perhaps, the traffic conditions, and would choose accordingly. This opens up a miscellaneous set of options in the simulator, which is intended for future work.
3.2. Global Traffic Status and Collapse Conditions
The previous subsection shows that detecting traffic thresholds is crucial to extract interesting conclusions about the arrangement of stations. The main parameter that indicates traffic saturation is the mean speed of vehicles.
Figure 5 shows the mean speed for different scenarios and the total number of cars. In general, the speed remains constant (meaning that there are no traffic jams or congestion in the simulations) until a certain threshold is reached. Once again, this is more abrupt for the case of having only one big station: here, EVs congest the surroundings of the station and cause a malfunction of the whole city.
Considering these conditions, we conclude that the critical parameters for the small city presented here are of the drivable cells occupied by vehicles and of EVs. By investigating these values in greater detail, we can draw the main conclusions, which should be predicted by urban designers to prevent these problematic situations. This is analyzed in the next subsection.
3.3. EV State Evolution
Mean values for the whole simulation give a big picture of the traffic and charging behavior, but it does not allow us to identify the reasons for an abnormal situation when approaching the threshold values that begin to collapse the city. Observing the evolution of the vehicle’s states sheds light on the reasons for this situation and its consequences. An interesting case is shown in
Figure 6 for the critical parameters
,
for two different station arrangements. We can distinguish 3 or 4 periods along the complete 1000 min simulation.
The first transient period (from 0 to approximately minute 200) obviously depends on the initial conditions. In this case, the initial state for every car is declared as “At destination” (b), meaning that after some time (defined as a Gaussian distribution), it starts moving towards a new destination, entering the state “Driving towards destination” (a). Note that the other Gaussian distribution for the initial charge assumes that each car has enough energy to start moving.
In the second transient (200-500 min. approximately), the vehicles need to recharge their batteries. Once they arrive at the station, they take available chargers. In some minute, the number of vehicles in “Charging state” (e) reaches its maximum (in our simulations, only 72 chargers are available). This means that the following vehicles entering the station necessarily remain in “At station” state (b) waiting for a free charger (that is, in mode “queuing time”).
At the end of this second transient, we can find two very different behaviors: the case of distributed stations (
Figure 6, right) enters a stationary state, but the case of a central station (
Figure 6, left) continues its evolution. Evidently, the few chargers of this layout imply that most cars are forming a queue to wait for a free charger. However, and paradoxically, a simple visual inspection seems to indicate that the case of a central station is behaving better than the other one because the number of cars in “At station” state starts to decrease (around minute 700). However, the truth is very different when a detailed examination is carried out: during the second transient (around minute 300) the number of cars going towards a station grows far more for the central station than for the distributed case. The real reason for this difference is that a traffic jam starts to appear around the central station (see
Section 3.2 and
Figure 5), so it takes a long time for cars to reach the station.
After some minutes (min. 650), this results in the battery of some EVs being drained, clearly indicating that this disposition is not practical. The rest of the curves behave accordingly due to this traffic jam effect: the number of vehicles in “Driving towards destination” and “Driving towards station” states increases steadily and simply, because many vehicles are spending more and more time in the traffic jam; the number of vehicles at a destination tends to zero because they are waiting in the other states.
Finally, note that there is an important dispersion for the different trials in
Figure 6 (left) once the traffic jam is reached, because this condition depends on the initial condition and on the specific points at which traffic starts to increase.
This is also supported by the “heat map” of the traffic during the simulation.
Figure 7 shows three visual representations of the amount of traffic in the city during three instants for the case of a central station with the same conditions as those of
Figure 6. We can observe that the area surrounding the central station (located at the middle cross) is becoming increasingly congested, leading, in the end, to a global deadlock in the city.
4. Conclusions and Future Work
In this work, we present a cellular-automata and agent-based hybrid model and a simulation tool implementing this model to analyze the layout of CSs deployment in cities. Our model and tool allow the execution of microscopic traffic simulations, tracing the behavior of individual vehicles, which include EVs, ICEVs and CSs. Each EV is modeled as an agent incorporating a complex behavior which involves decisions about the following factors: the specific destination of each vehicle, the route to follow, how to navigate there, taking into account the local traffic state, when to drive to a CS to recharge the battery, and what station to drive to.
Our study has been useful to analyze the effects of different dispositions of CSs in a city. We have found very different consequences for a large central station versus a set of smaller distributed stations.
In the case of a single central station, we found that when the vehicle density exceeds a certain critical level, a large traffic jam is formed in a wide area around the central station. It collapses the city and causes an abrupt decrease in the mean speed of vehicles, even resulting in many EVs running out of battery. We also verified the effects of this collapse directly in the ‘heat map’ of the amount of traffic, in the evolution of the number of EVs at each state and in the high increase of the time spent by EVs driving towards a CS.
This situation does not happen for the distributed CS layout: here, the city reaches a stationary state in which the traffic is fluid. However, we also found that a distributed layout conveys the following drawback: the time spent by EVs queuing in the CS (that is, waiting for a free charger) gets higher than for a single large station when the vehicle density increases. The reason is an unbalanced distribution of EVs among the different CSs, so that some of them become saturated, while other ones still hold free chargers.
Our results are useful for urban designers, planners and policy makers. First, our study shows that the best arrangement of public CSs in a city is a distributed network, which avoids the serious traffic jam introduced by a single central station, which indeed affects severely the traffic conditions in the whole city. Second, it is very important to introduce smart routing measures to balance the distribution of EVs among the different CSs in the distributed network. These measures can range from placing information signs showing the occupancy level of each CS to implementing an intelligent traffic management system that directly interacts with the EVs or with their users.
The work presented here is a first step that opens up very promising perspectives for further research directions, which we plan to undertake in the future.
Firstly, our simulator can be used to investigate other ‘what if’ scenarios, involving, for example, other patterns in the CS layout, how the location of traffic attraction areas (e.g., large shopping malls) affects the CS network and traffic, or how specific massive events (such as large sport competitions) or emergency situations can affect the charging network and the traffic.
Secondly, the results presented in this work indicate that an extension of our tool to incorporate geographic data from a real city extracted using a geographic information system (GIS) is of great interest. This extended tool can be highly valuable for predicting the effects of the layout of the CS network on urban traffic in real cities and for studying the optimal deployment of these stations. To carry out such a study of a real city, it will be important to take also into account a lot of city peculiarities as additional inputs, such as the urban planning policy, the land use, the power grid and so on.
Thirdly, this model and tool can also be extended to build a traffic simulation framework to study the effects of intelligent traffic management algorithms on EV traffic and CS network. To this end, the agents of our model should be provided with inputs coming from a smart city traffic management system, including information such as the current occupancy level of each CS, the current traffic state in each part of the city, or even predictions of intelligent algorithms about the near future states. Our tool is thus a first step towards an intelligent traffic management simulation system for both electric and combustion vehicles and for CSs in smart cities.