Next Article in Journal
Experimental Evaluation of Conjugate Flaws on Rock Dynamic Fracturing
Next Article in Special Issue
Is a Carbon-Neutral Pathway in Road Transport Possible? A Case Study from Slovakia
Previous Article in Journal
Energy and Exergy Enhancement Study on PV Systems with Phase Change Material
Previous Article in Special Issue
Identifying a Country’s Freight Transport-Intensive Economic Sectors and Their Logistics Emissions—Method Development and Exemplary Evaluation with Austria
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Robustness of Battery Electric Bus Transit Networks under Charging Infrastructure Disruptions

1
Department of Civil Engineering, McMaster University, Hamilton, ON L8S 4L8, Canada
2
Department of Mathematics & Engineering Physics, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(4), 3642; https://doi.org/10.3390/su15043642
Submission received: 19 January 2023 / Revised: 6 February 2023 / Accepted: 10 February 2023 / Published: 16 February 2023

Abstract

:
The optimization of battery electric buses (BEBs) systems in transit is receiving considerable scholarly and practical attention. The practice is to minimize the total system cost to inform the optimal resource allocation. However, a minimization approach is insensitive to assessing and accommodating the robustness of BEB transit systems under disruption. This study evaluates the robustness of the BEB transit system under charging infrastructure disruption using complex network theory. The results of a mid-size multi-hub network indicate that the BEB system is robust against disruption if the disruption is resolved in a timely manner (within one hour). Furthermore, multi-charger charging stations have severe impacts on the system’s robustness. Overall, the BEB system robustness is more sensitive to the hourly number of buses charging at each station and the duration of the charging events.

1. Introduction

The relevance of incorporating electric drive technologies into public transit networks is expanding to mitigate transportation-related greenhouse gas (GHG) emissions and improve energy efficiency. This motivated several cities to adopt alternative-fuel buses, particularly electric buses (e-Buses), in the transit sector. E-Buses have long been seen as a viable option to substitute conventional internal combustion engine (ICE) buses due to their enhanced operation performance, noise reduction, and reduced operational costs [1,2]. Furthermore, e-Buses promise significant GHG reductions per passenger kilometer traveled, even if the electricity is not 100% carbon-free [3].
In general, e-Buses are categorized into three main types: Fuel Cell e-Bus (FCEB), which generates onboard electricity from fuel cells; Trolly e-Bus, which utilizes overhead wires that supply the electric motor with the energy continuously; and Battery e-Bus (BEB), which stores the energy on onboard batteries to power the electric motor [4]. The latter is considered more economically feasible as it enables the optimal utilization of the charging process during dwelling times. Moreover, compared to BEB, FCEB faces various practical challenges due to the lack of global technical regulations for hydrogen vehicles. Also, Trolley e-Buses require an extensive overhead-wiring infrastructure system [4].
That said, the utilization of BEBs in transit necessitates robust infrastructure planning and optimization, where each bus route must have a compatible charging scheme to ensure BEB’s adherence to the operating schedule. However, before phasing out fuel buses, it is vital for decision-makers and transit agencies to comprehend the implications of an all-electric bus fleet on the energy infrastructure [5,6,7,8].
Towards that end, previous BEB research has been carried out across three broad domains: (1) BEB components, including the motor, battery, and auxiliary systems that determine the energy consumption rates [9,10,11]. Typically, this research domain is tied to advancing BEB design, and its practical implications are geared toward bus manufacturers. (2) Charging infrastructure, including energy storage systems (ESS), charger’s capacity, the number of chargers, and their spatial distribution. This domain targets the optimal spatiotemporal allocation of BEB infrastructure to satisfy the energy demand/utilization [12,13]. These studies implement two approaches to charging BEBs: en-route and overnight. (3) Bus transit network modifications, which involve altering bus routes, location of bus stops, and transfer stations in favor of electrification [14,15].
Together, these domains are addressed using advanced optimization models, which aim at minimizing the total costs of adopting BEBs systems while mitigating their external impacts (e.g., on the utility grid and GHG emissions). Furthermore, previous work also considered the optimal operational feasibility, namely, operational features such as energy consumption, charging time, availability of chargers, and state of charge (SoC) [1,16,17]; economic features such as the costs of charging infrastructure, operation, maintenance, vehicle, and battery [18,19,20]; and environmental features including GHG emissions and the air quality [21,22].
Despite the plethora of well-established optimization models in the BEB literature, almost all these models are based on a minimization cost function, meaning that the BEB system configuration is designed without any redundancy. However, recent optimization approaches in BEB transit network configuration address this issue by accommodating the uncertainty related to transit operation and energy demand/supply. These uncertainties are often associated with energy consumption [23,24,25], arrival delay to charging stations [26], and power supply variability [27]. Acknowledging the systemic uncertainty is indeed essential to offer practical solutions for the adoption of BEBs in transit. However, the scope of these studies focuses on designing the optimal BEB transit system that accommodates these uncertainties yet not measuring the impacts of disruptions on the system’s performance.
As such, there is a clear gap in understanding and assessing the behavior of the BEB transit system during disruption (e.g., electricity outage or equipment malfunction), which has been recently documented [28]. Furthermore, the cascading impacts of disruption on the BEB transit network operation are still unknown. One can argue for the dire need to quantify the vulnerability of the BEB transit system to accurately design an optimal BEB system configuration.
Put another way, the works of [29,30,31] documented the vulnerability of transit networks (using ICE buses) to disruption. In this context, vulnerability indicates the ability of the BEB transit network to operate under the cascading failure resulting from a disruption in the network. In this context of BEB, additional disruptions could take place due to an electricity outage at any point in the utility grid or a malfunction in any of the system components (chargers).
Toward that end, this study aims to address two research questions:
(Q1) What is the robustness of BEB transit networks under disruptive events such as electricity outages or equipment malfunction?
(Q2) Which charging profile parameters have significant impacts on the robustness of BEB transit systems?
In the context of this study, we follow [32] in their definition of robustness is defined as the system’s ability to maintain the desired level of service despite disruption. In this study, transit system robustness refers to the ability of the transit system to deliver scheduled trips. As such, robustness is measured using a service frequency indicator (detailed in Section 3).
Overall, addressing these three research questions brings three substantial contributions to the BEB research:
-
First, we quantify the robustness of a BEB transit network using complex network theory.
-
Second, we analyze the sensitivity of the BEB service robustness to several operational parameters.
-
Lastly, this is the first attempt to evaluate the robustness of BEB transit system networks using complex network theory, offering some new insights into the design, planning, and optimization of BEB transit networks.
Following this introduction, a literature review of the current BEB optimization studies is discussed in Section 2. Section 3 presents the methodology, including the optimization model and the proposed complex network theory model. Section 4 describes the case study. Section 5 discusses the results, optimization, and robustness. Section 6 explains the sensitivity analysis results. While section seven discusses the findings and presents the concluding remarks.

2. Literature Review

The state-of-the-art for optimizing BEB transit system configurations, along with the associated infrastructure and the charging profile, focuses on optimizing the system costs. That is to minimize the total cost of ownership (TCO) associated with the implementation of BEB fleets in transit networks through the optimal allocation of resources. Previous studies related to BEB system configuration optimization are categorized into three groups; cost, utility impact, and GHG remissions, as listed in Table 1.
Cost optimization models, where the main aim is to minimize the TCO of the BEB system, considering the number of chargers, their locations, and the charging schemes that fulfill the transit schedules limits [1,8,17,33,34,38]. Furthermore, various studies developed optimization models for fast-charging infrastructure system that minimizes the system cost by reaching the optimal battery size for each bus and their charger-rated power [2,7,12,19,36].
In comparison, cost-utility optimization models aim to minimize the operational cost and the impact on the utility grid [37]. Moreover, spatial optimization models are proposed to quantify the number of chargers required for a BEB fleet to minimize construction, operational and maintenance costs [35]. At the same time, cost-emission optimization models aim at minimizing the BEB system cost and the life-cycle GHG emissions for both the operation stage and the end-of-life stage [22].
Another distinct stream of research accommodates the uncertainty of BEB operation through the utilization of Robust and Two-stage Stochastic optimization approaches), such as the works of [13,24,25,26,39,40,41,42]. This approach considers the uncertainty associated with several parameters, as detailed in Table 2. The advantages of the two-stage stochastic and the robust optimization models are their ability (1) to design the system for the worst-case scenario and (2) to accommodate the uncertainty distribution of several parameters at the same time. However, their limitations, with respect to the scope of this study, are the inability to assess the robustness of the BEB system under disruptive events. The latter is the main aim of the present paper.
Overall, the literature adopts optimization models for BEB relying on a minimization function of the system cost, utility impact, and component sizing. However, a minimization function is not ideal for BEB system configuration, especially when challenged with disruption events (e.g., power outage) [43]. Indeed, there is a significant degree of uncertainty associated with the vulnerability of the BEB transit systems and their robustness against disruptive events. As such, it is essential to assess the robustness of BEB systems to inform the optimal BEB system configuration.

3. Methodology

3.1. Problem Description

Infrastructure disruption for BEB transit systems can happen due to two main reasons: electricity outages and charging station malfunction. Although there are other types of transit operation disruptions (e.g., congestions, etc.), here we focus on BEB charging infrastructure disruptions. Figure 1 represents a hypothetical example of the charging scheme for two BEBs operating on the same route (inbound and outbound). The BEBs are served by an en-route charging station. A disruption before the scheduled charging event (e.g., an electricity outage in the charging station) will have varying impacts on the operation of this route.
If the charging station is disrupted (e.g., no power or charger malfunction), BUS 1 will not be able to complete the assigned trip, given that the battery State of Charge (SoC) is lower than the required energy to complete the trip. In comparison, BUS 2 will complete the assigned trip as their SoC is higher, yet the scheduled charging event will be shifted to another time (if available).
Therefore, BEB’s infrastructure disruption entails the assessment of the spatiotemporal utilization of the charging network, BEB’s SoC, and transit timetables. Given the interdependency of these three elements, several complex scenarios arise during disruption.
First, some BEBs will not be able to complete their assigned trips on the available SoC (immediate operation impact). Second, some BEBs will be able to complete their trips; however, they will not be able to charge in the rolling horizon (shifted operation impact). Third, the shifted charging/energy demand will be higher than the capacity of the utility grid (cumulative utility impact).
Therefore, assessing BEB transit system disruption should consider the spatiotemporal features of (1) transit operation (timetable and BEB SoC over time), (2) the charging scheme including schedule, location, and utilization time, and (3) the utility grid specification.
To address this issue, we (1) optimized a BEB transit system using the dominant optimization technique in the literature. (2) Extract the resultant BEB system configuration, including battery size, chargers’ rated power, spatial allocation, and the charging schedule. (3) Model the BEB transit system as a directed weighted network and subject the network to charging station (node) distribution. (4) Quantify the BEB system robustness by evaluating the cascading impacts on the BEB transit system. (5) Conduct sensitivity analysis to identify how the operational parameters influence the robustness of BEB transit systems. Figure 2 details these procedures.
It should be noted that the contribution of the present study is not claimed from the optimization model. The contribution is associated with the measurement of BEB system robustness against charging infrastructure disruption. The optimization model is used as input to feed the network disruption model. Second, the study is not aimed at offering a service rescheduling solution, and we treat the scheduled timetable as a hard constraint in the paper that must be satisfied.

3.2. Optimization Model: Mathematical Formulation

In this work, the optimization model satisfies the objective functions and constraints based on recent BEB literature. Overall, the model minimizes the total system cost by considering infrastructure, fleet, and operation. The model satisfies three common assumptions:
  • The charging process is carried out using both en-route charging during the recovery time as well as overnight at the depot;
  • The charger-rated powers for all stations are homogenous [6];
  • The battery sizes for all buses are homogeneous, enabling flexibility of operation [44];
  • The model maintains the current fleet size and the operational timetable [13];
  • The model accommodates the electricity time of use (ToU) tariff [24].
In particular, the model identifies the optimal number and locations of the charging stations, which are selected from a set of candidate stations ( I ). Furthermore, the model quantifies the optimal configuration of each charging station i I , (i.e., the number of chargers N i and the rated power P c h ). From a fleet perspective, the model chooses the optimal battery capacity for the fleet E b a t from a set of battery capacities A b a t . A list of abbreviations and notations used in the BEB system optimization model is presented in Table 3.
The objective function of the model is described in Equation (1). F s y s t e m denotes the overall BEB system annual cost, which includes charging station construction cost F c o n s , chargers cost F c h a r g e r s , fleet cost F f l e e t , maintenance cost F m a i n t , and operation cost F o p .
All cost parameters in Equation (1), except the operational cost, are calculated for the system lifespan. Therefore, these parameters are multiplied by an annualized factor related to the lifespan β and the discount rate μ 0 to annualize the cost of the system.
F s y s t e m = ( F c o n s + F c h a r g e r s + F f l e e t + F m a i n t ) ( μ 0 ( 1 + μ 0 ) β ( 1 + μ 0 ) β 1 ) + F o p
The five individual costs in Equation (1) are calculated using Equations (2)–(6). The construction cost is estimated in (Equation (2)), Where x i ,   i I is a binary decision variable that denotes whether a charging station will be installed in location i or not, and I is the set of candidates charging stations. The charger cost is considered in (Equation (3)) as a linear function of the charger-rated power with a constant [2], where N i is the number of chargers deployed in location i , P c h is the charger’s power (kW), C c h is the charger cost related to the charger power ($/kW), and C f c h is the fixed cost ($).
The fleet cost Is presented in Equation (4) in two parts; the cost of the battery, C b a t ($/kWh), which is related to the battery capacity E b a t (kWh), and the cost of the bus, C f b u s ($). The maintenance cost is calculated in (Equation (5)) as a percentage, r m , of the purchase cost of the infrastructure, chargers, and fleet costs. Lastly, the operational cost is the electricity cost, which is related to the ToU tariff C e ( t ) ($/kWh), as described in (Equation (6)). Where θ is the number of network workdays, T s is the timeslot duration (h), P c h is the charger-rated power (kW), and y b , j , i , t is a binary decision variable that denotes whether bus b B will charge after sub-trip j J b in location i I during timeslot t T or not.
F c o n s = i I C c o n s x i
F c h a r g e r s = i I N i ( P c h C c h + C f c h )
F f l e e t = b B ( C b a t E b a t + C f b u s )
F m a i n t = r m ( F c o n + F c h a r g e r s + F f l e e t )
F o p = θ T s [ b B t T i I j J b C e ( t ) P c h y b , j , i , t ]
The total annual system cost in (Equation (1)) will be minimized under a set of constraints. First, the battery capacity constraints are presented in Equations (7)–(9). In Equations (7) and (8), the arrival ( S b , j , i a r r )\departure ( S b , j , i d e p ) battery energy for the bus b at\from the candidate charging stations i after\before accomplishing the sub-trip j ( S b , j , i a r r ) should be greater\less than a minimum\maximum threshold ratio ( r m i n \ r m a x ) from the battery size E b a t . These constraints are applied to all buses, b B , and for all sub-trips, j J b , except for the departure for the first sub-trip. The battery capacity equals the maximum threshold as presented in (Equation (9)).
S b , j , i a r r   r m i n E b a t b B ,   j J b ,   i I  
S b , j , i d e p r m a x E b a t b B ,   j J b ,   i I  
S b , 1 d e p = r m a x E b a t b B  
Equation (10) dictates that the arrival battery capacity for any bus b from a sub-trip j in location i is equal to the departure battery energy from the previous location i mines the energy consumed during the sup-trip j from i to i ( E b , j , i c ) . The energy consumption rate for each sub-trip is taken as a linear function of the battery capacity. The energy consumption for each sub-trip is calculated by multiplying the driving distance d b , j of the sub-trip j by the energy consumption rate ( e b , j f + e b , j b a t E b a t ) and the sub-trip route factor w b , j which is related to the road and traffic conditions following [5].
S b , j , i a r r = S b , j , i d e p E b , j , i c b B ,   j J b ,   i & i I
E b , j , i c = d b , j ( e b , j f + e b , j b a t E b a t ) w b , j b B ,   j J b ,   i I
The charging strategy constraints in Equations (12)–(16) state that the departure battery energy S b , j + 1 , i d e p of bus b from location i for sub-trip j + 1 is equal to the summation of the arrival energy S b , j , i a r r from the previous sub-trip j and the charged energy during the recovery time R b , j , i , if it exists. The charged energy during the recovery time is related to the charger-rated power P c h , charger efficiency η c h , and the decision if the bus was charged or not y b , j , i , t . The charger decision of bus b after sub-trip j in location i during the timeslot t ( y b , j , i , t ) is set to zero during the operation time t R b , j , i , as mentioned in (Equation (13)). As the time is discretized, Constraints (14)–(16) are used to ensure the charging continuity of each charging event. These constraints are drawn from the work of [45]. While ρ b , j , i , t and σ b , j , i , t are auxiliary variables used to calculate the change of the charging state during the recovery time. In addition, Constraint (16) emphasizes that if the bus charges, the change from the state of not charging to charging is only one and the same for the change from charging to not charging to ensure the charging continuity.
S b , j + 1 , i d e p = S b , j , i a r r + t R b , j , i η c h T s P c h y b , j , i , t b B ,   j J b ,   i I
y b , j , i , t = 0 b B ,   j J b ,   i I ,   t R b , j , i
ρ b , j , i , t y b , j , i , t y b , j , i , t + 1 b B ,   j J b ,   i I ,   t R b , j , i
σ b , j , i , t     y b , j , i , t y b , j , i , t 1 b B ,   j J b ,   i I ,   t R b , j , i
t R b , j , i ρ b , j , i , t = t R b , j , i σ b , j , i , t 1 b B ,   j J b ,   i I ,   t R b , j , i
For the charging station constraints, (Equation (17)) ensures no charger deployment in location i without selecting the location as a charging station. In addition, it constrains the number of chargers deployed in location i by an upper limit N i m which is related to the available area in site i . In Equation (18), the number of buses charged in location i should be lower than the number of chargers available.
  N i   N i m   x i i I
b B y b , j , i , t N i i I ,   t T
For the variable’s types of constraints, Equations (19)–(25) emphasize that x i , y b , j , i , t ,   ρ b , j , i , t , and σ b , j , i , t are binary, N i is a non-negative integer, and P c h and E b a t are selected from predefined sets containing different finite predefined levels.
x i { 0 , 1 } i I
N i Z 0 + i I
P c h A c h
E b a t A b a t
y b , j , i , t { 0 , 1 } b B ,   j J b ,   i I ,   t R b , j , i
ρ b , j , i , t { 0 , 1 } b B ,   j J b ,   i I ,   t R b , j , i
σ b , j , i , t { 0 , 1 } b B ,   j J b ,   i I ,   t R b , j , i
The utilized optimization model for the BEB system design is formulated as follows:
M i n   ( 1 )  
s . t .   ( 7 25 )
The optimization model in this formulation is represented as an integer non-linear programming. The non-linearity exists in Equations (3), (6) and (12). The charger-rated power ( P c h ) is the common variable in all the non-linear terms. In addition, the model considers the homogeneity feature of the charger-rated power for the entire network ( P c h ). Therefore, the utilized optimization model could be linearized using a scenario-based optimization approach by solving the model several times using each level of the charger-rated power in A c h set (finite set), compare them and select the one that minimizes the objective function. In this case, the linear model will be solved | A c h | times, and it will be converted to an integer linear programming (ILP) model.
Several methods and algorithms have been developed to solve ILP models [46]. Branch and bound, cutting plane algorithm, dynamic programming, linear programming relaxation methods (e.g., Lagrangian relaxation), metaheuristics, and population-based evolutionary algorithms [47]. These methods can be used alone or in combination to solve ILP problems such as branch and price, branch and cut, and decomposition methods (e.g., Benders’ decomposition and column generation technique) [48]. A detailed description of each method can be found in [47].
Any commercial solver developed for these methods could handle the ILP model in a reasonable computational time (e.g., Gurobi and CPLEX). As such, and based on the formulation of the proposed model, the model is solved using the Gurobi solver.

3.3. Complex Network Representation for BEB Transit System

The BEB transit system consists of e-buses operating several tips that connect between stops/stations on predefined routes, as shown in Figure 3a. The BEB transit system is represented through a Complex Network approach. As such, the system is modeled as a directed weighted network that consists of nodes and links, as shown in Figure 3b.
The charging stations (output of the optimization) represent nodes. In comparison, links are represented by the number of BEBs connected to each charging station. Using this classification, the Degree Centrality measure ( C D ) in Equation (26) [49] represents the total number of bus links connected to each charging station in the network (undirected network), which is divided into in-degree (Equation (27)) and out-degree (Equation (28)) (directed network). Where, a i j and a j i are the number of buses going to and from j , respectively. The total number of buses operating on a link between i and j is a i j = a i j + a j i . J is the set of stations connected to j .
Degree   Centrality   C D C D ( j ) = i = 1 J a i j
In - Degree   Centrality   ( C D i n ) C D i n ( j ) = j = 1 J a i j
Out - Degree   Centrality   ( C D o u t ) C D o u t ( j ) = j = 1 J a j i
The robustness of the BEB system is evaluated under charging infrastructure disruption using the service frequency indicator R t f r e q in Equation (29). It estimates, during the duration of the disruption, the losses in transit trips due to the disruption. In other words, BEBs that will not be able to charge during disruption and their SoC cannot satisfy the next trip will be deemed unavailable, and their assigned trips will be canceled.
Numerically, it represents the ratio between the frequency of bus trips (for each time t T ) operating on the network after disruption ( N t D ) to the total daily frequency ( N T o t a l ).
R t f r e q = N t D N T o t a l × 100     t T
It should be noted that: (1) BEBs are removed from operation if the available battery SoC is insufficient to complete their trips during the disruption time. In case SoC values are sufficient, no BEBs are removed from the network. (2) We did not re-distribute BEBs to other charging stations outside their routes, as the reallocation will also result in canceling the assigned trips.
Furthermore, two disruption scenarios are considered for the charging stations (nodes). First, we assume a node disruption that will be resolved within the next hour. This is referred to as hourly disruption. In this scenario, each charging station is disrupted for one hour, and the impacts of this disruption are assessed for each charging station. Second, we assume a node disruption that will be resolved by the end of the operation day. This is referred to as daily disruption. During the daily disruption, the charging station is assumed to be out of service from the disruption time till the next day of operation. In this case, the impact of the daily disruption during each operation hour is also assessed for each charging station independently. For each scenario, the two robustness indicators are quantified, and the impact of each charging station on the system’s robustness is assessed.
A methodological flowchart summarizing the methods used in this research is shown in Figure 4.

4. Case Study

The Guelph bus transit network is selected as the case study. Located in Guelph City, Ontario, Canada, this medium-size multi-hubs network operates a fleet of 55 buses that travel 5,144,238 km (478,120 trips) annually while transporting around 110,000 passengers daily.
The network dataset is collected for weekday operation from the REMIX platform [50] and timetable data [51]. Both are based on the general transit feed specification (GTFS) data of Guelph. A brief description of Guelph’s operation data is illustrated in Table 4.
Guelph network consists of 23 bus routes operating through 506 bus stops/stations (see Figure 5a). Using the longitude and latitude data for each bus stop/station, we estimated the distances between each pair of stops/stations; then, we estimated the consumed energy for each bus trip using the calibrated model [5,9]. Stations that serve more than one route (e.g., transfer and end stations) are identified as candidate locations for charging stations “candidate charging stations”. This follows previous research, such as the work of [12,16,37,52]. As a result, we obtained 19 candidate locations for charging stations, as shown in Figure 5b.

5. Results

5.1. Optimal BEB System Configuration

The model is coded in MATLAB and solved using the GUROBI solver. The model runs on a personal computer with Intel® Core i5, 16 GB Ram, and a 4.20 GHz CPU. Overall, there are 67,011 variables and 82,703 constraints in the study. The optimal solution was reached in 4 min and 26.17 s.
The results indicate that BEBs could be implemented for Guelph transit. The system configuration includes 55 BEBs, each with a battery of 100 kWh (homogonous BEB fleet), with a total energy demand of 28,853.33 kWh/day. In addition, out of the 19 candidate charging stations, seven charging stations are required for the system with varying rated power and number of chargers (Table 5). The total annual system cost is $4,840,277.81, distributed as detailed in Figure 6.
The optimization model allocated seven charging stations out of the 19 candidate stations (Table 5). Charging station #1 includes two chargers, and charging station #4 contains four chargers. Therefore, the total number of chargers in the network is 11 chargers. The spatial distribution of the charging stations is depicted in Figure 7. Given that a homogonous charger power is assumed, all chargers have a rated power of 400 kW.
In particular, Table 5 details the locations of the charging stations, the number of chargers at each station, and the routes served by each charging station. It is observable from the data that Charging station #12, located at the depot, is shared between all the buses; hence it serves all routes at night. In comparison, Charging stations #13 and #15 serve fewer bus routes.
Collectively, the BEB fleet utilizes the charging infrastructure for 2656 min (Table 6—the sum of column 4), and the maximum charging duration for a bus is 37 min/hour. While the total number of charging events is 1104 events/day (Table 6—the sum of column 7), and the maximum number of hourly charging events per charger is seven events/hour (Table 6—column 9).
The seven charging stations satisfy the charging demands of 23 bus routes (55 buses). We listed the hourly charging flow for each operation hour in Table 6. The results also indicate a significant variation in the utilization rates of the charging stations. This is attributed to the varying spatiotemporal energy demand induced by the transit operation timetable and the varied BEB energy consumption. The two energy demand peaks could be observed (shaded cells in Table 6) from 11:00 AM to 4:00 PM and from 7:00 PM to 9:00 PM.
Furthermore, the charging utilization rates (number of charging events) per route are depicted in Table 7. As seen in the last column, the number of charging events varies significantly per route. This is attributed to the route features, including length, headway, and energy consumption.
Along the same lines, the utilization rate (time of use) per charging station varies significantly, as depicted in Table 8. The number of buses operating at each hour and the BEBs’ energy demand overtime for each charging station is shown in Appendix A.

5.2. BEB System Robustness

The results of the service frequency robustness indicator ( R t f r e q ) due to the hourly and daily disruption scenarios reported in Figure 8 and Figure 9, and the numerical results are available in Appendix B. Each line in the graphs should be interpreted as the impact of the disruption to a charging station on the entire network. Please note that the operation ends at 1:00 AM; however, Charging Station # 12 (Depot) operates overnight. Hence the x-axis is fixed to 24 h.
The hourly disruption results (Figure 8) show a relatively small fluctuation in the service frequency robustness indicator associated with the disruption of each charging station. The trend is almost similar over time. In other words, an hourly disruption to any charging station would have the same impact on the service. The service is likely to operate at 97–100% of its original frequency, indicating a robust operation.
In comparison, for multi-charger stations (i.e., charging stations #1 and #4), the results show that the impact of an hourly electricity outage (disruption to all chargers) is significant compared to equipment malfunction (disruption to single chargers). However, the magnitude of impact is still marginal, with robustness values between 97.9% and 97.1% for chargers #1 and #4, respectively.
A daily disruption (Figure 9) to any charging station will reduce the frequency of the service significantly (up to 57%). Like the hourly disruption, multi-charger stations severely impact the service in the case of an electricity outage. It should be noted that each dot in Figure 9 represents the daily impact on the service frequency due to the disruption to the charger at the given time. It also assumes that the disruption will last until the end of the operation day. In other words, each dot is one disruption scenario.

5.3. Sensitivity of BEB Robustness: A Discussion

A Sensitivity Analysis is performed to ascertain how the hourly frequency robustness of the BEB transit system is affected by the BEB charging schedule parameters. These parameters include (a) the number of BEBs charging at each charging station during each hour; (b) the number of charging events during each hour per station; (c) the charging duration, which represents the time that each BEB requires to charge its batteries at each hour; and (d) energy demand represents the total energy demanded at each charging station during each hour.
This study uses a global variance-based sensitivity analysis method, Sobol Indices [53]. Sobol Indices consider the entire input domain and provide a means of gauging the interactions between the parameters. It calculates the robustness variance and decomposes it into input parameter contributions. Two Sobol indices are studied: the first-order effect index S i , which represents the effect of X i alone, and the total-order effect index S T i , considers the interactions of X i with other parameters as follows:
S i = V [ E ( Y X i ) ] V ( Y )
S T i = 1 V [ E ( Y / X i ) ] V ( Y )
where V [ . ] and E ( . ) are the variance and expected values, respectively. The Sobol sampling [54] and Saltelli estimator [55] are utilized to estimate these two indices (Figure 10).
The sensitivity results (Figure 10) show that the number of BEBs charging at each station has the greatest impact on the robustness of the BEB transit system for both first-order and total effects. The number of charging events and their duration have the second and third highest impacts on the BEB system hourly frequency robustness, respectively. In contrast, the energy demand has the lowest influence on robustness.
These results indicate that: First, the charger’s power is less impactful, from a service robustness perspective, compared to the number of charging events. In other words, optimizing the charging process to reduce the number of charging events is critical. Second, allocating multiple chargers at the same location severely impacts service robustness. We recommend spatially spreading the charging infrastructure to distribute the risk of having multiple chargers at the same location. This is critical to the depot charging concept.

6. Conclusions

Current BEB infrastructure optimization models adopted in the literature rely on a minimization function of the total system cost, utility impact, and component sizing. However, a minimization function is not optimal for BEB system configuration, particularly when challenged with disruptive events. Furthermore, the cascading impacts of disruptions on the BEB network are still unidentified [28]. As such, there is a clear gap in understanding and assessing the behavior of the BEB transit system during disruption.
This paper answers two primary research questions: First, what is the robustness of the BEB transit system during disruptive events such as electricity outages or equipment malfunction? Second, which charging process parameters significantly affect the robustness of the BEB transit system?
To address these questions, we developed a BEB system configurations optimization model applied to a Guelph bus transit network as a case study. The resultant BEB system configuration is extracted, including the battery size, chargers’ spatial distribution, chargers’ power, and charging schedule. The extracted results are used as input in a complex network, modeled as a directed weighted graph. The resultant charging stations are subjected to disruptive events representing hourly and daily (discrete events) electricity outages and/or equipment malfunction. The cascading impacts on the BEB transit system are evaluated to quantify the robustness of the bus transit network through a service frequency indicator. Furthermore, the sensitivity of the BEB transit system robustness towards the charging schedule parameters is analyzed.
The results show that the hourly disruption slightly impacts the BEB service frequency robustness. Still, the service can operate with more than 97% of its original frequency, which indicates a robust operation. In contrast, the daily disruption will diminish the service frequency by up to 57% of the total service frequency.
The sensitivity results show that the charging schedule has a noticeable impact on the robustness of the BEB system. The number of charging events and their duration significantly affect the BEB system’s robustness. Therefore, reducing the number of charging events is critical to attaining a robust operation. Besides, the multiple chargers assigned to a single location have a deleterious effect on the robustness of the BEB service.
Overall, we recommend transit operators reduce the number of charging events per charger. This could be addressed by selecting a larger battery size for BEBs or increasing the number/location of charging stations. Risk disruption is equally critical, which could be achieved by spreading the charging infrastructure to distribute the risk of having multiple chargers at the same location.
For scholars, it is worth noting that we presented a reactive robustness assessment of the BEB transit system. Although valid, it would be more beneficial to proactively address the BEB robustness in the optimization process by integrating robustness thresholds in the problem formulation. This item should be addressed in future research studies.

Author Contributions

Conceptualization, M.M. and H.A.; methodology, H.A. and A.F.; software, H.A. and A.F.; validation, M.M.; formal analysis, H.A. and A.F.; investigation, H.A. and A.F.; resources, M.M.; data curation, H.A.; writing—original draft preparation, H.A. and A.F..; writing—review and editing, M.M.; visualization, M.M.; supervision, M.M.; project administration, M.M.; funding acquisition, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to acknowledge support from the Natural Sciences and Engineering Research Council of Canada (NSERC) Grant No: RGPIN-2018-05994.

Data Availability Statement

The data will be made available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. BEB System Configuration

Table A1. Number of buses operating each hour.
Table A1. Number of buses operating each hour.
Route Name1A1B2A3A3B4567891011121415162050565758GCTotal
Fleet Size32411314221121111331222255
5:00 AM001100100011100010000007
6:00 AM3233314221121111330011140
7:00 AM2235313121101110320000032
8:00 AM2234313121111110330002237
9:00 AM2135314221111111231000036
10:00 AM2234314121101111320002136
11:00 AM2234314221121111321222144
12:00 PM2234314221111111331222144
1:00 PM2234314221111111321111140
2:00 PM2234314221111111321212142
3:00 PM2235314221121111321122245
4:00 PM1236314221111111331222246
5:00 PM0235303020101000330000026
6:00 PM2133312121100110320000128
7:00 PM3234314221121111331222146
8:00 PM2236314221121111331212045
9:00 PM2234314121121111331222043
10:00 PM2232314221101111320121036
11:00 PM2233314221121111221100036
12:00 AM2232314221121111320000034
1:00 AM0120202001000001100010011
2:00 AM000010111000000002000006
3:00 AM002000011001000001000006
4:00 AM111001000000110000000006
Total3838637860197533401920242019181657501118192314772
Table A2. BEBs energy demand overtime for each charging station (kW at each hour).
Table A2. BEBs energy demand overtime for each charging station (kW at each hour).
Charging Station14712131519MinimumMaximumTotal
5:00 AM2932931471471471475871475873227
6:00 AM121312135275275275271680527168011,600
7:00 AM867867390390390390118739011878307
8:00 AM117311735205205205201600520160011,360
9:00 AM108010804874874874871453487145310,440
10:00 AM109310934804804804801440480144010,453
11:00 AM194719478438438438432653843265319,093
12:00 PM170717077577577577572427757242717,200
1:00 PM152015206676676676672107667210714,800
2:00 PM172017207707707707702427770242717,080
3:00 PM173317337807807807802507780250717,787
4:00 PM189318938208208208202667820266718,880
5:00 PM787787370370370370109337010937613
6:00 PM800800360360360360108036010807720
7:00 PM180018007877877877872480787248017,800
8:00 PM197319738808808808802867880286719,800
9:00 PM184018408238238238232653823265318,587
10:00 PM149314936436436436432067643206714,520
11:00 PM140014006336336336331960633196013,600
12:00 AM134713475835835835831840583184012,933
1:00 AM4404402172172172177332177334533
2:00 AM2402401201201201204801204802640
3:00 AM18718793939393373933732053
4:00 AM3073071531531531536131536133373
Minimum1871879393939337393
Maximum197319738808808808802867 2867
Total28,85328,85312,85012,85012,85012,85040,973 285,400

Appendix B. BEB Robustness

Table A3. Service frequency under charging station hourly disruption.
Table A3. Service frequency under charging station hourly disruption.
Station ID11-11-244-14-24-34-4712131519
Time of Disruption
5:00 AM100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%
6:00 AM98.64%99.46%99.18%97.68%99.05%99.32%99.59%99.73%99.25%99.86%99.73%99.73%100.00%
7:00 AM99.18%99.32%99.86%97.82%98.77%99.86%99.46%99.73%99.39%100.00%99.73%99.86%100.00%
8:00 AM99.05%99.46%99.59%97.89%98.84%99.46%99.73%99.86%99.39%100.00%99.86%99.86%99.32%
9:00 AM99.05%99.46%99.59%97.75%99.11%99.32%99.59%99.73%99.46%100.00%99.73%99.73%99.80%
10:00 AM98.77%99.05%99.73%97.28%99.18%99.18%99.32%99.59%99.52%100.00%99.86%99.86%99.52%
11:00 AM98.09%98.98%99.11%97.82%99.18%99.46%99.59%99.59%99.32%100.00%99.73%99.86%98.84%
12:00 PM98.09%98.77%99.32%98.02%99.39%99.46%99.46%99.73%99.52%100.00%99.73%99.86%98.64%
1:00 PM98.77%99.18%99.59%97.21%98.98%99.46%99.32%99.46%99.46%100.00%99.73%99.86%99.39%
2:00 PM98.57%98.98%99.59%98.02%99.39%99.32%99.46%99.86%99.39%100.00%99.73%99.86%99.05%
3:00 PM98.43%99.18%99.25%97.41%99.25%99.11%99.46%99.59%99.11%100.00%99.73%99.73%98.98%
4:00 PM97.96%98.84%99.11%97.14%99.25%99.18%99.39%99.32%99.25%100.00%99.73%99.86%98.77%
5:00 PM99.18%99.32%99.86%98.57%99.11%99.66%99.80%100.00%99.46%100.00%99.86%99.73%100.00%
6:00 PM99.32%99.32%100.00%98.23%99.32%99.59%99.46%99.86%99.59%100.00%99.86%99.86%99.86%
7:00 PM97.89%98.77%99.11%97.41%99.18%99.32%99.46%99.46%99.32%99.73%99.86%99.86%98.57%
8:00 PM98.16%98.91%99.25%97.28%98.91%99.46%99.18%99.73%99.66%99.86%99.73%99.73%98.91%
9:00 PM98.16%98.77%99.39%97.68%99.05%99.73%99.46%99.46%99.52%99.93%99.73%99.86%98.91%
10:00 PM98.30%98.84%99.46%97.41%99.05%99.32%99.46%99.59%99.39%100.00%99.86%99.86%99.59%
11:00 PM98.64%99.05%99.59%97.96%98.77%99.86%99.59%99.73%99.52%99.93%99.86%99.73%100.00%
12:00 AM99.18%99.52%99.66%98.84%99.52%99.80%99.66%99.86%99.86%99.66%99.93%99.93%100.00%
1:00 AM100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%100.00%
2:00 AMNANANANANANANANANANANANANA
3:00 AMNANANANANANANANANANANANANA
4:00 AMNANANANANANANANANANANANANA
5:00 AM–next dayNANANANANANANANANANANANANA

References

  1. Xiong, J.; Bai, C.; Guan, W.; Chen, Y.; Xu, Y. Mixed Optimization on Vehicle Scheduling and Recharge Scheduling of Plug-In Electric Buses with Consideration of Partial Recharge. J. Transp. Eng. Part A Syst. 2022, 148, 04021111. [Google Scholar] [CrossRef]
  2. He, Y.; Song, Z.; Liu, Z. Fast-charging station deployment for battery electric bus systems considering electricity demand charges. Sustain. Cities Soc. 2019, 48, 101530. [Google Scholar] [CrossRef] [Green Version]
  3. Soukhov, A.; Mohamed, M. Occupancy and GHG emissions: Thresholds for disruptive transportation modes and emerging technologies. Transp. Res. Part D Transp. Environ. 2022, 102, 103127. [Google Scholar] [CrossRef]
  4. Mahmoud, M.; Garnett, R.; Ferguson, M.; Kanaroglou, P. Electric buses: A review of alternative powertrains. Renew. Sust. Energ. Rev. 2016, 62, 673–684. [Google Scholar] [CrossRef]
  5. Foda, A.; Mohamed, M.; Bakr, M. Dynamic Surrogate Trip-Level Energy Model for Electric Bus Transit System Optimization. Transp. Res. Rec. J. Transp. Res. Board 2022, 2677, 513–528. [Google Scholar] [CrossRef]
  6. El-Taweel, N.A.; Farag, H.E.Z.; Mohamed, M. Integrated Utility-Transit Model for Optimal Configuration of Battery Electric Bus Systems. IEEE Syst. J. 2020, 14, 738–748. [Google Scholar] [CrossRef]
  7. Liu, Z.; Song, Z.; He, Y. Economic Analysis of On-Route Fast Charging for Battery Electric Buses: Case Study in Utah. Transp. Res. Rec. J. Transp. Res. Board 2019, 2673, 119–130. [Google Scholar] [CrossRef]
  8. Rogge, M.; van der Hurk, E.; Larsen, A.; Sauer, D.U. Electric bus fleet size and mix problem with optimization of charging infrastructure. Appl. Energy 2018, 211, 282–295. [Google Scholar] [CrossRef] [Green Version]
  9. Abdelaty, H.; Al-Obaidi, A.; Mohamed, M.; Farag, H.E.Z. Machine learning prediction models for battery-electric bus energy consumption in transit. Transp. Res. Part D-Transp. Environ. 2021, 96, 102868. [Google Scholar] [CrossRef]
  10. Pamuła, T.; Pamuła, W. Estimation of the Energy Consumption of Battery Electric Buses for Public Transport Networks Using Real-World Data and Deep Learning. Energies 2020, 13, 2340. [Google Scholar] [CrossRef]
  11. Vepsäläinen, J.; Kivekäs, K.; Otto, K.; Lajunen, A.; Tammi, K. Development and validation of energy demand uncertainty model for electric city buses. Transp. Res. Part D Transp. Environ. 2018, 63, 347–361. [Google Scholar] [CrossRef]
  12. Liu, X.; Qu, X.; Ma, X. Optimizing electric bus charging infrastructure considering power matching and seasonality. Transp. Res. Part D Transp. Environ. 2021, 100, 103057. [Google Scholar] [CrossRef]
  13. An, K. Battery electric bus infrastructure planning under demand uncertainty. Transp. Res. Part C Emerg. Technol. 2020, 111, 572–587. [Google Scholar] [CrossRef]
  14. Liu, K.; Yamamoto, T.; Morikawa, T. Impact of road gradient on energy consumption of electric vehicles. Transp. Res. Part D: Transp. Environ. 2017, 54, 74–81. [Google Scholar] [CrossRef]
  15. Perrotta, D.; Macedo, J.L.; Rossetti, R.J.F.; Sousa, J.F.d.; Kokkinogenis, Z.; Ribeiro, B.; Afonso, J.L. Route Planning for Electric Buses: A Case Study in Oporto. Procedia-Soc. Behav. Sci. 2014, 111, 1004–1014. [Google Scholar] [CrossRef] [Green Version]
  16. Wei, R.; Liu, X.; Ou, Y.; Kiavash Fayyaz, S. Optimizing the spatio-temporal deployment of battery electric bus system. J. Transp. Geogr. 2018, 68, 160–168. [Google Scholar] [CrossRef]
  17. Uslu, T.; Kaya, O. Location and capacity decisions for electric bus charging stations considering waiting times. Transp. Res. Part D Transp. Environ. 2021, 90, 102645. [Google Scholar] [CrossRef]
  18. Rupp, M.; Rieke, C.; Handschuh, N.; Kuperjans, I. Economic and ecological optimization of electric bus charging considering variable electricity prices and CO2eq intensities. Transp. Res. Part D Transp. Environ. 2020, 81, 102293. [Google Scholar] [CrossRef]
  19. Kunith, A.; Mendelevitch, R.; Goehlich, D. Electrification of a city bus network—An optimization model for cost-effective placing of charging infrastructure and battery sizing of fast-charging electric bus systems. Int. J. Sustain. Transp. 2017, 11, 707–720. [Google Scholar] [CrossRef] [Green Version]
  20. Gormez, M.A.; Topcu, A.; Sozer, Y. Cost Optimization of an Opportunity Charging Bus Network. In Proceedings of the 2018 IEEE Energy Conversion Congress and Exposition (ECCE), Portland, OR, USA, 23–27 September 2018; pp. 2782–2789. [Google Scholar]
  21. Zhou, Y.; Liu, X.C.; Wei, R.; Golub, A. Bi-Objective Optimization for Battery Electric Bus Deployment Considering Cost and Environmental Equity. IEEE Trans. Intell. Transp. Syst. 2021, 22, 2487–2497. [Google Scholar] [CrossRef]
  22. Bi, Z.; Keoleian, G.A.; Ersal, T. Wireless charger deployment for an electric bus network: A multi-objective life cycle optimization. Appl. Energy 2018, 225, 1090–1101. [Google Scholar] [CrossRef]
  23. Alwesabi, Y.; Avishan, F.; Yanıkoğlu, İ.; Liu, Z.; Wang, Y. Robust strategic planning of dynamic wireless charging infrastructure for electric buses. Appl. Energy 2022, 307, 118243. [Google Scholar] [CrossRef]
  24. Liu, K.; Gao, H.; Wang, Y.; Feng, T.; Li, C. Robust charging strategies for electric bus fleets under energy consumption uncertainty. Transp. Res. Part D Transp. Environ. 2022, 104, 103215. [Google Scholar] [CrossRef]
  25. Liu, Z.; Song, Z.; He, Y. Planning of Fast-Charging Stations for a Battery Electric Bus System under Energy Consumption Uncertainty. Transp. Res. Rec. J. Transp. Res. Board 2018, 2672, 96–107. [Google Scholar] [CrossRef] [Green Version]
  26. Zheng, F.; Wang, Z.; Liu, M. Overnight charging scheduling of battery electric buses with uncertain charging time. Oper. Res. 2022, 22, 4865–4903. [Google Scholar] [CrossRef]
  27. Iliopoulou, C.; Kepaptsoglou, K. Robust electric transit route network design problem (RE-TRNDP) with delay considerations: Model and application. Transp. Res. Part C Emerg. Technol. 2021, 129, 103255. [Google Scholar] [CrossRef]
  28. Eldeeb, G.; Mohamed, M. Transit electrification state of the art: A machine-learning based text mining approach. Transp. Res. Part D Transp. Environ. 2022, 111, 103446. [Google Scholar] [CrossRef]
  29. Xu, Z.; Zhang, C.; Chopra, S.S. Robustness assessment of public bus transit system with a response-integrated approach for a resilient public transport system in Hong Kong. Transp. B Transp. Dyn. 2022, 1–15. [Google Scholar] [CrossRef]
  30. Jia, T.; Liu, W.; Liu, X. A cross-city exploratory analysis of the robustness of bus transit networks using open-source data. Phys. A Stat. Mech. Its Appl. 2021, 580, 126133. [Google Scholar] [CrossRef]
  31. Abdelaty, H.; Mohamed, M.; Ezzeldin, M.; El-Dakhakhni, W. Quantifying and classifying the robustness of bus transit networks. Transp. A Transp. Sci. 2020, 16, 1176–1216. [Google Scholar] [CrossRef]
  32. Barabási, A.L.; PÃ3sfai, M.Ã. Network Science; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
  33. Li, M.; Tang, P.; Lin, X.; He, F. Multistage planning of electric transit charging facilities under build-operate-transfer model. Transp. Res. Part D Transp. Environ. 2022, 102, 103118. [Google Scholar] [CrossRef]
  34. Benoliel, P.; Jenn, A.; Tal, G. Examining energy uncertainty in battery bus deployments for transit agencies in California. Transp. Res. Part D Transp. Environ. 2021, 98, 102963. [Google Scholar] [CrossRef]
  35. Wu, X.; Feng, Q.; Bai, C.; Lai, C.S.; Jia, Y.; Lai, L.L. A novel fast-charging stations locational planning model for electric bus transit system. Energy 2021, 224, 106. [Google Scholar] [CrossRef]
  36. Lotfi, M.; Pereira, P.; Paterakis, N.; Gabbar, H.A.; Catalão, J.P.S. Optimizing Charging Infrastructures of Electric Bus Routes to Minimize Total Ownership Cost. In Proceedings of the 2020 IEEE International Conference on Environment and Electrical Engineering and 2020 IEEE Industrial and Commercial Power Systems Europe (EEEIC/I&CPS Europe), Madrid, Spain, 9–12 June 2020; pp. 1–6. [Google Scholar]
  37. Lin, Y.; Zhang, K.; Shen, Z.-J.M.; Ye, B.; Miao, L. Multistage large-scale charging station planning for electric buses considering transportation network and power grid. Transp. Res. Part C Emerg. Technol. 2019, 107, 423–443. [Google Scholar] [CrossRef]
  38. Wang, Y.; Huang, Y.; Xu, J.; Barclay, N. Optimal recharging scheduling for urban electric buses: A case study in Davis. Transp. Res. Part E Logist. Transp. Rev. 2017, 100, 115–132. [Google Scholar] [CrossRef]
  39. Bie, Y.; Ji, J.; Wang, X.; Qu, X. Optimization of electric bus scheduling considering stochastic volatilities in trip travel time and energy consumption. Comput.-Aided Civ. Infrastruct. Eng. 2021, 36, 1530–1548. [Google Scholar] [CrossRef]
  40. Zhou, Y.; Wang, H.; Wang, Y.; Li, R. Robust optimization for integrated planning of electric-bus charger deployment and charging scheduling. Transp. Res. Part D Transp. Environ. 2022, 110, 103410. [Google Scholar] [CrossRef]
  41. Jiang, M.; Zhang, Y.; Zhang, Y.; Li, W. Optimal Electric Bus Scheduling under Travel Time Uncertainty: A Robust Model and Solution Method. J. Adv. Transp. 2021, 2021, 1–19. [Google Scholar] [CrossRef]
  42. Hu, H.; Du, B.; Liu, W.; Perez, P. A joint optimisation model for charger locating and electric bus charging scheduling considering opportunity fast charging and uncertainties. Transp. Res. Part C Emerg. Technol. 2022, 141, 103732. [Google Scholar] [CrossRef]
  43. Tessler, M.E.; Traut, E.J. Hurricane resiliency methods for the New York City electric bus fleet. Transp. Res. Part D Transp. Environ. 2022, 105, 103255. [Google Scholar] [CrossRef]
  44. Wang, Y.; Liao, F.; Lu, C. Integrated optimization of charger deployment and fleet scheduling for battery electric buses. Transp. Res. Part D Transp. Environ. 2022, 109, 103382. [Google Scholar] [CrossRef]
  45. Liu, K.; Gao, H.; Liang, Z.; Zhao, M.; Li, C. Optimal charging strategy for large-scale electric buses considering resource constraints. Transp. Res. Part D Transp. Environ. 2021, 99, 103009. [Google Scholar] [CrossRef]
  46. Darby-Dowman, K.; Wilson, J.M. Developments in linear and integer programming. J. Oper. Res. Soc. 2017, 53, 1065–1071. [Google Scholar] [CrossRef] [Green Version]
  47. Genova, K.; Guliashki, V. Linear integer programming methodsand approaches—A survey. Cybern. Inf. Technol. 2011, 11, 3–25. [Google Scholar]
  48. Kuma, P.H.; Mageshvaran, R. Methods And Solvers Used For Solving Mixed Integer Linear Programming And Mixed Nonlinear Programming Problems: A Review. Int. J. Sci. Technol. Res. 2020, 9, 1872–1882. [Google Scholar]
  49. Albert, R.; Barabási, A.-L. Statistical mechanics of complex networks. Rev. Mod. Phys. 2002, 74, 47–97. [Google Scholar] [CrossRef] [Green Version]
  50. Remix. Base Network_Hamilton Bus Transit Network. Available online: https://platform.remix.com/?explore=true&latlng=43.13105,-79.68847,8.061&layer=simple (accessed on 5 January 2023).
  51. Guelph_Schedule. Available online: https://transitfeeds.com/p/guelph-transit/53 (accessed on 23 December 2022).
  52. Lin, Y.; Zhang, K.; Shen, Z.-J.M.; Miao, L. Charging Network Planning for Electric Bus Cities: A Case Study of Shenzhen, China. Sustainability 2019, 11, 4713. [Google Scholar] [CrossRef] [Green Version]
  53. Sobol, I.M. Sensitivity Estimates for Nonlinear Mathematical Models. Math. Model. Comput. Exp. 1993, 4, 407–414. [Google Scholar]
  54. Sobol, I.M. On the distribution of points in a cube and the approximate evaluation of integrals. USSR Comput. Math. Math. Phys. 1967, 7, 86–112. [Google Scholar] [CrossRef]
  55. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Sensitivity Analysis: From Theory to Practice. In Global Sensitivity Analysis. The Primer; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2007; Volume 1, pp. 237–275. [Google Scholar]
Figure 1. A hypothetical example of the uncertainty in the BEB system’s vulnerability.
Figure 1. A hypothetical example of the uncertainty in the BEB system’s vulnerability.
Sustainability 15 03642 g001
Figure 2. Robustness Assessment Process.
Figure 2. Robustness Assessment Process.
Sustainability 15 03642 g002
Figure 3. A hypothetical example of the BEB representation in complex network theory. (a) BEB transit system; (b) BEB in complex network theory. Black dots represent bus stops/stations in the network; red lines represent bus routes; white circles represent the charging stations (nodes); dashed blue lines represent the directed bus trips (links).
Figure 3. A hypothetical example of the BEB representation in complex network theory. (a) BEB transit system; (b) BEB in complex network theory. Black dots represent bus stops/stations in the network; red lines represent bus routes; white circles represent the charging stations (nodes); dashed blue lines represent the directed bus trips (links).
Sustainability 15 03642 g003
Figure 4. Methodological flowchart.
Figure 4. Methodological flowchart.
Sustainability 15 03642 g004
Figure 5. The distribution of bus stops/stations and candidate charging stations in the Guelph bus network. (a) Bus stops/stations in the network (n = 506); (b) Candidate charging stations (n = 19).
Figure 5. The distribution of bus stops/stations and candidate charging stations in the Guelph bus network. (a) Bus stops/stations in the network (n = 506); (b) Candidate charging stations (n = 19).
Sustainability 15 03642 g005aSustainability 15 03642 g005b
Figure 6. The total annual costs of implementing the BEB system in Guelph city.
Figure 6. The total annual costs of implementing the BEB system in Guelph city.
Sustainability 15 03642 g006
Figure 7. Charging stations’ locations with ID.
Figure 7. Charging stations’ locations with ID.
Sustainability 15 03642 g007
Figure 8. Service frequency reduction due to hourly disruption to each charging station.
Figure 8. Service frequency reduction due to hourly disruption to each charging station.
Sustainability 15 03642 g008
Figure 9. Service frequency percentage due to daily disruption to each charging station.
Figure 9. Service frequency percentage due to daily disruption to each charging station.
Sustainability 15 03642 g009
Figure 10. Sensitivity analysis results (first- and total-order effects).
Figure 10. Sensitivity analysis results (first- and total-order effects).
Sustainability 15 03642 g010
Table 1. Recent BEB system configuration optimization studies.
Table 1. Recent BEB system configuration optimization studies.
AuthorObjective Functions (Minimize)BEB System Configuration
CostUtility
Impact
GHG
Emission
Li [33] Number of chargers
Xiong [1] Number of chargers
Benoliel [34] Number of chargers & Fleet size
Liu [12] Number of chargers & Chargers’ power
Uslu and Kaya [17] Number of chargers
Wu [35] Number of chargers & Fleet size
El-Taweel [6] Battery capacity, number of chargers & power
Lotfi [36] Battery capacity & Chargers’ power
He [2] Battery capacity & Chargers’ power
Lin [37]Number of chargers
Lin [37] Number of chargers
Liu [7] Battery capacity & Chargers’ power
Bi [22] Battery capacity & Number of chargers
Rogge [8] Number of chargers & Fleet size
Kunith [19] Battery capacity
Wang [38] Number of chargers
Table 2. BEB optimization studies with uncertainty.
Table 2. BEB optimization studies with uncertainty.
AuthorRobust and Two-Stage Stochastic Optimization
Energy
Consumption
Travel TimeCharging TimePassenger Load (Mass)Charging
Demand
Liu [24]
Zheng [26]
Zhou [40]
Bie [39]
Hu [42]
Jiang [41]
An [13]
Liu [25]
Table 3. Abbreviations and notations summary.
Table 3. Abbreviations and notations summary.
AbbreviationDescriptionIndexDescription
GHGGreenhouse gas i ,   i Index of candidate charging station location
e-BusesElectric buses b Index of BEBs
ICEInternal combustion engine j Index of sub-trips (the distance between each two consecutive candidate locations)
FCEBFuel cell electric bustIndex of timeslots
BEBBattery electric busVariablesDescription
ESSEnergy storage system F c o n s Cost of charging stations constructions ($)
SoCState of charge F c h a r g e r s Cost of chargers ($)
TCOThe total cost of ownership F f l e e t Cost of bus fleet, including battery ($)
ToUTime-of-use F m a i n t Cost of maintenance ($)
SetsDescription F o p Cost of electricity ($)
I Set of candidates charging stations locations F s y s t e m Total annual system cost
B Set of BEBs S b , j , i a r r Arrival battery energy of bus b after sub-trip j at charging station i (kWh)
J b Set of sub-trips of bus b S b , j , i d e p Departure battery energy of bus b before sub-trip j from charging station i (kWh)
A c h Set of charger-rated powers E b , j , i c Energy consumption of bus b during sub-trip j after departure from charging station i
A b a t Set of battery-rated energies
Decision
variables
Description
x i Binary   decision   variable ,   x i = 1 if the location i is   selected   to   build   a   charging   station ,   x i = 0 otherwise
N i A   non - negative   integer   indicates   the   number   of   the   chargers   deployed   in   location   i I
P c h Charger-rated power (kW)
E b a t BEBs battery-rated energy for all buses (kWh)
y b , j , i , t Binary decision variable, y b , j , i , t = 1 if the bus b charged at timeslot t in charging station i after sub-trip j, y b , j , i , t = 0 otherwise
ρ b , j , i , t Binary   variable ,   ρ b , j , i , t = 1 if the charged state of bus b after sub-trip j in the charging station i for the current timeslot t and the later timeslot t + 1 changes, ρ b , j , i , t = 0 otherwise
σ b , j , i , t Binary   variable ,   σ b , j , i , t = 1 if the charged state of bus b after sub-trip j in the charging station i for the current timeslot t and the previous timeslot t − 1 changes, σ b , j , i , t = 0 otherwise
ParametersDescriptionParametersDescription
C c o n s Cost of construction of candidate charging station i ($) e b , j f The energy consumption rate of bus b during sub-trip j (kWh/km)
C c h Charger cost that is related to its power ($/kW) e b , j b a t The energy consumption rate of bus b during sub-trip j caused by one unit of BEB battery size for 1 km distance (1/km)
C f c h charger fixed cost ($/unit) r m Maintenance cost percentage from purchase cost (%)
C b a t Battery cost ($/kWh) C e ( t ) Electricity rate in timeslot t depending on ToU ($/kWh)
C f b u s Bus cost without battery ($/bus) θ Number of workdays
T s Timeslot duration (min) μ 0 Discount rate (%)
w b , j Route factor of bus b in sub-trip j β Lifespan (years)
N i m Maximum limit of the number of chargers in charging station i r m i n Minimum limit (%)
R b , j , i Recovery time set of bus b after sub-trip j at charging station i r m a x Maximum limit (%)
η c h Charger efficiency (%) d b , j Length of sub-trip j of bus b (km)
Table 4. Timetable data for each bus route in the Guelph bus transit network.
Table 4. Timetable data for each bus route in the Guelph bus transit network.
Route ID T r s t a r t T r e n d T r t r i p ( min )   * Av. Headway Time (min)l(r)
(km)
Ntrips,r (#)
1A-College Edinburgh5:4524:15143018.01676
2A-West Loop Clockwise5:4524:15263035.35476
3A-East Loop Clockwise6:0022:30431324.528148
1B- College Edinburgh5:4524:15153018.91076
3B-East Loop5:4524:15243035.85664
4-York5:4524:15263010.45338
5-South Gordon5:4524:15553030.12476
6-Harvard Ironwood5:4524:15263014.58976
7-Kortright Downe5:4524:15263019.06576
8-Stone Road Mall5:4524:1526309.50438
9-West End Community Centre5:4524:15263011.18938
10-Imperia5:4524:15283010.50038
11-Willow West5:4524:15283010.13038
12-General Hospital5:4524:15263010.28938
14-Grange5:4524:1525309.75538
15-College Ave W5:4524:15263013.95176
16-Route 165:4524:15543032.11776
20-Northwest Industrial5:4524:15513029.74876
50-Route 508:0021:4013205.17842
56-Route 567:4521:4573011.44058
57-Route 577:4522:2515208.30590
58-Route 587:4521:5015458.52382
Gordon Corridor7:4519:00243614.91336
* The data represent inbound and outbound trips. For each route (r): T r s t a r t is the beginning of the operation time; T r e n d is the end of the operation time; T r t r i p is the time duration for each bus trip; l(r) is the total distance, and Ntrips,r is the total number of daily trips for route r .
Table 5. Configuration of the charging stations and the corresponding bus routes.
Table 5. Configuration of the charging stations and the corresponding bus routes.
Charging
Station ID
Station NameNumber of Chargers per Station (#)Bus Routes
1UC South Loop2 (1-1, 1-2)1A, 1B, 5, 6, 7, 15, 57
4GCS East4 (4-1, 4-2, 4-3, 4-4)2A, 3B, 4, 5, 8, 9, 10, 11, 12, 14, 16, 20
7Gordon St. at Harvard Rd.13A, 6, 7
12Depot11A, 1B, 2A, 3A, 3B, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 15, 16, 20, 50, 56, 57, 58, GC
13Goodwin Dr. at Ray Cres.12A, 5
15Woodlawn at Wal-Mart13B
19University150, 56, 57, 58, GC
Table 6. The hourly charging utilization.
Table 6. The hourly charging utilization.
HourNumber of Utilized Charging Stations (#)Charging Stations IDsCharging Duration (Minute)Number of Charging Events (#)Number of Buses at Each Charging Station (#)
Sum.Min.Max.Sum.Min.Max.
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)
5:00 AM1124444447117
6:00 AM61, 4, 7, 12, 13, 15180834571640
7:00 AM51, 4, 7, 13, 15130426381532
8:00 AM61, 4, 7, 13, 15, 19176634511637
9:00 AM61, 4, 7, 13, 15, 19160238461636
10:00 AM61, 4, 7, 13, 15, 19164632511636
11:00 AM61, 4, 7, 13, 15, 192921242681644
12:00 PM61, 4, 7, 13, 15, 19256638611644
1:00 PM61, 4, 7, 13, 15, 192281240571740
2:00 PM61, 4, 7, 13, 15, 192581242591642
3:00 PM61, 4, 7, 13, 15, 192601246671645
4:00 PM61, 4, 7, 13, 15, 192841242751646
5:00 PM51, 4, 7, 13, 15118228321526
6:00 PM51, 4, 7, 13, 15120224371628
7:00 PM71, 4, 7, 12, 13, 15, 192701240751746
8:00 PM71, 4, 7, 12, 13, 15, 192961240712745
9:00 PM71, 4, 7, 12, 13, 15, 192761246651643
10:00 PM71, 4, 7, 12, 13, 15, 192201044581636
11:00 PM61, 4, 7, 12, 13, 152101242471636
12:00 AM61, 4, 7, 12, 13, 15202842511634
1:00 AM31, 4, 1266244131311
2:00 AM1123636366126
3:00 AM1122828286126
4:00 AM1124646466116
Table 7. Hourly chagrining events per route.
Table 7. Hourly chagrining events per route.
Route Name1A1B2A3A3B4567891011121415162050565758GCTotal
Fleet Size32411314221121111331222255
5:00 AM001100100011100010000007
6:00 AM3433525462222221330011157
7:00 AM2235415131101120420000038
8:00 AM2244525461111110430002251
9:00 AM2145616242111112231000046
10:00 AM2344526242202211420002151
11:00 AM2344516351232223423235168
12:00 PM2444426242212212431242161
1:00 PM2344627232112123421112157
2:00 PM2344625342212121421213259
3:00 PM2245616561222221422125267
4:00 PM1446626462112225432242475
5:00 PM0245503050101000330000032
6:00 PM2143622141200110420000137
7:00 PM3344627472232122431254275
8:00 PM2447626342222223433224071
9:00 PM2444615242222122431345065
10:00 PM2442626451202222420132058
11:00 PM2233625332122112321100047
12:00 AM2342526451221113520000051
1:00 AM0120203001000001200010013
2:00 AM000010111000000002000006
3:00 AM002000011001000001000006
4:00 AM111001000000110000000006
Total3856797910732108559030322633272934755017193139181104
Values in the color-coded cells represent the number of hourly charging events for each route, and the last column and row represent the total.
Table 8. Hourly chagrining duration (minutes) per charger.
Table 8. Hourly chagrining duration (minutes) per charger.
Charger ID1-11-24-14-24-34-4712131519Total
5:00 AM00000004400044
6:00 AM2820301610816348100180
7:00 AM2242412101026012100130
8:00 AM2416341814630014128176
9:00 AM248381881020020122160
10:00 AM328221618162201686164
11:00 AM423642281820400221232292
12:00 PM36223418182036028638256
1:00 PM401632141622400221214228
2:00 PM421236281816420281224258
3:00 PM361636221816360221246260
4:00 PM423430261818400281236284
5:00 PM12228161402806120118
6:00 PM24024161041801482120
7:00 PM4028402216222022221226270
8:00 PM4024381626183828261230296
9:00 PM3820461618163024221234276
10:00 PM4418241416182810261210220
11:00 PM281242161612203418120210
12:00 AM36182422141484212120202
1:00 AM201064004400066
2:00 AM00000003600036
3:00 AM00000002800028
4:00 AM00000004600046
Total6323146343603002665383923662103084320
Values in the color-coded cells represent the duration (min) of charging events at each charger, and the last column and row represent the total.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abdelaty, H.; Foda, A.; Mohamed, M. The Robustness of Battery Electric Bus Transit Networks under Charging Infrastructure Disruptions. Sustainability 2023, 15, 3642. https://doi.org/10.3390/su15043642

AMA Style

Abdelaty H, Foda A, Mohamed M. The Robustness of Battery Electric Bus Transit Networks under Charging Infrastructure Disruptions. Sustainability. 2023; 15(4):3642. https://doi.org/10.3390/su15043642

Chicago/Turabian Style

Abdelaty, Hatem, Ahmed Foda, and Moataz Mohamed. 2023. "The Robustness of Battery Electric Bus Transit Networks under Charging Infrastructure Disruptions" Sustainability 15, no. 4: 3642. https://doi.org/10.3390/su15043642

APA Style

Abdelaty, H., Foda, A., & Mohamed, M. (2023). The Robustness of Battery Electric Bus Transit Networks under Charging Infrastructure Disruptions. Sustainability, 15(4), 3642. https://doi.org/10.3390/su15043642

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop