Research Article  Open Access
Optimizing Route for Hazardous Materials Logistics Based on Hybrid Ant Colony Algorithm
Abstract
Optimizing Route for Hazardous Materials Logistics (ORHML) belongs to a class of problems referred to as NPHard, and a strict constraint of it makes it harder to solve. In order to dealing with ORHML, an improved hybrid ant colony algorithm (HACA) was devised. To achieve the purpose of balancing risk and cost for route based on the principle of ACA that used to solve TSP, the improved HACA was designed. Considering the capacity of road network and the maximum expected risk limits, a route optimization model to minimize the total cost is established based on network flow theory. Improvement on route construction rule and pheromone updating rule was adopted on the basis of the former algorithm. An example was analyzed to demonstrate the correctness of the application. It is proved that improved HACA is efficient and feasible in solving ORHML.
1. Introduction
Hazardous materials, which have different physical and chemical properties, have high risk during transportation, as a series of problems may arise in this process. Route optimization is a complex combinatorial optimization problem, which is a typical NPcomplete problem and difficult to come up with a direct answer. It is a practical problem in urgent need of solution in which we can find the optimal plan under the restrictions quickly, accurately, safely, and economically.
Optimizing Route for Hazardous Materials Logistics (ORHML) can be described as follows. Given a set of hazardous materials and an underlying network consisting of a number of nodes and capacitated arcs, we wish to find an optimal routing plan to ship the hazardous materials through the network at lowest cost without violating the capacity limits. ORHML models also appear as subproblems in more complicated models, such as distribution system design and capacitated network design.
ORHML has attracted the attention of many OR researchers. Kara et al. [1] presented several route planning models of road. Verma and Verter [2] gave a number of route planning models of rail. Iakovou [3] provided route planning models of water. MillerHooks [4] modeled ORHML as a path selection problem in a stochastic timevarying network. Dell’Olmo et al. [5] finding a number of spatially dissimilar paths between an origin and a destination can also be considered in this area. Jin and Batta [6] presented a risk model that took into account the dependency to the impedances of preceding road segments. Erkut and Verter [7] proposed a collection of edges in place of an origindestination route for a hazmat shipment, where travel on this path can be viewed as a probabilistic experiment. It considered minimizing for a given OD pair in a hazmat transport network is a shortest path problem which is solved easily for even large networks. Their work also pointed out that this approximation is likely to result in a very small error in measuring the incidence probability along a hazmat transport route [8].
In general, in spite of their more realistic assumption, most of the exact versions of risk models have some puzzling properties and these models may not be suitable for hazmat transportation planning. We suggest that researchers and practitioners consider the properties of the risk models carefully.
Successful ant colony algorithms have been developed for several combinatorial optimization problems, such as Traveling Salesman Problem (TSP) and Vehicle Routing Problem (VRP).
Based on the successful application of ant colony algorithm to TSP, a new ACA was adopted in solving ORHML. A collaborative metaheuristic approach is inspired by the foraging behavior of real colonies of ants. The basic ACA application is a large number of simple artificial agents that are able to build good solutions for hard combinatorial optimization problems through lowlevel based communications. Real ants cooperate in their search for food by depositing chemical traces (pheromones) on the floor.
2. Model Analysis
ORHML optimization model can be demonstrated as follows: where represents the actual road network, in which the set of nodes represents a road intersection, and arc set indicates the connection between the intersection of sections. If denotes the sections, then obviously and . represent a different flow of hazardous materials set ; represents the hazardous materials . represents hazardous materials using the arc ; represents the unit cost when hazardous materials use arc ; represents the unit risk when hazardous materials use arc . represents the capacity constraints of the road. represents the nodearc incidence matrix as the parameter [9].
represents the supply amount or demand amount in point of hazardous materials , means that the supply, which means that the demand.
Based on the risk analysis theory, the risk of hazardous materials in can be calculated by the following formula [10]: where represents the traffic accident of hazardous materials in ; represents the leak accident of hazardous materials in ; represents the accident loss of hazardous materials in . = {casualties along the road, casualties in the vehicle, property damage, traffic interruption, evacuation, environmental damage}. represents the risk of loss caused by the harm .
(1) Casualties along the road : where represents the length of the road; represents the radius of the influence area for the accident disaster; represents the population density along the road; represents the value of the loss of life.
(2) Casualties in the vehicle : where represents the traffic flow of bus; represents the traffic flow of car; represents the patronage of bus; represents the patronage of car.
(3) Property damage : where represents the loss value of bus; represents the loss value of car.
(4) Traffic interruption : where represents the traffic disruption time; represents per unit time; represents the value of loss time.
(5) Evacuation : where represents average evacuation costs.
(6) Environmental damage : where represents value of environmental damage of unit area.
3. Solving Strategies
HACA goal is to find the shortest tour. In HACA ants build tours in parallel, where is a parameter. Each ant is randomly assigned to a starting node and has to build a solution, that is, a complete tour. A tour is built node by node: each ant iteratively adds new nodes until all nodes have been visited. Improvement process of solution strategy is present in the algorithm; this strategy can transform solution of the problem to a single ant path. Given multiple source points, ink points, and the corresponding flow, construct an ant colony; in the initial state, the ants are randomly distributed in each source point; the ants’ next point is selected at random. Updated according to the flow from the source point to sink point as weight. When the ant went to the source point and the corresponding sink point, place the ants to another source point, begin the next search process until all sources and sinks were finished, and then complete a cycle of the ant. As can be seen in Figure 1 [11].
When ant is located in node , it chooses the next node probabilistically in the set of feasible nodes according to .
In the HACA, original version formula for is
The number of next nodes which ant in node could select is represented by . In basic ant colony algorithm, , denotes the number of nodes in set . The value of should be adjusted dynamically:
The sum of ants is represented by ; there are trails from node ; the number of ants that pass node is , and is distributing in trails with number :
Using for reference the idea of ant algorithm of Dorigo, the meliorative route construction rule can be expressed as follows , called saving; . Variable induced considering the time constraint.
The algorithm constructs two ant colonies to realize the goal of optimization: minimize the cost and risk of the hazardous materials transportation. These two ant colonies are the HACARISK and HACACOST, which are used to determine the optimal route. In consideration of the two together, we may find better solutions to the problem. After getting the improved solution, the overall pheromone update would be used to get the exchange of information concerning the advantages and disadvantages of the solution [12].
The HACARISK optimization process is to get the feasible solution which uses one less risk than . HACARISK process is used to determine a shorter path under the circumstances. Once a new better feasible solution has been figured out, should be updated. These two processes will come to an end if a solution with the use of fewer cost comes up. At the same time, a new ant colony can be built on the basis of this new solution, which restarts a new interactive optimization process.
If the least cost has been figured out by this algorithm, the HACARISK optimization process can be stopped, and then we shift to the HACACOST optimization process.
In the HACACOST we introduce an integer vector ; means the times that the point has not been covered by the solution. This parameter should be updated correspondingly with the optimal solution. Whenever a new solution has been found, the overall pheromone update should be applied to reset .
The HACACOST optimization process is similar to the traditional HACA optimization process which is to optimize the utilization of the load capacity and volume of RISKs.
The optimizing process can be divided into two parts: first, the ants move between the various points in search of the optimal solution; second, decompose the optimal solution into subsets according to the location of the original node, then allocate goods to the RISKs.
The movement of ants in the HACACOST and HACARISK is similar to each other. At the beginning of the algorithm, ants are randomly distributed among all the nodes. From the initial node, the ants determine the next move within the scope of available points, on the basis of the probability and the restrictions on RISK load capacity and volume. The time and cost of all nodes to the virtual origin node are 0, and no direct connection between the virtual origin nodes is allowed. Once the ants arrive the terminal node, thus begins the construction of the next subset. After traversing all the nodes, the ants return to the initial node to form a loop, that is, feasible solution.
In the HACACOST optimization process, a feasible solution is formed through the collaboration of all the ant colonies. Compare with , if better, update with the value of .
When the above process was being completed, the obtained solution may still be incomplete (and may have missed points), so the insertion process is needed as further optimization. Insertion process is mainly to deal with the point that has not been included in the existing solution. The process determines the optimal feasible insertion position (the shortest path) for each point until the plug cannot find a viable location [13].
4. Example and Analysis
As can be seen in Figure 2, there are two kinds of hazardous materials (flammable liquid) and (poisonous gas) : , from note A to note H, flow is 60 units; , from note B to note G, flow is 20 units. That is, b(A) = 60, b(H) = −60, b(B) = 20, and b(G) = −20. The related data values are shown in Table 1.

Set freight = 0.5 yuan/km/unit, yuan/person, persons/vehicle, persons/vehicle, yuan, yuan, h, h, yuan, and yuan/h.
, km; , km. The probabilities of hazard losses can be seen in Table 2.

Total expected risk is shown in Table 3.

Set Fright = 0.5 yuan/km/unit. The cost can be shown in Table 4.

In view of complexity of this problem, the whole level optimization solutions are obtained based on the random searching and evolution process of the improved ant colony algorithm. The solution can be seen in Figure 3.
5. Summary
Solving the ORHML is to determine route plans that can make a minimum of risk and safety under the premise of restrictions of the placements. That is, to seek a safe routing plan that could obtain the lowest cost with the conditions has been known.
A new ant colony optimization based approach to solve ORHML was introduced. In particular, the algorithm has been designed to solve ORHML with a balance between the safe and the cost. Pinpointing the characteristics of this problem, our algorithm introduces a new methodology for optimizing multiple objective functions. We consider the optimization of the safe and the cost at the same time. This paper analyzes the differences of the ant colony algorithms in solving the restricted ORHML. In order to deal with ORHML, an improved HACA was devised. Pinpointing the characteristics of this problem, we consider the optimization of cost and risk together. The basic idea is to coordinate the activity of different ant colonies, each of them optimizing a different objective. These colonies work by using independent pheromone trails but they collaborate by exchanging information. Furthermore, the integrated use of HACASAFE and HACACOST has been applied as the improvement of the solving strategy.
Finally, the feasibility and effectiveness of this method has been scrutinized with practical examples. From the result on the test problem, we can conclude that the model and the heuristic procedures are quite successful in solving ORHML.
Acknowledgment
The research is supported by Basic Scientific Research Funding of Beijing Jiaotong University (Project name: Collaborative optimization for hazardous materials transportation route choice and logistics center location).
References
 B. Y. Kara, E. Erkut, and V. Verter, “Accurate calculation of hazardous materials transport risks,” Operations Research Letters, vol. 31, no. 4, pp. 285–292, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Verma and V. Verter, “Railroad transportation of dangerous goods: population exposure to airborne toxins,” Computers and Operations Research, vol. 34, no. 5, pp. 1287–1303, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. T. Iakovou, “An interactive multiobjective model for the strategic maritime transportation of petroleum products: risk analysis and routing,” Safety Science, vol. 39, no. 12, pp. 19–29, 2001. View at: Publisher Site  Google Scholar
 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  Zentralblatt MATH  MathSciNet
 P. Dell'Olmo, M. Gentili, and A. Scozzari, “On finding dissimilar Paretooptimal paths,” European Journal of Operational Research, vol. 162, no. 1, pp. 70–82, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. Jin and R. Batta, “Objectives derived from viewing Hazmat shipments as a sequence of independent Bernoulli trials,” Transportation Science, vol. 31, no. 3, pp. 252–261, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Erkut and V. Verter, “Modeling of transport risk for hazardous materials,” Operations Research, vol. 46, no. 5, pp. 625–642, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Erkut, S. A. Tjandra, and V. Verter, “Chapter 9 hazardous materials transportation,” Handbooks in Operations Research and Management Science, vol. 14, pp. 539–621, 2007. View at: Publisher Site  Google Scholar
 V. Marianov and C. ReVelle, “Linear, nonapproximated models for optimal routing in hazardous environments,” Journal of the Operational Research Society, vol. 49, no. 1, pp. 157–164, 1998. View at: Google Scholar  Zentralblatt MATH
 V. Verter and B. Y. Kara, “A GISbased framework for hazardous materials transport risk assessment,” Risk Analysis, vol. 21, no. 6, pp. 1109–1120, 2001. View at: Publisher Site  Google Scholar
 J. E. Bell and P. R. McMullen, “Ant colony optimization techniques for the vehicle routing problem,” Advanced Engineering Informatics, vol. 18, no. 1, pp. 41–48, 2004. View at: Publisher Site  Google Scholar
 S. Ghafurian and N. Javadian, “An ant colony algorithm for solving fixed destination multidepot multiple traveling salesmen problems,” Applied Soft Computing Journal, vol. 11, no. 1, pp. 1256–1262, 2011. View at: Publisher Site  Google Scholar
 M. H. Afshar, “A parameter free continuous ant colony optimization algorithm for the optimal design of storm sewer networks: constrained and unconstrained approach,” Advances in Engineering Software, vol. 41, no. 2, pp. 188–195, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2013 Haixing Wang 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.