Next Article in Journal
Optimized Dissolved Oxygen Fuzzy Control for Recombinant Escherichia coli Cultivations
Next Article in Special Issue
Metaheuristics in the Humanitarian Supply Chain
Previous Article in Journal / Special Issue
Feature Selection for High-Dimensional Datasets through a Novel Artificial Bee Colony Framework
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Travel Time Reliability-Based Rescue Resource Scheduling for Accidents Concerning Transport of Dangerous Goods by Rail

School of Traffic and Transportation Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
*
Author to whom correspondence should be addressed.
Algorithms 2021, 14(11), 325; https://doi.org/10.3390/a14110325
Submission received: 13 October 2021 / Revised: 1 November 2021 / Accepted: 2 November 2021 / Published: 5 November 2021
(This article belongs to the Special Issue Metaheuristic Algorithms in Optimization and Applications 2021)

Abstract

:
The characteristics of railway dangerous goods accidents are very complex. The rescue of railway dangerous goods accidents should consider the timeliness of rescue, the uncertainty of traffic environment and the diversity of rescue resources. Thus, the purpose of this paper is to confront the rescue resources scheduling problem of railway dangerous goods accident by considering factors such as rescue capacity, rescue demand and response time. Based on the analysis of travel time and reliability for rescue route, a multi-objective scheduling model of rescue resources based on travel time reliability is constructed in order to minimize the total arrival time of rescue resources and to maximize total reliability. The proposed model is more reliable than the traditional model due to the consideration of travel time reliability of rescue routes. Moreover, a two-stage algorithm is designed to solve this problem. A multi-path algorithm with bound constraints is used to obtain the set of feasible rescue routes in the first stage, and the NSGA-II algorithm is used to determine the scheduling of rescue resources for each rescue center. Finally, the two-stage algorithm is tested on a regional road network, and the results show that the designed two-stage algorithm is valid for solving the rescue resource scheduling problem of dangerous goods accidents and is able to obtain the rescue resource scheduling scheme in a short period of time.

1. Introduction

Dangerous goods include the substances and articles the carriage of which is prohibited by the Regulation Concerning the International Carriage of Dangerous Goods by Rail (RID) or authorized only under the conditions prescribed therein [1]. Specific regulations are in place in relation to the transport and storage of dangerous goods, as well as employee, consumer and environmental protection, in order to prevent accidents during carriage of such goods.
Transport of dangerous goods is regulated in order to prevent, as much as possible, accidents related to persons or property and damage to the environment, the means of transport employed or to other goods. The transportation of dangerous goods (TDG) includes activities related to the movement of dangerous goods from their place of manufacture or storage to their destination with the preparation of cargo, packaging, vehicles and crew, reception of goods, carrying out cargo operations and short-term storage of goods at all stages of their transfer [2].
Currently, a large quantity of different types of dangerous goods is transported by road, rail, inland waterways and maritime, pipelines and air. Among them, railway transportation undertakes the main traffic volume. The influence of random factors and events can result in an accident resulting in leakage, combustion or explosion of hazardous substances during loading and unloading, transportation and storage procedures. These types of incidents not only threaten the safety of rail transport but also the life, environment and property [3]. The most common places where dangerous goods accidents happen are dangerous goods handling locations of petrochemical and logistics enterprises, railway stations and lines along the railway. Moreover, these locations have many internal equipment and facilities, and their external rescue environment is very complex. Therefore, once an accident occurs, it is necessary to carry out emergency rescue in a timely manner; otherwise, it will cause serious consequences.
There has been a substantial amount of investigations on emergency resource scheduling or its related problems. Emergency resource scheduling is the scheduling optimization operation spanning across the categories of routing and emergency materials distribution.
In emergency relief research, reviews by Caunhye et al. [4], Galindo and Batta [5] and Zhou et al. [6] revealed that most of the models developed for emergency relief only permit scheduling of one type of resource, which is either expendable or non-expendable [7]. Özdamar et al. [8], Chang et al. [9], Balcik et al. [10] and Huang et al. [11] only considered resource scheduling related to expendable resources. In contrast, Lassiter et al. [12], Wex et al. [13] and Schryen et al. [14] only considered resource scheduling related to non-expendable resources in disaster emergencies.
Some researchers have also proposed models for multi-resource scheduling of emergency relief. Zheng and Gao [15] built a multi-resource emergency systems model with consumption rates of nonnegative and integrable functions based on the object of earliest emergency-start-time, and the corresponding algorithm was provided. Lee et al. [16] studied a scheduling problem with operations that required renewable as well as non-renewable resources and proposed a framework of heuristic procedures for solving this problem. Zhang et al. [17] built an emergency resource scheduling model, which included multiple suppliers with a variety of resources, a single accident site and some restrictions, and an adaptively mutate genetic algorithm is applied to figure out a superior solution. Shahparvari and Bodaghi [18] developed a mixed integer programming model to support tactical decision making in allocating emergency relief resources in the context of the Black Saturday bushfires. Furthermore, a possibilistic programming approach was employed to minimize the transportation disruption risk. In addition, Bodaghi et al. [19] analyzed several probabilistic scenarios to determine the most frequent emergency operation plan and the most persistent best compromised emergency operation plan. Li [20] constructed a resource scheduling model by considering influencing factors such as transportation time, cost and resource availability. Particle swarm optimization combined with Cuckoo search (PSO-CS) algorithm was applied to solve a resource scheduling process of multiple types and multiple classes. Bodaghi et al. [7] formulated the Multi-Resource Scheduling and Routing Problem (MRSRP) for emergency relief and developed a solution framework to effectively deliver expendable and non-expendable resources in emergency recovery operations. Tang and Sun [21] established two optimal models, respectively; the single objective model aimed at minimizing the time of emergency resource dispatch, and the multi objective model was for the sake of minimizing the time of emergency resource dispatch and the number of emergency rescue bases, and a voting analytic hierarchy process was used to solve the model. Wu and Xiang [22] constructed a capacitated vehicle routing problem model suitable for emergency resource scheduling process and designed a simulated annealing-genetic algorithm (SA-GA).
Many of the models developed for emergency resources scheduling are deterministic and had the minimization of travel time as one of the objectives; they ignored the rescue routing problem. Moreover, heuristic methods are always used to find a near-optimal due to the NP-hard nature of the problem. For example, the PSO-CS algorithm was applied to solve the resource scheduling process of multiple types and multiple classes by Li [20]. An adaptively mutate genetic algorithm (GA) was used to solve an emergency resource scheduling model and included multiple suppliers with a variety of resources by Zhang et al. [17]. SA-GA was proposed in order to confront the capacitated vehicle routing problem for emergency resource scheduling process by Wu and Xiang [22]. Wu et al. [23] used GA to find out the dynamic emergency distribution path. Zidi et al. [24] and Mguis et al. [25] proposed a multi-agents approach by using a genetic algorithm for scheduling vehicle routing and local search for the management of an eventual event. Moreover, Liu et al. [26] designed a quick heuristic algorithm in order to obtain a fleet dispatching plan. Ceselli et al. [27] proposed an algorithm based on the branch-and-cut-and-price paradigm in order to deal with the generalized location and distribution problem. Tang et al. [28] used the greedy algorithm to optimize the distribution path. Zhao et al. [29] proposed a two-stage shortest path algorithm, which is composed of K-paths algorithm and shuffled frog leaping algorithm, in order to solve the dynamic paths planning problem of emergency vehicles. Özdamar and Demir [30] designed a multi-level clustering algorithm for coordinating vehicle routing in large-scale post-disaster distribution and evacuation activities. He et al. [31] adopted the K-means clustering algorithm in order to obtain some of the local distribution centers and used the PSO algorithm to design the local optimal allocation routings of emergency relief vehicles. Vargas-Florez et al. [32] used the hierarchical ascending clustering approach in order to discover the reliable delivery routes. Gharib et al. [33] proposed the NSGA-II and multi-objective firefly algorithm in order to deal with the multiobjective vehicle routing model that was developed. Penna et al. [34] proposed a generic hybrid heuristic framework to solve the vehicle routing problem after a natural disaster.
All these related investigations have given us much inspiration when studying the rescue resource scheduling problem. For dangerous goods accidents of railways, the rescue resource scheduling should meet the following requirements:
  • The rescue resource scheduling problem is one of time urgency, and the rescue resources should arrive at the accident location in time. Since time is critical in railway dangerous goods accidents and short response time can help save more properties and lives and reduce the impact of accidents, the rescue work should be started at the first time of incidence. Therefore, the earlier the rescue resources are transported to the accident location, the better the rescue effect.
  • Railway dangerous goods accident has continuous expansibility, and emergency rescue has time window requirements. Expansion of accident is the main feature of railway dangerous goods accidents, and it very easily affects other equipment or facilities and causes accident expansion. The characteristics of accident expansion require that rescue resources must be transported to the accident location within the time window before accident expansion.
  • The traffic environment is uncertain, and rescue time reliability is an important quality performance measure of the rescue resource scheduling problem. Road traffic conditions are more complex and uncertain, which requires that the arrival time of rescue vehicles must be reliable.
  • Railway dangerous goods accidents require coordinated dispatching of multi-rescue resources. The characteristics of railway dangerous goods determine the multiplicity of accident consequences. Therefore, accident rescues of railway dangerous goods generally require multi-rescue resources of fire control, medical and health.
Therefore, in the case of a railway dangerous goods accident, transportation of rescue equipment (such as fire truck, fire sprinkler, fire water monitor, emergency lighting equipment, chemical rescue vehicle, special rescue vehicle, medical rescue vehicle, etc.) and resources (such as rescue medicine, air respirator, etc.) to the accident location is urgently required. At the same time, due to the influence of rescue capacity, location, road network, traffic environment and traffic flow of each rescue center, the appropriate rescue resource scheduling and routing plans need to be integrated and reasonably arranged.
The remainder of the paper is organized as follows: Some mathematical notations are defined in Section 2.1 and the travel time and reliability of a rescue route are analyzed in Section 2.2. The constraints and objective function are discussed, and then a multi-objective scheduling model of rescue resources is established in Section 2.3. A multi-path algorithm with a bound constraint as the first stage of the two-stage algorithm is proposed in Section 3.1, and the NSGA-II algorithm as the second step is designed in Section 3.2. The proposed algorithm is tested on a regional road network in Section 4. Finally, the conclusions are provided in Section 5.

2. The Reliability Based Rescue Resource Scheduling Model

In this section, a multi-objective optimization model of rescue resources is developed by considering the travel time and reliability of the rescue route. The travel time of links and intersection delay, the travel time and reliability of the rescue route are analyzed based on the mathematical notations defined in Section 2.1. Moreover, the constraints of resource demand, rescue capacity and flow balance, etc., are analyzed, and two objective functions are formulated in Section 2.3. The flow chart of this section is presented in Figure 1.

2.1. Notations

Variables and symbols in the model are defined as follows:

2.1.1. Sets and Indices

The sets and indices in this paper are listed in Table 1.

2.1.2. Parameters

The parameters of rescue demand and rescue center are listed in Table 2.

2.1.3. Decision Variables

The decision variables are listed in Table 3.

2.2. Analysis of Travel Time and Reliability

Due to the influence of various random events, such as weather conditions, traffic jams and traffic incidents, the supply of a traffic network may be uncertain, which varies the travel time of a rescue route randomly. The different positions of different rescue centers result in the differences in travel time and reliability of rescue routes. Therefore, it is necessary to analyze the travel time of a rescue route in rescue resource scheduling.

2.2.1. Link Travel Time and Intersection Delay

For a road section, its travel time is mainly affected by link traffic flow and intersection delay. The travel time t m n of rescue vehicles on a link mainly includes two parts: the travel time of vehicles passing through link e m n under free-flow conditions, t m n 0 , and the delay time on the link e m n due to vehicle speed fluctuation, ξ m n . Thus, the travel time of a link can be expressed as the following:
t m n = t m n 0 + ξ m n
where delay time ξ m n is a random variable that is related to the average speed of vehicles, vehicle density and time headway.
Rescue vehicles have priority when passing through the intersection, but they may face delay due to the influence of social vehicles. Intersection delay time d v can be divided into two parts—the dissipation time of vehicles driving in front of rescue vehicles and the time of rescue vehicles passing through the intersection [35]:
d v = d v 0 + ψ v
where d v 0 refers to the time when rescue vehicles pass through the intersection without delay, and ψ v represents the dissipation time of other vehicles driving in front of rescue vehicles, which is related to dissipation speed, queuing length and driving distance of the queuing vehicles. Since the arrival rate of traffic flow at the intersection is a random distribution, ψ v can be considered as a random variable.

2.2.2. Travel Time of a Rescue Route

A rescue route shall be a simple route, and its travel time can be expressed by the travel time of each link and intersection included in the route. The travel time of a rescue route k can be formulated as follows.
t 2 k r s = t 1 i r + e m n E t m n δ m n k r s + v V d v δ v k r s
Assuming that the travel time of each link and intersection is independent of each other, the expected value and standard deviation of the travel time for route k can be expressed as follows:
E ( t 2 k r s ) = e m n E t m n 0 δ m n k r s + v V d v 0 δ m n k r s + E ( t 1 i r ) + e m n E E ( ξ m n ) δ m n k r s + v V E ( ψ v ) δ v k r s
D ( t 2 k r s ) = D ( t 1 i r ) + e m n E D ( ξ m n ) δ m n k r s + v V D ( ψ v ) δ m n k r s
where E ( ) and D ( ) are the expected value and standard deviation of random variables, respectively.
According to the central limit theorem, the travel time of a rescue route can be considered to obey the normal distribution with the expected value of E ( t 2 k r s ) and standard deviation of D ( t 2 k r s ) [36,37].
t 2 k r s N ( E ( t 2 k r s ) , D ( t 2 k r s ) 2 )
Meanwhile, the probability density function for the travel time of a rescue route is as follows.
Φ ( t ) = 1 2 π D ( t 2 k r s ) t exp ( ( t E ( t 2 k r s ) ) 2 2 D ( t 2 k r s ) 2 ) d t
5.
Route travel-time reliability
According to the travel characteristics of dangerous goods accident rescue, the travel time reliability of a rescue route is the probability that the rescue vehicle can reach the accident location within the maximum tolerable time, which can be expressed as follows.
R k r s = P { t 2 k r s R T i } = Φ ( R T i )
Generally, the more routes that can be selected between the rescue center and the accident location, the greater the possibility that the rescue vehicle reaches the accident location within the maximum tolerable time, and the higher the reliability. If there are multiple routes between the rescue center and the accident location, it can be considered as a parallel system, and the multi-route travel time reliability can be formulated as follows [38].
R r s = 1 k r s K r s ( 1 R k r s )

2.3. The Optimization Model

  • Analysis of constraints
The sum of emergency materials, equipment or personnel provided by each rescue center shall meet the demands of the accident location.
r R E M i r E D i
In addition, the number of emergency materials, equipment or personnel provided by each rescue center shall be less than or equal to its rescue capacity.
E M i r E N i r
The rescue center r can reach the accident location within the maximum tolerable time by selecting the rescue route k r s with confidence level α .
R k r s = P { t 2 k r s R T i } α
Constraint (13) ensures that route k r s must not have a loop.
v V δ v k r s 1
The flow balance constraint for each node in a transportation model network is as follows.
n V : e m n E δ m n k r s n V : e n m E δ n m k r s = { 1 , m = s 1 , m = r 0 , otherwise
2.
The objective function
For accident rescue, the arrival time of emergency materials, equipment and personnel and the travel-time reliability of the rescue route are the main concerns. Once a dangerous goods accident occurs, all rescue centers will respond at once. Therefore, the arrival time can be expressed by the sum of the expected value of the route travel time of emergency materials, equipment and personnel of each rescue center.
Z 1 = i I r R [ E M i r E ( t 2 k r s ) ]
In terms of reliability, the total reliability of various emergency materials, equipment and personnel arriving at the accident location is used to represent the reliability of the rescue centers relative to the accident location, which can be expressed as follows.
Z 2 = i I r R ( E M i r R k r s )
3.
The multi-objective rescue resource scheduling model
According to the above analysis, by considering the factors such as the rescue capacity, rescue demand and response time, the multi-objective scheduling model of rescue resources based on travel time reliability is constructed in order to minimize the total arrival time of rescue resources and to maximize total reliability.
min Z 1 = i I r R [ E M i r E ( t 2 k r s ) ] , max Z 2 = i I r R ( E M i r R k r s ) , s.t. r R E M i r E D i , i I , E M i r E N i r , i I , r R , R k r s = P { t 2 k r s R T i } α , i I , r R , k r s K r s , v V δ v k r s 1 , k r s K r s , n V : e m n E δ m n k r s n V : e n m E δ n m k r s = { 1 , m = s 1 , m = r 0 , otherwise , k r s K r s , E M i r N , δ m n k r s , δ v k r s 0 , 1 , i I , r R , k r s K r s

3. Solution Algorithms

The rescue resource scheduling problem based on travel time reliability needs to confront two problems: one is the vehicle routing problem from each rescue center to the location of accident under the requirements of reliability, and the other is the quantity allocation of rescue resources provided by each rescue center. In order to reduce the impact caused by the accident as soon as possible, the rescue center with the shortest travel time and highest reliability should be preferred in terms of dispatching rescue resources.
The time used for solving a rescue resource scheduling problem in a realistic scenario may influence the rescue. In order to improve calculation speeds, the solution algorithm for the multi-objective rescue resource scheduling model can be divided into two stages. The first stage is to obtain a set of feasible routes from each rescue center to the location where the accident occurs that meets reliability constraints, and the second stage is to determine the rescue route selected and the quantity of various rescue resources provided by each rescue center.

3.1. Obtaining the Set of Feasible Routes

The travel time of the rescue route should meet constraint (12); that is, the rescue route from the rescue unit to the location of accident is bounded. In addition, when determining the vehicle routing problem from the rescue center to the accident location, the travel time and reliability of the route should be considered simultaneously. Therefore, designing a multi-path algorithm with a bound constraint is necessary in order to obtain a set of feasible routes.
  • Determination of the boundary constraint
The rescue route of each rescue center should meet the constraint of the maximum tolerable time constraint (12). Therefore, according to the definition of travel time reliability of rescue vehicle route k r s (Equation (8)) and the probability density function for route travel time (Equation (7)), the reliability constraint (12) can be transformed into the following:
E ( t 2 k r s ) + D ( t 2 k r s ) inf { G | G = φ 1 ( α ) } R T i
where φ 1 ( ) is the inverse function of the distribution function of the standard normal distribution, and inf { G | G = φ 1 ( α ) } represents the infimum of function φ 1 ( α ) .
Thus, the boundary constraint of a rescue route can be determined by Equation (18).
2.
Multi-path algorithm
Dijkstra’s algorithm is an effective shortest path algorithm, which calculates many nodes during the process of solving the shortest path. If the label information of these nodes is retained, it is conducive to the calculation of multi-path.
The rescue route studied in this paper is a simple path with a boundary constraint. Therefore, it is necessary to meet the boundary constraints when selecting temporary label nodes. Thus, this paper provides the definition of qualified nodes firstly.
Definition 1.
Let t 2 p r m represent the travel time of route p r m from the rescue center r to a node v m , and v n is the adjacent node of v m . If v n is not on route p r m and the travel time from the rescue center r to a node n satisfies the following:
E ( t 2 p r m + t m n ) + D ( t 2 p r m + t m n ) inf { K | K = φ 1 ( α ) } R T i
then v n is called the qualified node for route p r m .
Definition 1 ensures that the generated route is a simple path and meets the bound constraint. Based on Dijkstra’s algorithm and the definition of qualified nodes, a multi-path algorithm can be designed by retaining the temporary information of nodes. For convenience, p k r s is used to represent the k r s -th route from the rescue center r to the accident location s, and T ( ) and P ( ) represent temporary label and permanent label, respectively.
The specific steps of the multi-path algorithm are as follows:
  • Step 1: Let the number of routes k r s = 1, the temporary label of other notes ( v m V \ r ) T ( v m ) = + , the departure time of the rescue center T ( r ) = t 1 i r , p k r s = r , and let node r be a permanent label.
  • Step 2: Select the smallest permanent label v i and check whether each adjacent node v j of v i is a qualified note according to Equation (18). If it is a qualified node, add the temporary label of v j and T ( v j ) = T ( v i ) + d i + t i j and proceed to Step 3. If there is no qualified node, let P ( v i ) = + , and proceed to Step 3.
  • Step 3: Select the node v n with the smallest expected value among all the temporary labels, and mark it as the permanent label; let P ( v n ) = min j { E [ T ( v j ) ] } , p k r s = p k r s n .
  • Step 4: If v n is the location of accident, add route p k r s to K k r s and mark the permanent label value of nodes on this route with + , k r s = k r s + 1 ; then, proceed to Step 3. Otherwise, proceed to Step 2.
For the routes in the calculated route set K k r s , the expected values of travel time and reliability are calculated according to Equations (7), (8) and (12). Then, the dominant path is eliminated, and the rescue route set between each rescue center r and the accident location s can be obtained.

3.2. Determining Rescue Route and Rescue Resource Supply

Since the model in this paper is a bi-objective programming model, the expected value and the reliability of the route travel time should be considered when determining rescue route selection and the supply of rescue resources in this stage. One of the widely used algorithms among metaheuristics proposed to solve multi-objective APP problems is the Pareto-based non-dominated sorting genetic algorithm called NSGA-II (Non-dominated Sorting Genetic Algorithm-II) [39]. Therefore, in this stage, the NSGA-II algorithm is used to determine the scheduling of rescue resources for each rescue center.

3.2.1. Chromosome Structure

In this stage, it is necessary to determine the rescue route of each rescue center and the supply of rescue resources by each rescue center.
Accordingly, the index of rescue routes k r s and supply quantity of rescue centers E M i r , which are the decision variables, are both positive integers. Therefore, integer encoding is adopted, and CH = ( EM , k ) is designed as a chromosome. For sub-chromosome EM, a positive integer I×R matrix is used to represent the supply quantity of rescue centers. A gene in this sub-chromosome represents the supply quantity ( E M i r ) of the rescue center r for rescue resource i. For sub-chromosome k, a vector with r elements is adopted to represent the index of rescue routes for each rescue center. Figure 2 shows a chromosome with I rescue resource types and R rescue center.

3.2.2. Initialization

In the process of initial population generation, the genes of sub-chromosomes EM need to meet constraints (10) and (11). In order to improve the quality of the initial population, rescue centers are sorted from the perspective of two objective functions.
Let o i r 1 denote the sequence number of the rescue center r sorted in descending order according to the minimum expected travel time of the rescue route for rescue resource i. Correspondingly, let o i r 2 represent the sequence number of the rescue center r sorted in ascending order according to the reliability of the rescue route for rescue resource i. Generally, the higher the amount of resources provided by the top rescue centers, the better the quality of chromosomes for this object. Let us set the comprehensive ranking value of rescue centers for rescue resource i as o i r and the weight of random sorting as ω i , respectively. Then, comprehensive ranking can be defined as follows.
o i r = ω i o i r 1 + ( 1 ω i ) o i r 2
For sub-chromosome k, the first stage of this algorithm has obtained the set of feasible routes between each rescue center and the accident location, and at this stage, we only need to select an appropriate route.
The process of generating a chromosome based on random sorting is as follows:
  • Step 1: Parameter setting: set R, I, k r s , o i r 1 , o i r 2 , E D i and E N i r . Let g = 1 and ω i = 0 .
  • Step 2: For g = 1 to popsize (the number of chromosomes), repeat Step 3 to Step 5.
  • Step 3: For i = 1 to I, calculate the comprehensive ranking value o i r and sort the rescue centers according to o i r in descending order.
  • Step 4: For r = 1 to R, arrange the supply of rescue resource i in the order of o i r : let the supply of the first rescue center be E M i r = min { E N i r , E D i } , E D i = E D i min { E N i r , E D i } and calculate the supply of the subsequent rescue center E M i r = min { E N i r , E D i } . If E D i = 0 , set ω i = ω i + 1 / ( p o p s i z e 1 ) and proceed to Step 3.
  • Step 5: For r = 1 to R, generate a random integer number   k r [ 1 ,   k r s ] as a gene of sub-chromosome k.
For example, if the demand for the rescue resources E D i (I = 3) is (50, 40, 40), the rescue capacity E N i r and the comprehensive ranking value o i r of the rescue centers (r = 4) are shown in Figure 3, and the sub-chromosome of the supply quantity for rescue resources will be generated according to the above chromosome generation method. For resource 1, according to the comprehensive ranking value o 1 r , the first place is rescue center 1, and its supply quantity can be determined as min {50, 20} = 20. The second place is rescue center 2, and it can be determined that the supply quantity is min {50–20, 10} = 10. Similarly, the third place is rescue center 4, and its supply quantity is min {50–20–10, 30} = 20. The demand of resource 1 has been met, and we can obtain a feasible sub-chromosome of the supply quantity by repeating this operation, as shown in Figure 3.
Obviously, the above chromosome generation method satisfies constraints (10) and (11).

3.2.3. Crossover Operator

The crossover operation for sub-chromosome EM is carried out by a linear combination of two parent chromosome genes. Let us randomly select two parent chromosomes and randomly select a rescue resource type for crossover according to the following equation:
E M 1 i r = λ E M 1 i r + ( 1 λ ) E M 2 i r , E M 2 i r = ( 1 λ ) E M 1 i r + λ E M 2 i r ,   λ ( 0 , 1 )
where E M 1 i r and E M 2 i r are the genes of the supply quantity for selected parents, and E M 1 i r and E M 2 i r are the genes of an offspring. Figure 4 shows the crossover operation of resource 1 for sub-chromosome EM.
For sub-chromosome k, generate a random integer number r [ 1 , R ] and exchange genes in location r of two parent sub-chromosomes. The crossover process for sub-chromosome k is illustrated in Figure 5.

3.2.4. Mutation Operator

For sub-chromosome EM, the mutation operation is carried out by regenerating the random sorting weight. Let us randomly select a rescue resource type i, generate a random sorting weight ω i , calculate the comprehensive ranking value o i r and regenerate the sub-chromosome EM of i according to the generation method based on random sorting. Assume that the resource 2 is selected for mutation. Figure 6 shows the mutation process for sub-chromosome EM.
For sub-chromosome k, let us randomly select a gene and regenerate this gene in the feasible region.

3.2.5. Selection Process

In this paper, non-dominated sorting and crowding distance-based fitness assignment strategies of NSGA-II algorithm are used to implement competitive selection [40].

4. Experimental Analysis

A regional road network is shown in Figure 7. The accident location s = 1 and points 2–7 are rescue centers. Table 4 shows the resource demand of the accident location and the corresponding capacity of each rescue center. The maximum tolerable transport time for rescue resources and the assembly time of rescue resources for each rescue center are shown in Table 5. The travel time of links is listed in Table 6. In addition, the time for rescue vehicles to pass through the intersection without delay d v 0 = 0.05 min and the dissipation time of other vehicles at the intersection are shown in Table 7.
Let us set the confidence level α = 0.9 by using the multi-path algorithm with the bound constraint; the experiment is then performed on DELL Vostro 3800-R1846 (Intel Core i3 4130 (3.4 GHz 8 GB memory)) for 10 times, and the average runtime is about 10.3 s. The set of feasible routes from each rescue center to the location where the accident occurs can be obtained, which is listed in Table 8, Table 9, Table 10 and Table 11. Then, the dominant path is eliminated, and the rescue route set between each rescue center r and the accident location can be obtained.
Based on the rescue route set obtained by the first stage, the NSGA-II algorithm is used to determine the scheduling of rescue resources for each rescue center in the second stage. The size of population is 20, the maximum number of iterations is 100, the crossover probability is P c = 0.5 and the mutation probability is P m = 0.3 . Similarly, the algorithm is executed 10 times in this stage, with an average runtime of about 5 s. The Pareto solutions are shown in Figure 8. In the obtained Pareto solution set, two schemes with the minimum expected total travel time (Scheme 1) and the maximum total reliability (Scheme 7) are selected for analysis, and the results of rescue resource scheduling for the two schemes are shown in Table 12. The object values of Scheme 1 are 2943 and 386.84, and the object values of Scheme 7 are 3018.5 and 388.05. Compared with Scheme 1, the reliability of Scheme 7 is improved, but the total expected travel time increased. Rescue resource 1 of Scheme 7 needs to dispatch resources from three centers due to the capacity limitation of rescue center 4, while Scheme 1 only needs to dispatch resources from two centers.
By using the method of information entropy, synthesized objective values can be obtained, and the synthesized objective values of Scheme 2 and Scheme 3 exceed 0.9. It can be observed that Scheme 2 is relatively better in two objectives (2945, 387.92), and the rescue resource scheduling of Scheme 2 is shown in Figure 9.

5. Conclusions

The accident rescue of railway dangerous goods accident has the characteristic of time urgency, the accident itself has continuous expansion, the traffic environment has uncertainty and the accident rescue needs multiple rescue resources. Therefore, it is necessary to integrate and reasonably arrange appropriate rescue resource scheduling and routing plans. According to the analysis of travel time and reliability for rescue routes, considering the factors such as rescue capacity, rescue demand and response time, a multi-objective scheduling model of rescue resources based on travel time reliability is constructed in order to minimize the total arrival time of rescue resources and to maximize total reliability. Considering travel time and reliability of the rescue route, the proposed model is more reliable than the traditional model.
Furthermore, a two-stage algorithm is designed to solve this problem. In the first stage, a multi-path algorithm with a bound constraint is used to obtain the set of feasible routes from each rescue center relative to the accident location. In the second stage, the NSGA-II algorithm is used to determine the scheduling of rescue resources for each rescue center. Finally, a regional road network is collected to test multi-objective scheduling and algorithm. The results show that the designed two-stage algorithm can simplify problem solving and effectively improve solving speed, and it is valid for solving the rescue resource scheduling problem of dangerous goods accidents.
In order to reduce emergency response time when an accident occurs, the location in relation to transport and storage of dangerous goods should formulate an emergency plan. At the same time, emergency personnel should be trained to improve their emergency safety technical knowledge and master certain abilities in order to prevent and deal with accidents. Emergency drills should be organized regularly to improve risk awareness in dealing with emergencies and abilities related to emergency response. Furthermore, each rescue unit should reserve sufficient materials required for railway dangerous goods accidents within the emergency scope and plan the rescue route in advance.
In this paper, we assumed that the demand for the rescue resources is known. In reality, the running time may be uncertain due to the complexity of dangerous goods accidents. Moreover, travel times are time-varying and stochastic in transportation because of the effects of traffic flow and weather. Therefore, future research can consider the following:
  • Considering the uncertainty of the demand for rescue resources to design a model and algorithm of rescue resource scheduling problem;
  • Considering the travel time of links in stochastic, time-varying transportation networks to research the rescue resource scheduling problem;
  • Considering the priority of rescue equipment and resources to research on the rescue resource scheduling problem.

Author Contributions

The contributions of the authors were distributed in the following manner: Conceptualization, L.L.; methodology, X.Y. and L.L.; software, X.Y.; validation, L.L.; formal analysis, X.Y.; investigation, L.L.; data curation, X.Y.; writing—original draft preparation, X.Y.; writing—review and editing, L.L; supervision, X.Y.; project administration, X.Y.; funding acquisition, X.Y. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (Grant No. 71761024), the Natural Science Foundation of Gansu Province of China (No. 21JR1RA236) and Double-First Class Major Research Programs, Educational Department of Gansu Province (No. GSSYLXM-04).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Intergovernmental Organisation for International Carriage by Rail. Convention concerning International Carriage by Rail (COTIF). Appendix C—Regulations Concerning the International Carriage of Dangerous Goods by Rail (RID). Available online: https://otif.org/fileadmin/new/3-Reference-Text/3B-RID/RID_2021_e_01_July_2021.pdf (accessed on 1 July 2021).
  2. UIC. International Union of Railways. Dangerous Goods. Available online: https://uic.org/freight/dangerous-goods/ (accessed on 1 November 2021).
  3. Galieriková, A.; Sosedová, J.; Dávid, A.; Bariak, M. Transportation of dangerous goods by rail. In Proceedings of the Horizons of Railway Transport 2018—Proceedings of the 10th Year of International Scientific Conference, Strečno, Slovakia, 11–12 October 2018; Gašparik, J., Čamaj, J., Mašek, J., Zitrický, V., Eds.; EDP Sciences: Les Ulis, France, 2018; Volume 1. [Google Scholar]
  4. Caunhye, A.M.; Nie, X.; Pokharel, S. Optimization models in emergency logistics: A literature review. Soc. Econ. Plan. Sci. 2012, 46, 4–13. [Google Scholar] [CrossRef]
  5. Galindo, G.; Batta, R. Review of recent developments in OR/MS research in disaster operations management. Eur. J. Oper. Res. 2013, 230, 201–211. [Google Scholar] [CrossRef]
  6. Zhou, L.; Wu, X.H.; Xu, Z.S.; Fujita, H. Emergency decision making for natural disasters: An overview. Int. J. Disaster Risk Reduct. 2018, 27, 567–576. [Google Scholar] [CrossRef]
  7. Bodaghi, B.; Shahparvari, S.; Fadaki, M.; Lau, K.H.; Ekambaram, P.; Chhetri, P. Multi-resource scheduling and routing for emergency recovery operations. Int. J. Disaster Risk Reduct. 2020, 50, 101780. [Google Scholar] [CrossRef]
  8. Özdamar, L.; Ekinci, E.; Küçükyazici, B. Emergency logistics planning in natural disasters. Ann. Oper. Res. 2004, 129, 217–245. [Google Scholar] [CrossRef]
  9. Chang, M.-S.; Tseng, Y.-L.; Chen, J.-W. A scenario planning approach for the flood emergency logistics preparation problem under uncertainty. Transp. Res. Part E Logist. Transp. Rev. 2007, 43, 737–754. [Google Scholar] [CrossRef]
  10. Balcik, B.; Beamon, B.M.; Smilowitz, K. Last mile distribution in humanitarian relief. J. Intell. Transp. Syst. 2008, 12, 51–63. [Google Scholar] [CrossRef]
  11. Huang, K.; Jiang, Y.P.; Yuan, Y.F.; Zhao, L.D. Modeling multiple humanitarian objectives in emergency response to large-scale disasters. Transport. Res. Part E Logist. Transp. Rev. 2015, 75, 1–17. [Google Scholar] [CrossRef]
  12. Lassiter, K.; Khademi, A.; Taaffe, K.M. A robust optimization approach to volunteer management in humanitarian crises. Int. J. Prod. Econ. 2015, 163, 97–111. [Google Scholar] [CrossRef]
  13. Wex, F.; Schryen, G.; Feuerriegel, S.; Neumann, D. Emergency response in natural disaster management: Allocation and scheduling of rescue units. Eur. J. Oper. Res. 2014, 235, 697–708. [Google Scholar] [CrossRef]
  14. Schryen, G.; Rauchecker, G.; Comes, T. Resource planning in disaster response. Bus. Inf. Syst. Eng. 2015, 57, 243–259. [Google Scholar] [CrossRef]
  15. Zheng, H.; Gao, Y. Dispatching model and algorithm for multi-resource emergency systems. J. Univ. Shanghai Sci. Technol. 2013, 35, 415–419. [Google Scholar]
  16. Lee, K.; Lei, L.; Pinedo, M.; Wang, S. Operations scheduling with multiple resources and transportation considerations. Int. J. Prod. Res. 2013, 51, 7071–7090. [Google Scholar] [CrossRef]
  17. Zhang, L.; Lin, Y.; Yang, G.; Chang, H. Emergency resources scheduling based on adaptively mutate genetic algorithm. Comput. Hum. Behav. 2011, 27, 1493–1498. [Google Scholar] [CrossRef]
  18. Shahparvari, S.; Bodaghi, B. Risk reduction for distribution of the perishable rescue items: A possibilistic programming approach. Int. J. Disaster Risk Reduct. 2018, 31, 886–901. [Google Scholar] [CrossRef]
  19. Bodaghi, B.; Palaneeswaran, E.; Shahparvari, S.; Mohammadi, M. Probabilistic allocation and scheduling of multiple resources for emergency operations: A victorian bushfire case study. Comput. Environ. Urban Syst. 2020, 81, 101479. [Google Scholar] [CrossRef]
  20. Li, C. Optimization of emergency resource scheduling in serial terrorist attacks based on PSO-CS algorithm. J. Appl. Secur. Res. 2021, 16, 526–537. [Google Scholar] [CrossRef]
  21. Tang, Z.; Sun, J. Multi objective optimization of railway emergency rescue resource allocation and decision. Int. J. Syst. Assur. Eng. Manag. 2018, 9, 696–702. [Google Scholar] [CrossRef]
  22. Wu, L.; Xiang, X. Emergency resource scheduling during COVID-19. In Proceedings of the 2021 International Conference on Environmental and Engineering Management (EEM 2021), Changsha, China, 23–25 April 2021; Li, J., Wu, J., Chen, K., Eds.; EDP Sciences: Les Ulis, France, 2021; Volume 253. [Google Scholar]
  23. Wu, X.; Cao, Y.; Xiao, Y.; Guo, J. Finding of urban rainstorm and waterlogging disasters based on microblogging data and the location-routing problem model of urban emergency logistics. Ann. Oper. Res. 2020, 290, 865–896. [Google Scholar] [CrossRef]
  24. Zidi, K.; Mguis, F.; Borne, P.; Ghedira, K. Distributed genetic algorithm for disaster relief planning. Int. J. Comput. Commun. Contr. 2013, 8, 769–783. [Google Scholar] [CrossRef] [Green Version]
  25. Mguis, F.; Zidi, K.; Ghedira, K.; Borne, P. Distributed approach for vehicle routing problem in disaster case. IFAC Proc. 2012, 45, 353–359. [Google Scholar] [CrossRef]
  26. Liu, D.; Han, J.; Zhu, J. Vehicle routing for medical supplies in large-scale emergencies. Lect. Notes Oper. Res. 2007, 8, 412–419. [Google Scholar]
  27. Ceselli, A.; Righini, G.; Tresoldi, E. Combined location and routing problems for drug distribution. Discrete Appl. Math. 2014, 165, 130145. [Google Scholar] [CrossRef]
  28. Tang, P.; Zhu, S. Emergency logistics distribution VRP based on differential evolution optimization. J. Softw. Eng. 2015, 9, 105–113. [Google Scholar] [CrossRef] [Green Version]
  29. Zhao, J.D.; Guo, Y.J.; Duan, X.H. Dynamic path planning of emergency vehicles based on travel time prediction. J. Adv. Transp. 2017, 2017, 9184891. [Google Scholar] [CrossRef]
  30. Özdamar, L.; Demir, O. A hierarchical clustering and routing procedure for large scale disaster relief logistics planning. Transp. Res. Part E Logist. Transp. Rev. 2012, 48, 591–602. [Google Scholar] [CrossRef]
  31. He, Y.; Wen, J.; Huang, M. Study on emergency relief VRP based on clustering and PSO. In Proceedings of the 2015 11th International Conference on Computational Intelligence and Security (CIS), Shenzhen, China, 19–20 December 2015. [Google Scholar]
  32. Vargas-Florez, J.; Medina-Rodríguez, R.; Alva-Cabrera, R. A clustering optimization approach for disaster relief delivery: A case study in lima-Peru. In Proceedings of the 3rd Annual International Symposium on Information Management and Big Data, Cusco, Peru, 1–3 September 2016; Lossio-Ventura, J.A., Alatrista-Salas, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  33. Gharib, Z.; Bozorgi-Amiri, A.; Tavakkoli-Moghaddam, R.; Najafi, E. A cluster-based emergency vehicle routing problem in disaster with reliability. Sci. Iran. 2018, 25, 2312–2330. [Google Scholar] [CrossRef] [Green Version]
  34. Penna, P.H.V.; Santos, A.C.; Prins, C. Vehicle routing problems for last mile distribution after major disaster. J. Oper. Res. Soc. 2018, 69, 1254–1268. [Google Scholar] [CrossRef]
  35. Bai, W. Research on Intersection Signal Switching Model Under Emergency Situation. Master’s Thesis, Southwest Jiaotong University, Chengdu, China, May 2012. [Google Scholar]
  36. Lo, H.K.; Luo, X.W.; Siu, B.W.Y. Degradable transport network: Travel time budget of travelers with heterogeneous risk aversion. Transp. Res. Part B 2006, 40, 792–806. [Google Scholar] [CrossRef]
  37. Siu, B.W.Y.; Lo, H.K. Doubly uncertain transportation network: Degradable capacity and stochastic demand. Eur. J. Oper. Res. 2008, 191, 166–181. [Google Scholar] [CrossRef]
  38. Uno, N.; Iida, Y.; Kawaratani, S. Effects of dynamic information system on travel time reliability of road network. Traffic Transp. Stud. 2002, 2020, 911–918. [Google Scholar]
  39. Liu, L.F.; Yang, X.F. Multi-objective aggregate production planning for multiple products: A local search-based genetic algorithm optimization approach. Int. J. Comput. Intell. Syst. 2021, 14, 156. [Google Scholar]
  40. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multi objective genetic algorithm: NSGA-II. IEEE Trans. Evolut. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The flow chart of Section 2.
Figure 1. The flow chart of Section 2.
Algorithms 14 00325 g001
Figure 2. A simple example of a chromosome.
Figure 2. A simple example of a chromosome.
Algorithms 14 00325 g002
Figure 3. Illustration of sub-chromosome EM generation.
Figure 3. Illustration of sub-chromosome EM generation.
Algorithms 14 00325 g003
Figure 4. Illustration of crossover for sub-chromosome EM.
Figure 4. Illustration of crossover for sub-chromosome EM.
Algorithms 14 00325 g004
Figure 5. Illustration of crossover for sub-chromosome k.
Figure 5. Illustration of crossover for sub-chromosome k.
Algorithms 14 00325 g005
Figure 6. Illustration of mutation for sub-chromosome EM.
Figure 6. Illustration of mutation for sub-chromosome EM.
Algorithms 14 00325 g006
Figure 7. Schematic diagram of a regional road network.
Figure 7. Schematic diagram of a regional road network.
Algorithms 14 00325 g007
Figure 8. The distribution of obtained Pareto solutions.
Figure 8. The distribution of obtained Pareto solutions.
Algorithms 14 00325 g008
Figure 9. The rescue resource scheduling of Scheme 2.
Figure 9. The rescue resource scheduling of Scheme 2.
Algorithms 14 00325 g009
Table 1. Sets and indices in the model.
Table 1. Sets and indices in the model.
Sets/IndicesDescription
ISet of all resource types
i Index   of   rescue   resource   types ,   i I
RSet of rescue centers
r Index   of   rescue   centers ,   j J
sThe accident location
ESet of links
e m n Index   of   links ,   e m n E ,   m , m V , m n
VSet of intersections
v Index   of   intersections , v V
K r s Set of rescue routes between rescue center r and accident location s
Table 2. Parameters of rescue demand and rescue center.
Table 2. Parameters of rescue demand and rescue center.
ParametersDescription
E D i The demand for a rescue resource type i
R T i Maximum tolerable transport time for a resource type i
E N i r Capacity of rescue center r for a resource type i
t 1 i r Assembly time of rescue resource i of rescue center r
t m n Travel   time   of   link   e m n
d v Delay of intersection v
Table 3. Decision variables.
Table 3. Decision variables.
Parameters Description
E M i r The supply quantity of rescue center r for resource type i
k r s Index of routes between rescue center r and accident location s
t 2 k r s Travel   time   of   route   k r s from the rescue center r to the accident location s
R k r s Travel   time   reliability   of   route   k r s from the rescue center r to the accident location s
δ m n k r s A   binary   variable .   δ m n k r s = 1 means   link   e m n   on   route   k r s ,   otherwise   δ m n k r s = 0
δ v k r s A   binary   variable .   δ v k r s   = 1   means   intersection   v   on   route   k r s ,   otherwise   δ v k r s = 0
Table 4. Demand for rescue resources and rescue capacity of rescue centers.
Table 4. Demand for rescue resources and rescue capacity of rescue centers.
Rescue Resource i E D i E N i r
234567
11208060500070
21007050606000
380005050500
490500007050
Table 5. The maximum tolerable transport time and the assembly time of each rescue center.
Table 5. The maximum tolerable transport time and the assembly time of each rescue center.
Rescue Resource i R T i / min t 1 i r / min   ( E ( ) , D ( ) )
234567
115(3.5, 1.5)(3.5, 1.2)(4, 1)n.a.n.a.(3.5, 1.5)
215(4, 1.5)(4, 1.5)(5, 1.2)(4.5, 1.5)n.a.n.a.
310n.a.n.a.(2, 0.8)(1.5, 0.5)(1.5, 0.9)n.a.
410(1.5, 0.8)n.a.n.a.n.a.(2, 1.2)(1.5, 1)
Table 6. The travel time of links.
Table 6. The travel time of links.
m,n t m n 0 / min ξ m n / min   ( E ( ) , D ( ) ) m,n t m n 0 / min ξ m n / min   ( E ( ) , D ( ) )
2, 90.6(0.5, 0.05)7, 260.8(0.5, 0.1)
8, 90.5(0.3, 0.05)7, 240.4(0.1, 0.05)
6, 80.8(0.5, 0.1)7, 350.8(0.5, 0.1)
6, 100.6(0.5, 0.05)4, 280.8(0.5, 0.1)
2, 120.7(0.5, 0.08)28, 291.0(0.6, 0.1)
9, 110.6(0.5, 0.05)29, 300.5(0.3, 0.05)
11, 120.7(0.5, 0.08)30, 310.7(0.5, 0.08)
10, 111.2(0.7, 0.1)31, 320.3(0.1, 0.05)
10, 160.8(0.5, 0.1)3, 320.4(0.1, 0.05)
12, 170.9(0.55, 0.1)3, 330.7(0.5, 0.08)
16, 172.0(0.8, 0.1)4, 370.5(0.3, 0.05)
16, 210.7(0.5, 0.08)37, 380.7(0.5, 0.08)
17, 220.8(0.5, 0.1)28, 380.8(0.5, 0.1)
17, 181.0(0.6, 0.1)38, 391.0(0.6, 0.1)
15, 160.7(0.5, 0.08)38, 480.5(0.3, 0.05)
15, 201.0(0.6, 0.1)5, 480.4(0.1, 0.05)
18, 230.8(0.5, 0.1)5, 490.6(0.5, 0.05)
20, 210.8(0.5, 0.1)39, 490.4(0.1, 0.05)
1, 210.6(0.5, 0.05)39, 400.5(0.3, 0.05)
1, 220.8(0.5, 0.1)49, 500.5(0.3, 0.05)
22, 231.0(0.6, 0.1)29, 390.4(0.1, 0.05)
23, 240.8(0.5, 0.1)31, 400.8(0.5, 0.1)
20, 280.5(0.3, 0.05)40, 500.4(0.1, 0.05)
1, 300.7(0.5, 0.08)40, 410.3(0.1, 0.05)
22, 250.3(0.1, 0.05)50, 510.4(0.1, 0.05)
25, 320.5(0.3, 0.05)41, 510.3(0.1, 0.05)
25, 261.0(0.6, 0.1)32, 410.9(0.55, 0.1)
26, 330.6(0.5, 0.05)33, 340.2(0.1, 0.05)
23, 260.4(0.1, 0.05)34, 350.7(0.5, 0.08)
37, 470.5(0.3, 0.05)50, 610.7(0.5, 0.08)
47, 480.7(0.5, 0.08)51, 520.3(0.1, 0.05)
48, 500.4(0.1, 0.05)51, 420.8(0.5, 0.1)
49, 600.7(0.5, 0.08)42, 430.6(0.5, 0.08)
Table 7. The dissipation time of other vehicles at the intersection.
Table 7. The dissipation time of other vehicles at the intersection.
v ψ v / min   ( E ( ) , D ( ) ) v ψ v / min   ( E ( ) , D ( ) ) v ψ v / min   ( E ( ) , D ( ) )
6(0.15, 0.05)22(0.15, 0.05)35(0.1, 0.03)
8(0.1, 0.03)23(0.1, 0.03)37(0.1, 0.03)
9(0.1, 0.03)24(0.1, 0.03)38(0.15, 0.05)
10(0.2, 0.05)25(0.15, 0.05)39(0.15, 0.05)
11(0.15, 0.05)26(0.1, 0.03)40(0.15, 0.05)
12(0.2, 0.05)28(0.15, 0.05)41(0.1, 0.03)
15(0.1, 0.03)29(0.15, 0.05)47(0.15, 0.05)
16(0.15, 0.05)30(0.15, 0.05)48(0.15, 0.05)
17(0.15, 0.05)31(0.2, 0.05)49(0.25, 0.15)
18(0.1, 0.03)32(0.25, 0.15)50(0.15, 0.05)
20(0.1, 0.03)33(0.1, 0.03)51(0.15, 0.05)
21(0.15, 0.05)34(0.1, 0.03)52(0.1, 0.03)
Table 8. Feasible routes and their reliability for rescue resource 1.
Table 8. Feasible routes and their reliability for rescue resource 1.
Rescue Center (r)Route t 2 k r s / min   ( E ( ) , D ( ) ) R k r s
22-12-17-22-1(9.4, 2.03)0.994
2-9-11-12-17-22-1(11.95, 2.21)0.916
2-9-11-10-16-21-1(11.7, 2.16)0.937
2-12-17-16-21-1(11.6, 2.11)0.947
33-32-25-22-1(7.2, 1.7)1
3-32-31-30-1(7.55, 1.71)1
3-33-26-23-22-1(9.85, 1.72)0.999
3-33-26-25-22-1(9.8, 1.74)0.999
44-28-20-21-1(9.05, 1.43)1
4-28-29-30-1(9.5, 1.47)1
4-37-38-39-29-30-1(11.05, 1.64)0.992
4-37-38-28-20-21-1(11.4, 1.64)0.986
77-26-25-22-1 (8.65, 1.98)0.999
7-24-23-22-1(8.7, 1.96)0.999
7-26-25-32-31-30-1(11.1, 2.29)0.956
7-35-34-33-32-25-22-1(11.65, 2.4)0.919
Table 9. Feasible routes and their reliability for rescue resource 2.
Table 9. Feasible routes and their reliability for rescue resource 2.
Rescue Center (r)Route t 2 k r s / min   ( E ( ) , D ( ) ) R k r s
22-12-17-22-1(9.9, 2.03)0.994
2-9-11-10-16-21-1(12.2, 2.16)0.903
2-12-17-16-21-1(12.1, 2.11)0.915
33-32-25-22-1(7.6, 1.9)1
3-32-31-30-1(7.95, 1.91)1
3-33-26-23-22-1(10.35, 2.02)0.989
3-33-26-25-22-1(10.3, 2.04)0.989
44-28-20-21-1(10.05, 1.63)0.999
4-28-29-30-1(10.5, 1.67)0.996
4-37-38-39-29-30-1(12.05, 1.84)0.946
4-37-38-28-20-21-1(12.4, 1.84)0.921
55-49-39-29-30-1(9.1, 2.08)0.998
5-48-38-29-20-21-1(11.25, 2.13)0.961
5-49-50-40-31-30-1(11.85, 2.26)0.918
5-48-38-28-29-30-1(11.7, 2.18)0.935
Table 10. Feasible routes and their reliability for rescue resource 3.
Table 10. Feasible routes and their reliability for rescue resource 3.
Rescue Center (r)Route t 2 k r s / min   ( E ( ) , D ( ) ) R k r s
44-28-20-21-1(7.05, 1.23)0.992
4-28-29-30-1(7.49, 1.27)0.976
55-49-39-29-30-1(6.5, 1.08)0.999
5-48-38-29-20-21-1(8.25, 1.13)0.939
66-10-16-21-1(6.85, 1.33)0.991
Table 11. Feasible routes and their reliability for rescue resource 4.
Table 11. Feasible routes and their reliability for rescue resource 4.
Rescue Center (r)Route t 2 k r s / min   ( E ( ) , D ( ) ) R k r s
22-12-17-22-1(7.4, 1.33)0.975
66-10-16-21-1(7.35, 1.63)0.948
77-26-25-22-1(6.65, 1.48)0.988
7-24-23-22-1(6.7, 1.46)0.988
Table 12. The rescue resource scheduling of Scheme 1 and Scheme 7.
Table 12. The rescue resource scheduling of Scheme 1 and Scheme 7.
Rescue Resource i E D i E M i r
Scheme 1Scheme 7
356723457
112060006006050010
21005050000505000
3800503000030500
4900040504000050
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, L.; Yang, X. Travel Time Reliability-Based Rescue Resource Scheduling for Accidents Concerning Transport of Dangerous Goods by Rail. Algorithms 2021, 14, 325. https://doi.org/10.3390/a14110325

AMA Style

Liu L, Yang X. Travel Time Reliability-Based Rescue Resource Scheduling for Accidents Concerning Transport of Dangerous Goods by Rail. Algorithms. 2021; 14(11):325. https://doi.org/10.3390/a14110325

Chicago/Turabian Style

Liu, Lanfen, and Xinfeng Yang. 2021. "Travel Time Reliability-Based Rescue Resource Scheduling for Accidents Concerning Transport of Dangerous Goods by Rail" Algorithms 14, no. 11: 325. https://doi.org/10.3390/a14110325

APA Style

Liu, L., & Yang, X. (2021). Travel Time Reliability-Based Rescue Resource Scheduling for Accidents Concerning Transport of Dangerous Goods by Rail. Algorithms, 14(11), 325. https://doi.org/10.3390/a14110325

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop