Optimal Road Toll Design from the Perspective of Sustainable Development

This paper investigates the optimal road toll design problem from the perspective of sustainable development at network-wide level. In this paper the sustainable development level of transportation system is quantitatively described with the total vehicular emission, total fuel consumption, and total travel time in the network. In order to simultaneously consider the impacts of all these three indicators on sustainability of transportation system, we integrate them into a sustainable development index (SD-index) by a linear combination, and then we establish the corresponding bilevel optimization model. The upper level problem is the network toll design problem to maximize the SD-index from the viewpoint of traffic managers, and the lower level problem is to depict travelers’ route choice behavior under a certain road toll scheme. Finally, a combined genetic algorithm and gradient projection algorithm (GA-GP) is used to solve the bilevel model, in which the GP algorithm solves the traffic assignment problem with road toll scheme in the lower level. In order to verify the proposed model and algorithm, we take the Nguyen-Dupuis network for the numerical example, and the computing results show that the model and algorithm are effective and efficient.


Introduction
Traffic congestion is a global problem that almost everyone living in big cities has to face, which causes a great waste of time for the urban commuters.At the same time, the rapid development of transportation, especially the current unrestrictive use of private cars, also brings about the severe air pollution and energy crisis problem.It is well known that all these traffic problems seriously restrict the further development for the cities all over the world.And therefore, there is clearly an urgent need for effective measures and policies to combat all these issues caused by the unsustainable transportation systems and then to develop a sustainable, low-carbon, and energy-saving urban transportation systems.
Road toll pricing, which is a type of traffic demand management measure, is widely recognized as a useful tool for alleviating traffic congestion and reducing vehicular emissions [1][2][3][4][5][6].Recent developments in information and communication technologies have made the implementation of road pricing schemes easier.There are some wellknown successful examples of electronic road pricing all over the world, including the congestion charging schemes in California, Singapore, and London [7,8].
There is a substantial body of literatures on road toll pricing, and the policy objectives of these works can be generally summarized as follows: minimizing total travel time of transportation networks with fixed travel demand [1-3, 9, 10], maximizing total social welfare of all travelers in the network with elastic travel demand [11][12][13][14], maximizing the network capacity [15], or maximizing the network travel time reliability [16], and so on.For recent reviews on road tolls, see [17,18].However, as the issues of society, environment, and energy caused by traffic congestion become worse and worse, improving the sustainability of the urban transportation systems is nowadays the urgent concern of transportation network management.In this paper, we attempt to investigate the optimal road toll design problem from the perspective of sustainable development of transportation system.In order to quantitatively measure the sustainability of transportation system with different network toll schemes, some quantitative indicators must be adopted.In general, the sustainable development level of transportation system can be described from the aspects of transportation network efficiency and environmental pollution, as well as the energy consumption.Consequently, in this paper we take the network-wide total travel time, total vehicular emission, and total fuel consumption as the measuring indicators of the sustainable development level of transportation system, and furthermore, in order to simultaneously consider the impacts of all these three indicators on the sustainability of transportation system, we integrate all these three sustainability indicators into a sustainable development index (SD-index) by a linear combination.And thus, the optimization objective in the road toll design problem is to maximize the SD-index of the urban transportation system.
The remaining sections are organized as follows.Section 2 presents the bilevel programming model of optimal road toll design problem from the perspective of sustainable development.Section 3 introduces the solution approach, and in this section we adopt a combined GA-GP algorithm to solve the bilevel optimization model.In order to verify the proposed model and algorithm, the Nguyen-Dupuis network is taken as a numerical example in Section 4 and conclusions are provided in Section 5.

Bilevel Programming Model
Since the sustainable development index (SD-index) is determined by the total weighted sum of the emitted vehicular emission and the fuel consumption, as well as the total travel time in the network, and, furthermore, these three sustainability indicators can be largely affected by the link flows, thus we must figure out the relationships between these sustainability indicators and link flows before getting into the essential investigation; in other words, we must figure out how these sustainability indicators change with respect to the link flow volumes in the network.In fact, there are many investigations having focused on deriving the desirable function relationships to date, and some analytical function expressions are obtained based on the realistic survey data [4,5,19,20].In this paper, we adopted the emission function used in TRANSYT-7F [4,5,19] and the fuel consumption function derived in [20], as well as the commonly used BPR function.
The optimal road toll design problem at network-wide level from the perspective of sustainable development with fixed demand static user equilibrium flow constraint is formulated as follows: where x() is the equilibrium flow pattern defined by the following fixed demand static user equilibrium problem with a certain road toll scheme imposed on the network: The notations are defined in Notations section.
In this bilevel model, the upper level describes leader or policy problem and the lower-level model represents follower or user's behavioral problem.In the optimal road toll design problem from the perspective of sustainable development, the upper-level problem is to determine an optimal decision for road toll scheme to make the sustainable development index (i.e., SD-index) maximum.The lower-level problem represents a user equilibrium assignment problem that describes travelers' path-choosing behavior under a certain road toll scheme.
The proposed bilevel model is generally difficult to solve because evaluation of the upper-level objective functions requires solving the lower-level subprogram.Besides, since this bilevel model is a NP-hard problem and the upper level is nonconvex and nondifferentiable in x which is defined by the lower level model, therefore, the existing exact solution algorithms are not available for solving this bilevel model.Consequently, in this paper we adopt a combined genetic algorithm and gradient projection algorithm (GA-GP) to solve the proposed model; to be specific, the traffic assignment problem with road toll in the lower level is solved by GP algorithm while the toll design problem in the upper level is solved by GA.

Solution Approach
Genetic algorithm (GA) is one of the most well-known search heuristics for solving optimization problems [21][22][23][24].GA uses techniques such as inheritance, mutation, selection, and crossover which are inspired by natural evolution and can search for optimal or near-optimal solutions for an optimization problem over the search domain.GA is implemented as a computer simulation in which a population of chromosomes of candidate solutions to an optimization problem evolves toward better solutions.Solutions can be represented in binary or real-coded.The evolution usually starts from a population of randomly generated individuals and happens in generations.In each generation, the fitness value of every individual in the population is evaluated; multiple individuals are randomly selected from the current population (based on their fitness value) and modified (recombined and possibly randomly mutated) to form a new population.The new population is then used in the next iteration of the algorithm.Commonly, the algorithm terminates when either a maximum number of generations have been produced, or a satisfactory fitness level has been reached for the population.
In this paper the chromosomes of the GA are designed in this way: all the tolled links on the network are successively numbered, and each gene in one chromosome represents the toll level charged on the corresponding link.For the chromosomes in the initial generation, all their genes are randomly generated between 0 and .And meanwhile, since the optimal toll design variables in the upper-level subprogram are real, we adopt the real-coded representation to represent the design variables   , ∀ ∈  with a length equal to the number of tolled links in the network.The value of each gene represents the link toll level, which is constrained by the upper and lower bound, that is,  and 0 in this paper.In the following, a brief description of the GA implementation procedure is provided.The implementation procedure of GA is shown in Figure 1.
Step 1 (Initialization).A feasible search space should be first defined, and all the chromosomes will be operated and evaluated in this constrained space.Set the maximum number of generations as , population size as , crossover probability as   , and mutation probability as   .Then set  = 0, and create the original population [Pop(0)] consisting of  chromosomes, which can be denoted by [(0)] 1 , [(0)] 2 , . . ., [(0)]  , . . ., [(0)]  .It is noteworthy that these chromosomes are randomly generated from the feasible region, the length of each chromosome equals the number of the tolled links in the network, and each gene of the chromosomes represents the corresponding link toll level which lies in the interval between 0 and .
Step 2. For each chromosome [()]  ( = 1, 2, . . ., ) in the current population [Pop(t)], solve the lower level program by Gradient Projection (GP) algorithm and evaluate the fitness value with the sustainable development index (SDindex), which is the objective function of the upper level program.
Step 3. If  > , then stop and output the optimal chromosome; else set  :=  + 1, and then go to next step (Step 4).
Step 4 (Reproduction).In this process the chromosomes reproduce according to their fitness values; that is, the bigger the SD-index value is, the higher the chance becomes that the chromosome will be chosen as the offspring chromosome.And meanwhile, a simpler tournament selection is adopted to decide whether or not a chromosome can reproduce depending on SD-index value.The procedure of tournament selection is that   ⋅  chromosomes with maximum SDindex values are added into the population and   ⋅  chromosomes with minimum SD-index values are discarded from the population.The resulting population has the same size with the original one.After the selection, all chromosomes are completely put in the mating pool.
Step 5 (Crossover).Since the chromosomes in the toll design problem are coded using a real-coded representation, the arithmetic crossover is used.This method randomly divides the  chromosomes into /2 pairs, and each pair yields two new chromosomes via a linear combination.Assume that both  and   are selected and  is a random number chosen from 0 to 1.If  >   , then the following crossover operations for y and y  are performed: where () and (  ) are, respectively, the fitness values of chromosomes  and   , while  and   are the resulting children chromosomes, and  ∈ [0, 1] is a random number determining the crossover grade of chromosomes  and   .If  <   , no crossover operation is performed.
Step 6 (Mutation).With a lower probability, randomly choose some genes from all the chromosomes in current population, and then modify value of these genes by a pseudorandom number between 0 and .In this paper the multiple mutation is adopted and only   ⋅  random chromosomes in the population [Pop(t)] are chosen to be mutated.The formula of mutation operation for a selected  is given by where  is a small positive constant and  is a random perturbation vector to produce small disturbances on .If the resulting chromosome is outside the feasible region, the original one is retained.This process may also generate some new chromosomes.
Step 7.After the genetic operations of reproduction, crossover, and mutation, a new offspring population will be generated.And then go back to Step 2.

Numerical Example
In this section, an example demonstrates the performance of the proposed model and algorithm for the optimal toll design problem.The network used in the numerical experiment is depicted in Figure 2, which was also used for numerical examples in [25].Initialize: t = 0 No, t = t + 1 Set the GA parameters (T, N, P c , and P m ) Figure 1: The implementation procedure of genetic algorithm.given in Table 1.The link travel time function is the standard BPR function: where   ,  0  , V  , and   are the travel time, free-flow travel time, flow, and capacity of link , respectively.All the link related parameters are provided in Table 2.In this paper, Matlab GA Toolbox is applied to solve the bilevel programming model, with the GA parameters set as follows:  = 100,  = 20,   = 0.65, and   = 0.01.
The optimal toll scheme and the equilibrium link flow pattern under UE, as well as the corresponding link travel times, are presented in Table 3. Besides, in order to compare the results of the toll scenario with the no-toll scenario and analyze the impacts of road toll on SD-index, we also provide the UE solution of no-toll scenario in Table 3.In addition, under UE with the optimal road toll scheme the used equilibrium paths between each OD pair and the corresponding path flows, as well as the resultant path costs (in fact, the path cost is generalized path travel cost which consists of the path travel time and the total toll charged on that path), are given in Table 4.And meanwhile, the used paths and the equilibrium path flows under UE without road toll are also presented in Table 4.
As shown in Table 3, when there is no tolled link in the network (i.e., no-toll scenario), it gives and with this optimal toll scheme imposed in the network, the SD-index is about 0.802, with an increasing ratio of 15.56% compared to that in the no-toll scenario.
It can be seen from Table 4 that when the network flow pattern reaches UE, the total travel costs on all the used routes between the same OD pairs are equal to and less than the costs on all the unused routes, which is entirely consistent with the Wardrop UE principle.Besides, from Table 4 we can also find that when the optimal toll charge scheme is imposed in the network, there are 3 equilibrium routes used between OD pair (1, 3) under UE, while there are 2 equilibrium routes used between the same OD pairs in the no-toll scenario.In other words, it seems that the network flows under UE can be more evenly distributed throughout the whole network, and thus the traffic congestion in the road can be mitigated if the road toll scheme is properly designed.In fact, the V/C ratios in links 5, 9, 10, and 11 are all greater than 1 under UE in the no-toll scenario, while the V/C ratios are all smaller than 1 under UE with the optimal road toll scheme imposed on the network, which shows that the appropriate toll charge scheme can not only increase the SD-index of transportation system, but also alleviate the road congestion simultaneously.

Conclusions
In this paper, the optimal road toll design problem (TDP) from the perspective of sustainable development at networkwide level has been studied.The TDP is formulated as a bilevel programming model, and, specifically, the upper level model is to determine the optimal road toll scheme with the aim of maximizing sustainability level of transportation system and the lower level model is to depict the users' travel behaviors under a certain network toll scheme.Moreover, this paper integrates three measurable sustainability indicators of transportation system, that is, total travel time, total fuel consumption, and total vehicular emissions into a sustainable development index (SD-index), which is the final optimization objective of the bilevel programming model.Besides, the combined genetic algorithm and gradient projection algorithm (GA-GP) is adopted to solve the bilevel model, in which the gradient projection (GP) algorithm is mainly used to solve the lower level traffic assignment problem with a certain network toll scheme.In order to verify the proposed model and algorithm, we take the Nguyen-Dupuis network for the numerical example, and the computing results show that the model and algorithm are feasible and efficient and can be taken as a reference for the related research.Further investigations will concentrate on the toll design problem when there are multiple user groups, as well as the impacts of demand uncertainty on the optimal toll scheme.

A:
The set of links in the network R: The set of origins S: The set of destinations K  : The set of paths between OD pair (, ), ∀ ∈ ,  ∈  f: The vector of path flows,  = [   ], ∀ ∈ ,  ∈  x: The vector of link flows,  = [  ], ∀ ∈  : Thevectorofroadtolllevel, = [  ], ∀ ∈  l: Thevectoroflinklength, = [  ], ∀ ∈  U: The vector of upper bound of road toll charged on links  = [  ], ∀ ∈  : Digit number of the function value (, ), and 10 −1 ⩽ (, ) ⩽ 10  t: The vector of link travel times,  = [  (  ,   )], ∀ ∈    : The emission per vehicle on link    : 100 km fuel consumption per vehicle on link    , : The link-path incidence,   , = 1 if link  is on the th route connecting origin  and destination , and   , = 0 otherwise : Value of time, which is assumed to be 0.5 for all travelers in this paper , , : The monetary conversion factors.

Figure 2 :
Figure 2: Nguyen-Dupuis network used in the example.

Table 2 :
Link related parameters of the example network (min, veh/hr, km).