Vehicle Routing Problem for Collaborative Multidepot Petrol Replenishment under Emergency Conditions

In recent years, emergency events have affected urban distribution with increasing frequency. For example, the 2019 novel coronavirus has caused a considerable impact on the supply guarantee of important urban production and livingmaterials, such as petrol and daily necessities. On this basis, this study establishes a dual-objective mixed-integer linear programming model to formulate and solve the cooperative multidepot petrol emergency distribution vehicle routing optimization problem with multicompartment vehicle sharing and time window coordination. As a method to solve the model, genetic variation of multiobjective particle swarm optimization algorithm is considered. /e effectiveness of the proposed method is analyzed and verified by first using a small-scale example and then investigating a regional multidepot petrol distribution network in Chongqing, China. Cooperation between petrol depots in the distribution network, customer clustering, multicompartment vehicle sharing, time window coordination, and vehicle routing optimization under partial road blocking conditions can significantly reduce the total operation cost and shorten the total delivery time. Meanwhile, usage of distribution trucks is optimized in the distribution network, that is, usage of singleand double-compartment trucks is reduced while that of three-compartment trucks is increased. /is approach provides theoretical support for relevant government departments to improve the guarantee capability of important materials in emergencies and for relevant enterprises to improve the efficiency of emergency distribution.


Introduction
e collaborative multidepot petrol distribution vehicle routing problem under emergency conditions (CMPDVRPE) is an extension of the multidepot petrol station replenishment problem (MPSRP) [1]. Factors such as blocked roads under emergency conditions, cooperation among petrol depots (PDs), and distribution resources sharing are considered. CMPDVRPE optimizes vehicle routing arrangement in emergency environments through multicompartment vehicle sharing and demand time window coordination (TWC) and then reduces the total network operating cost and shortens the total delivery time. As early as 1995, multicompartment vehicle has been used to study oil distribution [2], and it is widely used in oil supply. In recent years, various emergencies, such as earthquakes and traffic accidents, have been occurring all over the world and cause road traffic blocks or closures and affect the use of several road lines. In particular, the 2019 novel coronavirus  outbreak in 2020 has spread around the world, causing road closures in many cities and communities and severely affecting the transportation and supply of important production and daily necessities. Such cases show the importance of the efficient and timely supply of petrol, an important material for production and daily life.
In the existing petrol distribution network, each depot is only responsible for specific petrol stations (PSs) in the region. PDs are independent of each other and lack coordination and sharing of distribution business and resources. In an emergency, roadblocks cause detours for several distribution trucks, resulting in low distribution efficiency and timeliness of the entire network. Furthermore, in such cases, petrol is an important and necessary material, such as the replenishment demand of vehicles for disaster relief and medical treatment, and thus, ensuring its efficient and timely supply is of considerable significance for rescue and recovery. erefore, cooperation between regional PDs must be enhanced and the routing arrangement of distribution trucks should be optimized through resource sharing and business coordination of multidepots on the premise of partial road blocking. is cooperation can effectively reduce the operation cost of the entire petrol distribution network, improve the efficiency of emergency petrol distribution in the region, shorten the total distribution time in the region, and ensure the timely supply of important production and living materials such as petrol.
In this study, multicompartment truck sharing (TS), TWC, truck route detour, and a cooperation mechanism are integrated into the traditional MPSRP as CMPDVRPE. An optimal mathematical model is established to minimize the total operating cost and total delivery time to optimize the CMPDVRPE and get good results. An improved multiobjective particle swarm optimization (MOPSO) algorithm considering genetic variation (GV) is designed to achieve the near-optimal solution. Correlation results before and after optimization are compared and analyzed, and the solution of CMPDVRPE can improve the efficiency and rationality of vehicle routing arrangement of petrol distribution networks. First, petrol products have special distribution natures, such as product diversification, and cannot be mixed. Second, multidepot cooperation, multicompartment vehicle use, different TWC mechanisms, and vehicle routing optimization are comprehensively considered to propel the sustainable development of vehicle routing problem (VRP) theory [3,4] and emergency urban transportation system. e remaining parts of this paper are organized as follows. In Section 2, relevant studies on petrol distribution optimization considering multidepot cooperation and vehicle routing optimization are reviewed. In Section 3, a practical example of CMPDVRPE is presented and a mathematical model is established using notations and definitions to minimize total operating costs and total delivery time. In Section 4, an improved heuristic algorithm is introduced to solve CMPDVRPE. In Section 5, a small-scale example and a case study in Chongqing City in China are conducted to verify the applicability of the proposed methodology. In Section 6, the conclusions and future directions are provided.

Literature Review
In the past few decades, relatively little research has been carried out on petrol replenishment in academic circles. Existing research focuses on the distribution of petrol in a single depot and the variations of PS replenishment problem. For example, the multiperiod, time window (TW), trip packing, and multidepot with TWS are separately considered for PSRP [1,[5][6][7], which with its related problems attract increasing research focus for their practical importance. Popović et al. [8] proposed a heuristic algorithm based on variable neighborhood search to solve the multiproduct and multiperiod inventory routing problem (IRP) for multicompartment homogeneous vehicles. eir method proved superior to other optimization methods. Vidović et al. [9] developed a mixed-integer programming model and a heuristic method to observe its effects on multiproduct multiperiod IRP in fuel supply. Wang et al. [10] proposed a mathematical model considering petrol trucks returning multiple times to a depot, solved using a heuristic algorithm according to a local branch-and-bound search with a Tabu list and the Metropolis acceptance criterion. Wang et al. [11] constructed an adaptive large neighborhood search (ALNS) to solve the fuel replenishment problem (FRP) and performed sensitivity analysis on different features, including the number of vehicles, products, and vehicle compartments and capacities. e above studies directly relate to the replenishment problem of PSs. However, the models mainly considered the multicompartment vehicle transportation and time window assignment (TWA). Huang [12] proposed the Tabu search to solve an advanced capacitated location routing problem in a distribution network with multiple pickup and delivery routes. Derigs et al. [13] introduced a formal model, an integer programming formula, and reference set of 200 examples for a class of delimited general VRP with compartments and presented a set of heuristic component solvers. Lahyani et al. [14] optimized a rich multiproduct, multiperiod, and multicompartment vehicle route for the collection of olive oil in Tunisia using mathematical formulas, especially a precise branch-and-cut algorithm. Coelho and Laporte [15] defined and compared four main categories of the multicompartment delivery problem (MCDP). e study proposed formulas and models for specific MCDP cases and versions and then described a branch-and-cut algorithm that works for all variants. Ostermeier and Hübner [16] identified vehicle-related costs in empirical data collection and used a large neighborhood to solve grocery delivery for an extended multicompartment vehicle routing problem (MCVRP). Qi et al. [17] proposed a VRP method for large-scale TW based on spatiotemporal partitioning using a genetic algorithm (GA) to cluster largescale customers with K-medoid. Moreira et al. [18] presented a new method to solve the time window assignment vehicle routing problem (TWAVRP) where TWS is defined for multiple product segments. A mathematical heuristic method based on fixed optimization is used to solve the twostage stochastic optimization problem. Martins et al. [19] extended research on MCVRP by tackling a multiperiod environment with a product-oriented TWA and proposed an adaptive large neighborhood search as a solution. Eshtehadi et al. [20] studied the VRP with multicompartment vehicles operating from a single depot to visit customers within the chosen time period by minimizing major operational costs.
One part of CMPDVRPE literature focuses on the cooperation in distribution networks, which largely influences the arrangement of vehicle routes and the entire network efficiency. Wang et al. [21] established an intraregional oil distribution model with allocation quantity and route as decision variables to extend the multidepot half-open VRP with TWS and then used GA as a solution. Wang et al. [22] used a mixed-integer linear programming model to minimize the total operating cost of the nonempty two-echelon heterogeneous cooperative logistics network alliance and the GA-particle swarm optimization (PSO) algorithm to reallocate customer clustering units. Wang et al. [23] established a biobjective programming model to optimize the total operation routing cost and the total number of delivery vehicles for collaborative multidepot VRP with TWA (CMDVRPTWA). e model was solved using a hybrid heuristic algorithm consisting of K-means clustering, Clark-Wright (CW) saving algorithm, and an extended nondominated sorting genetic algorithm-II (E-NSGA-II).
Another part of the CMPDVRPE literature focuses on the VRP in emergency distributions. Sheu [24] presented a hybrid fuzzy clustering-optimization approach to the emergency logistics codistribution response to the urgent relief demands in the crucial rescue period. Zhang and Xiong [25] studied the routing optimization problem of grain emergency vehicle scheduling with three objectives and presented a hybrid algorithm as a solution based on combining artificial immune and ant colony optimization (ACO) algorithms. Huizing et al. [26] described a mixed-integer linear program and several increasingly refined heuristics for the optimization problem of timetabling jobs and moving responders over a discrete network. A large set of benchmark instances, both from real-life case study data and from a generator, was created for the study.
Generally speaking, the mathematical model of VRP is hard NP, which requires solutions using appropriate algorithms. Kuo et al. [27] proposed a hybrid PSO with GA (HPSOGA) for solving capacitated VRP with fuzzy demand (CVRPFD). Wang et al. [28] minimized the total cost of the two-echelon logistics distribution network using a hybrid extended PSO and GA (EPSO-GA), which combines the merits of global and local search capabilities. Zhou et al. [29] proposed two partheno GA (PGA) and adapted PSO and one state-of-the-art method to solve the multiple traveling salesman problem (MTSP). e above studies tackle plentiful MPSRP aspects but suffer from the following issues. (1) e vehicle routing optimized design procedure rarely considers the cooperation among PDs by regional partitioning. (2) Minimal attention is paid to distribution TS, multicompartment truck application, roadblocks, and transship transportation in a collaborative multidepot optimization network. (3) Single intelligent algorithm and heuristic approach are difficult to apply directly to a specific scale of CMPDVRPE with numerous PSs.
Combining observations in Table 1, the main contributions of the present study lie in the following aspects. (1) A cooperation mechanism is proposed based on the regional partitioning method and a collaborative multidepot PS replenishment VRP with emergency conditions is constructed.
(2) A mixed-integer linear programming model is established on the basis of the minimum total operating cost and total delivery time for CMPDVRPE. (3) An improved MOPSO algorithm is designed to effectively address the optimization model. (4) Finally, a small-scale example and a real-world case study are used to assess the applicability of the proposed model and approach and then compare the costs and delivery time before and after the vehicle routing optimization. In addition, this study lays a foundation and strategy for optimizing the VRP in emergencies and is conducive to improving its application range.

Problem Statement and Model Formulation
3.1. Problem Statement. CMPDVRPE integrates the problems of cooperative VRP, TS, roadblocks, and TWAs. Figure 1 illustrates a noncooperative petrol distribution network in which PDs operate independently and serve only their own customers. Roads blocked in an emergency can severely affect the timely supply of petrol. As such, long-haul deliveries are inevitable. A single PD also needs a large number of delivery petrol tankers to meet the different demands and TW requirements of its customers. is case results in a substantial decrease in distribution efficiency and increase in distribution costs. Given the numerous crosstransportations, delivery times dramatically increase and cannot be controlled. Figure 2 displays an optimized petrol distribution network with cooperation, vehicle routing optimization, roadblocks, TS, and TWAs. In this network, large tankers transfer petrol across depots and are responsible for petrol delivery to many or a single PS. is cooperation leads to more efficient and reasonable distribution routes than the noncooperative distribution network. e delivery times, required number of petrol distribution trucks, and service scope of each PD are likewise reduced.
Based on the regional partitioning method, the transport time between PDs and stations is used to cluster corresponding PSs.
en appropriate distribution trucks are selected and distribution routes are arranged according to the demand of different petrol, demand TW, and roadblocks to each PS.
Seven assumptions underline the corresponding mathematical model. (1) In a relatively short working period, each PS only generates one distribution order for petrol demand. For each kind of petrol, the demand does not exceed the maximum loading capacity of most distribution truck compartments. (2) Each customer (PS) can only be served once in a working period, that is, only one distribution truck is used for its delivery service. (3) e petrol transfer and distribution trucks have constant transportation speed. (4) During a working period, distribution trucks may be dispatched repeatedly if time permits, regardless of the transfer time across PDs. (5) e loading and unloading service times of PDs and stations are related not to the petrol type but to the truck type and operational quantity. (6) Every 10 minutes is the time unit. (7) No consideration is given to the restriction of road geometry on multicompartment vehicles.

Model Formulation.
e proposed model is mathematically formulated as an optimization problem to minimize the total cost and delivery time when each PD is assigned to serve a group of PSs and VRP with different Journal of Advanced Transportation trucks and TWS [20,33]. Tables 2 and 3 list the related notations and definitions in the CMPDVRPE optimization and the judgement variables, respectively.
CMPDVRPE is formulated as a mixed-integer linear programming model to minimize the total cost and delivery time. e cost function contains four components, namely, C 1 , C 2 , C 3 , and C 4 , which are described below.
Equation (1) shows the formulation for C 1 , which denotes the transport cost of tankers that transfer petrol between PDs during a working period: Equation (2) shows the formulation for C 2 , which denotes the transportation costs of petrol distribution trucks during a working period, including transportation costs from the PD to the PS and from PS to PS:  Equation (3) shows the formulation for C 3 , which denotes the penalty cost of the distribution truck arriving early or late at the PS, that is, the impact of the distribution truck arrival time on the cost under the premise of fixed service TW at the PS: Equation (4) shows the formulation for C 4 , which denotes the maintenance cost of all petrol tankers in a working period, and the fixed cost of the PD after considering the incentive of the government competent authority or leader of the cooperative alliance: e optimization model of CMPDVRPE is defined as follows: subject to i∈D∪S j∈S   Journal of Advanced Transportation 5 p∈P k∈K j∈S p∈P Loading capacity of petrol tanker, o ∈ O Q k e loading capacity of single compartment of truck k, k ∈ K t D gh e transport time from petrol depot g to h, g, h ∈ D t D∪S ij e transport time from the depot to the petrol station or from the petrol station to another, e service time window of petrol station, j ∈ J t ik e departure time of distribution truck k from petrol depot i, i ∈ D, k ∈ K tr ik e time of distribution truck k returning to the petrol depot i or arriving at the petrol station i, i ∈ D ∪ S, k ∈ K t ijk e driving time of distribution truck k from petrol depot i to petrol station j or petrol station i to j, i ∈ D ∪ S, j ∈ S, k ∈ K T e maximum travel time allowed for each distribution truck μ e e penalty cost per unit of time for distribution trucks arriving early μ d e penalty cost per unit of time for distribution trucks arriving late M o e annual maintenance cost of a petrol tanker, o ∈ O M k e annual maintenance cost of distribution truck k, k ∈ K q ghp During a working period, the amount of petrol p is transferred from depot g to h, g, h ∈ D, p ∈ P q jp During a working period, the demand of petrol station j for petrol p, j ∈ J, p ∈ P q kp j During a working period, the amount of petrol p delivered to the petrol station j by the distribution truck k, j ∈ S, k ∈ K, p ∈ P C kp Integer variable, the number of compartments of distribution truck k loaded with petrol p, k ∈ K, p ∈ P Q d Service capacity of depot d during a working period, d ∈ D FI i e fixed cost of petrol depot i in a working period, that is, when the petrol depot i joins the cooperative network, the fixed cost can be covered by the rewards given to the partners by the government authorities or the cooperative alliance leaders, i ∈ D N i e number of distribution trucks at petrol depot i serving petrol stations during a working period, i ∈ D N k e number of compartments for distribution trucks, k ∈ K N w e number of working periods per year When the petrol tanker o transfers the petrol p from the petrol depot g to h, the variable is 1; otherwise, it is 0, When the distribution truck k starts from the petrol depot d and passes through the route(i, j), the variable is 1; When the serving petrol depot of petrol station j is adjusted from g to h after the petrol distribution network optimization, the variable is 1; otherwise, it is 0, g, h ∈ D, j ∈ S CO i When depot i agrees to join the cooperative petrol distribution network, the variable is i∈D ∪ S j∈S i,j∈s Equations 5 and (6) show the objective functions that minimize the total cost and total delivery time of the regional petrol distribution network, respectively. Constraint (7) stipulates that each distribution truck is required to serve at most one PD and one distribution route. Constraint (8) ensures flow conservation. Constraint (9) ensures that only one kind of petrol can be loaded into each distribution truck compartment, that is, the petrol cannot be mixed. Constraint (10) ensures that the number of compartments of the distribution truck k loaded with petrol p is less than the total number of compartments. Constraint (11) indicates that the total amount of all distribution trucks that deliver to PS j is equal to the station demand for petrol p. Constraint (12) ensures that the amount of petrol p distributed by the distribution truck k to PS j is not greater than the total truckload. Constraint (13) ensures that the distribution truck does not deliver more petrol p than it loads. Constraint (14) indicates that the amount of petrol p delivered by the distribution truck to PS j should not be greater than the smaller value between its maximum loading and the demand of the PS j for petrol p. Constraint (15) ensures that the PS demand for a product can only be delivered by one distribution truck. Constraint (16) stipulates that the distribution truck k is required to deliver no more petrol than its total loading capacity to all the stations visited. Constraint (17) ensures that the amount of petrol transferred between depots is not greater than the loading capacity of tankers. Constraint (18) indicates that the amount of petrol p transferred from PD g to h is equal to the demand for p from the newly added service stations after optimization of the cooperative network. Constraints (19) and (21) specify the service time of the PD and the maximum duration constraint of the distribution truck travel. Constraints (22) and (23) constrain the time of distribution truck serving PS. Constraint (24) represents the elimination of subloop constraints. Equations (25) and (30) are judgement and decision variables.

Solution Methodology
4.1. PSO. As a method of evolutionary computing, PSO was proposed in 1995 by Eberhart, an American electrical engineer, and Kenndey, a social psychologist. By observing the foraging behavior of birds, the authors believed that at first, the birds do not know where the food but come closer and closer to the food through a kind of information exchange.
is information includes the adaptability of each bird to estimate its own position according to certain rules, ability of each bird to remember its best position, and the best position found by all the birds in the flock. e best position found by the bird group is called individual optimal "pbest," while the best position found by the whole bird group is called global optimal "gbest" [34].
In PSO, the entire population is called particle swarm, and each individual in the population is called a particle. e target is a search space, and the position of each particle is a potentially feasible solution. Each particle in space adjusts its flight path based on personal and group experience to find the optimal solution.
Suppose the population is composed of n particles, and in a D-dimensional target search space, each particle can be regarded as a point in space, then x i � x i1 , x i2 , ..., x iD represents the D-dimensional vector of particle i(i � 1, 2, ..., m). A particle flies in space at a certain speed, as v i � v i1 , v i2 , ..., v iD . Generally, the fitness represented by the objective function is used to calculate the current particle fitness value, to judge the advantages and disadvantages of the particle position. e local optimal particle position is pbest � pbest i1 , pbest i2 , ..., pbest iD . At the same time, the best adaptive value experienced by the entire particle swarm is the global optimal value, gbest � gbest i1 ; gbest i2 , ..., gbest iD }.
e particle optimizes the search process iteratively, and for each of its generation, the d(1 ≤ d ≤ D) dimensional position and velocity are updated by the following equation iteration [34]: where c 1 , c 2 are the learning coefficients, also known as the acceleration coefficients, and r 1 , r 2 � ran d[0, 1], namely, r 1 , r 2 are the random numbers between 0 and 1.

MOPSO.
e single-objective optimization problem only needs to obtain a single or a set of continuous optimal solutions, while the multiobjective optimization problem generally obtains a set of continuous solutions. For the basic PSO, all particles converge in one direction following the best particle (leader), and thus, the single-objective optimization problem can be solved. When solving the multiobjective optimization problem, pbest (individual optimal solution) and gbest (globally optimal solution) cannot be directly determined due to the lack of a single optimal solution, and thus, the basic PSO method cannot be directly adopted.
In fact, the solution of the multiobjective optimization problem is composed of a set of noninferior solutions. Different particles seek different "leaders" in the optimization. e noninferior solutions are stored in the "external container" (archive). Furthermore, by evaluating the "quality" of these solutions, the gbest is determined to realize the guidance of particle iterative updating. Clearly, the archived solution can represent the optimal subset of particles that the algorithm finds in each generation of particle swarm, and thus, considering the archived set as a candidate set for the gbest is appropriate [35].
MOPSO is a cyclic process. When the algorithm starts, the first-generation noninferior solution is found from the initialized particle swarm and stored in the external archive. Subsequently, the gbest is selected from the external container by comparison, and then, the particle swarm updates the speed and position through equations (31) and (32). en, the loop operation is repeated. Based on the generation of new nondominant solutions after each particle iteration, the external archive set size, and thereby the computation, increases after multiple iterations. erefore, considering the effect of computational complexity, limiting the size of external collections is necessary [35]. Once the number of solutions in the external archive collection reaches the upper limit set for the external container, the noninferior solutions in the external container are compared and a few are deleted according to certain rules.
In solving MOPSO, the evaluation criteria of particles, construction, and preservation of noninferior solutions, selection of optimal positions, and the processing of constraints are mainly used.

IMOPSO.
PSO is characterized by rapid convergence speed, but for a multiobjective optimization problem, too fast convergence speed may cause particles to fall into the local optimal. us, achieving the global optimal becomes difficult. is problem may be prevented by introducing the mutation operation of particles.
In the early stage of the algorithm, the particle needs to search the entire target space to ensure the escape from the local optimal state and enhance the global search capability; in the late stage, a quick convergence is necessary to find the optimal solution. us, this study carries out mutation operation on particles through mutation probability. At the beginning of the algorithm, mutation should be carried out on all particles, and with the increase of iteration, the number of particles undergoing mutation should be gradually reduced. According to the variation operation requirements, the variation probability of particle swarm after each update is set as follows [36]: where mp is the mutation probability, P it is the current iterations, and max it is the maximum number of iterations.
In reference to the nonuniform mutation operator in the GA, a random variable rand between 0 and 1 is assigned to each particle. If the random variable of particle i is less than the probability of variation, then random nonuniform variation operation is carried out on the D-dimensional particle position: In equation (34), λ is the system random parameter (λ≻1). During the variation operation, if the result is greater than the threshold set by the system, the variation result is set at the boundary of the threshold range.
In addition, to balance the global and local search capability of basic PSO, adding inertia weight is necessary [37], that is, the influence of inertia weight coefficient is considered when updating the velocity of particles. erefore, the new velocity updating formula of particles is as follows: Combined with the basic PSO algorithm and MOPSO, Figure 3 describes the procedure of the IMOPSO algorithm. e steps of the MOPSO algorithm considering GV are as follows: (1) Initialization: set the particle swarm size, parameter coefficient, threshold, maximum number of iterations, and initial velocity and initial position of each particle. (2) Use equations (33) and (34) to carry out variation on particle swarm.  Journal of Advanced Transportation (8) Update the gbest of particle swarm. In the external archive, the front-end particles in descending order of crowding distance are randomly selected as the gbest. (9) Update the particle velocity and position according to equations (31) and (32), and return to step (2). When the maximum number of iterations is reached, the loop is terminated.

Small-Scale
Example. e quality of the proposed heuristic algorithm in solving the problem of multidepot vehicle routing with emergency effects is further demonstrated by a small example analysis. In this case, the petrol distribution network consists of two PDs and 31 PSs in the region. Under noncooperation, PD P is responsible for P1, Start Initialization: set the particle swarm size, parameter coefficient, threshold, maximum number of iterations, and initial velocity and position of each particle Using equations (33) and (34) to take mutation operation on the particle swarm Calculate the fitness value for each particle  Journal of Advanced Transportation P2, . . ., P16 and other 16 PSs petrol distribution, while PD Q is responsible for Q1, Q2, . . ., Q15 and other 15 PSs petrol distribution. With cooperation according to the customer clustering and the distance between each PS and the PD, the PSs distributed by the two PDs are first repartitioned. In addition, distribution routes are blocked under emergency control. en, the two PDs can share the petrol distribution trucks. e customers and routes of the distribution truck can be decided according to the demand TW, demand characteristics (including petrol products and quantity demand), and transportation distance and detour route of each petrol station in one distribution task, as shown in Figure 2.
For the convenience of calculation, the following assumption is made for this example: the demand of petrol product of PSs is divided into three categories, namely, nos. 92, 95, and 98. e distribution trucks have three types: single compartment with a capacity of 2,000 gallons, double compartment with a capacity of 1,500 gallons, and three compartments with a capacity of 1,500 gallons. Each time unit represents 10 minutes, and the transportation cost per unit time is $50 for a single-compartment truck, $70 for a double-compartment truck, and $80 for a three-compartment truck. e fixed cost of single-compartment truck, double-compartment truck, and three-compartment truck is $20, $30, and $40, respectively. In addition, considering the customer demand TW, the penalty cost of the distribution truck arriving early per unit time is $10, and the penalty cost of the distribution truck arriving late per unit time is $15. e single assignment cost of the distribution truck is $20. For example, the vehicle route P ⟶ Q4 ⟶ P7 ⟶ P6 ⟶ P uses a three-compartment vehicle, the entire route vehicle needs to run 11 time units, transportation cost is $880, fixed cost is $40, and the vehicle assignment cost is $20. According to (5), the total cost of the route is $1,160. Table 4 shows that compared with the condition of noncooperation, the cooperative multidepot petrol distribution network has a significant decrease in the delivery time, number of distribution trucks, transportation cost, and total cost. Total delivery time under the cooperative condition dropped to 1080 minutes, which was 52% lower than that under the noncooperative condition. With cooperation, the number of trucks required for distribution and use of single-and double-compartment trucks considerably decreased whereas the use of three-compartment trucks increased to a certain extent. Overall, the number of needed distribution trucks declined from 31 to 14, fixed cost declined from $1470 to $760, and transportation cost declined from $15,180 to $7,800. In the case of joint distribution, the integrated use of trucks and the reasonable arrangement of vehicle routes need consideration, resulting in a penalty cost of $350 for trucks arriving early or being late. In general, with cooperation, the total cost of dual depots petrol distribution network is $8,910, which is significantly reduced by 46% compared with the noncooperative condition. In addition, IMOPSO performed 30 random operations, with the optimal result on the 23rd. e results show that considering TS and route optimization, the petrol distribution network can achieve significant cost savings, considerably reduce the total delivery time, and significantly reduce the number of distribution trucks. Table 4 also reveals two interesting phenomena. First, the cooperative distribution network considerably reduced the use of single-and double-compartment trucks and increased the use of three-compartment trucks.
is finding indicates that the use of distribution trucks is coordinated through the integration of customer requirements, and multiple PSs can be served by sending one distribution truck at a time. Second, the penalty cost has increased because the joint distribution comprehensively considers the customer TW and truck route arrangement and the detour of trucks under emergency road blocking. In certain cases, trucks arrive early or late, but the overall delivery time is reduced and the benefits are optimized. Of course, attention should also be paid to the impact of early or late arrival on customer satisfaction.

Large-Scale Example.
In this section, a numerical experiment is applied to a large-scale multidepot petrol distribution network in Chongqing, China. e urban petrol distribution network in Chongqing municipality, which is directly under the jurisdiction of western China, is selected as the experimental object to demonstrate the effectiveness of IMOPSO in solving the vehicle routing optimization for the emergency joint distribution of multidepot. In this case, the petrol distribution network consists of five PDs (PD1, PD2, . . ., PD5) and 86 PSs (PS1, PS2, . . ., PS86). Figure 4 shows the distribution of all PDs and PSs, and each PS has a fixed demand TW. Specifically, PD1 serves the PSs represented by the circle, PD2 serves the PSs represented by the rectangle, PD3 serves the PSs represented by the triangle, PD4 serves the PSs represented by the cross, and PD5 serves the PSs represented by the hexagon. Table 5 shows PSs initially served by different PDs.
As mentioned in the previous Section 3.2, the optimization objective is to minimize the total cost and time of the distribution network through the cooperation of different PDs, customer clustering, TW coordination, integration of distribution trucks, and vehicle route optimization under emergency road blocking. In this numerical experiment, a petrol tanker with a loading capacity of 5,000 gallons was selected, and the distribution trucks were single-compartment k1, double-compartment k2, and three-compartment k3. Relevant parameter values of the optimization model and its algorithm are adopted from the literature [11,23,33] as follows: FI 2 � 1300, FI 3 � 1500, FI 4 � 1000, FI 5 � 1200, population pop − size � 500, learning factor c 1 � c 2 � 2, maximum iteration number max − it � 100, and inertia weight iw � 0.9. In addition, the service time of PD is from 6:00 to 24:00. Table 6 shows the service TW of PS while Table 7 shows the demand for different petrol products of each station in the working period.
In this study, a working period can be regarded as a working day. Based on the participation of different cooperation subjects, five PDs can have 2 5 − 1 different cooperation alliances. First, the PSs that are responsible for distribution in each PD are reallocated through a clustering algorithm, and then the assignment and route optimization of distribution trucks are calculated using the IMOPSO algorithm. us, the optimized distribution cost, delivery time, and the number of different distribution trucks used per time in a working period are obtained. For the convenience of comparison and analysis, Table 8 and Figure 5 show the cost optimization results of the different cooperative alliances. Figure 5 shows that the total operating cost of all alliances decreased compared with that of individual PD. In particular, when all PDs are willing to cooperate and a grand alliance can be formed, the reduction in operating costs is significant even if local truck detours are caused by emergency road blocking. Compared with the initial network for the entire petrol distribution, the cost decreased by 65.59%. erefore, cooperation can benefit the operation costs reduction of multidepot petrol emergency distribution network.     Noncooperative  2240  12  15  4  1470  15180  -16650  Cooperative  1080  2  4  8  760  7800  350  8910  Table 9 and Figure 6 show that the total delivery time of all alliances decreased compared with that of individual PD. In particular, when all PDs are willing to cooperate and a grand alliance can be formed, the total delivery time savings are significant even though several PSs have longer delivery times due to detours. Compared with the initial network for the entire petrol distribution, the delivery time decreased by 56.59%. erefore, the optimization of vehicle route can benefit the saving of delivery time in multidepot petrol emergency distribution network. Table 10 shows the changes in distribution truck usage before and after optimization in different cooperative alliances, mainly including the change of demand truck type and corresponding quantity. Specifically, distribution trucks are shared through cooperation between PDs. Based on the service TW of different PSs and the quantity of petrol required, the usage and route arrangement of different truck types are reasonably planned. As a result, usage considerably decreased for single-compartment trucks (k1), decreased to a certain extent for double-compartment trucks (k2), and significantly increased for three-compartment trucks (k3). For example, for the grand alliance {PD1 PD2 PD3 PD4 PD5}, usage of single-compartment trucks (k1) significantly decreased from 27 trucks to 6 trucks per time, that of doublecompartment trucks (k2) decreased from 36 to 17 trucks per time, while that of three-compartment trucks (k3) increased from 27 to 35 trucks per time. e change in the number of trucks is mainly due to the change of petrol distribution  1289 from one-to-one assignment to the one-to-many vehicle route optimization between the PD and the PS to coordinate the station TW and the demand for petrol products. us, the use of single-and double-compartment trucks decreased and that of three-compartment trucks increased. Figures 7 and 8 compare the use of distribution trucks k1 and k2 before and after optimization. After optimization, usage of both k1 and k2 decreased to a certain extent, while that of k1 decreased more. Due to the realization of the coordination between the service TW of the PS and petrol demand, in most cases, a distribution truck from the PD has to serve multiple PSs, and thus, the use of single-compartment trucks decreased. Figure 9 shows that, relative to the decreased usage of optimized distribution trucks k1 and k2, the use of distribution truck k3 increased after optimization. is result is mainly due to two reasons, the demand of PSs for three types of petrol and the vast majority of distribution trucks needed to serve multiple PSs at a time.
A large petrol distribution network consisting of five PDs and 86 PSs is relatively complex. erefore, for the optimized truck distribution routes, the distribution truck route network formed by PD3 is taken as an example, as shown in Table 11. Before cooperation, the PD was responsible for one-to-one distribution of petrol products to PSs and was only responsible for the distribution of PSs operated by itself. After cooperation, PD3 reduces the distribution of longdistance PSs through clustering, such as PS61 and PS62, and increases the distribution of PS31, PS32, and other PSs near the PD. Second, in most cases, a distribution truck serves multiple PSs at once. For example, in the vehicle route PD3 ⟶ PS21 ⟶ PS22 ⟶ PS31 ⟶ PD3, a three-compartment truck starts from PD3 in the PD and serves PS21, PS22, and PS31 at one time. In addition, considering possible road emergency blocking, detour routes for distribution trucks should be designed. In the vehicle route PD3 ⟶ PS33 ⟶ PS29 ⟶ PS25 ⟶ PS29 ⟶ PD3, if the road from PD3 to PS29 is blocked or closed due to emergency, then the distribution truck needs to detour PS29 back to PD3.

Implication.
e optimization of coordinated PS replenishment vehicle route based on regional partitioning and reasonable resource sharing promotes the sustainable development of petrol distribution and emergency energy supply system. rough the optimization of CMPDVRPE, the petrol distribution service area is reasonably divided, long-distance and cross-traffic are reduced, and the The demand quantity of k1 before cooperation (trucks/time) The demand quantity of k1 after cooperation (trucks/time) The demand quantity of k2 before cooperation (trucks/time) The demand quantity of k2 after cooperation (trucks/time)  The demand quantity of k3 before cooperation (trucks/time) The demand quantity of k3 after cooperation (trucks/time) PD3 ⟶ PS83 ⟶ PD3 * e petrol station passed by vehicles on a detour. distribution time is shortened. us, costs are minimized and additional benefits are provided to each PD. Effective vehicle routing arrangement and sharing matching strategy between vehicles and PSs are the important characteristics of this network. ese improvements greatly enhance energy and social resources, cost savings, and emergency response capabilities for PD operators and transportation management. Cooperation among logistics facilities plays an important role in optimizing the distribution in cases of emergency [42][43][44]. Further integration of transportation resource sharing and optimization of vehicle routing can further save on costs. In addition, traffic management policies that encourage joint distribution are also a sign of political will to achieve sustainable development in administrative areas [31,32]. As one of the major development factors, petrol distribution activities can be further organized under the coordination of vehicle routes, reducing the number of petrol distribution trucks. erefore, encouraging the formation of grand coalitions is a relevant method that can benefit not only PDs, but also society as a whole.

Conclusions
is study proposes an effective method to solve the CMPDVRPE optimization, which improves the cooperation of PDs and the efficient distribution vehicle routing optimization in emergencies. rough the cooperation between PDs, regional petrol joint distribution, optimization of distribution, and resource sharing can be formed. In the optimization, PS clustering mechanism, road blocking, and TS mechanism are considered. Comparison of the data before and after the cooperation shows that the total operating cost, total delivery time, and the total number of delivery vehicles are significantly decreased. e optimization model considers customer clustering, multicompartment TS, roadblocks, and TWC, thereby reducing the overall transport distance and the number of trucks used. Taking the regional petrol distribution network in Chongqing as an example, the application of the model and method is evaluated. A heuristic algorithm IMOPSO is proposed, and a case study of different scales of petrol distribution networks was carried out. e operating cost, delivery time, and the number of different types of trucks are compared before and after optimization.
In summary, the optimization of petrol emergency distribution vehicle routing is consistent with reality. e proposed optimization method is superior to the existing research in this field. On the basis of the analysis, the following conclusions are drawn. (1) rough customer clustering, multicompartment TS, vehicle routing optimization, and TWC, the regional petrol distribution network can considerably reduce the delivery time, number of distribution trucks required, and the total operating cost of the petrol distribution network. (2) Optimizing the PSs for each PD and sharing trucks when the TW demand allows can reduce traffic pressure in urban areas and its negative impact on the energy supply system and contribute to the sustainable development of urban traffic. (3) Timely and efficient supply of petrol is guaranteed through optimization of vehicle routing of petrol emergency distribution. e results of this study point to interesting research directions for the future. e following views can be considered. (1) is study only examines the cooperation between PD and PS in the secondary distribution of petrol.
us, the cooperation between the two sides can extend to the transport energy supply chain. (2) Consistent with most existing joint distribution literature, this study assumes a constant transportation speed of petrol distribution trucks. Future research can consider real-time urban traffic speed analysis to obtain more realistic results. (3) From the perspective of vehicle-road interaction, the influence of road geometry on the selection of multicompartment vehicle types can be considered in the future. (4) In the future, a dynamic CMPDVRPE model can be established by considering the spatiotemporal change of roadblocks or congestion.
Data Availability e service time window and demand quantity data used to support the findings of this study are available from the corresponding companies and administrative departments.

Conflicts of Interest
e authors declare that they have no conflicts of interest.