Charging Network Planning for Electric Bus Cities: A Case Study of Shenzhen, China

: In 2017, Shenzhen replaced all its buses with battery e-buses (electric buses) and has become the ﬁrst all-e-bus city in the world. Systematic planning of the supporting charging infrastructure for the electriﬁed bus transportation system is required. Considering the number of city e-buses and the land scarcity, large-scale bus charging stations were preferred and adopted by the city. Compared with other EVs (electric vehicles), e-buses have operational tasks and di ﬀ erent charging behavior. Since large-scale electricity-consuming stations will result in an intense burden on the power grid, it is necessary to consider both the transportation network and the power grid when planning the charging infrastructure. A cost-minimization model to jointly determine the deployment of bus charging stations and a grid connection scheme was put forward, which is essentially a three-fold assignment model. The problem was formulated as a mixed-integer second-order cone programming model, and a “No R” algorithm was proposed to improve the computational speed further. Computational studies, including a case study of Shenzhen, were implemented and the impacts of EV technology advancements on the cost and the infrastructure layout were also investigated.


100% Electric Bus City
Today, problems associated with internal combustion engine vehicles have become severe and include oil shortage and vehicle emission pollution [1,2]. PM 2.5 pollution concerns many Chinese people who are living in cities. PM 2.5, which is particle matter in air with a diameter no greater than 2.5 microns, can easily absorb toxic substances that are introduced into humans' pulmonary alveoli when we breathe them in. PM 2.5 comes from coal burning, vehicle emission, biomass burning, etc. In most cities in China, emissions from coal burning contribute the most to PM 2.5 levels. However, in big cities, such as Beijing, Shanghai, Guangzhou and Shenzhen, vehicle emission is the biggest contributor among the sources of PM 2.5, accounting for 45.0%, 29.2%, 21.7%, and 52.1% in those cities, respectively [3].
Bus transportation is an important mode of travel, and it can mitigate city traffic congestion effectively. Public transportation, including bus transportation, is also one of the main modes of travel. For example, in the central areas of Beijing, journeys on public transportation accounted for 72% among all the journeys in 2018 [4]. The electrification of public transportation vehicles can further alleviate the reliance on oil and reduce the emission levels [5,6]. Although the volume of buses is smaller than private cars, buses consume considerable amounts of fuel because of their long operation times. In addition, bus powertrain and working characteristics, such as frequent stopping and starting and the heavy loads transported, result in high fuel consumption. Diesel-powered buses readily emit more particles than gasoline-powered vehicles do, further aggravating the extent of PM 2.5 air pollution. For example, in Shenzhen, China, diesel-powered buses accounted for 0.5% of the city's vehicle fleet, but contributed 20% of vehicle emissions [7]. E-buses, powered by electricity and with zero emissions, are greatly advocated to replace traditional buses. In China, the central and local government has introduced a series of policies with massive financial subsidies to support bus electrification. Shenzhen, a city with a population of 12.5 million in Southern China, has become the first megacity that has achieved full bus electrification. By the end of 2017, the city replaced all the previous buses with 16,359 pure battery e-buses and has the largest e-bus fleet in the world [7]. Every bus running on the roads in Shenzhen is an electric bus [8]. It is estimated that an e-bus reduces energy consumption by 72.9% compared with a diesel bus. According to the Shenzhen Transport Commission, each year, the use of e-buses reduces fuel oil consumption by 345 thousand tons, and the e-bus fleet reduces emissions by 1.35 million tons of carbon dioxide locally. The overall environmental impact of e-bus use depends on the energy mix.
A large number of e-bus charging stations are to be constructed to support e-bus operation. According to the growth trend of automobiles and related policies for promoting EVs [9,10], the Chinese government has estimated that China will have more than 5 million EVs by the end of 2020 [11]. The volumes of different EV types are shown in Figure 1 (note: e-trucks include electric sanitation vehicles and logistics vehicles). As the figure shows, the number of e-buses is 200,000. With the development of EVs, the issue of the supporting charging infrastructure construction needs to be addressed urgently. The Chinese government has given priority to the charging infrastructure construction of public service vehicles, such as buses, taxies, and sanitation vehicles [11]. The charging station deployment plan in China from 2015 to 2020 is shown in Figure 2. Over these five years, China has planned to construct 3850 charging stations for e-buses to support public transportation electrification. travel. For example, in the central areas of Beijing, journeys on public transportation accounted for 72% among all the journeys in 2018 [4]. The electrification of public transportation vehicles can further alleviate the reliance on oil and reduce the emission levels [5,6]. Although the volume of buses is smaller than private cars, buses consume considerable amounts of fuel because of their long operation times. In addition, bus powertrain and working characteristics, such as frequent stopping and starting and the heavy loads transported, result in high fuel consumption. Diesel-powered buses readily emit more particles than gasoline-powered vehicles do, further aggravating the extent of PM 2.5 air pollution. For example, in Shenzhen, China, diesel-powered buses accounted for 0.5% of the city's vehicle fleet, but contributed 20% of vehicle emissions [7]. E-buses, powered by electricity and with zero emissions, are greatly advocated to replace traditional buses. In China, the central and local government has introduced a series of policies with massive financial subsidies to support bus electrification. Shenzhen, a city with a population of 12.5 million in Southern China, has become the first megacity that has achieved full bus electrification. By the end of 2017, the city replaced all the previous buses with 16,359 pure battery e-buses and has the largest e-bus fleet in the world [7]. Every bus running on the roads in Shenzhen is an electric bus [8]. It is estimated that an e-bus reduces energy consumption by 72.9% compared with a diesel bus. According to the Shenzhen Transport Commission, each year, the use of e-buses reduces fuel oil consumption by 345 thousand tons, and the e-bus fleet reduces emissions by 1.35 million tons of carbon dioxide locally. The overall environmental impact of e-bus use depends on the energy mix. A large number of e-bus charging stations are to be constructed to support e-bus operation. According to the growth trend of automobiles and related policies for promoting EVs [9,10], the Chinese government has estimated that China will have more than 5 million EVs by the end of 2020 [11]. The volumes of different EV types are shown in Figure 1 (note: e-trucks include electric sanitation vehicles and logistics vehicles). As the figure shows, the number of e-buses is 200,000. With the development of EVs, the issue of the supporting charging infrastructure construction needs to be addressed urgently. The Chinese government has given priority to the charging infrastructure construction of public service vehicles, such as buses, taxies, and sanitation vehicles [11]. The charging station deployment plan in China from 2015 to 2020 is shown in Figure 2. Over these five years, China has planned to construct 3850 charging stations for e-buses to support public transportation electrification.  travel. For example, in the central areas of Beijing, journeys on public transportation accounted for 72% among all the journeys in 2018 [4]. The electrification of public transportation vehicles can further alleviate the reliance on oil and reduce the emission levels [5,6]. Although the volume of buses is smaller than private cars, buses consume considerable amounts of fuel because of their long operation times. In addition, bus powertrain and working characteristics, such as frequent stopping and starting and the heavy loads transported, result in high fuel consumption. Diesel-powered buses readily emit more particles than gasoline-powered vehicles do, further aggravating the extent of PM 2.5 air pollution. For example, in Shenzhen, China, diesel-powered buses accounted for 0.5% of the city's vehicle fleet, but contributed 20% of vehicle emissions [7]. E-buses, powered by electricity and with zero emissions, are greatly advocated to replace traditional buses. In China, the central and local government has introduced a series of policies with massive financial subsidies to support bus electrification. Shenzhen, a city with a population of 12.5 million in Southern China, has become the first megacity that has achieved full bus electrification. By the end of 2017, the city replaced all the previous buses with 16,359 pure battery e-buses and has the largest e-bus fleet in the world [7]. Every bus running on the roads in Shenzhen is an electric bus [8]. It is estimated that an e-bus reduces energy consumption by 72.9% compared with a diesel bus. According to the Shenzhen Transport Commission, each year, the use of e-buses reduces fuel oil consumption by 345 thousand tons, and the e-bus fleet reduces emissions by 1.35 million tons of carbon dioxide locally. The overall environmental impact of e-bus use depends on the energy mix. A large number of e-bus charging stations are to be constructed to support e-bus operation. According to the growth trend of automobiles and related policies for promoting EVs [9,10], the Chinese government has estimated that China will have more than 5 million EVs by the end of 2020 [11]. The volumes of different EV types are shown in Figure 1 (note: e-trucks include electric sanitation vehicles and logistics vehicles). As the figure shows, the number of e-buses is 200,000. With the development of EVs, the issue of the supporting charging infrastructure construction needs to be addressed urgently. The Chinese government has given priority to the charging infrastructure construction of public service vehicles, such as buses, taxies, and sanitation vehicles [11]. The charging station deployment plan in China from 2015 to 2020 is shown in Figure 2. Over these five years, China has planned to construct 3850 charging stations for e-buses to support public transportation electrification.

Large-Scale E-Bus Charging Station Solution
Previously, bus charging stations were simply located in bus depots, which was a common practice and the most convenient arrangement considering the relatively small number of e-buses. However, with the advent of the all e-bus city, it would not be economic to transform every bus depot into a charging station. Moreover, the power grid will exceed its capacity in some areas when high-power charging stations connect to it, which is not permissible for grid operation. Another problem is that many depots are not owned by the government or bus companies but are temporally leased fields. Depot transformation requires multilateral negotiation. Nearby citizens are concerned that charging piles emit electromagnetic radiation when in use, which is harmful to human health. They are not willing to live in an environment with high levels of radiation. Therefore, the construction of charging facilities in bus depots is not supported by the public [12]. Due to these reasons, city planners are now seeking better locations for charging stations. The locations need not necessarily be bus depots and can be selected from anywhere in the city.
Systematic planning of the supporting charging infrastructure for the electrified bus transportation system is needed. Considering the land scarcity in cities, building large-scale integrated bus charging stations is a good solution. Shenzhen has planned to build 26 large-scale bus charging stations to serve the 16,359 e-buses, each of which will have multiple floors for bus charging, parking and maintenance. The design sketch of the first large-scale bus charging station in Shenzhen, Yueliangwan Charging Station, is shown in Figure 3 [13]. The station is now under construction and anticipated to be put to use in 2019. It has 11 floors, including 10 floors above ground and one floor underground. The site area of the station is 19,465 m 2 and the total building area is 98,478 m 2 . It has 660 e-bus parking spots and 237 chargers. Compared with charging stations located in bus depots, which only have ground floors, the design of multiple floors of large-scale charging stations highly utilizes the land resources of cities.

Large-Scale E-Bus Charging Station Solution
Previously, bus charging stations were simply located in bus depots, which was a common practice and the most convenient arrangement considering the relatively small number of e-buses. However, with the advent of the all e-bus city, it would not be economic to transform every bus depot into a charging station. Moreover, the power grid will exceed its capacity in some areas when highpower charging stations connect to it, which is not permissible for grid operation. Another problem is that many depots are not owned by the government or bus companies but are temporally leased fields. Depot transformation requires multilateral negotiation. Nearby citizens are concerned that charging piles emit electromagnetic radiation when in use, which is harmful to human health. They are not willing to live in an environment with high levels of radiation. Therefore, the construction of charging facilities in bus depots is not supported by the public [12]. Due to these reasons, city planners are now seeking better locations for charging stations. The locations need not necessarily be bus depots and can be selected from anywhere in the city.
Systematic planning of the supporting charging infrastructure for the electrified bus transportation system is needed. Considering the land scarcity in cities, building large-scale integrated bus charging stations is a good solution. Shenzhen has planned to build 26 large-scale bus charging stations to serve the 16,359 e-buses, each of which will have multiple floors for bus charging, parking and maintenance. The design sketch of the first large-scale bus charging station in Shenzhen, Yueliangwan Charging Station, is shown in Figure 3 [13]. The station is now under construction and anticipated to be put to use in 2019. It has 11 floors, including 10 floors above ground and one floor underground. The site area of the station is 19,465 m 2 and the total building area is 98,478 m 2 . It has 660 e-bus parking spots and 237 chargers. Compared with charging stations located in bus depots, which only have ground floors, the design of multiple floors of large-scale charging stations highly utilizes the land resources of cities. Bus charging stations should remain exclusive and not be shared with public charging stations for the following reasons: • Compared with other vehicle types, buses have strict operation tasks and their charging demand is rigid. The charging demand must be fulfilled at the time when needed, otherwise passenger transport tasks will be influenced, which is not permissible for bus operation. To eliminate the level of disturbance from personal EVs and taxies, bus charging stations are better to be appropriative. • In the case of public charging stations, these mainly serve small-size vehicles; therefore, the area of the charging spots is not usually large enough to accommodate buses, whose lengths can reach 9 to 12 m. • Since e-bus battery sizes are large and e-buses have operation tasks, a high charging speed is usually used to decrease charging time as the charging time takes away from the bus operation Bus charging stations should remain exclusive and not be shared with public charging stations for the following reasons: • Compared with other vehicle types, buses have strict operation tasks and their charging demand is rigid. The charging demand must be fulfilled at the time when needed, otherwise passenger transport tasks will be influenced, which is not permissible for bus operation. To eliminate the level of disturbance from personal EVs and taxies, bus charging stations are better to be appropriative.

•
In the case of public charging stations, these mainly serve small-size vehicles; therefore, the area of the charging spots is not usually large enough to accommodate buses, whose lengths can reach 9 to 12 m.

•
Since e-bus battery sizes are large and e-buses have operation tasks, a high charging speed is usually used to decrease charging time as the charging time takes away from the bus operation time.
The designed charging rate for public charging stations may not meet the charging requirement of e-buses. However, certain public charging stations that have extra capabilities to serve e-buses can be chosen as candidate bus charging stations as long as charging e-buses in the stations is differentiated from charging other social vehicles in both the charging areas and the charging service provided.
In this study, the charging stations which are planned are plug-in charging stations. Currently, there are three charging technologies: battery swapping [14,15], plug-in charging [16,17], and wireless charging [18]. Each has its advantages and disadvantages. Battery swapping was adopted early, such as for the e-buses with battery swapping in the 2008 Beijing Olympic Games. However, security concerns brought by frequent swapping and unrealized battery standardization hinder its popularization. For vehicle manufacturers, battery technology and the Battery Management System (BMS) are their core business and the source of profit. Batteries can even contribute 50% of the sale price in some EV models [19,20]. Therefore, vehicle manufacturers are reluctant to promote battery standardization because it means that this important part of the profit will be lost. For this reason, the swapping mode is not widely used. Better Place, the company once well-known for its battery swapping business model, was shut down in 2013 because few customers had battery-swapping needs. Wireless charging is the most promising technology as it allows charging while driving, and facilities for such charging have been installed on the Korea Advanced Institute of Science and Technology campus and in the Korean city of Gumi [18]. Utah State University also demonstrated dynamic wireless charging technology for an electric bus with a peak power of 25 kW at its test track in 2016 [21]. Although wireless charging is the most advanced, it is still not mature enough in terms of its technology and business model to be put into large-scale applications. Plug-in charging, which is more secure, technologically more mature, and more economic than the other two, currently dominates the EV charging market. All the e-buses in Shenzhen use plug-in charging to charge their batteries.
Compared with internal combustion engine vehicles, electric vehicles need to face the battery-life and battery replacement problems. A survey that investigated the EV battery replacement cost showed that the cost was high [22]. For one EV model, e5 450, produced by the auto manufacturer BYD, the sale price is 200 thousand yuan and the battery replacement cost is half of the price. For e-buses, the battery lifetime is quite limited because of the frequent usage of buses [23]. Currently, the battery lifetime is not consistent with the e-bus lifetime. In Shenzhen, e-buses can be in service for about 8 years [24], while their battery lifetime is shorter and is dependent on the recharge cycles. For the lithium iron phosphate batteries used in most e-buses [25], such as the e-bus model, K9, produced by BYD and currently running in Shenzhen, the batteries can recharge for about 4000 cycles [26]. According to this, the e-bus battery is estimated to have a lifetime of 4 to 5 years. Therefore, the batteries need to be replaced once during the e-bus' lifetime, which is costly [7]. The sale price of the e-buses of Shenzhen was 2 million yuan, and the subsidies from central and local government amounted to 1 million yuan per e-bus. For the remaining 1 million yuan, the battery costs were 0.35 million and the bus body costs were 0.65 million [27].

E-Bus Charging Infrastructure Planning
EV charging infrastructure planning has attracted great attention in recent years. The main research objects of the studies are general EVs. Regarding e-buses, they have distinct operation characteristics which should be considered when locating e-bus charging stations. Previous studies have investigated the charging infrastructure planning for e-buses, especially since the advent of wireless charging [18,[28][29][30]. Wireless charging facilities are installed at bus stops, which utilizes bus dwelling time at the stops, or are paved under the roads along bus routes, which can charge e-buses when buses are driving on the roads [31,32]. Since buses have fixed routes and high operation frequencies, wireless charging facilities along bus routes can be highly utilized. A path refueling location model (PRLM) locates facilities along the path to ensure that vehicles can finish their trips by considering the energy consumption on the road and the energy supply at the refueling stations [33][34][35]. The idea of PRLM has been applied in the e-bus charging facility planning based on wireless power transfer technology to optimize the location of wireless charging stops [36] or charging lanes [21,30]. Liu et al. [37] further considered the uncertainty of energy consumption and proposed a robust model to optimize the planning of bus charging stops. Bi et al. [38] proposed a multi-objective optimization framework that evaluated not only the minimization of costs but also newly incorporated the objectives of minimizing life cycle greenhouse gas emissions and energy consumption during the entire lifetime of a wireless charging bus system when deploying the wireless charging stops. Liu, Song and He [23] optimized the location of wireless charging lanes and the sizing of bus batteries. These studies did consider bus operation characteristics and buses can charge en route while they are in operation.
However, concerning plug-in charging, which this study focused on, buses cannot charge en route because plug-in charging requires bus stopping and charging takes considerable time, therefore, passengers on board would not have the patience to wait. Three plug-in charging levels are used in industry: Level 1, 2 and 3. Level 3 operates at the highest power level and is usually adopted in fast charging [39,40]. However, even under fast charging, it still requires long charging times because of the large battery capacity of e-buses. For example, the e-bus model, K9 of BYD Company, requires 3 h for a full charge under fast-charging mode and 6 h under slow-charging mode [26,41]. E-buses can only use plug-in charging service after finishing the passenger transport tasks. This study aimed at solving the plug-in charging station planning problem for e-buses according to bus operation and plug-in charging characteristics.
For other studies related to e-bus charging, charging stations are often assumed to sit in the depots or terminal stops. One category of the studies is the e-bus scheduling problem. Rogge et al. [42] analyzed bus charging scheduling and the impact on power grid. Rogge et al. [43] continued to focus on depot charging battery buses and solved the problem that covered the scheduling of battery buses, the fleet composition, and the optimization of charging infrastructure in a joint process. However, the optimization of the charging infrastructure only determined the number of chargers to be installed in the depots, which did not involve the issue of charging station location. Similarly, Ke et al. [44] proposed a model for simulating the operation and battery charging schedule of plug-in e-buses and determined the minimum construction cost of an all plug-in electric bus transportation system. Other focuses on e-buses were based on charging infrastructure comparison [45][46][47], the charging strategy interacting with the power grid [48], and the BMS [49][50][51][52].

Grid Consideration
The planning of large-scale e-bus charging stations must consider power grid factors, which is different from most facility location problems. For example, Shenzhen's e-bus fleet, with 16,359 e-buses, is the largest in the world at present, bigger than New York's, Los Angeles', New Jersey's, Chicago's and Toronto's e-bus fleets combined [7]. This huge e-bus fleet is only served by 26 charging stations, each of which is a large-scale electricity-consuming station. In China, the centralized large-scale charging stations usually connect to the exclusive transformer substations, which can avoid congestion with other loads in the distribution network. The access of large-scale electricity-consuming stations must conform to the requirements of the power grid, and on the other hand, it affects the operation of the power grid [53][54][55][56][57]. E-bus charging station construction is restricted by the power grid, such as grid capacity and grid safety [58]. Particularly with urban development, power utilization has been close to saturation in cities and the available power resources are limited. Charging puts stress on the power system, such as load surges [59][60][61]. Because the battery loaded on e-buses is large, to avoid the charging time taking too long, high charging power is usually adopted to charge e-buses rapidly. The intensive use of high-power electric equipment in centralized charging stations leads to an intense burden on the grid. Another consideration of the grid is from the cost perspective. According to our interview with charging service providers, different locations and sizing of charging stations bring different grid-related cost and the cost is high, with grid lines being particularly expensive.
Currently, much of the available literature has studied the influence of charging stations on the power grid, most of which is qualitative [16,53,54]. Only a few studies have investigated the connection to the power grid quantitatively [62,63]. However, the power grid needs to be considered quantitatively during planning, especially for large-scale charging stations and planning on the city level. In this study, the transportation system and power grid were combined and a detailed charging-station connection solution to the power grid was provided. More specifically, which charging station connects to which power grid node and how to allocate charging loads in the grid were determined.
The technical features of bus charging stations discussed in this study can be summarized as follows: (1) plug-in charging technology; (2) fast-charging mode; (3) large scale; (4) exclusive distribution network.
This study aimed at large-scale bus charging station planning, considering both the transportation network and the power grid. Besides the location and sizing of charging stations, bus charging assignment and a scheme connecting the charging stations to the power grid were also investigated. This research can help city planners construct an e-bus charging network in the city.
The main contributions of the research are listed below.
(1) The work specifically targeted large-scale plug-in charging station planning for e-buses according to bus operation characteristics. Since the operation mode of buses is significantly different from that of private cars and taxis, it was of merit to carry out the planning for bus charging stations. (2) Conduct interdisciplinary research to combine the transportation network with the power grid into the model to achieve the global optimization of the two systems. (3) Transform the model into mixed-integer second-order cone programming (MISOCP) and further propose a "No R" exact algorithm to improve the computational speed. (4) Implement a case study of the real-world transportation network in Shenzhen and study the impacts of EV technology advancements on the cost and the infrastructure layout. One major finding was that the e-bus driving range is the key factor to lower the cost of the bus charging system.
The rest of this study was organized as follows. Section 2 presents the large-scale charging station planning model for e-buses. The solution method is described in Section 3. Section 4 demonstrates the computational studies. Section 5 concludes the study and provides directions for future research.

E-Bus Charging Characteristics
The bus characteristics which should be considered in the model depend on the charging modes and facility types. In a typical bus operation, buses depart from depots and reach the initial stop, followed by sequential stops, and then reach the final stop. For wireless charging facility planning, bus routes and locations of stops must be considered. For plug-in charging, as discussed in this study, since charging requires extra operation and takes a considerable amount of time, passengers onboard would lack the patience to wait. Therefore, buses can only be charged after completing the trips. Namely, they cannot be charged en route while they are in operation, which is the main difference between plug-in charging and wireless charging behavior for e-buses. Three sites were considered in the plug-in charging station planning problem-bus depot, final stop, and initial stop (the initial stop becomes the final stop in the reverse direction)-which are the candidate origins of bus charging trips. Buses can travel directly from the bus depots to charging stations or from the terminal stop after finishing the passenger transport trips. The characteristic that charging trips always start from fixed points is unique for e-buses. This characteristic does not apply to other electric vehicles, such as taxies and private cars, because they usually do not start the charging trips from fixed points. This study considered the bus operational characteristics and plug-in charging behavior in the model.
There are two charging strategies that e-buses can adopt. The first one is to charge them if they are not in operation, even if the battery power is still high. The second strategy is to charge them when the power is low. Since most large-scale bus charging stations are not located at bus depots, it takes time to travel to charging stations. Bus companies usually adopt the second strategy to decrease charging trips, and buses undergo a full charge each time they go to charging stations. Table 1 summarized the notations used in the study, which included the sets, parameters and decision variables. The parameters in the table were exogenous, some of which will be explained later. the demand for chargers of bus line i ∈ I f j annual fixed costs of location j ∈ J, which could be the construction cost of a new station that includes the price of land, or the reconstruction cost of original bus depots c annual unit cost of chargers u driving range of an electric bus w i π j distance from π of bus line i to location j b transportation cost per kilometer, which is measured by the product of electricity consumed per kilometer and the market price of electricity α jm annual line construction cost to connect location j to power grid node m. H charging speed/the power of charger s l m original loads of node m s max the maximal loads of the power grid z mn z mn = r mn + ig mn denotes the impedance of branch mn where r mn and g mn denote the resistance and reactance, respectively t mn annual power loss time of branch mn β unit price of electricity l mn0 the original square of the magnitude of the complex current from node m to node n before charging stations connect to the power grid v m , v m voltage limit of node m, v is the square of the magnitude of the complex voltage l mn current limit of branch (m, n), l is the square of the magnitude of the complex current Decision variables S mn S mn = P mn + iQ mn denotes the sending-end power flow from node m to node n, where P mn and Q mn denote the real and reactive power flow, respectively s m s m = p m + iq m denotes the power injection of node m where p m and q m denote the real and reactive power injection, respectively v m the square of the voltage magnitude of node m l mn the square of the current magnitude of branch (m, n) The parameter h i refers to the total number of charging trips of e-buses in bus line i ∈ I per year. The measurement of this parameter depends on several factors: the driving range of e-buses, the daily operating range of passenger transport, and the safety range.
According to the specification, there are at least three types of e-bus driving range: constant-speed range, range under city conditions, and range under working conditions [26,41]. The constant-speed range is measured when the bus runs at a constant speed, which is an ideal condition. The range under city conditions is measured when the bus runs in the city, and therefore, faces traffic congestion and other city traffic conditions. The last one, the range under bus working conditions, is measured under a real working environment, such as frequent stopping and starting. Passenger volume, traffic congestion in the city, and other working conditions are also considered. For example, e-bus model K9 of BYD Company has a driving range of more than 500 km under 40 km/h constant driving speed, more than 300 km under city conditions, and 250 km under working conditions [26,41]. The driving range in this study refers to the range under working conditions.
For the safety range, buses do not fully discharge the batteries when operating and the batteries still have some power remaining before charging. This practice can protect batteries, improving their longevity. Another advantage of setting the safety range is that it can act as a buffer against the fluctuation of the operating range of passenger transport. For e-buses, the drivers estimate the power the next trip will consume before they embark on it. If the remaining power cannot support the next trip, such as 20 km for a single trip or 40 km for a round trip, they do not execute the passenger transport task. Setting the safety range can also avoid a situation in which buses stop halfway because of an unexpected sharp decline in battery power, which is strictly not allowed for customer service. Bus companies set a relatively wide safety range for e-bus operation. For example, e-buses serving Line 245 in Shenzhen do not operate passenger transport tasks when the remaining SOC (state of charge) is below 30% [64]. An SOC of 30% can support a 75-km journey. The distance associated with charging trips is included in the safety range. Although the exact distances of charging trips are unknown before the issue of assigning e-buses to a specific charging station is solved, the distance is usually not long because buses are assigned to nearby charging stations to save transportation costs generated by charging trips.
Assume e-buses go for charging every X day(s), The number of charging trips of e-buses in bus line i per day is λ i /X on average, therefore, Charging time y is mainly decided by the battery capacity C and the charging power H, In China, the government usually sets a car-to-pile ratio [11,65], which suggests the number of cars that one charging pile serves per day on average, such as 5:1. The ratio is set according to the charger working time U per day and charging time y required for a full battery charge. Suppose that a charger works 15 h per day, and a full charge requires 3 h, then the ratio = U/y = 5. When charging technology improves, the ratio adjusts with the shorter charging time. The two terms, charging pile and charger, mean the same in the study.
The demand for chargers of bus line i, For example, the driving range of e-buses is 220 km, the safety range is 50 km, and the daily operating range is 170 km. According to Formula (1), buses go for charging once per day, X = 1.
Suppose that a bus line i is run by λ i = 15 e-buses. The number of charging trips of the bus line per year is h i = 15/1 × 365 = 5475. The battery capacity of e-buses that currently run in Shenzhen is 324 kW·h and it takes 3 h to finish a full charge under fast-charging mode. Assume the chargers work U = 15 h per day. The demand for chargers for e-buses in this case is d i = (λ i y)/(XU) = (15 × 3)/(1 × 15) = 3.

The Planning Model
In addition to the problem of planning the location of charging stations, the sizing of facilities is also considered, namely, the number of chargers to be installed in each charging station. This depends on the number of buses that are assigned to the stations for charging. As stated in Section 1, large-scale bus charging stations are constrained by the power grid and, on the other hand, they affect grid operation. Regarding the power grid, two important factors were considered, grid line construction cost and power loss of the grid. To be specific, these include the cost of connecting charging stations to the power grid and the increased power loss after the connection.
The planning of plug-in bus charging stations involves three layers. The first layer is to decide where buses journey from to reach charging stations. The second layer is to choose from the candidate locations of charging stations and design their capacities. The sizing of charging stations simultaneously dictates the bus assignment. If more buses are assigned to the charging station, this means that the capacity will be larger. The third layer is to determine the connection scheme of charging stations to the power grid. The planning includes several optimization problems and is interdisciplinary. They can be solved separately; however, this cannot guarantee the system runs optimally most of the time. To achieve a global optimum, these three layers would be integrated into one optimization model. Essentially, it is a threefold assignment problem: (1) which points are assigned as the origins of bus charging trips; (2) which charging stations e-buses are assigned to; (3) which power grid nodes charging stations are assigned to. These three assignment issues are closely correlated.
(1) Objective Function The planning objective is to minimize the cost of facilities, transportation cost and the cost of the power grid. Min where the first and second terms, respectively, account for the fixed cost of charging stations and the variable cost of the chargers. The third term represents the transportation cost generated by charging trips. The coefficient "2" refers to the round trips for charging. The last two terms together account for the costs associated with the power grid, which are grid line construction cost and newly increased power loss cost after charging stations connect to the grid, respectively. In the objective function, the weighting factor of each term is 1. The factors can also be changed as needed. For example, the weighting factor of the third term, the transportation cost generated by charging trips, can be increased if the time cost is considered additionally. Charging trips not only generate transportation cost, but also occupy bus operation time, which affects bus companies' benefit. Increasing the weighting factor increases the charging trip cost. The weighting factor of the last term, power loss cost, can also be increased with more comprehensive consideration of grid operation. The power loads brought by the large-scale charging stations influence grid operation, which causes not only power loss, but also voltage deviation, load surges, and so on. They are all related to grid operation. Power loss can be quantified by cost, while other impacts are not easy to quantify. Increasing the weighting factor of the last term can include other impacts on the power grid. Other weighting factors can also be adjusted to meet the market requirement. The triangles in Figure 4 represent candidate charging stations in the city, and the other three shapes represent the candidate origins of charging trips, depots, initial stops and final stops. Some bus depots are also candidates for charging stations (where the triangles and the stars overlap in Figure 4). Take Bus Line 1 as an example. D 1 , S 1 and E 1 are the depot, initial stop and final stop of Line 1, respectively. The buses on Line 1 start from one of the origins, then travel to one of the candidate charging stations.
Constraint (6) suggests that buses on bus line i select one of the three candidate sites as the origin of the charging trips. Constraint (7) suggests that buses on bus line i are assigned to one charging station. Constraint (8) requires that if buses are assigned to a candidate charging station j, j must be built.  The triangles in Figure 4 represent candidate charging stations in the city, and the other three shapes represent the candidate origins of charging trips, depots, initial stops and final stops. Some bus depots are also candidates for charging stations (where the triangles and the stars overlap in Figure 4). Take Bus Line 1 as an example. D1, S1 and E1 are the depot, initial stop and final stop of Line 1, respectively. The buses on Line 1 start from one of the origins, then travel to one of the candidate charging stations.
≤ , ∀ , Constraint (6) suggests that buses on bus line select one of the three candidate sites as the origin of the charging trips. Constraint (7) suggests that buses on bus line are assigned to one charging station. Constraint (8) requires that if buses are assigned to a candidate charging station , must be built.
(b) Power network constraints Let = ( , ) denote a directed graph that represents a distribution network. Each node in represents a bus and each link in represents a line connecting the buses (the term bus in the power grid area is different from the vehicle bus in transportation). The distribution network is usually radial and the structure of graph will be a tree. The root of the tree, node 1, is a substation bus, which connects to the transmission network. The substation voltage is given and fixed. Each link connects an ordered pair of nodes, such as ( , ), where node lies on the unique path from node to node 1, as demonstrated in Figure 5. ( , ), and → are used interchangeably to refer to a link. For a complex number ∈ ℂ, denotes its conjugate. (b) Power network constraints Let G = (N, E) denote a directed graph that represents a distribution network. Each node in N represents a bus and each link in E represents a line connecting the buses (the term bus in the power grid area is different from the vehicle bus in transportation). The distribution network is usually radial and the structure of graph G will be a tree. The root of the tree, node 1, is a substation bus, which connects to the transmission network. The substation voltage is given and fixed. Each link connects an ordered pair of nodes, such as (m, n), where node n lies on the unique path from node m to node 1, as demonstrated in Figure 5. (m, n), mn and m → n are used interchangeably to refer to a link. For a complex number a ∈ C, a denotes its conjugate.  The triangles in Figure 4 represent candidate charging stations in the city, and the other three shapes represent the candidate origins of charging trips, depots, initial stops and final stops. Some bus depots are also candidates for charging stations (where the triangles and the stars overlap in Figure 4). Take Bus Line 1 as an example. D1, S1 and E1 are the depot, initial stop and final stop of Line 1, respectively. The buses on Line 1 start from one of the origins, then travel to one of the candidate charging stations.
≤ , ∀ , Constraint (6) suggests that buses on bus line select one of the three candidate sites as the origin of the charging trips. Constraint (7) suggests that buses on bus line are assigned to one charging station. Constraint (8) requires that if buses are assigned to a candidate charging station , must be built.
(b) Power network constraints Let = ( , ) denote a directed graph that represents a distribution network. Each node in represents a bus and each link in represents a line connecting the buses (the term bus in the power grid area is different from the vehicle bus in transportation). The distribution network is usually radial and the structure of graph will be a tree. The root of the tree, node 1, is a substation bus, which connects to the transmission network. The substation voltage is given and fixed. Each link connects an ordered pair of nodes, such as ( , ), where node lies on the unique path from node to node 1, as demonstrated in Figure 5. ( , ), and → are used interchangeably to refer to a link. For a complex number ∈ ℂ, denotes its conjugate. Given graph , the substation voltage and the impedance , the other variables in the power network are described by the branch flow equations as follows [66,67]. Given graph G, the substation voltage and the impedance z mn , the other variables in the power network are described by the branch flow equations as follows [66,67].
Constraints (9) to (12) are correlated and together calculate the voltages and power injections of grid nodes, currents and power flows of grid branches. In Constraint (9), k → m refers to the grid branch km, where node m lies on the unique path from node k to node 1. k : k → m refers to a node of k that conforms to the definition of k → m . For example, in the power network of Figure 6, suppose that node m is node 3, then node k would be node 11 and node 8. In Constraint (10), for branch k1, node k would be node 2, node 3 and node 4. In Constraint (11), "Re" means the real part of a complex number, and |·| means the modulus of a complex number. Constraint (12) indicates the relationship between power flows, currents and voltages.
(c) Coupled constraints of the two networks Selected charging stations connect to the power grid to obtain a power supply, as demonstrated in Figure 6. Besides the original loads, some grid nodes would, therefore, have more loads.
Constraints (9) to (12) are correlated and together calculate the voltages and power injections of grid nodes, currents and power flows of grid branches. In Constraint (9), → refers to the grid branch , where node lies on the unique path from node to node 1. : → refers to a node of that conforms to the definition of → . For example, in the power network of Figure 6, suppose that node is node 3, then node would be node 11 and node 8. In Constraint (10), for branch k1, node would be node 2, node 3 and node 4. In Constraint (11), "Re" means the real part of a complex number, and | • | means the modulus of a complex number. Constraint (12) indicates the relationship between power flows, currents and voltages.
(c) Coupled constraints of the two networks Selected charging stations connect to the power grid to obtain a power supply, as demonstrated in Figure 6. Besides the original loads, some grid nodes would, therefore, have more loads.
Constraint (13) suggests that a chosen charging station should connect to the power grid. Assume that charging piles' power factors are 1. Constraint (14) calculates the total power injection 0 ≤ l mn ≤ l mn , ∀(m, n) Constraint (13) suggests that a chosen charging station should connect to the power grid. Assume that charging piles' power factors are 1. Constraint (14) calculates the total power injection of each node after the charging stations connect to the grid. Constraint (15) indicates that the total loads cannot be greater than the grid's capacity. Constraint (16) suggests the fluctuation range of the voltage. Constraint (17) requires that the current does not exceed the limit.

Mixed-Integer Second-Order Cone Programming (MISOCP)
During the last two decades, MISOCP has been used to solve certain mixed-integer optimization problems with its high computational efficiency [68]. The planning model can be transformed into the MISOCP after some mathematical treatments. The nonlinear terms in the model were processed, including the division Constraint (12), the variable product of R πi Y ij in the objective function, and Y ij Φ jm in Constraint (14). Other terms in the objective function and constraints were linear with regard to the decision variables. Constraint (11) was also transformed into the form of real numbers in this subsection.
(1) Convexifying Constraint (12) Farivar and Low [66] showed that Constraint (12) in radial networks admits exact second-order cone relaxation, which relaxes the equality constraint to an inequality constraint, that is, "Exact" means that the constraint holds as an equality at the optimal solution, although it is relaxed. Hence, the relaxation does not influence the optimal solution of the original problem. S mn is the complex power, |S mn | 2 = |P mn + iQ mn | 2 = P mn 2 + Q mn 2 (23) Therefore, which can be equivalently transformed into the standard second-order cone constraint, as follows: For the detailed mathematical transformation, please refer to Part 1 of Supplementary Material 1.
(2) R πi Y ij and Y ij Φ jm To treat the product of variables, R πi Y ij and Y ij Φ jm , two auxiliary variables were introduced and several linear constraints were added. For R πi Y ij , let E πij replace R πi Y ij and add the four constraints.
When R πi , Y ij = (0, 0) or (0, 1) or (1, 0), the four constraints will limit E πij = 0, namely R πi Y ij will be equal to 0. The transformation is an equivalent transformation. For Y ij Φ jm , let W ijm replace Y ij Φ jm and add the constraints below.
v m − v n = 2Re(r mn P mn + g mn Q mn − ig mn P mn + ir mn Q mn ) − r mn 2 + g mn 2 2 l mn (36) v m − v n = 2(r mn P mn + g mn Q mn ) − r mn 2 + g mn 2 l mn (37) Thus, the original model is equivalent to In the next subsection, two algorithms based on MISOCP were further proposed for the model, which would help improve the computational speed.

Improved MISOCP
(1) Relaxation of R The first algorithm, or technique, to improve the computational speed is to relax R πi ∈ {0, 1} to 0 ≤ R πi ≤ 1 and this is called Relaxation of R for short. When R πi is relaxed, the product of R πi Y ij in the original model can, coincidently, be treated in the same way as that in Section 3.1. Constraints (26) to (29) can be applied to linearize the product of a binary variable and a continuous variable ranging from 0 to 1. The relaxation of R πi does not influence the MISOCP form of the model and can increase the computational speed significantly.
Generally, R πi only obtains the value of either zero or one, even though it is relaxed to be continuous from zero to one. This is because Y ij ∈ {0, 1} assigns buses of bus line i to a certain charging station j, and buses choose the origin which is closest to the assigned charging station to minimize the travel cost. Since distances from the candidate origin to the charging station, w i π j , are different, buses all select the one with the shortest distance rather than partially selecting and partially not selecting. Therefore, R πi will either be zero or one.
However, it is possible that some of the distances are equal and R πi can possess a value in the range 0 < R πi < 1 in this situation. There are multiple optimal solutions. If w i 1 j , w i 2 j , and w i 3 j are equal, R 1i , R 2i , and R 3i can have multiple values, as long as their summation is 1 and each one of them is valued from 0 to 1. For example, If two of the three distances are equal, such as w i 1 j = w i 2 j , and they are shorter than the third distance, R 3i will be 0 and the values of R 1i and R 2i can have multiple combinations in the optimal solutions. When the values of R πi are not binary in the solution, further comparison of the selected origins (R πi > 0 means π is selected, partially or all) is performed and one of them is reselected as the origin, according to additional preferences. Besides improving the computational speed, the advantage of the relaxation is that it reveals the information of other optimal origins. When R πi is forced to be binary, the information will be lost and the reselection with additional preferences cannot be operated.
If it is hoped that R πi is valued at either zero or one after relaxation, one method is to slightly tune some distance parameters, w i π j . The situation, 0 < R πi < 1 only occurs when the three distances from different origins of bus line i to the charging station j are equal, or two of the three distances are equal and shorter than the third one. Such origins with equal distances can be partially selected in the optimal solutions. The method is to compare these origins with additional preferences in advance and slightly increase the distance value(s) of the undesirable one(s) when setting the parameters. This ensures that R πi will obtain binary values and does not influence other variables and the optimal values.
(2) "No R" algorithm The second algorithm to improve computational speed is to preprocess the three-dimensional distance matrix w i π j to a two-dimensional matrix w ij . For a certain bus line i and candidate charging station j, buses of bus line i will always select the origin which is closest to the destination j. w ij is the shortest distance between such origin and charging station j. R πi is the decision variable about the selection of the charging trip origin. In this algorithm, because the origin of the charging trips is preselected according to the trip distances, there is no R πi and the algorithm is called "No R" for short.
First, compare w i 1 j , w i 2 j , and w i 3 j and then input the smallest one to w ij . In dimensionality reduction, the information of π is lost. When forming the matrix w ij , it is necessary to record π, because for the same i but different j, π is probably different. After the assignment of i to j is solved, retrieve the corresponding π of i. The formation of matrix w ij and retrieval of π can be realized by min( ) and find( ) functions within seconds by Matlab. This is equivalent to help the model solve the first-layer assignment problem and significantly reduces the computational loads.
For the "No R" algorithm, the MISOCP model in Section 3.1 was modified. The "No R" model is as follows: Min s.t. Since there is no variable R πi , there are no Constraints (6), (19), and (26)- (29) in the "No R" model. Figure 7 shows the flow chart of the "No R" algorithm. Since there is no variable , there are no Constraints (6), (19), and (26)-(29) in the "No R" model. Figure 7 shows the flow chart of the "No R" algorithm. (

3) Computational speed comparison
To test the performance of the algorithms, a series of instances of different problem scales was designed and their computational time was compared, as shown in Table 2. LN refers to the number of bus lines, SN refers to the number of candidate charging stations, and GN refers to the number of grid nodes. These three parameters determined the problem scales. For LN = , the candidate origins would be 3 . All instances were run on IBMX3850M2 with 16 2.40 GHz Intel(R) Xeon(R) CPUs and 110 GB memory. They were solved using CPLEX 12.8, called by MATLAB R2016b and YALMIP, running on CentOS release 6.5. All the solutions were exact optimal solutions. In Table 2, ∈ 0,1 used the original MISOCP in Section 3.1. As shown in the table, the improved MISOCPs, "Relaxation of R" and "No R", performed much better than the original MISOCP. For the improved MISOCPs, the "No R" algorithm performed much better and could solve large-scale instances. In the following section of computational studies, the "No R" algorithm was adopted for computation.  (

3) Computational speed comparison
To test the performance of the algorithms, a series of instances of different problem scales was designed and their computational time was compared, as shown in Table 2. LN refers to the number of bus lines, SN refers to the number of candidate charging stations, and GN refers to the number of grid nodes. These three parameters determined the problem scales. For LN = x, the candidate origins would be 3x. All instances were run on IBMX3850M2 with 16 2.40 GHz Intel(R) Xeon(R) CPUs and 110 GB memory. They were solved using CPLEX 12.8, called by MATLAB R2016b and YALMIP, running on CentOS release 6.5. All the solutions were exact optimal solutions. In Table 2, R πi ∈ {0, 1} used the original MISOCP in Section 3.1. As shown in the table, the improved MISOCPs, "Relaxation of R" and "No R", performed much better than the original MISOCP. For the improved MISOCPs, the "No R" algorithm performed much better and could solve large-scale instances. In the following section of computational studies, the "No R" algorithm was adopted for computation.

Computational Studies
In this section, numerical experiments and a case study of Shenzhen were implemented to demonstrate the effectiveness of the joint planning model and the solution method. The impacts of EV technology advancements were also studied. The solutions attained in the computational studies were global optimal solutions.

Numerical Experiments
This subsection demonstrated an e-bus charging station planning problem where 100 e-bus lines ran in the city, hence there were 300 candidate origins of charging trips (each bus line has three candidate origins of charging trips, bus depot, final stop and initial stop). It was assumed that there were 20 candidate charging stations for selection. Figure 8 shows the distribution of the 100-100-100-20 points in the transportation network. Since some bus lines shared the same bus depots, final stops or initial stops, some of the 300 candidate origins overlapped in the figure. In China, the centralized large-scale charging stations usually connect to the exclusive 10-kV transformer substations in practice. In other words, the charging station occupies an entire 10-kV distribution network, which can avoid congestion with other loads in the 10-kV network. Therefore, the problem of the access nodes and access loads in the power grid is not an issue in the 10-kV voltage class but does matter in the upper voltage class. The 110-kV distribution network of Zhang [62] was adopted as the upper voltage class. The e-bus parameters referred to the real operating buses in Shenzhen, China, where the driving range of e-buses is 250 km; the battery capacity is 324 kW·h [26,69]. The charging stations provided plug-in charging service for e-buses and adopted a fast-charging mode with a charging speed of 108 kW. A full battery charge for e-buses took around 3 h. For other parameter settings, please refer to Part 2 of Supplementary Material 1. In this section, numerical experiments and a case study of Shenzhen were implemented to demonstrate the effectiveness of the joint planning model and the solution method. The impacts of EV technology advancements were also studied. The solutions attained in the computational studies were global optimal solutions.

Numerical Experiments
This subsection demonstrated an e-bus charging station planning problem where 100 e-bus lines ran in the city, hence there were 300 candidate origins of charging trips (each bus line has three candidate origins of charging trips, bus depot, final stop and initial stop). It was assumed that there were 20 candidate charging stations for selection. Figure 8 shows the distribution of the 100-100-100-20 points in the transportation network. Since some bus lines shared the same bus depots, final stops or initial stops, some of the 300 candidate origins overlapped in the figure. In China, the centralized large-scale charging stations usually connect to the exclusive 10-kV transformer substations in practice. In other words, the charging station occupies an entire 10-kV distribution network, which can avoid congestion with other loads in the 10-kV network. Therefore, the problem of the access nodes and access loads in the power grid is not an issue in the 10-kV voltage class but does matter in the upper voltage class. The 110-kV distribution network of Zhang [62] was adopted as the upper voltage class. The e-bus parameters referred to the real operating buses in Shenzhen, China, where the driving range of e-buses is 250 km; the battery capacity is 324 kW·h [26,69]. The charging stations provided plug-in charging service for e-buses and adopted a fast-charging mode with a charging speed of 108 kW. A full battery charge for e-buses took around 3 h. For other parameter settings, please refer to Part 2 of Supplementary Material 1.

Joint Optimization
The joint optimization model, which combined the transportation network with the power grid, was applied to solve the e-bus charging station planning problem.
The planning result is shown in Figure 9. Charging station Nos. 9, 11 and 13 were selected and connected to grid node Nos. 2, 2 and 14, respectively. A total of 127, 151 and 67 chargers were installed in the charging stations, respectively (note that the number of chargers in the charging stations may not be integers in the optimal solutions. In that case, the amounts will be rounded up to integers). The bus charging trips are shown in green lines. Table 3 shows the detailed bus charging assignment

Joint Optimization
The joint optimization model, which combined the transportation network with the power grid, was applied to solve the e-bus charging station planning problem.
The planning result is shown in Figure 9. Charging station Nos. 9, 11 and 13 were selected and connected to grid node Nos. 2, 2 and 14, respectively. A total of 127, 151 and 67 chargers were installed in the charging stations, respectively (note that the number of chargers in the charging stations may not be integers in the optimal solutions. In that case, the amounts will be rounded up to integers). The bus charging trips are shown in green lines. Table 3 shows the detailed bus charging assignment scheme of the first 20 bus lines. In the table, Origin 1 stands for bus depot, 2 for the initial stop, 3 for the final stop, and CS stands for charging station. For example, e-buses of bus line 3 traveled from the bus depot to charging station No. 11 for charging. For the complete scheme of the 100 bus lines, please refer to Part 2.2 in Supplementary Material 1.    1  2  11  11  3  11  2  2  13  12  2  11  3  1  11  13  2  11  4  1  11  14  1  11  5  3  11  15  1  13  6  3  11  16  1  9  7  3  11  17  3  9  8  3  11  18  3  9  9  3  11  19  2  11  10  3  11  20  1  9 4.1.2. Joint Optimization vs. Separate Optimization The transportation system and power grid were jointly optimized in the planning model. In this subsection, optimizations of each were conducted, i.e., the two systems were optimized separately. For the individual optimizations, two optimization methods are used, transportation-first   1  2  11  11  3  11  2  2  13  12  2  11  3  1  11  13  2  11  4  1  11  14  1  11  5  3  11  15  1  13  6  3  11  16  1  9  7  3  11  17  3  9  8  3  11  18  3  9  9  3  11  19  2  11  10  3  11  20  1  9 4.1.2. Joint Optimization vs. Separate Optimization The transportation system and power grid were jointly optimized in the planning model. In this subsection, optimizations of each were conducted, i.e., the two systems were optimized separately. For the individual optimizations, two optimization methods are used, transportation-first optimization and grid-first optimization. Which optimization method to be used depends on which system it prioritizes. optimization and grid-first optimization. Which optimization method to be used depends on which system it prioritizes.
Step 1. Focus on the transportation network to optimize the location and sizing of charging stations and the bus charging assignment. The objective function is Part 1. The results acquired in this step are used as inputs in the next step.
Step 2. Optimize the power system to determine the connection of charging stations to power grid nodes. The objective function is Part 2 in the second step.
Step 1. Use Part 2 as the objective function to optimize the power system.
Step 2. Calculate Part 1. The results of the joint optimization are shown in the previous subsection. Figure 10 and Figure  11 show the results of transportation-first optimization and grid-first optimization, respectively. As shown in Figure 10, to reduce the transportation cost generated by charging trips, as many as nine charging stations were built. Based on the planning result in Step 1, Step 2 optimized the grid access scheme. For the grid-first optimization shown in Figure 11, only one charging station was built to minimize grid access cost, and all the e-buses were assigned to it for charging.
Step 1. Focus on the transportation network to optimize the location and sizing of charging stations and the bus charging assignment. The objective function is Part 1. The results acquired in this step are used as inputs in the next step.
Step 2. Optimize the power system to determine the connection of charging stations to power grid nodes. The objective function is Part 2 in the second step.
Step 1. Use Part 2 as the objective function to optimize the power system.
The results of the joint optimization are shown in the previous subsection. Figures 10 and 11 show the results of transportation-first optimization and grid-first optimization, respectively. As shown in Figure 10, to reduce the transportation cost generated by charging trips, as many as nine charging stations were built. Based on the planning result in Step 1, Step 2 optimized the grid access scheme. For the grid-first optimization shown in Figure 11, only one charging station was built to minimize grid access cost, and all the e-buses were assigned to it for charging. optimization and grid-first optimization. Which optimization method to be used depends on which system it prioritizes.
Step 1. Focus on the transportation network to optimize the location and sizing of charging stations and the bus charging assignment. The objective function is Part 1. The results acquired in this step are used as inputs in the next step.
Step 2. Optimize the power system to determine the connection of charging stations to power grid nodes. The objective function is Part 2 in the second step.
Step 1. Use Part 2 as the objective function to optimize the power system.
Step 2. Calculate Part 1. The results of the joint optimization are shown in the previous subsection. Figure 10 and Figure  11 show the results of transportation-first optimization and grid-first optimization, respectively. As shown in Figure 10, to reduce the transportation cost generated by charging trips, as many as nine charging stations were built. Based on the planning result in Step 1, Step 2 optimized the grid access scheme. For the grid-first optimization shown in Figure 11, only one charging station was built to minimize grid access cost, and all the e-buses were assigned to it for charging.

Comparison
The cost comparison of the three optimization methods is shown in Figure 12. The joint optimization method balanced the costs generated in Part 1 and 2 and achieved the minimal total cost of 8.41 million yuan. Transportation-first optimization and grid-first optimization achieved minimal cost for Part 1 and Part 2, respectively. However, the cost for the other part was huge, especially for the grid-first optimization. Joint optimization, by considering transportation and power network simultaneously, was the most cost-efficient, which saved 18% of the total cost compared with transportation-first optimization, and 33% compared with grid-first optimization.

Comparison
The cost comparison of the three optimization methods is shown in Figure 12. The joint optimization method balanced the costs generated in Part 1 and 2 and achieved the minimal total cost of 8.41 million yuan. Transportation-first optimization and grid-first optimization achieved minimal cost for Part 1 and Part 2, respectively. However, the cost for the other part was huge, especially for the grid-first optimization. Joint optimization, by considering transportation and power network simultaneously, was the most cost-efficient, which saved 18% of the total cost compared with transportation-first optimization, and 33% compared with grid-first optimization.

Comparison
The cost comparison of the three optimization methods is shown in Figure 12. The joint optimization method balanced the costs generated in Part 1 and 2 and achieved the minimal total cost of 8.41 million yuan. Transportation-first optimization and grid-first optimization achieved minimal cost for Part 1 and Part 2, respectively. However, the cost for the other part was huge, especially for the grid-first optimization. Joint optimization, by considering transportation and power network simultaneously, was the most cost-efficient, which saved 18% of the total cost compared with transportation-first optimization, and 33% compared with grid-first optimization.  The joint optimization model breaks the boundary between transportation and power systems and integrates them to achieve global optimization of the two systems, which the separate optimization methods cannot guarantee. The set of numerical experiments justified the effectiveness of this model.

Technological Considerations
EV technology will continue to develop in the near future, such as charging technology and driving range. EV technology influences the cost of e-bus charging infrastructure and the infrastructure layout. In this subsection, the impacts of the advancements in charging speed and driving range on the cost and the layout were studied.

Charging Speed
The current charging power of EV chargers in Shenzhen is 108 kW. In China, the maximal EV charging current can reach up to 1000 A and the charging power reaches 750 kW. In this subsection, the charging speed was increased from 1 time (108 kW) to 10 times (1080 kW) to see the cost changes and examine the layout robustness to progress in charging technology. Table 4 summarizes the cost and layout change with the increase in the charging speed. The cost of Term 2 in the objective function is the charger cost, which is directly related to the number of chargers, and the amount is decided by the charging power. Since the change of charging power only affects the second term of the objective function, the costs of other terms remain the same and are omitted in the table. The parameter settings of the first instance in the table were the same as in Section 4.1. Apart from the charging power, the other parameters of the ten instances were the same. For the layout, charging station Nos. 9, 11 and 13 were connected to grid node Nos. 2, 2 and 14, respectively. As shown in the table, with the increase in the charging power, the charger cost and total cost decreased. Figure 13 shows the decline curve of the cost. The total cost decreased with the cost of Term 2. The infrastructure layout, including charging stations selected and the connection scheme to the power grid, was identical for the ten instances, which indicated that the layout was very robust to the change of charging power.
It should be noted that the higher the charging power is, the heavier the thermal loads on batteries are, which results in an immediate effect in reducing the battery lifetime. Additional costs of battery replacement or maintenance may apply. It should be noted that the higher the charging power is, the heavier the thermal loads on batteries are, which results in an immediate effect in reducing the battery lifetime. Additional costs of battery replacement or maintenance may apply.

Driving Range
The current driving range of e-buses in Shenzhen is 250 km. In this subsection, the range was increased to 500 km, at an interval of 25 km. Two reasons explain why very long ranges were not set.
(1) The e-bus driving range is mainly decided by the battery capacity. To increase the range, heavier batteries will be loaded onto the e-buses and occupy more space on the carriages. E-buses must carry the batteries for every trip, which will consume large amounts of electricity and is, therefore, not desirable. An appropriate battery capacity would be preferable, rather than a larger battery. (2) Considering and comparing with the operation of traditional buses, the daily operating range of traditional buses is about 200 km on average. Therefore, although it is technically feasible to increase the driving range greatly, it is not economical and is unnecessary. It was assumed that the driving range is linearly dependent on the battery capacity, = × . Table 5 summarizes the cost and layout change with the increase in the driving range. The parameter settings of the first instance were the same as in Section 4.1. Figure 14 shows the decline curve of the costs. The total cost declined more rapidly than that shown in Figure 13. Although the technology was different, a rough comparison can be made by doubling the values of both parameters. When doubling the driving range, the total cost decreased to 5.24 million yuan (corresponding to a 38% decrease). The total cost was 7.33 million yuan when doubling the charging speed (corresponding to a 13% decrease). Even when the charging speed was multiplied by 10, the total cost was 6.48 million yuan, which was still higher than 5.24 million yuan. The improvement in the driving range could reduce the cost more than the advances in the charging technology.

Driving Range
The current driving range of e-buses in Shenzhen is 250 km. In this subsection, the range was increased to 500 km, at an interval of 25 km. Two reasons explain why very long ranges were not set. (1) The e-bus driving range is mainly decided by the battery capacity. To increase the range, heavier batteries will be loaded onto the e-buses and occupy more space on the carriages. E-buses must carry the batteries for every trip, which will consume large amounts of electricity and is, therefore, not desirable. An appropriate battery capacity would be preferable, rather than a larger battery.
(2) Considering and comparing with the operation of traditional buses, the daily operating range of traditional buses is about 200 km on average. Therefore, although it is technically feasible to increase the driving range greatly, it is not economical and is unnecessary. It was assumed that the driving range is linearly dependent on the battery capacity, u = 250 324 × C. Table 5 summarizes the cost and layout change with the increase in the driving range. The parameter settings of the first instance were the same as in Section 4.1. Figure 14 shows the decline curve of the costs. The total cost declined more rapidly than that shown in Figure 13. Although the technology was different, a rough comparison can be made by doubling the values of both parameters. When doubling the driving range, the total cost decreased to 5.24 million yuan (corresponding to a 38% decrease). The total cost was 7.33 million yuan when doubling the charging speed (corresponding to a 13% decrease). Even when the charging speed was multiplied by 10, the total cost was 6.48 million yuan, which was still higher than 5.24 million yuan. The improvement in the driving range could reduce the cost more than the advances in the charging technology.
The optimal layout solution may change with the parameters in two cases. The first one would be to open new facilities based on the original layout or close some existing ones. The second one would be to reselect a new group of facilities, which would be quite different from the original ones and, therefore, generate a huge cost compared with the first case. Table 5 shows that the layout was the same across the range from 250 km to 325 km. Since the range of 350 km, stations Nos. 9 and 11 were selected and both were connected to grid node No. 2. The layout difference was station No. 13. Although the layout solution changed when the driving range increased, the change was relatively small in the first case. Additional numerical experiments were conducted to find over what range the layout would continue to change. The results show that in the range 350-1075 km, the layout remained the same. Since the range of 1100 km, only station No. 11 was selected, and it was connected to grid node No. 2.   The optimal layout solution may change with the parameters in two cases. The first one would be to open new facilities based on the original layout or close some existing ones. The second one would be to reselect a new group of facilities, which would be quite different from the original ones and, therefore, generate a huge cost compared with the first case. Table 5 shows that the layout was the same across the range from 250 km to 325 km. Since the range of 350 km, stations Nos. 9 and 11 were selected and both were connected to grid node No. 2. The layout difference was station No. 13. Although the layout solution changed when the driving range increased, the change was relatively small in the first case. Additional numerical experiments were conducted to find over what range the layout would continue to change. The results show that in the range 350-1075 km, the layout remained the same. Since the range of 1100 km, only station No. 11 was selected, and it was connected to grid node No. 2.
Note that in Figure 14, for the range of 350 km, the cost of Term 3 increased suddenly, which was because the number of charging stations decreased from 3 to 2 and e-buses only had two charging options. Therefore, some e-buses would have longer charging trip distances than they had previously. However, the reduction in the facility costs, such as in Term 1 and Term 4, was larger than the increased cost of Term 3. Hence, the total cost still decreased.

Analyses and Findings
As indicated in Section 2.1, the value of the charging power only affected one input parameter of the model, the demand for chargers , which was reflected in the second term of the objective function. The cost of Term 2 only accounted for about 10% of total cost on average. The influence Note that in Figure 14, for the range of 350 km, the cost of Term 3 increased suddenly, which was because the number of charging stations decreased from 3 to 2 and e-buses only had two charging options. Therefore, some e-buses would have longer charging trip distances than they had previously. However, the reduction in the facility costs, such as in Term 1 and Term 4, was larger than the increased cost of Term 3. Hence, the total cost still decreased.

Analyses and Findings
As indicated in Section 2.1, the value of the charging power only affected one input parameter of the model, the demand for chargers d i , which was reflected in the second term of the objective function. The cost of Term 2 only accounted for about 10% of total cost on average. The influence exerted by the charging power was relatively small; hence, the layout did not change and the total cost did not decrease too much. The change in the driving range influenced (1) the demand for chargers by the change of X and increased battery capacity C; (2) the number of charging trips, further to the transportation cost of Term 3. Therefore, the optimization results changed more significantly with the increase in the driving range.
In summary, from the cost perspective, the advancements in both charging speed and driving range can reduce the cost, as expected. However, increasing the driving range can reduce the cost more than improving the charging speed. To maximize cost savings of the e-bus charging infrastructure, the R&D priority should be on increasing the driving range of e-buses. From the layout perspective, the model is remarkably robust for a wide range of charging speed values.

Real-World Transportation Network in Shenzhen
In this subsection, a case study was implemented to demonstrate the application of the proposed method in real-world settings and to examine whether the "No R" algorithm could solve large-scale problems.
The real-world transportation network of Shenzhen, China, was adopted. As mentioned in the Introduction, Shenzhen is the first megacity in the world to have electrified its bus transportation system. Since bus transportation is one of the main transport modes in China and is more popular than in most European and American countries, the scale of the bus transportation network is large. The operation data of a bus company, Shenzhen Bus Group, was adopted. It is the largest among the bus companies in Shenzhen, which accounts for about 1/3 of bus business. At the time of analysis, it operated 333 bus lines (data updated on Nov. 2018) and owned 5698 e-buses. The real operational data of the company contained the initial stops and final stops of the 333 bus lines. It was assumed that this bus company had 30 candidate e-bus charging stations. Figure 15 shows the distribution of the initial stops, final stops and candidate charging stations in Shenzhen. Python was used to inquire as to the longitude and latitude of the initial stops, final stops, and candidate charging stations. The API (application program interface) of Baidu Map was called to inquire the real driving distances w i 2 j and w i 3 j from initial stops and final stops to charging stations, respectively. The distance matrixes are attached in Supplementary Material 2. While the distance data from depot to charging stations w i 1 j was virtual, it was generated through a given rule, w i 1 j − w i 2 j ≤ 3 km, considering that the initial stops and depots are usually close geographically. The power network is the same as that in Section 4.1. The other parameter settings of the case study are attached in Part 3 of Supplementary Material 1.
In summary, from the cost perspective, the advancements in both charging speed and driving range can reduce the cost, as expected. However, increasing the driving range can reduce the cost more than improving the charging speed. To maximize cost savings of the e-bus charging infrastructure, the R&D priority should be on increasing the driving range of e-buses. From the layout perspective, the model is remarkably robust for a wide range of charging speed values.

Real-World Transportation Network in Shenzhen
In this subsection, a case study was implemented to demonstrate the application of the proposed method in real-world settings and to examine whether the "No R" algorithm could solve large-scale problems.
The real-world transportation network of Shenzhen, China, was adopted. As mentioned in the Introduction, Shenzhen is the first megacity in the world to have electrified its bus transportation system. Since bus transportation is one of the main transport modes in China and is more popular than in most European and American countries, the scale of the bus transportation network is large. The operation data of a bus company, Shenzhen Bus Group, was adopted. It is the largest among the bus companies in Shenzhen, which accounts for about 1/3 of bus business. At the time of analysis, it operated 333 bus lines (data updated on Nov. 2018) and owned 5698 e-buses. The real operational data of the company contained the initial stops and final stops of the 333 bus lines. It was assumed that this bus company had 30 candidate e-bus charging stations. Figure 15 shows the distribution of the initial stops, final stops and candidate charging stations in Shenzhen. Python was used to inquire as to the longitude and latitude of the initial stops, final stops, and candidate charging stations. The API (application program interface) of Baidu Map was called to inquire the real driving distances and from initial stops and final stops to charging stations, respectively. The distance matrixes are attached in Supplementary Material 2. While the distance data from depot to charging stations was virtual, it was generated through a given rule, | − | ≤ 3 km, considering that the initial stops and depots are usually close geographically. The power network is the same as that in Section 4.1. The other parameter settings of the case study are attached in Part 3 of Supplementary Material 1.  The planning results are shown in Table 6. Eleven charging stations were built to serve the 5698 e-buses of the company. The assignment scheme of buses to charging stations is attached in Part 3.2 of Supplementary Material 1.