Research Article  Open Access
Zhengfeng Huang, "Proactive Traffic Information Control in Emergency Evacuation Network", Discrete Dynamics in Nature and Society, vol. 2015, Article ID 958052, 10 pages, 2015. https://doi.org/10.1155/2015/958052
Proactive Traffic Information Control in Emergency Evacuation Network
Abstract
Traffic demand in emergency evacuation is usually too large to be effectively managed with reactive traffic information control methods. These methods adapt to the road traffic passively by publishing realtime information without consideration of the routing behavior feedback produced by evacuees. Other remedy measures have to be prepared in case of nonrecurring congestion under these methods. To use the network capacity fully to mitigate nearfuture evacuation traffic congestion, we propose proactive traffic information control (PTIC) model. Based on the mechanism between information and routing behavior feedback, this model can change the route choice of evacuees in advance by dissipating strategic traffic information. Generally, the nearfuture traffic condition is difficult to accurately predict because it is uncertain in evacuation. Assume that the value of traffic information obeys certain distribution within a range, and then realtime traffic information may reflect the mostlikely nearfuture traffic condition. Unlike the realtime information, the proactive traffic information is a selection within the range to achieve a desired level of the road network performance index (total system travel time). In the aspect of the solution algorithm, differential equilibrium decomposed optimization (DEDO) is proposed to compare with other heuristic methods. A field study on a road network around a large stadium is used to validate the PTIC.
1. Introduction
For the natural or manmade disasters, the surrounding road network should be fully utilized to evacuate the affected people. Otherwise, stochastic evacuation flows may induce nonrecurring traffic congestion. Therefore, optimizing the operation of evacuation flows to improve emergency evacuation efficiency is one of the main goals of the disaster relief. From the perspective of behavior, how to adjust evacuees’ routechoice result to mitigate congestion in advance is a critical point. In other words, we should consider evacuees’ routing behavior feedback in the evacuation management.
Compared to general transportation network management, routing behavior feedback is less considered in the evacuation condition. Unlike the widely applied cases in the former condition, such as street capacity optimization [1] and signal optimization [2], only few studies are found in the latter condition, that is, Xie and Turnquist [3] in the planning of lane reversals and intersection crossings mitigation, Chiu and Mirchandani [4], and Paz and Peeta [5] in the dynamic choice of informationpublishing routes.
Traffic information control can be a feasible evacuation management. However, decision makers usually passively published realtime information to evacuees whose routechoice result is not being optimized. Thus, the performance index of evacuation network can be improved only to a limited extent. Generally, the realtime information can solely reflect the nearfuture traffic condition with the maximum probability because this condition is uncertain in evacuation. At the current time interval, the value of nearfuture traffic information may emerge to be a probabilistic distribution within a range. Selecting any of them may lead to varied routechoice results of evacuees; thus, we need to develop a method to choose the appropriate value instead of the realtime value as the information provision to improve evacuation efficiency.
To design traffic information considering behavior feedback, the mapping between traffic information and evacuation flows should be addressed. Under specific information, the mapping is generated by two intermediate processes, the route choice and traffic loading. This is also the routing behavior feedback of evacuees, through which the decision maker can design the information to optimize the system performance. We define such optimization as proactive traffic information control (PTIC). Thus, this type of information optimization is akin to network design problem, which belongs to Stackelberg game. Because of the nonclosedform of trafficloading simulation steps involved in this application, the general operations research methods used for nonlinear planning cannot be applied. The useful methods are heuristic, such as equilibrium decomposed optimization algorithm [6], linear approximation algorithm [7], and intelligent algorithm [8]. The last algorithm can have the potency of global convergence because the searching process is akin to enumeration. However, it is imperfect for its low efficiency. The first method is a numerical method which uses onedimensional search to reduce the solution range, but the global solution may be ignored at the initial iterations with large step size and also we cannot ensure the convergence of the solutions in the whole operation. The main idea of the second method is to use Taylor expansion to change the objective to a linear function to solve drop direction and golden selection method to solve step size, such as FrankWolfe (FW) algorithm. The only note is that the derivative or difference of objective with respect to decision variable should be provided in the generated subproblem. However, this method may converge to a local optimum near the initial solution because the drop direction and step size are too accurate to jump out of the local trap.
To better mitigate nearfuture congestion for emergency evacuation network, we propose PTIC to provide the traffic information strategically according to routing behavior feedback. The minimal total system travel time (TSTT) is set as the objective, which is composed of total evacuation time, arrival time penalty, and congestion cost. A corresponding differential equilibrium decomposition optimization (DEDO) algorithm is proposed and demonstrated to be advantageous in the accuracy.
In what follows, the PTIC model and solving steps of DEDO are introduced. Then, the theoretical method is applied in the field case of evacuation network around Nanjing Olympic Stadium in China.
2. Method
2.1. PTIC Model
The following constraints (1–8) and the minimal TSTT objective (9) are designed to formulate PTIC model. The constraints are in accord with oneshot stochastic dynamic traffic assignment, composed of stochastic route choice and oneshot traffic loading. The reason for adopting a nonequilibrium model is that no historic travel experience can help evacuees make a game decision; the reason for adopting stochastic model is that evacuees cannot know the a priori flow condition in the nonrecurring event.
2.1.1. RouteChoice Utility
There are some exogenous routerelated variables influencing evacuees on the route choice. The determination of exogenous variables mainly depends on two criteria (evacuation time and information reliability). As for evacuation time, it is no doubt that people prefer to choose shorter routes. Besides, Wardman et al. [9] concluded, using a stated preference (SP) survey, that delay time is more useful than travel time in the route choice. With respect to information reliability, risk aversion will emerge in travelers’ mind when unreliable information is provided. Hua et al. [10] showed that the disseminated delaytime reliability is a significant factor influencing route choice. Therefore, exogenous variables of route length, delay time, and disseminated delaytime reliability are chosen as the routechoice explanatory variables. For easy recognition, these three explanatory variables corresponding to route at an informationupdating interval are represented as (m), (min), and (%) in sequence. The utility value of route at the informationupdating interval can be calculated by
In formula (1), the evacuees are assumed to be homogeneous, so use constant to represent the utility value of endogenous attributes of evacuees. We can lose this term “commensurable and can be removed” in the discrete choice model. To consider the overlap line problem, the Clogit model [11], a type of multinomial logit model, is usually adopted as the discrete choice model because of its simple form and easy calibration. The variable is the commonality factor in the utility of Clogit model. The variable is the proportional weight (often represented as a length ratio) of the overlap section between path and another route of the same OD pair. The set includes all the overlap sections related to route . The variable is the number of routes, connecting the same OD pair with , which share the overlap section . The explanatoryvariable coefficients , , , and need to be calibrated using SP survey.
2.1.2. Analytical RouteChoice Probability Function
According to the Clogit model, the analytical routechoice probability for choosing route from origin cell to destination at an informationupdating interval can be determined by where is the route set from cell to destination. Because the destination of evacuation network is normally set as a virtual node, the index for destination can be eliminated to reduce notation complexity. The destination is known implicitly through the route index.
2.1.3. OneShot Traffic Loading
Daganzo’s [12] cell transmission model (CTM) can load traffic dynamically and successfully embedded in simulation software VISTA. But this basic cell cannot extract congestion information. In the following, we propose an enhanced CTM to overcome this shortcoming. The following variables at each (loading) interval are defined a priori:: set of cells succeeding cell , : set of cells preceding cell , : maximal number of vehicles that can flow into or out of cell at interval , : maximal number of vehicles that can be present inside a cell at interval , : demand from source cell at interval , : vehicle occupancy of cell at interval , : flow outflowing from cell at interval , : flow from cell to cell at interval .It is noted that the unit time interval for traffic loading is shorter than the informationupdating time interval.
Assuming the timevarying demands and routechoice ratios are provided, we can turn to the procedure of the enhanced CTM (shown in Figure 1) to simulate the flow. Regarding the proposed varied length cell, each link in our model is divided into three types of cells composed of 1 tail cell, 1 head cell, and a few approach cells (its number is based on the number of stream directions). Head cell and tail cell are used to capture the intersection queue length and traffic spillback separately. The latter two indexes could serve the representation of congestion cost in the PTIC objective. For example, we can let the length of approach cell be equal to unit time first and set the length of head cell to be equal to a threshold value required by decision maker to penalize the intersection queue and give the remaining length to tail cell to check the queue spillback for penalization.
Assuming backward propagation owns freeflow speed, the following procedure can be used to simulate vehicle occupancy of cell for interval , based on previous vehicle occupancies ( is set to be 0).
(1) Formulations for Sending and Receiving Flows. The formulae for receiving and sending flows of cell can be described in where is the needed number of time intervals used to pass through cell with the freeflow speed. It depends on the type that cell belongs to.
(2) Transmission Flow Formulation. We assume cell is a tail cell here. The transmission flow formula for flow from tail cell to head cell is expressed in
The transmission flow formula for flow from head cell to approach cell is denoted by formula (5). is the proportions of vehicle occupancy that can go to cell . It is an implicit function with respect to routechoice probability . If equals zero then the corresponding term inside the brace of formula (5) should be deleted:
The transmission flow formula for flow from approach cell to tail cell is denoted by where equals which represents the fraction of sending flow in all the sending flows at preceding cells of cell at interval .
During each time interval, both transmission flows in formulae (5) and (6) should be rounded to ensure that the flows and vehicle occupancies are integers.
(3) Vehicle Occupancy Formulation. The following formula is used to update the vehicle occupancy of cell for the next interval:
2.1.4. Objective Formulation
The objective is designed in formula (8) and is composed of three terms. The first term is the total evacuation time. In the enhanced CTM, the total evacuation time is equal to total vehicle occupancies during the evacuation when veh time interval is the unit. The second term is the arrival time penalty, which means each consumed time interval should be penalized once the vehicle departs the origin. Why is penalty added to the cost of vehicles according to its arriving time at destination? One reason is to quicken the step of emergency. Another advantage is that the earlier leaving time will be chosen in advance when two routes with the same cost of different departure time need to be compared. The third term is the congestion cost, which penalizes the consumed time in the long queue and spillback: where is the set of all the cells, is the set of all the head and tail cells, is the set of origin cells, and is the dummy destination cell. is the set of all time intervals and and are the penalties for the arrival time and oversaturated vehicles separately. The traffic state judgment of cell at time interval can roughly be the following formula: It is obvious that if long queue or queue spillback occurs at cell , corresponding penalty is added to decrease the dynamic routes and their vehicles passing through this cell at this time.
The vector in formula (8) represents the decision variable vector composed of information variables shown in formula (1). In the proactive traffic information model, can affect the objective indirectly. Its influence process is shown as follows: information variables analytical routechoice probability simulated vehicle occupancy TSTT. To guarantee the mapping between the information and the evacuation flows is unique, we set single decision variable at each origin point per informationupdating interval. The problem about how to determine decision variables from the selected information should be solved according to the specific application.
2.2. Solving Steps of DEDO
The logical framework of solving PTIC in emergency evacuation is shown in Figure 2. The inputs are the physical attributes of the evacuation network, timevarying evacuation demand, the explanatoryvariable coefficients of routechoice utility, and the timevarying information boundary values. The output is the published strategic information. Two main modules (oneshot stochastic dynamic traffic assignment and DEDO procedure) shown in Figure 2 are elaborated in the following subsections.
Oneshot stochastic dynamic traffic assignment module is actually within the module of DEDO algorithm. It is composed of stochastic route choice and dynamic traffic loading. The updated information of DEDO should be the preparedness of utility, which can immediately be the input of this module. At each informationupdating interval of the internal process, Clogit function is applied to obtain the analytical routechoice ratio; and then a oneshot CTM simulation is performed. Step by step, traffic loadings at the following informationupdating intervals are calculated. This module provides the data for the corresponding objective calculation of the DEDO algorithm.
The DEDO is proposed by combining the EDO and FW algorithms. The EDO is set as the outside layer algorithm. For every certain number (equal to the decision variable number) of iterations, we would choose the current minimal boundary solution of EDO as the initial input of FW. We could solve its surrounding local optimum by the inside layer FW algorithm then. When the algorithm of EDO is over, we compare all these local solutions to obtain the minimal result as the optimum. The EDO can go through diverse feasible solutions as they are narrowed gradually from the initial maximal boundary range. The effectiveness of FW is related to initial input. The global optimum is more likely to be found using FW with diverse initial inputs than uniform ones. Therefore, the advantages of both EDO and FW are utilized to obtain a more accurate result which is expected to avoid the local optima. In the drop direction searching process of FW, difference instead of deviation is prepared because the nonclosedform of the problem (due to the embedding of the simulation process) prohibits the solving of the deviation. Therefore, we name this method as differential decomposed equilibrium optimization method.
The steps of DEDO are as follows.
Step 1 (initialization). Select vectors and representing, respectively, the initial lower and upper bounds on information variable vector, where is the origin cell and is the informationupdating time interval. Simulate the vehicle occupancy vectors and . Evaluate the objective value for each vector, and . Set .
Step 2 (FW algorithm).
Step 2.1. Set . Let be the feasible solution.
Step 2.2. Calculate the difference of (and ) with respect to through flow simulation; obtain the form (and ). Find auxiliary feasible solution such that the direction from to provides the maximum feasible drop in the following linearization, which has already deleted the constant terms in the initial objective function:
where is the set of informationupdating intervals and and are the elements in and .
Step 2.3. Use golden selection method to find the associated optimal step size such that
and update the current solution by .
Step 2.4. If the gap between two successive solutions is less than , then set and turn to Step 3 to continue the algorithm of EDO; otherwise, turn back to Step 2.2.
Step 3 (golden selection). For each and , the search procedure specifies the following two points in the interval : Let other elements in (and ) be assigned to other elements in (and ). Perform flow simulation to obtain the vehicle occupancy vectors and . Evaluate the objective value for each vector, and . If , then let ; otherwise, let .
Step 4 (stopping criterion). If , output optimized objective and the corresponding solution y; otherwise, turn back to Step 2.
3. Method Application
3.1. Input
Input the evacuation network, evacuation demand, analytical routechoice probability function, and initialized information in sequence.
3.1.1. Evacuation Network and Demand
The Nanjing Olympic Center (Figure 3(a)) in China, with a capacity of 87,000 people, is assumed to be under the threat of a terrorist incident. People in the stadium should evacuate the stadium and nearby parking lots. They are not safe until they have escaped the described road network that surrounds the Olympic Stadium. The roads within the network constitute a closed area that allows only a northeast bound flow. Parking lots, located on the north side, south side, and east side (two lots on the east side) of the Olympic Center, are the origins for the evacuation. The VMSs are set up at origins and . The cell structure containing a dummy destination (with unlimited maximal vehicle occupancy and flow capacity) is shown in Figure 3(b). Its corresponding evacuation traffic demand and cell characteristics are shown in Table 1.

(a)
(b)
3.1.2. Analytical RouteChoice Probability Function and Information Initialization
Considering the limitation of evacuees’ information analysis capability, at most three routes are assumed available for travelers to choose from at each origin cell. These chosen routes () are spaceconstrained candidate routes in this study. The comparison factors are route length, number of intersections, number of lanes, and number of left turns, in that order. The routes with lower factor values are more likely to be chosen for the set.
The candidate route sets of the origin cells are listed in Table 2. The third column shows the resolved commonality factor.

We implement routechoice SP survey under aforementioned evacuationscenario assumption. Using the survey data, specify the utility coefficients by maximum likelihood method. The probability of a route chosen by evacuees at interval is formulated as follows:
Only one information factor is allowed to be selected as the decision variable at each origin per informationupdating interval. The route feature and variable range are the criteria in this case.(i)The route of shortest length from each origin cell is chosen as the route for decision variable setting, because more variations occurring in the shortest length route can create a broader range as the candidate of the strategic information variable.(ii)When the coefficients are close, the information factors with broader ranges are chosen as the decision variable because they allow adjusting routechoice probability greatly.
Once the decision variable is determined, other information variables could be assigned with realtime values or the other corresponding value depending on the relation among the information values.
In this case, for each origin cell, the delay time of shortestlength route is set as the decision variable according to the criterion description. The delays in other routes are realtime values obtained from the flow simulation, which could occur most probably in the nearfuture traffic conditions. The disseminated delaytime reliabilities are determined based on the nearfuture delaytime distribution. Each probability of the disseminated delay time is assigned to its reliability in the application. To be convenient, we assume the delay time obeys the normal distribution with the mean of the realtime delay and standard deviation of . The realtime delay is undoubtedly the most probable value, which is in accord with the previous assumption. The disseminated delaytime reliability can be represented as follows: where the function is the probability density function of standard normal distribution. is the divided number of delaytime range, which is set as 10 in the following application. It is used to ensure that the reliability is roughly equal to probabilistic density multiplying unit range of the divided delay time. Besides, is set to be 1; the range of strategic delay time is set to be mins.
3.2. Test
The advantage of the application of the DEDO in creating an optimal traffic information strategy is verified. This strategy outperforming reactive information control is indicated. In addition, the influence of informationupdating frequency and demand level are analyzed.
3.2.1. Effectiveness of the DEDO
is set as the end criteria of DEDO algorithm. The corresponding solution convergence of the DEDO is shown in Figure 4. When EDO algorithm completes an updating of all the variables per loop, we use the current minimal boundary solution of the acquired two boundary solutions as the initial point to obtain the corresponding FW result. The minimal boundary solution convergence of EDO showed that the single EDO method may not achieve the minimal TSTT even compared with the historical iterations. However, the FW can contribute to obtaining lower TSTT when taking the intermediate EDO minimal boundary solution as the initial solution. Varied initial solution can achieve varied result using FW; thus, the EDO combined with FW can complement each other. The minimal TSTT selected from all the FW solutions is 122,500 veh time intervals, with the variable solution listed in Table 3.

We use the genetic algorithm (GA) as another heuristic algorithm to compare with DEDO. The GA is configured to have an initial population of 40, a crossover probability of 0.9, and mutation probability of 0.1. We set the maximal allowed tolerance generation number of unchanged fitness to be 5, which is the end criterion.
Seeing from Figure 4, the DEDO achieves a better solution compared with the GA. The primary advantage of DEDO method over the GA method is that the accuracy result is obtained. Although the GA converges with few generations, the calculation time is not saved significantly with nearly 40 times of flow simulation in each generation. To avoid the premature convergence, some improved GA (such as micro GA) may be available in finding the optimal solution, but the calculation time will also increase. Therefore, we use the DEDO algorithm to analyze the following problems.
3.2.2. Advantage of the Proactive Traffic Information Compared with Reactive Traffic Information
If the realtime delays are set as the information, the ISTT would be 138,212 veh time intervals, which is higher than that using the proactive values. The realtime delay information has not considered the routechoice feedback influenced by congestion. However, the PTIC can mitigate congestion effectively and the eventual evacuation network will be more safe and reliable compared to the result of reaction traffic information control.
3.2.3. Correlation of InformationUpdating Frequency and TSTT
The needed time interval per informationupdating gap at each origin cell is changed from 10 to 50; then, we calculate TSTT for each informationupdating frequency combination during the change process. The results are listed in Figure 5. Although few unusual situations exist where the increase of informationupdating frequency creates higher TSTT, the majority of TSTT trends can roughly show that the higher frequency has more opportunity to create a lower TSTT. This implies that an increase in frequency has a positive effect on evacuation efficiency improvement.
Before the recommendation of informationupdating frequency, we classify the frequency group into three types according to the closeness of TSTT. High frequency is normally recommended; however, some inflection points of frequency may be better choices. There are two reasons why overhigh frequency is not a better choice. First, it is lack of room for evacuation efficiency improvement in high frequency range as Figure 5 shows: the TSTT improves rapidly at the low frequency range (group 3), while it improves more slowly at a high frequency range (group 1). Second, it is subject to the constraints of information collection technology, which does not allow a short updating time. Therefore, the frequency range in group 2 may be a recommendation for decision makers.
3.2.4. Correlation of Demand Matrix Multiplication Factor and Average Cost per Evacuee
As shown in Figure 6, when the demand matrix multiplication factor is below 1, average cost per evacuee increases slowly. When the factor exceeds 1, average cost per evacuee increases dramatically with demand matrix multiplication factor and the rate of change is almost the same. This indicates that network operation using PTIC may be not effective in the high evacuation demand range. We can set the demand assumed in this case as the boundary in which the PTIC works. When the average cost per evacuee increases sharply, we should also depend on other methods such as the design of contraflow lanes [13, 14] and evacuation schedule [15] to control the evacuation congestion.
4. Conclusion and Discussion
Considering the uncertainty of the evacuation flows, a PTIC method for network evacuation is constructed, wherein the information value is selected from a range instead of adopting realtime value directly. With the preparedness of routechoice utility investigation, we adopt the oneshot stochastic dynamic traffic assignment as the constraint and set the minimal total system travel time as the objective to construct the framework model, wherein congestion cost is embedded to mitigate oversaturated traffic. DEDO is proposed to be the algorithm to solve the proactive traffic information. Compared to other network design methods, DEDO can achieve a more accurate result. The result shows that PTIC can obtain lower objective value than reactive traffic information control.
We extend the study to the influence of informationupdating frequency and demand level. An increase in informationupdating frequency is welcome. However, the frequency of inflection point is advised to be chosen as the application value instead of high technologycost frequency. With respect to the demand influence, average cost per evacuee increases rapidly in higher demand range. Thus, the application of PTIC is effective under a specific demand level.
Three points should be noticed in the future research. (1) The evacuation of excessive demands is beyond the application scope of this method. How to combine PTIC and other traffic management measures to deal with excessive demands is a significant study direction. (2) The candidate route set of the case is generated a priori, because dynamically updating the routes is a challenge. The technology of column generation applied in static traffic assignment or numberlimited dynamic routes proposed by Chiu and Mirchandani [4] may be useful. (3) Before making this method a promising solution to meet the emergency evacuation need, the distribution function of the fluctuated delay time should be investigated with field delay data.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by the National Natural Science Foundation of China (no. 51408321). The author appreciates the Wisconsin Traffic Operations and Safety Laboratory. The author appreciates Professor Gang Ren, Professor Bin Ran, Professor Haixu Liu, Dr. Yang Cheng, and Dr. Peter J. Jin for their professional advices.
References
 G. A. Davis, “Exact local solution of the continuous network design problem via stochastic user equilibrium assignment,” Transportation Research Part B, vol. 28, no. 1, pp. 61–75, 1994. View at: Publisher Site  Google Scholar
 E. Cascetta, M. Gallo, and B. Montella, “Models and algorithms for the optimization of signal settings on urban networks with stochastic assignment models,” Annals of Operations Research, vol. 144, pp. 301–328, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Xie and M. A. Turnquist, “Integrated evacuation network optimization and emergency vehicle assignment,” Transportation Research Record: Journal of the Transportation Research Board, no. 2091, pp. 79–90, 2009. View at: Publisher Site  Google Scholar
 Y.C. Chiu and P. B. Mirchandani, “Online behaviorrobust feedback information routing strategy for mass evacuation,” IEEE Transactions on Intelligent Transportation Systems, vol. 9, no. 2, pp. 264–274, 2008. View at: Publisher Site  Google Scholar
 A. Paz and S. Peeta, “Behaviorconsistent realtime traffic routing under information provision,” Transportation Research Part C: Emerging Technologies, vol. 17, no. 6, pp. 642–661, 2009. View at: Publisher Site  Google Scholar
 C. Suwansirikul, T. L. Friesz, and R. L. Tobin, “Equilibrium decomposed optimization: a heuristic for the continuous equilibrium network design problem,” Transportation Science, vol. 21, no. 4, pp. 254–263, 1987. View at: Publisher Site  Google Scholar
 A. H. F. Chow and H. K. Lo, “Sensitivity analysis of signal control with physical queuing: delay derivatives and an application,” Transportation Research Part B: Methodological, vol. 41, no. 4, pp. 462–477, 2007. View at: Publisher Site  Google Scholar
 Y. Yin, “Geneticalgorithmsbased approach for bilevel programming models,” Journal of Transportation Engineering, vol. 126, no. 2, pp. 115–119, 2000. View at: Publisher Site  Google Scholar
 M. Wardman, P. W. Bonsall, and J. D. Shires, “Driver response to variable message signs: a stated preference investigation,” Transportation Research Part C: Emerging Technologies, vol. 5, no. 6, pp. 389–405, 1997. View at: Publisher Site  Google Scholar
 H. Hua, Y.F. Gao, and X.G. Yang, “Probabilistic traffic forecast method based on comprehensive transport information platform,” Journal of Traffic and Transportation Engineering, vol. 9, no. 3, pp. 122–126, 2009. View at: Google Scholar
 E. Cascetta, A. Nuzzolo, F. Russo, and A. Vitetta, “A modified logit route choice model overcoming path overlapping problems: specification and some calibration results for interurban networks,” in Proceedings of the 13th International Symposium on Transportation and Traffic Theory, Lyon, France, 1996. View at: Google Scholar
 C. F. Daganzo, “The cell transmission model: a dynamic representation of highway traffic consistent with the hydrodynamic theory,” Transportation Research Part B: Methodological, vol. 28, no. 4, pp. 269–287, 1994. View at: Publisher Site  Google Scholar
 B. Wolshon, “Onewayout: contraflow freeway operation for hurricane evacuation,” Natural Hazards Review, vol. 2, no. 3, pp. 105–112, 2001. View at: Publisher Site  Google Scholar
 H. Tuydes and A. Ziliaskopoulos, “Tabubased heuristic approach for optimization of network evacuation contraflow,” Transportation Research Record, vol. 1964, pp. 157–168, 2006. View at: Google Scholar
 H. Sbayti and H. S. Mahmassani, “Optimal scheduling of evacuation operations,” Transportation Research Record: Journal of the Transportation Research Board, vol. 1964, pp. 238–246, 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Zhengfeng Huang. 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.