# Flexible Emergency Vehicle Network Design considering Stochastic Demands and Inverse-Direction Lanes

Authors: Hua Wang; Ling Xiao; Zhang Chen
Journal: Journal of Advanced Transportation
 (2018)

---

## Abstract

We study transportation network design with stochastic demands and emergency vehicle (EV) lanes. Different from previous studies, this paper considers two groups of users, auto and EV travelers, whose road access rights are differentiated in the network, and addresses the value of incorporating inverse-direction lanes in network design. We formulate the problem as a bilevel optimization model, where the upper-level model aims to determine the optimal design of EV lanes and the lower-level model uses the user equilibrium principle to forecast the route choice of road users. A simulation-based genetic algorithm is proposed to solve the model. With numerical experiments, we demonstrate the value of deploying inverse-direction EV lanes and the computational efficiency of the proposed algorithm. We reach an intriguing finding that both regular and EV lane users can benefit from building EV lanes.


---

## Body

## 1. Introduction

Emergency incidents, such as traffic accidents, mere spasmodic diseases, and fire disasters, have often been observed in metropolises with a high population density [1, 2], which generates an increasing travel demand of emergency vehicles (EVs) (e.g., ambulances and fire engines). Currently, transportation emergency evacuation and rescue have attracted greater attention from transportation researchers and practitioners than before. For example, So and Daganzo [3] proposed an inner-first-out control strategy to manage evacuation routes, and Daganzo and So [4] extended such a nonanticipative and adaptive control strategy to manage traffic networks in real-time emergency evacuations. Meanwhile, with advances in intelligent transport techniques, fast development in building various emergency evacuation systems is witnessed, for example, the corridor-based and region-wide emergency evacuation systems for Washington D.C. [5] and Baltimore’s multimodal evacuation system [6]. In addition, the Transportation Research Board [7] made a tentative investigation to accommodate emergency evacuation into transit service systems.Effective strategic-level or tactical-level planning is one of the important prerequisites to ensure the success of operational-level emergency evacuation management [8]. For example, transportation authorities often reserve proper emergency lanes on highway networks to promote travel efficiency of emergency vehicles in accident rescues. As a result, a few interesting questions arise. Can we design similar emergency lanes in urban transportation networks? Also, how can we achieve successful planning of emergency lanes in urban transportation systems with limited road infrastructure resources? Moreover, complex traffic conditions in urban transportation networks, such as demand uncertainties and complicated user behaviors, lead to new challenges for the planning of emergency lanes [9, 10]. In this study, we aim to present feasible answers from the perspective of transportation network design. A running example is presented as follows to illustrate the problem of designing emergency vehicle (EV) lanes.


### 1.1. Running Example

Figure1(a) shows a small urban transportation network with three nodes and five directed links. Link 1 has two lanes and each of links 2–5 is deployed with a respective single lane. All lanes have an identical service capacity of 1,000 vehicles/hour. Two origin-destination (OD) pairs, A→C and C→A, are taken into account. There exist two classes of users, auto commuters and EV travelers, in the network. Both classes of network users are assumed to follow Wardrop’s first equilibrium principle that each of the road users expects to find her/his shortest path. The auto travel demands for two OD pairs are, respectively, qA→Cp = 2,000 and qC→Ap = 1,000. The EV travel demands for two OD pairs are set as qA→Ce = 50 and qC→Ae = 60. To improve the travel efficiency of EVs, the network authority attempts to build an EV lane on Link 1. As shown in Figure 1(b), the original Link 1 with two lanes would be then decomposed into a new regular road with one lane (Link 1) and a new EV road with one lane (Link 6). Let us define a link set A≔1,2,…,6. A linear link travel time function is adopted here to evaluate the network performance:(1)tavap,vae=ta0+ta0vap+vaeca,a∈A,where ta0 is the flow-independent travel time of link a,a∈A; vap and vae, respectively, represent link flows of auto users and EV travelers; and the aggregate link flow va=vap+vae. In this example, we set t20=t30=5mins, t40=t50=7mins, and t10=t60=10mins.Figure 1
A three-node network for the running example.
(a)
Original transportation network (b)
New network with an EV laneWe consider and compare the following three EV lane design schemes.Scheme 1.
Take no action of setting EV lanes.Scheme 2.
Build an EV lane to connect the nodesA and C (Link 6 in Figure 1).Scheme 3.
Build an inverse-direction EV lane to accessC→A by reversing the direction of Link 6.Table1 gives the aggregate link flows at biclass user equilibrium and total travel time (TTT) for auto and EV users.Table 1
Results of three design schemes in the running example.
Design schemes Equilibrium link flows TTT for auto users (mins) TTT for EV users (mins) v 1 v 2 v 3 v 4 v 5 v 6 Scheme1 26600 / 17 8250 / 17 1060 8250 / 17 1060 – 60,367.06 2,374.38 Scheme2 13000 / 11 9000 / 11 1060 9000 / 11 1060 50 68,356.36 2,008.20 Scheme3 13300 / 11 9250 / 11 1000 9250 / 11 1000 60 68,181.82 1,740.55As shown in Table1, Scheme 2 sets an EV link and saves about 15.4% golden rescue time, although it makes auto travelers suffer uncomfortable experiences in a more congested network. Compared to Scheme 2, Scheme 3 not only remarkably improves the travel efficiency of EVs but also alleviates the traffic congestion of auto commuters through flexibly incorporating an inverse-direction EV lane. In this sense, Scheme 3 better balances the two design objectives of TTTs associated with auto and EV users. The results of network performance improvements clearly demonstrate the significance of building EV lanes and reveal the value of designing inverse-direction EV lanes in EV network planning.The running example clearly shows that the investigated problem of setting EV lanes might fall into the category of classical discrete transportation network design problems (NDPs), which can be defined as a new link/lane addition that maximizes network performance subject to a series of side constraints (e.g., investment budget constraint and equilibrium constraints of the traffic pattern that characterizes travelers’ route choice behaviors). However, the NDP that optimizes EV lane setting (EV-NDP) also exhibits a few new features and disparities from traditional NDPs (explained later in the next subsection). This study is devoted to proposing an effective approach for the EV-NDP that accounts for more complex and realistic traffic conditions, including application of flexible inverse-direction EV lanes and travel demand uncertainties.


### 1.2. Literature Review

The NDP typically aims to optimize transportation network performance via link capacity expansions (either current road maintenance/improvement or new road addition), while taking into account travelers’ route choice behaviors and financial constraints. After the seminal works of Leblanc [11] and Abdulaal and LeBlanc [12], we have witnessed a blooming of NDPs on model formulations, solution algorithms, and applications over the past half century. Regarding early achievements made in the last century, interested readers are referred to two surveys, Magnanti and Wong [13] and Yang and Bell [14], for a comprehensive literature review and detailed discussions. Meng and Wang [15] presented a structured overview of the NDP literature from 1999 onward and summarized a series of new developments and advances that consider complicated traffic conditions, including uncertainty parameters, traffic dynamics, and multiple optimization objectives. We briefly review the recent achievements of the NDP literature in the following four NDP categories.(1  ) Stochastic NDPs. It is well recognized that network uncertainties with respect to demand fluctuations and road capacity degradation inevitably and inherently exist on transportation networks. Recently, more researchers and practitioners have pointed out that network uncertainty would have significant impacts on network planning and other transportation management decisions, and ignoring this factor might result in suboptimal or even misleading optimization schemes [8, 16, 17]. Chen et al. [18] provided a comprehensive review of stochastic NDPs on various model formulations (e.g., expected-value model, mean-variance model, chance-constrained model, probability model, and min-max model) and proposed an attractive simulation-based genetic algorithm. To properly characterize the effect of network stochasticity, simulation-based sampling (e.g., Monte-Carlo simulation) approximation technique is often used to estimate probabilistic statistics, including expectation and variance, of various system performance measures [8, 18–24]. Alternatively, probabilistic or reliability-based user equilibriums are proposed to incorporate the effect of network uncertainties into travelers’ route choice behaviors in the lower-level traffic assignment problem and similar bilevel stochastic NDP models are then developed [25, 26]. Some stochastic NDP models had also integrated other traffic factors and/or optimization targets (e.g., sustainability, traffic dynamics, and land use). For example, Waller and Ziliaskopoulos [16] formulated a chance-constrained NDP model that takes into account both uncertain demands and dynamic traffic flows. Wang et al. [26] incorporated multiple indices of sustainability, including traffic emission and social equity, into a stochastic NDP.(2  ) Dynamic and Multiperiod NDPs. Time dimension is another important factor in transportation network design. In general, three scales of time dimension are considered: short-term/real-time (seconds), medium-term (days), and long-term (years). Overall, short-term traffic dynamics are attributed to the description of real-time parameter fluctuation in a transportation system, and long-term network decisions across multiple phases help to assess the value of a network improvement plan over time. Specifically, dynamic NDP models have been developed to properly describe real-time traffic dynamics and/or other unsteady-state traffic conditions (e.g., traffic shockwaves propagation, build-up process of queues) [27]. For example, Waller and Ziliaskopoulos [16] proposed a two-stage linear programming NDP model where a system-optimal dynamic traffic assignment and cell transmission loading technique were used to depict dynamic traffic propagation on the network. Karoonsoontawong and Waller [27] presented a robust dynamic NDP model with the objective of minimizing the weighted sum of the expected total travel time and expected risk. In detail, multiperiod NDP formulations were proposed to optimize a design scheme where construction and maintenance of transport infrastructures took place over a long-term horizon. Lo and Szeto [28] proposed the first multiperiod NDP model to examine tradeoffs among several stakeholders: travelers, private toll road operators, and the government. Ukkusuri and Patil [29] made an extension to formulate a multiperiod NDP considering demand uncertainty and demand elasticity, and they found that a multiperiod optimization model could generate more favorable decisions than the conventional single-stage NDP formulation. A more comprehensive literature review and detailed discussions on dynamic and multiperiod NDPs can be found in Meng and Wang [15].(3  ) Multiobjective NDPs. As reported by Yang and Bell [14] and Chen et al. [20], researchers often need to consider multiple design objectives in the NDP, for instance, travel efficiency, emission, social equity, user experience, construction cost, and system reliability. There is a great deal of literature on multiobjective NDPs. For example, Meng and Yang [30] discussed a biobjective NDP regarding travel efficiency and social equity, and Chen and Yang [23] made an extension to analyze the NDP incorporating social equity issue in stochastic networks and developed a chance-constrained optimization model. Wang et al. [26] investigated a chance-constrained multiobjective NDP considering social equity, emission reduction, and network reserve capacity maximization. In general, Pareto optimization method (e.g., [19–21, 31]) and weighted-sum approach (e.g., [8, 22]) are widely applied for multiobjective NDP formulations. For more information on deterministic multiobjective NDPs, refer to review reports of Yang and Bell [14] and Farahani et al. [32], and more literature on stochastic multiobjective NDPs can be found in Chen et al. [18].(4  ) Global Solution Methods. Solving a bilevel NDP model has been seen as a challenge all the time. Global solution algorithms have increasingly received attention in the last decade. For example, based on the Multivariate Taylor Series, Wang and Lo [33] transformed the classical NDP model into equivalent mixed-integer linear programming (MILP), in which a global-optimal solution can be guaranteed. In spirit of the transformed MILP, Luathep et al. [34] made an extension to solve a mixed NDP with both continuous and discrete decision variables. Motivated by a similar idea, Liu and Wang [35] developed an equivalent MILP model for the NDP that took into account logit-based stochastic user equilibrium. The above-mentioned MILP models were all solved by the CPLEX solver. In addition, according to the relationship between user equilibrium and system optimum, Wang et al. [36] presented a system optimum-relaxation based method to solve the discrete NDPs.The EV-NDP in this study exhibits some new features and differs from existing NDPs in the following aspects:(1)
The investigated EV-NDP involves two classes of users: auto commuters and EV travelers. The above literature review indicates that the issue of multiclass users has received rather limited attention. The only NDP work that considered multiclass users was carried out by Dimitriou et al. [37]. Incorporating multiclass users leads to challenges in two aspects, especially when considering uncertain demands: (i) the lower-level traffic assignment problem needs to characterize route choice behaviors of multiple classes of users and resultant equilibrium traffic patterns may not be unique and (ii) the presence of multiclass users probably gives rise to conflict design objectives. It is indeed a challenge to manage tradeoffs among different design objectives associated with multiple classes of users on stochastic networks.(2)
The EV-NDP has a different network reconstruction principle from traditional discrete NDPs. Traditional discrete NDP determines link/lane additions to expand network capacities [14]. The investment is directed to improve link capacity, while network structure accessibility (i.e., topological structure) might keep unchanged. However, in EV-NDP, the total capacity remains unchanged, but certain capacity is removed from original regular links to form new EV links (hereafter, EV links are called derivative links). Take the case in Figure 1 as an example. EV Link 6 is a derivative of Link 1, and setting the EV link results in a capacity reduction of the original Link 1. This can be regarded as link decomposition and capacity reassignment/redeployment. For an EV-NDP scheme, the construction budget is invested to build supporting infrastructures (e.g., pilot signals, guide boards, and communication facilities) that control and manage various vehicles, rather than directly expanding link capacities.(3)
Road access rights are differentiated for two classes of users. As mentioned above, conventional NDPs consider homogeneous commuters and every traveler can freely utilize any link if accessible. This condition is also satisfied in Dimitriou et al. [37], in which multiclass users are taken into account. However, the situation in EV-NDPs differs. EVs can freely use all roads without any restraints, but private cars are not allowed to run on EV lanes. Such differentiated road access rights lead to a challenge in formulating a biclass user equilibrium that describes different travelers’ route choice behaviors.(4)
Inverse-direction EV lanes are introduced for EV-NDPs to fulfill favorable decision schemes. Given the fact that EV trips are independently and separably regulated by the traffic control center during emergency periods, we could have a chance to set inverse-direction EV lanes. This interesting concept makes network design scheme more flexible and attractive.With the above stated disparities, we need to develop new EV-NDP models and propose pertinent solution methods.


### 1.3. Objective and Contributions

The objective of this paper is to propose an effective approach for the EV-NDP with demand uncertainty. To achieve this objective, we firstly propose a network reconstruction method to build a new unified network where both auto and EV users are properly accommodated. We then develop a bilevel optimization model for the EV-NDP on stochastic networks, where the upper level aims to minimize the weighted sum of the expected total travel time of two classes of users and the lower level characterizes route choice behaviors of auto and EV travelers under each demand realization.This study makes the following substantial contributions to the literature:(1)
Our work investigates a new NDP in urban transportation networks that exhibit several disparities from traditional NDPs. To the best of our knowledge, this work addresses a new EV network planning problem and is the first to formulate NDP with multiclass users, which substantially bridges the gap of the existing literature.(2)
We develop a bilevel optimization EV-NDP model considering uncertain demands to find the tradeoff between travel efficiencies for auto and EV travelers. The proposed methods can help transportation authorities make tangible EV network decisions.(3)
We introduce a new concept of inverse-direction EV lane in the EV-NDP. Incorporating inverse-direction EV lanes could promote more flexible and attractive decisions.(4)
We propose a simulation-based genetic algorithm (GA) to solve the developed EV-NDP model and further evaluate the model and computational efficiency of the solution algorithm through numerical experiments. The evaluation results help us better understand the importance of considering inverse-direction EV lanes in EV network planning and the value of simulation-based GA in solving a bilevel stochastic programming model.With all the above explanations, we shall devote Section2 to defining the problem and formulating it as a bilevel optimization model. We then present in Section 3 a simulation-based genetic algorithm to solve the EV-NDP model. Numerical examples are provided in Section 4. The final section, Section 5, concludes the paper.


## 1.1. Running Example

Figure1(a) shows a small urban transportation network with three nodes and five directed links. Link 1 has two lanes and each of links 2–5 is deployed with a respective single lane. All lanes have an identical service capacity of 1,000 vehicles/hour. Two origin-destination (OD) pairs, A→C and C→A, are taken into account. There exist two classes of users, auto commuters and EV travelers, in the network. Both classes of network users are assumed to follow Wardrop’s first equilibrium principle that each of the road users expects to find her/his shortest path. The auto travel demands for two OD pairs are, respectively, qA→Cp = 2,000 and qC→Ap = 1,000. The EV travel demands for two OD pairs are set as qA→Ce = 50 and qC→Ae = 60. To improve the travel efficiency of EVs, the network authority attempts to build an EV lane on Link 1. As shown in Figure 1(b), the original Link 1 with two lanes would be then decomposed into a new regular road with one lane (Link 1) and a new EV road with one lane (Link 6). Let us define a link set A≔1,2,…,6. A linear link travel time function is adopted here to evaluate the network performance:(1)tavap,vae=ta0+ta0vap+vaeca,a∈A,where ta0 is the flow-independent travel time of link a,a∈A; vap and vae, respectively, represent link flows of auto users and EV travelers; and the aggregate link flow va=vap+vae. In this example, we set t20=t30=5mins, t40=t50=7mins, and t10=t60=10mins.Figure 1
A three-node network for the running example.
(a)
Original transportation network (b)
New network with an EV laneWe consider and compare the following three EV lane design schemes.Scheme 1.
Take no action of setting EV lanes.Scheme 2.
Build an EV lane to connect the nodesA and C (Link 6 in Figure 1).Scheme 3.
Build an inverse-direction EV lane to accessC→A by reversing the direction of Link 6.Table1 gives the aggregate link flows at biclass user equilibrium and total travel time (TTT) for auto and EV users.Table 1
Results of three design schemes in the running example.
Design schemes Equilibrium link flows TTT for auto users (mins) TTT for EV users (mins) v 1 v 2 v 3 v 4 v 5 v 6 Scheme1 26600 / 17 8250 / 17 1060 8250 / 17 1060 – 60,367.06 2,374.38 Scheme2 13000 / 11 9000 / 11 1060 9000 / 11 1060 50 68,356.36 2,008.20 Scheme3 13300 / 11 9250 / 11 1000 9250 / 11 1000 60 68,181.82 1,740.55As shown in Table1, Scheme 2 sets an EV link and saves about 15.4% golden rescue time, although it makes auto travelers suffer uncomfortable experiences in a more congested network. Compared to Scheme 2, Scheme 3 not only remarkably improves the travel efficiency of EVs but also alleviates the traffic congestion of auto commuters through flexibly incorporating an inverse-direction EV lane. In this sense, Scheme 3 better balances the two design objectives of TTTs associated with auto and EV users. The results of network performance improvements clearly demonstrate the significance of building EV lanes and reveal the value of designing inverse-direction EV lanes in EV network planning.The running example clearly shows that the investigated problem of setting EV lanes might fall into the category of classical discrete transportation network design problems (NDPs), which can be defined as a new link/lane addition that maximizes network performance subject to a series of side constraints (e.g., investment budget constraint and equilibrium constraints of the traffic pattern that characterizes travelers’ route choice behaviors). However, the NDP that optimizes EV lane setting (EV-NDP) also exhibits a few new features and disparities from traditional NDPs (explained later in the next subsection). This study is devoted to proposing an effective approach for the EV-NDP that accounts for more complex and realistic traffic conditions, including application of flexible inverse-direction EV lanes and travel demand uncertainties.


## 1.2. Literature Review

The NDP typically aims to optimize transportation network performance via link capacity expansions (either current road maintenance/improvement or new road addition), while taking into account travelers’ route choice behaviors and financial constraints. After the seminal works of Leblanc [11] and Abdulaal and LeBlanc [12], we have witnessed a blooming of NDPs on model formulations, solution algorithms, and applications over the past half century. Regarding early achievements made in the last century, interested readers are referred to two surveys, Magnanti and Wong [13] and Yang and Bell [14], for a comprehensive literature review and detailed discussions. Meng and Wang [15] presented a structured overview of the NDP literature from 1999 onward and summarized a series of new developments and advances that consider complicated traffic conditions, including uncertainty parameters, traffic dynamics, and multiple optimization objectives. We briefly review the recent achievements of the NDP literature in the following four NDP categories.(1  ) Stochastic NDPs. It is well recognized that network uncertainties with respect to demand fluctuations and road capacity degradation inevitably and inherently exist on transportation networks. Recently, more researchers and practitioners have pointed out that network uncertainty would have significant impacts on network planning and other transportation management decisions, and ignoring this factor might result in suboptimal or even misleading optimization schemes [8, 16, 17]. Chen et al. [18] provided a comprehensive review of stochastic NDPs on various model formulations (e.g., expected-value model, mean-variance model, chance-constrained model, probability model, and min-max model) and proposed an attractive simulation-based genetic algorithm. To properly characterize the effect of network stochasticity, simulation-based sampling (e.g., Monte-Carlo simulation) approximation technique is often used to estimate probabilistic statistics, including expectation and variance, of various system performance measures [8, 18–24]. Alternatively, probabilistic or reliability-based user equilibriums are proposed to incorporate the effect of network uncertainties into travelers’ route choice behaviors in the lower-level traffic assignment problem and similar bilevel stochastic NDP models are then developed [25, 26]. Some stochastic NDP models had also integrated other traffic factors and/or optimization targets (e.g., sustainability, traffic dynamics, and land use). For example, Waller and Ziliaskopoulos [16] formulated a chance-constrained NDP model that takes into account both uncertain demands and dynamic traffic flows. Wang et al. [26] incorporated multiple indices of sustainability, including traffic emission and social equity, into a stochastic NDP.(2  ) Dynamic and Multiperiod NDPs. Time dimension is another important factor in transportation network design. In general, three scales of time dimension are considered: short-term/real-time (seconds), medium-term (days), and long-term (years). Overall, short-term traffic dynamics are attributed to the description of real-time parameter fluctuation in a transportation system, and long-term network decisions across multiple phases help to assess the value of a network improvement plan over time. Specifically, dynamic NDP models have been developed to properly describe real-time traffic dynamics and/or other unsteady-state traffic conditions (e.g., traffic shockwaves propagation, build-up process of queues) [27]. For example, Waller and Ziliaskopoulos [16] proposed a two-stage linear programming NDP model where a system-optimal dynamic traffic assignment and cell transmission loading technique were used to depict dynamic traffic propagation on the network. Karoonsoontawong and Waller [27] presented a robust dynamic NDP model with the objective of minimizing the weighted sum of the expected total travel time and expected risk. In detail, multiperiod NDP formulations were proposed to optimize a design scheme where construction and maintenance of transport infrastructures took place over a long-term horizon. Lo and Szeto [28] proposed the first multiperiod NDP model to examine tradeoffs among several stakeholders: travelers, private toll road operators, and the government. Ukkusuri and Patil [29] made an extension to formulate a multiperiod NDP considering demand uncertainty and demand elasticity, and they found that a multiperiod optimization model could generate more favorable decisions than the conventional single-stage NDP formulation. A more comprehensive literature review and detailed discussions on dynamic and multiperiod NDPs can be found in Meng and Wang [15].(3  ) Multiobjective NDPs. As reported by Yang and Bell [14] and Chen et al. [20], researchers often need to consider multiple design objectives in the NDP, for instance, travel efficiency, emission, social equity, user experience, construction cost, and system reliability. There is a great deal of literature on multiobjective NDPs. For example, Meng and Yang [30] discussed a biobjective NDP regarding travel efficiency and social equity, and Chen and Yang [23] made an extension to analyze the NDP incorporating social equity issue in stochastic networks and developed a chance-constrained optimization model. Wang et al. [26] investigated a chance-constrained multiobjective NDP considering social equity, emission reduction, and network reserve capacity maximization. In general, Pareto optimization method (e.g., [19–21, 31]) and weighted-sum approach (e.g., [8, 22]) are widely applied for multiobjective NDP formulations. For more information on deterministic multiobjective NDPs, refer to review reports of Yang and Bell [14] and Farahani et al. [32], and more literature on stochastic multiobjective NDPs can be found in Chen et al. [18].(4  ) Global Solution Methods. Solving a bilevel NDP model has been seen as a challenge all the time. Global solution algorithms have increasingly received attention in the last decade. For example, based on the Multivariate Taylor Series, Wang and Lo [33] transformed the classical NDP model into equivalent mixed-integer linear programming (MILP), in which a global-optimal solution can be guaranteed. In spirit of the transformed MILP, Luathep et al. [34] made an extension to solve a mixed NDP with both continuous and discrete decision variables. Motivated by a similar idea, Liu and Wang [35] developed an equivalent MILP model for the NDP that took into account logit-based stochastic user equilibrium. The above-mentioned MILP models were all solved by the CPLEX solver. In addition, according to the relationship between user equilibrium and system optimum, Wang et al. [36] presented a system optimum-relaxation based method to solve the discrete NDPs.The EV-NDP in this study exhibits some new features and differs from existing NDPs in the following aspects:(1)
The investigated EV-NDP involves two classes of users: auto commuters and EV travelers. The above literature review indicates that the issue of multiclass users has received rather limited attention. The only NDP work that considered multiclass users was carried out by Dimitriou et al. [37]. Incorporating multiclass users leads to challenges in two aspects, especially when considering uncertain demands: (i) the lower-level traffic assignment problem needs to characterize route choice behaviors of multiple classes of users and resultant equilibrium traffic patterns may not be unique and (ii) the presence of multiclass users probably gives rise to conflict design objectives. It is indeed a challenge to manage tradeoffs among different design objectives associated with multiple classes of users on stochastic networks.(2)
The EV-NDP has a different network reconstruction principle from traditional discrete NDPs. Traditional discrete NDP determines link/lane additions to expand network capacities [14]. The investment is directed to improve link capacity, while network structure accessibility (i.e., topological structure) might keep unchanged. However, in EV-NDP, the total capacity remains unchanged, but certain capacity is removed from original regular links to form new EV links (hereafter, EV links are called derivative links). Take the case in Figure 1 as an example. EV Link 6 is a derivative of Link 1, and setting the EV link results in a capacity reduction of the original Link 1. This can be regarded as link decomposition and capacity reassignment/redeployment. For an EV-NDP scheme, the construction budget is invested to build supporting infrastructures (e.g., pilot signals, guide boards, and communication facilities) that control and manage various vehicles, rather than directly expanding link capacities.(3)
Road access rights are differentiated for two classes of users. As mentioned above, conventional NDPs consider homogeneous commuters and every traveler can freely utilize any link if accessible. This condition is also satisfied in Dimitriou et al. [37], in which multiclass users are taken into account. However, the situation in EV-NDPs differs. EVs can freely use all roads without any restraints, but private cars are not allowed to run on EV lanes. Such differentiated road access rights lead to a challenge in formulating a biclass user equilibrium that describes different travelers’ route choice behaviors.(4)
Inverse-direction EV lanes are introduced for EV-NDPs to fulfill favorable decision schemes. Given the fact that EV trips are independently and separably regulated by the traffic control center during emergency periods, we could have a chance to set inverse-direction EV lanes. This interesting concept makes network design scheme more flexible and attractive.With the above stated disparities, we need to develop new EV-NDP models and propose pertinent solution methods.


## 1.3. Objective and Contributions

The objective of this paper is to propose an effective approach for the EV-NDP with demand uncertainty. To achieve this objective, we firstly propose a network reconstruction method to build a new unified network where both auto and EV users are properly accommodated. We then develop a bilevel optimization model for the EV-NDP on stochastic networks, where the upper level aims to minimize the weighted sum of the expected total travel time of two classes of users and the lower level characterizes route choice behaviors of auto and EV travelers under each demand realization.This study makes the following substantial contributions to the literature:(1)
Our work investigates a new NDP in urban transportation networks that exhibit several disparities from traditional NDPs. To the best of our knowledge, this work addresses a new EV network planning problem and is the first to formulate NDP with multiclass users, which substantially bridges the gap of the existing literature.(2)
We develop a bilevel optimization EV-NDP model considering uncertain demands to find the tradeoff between travel efficiencies for auto and EV travelers. The proposed methods can help transportation authorities make tangible EV network decisions.(3)
We introduce a new concept of inverse-direction EV lane in the EV-NDP. Incorporating inverse-direction EV lanes could promote more flexible and attractive decisions.(4)
We propose a simulation-based genetic algorithm (GA) to solve the developed EV-NDP model and further evaluate the model and computational efficiency of the solution algorithm through numerical experiments. The evaluation results help us better understand the importance of considering inverse-direction EV lanes in EV network planning and the value of simulation-based GA in solving a bilevel stochastic programming model.With all the above explanations, we shall devote Section2 to defining the problem and formulating it as a bilevel optimization model. We then present in Section 3 a simulation-based genetic algorithm to solve the EV-NDP model. Numerical examples are provided in Section 4. The final section, Section 5, concludes the paper.


## 2. Problem Description and Model Formulation

We first define in this section the EV-NDP. Then, we explain the network reconstruction and formulate the network design model of setting EV lanes (the EV-NDP model).


### 2.1. Problem Description

Consider an urban transportation networkG=(N,A), where N and A are the sets of nodes and links, respectively. There exist two classes of users in the network: auto travelers and EV users. To guarantee the travel efficiency of EVs, the network authority designs an EV-NDP scheme that deploys exclusive EV lanes on the selected roads. To better evaluate the network performance, we take into account random travel demands for both auto and EV users.The network planner not only aims to guarantee the travel efficiency of EVs, but also makes an attempt to improve the travel efficiency of auto travelers. The design objective function is thus to minimize the weighted sum of the expected total travel time of two classes of users on the network. LetA¯ denote the set of candidate links to be deployed with EV lanes. The decision variables can be defined as(2)x^a≔1if  an  EV  lane  in  the  same  direction  of  linkais  deployed0otherwise∀a∈A¯,x~a≔1if  an  EV  lane  in  the  inverse  direction  of  linkais  deployed0otherwise∀a∈A¯.


### 2.2. Notations and Assumptions

The following notations are used throughout the paper.N: the set of nodes in a given urban transportation network G;A: the set of original links in the urban transportation network G;A-: the set of candidate physical links considered in the EV-NDP scheme, A-⊆A;A^: the set of normal-direction EV links, A^⊆A, A^=A-;A~: the set of inverse-direction EV links, A~⊆A, A~=A-;A′: the set of links in the reconstructed unified network, A′≔A∪A^∪A~;W: the set of OD pairs in the transportation network;K: the set of user classes where K≔p,e, where p and e represent private car and emergency vehicle, respectively;Rwk: the set of travel paths associated with kth, k∈K users and OD pair w∈W;Ω: the set of demand realizations;ϖ: a specific demand realization ϖ≔(ϖp,ϖe)∈Ω;S: the set of scenarios to realize the random travel demands;vak,ϖ: traffic flow of kth, k∈K users on link a under demand realization ϖ∈Ω;vaϖ: the aggregate traffic flow on link a under demand realization ϖ∈Ω;fr,wk,ϖ: traffic flow of kth, k∈K users on travel route r, r∈Rwp for OD pair w, w∈W under demand realization ϖ∈Ω;taϖ(·): travel time on link a, a∈A under demand realization ϖ∈Ω;tak,ϖ(·): the link travel time for the kth, k∈K users under demand realization ϖ∈Ω;δa,rk,w: link-path incidence; δa,rk,w=1 if link a is on the travel route r, r∈Rwp for the kth, k∈K users’ trips of OD pair w, w∈W, and 0 otherwise;qwk,ϖ: travel demand associated with kth, k∈K users and OD pair w, w∈W under demand realization ϖ∈Ω;Prϖ: the probability that the demand ϖ∈Ω is realized;ρ: the weighted parameter for the expected total travel time of auto travelers, ρ∈[0,1];B: the total budget for the EV network design;Ca0: the unit lane capacity of physical link a, a∈A;La0: the initial number of lanes of link a, a∈A (before setting EV lanes);Ca: the capacity of link a, a∈A∪A^∪A~;ta0: free-flow travel time of link a, a∈A;αa, βa: two coefficients in the travel time function;ba(·): the construction cost function with respect to decision variables x^a and x~a and the parameter γa.To facilitate model formulation, we make three assumptions below.Assumption 4.
Link travel time functiontaϖ(vaϖ) is a differentiable strictly increasing function with respect to the aggregate link flow vaϖ; that is, dta(va)/dva>0.Assumption 5.
The network design problem takes into account uncertain peak-hour travel demands for both auto and EV users. The random demands for the two classes of users are assumed to follow known probability distributions.Assumption 6.
The EV lanes are assumed to open during the peak-hour period and private cars are not allowed to use any EV lanes. The network planner evaluates the network performance with constantly open EV lanes.Assumption6 implies EVs can either make use of EV lanes or move on regular lanes that eventually generate mixed traffic. We only study an extreme network state with constantly open EV lanes. One may argue that some individuals cannot tolerate building exclusive EV lanes to compete against limited road resources, which may lead to a more congested network. This concern could be greatly alleviated under traffic control management. When the traffic control center receives callings of EVs, private cars are not allowed to use EV lanes; otherwise, the EV lanes are fully free for use. Actually, an interesting extension could be made in the future to simultaneously optimize the EV-NDP scheme and working state setting of EV lanes. The working state setting is to determine an optimal active time of EV lanes (i.e., an optimal proportion to the entire peak-hour period).


### 2.3. Network Reconstruction

Partitioning a regular road generates two independent derivative products, an EV link and a new regular link with reduced capacity (one may have no capacity), and as a result, setting EV lanes changes the network’s topological structures for auto and EV users. It is not surprising that the access rights of private cars might be restricted due to the implementation of EV-NDP schemes. On the one hand, private cars cannot use newly added EV links. On the other hand, private cars that pass new regular links with possibly reduced capacities would face more congested traffic and take longer travel time. In turn, EVs enjoy more road infrastructures and more flexible routing choices due to the newly added EV links on the network. In addition, incorporating inverse-direction EV lanes would make vehicle accessibilities more complicated.To properly address the distinct vehicle accessibilities, we build a unified urban transportation network that can better accommodate two classes of users. In detail, for each candidate linka, a∈A¯ in the network, we add two dummy EV links, one normal-direction EV link a^ and one inverse-direction EV link a~, and map them together as a triplet a,a^,a~. As shown in Figure 2, let us build an EV lane on the candidate Link 1. A corresponding triplet 1,6,7 is thus constructed and A^=6 and A~=7. We then can see that auto drivers can use the roads in the subnetwork G=(N,A), where A=1,2,3,4,5, and EVs commute on the network G=(N,A′), where A′=A∪6,7. Note that if an EV lane is set, Link 1 drops one-lane capacity. Evidently, the unified network is general enough to characterize the differentiated vehicle accessibilities of two classes of users.Figure 2
An illustrative example to reconstruct the network.
(a)
The original network (b)
The reconstructed networkWe next need to identify the association between decision variables and network parameters. We now focus on determining link capacities. Define the EV-NDP scheme asx≔(x^,x~), where x^≔(x^a,a∈A^) and x~≔(x~a,a∈A~). Under a given EV-NDP scheme x, the capacities of links contained in every triplet 〈a,a^,a~〉 can be determined by(3)Ca=Ca0La0-x^a^-x~a~,a∈A¯,Ca′=Ca0x^a^,a^∈A^,Ca~=Ca0x~a~,a~∈A~,where La0 is the number of original lanes of physical link a before setting EV lanes. Here, the unit lane capacity is evenly divided according to the lane number and total link capacity. In practice, on each selected road, at most one EV lane (either normal-direction or inverse-direction EV lane) could be built. We thus have the following feasibility constraint:(4)x^a^+x~a~≤1,for  each  tripleta,a^,a~.To distinguish the access rights of EV linksa∈A^∪A~ for two types of users, we need to redefine the link travel cost functions for the two classes of users:(5)tae,ϖ=taϖvaϖ,for  EV  users,a∈A^∪A~,tap,ϖ=taϖvaϖ+x^aτ¯,for  auto  users,a∈A^,tap,ϖ=taϖvaϖ+x~aτ¯,for  auto  users,a∈A~,where tae,ϖ and tap,ϖ are, respectively, the redefined link travel cost functions of EV users and auto users under demand realization ϖ∈Ω; τ¯ is a sufficiently large road access charge to prohibit auto travelers from using EV links. The different access rights on EV links of the two classes of users could be easily verified in (5).


### 2.4. The EV-NDP Model

In this paper, we aim to investigate an EV-NDP in the strategic planning level, which takes into account long-term uncertainties of auto and EV demands. In spirit of the classical NDP modeling framework, the EV-NDP can be formulated as a bilevel programming model. The upper-level problem aims to minimize the weighted sum of the expected total travel time (ETTT) of auto and EV travelers subject to a series of constraints, while the lower-level problem is to characterize the travelers’ route choice behaviors under given demand realizations. We thus formulate the upper-level problem as follows:(6)gx=minxρ∑∀ϖPr⁡ϖ∑∀avap,ϖtap,ϖvaϖ,x^a,x~a+1-ρ∑∀ϖPr⁡ϖ∑∀avae,ϖtae,ϖvaϖ,x^a,x~a(7)s.t.4,∑a∈A^bax^a,γa+∑a∈A~bax~a,γa≤B,(8)x^a^,x~a~∈0,1,∀a^∈A^,a~∈A~,where ρ∈[0,1] is the weighted parameter. The auxiliary variable vaϖ is the aggregate link flow given that the demand ϖ∈Ω is realized, which is further specified as(9)vaϖ=vap,ϖ+ηvae,ϖ,∀a∈A,where η is the coefficient to convert the EV traffic volumes to passenger car equivalent (PCU) traffic volumes (empirically, η≥1.0).The traffic flowsvap,ϖ and vae,ϖ to the upper level are obtained by solving a biclass user equilibrium in the following lower-level model for each demand realization ϖ∈Ω:(10)minvap,ϖ,vae,ϖ∑∀a∫0vaϖtaϖw,x^a,x~adws.t.∑r∈Rwkfr,wk,ϖ=qwk,ϖ,w∈W,k∈K,∀ϖ,vak,ϖ=∑w∈W∑r∈Rwkfr,wk,ϖδa,rk,w,a∈A′,k∈K,∀ϖ,fr,wk,ϖ≥0,r∈Rwk,k∈K,w∈W,∀ϖ,and the linear constraints (9) of the aggregate link flows.With the strictly increasing link travel time function (Assumption4), the biclass user equilibrium model in the lower level is a strictly convex minimization problem. That is, under given EV-NDP scheme x and demand realization ϖ, solving the user equilibrium generates a unique traffic flow pattern (vap,ϖ,vae,ϖ).As the EV travel demand for each OD pair is far less than auto commuting demand, the magnitude of the ETTT of auto travelers would be significantly greater than that of EV users (see the third case of 68,181.82 versus 1,740.55 in the illustrative example). The huge difference in magnitude leads to difficulty in analyzing the impact of the weighted parameter setting in designing desirable EV-NDP schemes. To better examine the influence of the weighted parameter on balancing the two network performance measures, we introduce two performance benchmarks and redefine a normalized design objective function. Letgp,0 and ge,0, respectively, represent the ETTTs of auto and EV travelers without implementing any EV-NDP schemes, which can be easily obtained by solving the lower-level user equilibrium model (10). Based on the two performance benchmarks gp,0 and ge,0, we thus can rewrite the normalized design objective function g(x) as(11)minxρgp,0∑∀ϖPrϖ∑∀avap,ϖtap,ϖvaϖ,x^a,x~a+1-ρge,0∑∀ϖPrϖ∑∀avae,ϖtae,ϖvaϖ,x^a,x~a.Therefore, the EV-NDP model can be expressed by (11) and is subject to the feasible solution space defined by constraints (4) and (7)–(10).We further give a brief remark on designing the two aforementioned EV-NDP schemes: the scheme with pure normal-direction EV lanes and the scheme with flexible inverse-direction EV lanes. It is not difficult to find that the former scheme is a special case of the latter one. Setting all decision variablesx~=0 of an inverse-direction EV-NDP scheme generates a pure normal-direction EV-NDP scheme.


## 2.1. Problem Description

Consider an urban transportation networkG=(N,A), where N and A are the sets of nodes and links, respectively. There exist two classes of users in the network: auto travelers and EV users. To guarantee the travel efficiency of EVs, the network authority designs an EV-NDP scheme that deploys exclusive EV lanes on the selected roads. To better evaluate the network performance, we take into account random travel demands for both auto and EV users.The network planner not only aims to guarantee the travel efficiency of EVs, but also makes an attempt to improve the travel efficiency of auto travelers. The design objective function is thus to minimize the weighted sum of the expected total travel time of two classes of users on the network. LetA¯ denote the set of candidate links to be deployed with EV lanes. The decision variables can be defined as(2)x^a≔1if  an  EV  lane  in  the  same  direction  of  linkais  deployed0otherwise∀a∈A¯,x~a≔1if  an  EV  lane  in  the  inverse  direction  of  linkais  deployed0otherwise∀a∈A¯.


## 2.2. Notations and Assumptions

The following notations are used throughout the paper.N: the set of nodes in a given urban transportation network G;A: the set of original links in the urban transportation network G;A-: the set of candidate physical links considered in the EV-NDP scheme, A-⊆A;A^: the set of normal-direction EV links, A^⊆A, A^=A-;A~: the set of inverse-direction EV links, A~⊆A, A~=A-;A′: the set of links in the reconstructed unified network, A′≔A∪A^∪A~;W: the set of OD pairs in the transportation network;K: the set of user classes where K≔p,e, where p and e represent private car and emergency vehicle, respectively;Rwk: the set of travel paths associated with kth, k∈K users and OD pair w∈W;Ω: the set of demand realizations;ϖ: a specific demand realization ϖ≔(ϖp,ϖe)∈Ω;S: the set of scenarios to realize the random travel demands;vak,ϖ: traffic flow of kth, k∈K users on link a under demand realization ϖ∈Ω;vaϖ: the aggregate traffic flow on link a under demand realization ϖ∈Ω;fr,wk,ϖ: traffic flow of kth, k∈K users on travel route r, r∈Rwp for OD pair w, w∈W under demand realization ϖ∈Ω;taϖ(·): travel time on link a, a∈A under demand realization ϖ∈Ω;tak,ϖ(·): the link travel time for the kth, k∈K users under demand realization ϖ∈Ω;δa,rk,w: link-path incidence; δa,rk,w=1 if link a is on the travel route r, r∈Rwp for the kth, k∈K users’ trips of OD pair w, w∈W, and 0 otherwise;qwk,ϖ: travel demand associated with kth, k∈K users and OD pair w, w∈W under demand realization ϖ∈Ω;Prϖ: the probability that the demand ϖ∈Ω is realized;ρ: the weighted parameter for the expected total travel time of auto travelers, ρ∈[0,1];B: the total budget for the EV network design;Ca0: the unit lane capacity of physical link a, a∈A;La0: the initial number of lanes of link a, a∈A (before setting EV lanes);Ca: the capacity of link a, a∈A∪A^∪A~;ta0: free-flow travel time of link a, a∈A;αa, βa: two coefficients in the travel time function;ba(·): the construction cost function with respect to decision variables x^a and x~a and the parameter γa.To facilitate model formulation, we make three assumptions below.Assumption 4.
Link travel time functiontaϖ(vaϖ) is a differentiable strictly increasing function with respect to the aggregate link flow vaϖ; that is, dta(va)/dva>0.Assumption 5.
The network design problem takes into account uncertain peak-hour travel demands for both auto and EV users. The random demands for the two classes of users are assumed to follow known probability distributions.Assumption 6.
The EV lanes are assumed to open during the peak-hour period and private cars are not allowed to use any EV lanes. The network planner evaluates the network performance with constantly open EV lanes.Assumption6 implies EVs can either make use of EV lanes or move on regular lanes that eventually generate mixed traffic. We only study an extreme network state with constantly open EV lanes. One may argue that some individuals cannot tolerate building exclusive EV lanes to compete against limited road resources, which may lead to a more congested network. This concern could be greatly alleviated under traffic control management. When the traffic control center receives callings of EVs, private cars are not allowed to use EV lanes; otherwise, the EV lanes are fully free for use. Actually, an interesting extension could be made in the future to simultaneously optimize the EV-NDP scheme and working state setting of EV lanes. The working state setting is to determine an optimal active time of EV lanes (i.e., an optimal proportion to the entire peak-hour period).


## 2.3. Network Reconstruction

Partitioning a regular road generates two independent derivative products, an EV link and a new regular link with reduced capacity (one may have no capacity), and as a result, setting EV lanes changes the network’s topological structures for auto and EV users. It is not surprising that the access rights of private cars might be restricted due to the implementation of EV-NDP schemes. On the one hand, private cars cannot use newly added EV links. On the other hand, private cars that pass new regular links with possibly reduced capacities would face more congested traffic and take longer travel time. In turn, EVs enjoy more road infrastructures and more flexible routing choices due to the newly added EV links on the network. In addition, incorporating inverse-direction EV lanes would make vehicle accessibilities more complicated.To properly address the distinct vehicle accessibilities, we build a unified urban transportation network that can better accommodate two classes of users. In detail, for each candidate linka, a∈A¯ in the network, we add two dummy EV links, one normal-direction EV link a^ and one inverse-direction EV link a~, and map them together as a triplet a,a^,a~. As shown in Figure 2, let us build an EV lane on the candidate Link 1. A corresponding triplet 1,6,7 is thus constructed and A^=6 and A~=7. We then can see that auto drivers can use the roads in the subnetwork G=(N,A), where A=1,2,3,4,5, and EVs commute on the network G=(N,A′), where A′=A∪6,7. Note that if an EV lane is set, Link 1 drops one-lane capacity. Evidently, the unified network is general enough to characterize the differentiated vehicle accessibilities of two classes of users.Figure 2
An illustrative example to reconstruct the network.
(a)
The original network (b)
The reconstructed networkWe next need to identify the association between decision variables and network parameters. We now focus on determining link capacities. Define the EV-NDP scheme asx≔(x^,x~), where x^≔(x^a,a∈A^) and x~≔(x~a,a∈A~). Under a given EV-NDP scheme x, the capacities of links contained in every triplet 〈a,a^,a~〉 can be determined by(3)Ca=Ca0La0-x^a^-x~a~,a∈A¯,Ca′=Ca0x^a^,a^∈A^,Ca~=Ca0x~a~,a~∈A~,where La0 is the number of original lanes of physical link a before setting EV lanes. Here, the unit lane capacity is evenly divided according to the lane number and total link capacity. In practice, on each selected road, at most one EV lane (either normal-direction or inverse-direction EV lane) could be built. We thus have the following feasibility constraint:(4)x^a^+x~a~≤1,for  each  tripleta,a^,a~.To distinguish the access rights of EV linksa∈A^∪A~ for two types of users, we need to redefine the link travel cost functions for the two classes of users:(5)tae,ϖ=taϖvaϖ,for  EV  users,a∈A^∪A~,tap,ϖ=taϖvaϖ+x^aτ¯,for  auto  users,a∈A^,tap,ϖ=taϖvaϖ+x~aτ¯,for  auto  users,a∈A~,where tae,ϖ and tap,ϖ are, respectively, the redefined link travel cost functions of EV users and auto users under demand realization ϖ∈Ω; τ¯ is a sufficiently large road access charge to prohibit auto travelers from using EV links. The different access rights on EV links of the two classes of users could be easily verified in (5).


## 2.4. The EV-NDP Model

In this paper, we aim to investigate an EV-NDP in the strategic planning level, which takes into account long-term uncertainties of auto and EV demands. In spirit of the classical NDP modeling framework, the EV-NDP can be formulated as a bilevel programming model. The upper-level problem aims to minimize the weighted sum of the expected total travel time (ETTT) of auto and EV travelers subject to a series of constraints, while the lower-level problem is to characterize the travelers’ route choice behaviors under given demand realizations. We thus formulate the upper-level problem as follows:(6)gx=minxρ∑∀ϖPr⁡ϖ∑∀avap,ϖtap,ϖvaϖ,x^a,x~a+1-ρ∑∀ϖPr⁡ϖ∑∀avae,ϖtae,ϖvaϖ,x^a,x~a(7)s.t.4,∑a∈A^bax^a,γa+∑a∈A~bax~a,γa≤B,(8)x^a^,x~a~∈0,1,∀a^∈A^,a~∈A~,where ρ∈[0,1] is the weighted parameter. The auxiliary variable vaϖ is the aggregate link flow given that the demand ϖ∈Ω is realized, which is further specified as(9)vaϖ=vap,ϖ+ηvae,ϖ,∀a∈A,where η is the coefficient to convert the EV traffic volumes to passenger car equivalent (PCU) traffic volumes (empirically, η≥1.0).The traffic flowsvap,ϖ and vae,ϖ to the upper level are obtained by solving a biclass user equilibrium in the following lower-level model for each demand realization ϖ∈Ω:(10)minvap,ϖ,vae,ϖ∑∀a∫0vaϖtaϖw,x^a,x~adws.t.∑r∈Rwkfr,wk,ϖ=qwk,ϖ,w∈W,k∈K,∀ϖ,vak,ϖ=∑w∈W∑r∈Rwkfr,wk,ϖδa,rk,w,a∈A′,k∈K,∀ϖ,fr,wk,ϖ≥0,r∈Rwk,k∈K,w∈W,∀ϖ,and the linear constraints (9) of the aggregate link flows.With the strictly increasing link travel time function (Assumption4), the biclass user equilibrium model in the lower level is a strictly convex minimization problem. That is, under given EV-NDP scheme x and demand realization ϖ, solving the user equilibrium generates a unique traffic flow pattern (vap,ϖ,vae,ϖ).As the EV travel demand for each OD pair is far less than auto commuting demand, the magnitude of the ETTT of auto travelers would be significantly greater than that of EV users (see the third case of 68,181.82 versus 1,740.55 in the illustrative example). The huge difference in magnitude leads to difficulty in analyzing the impact of the weighted parameter setting in designing desirable EV-NDP schemes. To better examine the influence of the weighted parameter on balancing the two network performance measures, we introduce two performance benchmarks and redefine a normalized design objective function. Letgp,0 and ge,0, respectively, represent the ETTTs of auto and EV travelers without implementing any EV-NDP schemes, which can be easily obtained by solving the lower-level user equilibrium model (10). Based on the two performance benchmarks gp,0 and ge,0, we thus can rewrite the normalized design objective function g(x) as(11)minxρgp,0∑∀ϖPrϖ∑∀avap,ϖtap,ϖvaϖ,x^a,x~a+1-ρge,0∑∀ϖPrϖ∑∀avae,ϖtae,ϖvaϖ,x^a,x~a.Therefore, the EV-NDP model can be expressed by (11) and is subject to the feasible solution space defined by constraints (4) and (7)–(10).We further give a brief remark on designing the two aforementioned EV-NDP schemes: the scheme with pure normal-direction EV lanes and the scheme with flexible inverse-direction EV lanes. It is not difficult to find that the former scheme is a special case of the latter one. Setting all decision variablesx~=0 of an inverse-direction EV-NDP scheme generates a pure normal-direction EV-NDP scheme.


## 3. Solution Algorithms

As the proposed bilevel programming model is nonlinear and nonconvex in nature, it would be rather difficult to solve the EV-NDP model using conventional solution methods. Moreover, the incorporation of demand uncertainties results in a higher computational complexity, which makes the model more intractable to be solved. At early stages, a few greedy/descent search algorithms were developed to solve the NDPs with deterministic demands, including the Hooke–Jeeves algorithm [12], sequential quadratic programming method [38], augmented Lagrangian algorithm [39], and gradient/conjugate gradient/quasi-Newton projection methods [40]. In the last decade, more and more researchers have paid attention to global optimization methods, including mixed-integer linear programming (MILP) method [33–35], benders decomposition algorithm [41], and system optimum-relaxation based method [36].However, the above-mentioned solution algorithms are inapplicable for solving the EV-NDP model. On the one hand, we can hardly use descent search algorithms to solve the EV-NDP model due to two reasons: first, the search space is large and highly complex; second, it is intractable to obtain effective descent gradients for a discrete bilevel optimization model, especially in the presence of stochastic demands. On the other hand, incorporating demand randomness requires us to address a large number of demand realizations and consequently makes it almost impossible to design a global optimization method.Fortunately, metaheuristic solution algorithms can provide us with alternative options. Previous studies clearly indicated that the genetic algorithm (GA) is one of the most prevailing and recommendable tools to solve NDP models [18, 42]. For example, Chen and Yang [23] proposed a simulation-based GA procedure to solve expected-value and chance-constraint NDP models that considered both spatial equity and demand uncertainty. Chootinan et al. [24] presented a bilevel capacity reliability-based network design model, where the upper-level model was solved by Monte-Carlo simulation-based GA and probit-based stochastic user equilibrium in the lower level was solved by the method of successive averages. Chen et al. [20] and Chen and Xu [19] applied the GA to solve multiobjective NDPs with demand uncertainties. Ukkusuri et al. [8] developed a GA to solve a robust NDP, in which the design objective was to minimize the weighted sum of the expected total travel time and related standard variance. Sharma et al. [21, 22] applied the GA to solve similar mean-variance form robust NDPs with random demands. More information about the GA in solving stochastic NDPs can be found in Chen et al. [18].In this paper, we also design a simulation-based GA to solve the proposed EV-NDP model. Figure3 presents the explicit framework of the proposed GA. It works with four submodules: initialization module, simulation module, traffic assignment module, and GA operation module. Each of them is elaborated below:(1)
Initialization module: this module translates the decision variables to GA codes, sets input parameters of experiments, and randomly generates an initial population.(2)
Simulation module: it is used to produce a series of demand realizations according to given demand probability distributions. Next, combined with traffic patterns obtained by solving the lower-level user equilibrium under each demand realization, we use Monte-Carlo simulation to estimate the design objective function using sample average approximation (SAA). Interested readers can refer to Sharma et al. [22] for other attractive sampling approximation techniques, for example, randomized quasi-Monte-Carlo sampling and single-point approximation.(3)
Traffic assignment module: it characterizes the users’ route choice behaviors associated with the given demand realization and EV-NDP scheme. The widely applied Frank-Wolfe algorithm is used to solve the lower-level traffic assignment problem.(4)
GA operation module: this module aims to generate the next trial solution via various genetic operations. Applying a series of GA operations, including reproduction, crossover, and mutation, produces better offspring. More details could be found in Ukkusuri et al. [8].Figure 3
The simulation-based genetic algorithm for the EV-NDP model.Algorithm1 presents the detailed procedure to solve the developed EV-NDP model.Algorithm 1:Simulation-based genetic algorithm.
Input: The EV-NDP model, network parameters, sample size MS, population size MP, generation size MG, crossover rate pc, mutation rate pm, construction budget B, and etc. Output: the optimal objective value (ETTT) and optimal solution x∗. ( 1 ) Code variables and initialize the population; ( 2 )  for  every generation  i=1  to  i=MG  do ( 3 )  for  every population  j=1  to  j=MP  do ( 4 ) Decode the population and update link capacities; ( 5 )  for  every demand sample  ϖ=1  to  ϖ=MS  do ( 6 ) Generate random travel demands for two classes of users; ( 7 ) Solve the bi-class user equilibrium traffic assignment with demand realization ϖ by Frank-Wolfe Algorithm and consequently obtain the link flow patterns; ( 8 ) Compute the design objective function; ( 9 ) Compute the fitness function; ( 10 ) Generate better offspring by genetic operations, including reproduction, crossover and mutation. ( 11 ) End the algorithm and obtain the solution.In solving the EV-NDP model, we take an improved strategy of population initialization to filter out infeasible populations that violate the feasibility constraint (4) and budget constraint (7). Furthermore, a population updating strategy (e.g., reinitializing population frequently) can be used to avoid a local optimum as far as possible.


## 4. Numerical Examples

In this section, numerical examples are provided to evaluate the developed EV-NDP model and the computation efficiency of the proposed algorithm. The programming code is compiled by VC++ and run on Windows 7 system with the following attributes: Intel Core i5-2400 3.1GHz×2 and 4 GB RAM.As shown in Figure4, a four-node urban transportation network with nine links is used for the numerical experiments. Link parameter settings including the number of lanes, unit lane capacity, and free-flow travel time are provided in Table 2.Table 2
The detailed link parameter setting.
Link number 1 2 3 4 5 6 7 8 9 C a 0 1200 1000 1000 800 800 1200 1000 1000 1200 L a 0 3 2 2 2 2 3 2 2 3 t a 0 10 8 8 4 4 6 5 5 12Figure 4
The used transportation network.Six OD pairs are taken into account. For each OD pair, the auto and EV travel demands are all assumed to follow given normal distributions. Detailed parameters of respective demand distributions are given in Table3.Table 3
Travel demand distributions for two classes of users.
OD pair 1 → 2 1 → 4 2 → 1 2 → 4 4 → 1 4 → 2 Demand for cars ( 1000,10 0 2 ) ( 2000,10 0 2 ) ( 1000,10 0 2 ) ( 2200,10 0 2 ) ( 1800,10 0 2 ) ( 2100,10 0 2 ) Demand for EVs ( 50 , 5 2 ) ( 50 , 5 2 ) ( 50 , 5 2 ) ( 50 , 5 2 ) ( 50 , 5 2 ) ( 50 , 5 2 ) Note. Numerics (a,b2) are the mean and variance of a given normal distribution for each OD pair.The link travel time function (BPR function) is given below:(12)taϖvaϖ=ta01+0.15vaϖCa4,where vaϖ=vap,ϖ+1.5vae,ϖ.A linear function is used to determine the construction cost of setting an EV lane:(13)ba=10.0+10.0ta0x^a,a∈A^,10.0+10.0ta0x~a,a∈A~.Other parameters used in the GA procedure are listed in Table4.Table 4
Other parameters used in the numerical experiments.
Population size M P = 10 for the normal-direction EV-NDP scheme M P = 30 for the inverse-direction EV-NDP scheme Generation size M G = 50 for the normal-direction EV-NDP scheme M G = 80 for the inverse-direction EV-NDP scheme  Probability of crossover 0.75 Total investment budget B = 400 Probability of mutation 0.1 Sample size M S = 500For both normal-direction and inverse-direction EV-NDP schemes, we have tested the experiments by changing the weighted parameter in the design objective function (an equal increment of 0.1 is set). In addition to applying the simulation-based GA, we solve the normal-direction EV-NDP schemes by the enumeration method in order to compare the solution quality and computational efficiency. Solving the lower-level traffic assignment problem of the normal-direction EV-NDP generates the benchmark measures wheregp,0 = 113,779.42 and ge,0 = 3,395.93 (note that solving the user equilibrium model of the inverse-direction EV-NDP yields gp,0 = 113,777.49 and ge,0 = 3,395.63; the small differences are partly due to computational accuracy and partly resulted from two different feasible path sets between the normal-direction and inverse-direction networks).


### 4.1. Computational Results

Table5 shows the outcomes (including the optimal solution—the “optimal solution” obtained by the simulation-based GA may not be global-optimal but only local-optimal—ETTTs for two classes of users, the objective value, and computational time) of the normal-direction EV-NDP schemes. Table 6 provides the computed results of the inverse-direction EV-NDP schemes. As shown in Table 5, for the normal-direction EV-NDPs, the simulation-based GA and enumeration method generate identical results except for the case with ρ=0.1. In what follows, we make further discussions on the computed results and expect to gain some insightful findings.Table 5
The outcomes for the normal-direction EV-NDP schemes.
Weighting factorρ Optimal solutionx^∗ ETTT for cars ETTT for EVs Objective valueg∗(x^∗) (%) Computational time (s) 0.00 ( 0,0 , 1,0 , 1,1 , 1,0 , 1 ) 151,763.02 3,129.54 92.156 1,452.29 0.10 ( 0,0 , 1,0 , 1,0 , 1,1 , 1 ) 148,924.72 3,140.92 96.331 1,471.84 (0, 0, 1, 1, 1, 1, 1, 1, 0) 146,643.41 3,146.02 96.265 1,579.12 0.20 ( 0,0 , 1,0 , 0,0 , 0,1 , 1 ) 119,101.05 3,322.89 99.215 1,522.38 0.30 ( 0,0 , 1,0 , 0,0 , 0,1 , 0 ) 117,732.08 3,333.16 99.748 1,526.89 0.40 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,467.46 0.50 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,481.86 0.60 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,567.94 0.70 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,542.40 0.80 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,521.61 0.90 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,508.46 1.00 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,529.05 Note. For the experiment with ρ=0.1, the computed results solved by the enumeration method are underlined, while the results obtained by the simulation-based GA are presented without underline.Table 6
The outcomes for the inverse-direction EV-NDP schemes.
Weighting factorρ Optimal solutionx∗=(x^∗,x~∗) ETTT for cars ETTT for EVs Objective valueg∗(x^∗) (%) Computational time (s) 0.00 x ^ ∗ = 0, x~∗=(1,0,0,1,0,1,0,0,1) 116,502.41 2,884.32 84.935 25,512.50 0.10 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 86.548 25,836.17 0.20 x ^ ∗ = 0, x~∗=(1,1,0,0,0,1,0,0,1) 116,230.05 2,884.32 88.379 26,180.31 0.30 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 89.775 25,314.65 0.40 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 91.389 25,788.18 0.50 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 93.002 25,779.22 0.60 x ^ ∗ = 0, x~∗=(1,1,0,0,0,1,0,0,0) 114,740.04 2,922.68 94.932 25,747.76 0.70 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 96.229 25,739.35 0.80 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,0) 113,764.05 2,970.79 97.485 25,749.86 0.90 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,0) 113,764.05 2,970.79 98.736 26,324.82 1.00 x ^ ∗ = 0, x~∗=(1,0,0,0,0,0,0,0,0) 113,759.80 3,286.19 99.983 25,750.46


### 4.2. The Value of Incorporating Inverse-Direction EV Lanes

We first analyze the importance of incorporating inverse-direction EV lanes. The objective values of the normal-direction and inverse-direction EV-NDPs with different weighted parameter settings are illustrated in Figure5 for an easy comparison. It is not surprising to realize that, with a given weighted parameter, the scheme considering inverse-direction EV lanes performs better than the corresponding normal-direction EV-NDP in terms of the design objective value since the latter scheme is a special case of the inverse-direction EV-NDP. Take the experiment with ρ=0.2 as an example. We obtain a global-optimal objective value g∗=99.215% by solving the normal-direction EV-NDP scheme with the enumeration method. Solving the inverse-direction EV-NDP scheme by the simulation-based GA generates an optimal objective value g∗=88.379%. Although the solution of the inverse-direction EV-NDP scheme leads to a local optimum (the solution (x^∗=0,x~∗=(1,0,0,0,0,1,0,0,1)) yields a better objective value of g∗=88.162%), it still remarkably improves the travel efficiencies of both classes of users. It can be seen from Tables 5 and 6 that both classes of users benefit from incorporating inverse-direction EV lanes: the ETTT of private cars drops from 119,101.05 to 116,230.05 and that of EVs decreases from 3,322.89 to 2,884.32. In summary, incorporating inverse-direction EV lanes in the EV-NDP offers a chance to design more attractive EV-NDP schemes.Figure 5
Comparison between the normal-direction and inverse-direction EV-NDPs.


### 4.3. The Impact of Weighted Parameter Setting

As shown in Figure5, for both EV-NDP schemes, the design objective values climb steadily with an increasing weighted factor. As stated in the Introduction, the nature of EV-NDP is to reallocate the road capacity to each class of travelers and the total capacity on each link is fixed. Consequently, allocating part of the road capacity to EVs often makes private cars suffer a more congested traffic condition unless a paradox occurs (this indeed happens and will be discussed later). In this sense, balancing two conflicting objectives of minimizing ETTTs of private cars and EVs eventually creates, to a certain extent, a zero-sum game. As shown in Table 5, the ETTT for EVs falls monotonically when a higher priority is given to enhance EVs’ travel efficiency. On the contrary, the ETTT for private cars rises as the weighted parameter increases. A similar trend can be found in the outcomes of inverse-direction EV-NDP schemes, although it is not strictly monotone due to local-optimal solutions attained by the simulation-based GA. Overall, the developed weighted-sum EV-NDP model is capable of finding an appropriate tradeoff between two conflicting design objectives for two classes of users.


### 4.4. Solution Convergence, Quality, and Computational Efficiency

We now turn to analyze the solution convergence, quality, and computational efficiency of the proposed algorithm. Figure6 clearly displays the convergence of the proposed simulation-based GA in solving both the normal-direction and the inverse-direction EV-NDPs. With no doubt, the simulation-based GA is able to obtain local-optimal solutions.Figure 6
Examples to illustrate the convergence of the simulation-based GA.
(a)
Normal-direction EV-NDP withρ=0.1 (b)
Inverse-direction EV-NDP withρ=0.3Let us next examine the solution quality. The results of normal-direction EV-NDPs solved by the simulation-based GA and enumeration method are depicted in Figure7 for a straightforward comparison. It can be observed that the outcomes from the simulation-based GA are identical to those from the enumeration method, except for the case with ρ=0.1. Actually, for the experiment with ρ=0.1, the optimal objective values solved by the two methods are very close (96.331 versus 96.265). In general, the simulation-based GA can generate satisfactory solutions, although the tested experiments might be network-specific.Figure 7
Comparison of the results of normal-direction EV-NDPs solved by the simulation-based GA and enumeration method.We then focus on the computational efficiency of the proposed algorithm. It can be seen in Table5 that it takes about 25 minutes to solve the normal-direction EV-NDP model by the simulation-based GA (500=10×50 solutions are examined). The enumeration method consumes a bit longer computation time since all 512=29 combinatorial solutions are computed and compared. A heavy computational burden results from repeatedly solving the lower-level traffic assignment problem until all demand scenarios are realized. The larger the sample size is, the heavier the computational burden is. For each given solution x, it approximately takes 10 seconds to calculate the objective value for the inverse-direction EV-NDP model. Recall that more than 262,144 = 218 solutions need to be evaluated. Although the number of feasible solutions could be reduced according to the feasibility constraints (4) and (7), the enumeration searching space is still unbelievably large. It is expected to take more than hundreds of hours to solve the inverse-direction EV-NDP model by the enumeration method. Therefore, it is almost intractable to do that even for a small-scale network. Fortunately, the simulation-based GA offers an alternative to find satisfactory solutions for the inverse-direction EV-NDPs within 7.5 hours.


### 4.5. An Insightful Finding of the EV-NDP

Last but not least, we obtain an interesting insight into the inverse-direction EV-NDP withρ=1.0. Setting the weighted parameter ρ equal to 1.0 implies that the network planner is only concerned about the ETTT of private cars. To achieve such objective, it is always believed that the decision-maker needs to provide infrastructure resources to private cars as much as possible. Usually, building EV lanes diverts the resources away and would make auto users suffer a loss. However, a paradox occurs. As shown in Table 6, the optimal solution suggests constructing an inverse-direction EV lane on Link 1. As a result, both EVs and private cars benefit from this decision, where the ETTT for EVs decreases from 3,395.63 to 3,286.19 and that for private cars goes down slightly from 113,777.49 to 113,759.80. An inherent reason is that adding an inverse-direction EV lane makes Link 1 become more congested, and thus some auto users (OD pair 2→1) choose to switch to other travel paths/links. This consequently drives the user equilibrium approaching system optimum. It is not difficult to verify that the EV travelers of OD pair 1→2 benefit significantly from the new added EV lane.The lesson we have learned here is that the objectives to maximize the travel efficiencies of both classes of users are not always conflicting but could be reconcilable in some cases. Meanwhile, the occurrence of the above paradox has highlighted the value of incorporating inverse-direction EV lanes once again (such a win-win case does not happen in normal-direction EV-NDPs).


## 4.1. Computational Results

Table5 shows the outcomes (including the optimal solution—the “optimal solution” obtained by the simulation-based GA may not be global-optimal but only local-optimal—ETTTs for two classes of users, the objective value, and computational time) of the normal-direction EV-NDP schemes. Table 6 provides the computed results of the inverse-direction EV-NDP schemes. As shown in Table 5, for the normal-direction EV-NDPs, the simulation-based GA and enumeration method generate identical results except for the case with ρ=0.1. In what follows, we make further discussions on the computed results and expect to gain some insightful findings.Table 5
The outcomes for the normal-direction EV-NDP schemes.
Weighting factorρ Optimal solutionx^∗ ETTT for cars ETTT for EVs Objective valueg∗(x^∗) (%) Computational time (s) 0.00 ( 0,0 , 1,0 , 1,1 , 1,0 , 1 ) 151,763.02 3,129.54 92.156 1,452.29 0.10 ( 0,0 , 1,0 , 1,0 , 1,1 , 1 ) 148,924.72 3,140.92 96.331 1,471.84 (0, 0, 1, 1, 1, 1, 1, 1, 0) 146,643.41 3,146.02 96.265 1,579.12 0.20 ( 0,0 , 1,0 , 0,0 , 0,1 , 1 ) 119,101.05 3,322.89 99.215 1,522.38 0.30 ( 0,0 , 1,0 , 0,0 , 0,1 , 0 ) 117,732.08 3,333.16 99.748 1,526.89 0.40 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,467.46 0.50 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,481.86 0.60 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,567.94 0.70 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,542.40 0.80 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,521.61 0.90 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,508.46 1.00 ( 0,0 , 0,0 , 0,0 , 0,0 , 0 ) 113,779.42 3,395.93 100.00 1,529.05 Note. For the experiment with ρ=0.1, the computed results solved by the enumeration method are underlined, while the results obtained by the simulation-based GA are presented without underline.Table 6
The outcomes for the inverse-direction EV-NDP schemes.
Weighting factorρ Optimal solutionx∗=(x^∗,x~∗) ETTT for cars ETTT for EVs Objective valueg∗(x^∗) (%) Computational time (s) 0.00 x ^ ∗ = 0, x~∗=(1,0,0,1,0,1,0,0,1) 116,502.41 2,884.32 84.935 25,512.50 0.10 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 86.548 25,836.17 0.20 x ^ ∗ = 0, x~∗=(1,1,0,0,0,1,0,0,1) 116,230.05 2,884.32 88.379 26,180.31 0.30 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 89.775 25,314.65 0.40 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 91.389 25,788.18 0.50 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 93.002 25,779.22 0.60 x ^ ∗ = 0, x~∗=(1,1,0,0,0,1,0,0,0) 114,740.04 2,922.68 94.932 25,747.76 0.70 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,1) 114,996.07 2,884.32 96.229 25,739.35 0.80 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,0) 113,764.05 2,970.79 97.485 25,749.86 0.90 x ^ ∗ = 0, x~∗=(1,0,0,0,0,1,0,0,0) 113,764.05 2,970.79 98.736 26,324.82 1.00 x ^ ∗ = 0, x~∗=(1,0,0,0,0,0,0,0,0) 113,759.80 3,286.19 99.983 25,750.46


## 4.2. The Value of Incorporating Inverse-Direction EV Lanes

We first analyze the importance of incorporating inverse-direction EV lanes. The objective values of the normal-direction and inverse-direction EV-NDPs with different weighted parameter settings are illustrated in Figure5 for an easy comparison. It is not surprising to realize that, with a given weighted parameter, the scheme considering inverse-direction EV lanes performs better than the corresponding normal-direction EV-NDP in terms of the design objective value since the latter scheme is a special case of the inverse-direction EV-NDP. Take the experiment with ρ=0.2 as an example. We obtain a global-optimal objective value g∗=99.215% by solving the normal-direction EV-NDP scheme with the enumeration method. Solving the inverse-direction EV-NDP scheme by the simulation-based GA generates an optimal objective value g∗=88.379%. Although the solution of the inverse-direction EV-NDP scheme leads to a local optimum (the solution (x^∗=0,x~∗=(1,0,0,0,0,1,0,0,1)) yields a better objective value of g∗=88.162%), it still remarkably improves the travel efficiencies of both classes of users. It can be seen from Tables 5 and 6 that both classes of users benefit from incorporating inverse-direction EV lanes: the ETTT of private cars drops from 119,101.05 to 116,230.05 and that of EVs decreases from 3,322.89 to 2,884.32. In summary, incorporating inverse-direction EV lanes in the EV-NDP offers a chance to design more attractive EV-NDP schemes.Figure 5
Comparison between the normal-direction and inverse-direction EV-NDPs.


## 4.3. The Impact of Weighted Parameter Setting

As shown in Figure5, for both EV-NDP schemes, the design objective values climb steadily with an increasing weighted factor. As stated in the Introduction, the nature of EV-NDP is to reallocate the road capacity to each class of travelers and the total capacity on each link is fixed. Consequently, allocating part of the road capacity to EVs often makes private cars suffer a more congested traffic condition unless a paradox occurs (this indeed happens and will be discussed later). In this sense, balancing two conflicting objectives of minimizing ETTTs of private cars and EVs eventually creates, to a certain extent, a zero-sum game. As shown in Table 5, the ETTT for EVs falls monotonically when a higher priority is given to enhance EVs’ travel efficiency. On the contrary, the ETTT for private cars rises as the weighted parameter increases. A similar trend can be found in the outcomes of inverse-direction EV-NDP schemes, although it is not strictly monotone due to local-optimal solutions attained by the simulation-based GA. Overall, the developed weighted-sum EV-NDP model is capable of finding an appropriate tradeoff between two conflicting design objectives for two classes of users.


## 4.4. Solution Convergence, Quality, and Computational Efficiency

We now turn to analyze the solution convergence, quality, and computational efficiency of the proposed algorithm. Figure6 clearly displays the convergence of the proposed simulation-based GA in solving both the normal-direction and the inverse-direction EV-NDPs. With no doubt, the simulation-based GA is able to obtain local-optimal solutions.Figure 6
Examples to illustrate the convergence of the simulation-based GA.
(a)
Normal-direction EV-NDP withρ=0.1 (b)
Inverse-direction EV-NDP withρ=0.3Let us next examine the solution quality. The results of normal-direction EV-NDPs solved by the simulation-based GA and enumeration method are depicted in Figure7 for a straightforward comparison. It can be observed that the outcomes from the simulation-based GA are identical to those from the enumeration method, except for the case with ρ=0.1. Actually, for the experiment with ρ=0.1, the optimal objective values solved by the two methods are very close (96.331 versus 96.265). In general, the simulation-based GA can generate satisfactory solutions, although the tested experiments might be network-specific.Figure 7
Comparison of the results of normal-direction EV-NDPs solved by the simulation-based GA and enumeration method.We then focus on the computational efficiency of the proposed algorithm. It can be seen in Table5 that it takes about 25 minutes to solve the normal-direction EV-NDP model by the simulation-based GA (500=10×50 solutions are examined). The enumeration method consumes a bit longer computation time since all 512=29 combinatorial solutions are computed and compared. A heavy computational burden results from repeatedly solving the lower-level traffic assignment problem until all demand scenarios are realized. The larger the sample size is, the heavier the computational burden is. For each given solution x, it approximately takes 10 seconds to calculate the objective value for the inverse-direction EV-NDP model. Recall that more than 262,144 = 218 solutions need to be evaluated. Although the number of feasible solutions could be reduced according to the feasibility constraints (4) and (7), the enumeration searching space is still unbelievably large. It is expected to take more than hundreds of hours to solve the inverse-direction EV-NDP model by the enumeration method. Therefore, it is almost intractable to do that even for a small-scale network. Fortunately, the simulation-based GA offers an alternative to find satisfactory solutions for the inverse-direction EV-NDPs within 7.5 hours.


## 4.5. An Insightful Finding of the EV-NDP

Last but not least, we obtain an interesting insight into the inverse-direction EV-NDP withρ=1.0. Setting the weighted parameter ρ equal to 1.0 implies that the network planner is only concerned about the ETTT of private cars. To achieve such objective, it is always believed that the decision-maker needs to provide infrastructure resources to private cars as much as possible. Usually, building EV lanes diverts the resources away and would make auto users suffer a loss. However, a paradox occurs. As shown in Table 6, the optimal solution suggests constructing an inverse-direction EV lane on Link 1. As a result, both EVs and private cars benefit from this decision, where the ETTT for EVs decreases from 3,395.63 to 3,286.19 and that for private cars goes down slightly from 113,777.49 to 113,759.80. An inherent reason is that adding an inverse-direction EV lane makes Link 1 become more congested, and thus some auto users (OD pair 2→1) choose to switch to other travel paths/links. This consequently drives the user equilibrium approaching system optimum. It is not difficult to verify that the EV travelers of OD pair 1→2 benefit significantly from the new added EV lane.The lesson we have learned here is that the objectives to maximize the travel efficiencies of both classes of users are not always conflicting but could be reconcilable in some cases. Meanwhile, the occurrence of the above paradox has highlighted the value of incorporating inverse-direction EV lanes once again (such a win-win case does not happen in normal-direction EV-NDPs).


## 5. Concluding Remarks

In this paper, we studied a new network planning problem to build exclusive EV lanes in an urban transportation network, which falls into the category of discrete NDPs. The investigated EV-NDP differs from the classical NDPs in several aspects: it involves two classes of network users and has a distinct network reconstruction principle; auto and EV travelers have differentiated road access rights on the network. To make more tangible and practical planning decisions, random travel demands that follow known probability distributions are taken into account in the EV-NDP. Moreover, inverse-direction EV lanes are incorporated to promote flexible and attractive EV-NDP schemes.A unified network framework is introduced to properly describe the differentiated road access rights of private cars and EVs. We then formulate a bilevel programming EV-NDP model. In the upper level, the network planner intends to minimize the weighted sum of the ETTTs of auto and EV travelers by setting either normal-direction or inverse-direction EV lanes on selected roads. In the lower level, a biclass user equilibrium model is used to characterize travelers’ route choice behaviors under given demand realizations. To handle the stochasticity of travel demands, a simulation-based GA is proposed to solve the developed EV-NDP model. We further evaluate the developed EV-NDP model and the proposed solution algorithm through numerical examples.We arrive at a bunch of useful and interesting observations: incorporating inverse-direction EV lanes could help to design more favorable EV-NDP schemes; the simulation-based GA is capable of finding satisfactory solutions and guaranteeing a high solution quality with an acceptable computational burden; an interesting paradox is found where building EV lanes might contribute to reduced traffic congestion for both auto and EV travelers. This observation tells us that the intuitively conflicting objectives of the two classes of users could be reconcilable in some cases.This research makes an effort to fill in the gap in the literature by developing a bilevel optimization model for the EV-NDP. We can further make some interesting extensions. Firstly, a joint optimization of the EV-NDP scheme and working state setting of EV lanes can be considered so that more reasonable planning decisions could be obtained. Secondly, we analyze a static analytical framework of the EV-NDP in this paper. Extension to a dynamic environment is a more appealing research direction.



---

*Source: 2524838-2018-03-08.xml*