Research Article  Open Access
Zhen Zhou, Ada Che, Feng Chu, Chengbin Chu, "Model and Method for Multiobjective TimeDependent Hazardous Material Transportation", Mathematical Problems in Engineering, vol. 2014, Article ID 832814, 11 pages, 2014. https://doi.org/10.1155/2014/832814
Model and Method for Multiobjective TimeDependent Hazardous Material Transportation
Abstract
In most of the hazardous material transportation problems, risk factors are assumed to be constant, which ignores the fact that they can vary with time throughout the day. In this paper, we deal with a novel timedependent hazardous material transportation problem via lane reservation, in which the dynamic nature of transportation risk in the reallife traffic environment is taken into account. We first develop a multiobjective mixed integer programming (MIP) model with two conflicting objectives: minimizing the impact on the normal traffic resulting from lane reservation and minimizing the total transportation risk. We then present a cutandsolve based constraint method to solve this model. Computational results indicate that our method outperforms the constraint method based on optimization software package CPLEX.
1. Introduction
Hazardous materials are a kind of goods with physical, chemical, and biological properties, which could cause lots of accidents, such as flame, explosion, and leak, in the course of production, storage, and transportation. With the development of industry, more and more hazardous substances, including raw materials, intermediate, final products, and wastes, are produced and moved daily through the transportation network in different transportation modes, such as road, rail, water, air, and pipeline. Their transportation risk imposed on environment and human is widening and deepening year by year. Hazardous material transportation has become a serious problem worldwide and attracts many researchers’ attentions in the related field.
It has been pointed out in much relevant literature that the essential objective of hazardous material transportation is to minimize the transportation risk due to the nature of this problem. As we know, the selection of routes in a network for hazardous material transportation can affect its risk factors, such as the probability of hazardous material accidents and the risk exposure to the surrounding population and environment. Therefore, appropriate routing decisions are very important for hazardous material transportation management. In the last couple of decades, various applications of operations research models to hazardous material transportation have focused on risk reduction and fruitful achievements in this area have been published; please see [1] for details. However, in almost all of the hazardous material transportation problems, the transportation risk is considered to be timeinvariant. That is to say, the risk of a road segment in a transportation network is assumed to be constant, which fails to capture the dynamic nature of the reallife traffic environment. In real life, risk on road segments is timedependent on population density subject to timeofday variation, peak and offpeak periods, various weather conditions, and so on. The timedependent risk is one of the important features of hazardous material transportation. The timedependent transportation problem is to decide the path for each shipment and its starting time so as to minimize the transportation risk.
An important branch of the timedependent transportation concerns hazardous materials. Timedependent transportation problems can be distinguished into deterministic and stochastic settings. The dynamic characteristic of transportation networks usually depends on one or more traffic quantities [2], such as travel time, link volume, and queue length. For deterministic timedependent transportation problems, part or all of the traffic quantities are assumed to be variant but all of the traffic quantities are known for a road segment. For stochastic timedependent transportation problems, the traffic quantities are considered as random variables with timedependent distribution functions, such as in [2–8]. For timedependent hazardous material transportation, traffic quantities are usually travel time and transportation risk. According to optimization criteria, timedependent transportation problems can be divided into singleobjective ones and multiobjective ones.
For the hazardous material transportation routing problem, a lot of works have concentrated on timeinvariant risk and travel time. However, timedependent hazardous material transportation problems have not been widely studied and only a few related publications can be found in the literature. Jia et al. [9] investigated a hazardous material transportation problem with deterministic timedependent risk for minimizing the transportation risk. The proposed model guaranteed the minimum distance between hazardous material shipments at any time. They transformed the considered problem into a number of timedependent shortest path problems for each truck and proposed an iterative heuristic. Erkut and Alp [10] proposed an integrated routing and scheduling problem for hazardous material transportation in a network with stochastic timedependent accident probability, population exposure, and travel time. The model aimed to minimize transportation risk while imposing a constraint on the total travel time of the shipment. Meng et al. [11] examined a similar problem with multiple objectives. They transformed it into a timedependent multiobjective shortest path problem subject to three types of time constraints. A dynamic programming approach was constructed to solve the problem. Both of the methods of references [10, 11] were pseudopolynomial. Nozick et al. [12] developed an approach for a hazardous material routing and scheduling problem with deterministic timedependent risk. But their approach could not guarantee generating all Paretooptimal paths. Chang et al. [13] proposed an effective algorithm for finding a path in stochastic timedependent networks that could address multiple optimization criteria. In their work, travel time, transportation risk, and other traffic quantities along paths were random variables. However, the performance of the algorithm was sensitive to some parameters and the computational burden increased with the number of the dominated paths. MillerHooks and Mahmassani [14] proposed an optimal routing algorithm for a single hazardous material shipment in a stochastic timedependent network, where the travel time and risk followed timedependent normal distribution functions. They presented several procedures for determining the best compromise path to minimize the total travel time and risk (population exposure).
Lane reservation is considered as a flexible and economic strategy for special events or situations, like sport games and emergencies. As stated in [15], it has been successfully applied to the Olympic Games in 2000 in Sydney [16] and in 2004 in Athens [17]. The principle of lane reservation is to reserve lanes on some road segments and/or in some specific time periods in a transportation network. Only special transportation tasks are allowed to pass through them. A probable consequence of this strategy is that it may worsen traffic congestions on general lanes, which is also an important feature of the lane reservation problem. Research works about lane reservation involving simulation tools and statistical methods can be found in the literature [18–20]. Some recent works focused on minimizing the traffic impact by linear programming models and optimal methods for various applications, such as automated truck, large sportive events, and hazardous material transportation [15, 21–24]. Our previous work [15] investigated a hazardous material transportation problem via lane reservation, in which the risk along a road segment was assumed to be constant. The goal was to minimize the transportation risk and the impact on normal traffic due to lane reservation. Remarkably, it has been shown in [15] that the lane reservation strategy can greatly reduce the transportation risk at a reasonable cost of its traffic impact. The problem proposed in this paper, in which the transportation risk is considered to be deterministic and timedependent, is more realistic. To the best of our knowledge, this is the first work for timedependent hazardous material transportation via lane reservation.
In this paper, we investigate a novel multiobjective hazardous material transportation problem via lane reservation with a deterministic timedependent transportation risk. As we know, the factors of the transportation risk generally include the hazardous material accident probability and the population exposure to the accidents. The accident probability estimation is influenced by the nature of roads, characteristics of the trucks, transportation environment, driver conditions [25], and so forth. Estimating the accident probability is a complicated and difficult work. For simplification, the probability of an accident is regarded to be timeinvariant in this paper. Population exposure is determined by population density and area. In real life, the population density along a road segment strongly depends on time and space. As we know, the population density in hospitals, schools, factories, and so on in day time is greater than that at night, and the opposite happens in residential areas. In this work, the accident probability on reserved lanes is assumed to be constant and the population exposure is assumed to be timedependent. Therefore, the transportation risk varies with time and space. This work is motivated by the dynamic characteristic of risk and it is a natural extension of our previous work.
The contributions of this work can be summarized as follows. Firstly, we propose a novel multiobjective MIP model for the timedependent hazardous material transportation problem via lane reservation. Secondly, we develop a cutandsolve based constraint method for the considered problem. Computational results show the effectiveness of the method.
The remainder of this paper is organized as follows. Section 2 describes the timedependent hazardous material transportation problem via lane reservation and a new multiobjective model is formulated with the objectives of minimizing the total impact on the normal traffic and the total transportation risk. In Section 3, some properties are analyzed to reduce the search space of solutions and a cutandsolve based constraint method is developed. Section 4 reports the computational results of experiments. Section 5 concludes the paper and discusses future research directions.
2. Problem Formulation
2.1. Problem Description
Consider a bidirectional graph , where is the set of nodes and is the set of arcs. The graph represents the transportation network in which the vehicles carrying hazardous materials are allowed to move. Arc denotes a road segment from node to node . kinds of hazardous materials will be transported from origin nodes to destination nodes .
Before formulation of the problem, we make the following nonrestrictive assumptions. There are at least two lanes on a road segment such that one lane is allowed to be reserved. From the point of view of transport safety, any path for hazardous material shipments consists of only reserved lanes. Vehicles with hazardous materials travel on the reserved lanes without congestion. Consequently, travel time on reserved lanes is timeinvariant. Any two hazardous material shipments on the same road segment must maintain a minimum time interval, called safety time interval. The accident probability on a reserved lane is constant and hazardous material accidents happen independently.
In the network, nonreserved lanes are called general lanes, as shown in Figure 1. From assumption , the travel time on a reserved lane is timeinvariant throughout the day. Nevertheless, the risk on a road segment should be timedependent because the population exposure varies with time in nature. The population exposure on each arc at time period , denoted as , depends on the departure time from node . Without loss of generality, set as the beginning time of the first period. In a day, there are usually only several time periods [26]. So travel time on the reserved lane is far less than the length of a time period; that is, . The problem is to choose lanes on the existing network to be reserved, select the path for each hazardous material shipment, and decide the travel time period for each shipment on its selected path. The objective of the problem is to seek the best tradeoff for minimizing the total traffic impact on the normal traffic and the total transportation risk.
2.2. Problem Definition
To model our problem, we introduce the following notation.
Sets and Parameters : The set of hazardous material shipments : The set of origin nodes : The set of destination nodes : The set of time periods : The safety time interval between any two shipments which pass the same arc : The travel time on the reserved lane on arc : The impact on the normal traffic due to the reserved lane on arc : The accident probability of hazardous material on a reserved lane on arc : The population exposure along arc at time period .
Decision Variables is the arrival or departure time of shipment at node . Note that if shipment does not pass node , , ,
2.3. Mathematical Model
The mathematical model for the hazardous material transportation problem via lane reservation in a timedependent network is presented by constraints (3)–(17). Consider the following:
: subject to where is a very large positive number.
Objective (1) is to minimize the total impact on the normal traffic. The impact can be considered as the increase in travel time on the general lane(s) because of the lane reservation strategy, which is proportional to the travel time on the general lane(s) of arc before lane reservation and inversely proportional to the total number of lanes on arc [13]. We use the same formula in [21] to evaluate the impact , where is the travel time on the general lane(s) of arc and is the total number of lanes in the road segment. Statistical results from [27] revealed that the travel time on the general lanes increased by about 53% after reserving one of three lanes on A1 highway in Paris. This figure is very close to the theoretical test result (50%) obtained by the computational experiment done in [22] using the above formula to estimate the impact. If no lane is reserved, then . Objective (2) is to minimize the total transportation risk. In this paper, the risk is measured in a traditional way, and it is measured by multiplying the probability of the hazardous material accident by its consequence. See [1] for more details. Constraint (3) (resp., (4)) means that, for shipment , there is one and only one path departing from the source node (resp., arriving at the destination node ) during one and only one time period. Constraint (5) ensures the flow conservation constraint for node in on space and time. It represents that if shipment arrives at a node via a reserved lane during time period , it must also depart from via a reserved lane during time period . Constraints (3), (4), and (5) together ensure that there is one and only one path for each shipment from its origin to its corresponding destination during one and only one time period. Constraint (6) guarantees that no shipment would pass through arc during any time period if no lane on the arc has been reserved. Constraints (7) and (8) mean that if shipment passes through the reserved lane on arc , then its travel time is . Constraint (9) means that there is exactly one time period for on any node . In a feasible solution, shipment passes through arc if and only if two conditions are satisfied: and . Constraints (9) and (10) imply that departure time should be located within one and only one time period and on exactly one arc. Constraints (11) and (12) guarantee that if two or more shipments pass the same reserved lane, then the safety time interval between any two shipments must be satisfied. Constraints (13)–(17) specify the restriction on the decision variables.
3. Solution Algorithm
In this section, an improved algorithm, called the cutandsolve based constraint method, is developed for solving the multiobjective model. There are several techniques to solve a multiobjective problem in the literature, such as the weightedsum method, the constraint method, the goal attainment approach, and metaheuristics [28]. With the constraint method, a multiobjective problem can be converted into a series of singleobjective problems to obtain Paretooptimal solutions. Moreover, this method can alleviate the difficulties in setting up an appropriate weight vector faced by the weighted sum method and the goal attainment approach. Therefore, in this paper, the constraint method is used to solve the considered problem.
3.1. Constraint Method
Introduced by Haimes et al. [29], the constraint method is based on a scalarization where only one objective function (commonly, it may be the most preferred or primary one) is minimized while the others are bounded by some allowable values , , and added to the original model as constraints. The original multiobjective optimization problem is converted into a series of singleobjective optimization problems as given below:
: subject to where is the feasible solution space.
According to [30], for a given vector of , optimal solutions of are weakly Paretooptimal solutions. For solving problem , it is necessary to determine the range of that can be defined by ideal and nadir points of the original problem [31]. The ideal and nadir points can be obtained from the solutions of subproblems using the procedure described below.
The proposed mathematical model in Section 2 has conflicting objective functions (1) and (2) subject to constraints (3)–(17). From assumption , any path for hazardous material shipments contains only exclusive reserved lanes. Therefore, in this paper, the objective of minimizing the traffic impact is considered as the main objective function. In this way, the original multiobjective model can be converted into a singleobjective one, in which the traffic impact is minimized, while the transportation risk subject to forms a new constraint and is added to the original model.
The procedure of the constraint method is described as follows [15, 29].
Step 1. Transform problem into problem , which aims at minimizing objective function (1) subject to constraints (3)–(17) and
Step 2. Determine the ideal point () and the nadir point ().
The ideal point corresponds to lower limits of , where subject to constraints (3)–(17). The nadir point indicates upper limits of and , respectively, where subject to constraints (3)–(17) and . And subject to constraints (3)–(17) and . That is, the values of the objective function are bounded by .
Step 3. Fix the values of .
Determine the range of :
The range is divided into equal intervals with points, called equidistant grid points. The value of in constraint (20) is defined by these grid points with the following formula:
Step 4. Repeat to solve problem with and obtain Paretooptimal solutions.
3.2. CutandSolve Method
With the constraint method, the considered problem is transformed into a series of singleobjective MIP problems . In this section, we propose an exact and iterative approach, called cutandsolve (CS), to solve . The CS method was first presented by Climer and Zhang and was proved to be efficient for the traveling salesman problem [32]. It is a special branchandbound search strategy but can avoid making wrong choices in depthfirst branchandbound. Given an integer programming problem (IP) with the minimization of an objective function, at the th iteration of the branching tree of CS, there are only two nodes, corresponding to sparse problem () and remaining problem (), respectively [33]. The can be optimally solved in a reasonable amount of computation time because it is relatively small. The optimal solution of , if it exists, provides an upper bound of the IP, denoted as . And if it is “good” enough, then the best upper bound of the IP found so far, denoted as , is updated. Meanwhile, a lower bound of the IP, denoted as , can be obtained by solving the linearly relaxed . When the best upper bound found so far is smaller than or equal to , a global optimal solution is obtained and the CS iteration stops. Otherwise, the current is further divided into a new and a new by adding a piercing cut for the next iteration. The process continues until an optimal solution to the original IP is obtained.
(1) Preprocessing. Before using the cutandsolve method, we analyze properties of the model in order to reduce the search space. If some variables can be fixed in advance, the reduction in the number of decision variables will possibly speed up the CS process.
As defined in Section 2, means that shipment passes through the reserved lane on arc when time occurs at time period ; means that the arriving time of shipment at node occurs at time period and otherwise . We have the following property.
Property 1. If , then , , , , .
Proof. Note that if , then . It implies that . That is, . For arc , we can deduce that . Note that . We thus have , which implies two cases. In one case, , and, in the other case, . That is to say, either or . It follows that .
From Property 1, constraint (23) will be added to problem to obtain problem , which aims at minimizing objective function (1) subject to constraints (3)–(17), (20) and
(2) Definition of Piercing Cut, Remaining Problem, and Sparse Problem. A key issue for the cutandsolve method is to find an appropriate piercing cut that separates the current remaining problem into a new sparse problem and a new remaining problem. The reasons are the following. If the solution space of is too small, the optimal solution is not “good” enough to update the upper bound; if the solution space of is too large, it will take too much time to obtain an optimal solution. Another important factor which influences the efficiency of the cutandsolve method is that tight should be provided in each iteration. For an integer programming model, is usually obtained by solving a linear relaxation problem of . Note that means that shipment passes through the reserved lane on arc during time period and, otherwise, . If are not relaxed, our preliminary simulation experiments show that most of them take the values of zero, which means that very few reserved lanes are passed by the shipments at any time period. The obtained in this case is usually very “bad.” Therefore, should be relaxed to obtain more subpaths. Similar observation can be also found for . Another important decision variable means whether there is a reserved lane on arc . Our preliminary simulation experiments show that if are relaxed the same as , a considerable number of take the values greater than zero. It means that there may exist a reserved lane on the corresponding arc. However, the experiments also show that if are not relaxed, a better can be obtained. Therefore, are not relaxed as continous variables in this paper. Similar observation can be also found for . In a word, in order to obtain an improved , we apply the partial integrality strategy for , in which only and are relaxed as continuous variables while the integrality of and is maintained.
Climer and Zhang introduced a general procedure for generating piercing cuts based on reduced cost from an optimal solution of linear relaxation problem [32]. In their work, was defined as a set including the decision variables with large reduced cost. But the general procedure is not appropriate to MIP because it has been shown by our preliminary experimental results that the lower bound of the proposed problem obtained by the linearly relaxed is not good enough and the reduced costs of decision variables are often missing. Reference [24] proposed a new piercing cut technique for MIP using the number of “critical links.” Since the value of may be fractional in an optimal solution of the relaxed , , there may be multiple paths for some shipments. The link with the greatest value of is called “critical link.” For more details of the piercing cut based on “critical link,” please see [24]. The considered problem in this paper is different from that in [24]. The piercing cut in [24] cannot be used directly for our problem. A new piercing cut based on “critical link” is proposed as follows.
Let represent the set of shipments which contains more than one path in the fractional solution and let represent the first node where multipath appears for shipment . Define as the set of the most potential arcs for all shipments in , which implies that the arcs in are very likely to be selected in the final optimal solution of . We have where refers to the arc with the largest value among all the arcs outgoing from node .
In [32], the piercing cut is considered as a combination of some decision variables in a certain set, called . In this paper, refers to the set of these decision variables , where is a shipment which has multiple paths in the fractional solution and () is the most potential arcs for shipment . Set is defined as follows:
The piercing cut () is defined as follows: where is a given integer in .
Accordingly, the additional constraint associated with can be written as follows:
As mentioned above, and are generated by adding constraints (26) and (27) to , respectively, that is, : objective function (1) subject to constraints (3)–(17), (20), (23), (26), and : objective function (1) subject to constraints (3)–(17), (20), (23), (27), and (28).
It is worth pointing out that when , the original problem is considered as CP_{1}. Hence, constraint (28) should be removed for RP_{1} and SP_{1}.
The process of the proposed constraint is illustrated in Figure 2.
4. Computational Results
To evaluate the efficiency of the proposed method, 155 instances (31 sets × 5 instances) are randomly generated. The proposed algorithm is coded in . The computational experiments are carried out on an HP PC with 3.10 GHz Intel Core processor and 4 GB RAM under Windows 7 environment. The CPLEX MIP solver (version 12.5) under default settings is used to solve . It is allowed to run until problem is solved to optimality.
The transportation network in this work is generated according to the random network topology generator proposed by Waxman [34]. The nodes are randomly and uniformly generated in the plane , while arc is generated by a probability function , where and are the Euclidean distance between and and the maximum distance between any two nodes in the graph, respectively, and . The origin and destination pairs are randomly selected from the set of nodes. Parameter is set from 1 to 5 because in real life the number of time periods is not very large [26]. Let be the accident probability of hazardous material happening on the general lane(s) of arc . Note that . Our preliminary experimental results also show that the values of the input parameters , , , , , , and have little effect on the performance of the proposed algorithm for a considerable number of instances and they are generated according to the following ways. The data of parameters are set with the similar way of [15]: and , where is a uniform distribution; ; , whose unit is 10^{−7}; is generated by , whose unit is 10^{4}; ; and is generated by . The number of iterations of constraint method, , is set to 20.
In this paper, represents the average degree of graph , where and are its number of nodes and arcs, respectively [35, 36]. The average degree of graph is defined as its number of arcs per node, which implies the density of the graph.
Table 1 summarizes the computational results on the randomly generated instances with , , and . The total computation time for an instance represents its total running time for obtaining 21 solutions in the constraint method. Columns and represent the average computation time (in CPU seconds) of five instances for each set by the proposed constraint method, in which singleobjective problem is solved by the optimization software package CPLEX and the CS based method, respectively. We can observe from Table 1 that the proposed CS based method is more efficient than CPLEX and both of the total computation times moderately increase with the number of nodes. It is worth pointing out that the trends of two curves of and are almost the same. But increases with the number of nodes more slightly than in Figure 3. For example, for set 1 with , is only 1.261 times as much as , whereas, for set 9 with , the ratio / is 2.223 times.

Table 2 summarizes the computational results on the randomly generated instances with , , and . It can be observed from Table 2 that the computation times of CPLEX and the proposed method drastically increase with the number of nodes, but the latter increases more slightly than the former. Take sets 10 and 11 for example; given the number of shipments and time periods, the computation times and for set 11 are 6.636 and 4.440 times as much as those for set 10, respectively. From Tables 1 and 2, we can also observe that the CS based method is more efficient for the sets with than for the sets with . For example, given the number of nodes and time periods, for set 3 is only 1.476, while it increases to 3.573 for set 12. When the number of nodes increases to 60 in Table 2, CPLEX cannot find 21 optimal solutions of within 36000 s, but the proposed method can do it well for sets with up to 80 nodes. We set a threshold for the computation time for each problem to 36000 s/21. When the threshold is reached and a problem is not solved to optimality, CPLEX is terminated. For sets 14–16, although CPLEX fails to find all optimal solutions, it can provide lower bounds and upper bounds of which are not solved to optimality when CPLEX is terminated. The gap between the lower bound and upper bound is denoted by (upper bound − lower bound)/upper bound and it implies the extent of optimality. GAP in Table 2 denotes the average gap for problems that are not solved to optimality. The average gaps of sets 14–16 are 5.309%, 3.217%, and 3.949%. The gaps are relatively small, which means that the obtained solutions might be relatively close to optimality for sets 14–16. In Table 2, given that , the computation time for set 16 is nearly up to 36000 s so that this set can be considered as one of the largestscale problems which can be solved in reasonable time.

Table 3 shows the total computation times of the proposed method with , , and different . We can observe from Table 3 that, for a given number of , the computation time increases quickly with the value of . For example, the computation times for 100 nodes with (set 30) and (set 31) are 11.329 and 24.237 times as much as that with (set 29), respectively. As shown in Tables 1 and 3, the computation time also increases with the average degree of network. When , given the number of nodes and shipments, the computation time with is more than that with . For example, the computation time for set 24 is 2.645 times as much as that for set 5.

5. Conclusion
This paper investigated the timedependent lane reservation problem for hazardous material transportation, in which the transportation risk varies with time throughout the day. A new multiobjective mixed integer programming model was first proposed for the considered problem. Then a cutandsolve method based constraint method was developed. For improving the cutandsolve method, a specific property of the considered problem was explored. Computational results showed that our method outperforms the optimization software package CPLEX.
For the considered problem, one of the future research directions is to study properties of the model and attempt to reduce the search space of solution. Exploiting efficient piercing cuts is another way to improve the efficiency of the cutandsolve method. Finally, developing high quality heuristics is a future direction for largescale timedependent hazardous material transportation problems via lane reservation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (71071129 and 71471145), the Natural Science Foundation of Shaanxi Province, China (2012JM9001), Humanities, Social Sciences and Management Innovation Foundation of Northwestern Polytechnical University (RW201301), 1000 Plan for Foreign Talent (WQ20123400070), Foreign Experts HighEnd Projects (GDW20123400148), and the Ministry of Education Fund for Doctoral Student Newcomer Awards of China.
References
 E. Erkut, S. A. Tjandra, and V. Verter, “Hazardous materials transportation,” in Handbook in OR&MS, vol. chapter 9, pp. 539–621, Elsevier, Amsterdam, The Netherlands, 2007. View at: Google Scholar
 S. Gao and I. Chabini, “Optimal routing policy problems in stochastic timedependent networks,” Transportation Research Part B: Methodological, vol. 40, no. 2, pp. 93–122, 2006. View at: Publisher Site  Google Scholar
 E. D. MillerHooks and H. S. Mahmassani, “Least possible time paths in stochastic, timevarying networks,” Computers & Operations Research, vol. 25, no. 12, pp. 1107–1125, 1998. View at: Publisher Site  Google Scholar
 E. D. MillerHooks and H. S. Mahmassani, “Least expected time paths in stochastic, timevarying transportation networks,” Transportation Science, vol. 34, no. 2, pp. 198–215, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. MillerHooks, “Adaptive leastexpected time paths in stochastic, timevarying transportation and data networks,” Networks, vol. 37, no. 1, pp. 35–52, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 E. MillerHooks and H. Mahmassani, “Path comparisons for a priori and timeadaptive decisions in stochastic, timevarying networks,” European Journal of Operational Research, vol. 146, no. 1, pp. 67–82, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Opasanon and E. MillerHooks, “Multicriteria adaptive paths in stochastic, timevarying networks,” European Journal of Operational Research, vol. 173, no. 1, pp. 72–91, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 S. Wang, M. Wang, and J. Hu, “Shortest path problem in stochastic timevarying transportation networks considering risk measurement,” International Journal of Advancements in Computing Technology, vol. 3, no. 6, pp. 169–175, 2011. View at: Publisher Site  Google Scholar
 H. Jia, L. Zhang, M. Duan, and G. Fu, “A timedependent flow model for Hazmat transportation routing,” in Proceedings of the IEEE International Conference on Industrial Engineering and Engineering Management (IEEM '10), pp. 1553–1557, Macau, China, December 2010. View at: Publisher Site  Google Scholar
 E. Erkut and O. Alp, “Integrated routing and scheduling of hazmat trucks with stops en route,” Transportation Science, vol. 41, no. 1, pp. 107–122, 2007. View at: Publisher Site  Google Scholar
 Q. Meng, D.H. Lee, and R. L. Cheu, “Multiobjective vehicle routing and scheduling problem with time window constraints in hazardous material transportation,” Journal of Transportation Engineering, vol. 131, no. 9, pp. 699–707, 2005. View at: Publisher Site  Google Scholar
 L. K. Nozick, G. F. List, and M. A. Turnquist, “Integrated routing and scheduling in hazardous materials transportation,” Transportation Science, vol. 31, no. 3, pp. 200–215, 1997. View at: Publisher Site  Google Scholar
 T.S. Chang, L. K. Nozick, and M. A. Turnquist, “Multiobjective path finding in stochastic dynamic networks, with application to routing hazardous materials shipments,” Transportation Science, vol. 39, no. 3, pp. 383–399, 2005. View at: Publisher Site  Google Scholar
 E. MillerHooks and H. S. Mahmassani, “Optimal routing of hazardous materials in stochastic, timevarying transportation networks,” Transportation Research Record, no. 1645, pp. 143–151, 1998. View at: Google Scholar
 Z. Zhou, F. Chu, A. Che, and M. Zhou, “εConstraint and fuzzy logicbased optimization of hazardous material transportation via lane reservation,” IEEE Transactions on Intelligent Transportation Systems, vol. 14, no. 2, pp. 847–857, 2013. View at: Publisher Site  Google Scholar
 J. Black, “Strategic transport planning, demand analysis of transport infrastructure and transport services for the 27th summer, olympiad held in Sydney, Australia, 2000,” Journal of Transportation Systems Engineering and Information, vol. 2, no. 2, pp. 14–30, 2004. View at: Google Scholar
 E. Zagorianakos, “Athens 2004 olympic games’ Transportation Plan: a missed opportunity for Strategic Environmental Assessment (SEA) integration,” Journal of Transport Geography, vol. 12, no. 2, pp. 115–125, 2004. View at: Publisher Site  Google Scholar
 A. F. Abdelghany, K. F. Abdelghany, H. S. Mahmassani, and P. M. Murray, “Dynamic traffic assignment in design and evaluation of highoccupancy toll lanes,” Transportation Research Record, no. 1733, pp. 39–48, 2000. View at: Google Scholar
 P. M. Murray, H. S. Mahmassani, and K. F. Abdelghany, “Methodology for assessing highoccupancy tolllane usage and network performance,” Transportation Research Record, vol. 1765, pp. 8–15, 2001. View at: Google Scholar
 A. S. Shalaby, “Simulating performance impacts of bus lanes and supporting measures,” Journal of Transportation Engineering, vol. 125, no. 5, pp. 390–397, 1999. View at: Publisher Site  Google Scholar
 Y. Wu, C. Chu, F. Chu, and N. Wu, “Heuristic for lane reservation problem in time constrained transportation,” in Proceedings of the IEEE International Conference on Automation Science and Engineering (CASE '09), pp. 543–548, Bangalore, India, August 2009. View at: Publisher Site  Google Scholar
 Y. Fang, F. Chu, S. Mammar, and A. Che, “An optimal algorithm for automated truck freight transportation via lane reservation strategy,” Transportation Research Part C: Emerging Technologies, vol. 26, pp. 170–183, 2013. View at: Publisher Site  Google Scholar
 Y. Fang, F. Chu, S. Mammar, and M. Zhou, “Optimal lane reservation in transportation network,” IEEE Transactions on Intelligent Transportation Systems, vol. 13, no. 2, pp. 482–491, 2012. View at: Publisher Site  Google Scholar
 Y. Fang, F. Chu, S. Mammar, and A. Che, “A cutandsolvebased algorithm for optimal lane reservation with dynamic link travel times,” International Journal of Production Research, vol. 52, no. 4, pp. 1003–1015, 2014. View at: Publisher Site  Google Scholar
 Y. Qiao, N. Keren, and M. S. Mannan, “Utilization of accident databases and fuzzy sets to estimate frequency of HazMat transport accidents,” Journal of Hazardous Materials, vol. 167, no. 1–3, pp. 374–382, 2009. View at: Publisher Site  Google Scholar
 H. L. Khoo, L. E. Teoh, and Q. Meng, “A biobjective optimization approach for exclusive bus lane selection and scheduling design,” Engineering Optimization, vol. 46, no. 7, pp. 987–1007, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 J. Princeton and S. Cohens, “Impact of a dedicated lane on the capacity and the level of service of an urban motorway,” Procedia—Social and Behavioral Sciences, vol. 16, pp. 196–206, 2011. View at: Publisher Site  Google Scholar
 V. T'kindt and J.C. Billaut, “Multicriteria scheduling problems: a survey,” RAIRO Operations Research, vol. 35, no. 2, pp. 143–163, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Y. Haimes, L. S. Lasdon, and D. A. Wismer, “On a bicriterion formulation of the problems of integrated system identification and system optimization,” IEEE Transactions on Systems, Man and Cybernetics, vol. SMC1, no. 3, pp. 296–297, 1971. View at: Publisher Site  Google Scholar  MathSciNet
 V. Chankong and Y. Haimes, Multiobjective Decision Making Theory and Methodology, Elsevier Science, New York, NY, USA, 1983.
 J.F. Bérubé, M. Gendreau, and J.Y. Potvin, “An exact εconstraint method for biobjective combinatorial optimization problems: application to the traveling salesman problem with profits,” European Journal of Operational Research, vol. 194, no. 1, pp. 39–50, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. Climer and W. Zhang, “Cutandsolve: an iterative search strategy for combinatorial optimization problems,” Artificial Intelligence, vol. 170, no. 89, pp. 714–738, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Yang, F. Chu, and H. Chen, “A cutandsolve based algorithm for the singlesource capacitated facility location problem,” European Journal of Operational Research, vol. 221, no. 3, pp. 521–532, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 B. M. Waxman, “Routing of multipoint connections,” IEEE Journal on Selected Areas in Communications, vol. 6, no. 9, pp. 1617–1622, 1988. View at: Publisher Site  Google Scholar
 J. Li, F. Chu, and C. Prins, “Lower and upper bounds for a capacitated plant location problem with multicommodity flow,” Computers & Operations Research, vol. 36, no. 11, pp. 3019–3030, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 R. Diestel, Graph Theory, Springer, Berlin, Germany, 2005. View at: MathSciNet
Copyright
Copyright © 2014 Zhen Zhou et al. 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.