Real-Time Vehicle Routing for Repairing Damaged Infrastructures Due to Natural Disasters

We address the task of repairing damaged infrastructures as a series of multidepot vehicle-routing problems with time windows in a time-rolling frame. The network size of the tackled problems changes from time to time, as new disaster nodes will be added to and serviced disaster nodes will be deleted from the current network. In addition, an inaccessible disaster node would become accessible when one of its adjacent disaster nodes has been repaired. By the “take-and-conquer” strategy, the repair sequence of the disaster nodes in the affected area can be suitably scheduled. Thirteen instances were tested with our proposed heuristic, that is, Chen et al.’s approach. For comparison, Hsueh et al.’s approach 2008 with necessary modification was also tested. The results show that Chen et al.’s approach performs slightly better for larger size networks in terms of objective value.


Introduction
Taiwan is located at the active arc-continent collision region between the Luzon arc of the Philippine Sea.Due to the collision of the Eurasian plate and the Philippine plate, Taiwan is at high risk of earthquakes.A 1999 Chi-Chi earthquake with a Richter magnitude of 7.3 1 resulted in 2415 deaths, the collapse of 40,845 buildings, and damage to 44 major roads 2 .In addition, Taiwan also suffers from typhoons each year.Typhoons such as Kalmaegi and Sinlaku in the summer of 2008 and Morakot in the summer of 2009 cause severe fatalities and uncountable damage to agricultural products and transportation infrastructures.
The Federal Emergency Management Agency 3 published Guide for All-Hazard Emergency Operations Planning which is to develop, in partnership with state and local governments, a national emergency management system that is comprehensive, risk based, and all-hazard in approach.Crucial to this system are emergency operations plans EOPs , which describe who will do what, as well as when, with what resources, and by what authority-before, during, and immediately after an emergency.It is apparent that all-hazard emergency operations are interdisciplinary in nature and indeed too broad to be handled by staff in a specialized field.Therefore, to focus, this paper will only address the sequence for repairing damaged infrastructures due to natural disasters.
When a severe natural disaster occurs, transportation infrastructures either areas, roadways, or both in the affected area are often seriously damaged, which makes evacuation and rescue operations difficult.The number of fatalities and property losses can increase if damaged transportation infrastructures are not repaired in a reasonably short period of time.
To accelerate the completion of repairs, a systematic approach is needed.In this research, we address the task of repairing damaged infrastructures, named the real-time vehicle-routing problem with time windows, as a series of vehicle-routing problems with time windows in a time-rolling frame.
The rest of the paper is organized as follows.Section 2 conducts a literature review.Section 3 describes the problem.Section 4 formulates a multidepot vehicle-routing problem with time windows VRPTWs as a mixed integer programming model.Section 5 elaborates the solution algorithm, which is demonstrated with numerical examples in Section 6. Section 7 concludes the paper.

Literature Review
Scheduling the repair sequence of damaged infrastructures is only part of emergency operations plans.The related literature is relatively scarce, and moreover, the approaches adopted to tackle this type of problems are pretty diverse, covering both simulation and optimization techniques, as well as varying from single objective to multiobjective, or even from one level to bilevel optimization.Literature review is, thus, briefly conducted in the following.
Fiedrich et al. 4 suggest a resource allocation decision supporting system for postearthquake search-and-rescue SAR missions.Both simulation and analytical modeling techniques are employed in the system.The system aims to minimize the number of fatalities and injuries and optimally allocate the available resources.To solve the problem, three methods were compared: simulated annealing, tabu search, and hill climbing.Simulated annealing outperformed the other two methods.Chang 5 presents the damaged infrastructure repair problem as a multiple objective vehicle-routing problem with time windows.The objectives were to minimize the total travel time of all construction teams, the total rescue operation time, as well as the risk of construction.Two different approaches were used to solve the problem.The first approach adopts the weighted sum of all three objectives as a single objective, while the second approach uses the "maxmin" operator.
Barbaroso glu et al. 6 develop a mathematical model for helicopter mission planning during a disaster relief operation.The decisions inherent in the problem decompose hierarchically into two subproblems, where tactical decisions are made at the top level, and the operational routing and loading decisions are made at the base level.Consistency beaten the decomposed problems is achieved with an interactive coordination procedure which transfers anticipated information from the base level to improve the top level decisions.The existence of conflicting multiple objectives in this hierarchical structure requires the development of a multicriteria analysis, and an interactive procedure is designed with the top level decision maker to assess the preference of alternative nondominated solutions.Lee 7 uses two models to tackle the damaged infrastructure repair problem.The first model determines the sequence for repairing disaster spots by a fuzzy ranking method and assumes that each disaster spot can only be handled by one construction team.The second model assumes that the travel times are deterministic and allows each disaster spot to be handled by one or more construction teams.Özdamar et al. 8 propose a mathematical model that integrates the multicommodity network flow problem and the vehicle-routing problem for emergency logistics planning in natural disasters.The model is readily decomposed into two multicommodity network flow problems, the first one being linear for conventional commodities and the second integer for vehicle flows .In the solution approach, these submodels are coupled with relaxed arc capacity constraints using Lagrangean relaxation.The convergence of the proposed algorithms is tested on small test instances as well as on an earthquake scenario of realistic size.
To consider both the overall efficiency and total completion time, Liao 9 views the damaged infrastructure repairing problem as a bilevel programming model.The upperlevel subproblem maximizes the overall efficiency of the repair work, whereas the lowerlevel subproblem minimizes the total completion time subject to the time constraint for each disaster spot.
Pettit and Beresford 10 examine the processes involved in emergency relief and highlight the relationships between the participating bodies, that is, military and nonmilitary involvement.A composite model is proposed, which incorporates Jennnings et al.'s 11 model as well as the military/nonmilitary dimension, and it appears to be robust and workable in a range of geopolitical and operation circumstances.
Hsueh et al. 12 optimize the rescue operation using a time-constrained vehiclerouting approach.In their scheme, a fleet of vehicles is dispatched to simultaneously deliver emergency materials and medical equipment to the affected areas and to pick up materials and equipment that are no longer needed.The time-rolling horizon approach is adopted to take care of demands that are either known in advance or that appear in real time.The severity of each disaster node is reflected by the associated time window the tighter the time window, the more severe the disaster node .In addition, it is also assumed that there is no need for the dispatched vehicles to return to the depot.A two-stage solution heuristic is proposed.In the first stage, the minimum cost insertion method is used to find an initial solution; in the second stage, Or-opt method is adopted.
Wang 13 addresses the damaged infrastructure repair problem as a multidepot vehicle-routing problem that adapts fuzzy repair times.A time rolling horizon approach is used to solve the problem.At each specific time instant, only the part of the network consisting of accessible damaged infrastructure segments is treated.The damaged infrastructure segments shown in the partial network are optimally scheduled for repair by the construction teams.In other words, those inaccessible damaged infrastructure segments are not considered for repair at the moment.However, an inaccessible damaged infrastructure segments may become accessible later when its adjacent disaster nodes are repaired either partially or fully.Once a damaged infrastructure segment becomes accessible, it will become visible in the updated network and can be scheduled for repair.At each specific instant, the model minimizes the total travel time.A two-stage solution algorithm is proposed.The first stage uses the nearest/farthest insertion method for generating an initial solution, while the second stage uses the 2-opt, swap, and or-opt methods for improving the initial solution.The disadvantage associated with this approach is that the disaster nodes not appeared in the earlier partial network cannot be considered, and the severity as well as time window constraint of these "invisible" disaster nodes is not suitably reflected in the very beginning of the entire optimization procedure.
Since disaster nodes in the affected area are often arranged in order and the number of construction teams available is usually smaller than the number of disaster nodes, the vehicle-routing approach is appropriate to address the damaged infrastructure repair problem.Therefore, this research treats the damaged infrastructure repair problem as a multidepot vehicle-routing problem with time windows and formulates it as a series of mixed integer programming problems.This is analogous to Hsueh et al.'s 12 approach but with necessary modifications, especially in scheduling the sequence for repairing the damaged infrastructures, which are largely, affected by the changeable link travel times and the completion time of repairing adjacent disaster nodes.

Problem Description
Natural disasters like earthquakes and typhoons often damage infrastructures that need to be repaired.The damaged infrastructure repair problem has the following special features.
1 Repair sequence must be complied between some disaster nodes.One disaster node may be unreachable if it is surrounded by adjacent disaster nodes.In other words, this disaster node can be repaired only when one of adjacent disaster nodes has been repaired.This feature is unique and has to be esteemed in scheduling the repair sequence of disaster nodes.
2 Disaster-related information is changeable.The information such as road travel times, the severity of disaster nodes, or the number of disaster nodes, and construction teams may change from time to time, implying the damaged infrastructure repair plan must be updated in a real-time fashion.
3 Construction teams do not have to return to associated depots when their assigned work has been completed.This assumption is mainly due to the consideration that additional disaster nodes may appear after the work of a construction team has been completed.Note that assigning a construction team on-the-spot to a new disaster node for repair is usually faster than assigning a construction team that is already back to depot.
All these special features need to be properly tackled by the following proposed approach.

VRP Network
To address the damaged infrastructure repair problem as a series of vehicle-routing problems VRPs , a VRP network must be constructed from the corresponding physical transportation network.To show how a physical network can be transformed into a corresponding VRP network, we use a numerical example for illustration.Suppose that we have the physical network of a disaster area as shown in Figure 1 a .One depot is marked in green diamonds; 9 damaged infrastructure segments are marked in gray circles in which the first 8 disaster nodes are identified before performing the construction work and disaster node 9 is reported later , and 15 intersections and/or landmarks are marked in white rectangles.The current link travel times are fixed and externally determined.
To transform a physical transportation network into its corresponding VRP network, the following two steps are performed.Step 1. Omit the intersections marked in white rectangles from the VRP network.
Step 2. Represent headquarters and damaged infrastructure segments as nodes in the VRP network and denote link travel times by the shortest distances between nodes i.e., measured from the physical transportation network for each pair of nodes in Figure 1 a .
The travel time on link 3 → 4 in the VRP network, as shown in Figure 1 b , is represented by the travel time on path 3 → 18 → 4 in Figure 1 a which can be obtained by adding the link travel times for links 3 → 18 and 18 → 4 .When the two nodes are disconnected due to natural disasters, then the corresponding link travel time is represented by a prohibitive large number.
In Figure 1 b , the VRP network is incomplete; that is, some nodes are not reachable and can only be accessed with prohibitive large link travel times.This special feature does not exist in traditional VRP problems, which implies that the solution algorithms for traditional VRP cannot be directly used.For instance, damaged infrastructure 4 cannot be repaired unless damaged infrastructure 3 or 2 has been repaired.This implies that the repair work must be performed in a one-by-one sequential manner.

Repair Sequence
To show how the sequence for repairing damaged infrastructure segments proceeds in natural disasters, we use the hypothetical network shown in Figure 1 b for illustration.At the beginning, this network consists of one depot numbered as "0" and eight disaster nodes numbered from "1" to "8" .A depot is where construction teams initially reside; disaster nodes represent damaged infrastructures to be repaired.Two types of links are used.A solid line denotes a normal link, that is, the corresponding link travel time is finite, and its end nodes are connected.A dotted line indicates that the corresponding link travel time is prohibitively high, and its end nodes are adjacent but cannot be accessed by passing through it, that is, an isolated link.The accessibility of one disaster node by the construction team is clearly dependent on whether there are one or more solid links associated with it.
With this updated travel time information, as shown in Figure 2 b , node 2 now acts as the new starting node for the dispatched vehicle in the network.Referring to Chen et al. 14 , this new starting node is called the critical node.The critical node plays an important role in distinguishing between a serviced disaster node and an unserviced disaster node.A critical node is defined as a node that is currently using a construction team or to where a construction team is heading.Critical nodes need to be identified instantly when a demand or travel time has changed in real time so that the route can be reconstructed.
Upon repairing node 2, nodes 1, 3, 4, 5, and 8 become accessible.We assume by some criterion that the construction team selects node 1 among five accessible nodes for repair.As shown in Figure 2 c , upon repairing node 1, nodes 3, 4, 5, and 8 become accessible.Of the four accessible nodes, the construction team next chooses node 3 for repair.Figure 2 d indicates that nodes 4, 5, and 8 become accessible after repairing node 3. Figure 2 e shows nodes 5, 7, and 8 become accessible when the construction work for node 4 is completed.In addition, disaster node 9 is identified now.Figures 2 f , 2 g , 2 h , 2 i , and 2 j indicate that nodes 5, 9, 7, 6, and 9 are repaired sequentially, and the entire construction work has been completed.Figure 2 k shows the entire sequence for repairing all disaster nodes.
From the above description, it is clear that a disaster node that is initially inaccessible by the construction team can become accessible if its adjacent disaster nodes are repaired.This description naturally motivates a "take-and-conquer" solution concept for the damaged infrastructure repair problem, which is essentially equivalent to solving a series of VRP problems in a time-rolling horizon.Each VRP network to be tackled at a specific instant is constructed based on updated information such as link travel times and real-time repairing requests from the disaster area.The detailed steps of the proposed solution algorithm will be discussed later.

Assumptions
The damaged infrastructure repair problem is very complex, and there is no common definition in the literature.For the purpose of our research, we make the following assumptions.must be restored by the emergency construction teams.In the second stage, full functioning of the infrastructure segments must be achieved.In this research, only the first stage is addressed.

Notations
The following notations are used in this paper.
2 Parameters and constants Δ: CPU time estimated for solving the problem, τ: time instant equal to the time for receiving new information plus the CPU time Δ required for solving the problem.Note the CPU time can be neglected without causing significant effect, given the relatively short solution times and long travel and repair times, c ij τ : estimated travel time on link i → j at time instant τ, Mathematical Problems in Engineering r i τ : estimated repair time for disaster node i at time instant τ, l i : latest time or upper bound for node i to be repaired, k i τ : specific construction team that leaves critical node i at time instant τ, W i τ : estimated weight associated with the violation of time window, reflecting the importance of the time window associated with disaster node i at time instant τ.
3 Node sets N 0 τ : set of depots at time instant τ, N c τ : set of critical nodes at time instant τ, N u τ : set of unserviced disaster nodes at time instant τ, N cu τ : set of critical and unserviced disaster nodes at time instant τ, N 0u τ : set of depots and unserviced disaster nodes at time instant τ, N 0cu τ : set of depots, critical and unserviced disaster nodes at time instant τ.
4 Sets of construction teams K i τ : set of construction teams available at depot i, i ∈ N 0 τ , at time instant τ; K i τ ⊆ K τ , K τ : set of the total construction teams available in the entire affected area at time instant τ.
5 Variables a i : arrival time at disaster node i, d i : departure time from disaster node i, d ik : departure time of vehicle k from depot i, i ∈ N 0 τ , x ijk : 0-1 decision variable; if construction team k leaves disaster node i toward disaster node j, then x ijk 1; otherwise, x ijk 0.

Model Formulation
The damaged infrastructure repair problem is addressed as a series of multidepot VRP problems in the time-rolling horizon.Whenever network information is updated, a new VRP problem must be formulated as a 0-1 mixed integer programming problem and solved accordingly.Referring to Hsueh et al. 2008 , the multidepot VRP problem at specific instant τ can be formulated as follows.

4.10
Objective function 4.1 minimizes the weighted sum of the total vehicle travel time and the penalty associated with violating the time window at all disaster nodes.The penalty incurred at node i is defined by the multiplication of the violation of time window and the estimated weight W i τ .If the lower limit of the time window at disaster node i is violated, that is, a i − l i > 0, then a penalty occurs; otherwise, the penalty vanishes.Note that the larger the value of weight W i τ , the higher the priority given to disaster node i.When all the assigned weights are very large, the final solution will be largely determined by the "soft" time-window constraints, rather than the total travel time, in the objective function.
Constraints 4.2 through 4.7 manage construction teams at disaster nodes.Equation 4.2 requires that only one construction team can leave from a critical disaster node or unserviced disaster node i once.Equation 4.3 denotes that only one construction team can arrive at unserviced disaster node j once.Constraints 4.2 and 4.3 together indicate that each disaster node can be serviced by exactly one construction team.Constraint 4.4 requires that one construction team that enters disaster node h earlier must leave that node later.Equations 4.5 -4.6 designate that each vehicle can leave the depot once at most.Constraint 4.5 requires that construction team k i which has arrived at or is approaching critical disaster node i, must also leave this disaster node once.Note that construction team k i is known at time τ.Constraint 4.6 denotes that each construction team k can be dispatched from depots at most once.Constraint 4.7 defines variable x ijk as a 0-1 integer variable.If construction team k leaves node i toward j, then x ijk 1; otherwise, x ijk 0. Constraints 4.8 through 4.10 define arrival and departure times.Constraint 4.8 defines the arrival time at node j by adding the departure time from node i to the link travel time between node i and node j.Constraint 4.9 defines the completion time of repair work at node i.If node i is an unserviced node, then the completion time of the repair work is defined by the arrival time of the construction team plus the service time; that is, d i a i s i .If node i is a critical node, there are two conditions to be considered.First, if disaster node i is not the last to be serviced, that is, j∈N 0 k∈K x ijk 0, then the completion time for the repair work is defined by formula d i max τ, a i s i .On the other hand, if node i is the last to be serviced, that is, j∈N 0 k∈K x ijk 1, then the departure time of the construction team is set as infinitive; that is, d i ∞.This is equivalent to requiring the construction team stay at the current node to wait for the next assignment after the repair work is finished.Constraint 4.10 defines the departure time for construction team k to leave the depot.

Solution Procedure
As mentioned, the damaged infrastructure repair problem is addressed as a series of multidepot VRP problems with time windows.The solution procedure continuously updates real-time information such as demands and link travel times and determines when a new VRP problem is to be formed.Each VRP problem is then solved using a two-stage solution algorithm in which an initial solution is generated in the first stage and further improved upon in the second stage.In the following sections, the framework of the solution procedure is introduced, followed by the initial solution construction and improvement.

Unified Framework
While repairing the damaged infrastructures, relevant information is continuously updated.To take this new information into consideration, the routing schedule of the construction teams must be reoptimized.Changeable information could include the link travel time, the number of construction teams, as well as the number of disaster nodes.
1 New link travel times: suppose that we have two nearby damaged bridges, denoted by two disaster nodes, along a freeway.If the nearest damaged bridge has not been repaired, then the second bridge cannot be easily reached.In this situation, the link travel time between them is very large and is symbolically represented by the total network travel time under normal conditions.However, when the near damaged bridge is repaired, the farther damaged bridge becomes accessible, and repair work can be scheduled.In this situation, the link travel time is reduced from a large number to a much smaller value and hence must be updated.This operational sequence is called a take-and-conquer strategy.In other words, whenever new information is generated the network data, and hence, the optimal schedule for repairing the damaged infrastructure segments will be changed accordingly.
2 new construction teams: new construction teams may be formed or borrowed/ leased from other administrative agencies.Since different number of construction teams may have different routing plans, the corresponding VRP problem needs to be reformulated and reoptimized.3 New disaster nodes: information about infrastructure damages is not perfect immediately following a natural disaster but will become more transparent over time.When new disaster nodes are taken into consideration, the repair plan must be reoptimized.
The framework for the entire solution procedure is depicted in Figure 3.The repair operation begins with a routing plan based on the available information.This routing plan can be changed, if necessary, to include updated information such as new link travel times, new disaster nodes, as well as new construction teams.The new routing plan can be obtained by solving the VRP problem, which must start routes from the critical nodes.A critical node may be thought of as a virtual depot for the subnetwork consisting of the remaining unserviced disaster nodes.This procedure for solving a series of multidepot VRP problems with time windows continues until no more unserviced disaster nodes exist.As in the literature, the time-constraint VRP problem is solved using a two-stage algorithm, that is, initial solution generation and initial solution improvement, which will be detailed in the following sections.

Initial Solution Generation
To generate initial routes, the minimum cost insertion method is adopted.A route starts from the depot and chooses the next disaster node with the minimum insertion cost.
The incurred insertion cost for unserviced disaster node k between two consecutive serviced disaster nodes i and j, in a "construction team" route, is computed as follows.Insertion cost cost of travel from disaster node i to the newly scheduled disaster node k cost of travel from the newly scheduled disaster node k to disaster node j − cost of travel between disaster node i and disaster node j on the "construction team" route .This process for choosing the next disaster node continues until all disaster nodes have been serviced see Figure 4 .The minimum cost insertion method can also be applied with minor modifications to situations at which multiple depots exist.

Mathematical Problems in Engineering
The steps for generating initial routes can be described as follows.
Step 1. Input data.Input link travel times, estimated repair time for each disaster node, and upper limit of time windows for each disaster node.
Identify critical disaster nodes to begin with initial route construction for the remaining unserviced disaster nodes.
Step 3. Construct initial VRP routes using the minimum insertion cost method.Substep 3.1.Select an unserviced disaster node one by one for each VRP route and compute its insertion cost including its added travel time and penalty incurred due to the violation of time windows at different positions in the route.Substep 3.2.Choose the route position with the minimum insertion cost.
Step 4. Insert the chosen disaster node into the best position of the route and update the link travel times for the subsequent disaster nodes accordingly.
Step 5.If no more unserviced nodes exist and the initial routes have been constructed, then return to the unified framework; otherwise, go back to Step 3.

Initial Solution Improvement
To improve the initial solution, two types of heuristics can be considered.The first type is known as local search in which nodes and/or links are exchanged within intraroutes or between interroutes.It is easy to implement; however, in some circumstances, the heuristic may be trapped in a local region.In addition, the repair sequence or precedent-successor relationship of disaster nodes makes the local search complicated.The second type, called metaheuristics, is capable of moving beyond local optima and hence suitable for difficult combinatorial problems such as vehicle-routing problems.A well-known metaheuristics is called ant colony optimization ACO which involves a collection of various ant-related algorithms.The ant system 15 was the first ACO made available to academia.Later, various modified algorithms were published, such as the elitist ant system 16 , the ant colony system ACS 17 , the max-min ant system 18 , the rank-based ant system 19 , the fast ant system 20 , and approximate nondeterministic tree-search procedures 21 .
In this paper, we adopt the ACS, to improve the initial solution, because it has been successfully employed in solving VRP problems.As an ant-related algorithm, the ACS is motivated by the behavior of ants seeking a path between their colony and a source of food.When one ant finds a good i.e., short path from the colony to a food source, other ants are more likely to follow that path, and positive feedback eventually leads all the ants following a single path.The idea of the ACS is to mimic this behavior with "simulated ants" walking around the graph representing the problem to solve.
The ACS consists of four major steps: initial condition, edge selection, pheromone update, and termination criterion.We describe them below. (

1) Initial Condition
Six parameters are used to set the initial values Table 1 . (

2) Link Selection Rule
To select the next link for a route to augment, the state transition rule is employed.Two rules, exploitation and exploration, can be used depending on the random number q generated.Constant between 0 and 1, 0 ≤ q 0 ≤ 1; we set q 0 0.9 If q ≤ q 0 , then the link with the higher pheromone value and the shorter length will be chosen as shown in 5.1 .

5.1
On the contrary, if q > q 0 , then an ant will move from node i to node j with the probability computed by the following formula: where, τ ij is the amount of pheromone on link i, j , η ij is the desirability or visibility of link i, j a priori knowledge, typically η ij 1/L ij , L ij is the length of link i, j , J k i is the set of unserviced nodes indexed as i by ant k, β represents the importance of link length relative to pheromone value β > 0 , q is a random variable uniformly distributed between 0 and 1, S is the next node as determined by 5.1 , S k is the next node to visit as determined by 5.2 .
The probability for link i, j being chosen will increase as the associated pheromone value gets higher and distance gets shorter.However, it is no guarantee that the link with the highest probability will be chosen for route expansion.
(3) Pheromone Update Two rules are used to update pheromone values.The local updating rule shown in 5.3 is used to update pheromones on all routes when all ants have moved to next node.
As shown in 5.4 , the global updating rule is used to increase the pheromone value on the shortest path only when all ants have completed their routes where ρ is the rate of pheromone evaporation in the local updating rule, 0 < ρ < 1, α is the rate of pheromone evaporation in the global updating rule, 0 < α < 1, Δτ ij is the amount of pheromone deposited, typically given by where Q is constant and L gb is the length of the shortest path.

(4) Termination Criterion
Four termination criteria are commonly used in the ACS method: the maximum number of iterations, the maximum number of consecutive iterations in which the objective value is not improved, the maximum execution time, or the acceptable objective value.In this paper, a maximum number of iterations after which the corresponding objective value cannot be improved are used as the termination criterion.The flow chart for the initial solution improvement is depicted in Figure 5.
According to Figure 5, the algorithmic steps of the ACS can be described as follows.
Step 1. Set parameters.Initialize the number of antsm, the rate of pheromone evaporation ρ in the local updating rule, the rate of pheromone evaporation α in the global updating rule, the parameter to represent the relative importance of link length to pheromone value β, and the initial pheromone value τ 0 .Set the initial solution as the temporary best solution.
Step 2. Dispatch all ants from the depot s .
Step 3.For each ant, select the most suitable construction team and send it to the next disaster node using the link selection rule.
Step 4. For each ant, perform the local rule to update pheromone values.
Step 5.If all ants have completed route construction, continue; otherwise, go back to Step 3.
Step 6. Determine if any constructed solution is superior to the current best solution.If not, make no change; otherwise, replace the current best solution and perform the global rule to update pheromone values.
Step 7. Determine if the termination criterion is satisfied.If so, then stop; the stage of initial solution improvement has completed.Otherwise, go back to Step 2.

Computational Experiments
We generate 13 numerical examples for testing.Our assumptions are described as follows.
1 There are two depots, indexed as 1 and 2, respectively.
2 For example 1, 5 construction teams are available by the time when disaster occurs.Among which 2 construction teams are located in depot 1 while the other 3 construction teams in depot 2.
3 For examples 2-13, 15 construction teams are available by the time when disaster occurs.Among which 7 construction teams are located in depot 1 while the other 8 construction teams in depot 2.
4 There are between 10 and 70 nodes, indexed from 3 on.
5 The recovered travel times for all pairs of two disaster nodes are randomly generated between 0.3 and 28 hours; hence, all test problems are complete networks after repairing.
6 The repair time for each disaster node is randomly generated from 2 to 90 hours.
7 The time window for each disaster node is externally determined.
8 The unit penalty for the violation of gold rescue time i.e., 72 hours is set as 10.
Using the above data two experiments are conducted.The first experiment adopts our model and solution algorithm called Chen et al.'s approach hereafter , while the second experiment chooses an outside source, that is, Hsueh et al.'s approach 2008 with necessary modifications.Hsueh et al.'s approach is different from ours in four aspects: 1 it is designed for single depot problem rather than multiple-depot problem; 2 it takes care of both delivery and pickup operations; 3 it assumes network is always connected and will not change throughout the analysis period; 4 there is no restriction on the sequence of visiting customer nodes.The obtained results from the test networks are summarized in Table 2.
Both Chen et al.'s and Hsueh et al.'s approaches can improve 12 of 13 test examples in terms of objective function value, with average improvement rates of 11.13% and 11.37%, respectively.As can be seen in Figure 6, Hsueh et al.'s approach seems to perform slightly better for smaller size networks in terms of objective values.The detailed information for the best results is included as Table 3.Note that for smaller size networks, each construction team can only service very few disaster nodes, usually 2∼3 disaster nodes; no apparent advantage can be concluded in favor of either approach.However, Chen et al.'s approach tends to outperform for larger size networks in terms of objective value.
The advantage of Chen et al.'s approach over Hsueh et al.'s approach on Problem 12 can be further justified by varying number of construction teams.As shown in Figure 7, Chen et al.'s approach is better than Hsueh et al.'s approach in terms of objective value for the number of construction teams ranging from 4 to 15.In a separate experiment see Table 4 , we also observed that Chen et al.'s approach is better than Hsueh et al.'s approach in terms of objective value for the majority of test problems when either depot 1 or depot 2 is closed.
In Table 2, the average computational CPU times for the two approaches are 20.64 and 5.94 seconds, respectively, implying that the Chen et al.'s algorithm is computationally intensive compared to the Hsueh et al.'s method.As shown in Figure 8, both approaches indicate that the larger the network size, the higher the CPU time.However, even with the largest test network, the problem can still be solved within 3 minutes using any method.These low CPU times are not a large concern for this type of problem.Based on the above discussions, it is concluded that both Chen et al.'s and Hsueh et al.'s approaches are suitable for scheduling rescue operations for smaller size networks.However, for larger size networks, Chen et al.'s approach is preferred due to lower objective values.In all cases, CPU time is not a big concern.

Conclusion
This research addresses the damaged infrastructure repair problem as a series of multidepot VRP problem with time windows and presents 13 hypothetical test examples.In addition to link travel times, other real-time information such as new disaster nodes as well as added construction teams can also be considered in the unified framework.The time-rolling horizon approach is appropriate for dealing with the damaged infrastructure repair problem, especially when some disaster nodes have precedentsubsequent sequential relationships.As precedent disaster nodes are repaired, subsequent disaster nodes become reachable and can be repaired.
To deal with this type of problem, we propose a heuristic, that is, Chen et al.'s, approach and include an outside source, that is, Hsueh et al.'s 2008 approach with necessary modifications for comparison.The former embeds ACS which can efficiently improve the initial solution generated by the minimum cost insertion method, with an average improvement rate of 11.13% on 13 numerical examples.
The experiments show that both Chen et al.'s and Hsueh et al.'s approaches are suitable for scheduling rescue operations.However, Chen et al.'s approach is slightly preferred due to lower objective values.In all cases, CPU time is not a big concern.

Suggestions
Some important issues are listed as follows.
1 Because each construction team is equipped with different machines and tools that may not meet the requirements of all disaster nodes, a heterogeneous VRP problem may be more appropriate.
3 If spatial information can be displayed on a computer screen, then the decision maker can make more reasonable decisions about the routing schedule of construction teams.For real-world applications, a screen-based decision support system using the best proposed method along with the geographic information system GIS is worth further exploring.This is because the geographical characteristics of disaster nodes as well as the damage conditions associated with these disaster nodes can be elaborated.

Figure 1 :
Figure 1: a Physical transportation network for the disaster area.b VRP network for the disaster area.

Figure 3 :
Figure 3: Unified framework for solution procedure.

Figure 4 :
Figure 4: Flow chart for initial solution construction.

. 1 α. 9 β 5 m
Rate of pheromone evaporation in local updating rule, 0 < ρ < 1; we set ρ 0Rate of pheromone evaporation in global updating rule, 0 < α < 1; we set α 0Parameter to control the importance of edge desirability, determined by the trial-and-error method, usually greater than 0; we set β Number of artificial ants used; we set m 10.Each ant simulates a solution consisting of a set of routes τ 0 Initial pheromone value, usually set as τ 0 NL nn −1 , where N denotes the number of links and L nn represents the route length computed by the nearest neighbor search Iter the maximum number of ACS iterations; we set Iter 30

Figure 5 :
Figure 5: Flow chart for route improvement of emergency vehicle Routing.

Table 2 :
The results for the test networks.

Table 3 :
The best results for the test networks.Problem no.Nodes Algorithm Total travel time Violation of time windows Objective Teams no.
Remark: CPU time is recorded for Microsoft Windows XP, INTEL Pentium 4 3.0 GHz, 1 GB RAM.

Table 4 :
The results for the test networks with a single open depot.: CPU time is recorded for Microsoft Windows XP, INTEL Pentium 4 3.0 GHz, 1 GB RAM. Remark