Research Article  Open Access
Optimal Road Toll Design from the Perspective of Sustainable Development
Abstract
This paper investigates the optimal road toll design problem from the perspective of sustainable development at networkwide 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 (SDindex) 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 SDindex 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 (GAGP) 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 NguyenDupuis network for the numerical example, and the computing results show that the model and algorithm are effective and efficient.
1. 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, lowcarbon, and energysaving 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–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–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 networkwide 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 (SDindex) by a linear combination. And thus, the optimization objective in the road toll design problem is to maximize the SDindex 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 GAGP algorithm to solve the bilevel optimization model. In order to verify the proposed model and algorithm, the NguyenDupuis network is taken as a numerical example in Section 4 and conclusions are provided in Section 5.
2. Bilevel Programming Model
Since the sustainable development index (SDindex) 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 TRANSYT7F [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 networkwide level from the perspective of sustainable development with fixed demand static user equilibrium flow constraint is formulated as follows: where 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 lowerlevel model represents follower or user’s behavioral problem. In the optimal road toll design problem from the perspective of sustainable development, the upperlevel problem is to determine an optimal decision for road toll scheme to make the sustainable development index (i.e., SDindex) maximum. The lowerlevel problem represents a user equilibrium assignment problem that describes travelers’ pathchoosing behavior under a certain road toll scheme.
The proposed bilevel model is generally difficult to solve because evaluation of the upperlevel objective functions requires solving the lowerlevel subprogram. Besides, since this bilevel model is a NPhard problem and the upper level is nonconvex and nondifferentiable in 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 (GAGP) 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.
3. Solution Approach
Genetic algorithm (GA) is one of the most wellknown search heuristics for solving optimization problems [21–24]. GA uses techniques such as inheritance, mutation, selection, and crossover which are inspired by natural evolution and can search for optimal or nearoptimal 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 realcoded. 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 upperlevel subprogram are real, we adopt the realcoded 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 , and create the original population [Pop(0)] consisting of chromosomes, which can be denoted by . 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 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 , 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 SDindex 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 SDindex value. The procedure of tournament selection is that chromosomes with maximum SDindex values are added into the population and chromosomes with minimum SDindex 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 realcoded representation, the arithmetic crossover is used. This method randomly divides the chromosomes into 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 and are performed: where and are, respectively, the fitness values of chromosomes and , while and are the resulting children chromosomes, and 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.
4. 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]. The road network consists of 13 nodes, 19 links, and 25 paths. There are 4 OD pairs; that is, , , , and . The travel demands for each OD pair are given in Table 1. The link travel time function is the standard BPR function: where , , , and are the travel time, freeflow 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: , , , and .
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 notoll scenario and analyze the impacts of road toll on SDindex, we also provide the UE solution of notoll 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., notoll scenario), it gives the result of deterministic user equilibrium solution without implementing road pricing scheme. In this situation, the SDindex is about 0.694, which indicates that the sustainable development level in the notoll scenario is relatively not very high and may increase to a higher level by imposing an appropriate road toll scheme. In this example the optimal road toll scheme is = 1.31 0.29 0.04 0.52 0.70 0.15 0.34 0.02 0.67 2.10 1.95 0.20 1.23 0.42 0.18 0.66 0.00 0.45 0.81, and the corresponding UE link flow is = 800 400 245 555 752 293 621 131 245 376 645 355 493 486 355 507 0 400 493, and with this optimal toll scheme imposed in the network, the SDindex is about 0.802, with an increasing ratio of 15.56% compared to that in the notoll 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 notoll 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 notoll 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 SDindex of transportation system, but also alleviate the road congestion simultaneously.
5. 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 (SDindex), which is the final optimization objective of the bilevel programming model. Besides, the combined genetic algorithm and gradient projection algorithm (GAGP) 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 NguyenDupuis 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.
Notations
:  The set of links in the network 
:  The set of origins 
:  The set of destinations 
:  The set of paths between OD pair 
:  The vector of path flows, 
:  The vector of link flows, 
:  The vector of road toll level, 
:  The vector of link length, 
:  The vector of upper bound of road toll charged on links 
:  Digit number of the function value , and 
:  The vector of link travel times, 
:  The emission per vehicle on link 
:  100 km fuel consumption per vehicle on link 
:  The linkpath incidence, if link is on the th route connecting origin and destination , and otherwise 
:  Value of time, which is assumed to be 0.5 for all travelers in this paper 
The monetary conversion factors. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work described in this paper was supported by Grants from the National Natural Science Foundation of China (nos. 51178110, 51378119), the Graduate Innovation Project of Jiangsu Province (No. KYLX 0178), the Ph.D. Programs Foundation of Ministry of Education of China (20120092110062) and the Six Talent Peaks Program of Jiangsu Province of China (JZ003). The authors are also grateful to the anonymous referees for providing useful comments and suggestions to improve the quality and clarity of the paper.
References
 H. Yan and W. H. K. Lam, “Optimal road tolls under conditions of queueing and congestion,” Transportation Research Part A: Policy and Practice, vol. 30, no. 5, pp. 319–332, 1996. View at: Publisher Site  Google Scholar
 H. Yang and H.J. Huang, “Principle of marginalcost pricing: How does it work in a general road network?” Transportation Research A, vol. 32, no. 1, pp. 45–54, 1998. View at: Publisher Site  Google Scholar
 H. Yang, “System optimum, stochastic user equilibrium, and optimal link tolls,” Transportation Science, vol. 33, no. 4, pp. 354–360, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Yin and S. Lawphongpanich, “Internalizing emission externality on road networks,” Transportation Research D, vol. 11, no. 4, pp. 292–301, 2006. View at: Publisher Site  Google Scholar
 L. Chen and H. Yang, “Managing congestion and emissions in road networks with tolls and rebates,” Transportation Research B, vol. 46, no. 8, pp. 933–948, 2012. View at: Publisher Site  Google Scholar
 Z.C. Li, W. H. K. Lam, S. C. Wong, and A. Sumalee, “Environmentally sustainable toll design for congested road networks with uncertain demand,” International Journal of Sustainable Transportation, vol. 6, no. 3, pp. 127–155, 2012. View at: Publisher Site  Google Scholar
 G. Santos, Road Pricing: Theory and Evidence, Elsevier, Amsterdam, The Netherlands, 2004.
 T. D. Hau, “Congestion charging mechanisms for roads, part II—case studies,” Transportmetrica, vol. 2, no. 2, pp. 117–152, 2006. View at: Google Scholar
 X. Ban and H. X. Liu, “A linknode discretetime dynamic second best toll pricing model with a relaxation solution algorithm,” Networks and Spatial Economics, vol. 9, no. 2, pp. 243–267, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Lou, Y. Yin, and S. Lawphongpanich, “Robust congestion pricing under boundedly rational user equilibrium,” Transportation Research Part B: Methodological, vol. 44, no. 1, pp. 15–28, 2010. View at: Publisher Site  Google Scholar
 H. Yang and Q. Meng, “Departure time, route choice and congestion toll in a queuing network with elastic demand,” Transportation Research B, vol. 32, no. 4, pp. 247–260, 1998. View at: Publisher Site  Google Scholar
 E. T. Verhoef, “Secondbest congestion pricing in general networks: heuristic algorithms for finding secondbest optimal toll levels and toll points,” Transportation Research B: Methodological, vol. 36, no. 8, pp. 707–729, 2002. View at: Publisher Site  Google Scholar
 X. Zhang and H. Yang, “The optimal cordonbased network congestion pricing problem,” Transportation Research B, vol. 38, no. 6, pp. 517–537, 2004. View at: Publisher Site  Google Scholar
 Z. Liu, Q. Meng, and S. Wang, “Speedbased toll design for cordonbased congestion pricing scheme,” Transportation Research Part C: Emerging Technologies, vol. 31, pp. 83–98, 2013. View at: Publisher Site  Google Scholar
 X. Zhang and B. van Wee, “Enhancing transportation network capacity by congestion pricing with simultaneous toll location and toll level optimization,” Engineering Optimization, vol. 44, no. 4, pp. 477–488, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 H. Li, M. C. J. Bliemer, and P. H. L. Bovy, “Network reliabilitybased optimal toll design,” Journal of Advanced Transportation, vol. 42, no. 3, pp. 311–332, 2008. View at: Publisher Site  Google Scholar
 H. Yang and H. J. Huang, Mathematical and Economic Theory of Road Pricing, Elsevier, Amsterdam, The Netherlands, 2005.
 A. de Palma and R. Lindsey, “Traffic congestion pricing methodologies and technologies,” Transportation Research C, vol. 19, no. 6, pp. 1377–1399, 2011. View at: Publisher Site  Google Scholar
 C. E. Wallace, K. G. Courage, M. A. Hadi, and A. C. Gan, TRANSYT7F User’s Guide, University of Florida, Gainesville, Fla, USA, 1998.
 L. Xu and G.Z. Cheng, “Setting minimum vehicle speed limit on freeway based on speed scattering and economic speed,” Journal of Jilin University, vol. 40, no. 3, pp. 661–665, 2010. View at: Google Scholar
 D. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning, AddisonWesley, Reading, Mass, USA, 1989.
 T. Gen and R. Cheng, Genetic Algorithms and Engineering Design, New York, NY, USA, John Wiley & Sons, 1997.
 X.D. Xu and L. Cheng, “Bilevel programming model and hybrid solution algorithm for the configuration of oneway streets,” System Engineering Theory and Practice, vol. 29, no. 10, pp. 180–187, 2009. View at: Google Scholar
 T. Xu, H. Wei, and G. Hu, “Study on continuous network design problem using simulated annealing and genetic algorithm,” Expert Systems with Applications, vol. 36, no. 2, pp. 1322–1328, 2009. View at: Publisher Site  Google Scholar
 S. Nguyen and C. Dupuis, “An efficient method for computing traffic equilibria in networks with asymmetric transportation costs,” Transportation Science, vol. 18, no. 2, pp. 185–202, 1984. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Lin Cheng and Fei Han. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.