Next Article in Journal
File System Support for Privacy-Preserving Analysis and Forensics in Low-Bandwidth Edge Environments
Previous Article in Journal
Robust and Precise Matching Algorithm Combining Absent Color Indexing and Correlation Filter
Previous Article in Special Issue
Ship Classification Based on Improved Convolutional Neural Network Architecture for Intelligent Transport Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Customized Shuttle Transit Mode Responding to Spatiotemporal Inhomogeneous Demand in Super-Peak

School of Traffic and Transportation, Beiing Jiaotong University, No. 3 Shang Yuan Cun, Hai Dian District, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Information 2021, 12(10), 429; https://doi.org/10.3390/info12100429
Submission received: 13 August 2021 / Revised: 22 September 2021 / Accepted: 25 September 2021 / Published: 18 October 2021

Abstract

:
Instantaneous mega-traffic flow has long been one of the major challenges in the management of mega-cities. It is difficult for the public transportation system to cope directly with transient mega-capacity flows, and the uneven spatiotemporal distribution of demand is the main cause. To this end, this paper proposed a customized shuttle bus transportation model based on the “boarding-transfer-alighting” framework, with the goal of minimizing operational costs and maximizing service quality to address mega-transit demand with uneven spatiotemporal distribution. The fleet application is constructed by a pickup and delivery problem with time window and transfer (PDPTWT) model, and a heuristic algorithm based on Tabu Search and ALNS is proposed to solve the large-scale computational problem. Numerical tests show that the proposed algorithm has the same accuracy as the commercial solution software, but has a higher speed. When the demand size is 10, the proposed algorithm can save 24,000 times of time. In addition, 6 reality-based cases are presented, and the results demonstrate that the designed option can save 9.93% of fleet cost, reduce 45.27% of vehicle waiting time, and 33.05% of passenger waiting time relative to other existing customized bus modes when encountering instantaneous passenger flows with time and space imbalance.

1. Background

In megacities with complex travel demand structures, the transportation organization problems triggered by instantaneous peak passenger demand has troubled planners and operators for a long time. Despite the proved successfulness of fixed-route transit (FRT) in providing fixed and scheduled service, their performance is unsatisfying when dealing with irregular passenger flow with uneven spatiotemporal distribution. Their bad performance results from the stableness of service capacity which FRT provides, and the fact that it is difficult for FRT to do demand-responsive adjustments to operation lines. Thus, as the transportation capacity of FRT system always appears to be insufficient when coping with flow peaks, large, irregular passenger flow with uneven spatiotemporal distribution has always been the main cause of urban traffic problems under the current FRT system.
Under this background, Customized Shuttle Transit (CST) provides a new option for cities. CST is a user-oriented, demand-responsive, and individual customized transportation service, which generates operating routes and stops by collecting and calculating real-time requests. Thus, CST could operate efficiently in satisfying fluctuant and unpredictable requests, in dealing with peak flow, and providing services in the FRT-vacant period. In this way, the implementation of CST could alleviate peak traffic pressure, cover the unserviceable circumstances of the FRT, and provide an option to satisfy instantaneous peak passenger flow.
However, the distribution of requests submitted to CST system can be scattered, thus it might be difficult to integrate different customer demands into same bus routes. In this way, CST operators need to achieve balance between passenger satisfaction and operating cost. This is always a difficult balance to reach: improving passengers’ satisfaction always results in an increase in operational cost as well. Here is what always happens: in order to provide satisfactory services, the traveling time has to be shortened, as a result the number of passengers CST could serve per unit time will decrease, which will increase the average operating cost per customer, as operators has to spend more on vehicles and drivers, vice versa.
How to integrate transportation resources effectively to reach balance between service quality and operating cost has been a difficult problem. Transferring could help to alleviate this problem by splitting one service unit into one pick-up section and one delivery section and satisfy them separately, to avoid long-distance transregional travel of vehicles. However, new problems such as longer waiting time also emerges when transferring is introduced into the service mode. According to this, a new route planning system which could at the same time guarantee the service quality, keep the operational cost relatively low and provide adequate service capacity is in urgent need. To this end, this paper proposes a shuttle bus service that can gather and transport passengers efficiently. It can reduce operating costs and improve passenger satisfaction in facing the instantaneous mega-demand of spatiotemporal unevenness.
The contributions of this paper are as follows:
I.
Aiming at solving the problems caused by instantaneous super-peak passenger demand, this paper proposed a shuttle transit mode that organize passenger flow by a boarding-transferring-alighting framework. This framework, which targets on minimizing total operating costs and maximizing passenger satisfaction, can be well adapted to the needs of spatiotemporal uneven passenger flows;
II.
Considering time window, capacity, etc., constraints, a model based on PDPTWWT is proposed. Furthermore, a heuristic algorithm based on ALNS algorithm and Tabu-search algorithm is proposed to solve large scale problems. The proposed algorithm is proven outperform commercial solver;
III.
The efficiency of the proposed method is verified by comparing the metrics with other models based on different spatiotemporal homogeneity cases on the background of the actual scale case of Beijing.
The rest of this paper is organized as follows. Section 2 surveyed the works related to model establishment and algorithm construction for related problems. Section 3 constructed a model for shuttle bus route planning for instantaneous large passenger flow. A heuristic solution combined tabu search algorithm with ALNS algorithm was developed in Section 4. A set of numerical experiments were given in Section 5 to verify the validity and the effectiveness of the approach. Section 6 concluded the paper and discussed the further directions.

2. Literature Reviews

The study of customized bus system was carried out 60 years ago, while no project has landed until the 1970s. In 1973, Operators in Wisconsin, America has launched “Merrill-Go-Round” as the very first customized bus system, which is still being used at present. After that, more customized bus systems were launched. In America, ADA (Americans with Disabilities Act) paratransit bus service emerged in accordance with the ADA campaign of 1990. Dial-a-ride Transit was launched by Hopelink, which can be reserved through phone call. Other successful examples are OmniLink service, Peninsula Transit and Hampton Roads Transit (HRT) in America, Kan-Go, Flexible Transport System and Pocket Ride in Australia, and Winnipeg Transit in Canada.
The model of CB route planning problem (CBRPP) is a special case of PDPTW (pickup and delivery problem with time window) or PDPTWT (pickup and delivery problem with time window and transfers). The PDPT has been intensively studied over recent years.
In model establishment, Masson et al. constructed three mathematical models for the PDPS (The Pickup and Delivery Problem with Shuttle routes). Mitrović-Minić [1] proved the usefulness of transshipment points in the context of pickup and delivery problems with time windows. Deleplanque [2] introduced transshipment into traditional dial-a-ride problems (a variant of VRP problems) and solved it by an algorithm based on insertion techniques and constraints propagation [3].
In algorithm designs, exact algorithms are always used in cases. Sol (1994) discussed about general column generation techniques in solving pickup and delivery problems [4]. Masson et al. proposed a branch-and-price algorithm to solve the PDPS model they constructed [1]. Naccache et al. solved multi-PDPTWT firstly by a branch-and-bound algorithm, and then developed a hybrid adaptive large neighborhood search with improvement operations to solve large scale problems [5]. Parragh and Schmid modeled the underlying optimization problem of proposing specialized transportation systems to complement the traditional public transportation into a dial-a-ride problem, and proposed a hybrid column generation and large neighborhood search algorithm to solve it. Ropke et al. [6] introduced two new formulations to improve model performance for the PDPTW and dial-a-ride problem (DARP) problems, and developed branch-and-cut algorithms to solve instances with up to 8 vehicles and 96 requests [7]. However, the problem scales that exact algorithms can solve are less than 100 demands, and transferring are not allowed. Therefore, when facing large-scale problems, heuristic methods are always introduced.
Heuristic algorithms transcend in the ability of solving large-scale problems and are more often proposed to deal with more complicated matters. Among the researches in recent years, genetic algorithm (GA) [8,9,10], tabu search algorithm (TS) [11,12,13,14], ant colony optimization algorithm (ACO), simulated annealing algorithm (SA), Particle Swarm Optimization (PSO) [15,16,17,18], Adaptive large neighborhood search algorithm (ALNS) [5,6,19,20,21], and dynamic programming [22] are most frequently used.
GA has been studied by many scholars including Tan [8], Baker et al. [9] and Ombuki [10] on VRP and VRP-related problems. TS has been studied by Cordeau to solve dial-a-ride problem [11], by Joséto to solve open VRP problems in which vehicles do not have to go back to ending depot after delivering goods to customers [12], and by Lai [13] and Taillard [14] on PDPTW-related problems as well. ACO has been studied by Gambardella [15], Donati et al. [16], especially by Li to solve a multi-depot green vehicle routing problem (MDGVRP) [17] and by Zhang to solve the multi-objective vehicle routing problem [18]. SA has been studied by Bachem [23], Vincent [24], Chao [25], Chiang et al., and especially by Russell to study VRPTW [26]. PSO has been studied by Kachitvichyanukul to solve the GVRP-MDPDR (generalized multi-depot vehicle routing problem with multiple pickup and delivery requests) [27] and by Dridto to solve MDPDPTW (pickup and delivery problem with time windows and multiple depots) [28]. ALNS has been more frequently used in large scale PDPTW-related problems, which has been used by Ropke and Parragh [19], Ghilas [20], Côté [5], and Parragh and Schmid [6], and especially by Masson to solve PDPT (pickup and delivery problem with transfers) problems [21]. Dynamic programming has been studied by Ritzinger to solve VRP related problems as well [22].
The applicability of different heuristic algorithms in solving different PDPTW-related problems could be concluded from the researches above. GA’s performance is satisfying when problem is less complicated and problem scale is relatively small, but when problem scale grows, the searching procedure prolongs rapidly. ACO has strong robustness, and performs well in multi-objective optimization of PDPTW-related problem, but is relatively easy to sink into local optimum when required iteration becomes more. SA performs very well in small scale PDPTW-related problems with the advantage of strong robustness in initial stage, but its performance also declines when the scale of problem gets larger. TS and ALNS, however, with the advantage of low computational complexity, have been proved to be efficient in solving large scale PDPTW-related problems, especially ALNS is always used as an alternative to rigorous algorithms when problem scale gets larger. In conclusion, when it comes to large scale complicated single objective optimization problem, researchers always prefer TS or th eALNS algorithm.
It could also be concluded that among the PDPTW-related problems, traditional PDPTW draws most of the attentions, while PDPTW with transfers is very rarely mentioned, which is mostly solve by adaptive large neighborhood search (ALNS). Due to the spatiotemporal uncertainty and the short pickup and delivery time-windows, the transporting demand which generated by spatiotemporal homogeneity peak passenger flow is hard to satisfy. Thus, the service pattern of shuttle transit which is proposed to serve this demand might be different with the existing ones, with demand-responsive route planning functions, multiple flexible depots, and transferring allowed.

3. Customized Shuttle Transit for Instantaneous Large Passenger Flow Model Description (CSTILP)

To serve instantaneous large passenger flow, massive transportation resources is needed, and vehicles always need to travel through large acreage to satisfy every passenger. These make it difficult for transportation system to provide service for instantaneous large passenger flow. To get this problem solved, an operational system with transferring is needed, which could shrink proportion of transregional transportation by transferring and downsize the vehicle fleet.
The key point of the establishment of a transferring-allowed model is to avoid passengers or vehicles waiting long time at transfer point, and to guarantee the consistency of service after transferring has taken place.

3.1. Notations

3.1.1. Variables

We present the variables and the description of variables in Table 1.

3.1.2. Sets

We present the sets and the description of sets in Table 2.

3.1.3. Parameters

We present the parameters and the description of parameters in Table 3.

3.2. Model

We formally formulate the LIPPD as follows:
Objectives:
m i n f ( x ) = m i n ( C f + C o p + C p w t + C v w t + C q u a )  
subject to
j : ( i , j ) A x i j k j : ( j , i ) A x j i k = 0 ,   k K , i V N o D e ,
j : ( i , j ) A x i j k j : ( j , i ) A x j i k = z n k , k K ,   i V D e + , n N ,
i : ( i , j ) A x i j k i : ( j , i ) A x j i k = z n k , k K ,   j V D e ,   n N ,
k M z n + k = k M z n k = 1 ,   n N ,
j : ( n + , j ) A x n + j k = z n + k , k M ,   n N ,
j : ( n , j ) A x n j k = z n k ,   k M , n N ,
z n k + t T S , l M w t n l k t T S , l M w t n k l = z n + k , k K , n N ,
j V x d e k j k = 1 ,   k K p , d e k i s   t h e   s t a r t i n g   d e p o t   o f   k ,
p + D e + i V x i p + k = 1 , k K ,
i V , k M x i d e + k k d e , d e + D e + ,
i V x i t k 1 , k M , t T S ,
i V x i t k l M w t n k l 0 , k K , t T S , n N ,
i V x i t l k M w t n k l 0 , l K , t T S , n N ,
i V x i n l k M w t n k l 0 , l K , t T S , n N ,
i V x i n k l M w t n k l 0 , k K , t T S , n N ,
D i k + t i j D j k + o p j M · ( 1 x i j k ) ,   i V N   ,   j V , k K ,
D i k + t i d e k + k i V ( x i d e + k · l d e + ) + o p d e k + M · ( 1 x i d e k + k ) , i V N o D e , k K , d e + D e + , d e k + D e k + ,
D t k + t w p + t s o p t D t l M · ( n N w t n k l 1 ) , t T S ,   k , l K ,
D t k t w c + t s o p t D t l M · ( 1 n N w t n k l ) , t T S , k , l K ,
D n + k D t k M · ( 1   l M w t n k l ) , n N , t T S , k K ,
D t l D n l M · ( 1   k M w t n k l ) , n N , t T S , l K ,
D n + k D n l M · ( 1 t T S w t n k l ) ,   k , l K , n N ,
D n + k D n k M · t T S , l M w t n k l ,   k M , n N ,
y j k y i k q j M · ( 1 x i j k ) , i V N , j V ,
y j k y i k q j M · ( 1 x i j k ) , i V N , j V ,
y t k y i k + l M , n N ( w t n k l q n ) l M , n N ( w t n l k q n ) M · ( 1 x i t k ) , i V , k K , t T S ,
y t k y i k + l M , n N ( w t n k l q n ) l M , n N ( w t n l k q n ) M · ( x i t k 1 ) , i V , k K , t T S ,
i V , j V x i j k d i j d m a x k , k M
i V , j V x i j k t i j k t m a x k , k M ,
y d e k = 0 , d e D e ,
D d e k = 0 , d e D e ,
e i D i k l i , i V N ,
0 y i k Q k , k K , i V ,
x i j k { 0 , 1 } , i , j V ,
z n k { 0 , 1 } , n N ,
w t n k l { 0 , 1 } , n N .
Objective function (1) is composed of C f   ( fixed cost), C o p (operating cost), C p w t (passenger waiting cost), C v w t (vehicle waiting cost), and C q u a (service quality cost).
(1) C f   ( fixed cost)
The fixed cost refers to the cost of fixed equipment, including sales personnel salary, cost of vehicle maintenance, enterprise management, and vehicle depreciation. The calculation of total fixed cost is as shown in the formula:
C f = j V + k M f k x p j k .
(2) C o p (operating cost)
The operating cost is related to the total operating time and other factors. The total operating time cost is calculated as shown in the formula:
C o p = i V j V k M x i j k c i j .
(3) C p w t (passenger waiting cost)
Passenger waiting cost refers to the waiting time cost caused by time window limitations when passengers are waiting for follow-up passengers to get on the vehicle, which focus on the perspective of passengers. The total cost of passenger waiting time is calculated as follows:
C p w t = k K i V j V c p x i j k y i k W j k .
(4) C v w t (vehicle waiting cost)
Vehicle waiting time cost refers to the time cost caused by time window limitations, which focus on the perspective of vehicles. The total cost of vehicle waiting time is calculated as shown in the formula:
C p w t = k K i V c w k W i k .
(5) Service quality penalty
C q u a = C p i c k u p u n s a t i s + C d e l i v e r y u n s a t i s + C d t
In order to measure the service quality, we define three quality indexes. Research shows that most passengers prefer the vehicle to arrive at the mid-point of the pickup time window and the low limit of the delivery time window [29]. C p i c k u p u n s a t i s indicates the time gap between the expected pickup time and the actual ones, C d e l i v e r y u n s a t i s indicates the time gap between the expected delivery time and the actual ones. C d t is the difference between direct travel time and actual travel time which caused by necessary detouring distance.
C p i c k u p u n s a t i s = i N k K | e i + + l i + 2 D i k | · i V x p i k ,
C d e l i v e r y u n s a t i s = i N k K | e i D i k | · i V x p i k ,
C d t = i N k K | t s e r v e i p t i , i + n | · i V x p i k .
Equations (2)–(4) is the equilibrium constraint. Equation (5) ensures each section of every service unit is assigned to exactly one vehicle. Equations (6) and (7) ensures only if the vehicle passes certain section of the service unit, can the demand be served. Equation (8) ensures the continuity of service with transfers. Equation (9) ensures the starting point of every vehicle is corresponded to the starting station it belongs to. Equation (10) ensures each vehicle only goes back to one of the returning stations. Equation (11) ensures the total number of vehicles returned to every station is higher than the minimum value. Equation (12) ensures vehicles will not visit a same transfer point repeatedly in one service procedure. Equations (13)–(16) ensures the condition for the occurrence of transferring: only if vehicle k and vehicle l passes transfer point t in a same service procedure, can they exchange passengers in transfer point t, and only if demand n has been picked up by vehicle k, can it be transferred from vehicle k to vehicle l. Equations (17) and (18) defines the arrival time of each node. Equation (19) ensures the transfer waiting time of passengers is shorter than   t w p . Equation (20) ensures the transfer waiting time of vehicles is shorter than   t w c . Equations (21) and (22) ensures if   w t i k l = 1 ,   D i + k < D t k and   D t l < D i k , which defined the time sequence of the procedure of transferring. Equations (23) and (24) furtherly define the sequence of arrival. Equations (25) and (26) is car loading constraint which define how vehicle pickup and deliver customers when without passing transfer points. Equations (27) and (28) define how vehicle pickup and deliver customers when it comes to transfer. Equations (29) and (30) define the maximum operating time and maximum travelling distance of each vehicle.

3.3. Operational Mechanism

A five-step process which describes the mechanism of CSTILP system is designed, as shown in Figure 1.
Step 0: Request reservations are submitted before the operational process of vehicles starts. The information each reservation provides include the number of passengers, the pickup and delivery locations and the starting and ending time windows.
Step 1: Vehicles depart from their assigned starting depot vacantly, described by constraint (31).
Step 2: Vehicles pick up passengers. This could happen by passing the pickup location and pick passengers up directly which is described by constraints (6)~(7), or by accepting passengers from other vehicles at transfer locations, which was described by constraint (8).
Step 3: At the same time, vehicles deliver passengers into their destinations or transfer them to other vehicles. In the end, every picked-up passenger will be sent to their destinations, described by constraint (5).
Step 4: Vehicles go back to ending depots. The ending depot for each vehicle is not fixed, vehicles are assigned to ending depots according to the minimum returning number of vehicles for each depot, which is described by constraint (10).

4. ALNS-TS Algorithm

As accurate algorithm performs poorly when solving large-scale instances, introducing heuristic algorithm (hybrid ALNS-TS) is necessary.
Different from the standard ALNS-TS, hybrid ALNS-TS has the following characteristics:
(1) Destroying and repairing only two routes at a time. This is used to improve the performance of transfer operator. Transferring happens between only two routes (which constructed by two vehicles), thus it is easier for transferring to happen if the optimization process takes place on only two vehicles at a time, where the destroy operator and transfer operator are able to work corporately during whole optimization procedure. The advantage of this is that it allows pick-up and delivery sequence to be lined up in a way which could take place only if transferring is allowed. For example, if 10 passengers all boarded from the same place, and the destination of eight of them is the same (let us call it D 1 ), while the destination of the remaining two passengers is D 2 , which is very far from D 1 that makes it impossible for one vehicle to deliver passengers to their coordinate destinations without violating time windows. In this way, if the repair operator and the transfer operator did not work corporately in the optimization procedure, it would be very likely for the algorithm to assign another vehicle to serve this passenger, or refuse to provide service. In this way, second-best solutions without transferring would very likely be generated. Although the optimization process can further improve the solution, this unneeded “detouring” process reduce the efficiency of this algorithm by making the happening of transfers harder. Moreover, as transfer operator only work with the second half of repairing operator which inserts only delivery points instead of pick-up points, this second-best solution has to wait for another iteration to be fixed, which is likely not to happen due to the randomness of this algorithm.
However, by optimizing only two routes at a time, repairing operator, destroy operator and transfer operator could work corporately at the same time, which provides the maximum advantageous conditions for transferring to happen. For the previous example (10 passengers, the destination of 8 of them was D 1 , the destination of the rest two was D 2 ), if transfer operator was able to work together with repair operator in the optimization process, all the 10 passengers will be allowed to be picked up by a same vehicle. Then, this vehicle could transfer the two passengers with D 2 destination to other vehicles, maybe to vehicles which has passengers whose destinations have similar geographical location with the transferred ones. Considering that transferring is very helpful in serving instantaneous large passenger flows, we think it is necessary to put the happening of transferring in priority.
(2) Tabu list is introduced. The traditional ALNS is characterized by large amount of neighborhood solutions with low similarities, and it is unlikely for same neighborhood solutions be reached in different iterations. Thus, the risk of being sank in local optimum is of low possibility, thus seldom to be considered. However, as the improved ALNS algorithm shrinks the search domain, the risk of falling into local optimum cannot be neglected. Thus, the technique of tabu search is introduced in the optimization procedure.
(3) 2-opt destroy operators are introduced to enhance the possibility for the length of vehicle routes to change in the optimization process, in order to improve the performance of ALNS-TS algorithm.

4.1. Main Scheme of Hybrid ALNS-TS

4.1.1. Generating Initial Solution

The Clarke and Wright savings algorithm (C-W algorithm) was implemented to generate initial solution [30]. The initial solution was obtained by inserting pickup and delivery points into routes, and accepting the insertion option with maximum savings. The initial solution is feasible and without any transfer point.

4.1.2. ALNS-TS Algorithm

The fundamental optimization technique of ALNS-TS is to start from initial solution and search its adjacent domain for satisfactory solutions, where ALNS operators ( o p i + , o p i ), TS (tabu) operators, and Transfer operator ( T r s ) work in collaboration. ALNS operators can be characterized into destroy ( o p i ) and repair ( o p i + ) operators, which work together to obtain adjacent solution by destroying and repairing the current solution iteratively in order to improve it.   T r s works together with repair operator in order to add transfer points in route. T S operator delimits the searching procedure to avoid trapping into local optimization, rejecting moves to points already visited in the search space.
Based on the proposed ALNS-TS, the framework is depicted as follows.
1: Generating   I S by C-W algorithm
2: T S _ A L N S _ A l g o r i t h m ( I S ,   N )
3: B S = I S ,   C S = I S
4: N e i g C S = Null
5:  while the termination criterion is not satisfied
6:   Repeat:
7:     Select   destroy   operator   o p i   and   repair   operator   o p i   according   to   probabilities   p i , p i +
8:     C S = ( o p i ( C S ) ) ,
9:     C S = ( o p i + , T r s ( C S ) ) ,
10:     Add   C S   to   Neig C S
11:   Until: Termination criterion is satisfied
12:   Repeat:
13:     if   C S j is in tabu list
14:      if   f ( C S j ) <   f ( B S ) ,   B S = C S j ,   C S n e w = C S j
15:    else
16:      if   f ( C S j ) <   f ( B N S ) ,   B N S = C S j
17:    end
18:   Until: all neighborhood solutions have been searched.
19:    if   f ( B N S ) < f ( B S ) ,   B S = B N S , C S n e w = B N S
20:    C S = C S n e w
21:    Update     p i + , p i   according to adaptive Mechanism
22:   Update tabu list
23:  end while
24: return   B S , C S ,   p i + , p i
25:end
Where C S : Current-solution, C S : Neighborhood-solution ,   N e i g C S : Set of Neighborhood solutions,   I S : Initial-solution , B S : Best-solution, B N S : Best-solution in this neighborhood.

4.2. ALNS Components

ALNS iteratively destroys and repairs current solution in order to improve it. In this way, the components of a complete ALNS algorithm can be categorized into destroy operators ( o p i ) which remove request from current solution and repair operators ( o p i + ) which reinsert request back. In each iteration of improving current solution, there are   N e i N u m neighborhood solutions need to be generated. To generate each neighborhood solution, only one destroy operator and one repair operator is chosen and used to optimize the selected route.
To shrink the computational labor and avoid trapping into local optimization, we introduced modified version of standard ALNS operators and Tabu (TS) operator into this algorithm. Different from the standard version of ALNS operators, the modified version generates neighborhood solution by optimizing two routes at a time instead of optimizing all the given routes.

4.2.1. Destroy Operators

The destroy operator set   ( O P { o p 1 ,   o p 2 , ,   o p 5 } )   contains five operators, each of them belongs to one of these two subsets:   O P s w a p and   O P 2 o p t , O P = O P s w a p O P 2 o p t . One of the destroy operators in   O P is selected each time when generating one neighboring solution. The possibility of   o p i to be selected is determined by the corresponding possibility   p i .
(1) Swap Operator
The general working procedure of swap operator is as follows.
1: Procedure   of   removal   ( C S   , o p i ) {
2:    Randomly   choose   two   routes :   R a   and   R b .
3:    R a = R a ,   R b = R b ,   C S = C S
4:   Calculating the temporary variables.
5:    Selecting   the   requests   to   be   removed   according   to   the   removing   rule   of   o p i and the corresponding temporary variables.
6:    Remove   the   chosen   request   from   R a   and   R b
7:    Define   the   request   set   removed   from   R a   as   S e t a r ,   and   which   removed   from   R b   as   S e t b r .
8:    Return   C S   ,   R a ,   R b , S e t a r ,   S e t b r
9:}
The general framework of each operator is the same, the difference is displayed in the procedure of generating temporary variables and the rules they follow when removing requests. The corresponding rule to be followed of each operator is as follows.
I. Random Removal
Randomly remove   p r % of the requests from route a and b.
II. Worst Removal
The cost-savings by the removal of each request in route a and b are computed, then p w % of the requests are selected to be removed according to the probability which increases with the cost-savings of their removal.
III. Related Removal
This heuristic operates basing on the belief that “similar requests are more likely to be exchanged between routes and might lead to better solution.” We define the relatedness between request i and request j according to their attributes: location, starting time, end time, quantity of their customers, etc. According to the paper written by Ropke and Pisinger [19], We define the relatedness as
R L ( i , j ) = φ d i j + ζ | D i k D j l | + κ ( 1 | K i K j | m i n { | K i | , | K j | } )
The operating process of this operator is as follows:
  • Randomly chooses the first request to be removed, named   n a i .
  • Calculate the relatedness between   n a i and every request in route b.
  • List them in descending order according to the relatedness measure to form a list   L .
  • Choose requests to be removed. For each request n 0 , its possibility of being chosen is R L ( i , n 0 ) n N l R L ( i , n ) .
(2) 2-opt operator
The general working procedure of 2-opt operator is as follows.
1: Destroying   ( C S   , o p i ) {
2:    Randomly   choose   two   routes :   R a   and   R b .
3:    R a = R a ,   R b = R b , C S = C S
4:    Select   the   chosen   sec tion   R a [ i ,   j ]   and   R b [ k ,   l ] to be removed according to the corresponding rules of each operator. R a [ i ,   j ]   is   the   sec tion   between   position   i   and   position   j   in R a ,   i ,   j < l e n g t h ( R a ) ) ,   sec tion   R b   is   the   sec tion   between   position   k   and   position   l   in   R b ,   k ,   l < l e n g t h ( R b ) )
5:    Extract   R a [ i ,   j ] from R a and R b [ k ,   l ] from R b .
6:    Define   the   request   set   removed   from   R a   as   S e t a r   ,   and   which   removed   from   R b   as   S e t b r .
7:    Return C S , R a ,     R b , S e t a r ,   S e t b r ,   R a [ i ,   j ] ,   R b [ k ,   l ] ,   i ,   j ,   k ,   l
8:}
The general framework of each operator is the same, the difference between them is displayed in the procedure of the rule of choosing sections to be removed. The corresponding rule to be followed of each operator is as follows.
I. Random Removal
Randomly pick one piece of section from each route.
II. Best Removal
Generate two random numbers   ( r ,   t ) . Calculate the average cost saving of route   R a if a section with a length of   r is deleted iteratively until all possibilities have been calculated. Calculate the average cost saving of route   R b if a section with a length of   t is deleted iteratively until all possibilities have been calculated. Choose the option with maximum saving and delete the chosen section from the original chains.

4.2.2. Repair Operator

The repair operator set   ( O P + { o p 1 + ,   o p 2 + ,   o p 3 + } )   contains three operators. One of the repair operators in   O P + is selected each time when generating one neighboring solution. The possibility of   o p i + to be selected is determined by the corresponding possibility   p i + .
The general working procedure of repair operator o p k + is as follows.
1: R e p a i r i n g   ( Route   R a ,   R b ,   unoccupied   request   set   U R ,   S e t a r ,   S e t b r ) {
2:   Define   U R B = U R S e t a r ,   U R A = U R S e t b r
3:  for R a , R b
4:  do
5:    pick   one   request   R from   set   U R A ( R B ) following   R u l k
6:   do
7:     look   for   a   position   where   R can be inserted in Ra(Rb) according to the repairing rule of o p i +
8:    if at least one feasible solution exists
9:      accept   this   insertion   and   renew   R a ( R b ).
10:     URA = URAR, URB = URBR.
11:     break
12:    else
13:     continue
14:    until both chains have been processed.
15:    if there is not any feasible solution
16:     continue
17:   until all unoccupied requests have been processed.
18:    set C S ( R a ) = R a ,   C S ( R b ) = R b
19:    set   C S =   set   C S
20:   UR = URAURB
21:    return   set   C S ,   U R
22:}
The general framework of each operator is the same, the difference is displayed in the procedure of the rule of choosing and inserting requests. The corresponding rule to be followed of each operator is as follows.
I. Random Insertion
Random insertion heuristic is a serial computational process. At each iteration, unoccupied requests are randomly selected and randomly inserted in random positions of selected vehicle chain.
II. Greedy Insertion
Define the cost increase after the insertion of request i into route r as
c i , r a d d = c r o r i c i , r i n s
Define Cg as the minimum cost of the insertion of request i:
c g = m i n r C S c i , r a d d
where g is the best inserting position.
This operator works by iteratively assigns requests to the batch according to best inserting position g. Note that this operator also ensures removed requests not being inserted into the original route.
III. Regret-2 insertion
The operating process of this operator is as follows:
  • Define   α k i as the number k best insertion option for request i.
  • Define β 2 i as the difference between the optimal inserting solution and the suboptimal inserting solution.   β 2 i = c α 2 i c α 1 i
  • Choose the operating request   i o p by the following expression:
    i o p = arg m a x i N β 2 i
  • Iteratively assigns requests by descending order to the batch according to best inserting position. Note that this operator also ensures removed requests not being inserted into the original route.

4.2.3. Transfer ( T r s ) Operator

In the procedure of generating neighboring solutions, this is how   o p i + and   o p i work after they are chosen:
C S d e s t o r y i n g   c u r r e n t   s o l u t i o n     o p i , T r s   o p e r a t o r C S i n s e r t i n g   p i c k u p   p o i n t s                           o p i + ( 1 ) C S i n s e r t i n g   d e l i v e r y   p o i n t s o p i + ( 2 ) , T r s   o p e r a t o r C S ,
where   o p i deletes both the pickup section and the delivery section from vehicle routes. T s (tabu) operator determines whether and which transfer point will be removed.   o p i + ( 1 ) inserts the pickup section into the current solution. T r s operator works together with   o p i + ( 2 ) to generate solutions with different numbers of transfer points. In this way, after the process of repairing, two neighboring solutions with exactly same sequence of pickup points are generated, one with transfer points and one without, and only the better ones will be accepted.
Define the maximum number of transfer points in one route as   n t m . In this paper, the value of   n t m equals to 1.
The pseudocode of how this procedure works is as follows.
1: G e n e r a t i n g   q t h   C S ( o p i , o p i + ( 1 ) , o p i ( 2 ) , T r s   o p e r a t o r ,   R a , R b ) {
DELETIING   TS   POINTS
2:   if (number of Ts points in route   C S ( R x )   equals to or is larger than   n t m ,):
3:    Define the route which has transferring relationship with   R x as:   R x c o r r e s
4:    Replace R x with   Set S o l u t i o n N o T s ( x ) and   R x c o r r e s with   Set S o l u t i o n N o T s ( x c o r r e s )
5:   end
6:   optimizing   R x :   R x d e s t o r y i n g   n e i g h b o r h o o d o p i ,   T r s   o p e r a t o r R x i n s e r t i n g   p i c k u p   p o i n t s o p i + ( 1 ) R x
GENERATING   TSBESTSOLUTION
7:   Repeat:
8:    Repeat:
9:     judge whether transfer point   t can be inserted in   R x ( i ) .
10:     if insertion ( i , t , R x ) is feasible
11:       a d d = [ i ,   t ,   E a r l i e s t A r r i v a l T i m e ,   L a s t e s t L e a v i n g T i m e ] ,
12:       S e t T s ( R x ) =   S e t T s ( R x ) + a d d
13:     end
14:    Until all Transfer points have been checked
15:   Until all possible positions in route R x have been checked
16:   Determine where can transfer point be inserted in order to construct transferring relationship between   R a and   R b .
Describe the inserting option as:
F e a s i b l e T s = [ i , j ,   t ,   E a r l i e s t A r r i v a l T i m e ,   L a s t e s t L e a v i n g T i m e ] , where i denote the inserting position of transfer point t in R a , j denote the inserting position of transfer point t in R b .   S e t T s = S e t T s + F e a s i b l e T s
17:   Repeat:
18:    Determine which of the request in each route can be transferred by their time windows. Develop sets   Set a t s and   Set b t s . List all the combinations of request from two groups in set   T s O p t i o n s .
19:    Repeat:
20:     Judge whether this option is feasible. If it is, replace the original   R a , R a c o r r e s ,   R b , R b c o r r e s with the new ones. calculate the total cost.
21:    Until all options of transferring in   T s O p t i o n s have been calculated.
22:   Until all options of inserting transfer point in   S e t T s have been calculated.
23:   Pick the ones with minimum cost as   B e s t T s S o l u t i o n
24: GENERATING   NOTSSOLUTION
25:   Use the traversal methods which were discussed in 4.2.2 to repair R a and R b .
Replace the original R a and R b in   C S with R a and R b . Name the new solution as   N o T s S o l u t i o n and calculate the total cost.
26:    Set S o l u t i o n N o T s ( q ) = N o T s S o l u t i o n .
27:   if   Cos t N o T s S o l u t i o n > Cos t B e s t T s S o l u t i o n :
28:     Neig C S = Neig C S + N o T s S o l u t i o n .
29:   else
30:     Neig C S = Neig C S + B e s t T s S o l u t i o n .
33:end
34:}
Where, Set S o l u t i o n N o T s :a set of solutions with no transfer points which was generated through the optimizing process. Set S o l u t i o n N o T s is restored for usage, routes X in this set will be used to replace the current route when corresponding route X in current solution with transfers is being operated in another iteration. N S T :   N o T s S o l u t i o n .

4.3. Adaptive Mechanism

Adaptive mechanism records the performance   α and usage count θ of each destroy or repair operator and adjust their possibilities to be chosen according to their performance in past iterations.
The detailed mechanism is as follows.
Before the searching process starts, all destroy or repair operators are set to the same initial weight   ω . The ALNS searching process is divided into several segments, each segment contains   λ iterations. After   λ iterations, weights are updated according to the quality of solutions the corresponding operator reached, which are measured by four scenarios   ( α 1 > α 2 , α 3 ) .
In every iteration:
1.
β i = β i + α 1 , if a new best solution is reached by   o p i in this iteration.
2.
β i = β i + α 2 , if the solution   o p i reached in this iteration is better than the current solution.
3.
β i = β i + α 3 , if the solution   o p i reached in this iteration is worse than the current solution but is accepted.
4.
β i = β i , else.
After   λ iterations, renew corresponding weight of each operator:
ω i = { ω i ( 1 r ) + r β i θ i ,     θ i > 0 ω i ( 1 r ) ,   θ i = 0 .
For all the destroy or repair operator to be chosen, their possibility to be chosen is   ω i i ω i .

5. Numerical Experiment

In the numerical experiment, we proved the validity of the proposed ALNS-TS algorithm and the CSTILP model by comparing the performance of the ALNS-TS algorithm with that of GAMS (A commonly used commercial solver), based on a set of cases with incremental scale. Besides, the comparison can further prove efficiency of the proposed algorithm on relatively large-scale cases.
After that, we evaluated the performance of ALNS-TS algorithm on six different cases derived from realistic demands. These six cases differ in their temporal and spatial homogeneity and can be used to verify the adaptation of the proposed method in different scenarios.

5.1. Validation of the ALNS-TS Method

We developed instances with 4–10 service units which were solved by both GAMS and proposed ALNS-TS algorithm. The solving result GAMS presented were regarded as the standard solution, and the feasibility and efficiency of ALNS-TS algorithm is proven by comparing the solutions and operation time of these two kinds of solutions.
Figure 2 presents a comparison of the GAMS solutions and ALNS-TS solutions from the perspectives of the operation time. The vehicle routing solutions and operation time are also shown in Table 4, which shows that ALNS-TS algorithm obtains the same results as GAMS does; thus, the proposed method is feasible.
To prove the effectiveness of this algorithm, we discuss the operation time of these two methods. As the number of service units increases, the calculation times for GAMS grows exponentially, as is shown in Figure 2. When there are four service units, the calculation time is 4.55 s, but when the number grows from 4 to 5, the operation time increases by about 48 times to 168.20 s. The operation time reaches 294,506.23 s (approximately 82 h) when the number of service units reaches 10, which is regarded as the limit, that is at least 24,000 times that of the time cost by ALNS-TS algorithm.
Although the calculation time of ALNS-TS generally grows when the scale of service units gets larger, it shows a decrease when the number of service unit grows from 4 to 5, then shows an increase when the number of service unit reaches 7. This results from the fact that when number of service units decreases, it becomes harder for ALNS-TS to reach new feasible solutions, in this way the algorithm spends more time in searching for better solutions.
These observations demonstrate the effectiveness of the proposed method.

5.2. ALNS Computational Experiments

The method was coded in MATLAB and experiments were run on an E74830-CPU computer operated by Windows Server 2016. We discussed ALNS settings and used data based on real-life to evaluate the benefits of routing with transfer points.

5.2.1. Cases Design

To discuss the applicability of this algorithm, we developed six cases based on Beijing’s layout to carry out experimentation. As train stations and airports are important passenger distributing centers, the passenger flow of the instances used in this paper mainly arose from sources of stations, airports and residential area.
We developed six cases which were grouped into three case pairs with different spatiotemporal distribution characteristics to examine the spatiotemporal applicability of this algorithm. Besides, one of the two cases in a same case pair is with transfer option, the other is without transfer option, in order to examine whether improvement will be made and how much improvements will be achieved after transferring are implemented in operation.
The flow-generation points include: Beijingxi Railway Station, Beijingbei Railway Station, Beijingnan Railway Station, Beijingchaoyang Railway Station, Beijing Railway Station, airports: Beijing Daxing International Airport, Beijing Capital International Airport, and 13 different kinds of trip generation centers in residentials areas, known as housing estates attachments, other housing estates, community centers, dormitories, residential estates, universities, hotel, resort, apartment hotel, economy hotel, hotel affiliations, guest house, other hotels.
A request   n is composed of OD (origin-destination) pair with time windows. In the first four instances, we took airports or stations as request origin or destination, and took trip generation centers in residential areas as the other end of OD pair, to generate the spatial distribution of requests. After that, we used real-time timetable of airports and train stations to get the OD time windows of airports/stations end, and used distances between origin and destination to generate time windows of the other end. We divided passenger’s demand into 34 categories, including travelling from different airports or stations to residential area, travelling from residential area to different airports or stations, and travelling between different airports and stations. The proportion of each kind of request was decided by the actual number of passengers arrive with flights and trains. In the fifth and sixth instance, we used trip generation centers in residential areas to generate both origin and destination, whose one end of time windows were the same as the first and second instance, the other end was generated based on the travelling distance divided by traveling speed.
In the first two instances, the temporal distribution of request was comparatively even, which was based on real-time timetable from 23:00 to 7:00 the next day. However, the spatial distribution of which was centralized, where one end of the tour was airports or train stations, and the other end was randomly picked in residential area. In the second instance, transferring was allowed, while in first transferring was not allowed. The spatial distribution of the third and fourth instances was the same as the first and the second one, but in temporal distribution, service requests were distributed more densely from 5:00 to 7:00 to form super-peak. In the fourth instance, transferring was allowed, while in the third transferring was not allowed. The spatiotemporal distribution of both the fifth and sixth instance were scattered, with time windows generated randomly from 23:00 to 7:00 the next day, and requests picked up randomly from Beijing’s residential area. In the sixth instance, transferring was allowed, while in the fifth transferring was not allowed.
The settings of these six instances are shown in Table 5. Vehicles with different capacity and different fixed cost are assigned to each depot before the operational process starts.

5.2.2. Parameter Settings

For most of the parameters of the ALNS algorithm, we validate the values proposed by Ropke and Pisinger [19]. The parameter description and settings are shown in Table 6.

5.3. Research Results

The operational results of these six instances are shown in Table 7. The improvements in quality index of different case pairs after transferring was introduced are shown in Table 8. Figure 3 visualized the improvements in quality index after transferring is introduced by a bar chart.
It could be concluded from the experimental result that fixed cost reduced by 4.51% and the total cost reduced by 6.84% when transferring became allowed in cases with centralized passenger flow, and the number was 9.90% and 9.93% in super-centralized ones, while in cases with scattered passenger flow the fixed cost did not change and total cost reduced by 4.11%. This shows transferring can to a large extent reduce the operating cost of CB system, especially when the spatiotemporal distribution of passengers is centralized. Comparing the case with transfer option and the case without transfer option in all three case pairs, the travelling time of one passenger on average reduced after transferring is allowed, where the case with centralized passenger flow reduced by 7.84%, the case with super-centralized flow reduced by 9.95%, which was very satisfying, and the one with scattered passenger flow reduced by 5.80%. Both passenger waiting time and vehicle waiting time reduced when transferring became allowed in all three instance pairs, but the improvement appeared to be more significant only when the concentration degree got relatively high. In centralized cases numbers of passenger waiting time and vehicle waiting time were 0.05% and 18.69%, which were even less than the number in scattered case, that were 15.72% and 36.94%. Improvement in passenger waiting time and vehicle waiting time reached 33.05% and 45.27% when concentration ratio continued to increase. The number of vehicles in operation reduced by 9.40% after transferring became allowed in super-centralized flow case and reduced by 5.66% in centralized flow case, whereas in scattered ones the number stayed still. This might because when the requests are scattered, even if transferring could integrate some requests in other routes, it might be difficult to put all of the requests especially the ones with remote origins or destinations into other routes as well. Table 9 shows the details of the optimized solutions.
It could be concluded from the comparison of these three case pairs that transferring could enhance the performance of ALNS-TS when spatiotemporal distribution of requests is centralized. It could increase the indexes which indicate passengers’ satisfaction and reduce the operating cost at the same time. However, the enhancement brought by transferring is not significant in when dealing with instances with scattered spatiotemporal distribution of requests.

6. Conclusions

Aiming on making customized bus an option in dealing with instantaneous passenger flow, this paper proposed a CSTILP model to satisfy the traffic demands of these passengers. Modes such as multi-depots, multi-vehicle types and transferring are introduced in this model, which help reach a satisfying balance between passenger satisfaction and operating cost.
When the scale of problem grows, the combination explosion problem emerges, which limits the solving capacity of commercial solvers to a relatively small scale. In order to provide feasible solutions of large-scale problems, this paper proposed a matching algorithm which combines traditional ALNS and TS problem into a hybrid one according to the uniqueness of this problem.
In order to prove the feasibility and efficiency of the algorithm, seven small-scale instances and six large-scale ones were proposed. We proved the feasibility of the algorithm by comparing the one solved by GAMS and by our proposed method, and after that we ran six large-scale experiments.
Experiments show: (1) the proposed algorithm has the same accuracy as that of GAMS, a commercial solution software, but has a higher speed: when the demand size is set as 10, the proposed algorithm can save 24,000 times of time; (2) by comparing results of six reality-based cases, we found that the designed option could save 9.93% times of fleet cost and reduce 33.05% times of passenger waiting time relative to other existing customized bus modes when encountering instantaneous passenger flows with time and space imbalance.
Then we can conclude that (1) in dealing with cases of instantaneous super-peak passenger demand, comparing our proposed method with which without transferring, ours can at the same time shrink the operational cost and improve the service quality index from the perspective of both passengers and drivers; (2) our proposed method specializes in satisfying instantaneous passenger flow, it is proved to have better performance in cases where there exists group of people making transportation requests to similar geographical area; (3) our proposed algorithm could not only to a large extent shrink the operational time, but also reach relative optimal solution.
In this study, we only considered providing service for instantaneous large passenger flow with centralized spatiotemporal distribution, and took customized shuttle bus as the only transport means. Therefore, we will seek to expand the scope of application in future research, and integrate more means of transportation in our operational mode by cooperating our service mode with other public transports.

Author Contributions

Conceptualization, H.Z.; data curation, J.C.; formal analysis, J.C.; funding acquisition, J.C.; supervision, X.Z.; writing—original draft, H.Z.; writing—review & editing, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the National Science Foundation of China (XZ) (Grant No. U1734204), the Science and Technology Department of the China Railway Corporation (XZ) (Grant No. N2020x022) and Fundamental Research Funds for the Central Universities (Grant No. 2020JBZD006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Masson, R.; Ropke, S.; Lehuédé, F.; Péton, O. A branch-and-cut-and-price approach for the pickup and delivery problem with shuttle routes. Eur. J. Oper. Res. 2014, 236, 849–862. [Google Scholar] [CrossRef]
  2. Mitrović-Minić, S.; Laporte, G. The pickup and delivery problem with time windows and transshipment. INFOR 2006, 44, 217–227. [Google Scholar] [CrossRef]
  3. Deleplanque, S.; Quilliot, A. Dial-a-ride problem with time windows, transshipments, and dynamic transfer points. IFAC Proc. 2013, 46, 1256–1261. [Google Scholar] [CrossRef]
  4. Sol, M. Column generation techniques for pickup and delivery problems. Ph.D. Dissertation, Eindhoven University of Technology, Eindhoven, The Netherlands, 8 November 1994. [Google Scholar]
  5. Naccache, S.; Côté, J.; Coelho, L.C. The multi-pickup and delivery problem with time windows. Eur. J. Oper. Res. 2018, 296, 353–362. [Google Scholar] [CrossRef]
  6. Parragh, S.N.; Schmid, V. Hybrid column generation and large neighborhood search for the dial-a-ride problem. Comput. Oper. Res. 2013, 40, 490–497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ropke, S.; Cordeau, J.F.; Laporte, G. Models and branch-and-cut algorithms for pickup and delivery problems with time windows. Networks 2010, 49, 258–272. [Google Scholar] [CrossRef]
  8. Tan, K.C.; Lee, T.H.; Ou, K.; Lee, L.H. A messy genetic algorithm for the vehicle routing problem with time window constraints. In Proceedings of the 2001 Congress on Evolutionary Computation, Seoul, Korea, 27–30 May 2001; pp. 679–686. [Google Scholar]
  9. Wester; Dawn, V. A Genetic Algorithm for the Vehicle Routing Problem. Available online: https://trace.tennessee.edu/cgi/viewcontent.cgi?article=6358&context=utk_gradthes (accessed on 24 September 2021).
  10. Ombuki, B.; Ross, B.J.; Hanshar, F. Multi-objective genetic algorithms for vehicle routing problem with time windows. Appl. Intell. 2006, 24, 17–30. [Google Scholar] [CrossRef]
  11. Cordeau, J.F.; Laporte, G. A tabu search heuristic for the static multi-vehicle dial-a-ride problem. Transp. Res. Part B Methodol. 2003, 37, 579–594. [Google Scholar] [CrossRef] [Green Version]
  12. Brandão, J. A tabu search algorithm for the open vehicle routing problem-sciencedirect. Eur. J. Oper. Res. 2004, 157, 552–564. [Google Scholar] [CrossRef]
  13. Lai, D.S.; Demirag, O.C.; Leung, J.M. A tabu search heuristic for the heterogeneous vehicle routing problem on a multigraph. Transp. Res. E Logist. Transp. Rev. 2016, 86, 32–52. [Google Scholar] [CrossRef]
  14. Taillard, P.E. A parallel tabu search heuristic for the vehicle routing problem with time windows. Transp. Res. Part C Emerg. Technol. 1997, 5, 109–122. [Google Scholar]
  15. Gambardella, L.M.; Dorigo, M. Solving symmetric and asymmetric tsp by ant colonies. In Proceedings of the IEEE International Conference on Evolutionary Computation, Nagoya, Japan, 20–22 May 1996. [Google Scholar]
  16. Donati, A.V.; Montemanni, R.; Casagrande, N.; Rizzoli, A.E.; Gambardella, L.M. Time dependent vehicle routing problem with a multi ant colony system. Eur. J. Oper. Res. 2008, 185, 1174–1191. [Google Scholar] [CrossRef]
  17. Li, Y.; Soleimani, H.; Zohal, M. An improved ant colony optimization algorithm for the multi-depot green vehicle routing problem with multiple objectives. J. Cleaner Prod. 2019, 227, 1161–1172. [Google Scholar] [CrossRef]
  18. Zhang, H.; Zhang, Q.; Ma, L.; Zhang, Z.; Liu, Y. A hybrid ant colony optimization algorithm for a multi-objective vehicle routing problem with flexible time windows. Inf. Sci. 2019, 490, 166–190. [Google Scholar] [CrossRef]
  19. Ropke, S.; Pisinger, D. An adaptive large neighborhood search heuristic for the pickup and delivery problem with time windows. Transp. Sci. 2006, 40, 455–472. [Google Scholar] [CrossRef]
  20. Ghilas, V.; Demir, E.; Woensel, T.V. An adaptive large neighborhood search heuristic for the pickup and delivery problem with time windows and scheduled lines. Comput. Oper. Res. 2016, 12–30. [Google Scholar] [CrossRef] [Green Version]
  21. Masson, R.; Lehuede, F.; Peton, O. An adaptive large neighborhood search for the pickup and delivery problem with transfers. Transp. Sci. 2013, 47, 344–355. [Google Scholar] [CrossRef] [Green Version]
  22. Ritzinger, U.; Puchinger, J.; Hartl, R.F. Dynamic programming based metaheuristics for the dial-a-ride problem. Ann. Oper. Res. 2016, 236, 341–358. [Google Scholar] [CrossRef] [Green Version]
  23. Bachem, A.; Malich, M. The simulated trading heuristic for solving vehicle routing problems. Discrete Appl. Math. 1996, 65, 47–72. [Google Scholar] [CrossRef] [Green Version]
  24. Yu, V.F.; Lin, S.W.; Lee, W.; Ting, C.J. A simulated annealing heuristic for the capacitated location routing problem. Comput. Ind. Eng. 2010, 58, 288–299. [Google Scholar] [CrossRef]
  25. Chao, W.; Dong, M.; Fu, Z.; Sutherland, J.W. A parallel simulated annealing method for the vehicle routing problem with simultaneous pickup–delivery and time windows. Comput. Ind. Eng. 2015, 83, 111–122. [Google Scholar]
  26. Chiang, W.C.; Russell, R.A. Simulated annealing metaheuristics for the vehicle routing problem with time windows. Ann. Oper. Res. 1996, 63, 3–27. [Google Scholar] [CrossRef]
  27. Kachitvichyanukul, V.; Sombuntham, P.; Kunnapapdeelert, S. Two solution representations for solving multi-depot vehicle routing problem with multiple pickup and delivery requests via pso. Comput. Ind. Eng. 2015, 89, 125–136. [Google Scholar] [CrossRef]
  28. Dridi, I.H.A.; Alaa, E.B.; Borne, P.; Bouchriha, H. Optimisation of the multi-depots pick-up and delivery problems with time windows and multi-vehicles using pso algorithm. Int. J. Prod. Res. 2020, 58, 4201–4214. [Google Scholar] [CrossRef]
  29. Paquette, J.; Cordeau, J.F.; Laporte, G.; Pascoal, M.M. Combining multicriteria analysis and tabu search for dial-a-ride problems. Transp. Res. B Meth. 2013, 52, 1–16. [Google Scholar]
  30. Clarke, G.; Wright, J.W. Scheduling of vehicles from a central depot to a number of delivery points. Oper. Res. 1964, 12, 568–581. [Google Scholar] [CrossRef]
Figure 1. Illustration of the CSTILP design model. This figure demonstrates the main optimization content of the proposed model. There are four types of nodes in this figure, where the cycles represent origins of requests, the squares represent destination of requests, the rounded squares with car icon represents depots, and the bigger rounded square represents the transfer points. The request includes information about location and time window. The figure shows the designed operational process of CSTILP design model.
Figure 1. Illustration of the CSTILP design model. This figure demonstrates the main optimization content of the proposed model. There are four types of nodes in this figure, where the cycles represent origins of requests, the squares represent destination of requests, the rounded squares with car icon represents depots, and the bigger rounded square represents the transfer points. The request includes information about location and time window. The figure shows the designed operational process of CSTILP design model.
Information 12 00429 g001
Figure 2. Comparison of the operation time of two methods.
Figure 2. Comparison of the operation time of two methods.
Information 12 00429 g002
Figure 3. Improvements in indexes after transferring is introduced.
Figure 3. Improvements in indexes after transferring is introduced.
Information 12 00429 g003
Table 1. Names and Description of Variables.
Table 1. Names and Description of Variables.
Names of VariablesDescription
x i j k { 1 ,   if   vehicle   k   travels   from   location   i   to   location   j   0 , otherwise , ( i , j ) V
z n ± k { 1 ,   if   the   pickup   section ( n + )   or   delivery   section   ( n )   is   assigned   to   vehicle   k 0 , otherwise   ,   k K ,   n N
w t n k l { 1 ,   if   vehicle   k   and   l   transfered   request   n   a t   transfer   point   t 0 , otherwise   , k , l K ,   n N , t T S
D i k depicts   the   departure   time   of   vehicle   k ,   i V , k K
y i k specifies   the   load   of   vehicle   k   when   it   arrives   at   location   i ,   i V , k K
Table 2. Names and Description of Sets.
Table 2. Names and Description of Sets.
Names of SetsDescription
N The   set   of   service   units ,   N N + N .
q n For   each   service   unit   n ,   there   are   q n   passengers   need   to   be   transported   from   position   n +   to   position   n ,   and   q n = q n + = q n .
N + N + { n + | n N } is defined as the set of pickup locations.
N N { n | n N } as the set of delivery locations.
D e The   set   of   depots ,   D e D e + D e .
D e The   set   of   starting   depots ,   D e { d e 1 ,   d e 2 ,   } .   Every   vehicle   has   a   fixed   starting   depot ,   vehicles   in   set   M α   have   p α as their starting depot.
D e + The   set   of   ending   depots , D e + { d e 1 + , d e 2 + ,   } .   Number   of   vehicles   returning   to   ending   depot   d e p +   need   to   be   bigger   than   the   minimum   number   k p .  
D e k + The   set   of   ending   depots   which   vehicle   k   can   go   back   to ,   D e k + { d e k + 1 , d e k + 2 ,   } .
T s The   set   of   transfer   points ,   T s { t 1 , t 2 t n } .
V The   locations   vehicles   can   visit ,   V N + N T s D e D e +
V N V N N + N denotes the service locations.
V D e + V D e + N + N T s   D e +
V D e V D e N + N T s   D e
V N o D e V N o D e N + N T s
K The   vehicle   set .   K { k 1 , k 2 k n }
Table 3. Names and Description of Parameters.
Table 3. Names and Description of Parameters.
Names of ParametersDescription
c i j Travel   cos t   of   arc   ( i ,   j )
D i j Travel   time   of   arc   ( i ,   j )
q i Number of passengers to be picked up or delivered
F k Fixed   cos t   of   vehicle   k
e i Earliest   arrival   time   of   request   i
l i Latest   arrival   time   of   request   i
Q k Maximum number of passengers vehicle k could load
p o s i + Pickup   location   of   request   i
p o s i Delivery   location   of   request   i
o p i Operation   time   of   location   i
k p Minimum   number   of   vehicles   to   be   arrived   at   depot   p
t w p The maximum amount of time a passenger can wait at transfer points
t w c The maximum amount of time a vehicle can wait at transfer points
OPiThe operation time at service point i
t s o p t The operation time at transfer point t
d m a x k The maximum operating distance of vehicle k
t m a x k The maximum operating time of vehicle k
C f The total fixed cost
C o p The total operating cost
C p w t Passenger waiting cost
C v w t Vehicle waiting cost
C q u a Service quality penalty
C p i c k u p u n s a t i s The unsatisfaction caused by the time gap between the expected pickup time and the actual ones
C d e l i v e r y u n s a t i s The unsatisfaction caused by the time gap between expected delivery time and the actual ones
C d t The unsatisfaction caused by the difference between direct travel time and actual travel time
f k The fixed cost of vehicle k
c p The cost per unit time of passenger waiting
c w k The cost per unit of waiting time of vehicle k
W j k The waiting time of vehicle k at point j
Table 4. Comparison of the operation time of two methods.
Table 4. Comparison of the operation time of two methods.
Case IDNumber of
Service Units
Operating Time of Solving Method/SecVehicle Routing Solution
GAMSOur Proposed Method
144.5580.42 r o u t e 1 :   d e 1 - n 2 + - n 2 - n 1 + - n 1 - n 3 + - n 4 + - n 3 - n 4 d e 1
25168.245.42 r o u t e 1 :   d e 1 - n 1 + - n 1 - n 2 + - n 5 + - n 5 - n 2 - n 4 + - n 3 + - n 4 - n 3 d e 1
363027.810.64 r o u t e 1 :   d e 1 - n 2 + - n 2 - n 5 + - n 5 - n 3 + - n 6 + - n 6 - n 3 - n 4 + - n 4 - n 1 + - n 1 d e 1
4749,402.269.91 r o u t e 1 :   d e 1 - n 4 + - n 3 + - n 4 - n 5 + - n 3 t s 1 - n 1 - n 5 - n 2 + - n 2 + - n 2 d e 1
r o u t e 2 :   d e 1 - n 1 + - n 7 + - n 6 + t s 1 - n 6 - n 7 d e 1
5895,373.5810.27 r o u t e 1 :   d e 1 - n 5 + - n 4 + - n 5 - n 6 + - n 4 - n 8 + - n 7 + - n 8 - n 7 - n 1 + - n 3 + t s 1 - n 6 d e 1 +
r o u t e 2 :   d e 1 - n 2 + - n 2 t s 1 - n 1 - n 3 d e 1 +
69194,795.1612.54 r o u t e 1 :   d e 1 - n 2 + - n 3 + - n 2 - n 7 + - n 8 + - n 3 - n 8 - n 7 - n 1 + - n 5 + t s 1 - n 1 - n 6 d e 1 +
r o u t e 2 :   d e 1 - n 9 + - n 9 - n 4 + - n 4 - n 6 + t s 1 - n 5 d e 1 +
710294,506.2312.25 r o u t e 1 :   d e 1 - n 1 + - n 1 - n 7 + - n 2 + - n 8 + - n 6 + - n 5 + - n 5 - n 7 - n 6 - n 2 - n 9 + - n 10 + - n 4 + - n 9 - n 10 - n 3 + - n 4 - n 3   d e 1
Table 5. Instance Settings.
Table 5. Instance Settings.
Case IndexSpatial DistributionTemporal DistributionSpatiotemporal DistributionTransferring Allowed or NotNumber of Requests
1CentralizedScattered, from 23:00 to 7:00CentralizedNot Allowed488
2CentralizedScattered, from 23:00 to 7:00CentralizedAllowed488
3CentralizedCentralized, from 5:00 to 7:00Super-centralizedNot Allowed538
4CentralizedCentralized, from 5:00 to 7:00Super-centralizedAllowed540
5ScatteredScattered, from 23:00 to 7:00ScatteredNot Allowed396
6ScatteredScattered, from 23:00 to 7:00ScatteredAllowed396
Table 6. Parameters in algorithm.
Table 6. Parameters in algorithm.
NotationDescriptionValue
ω i Initial weight of operator i0.1
σ 1 Operator evaluation parameter 133
σ 2 Operator evaluation parameter 29
σ 3 Operator evaluation parameter 313
r Operator weight adjustment parameter0.1
λ Length of iteration fragments100
i t e r a t i o n Max iteration2000
φ Similarity parameter 19
ζ Similarity parameter 23
κ Similarity parameter 32
t w p Passenger waiting time5 min
t w c Vehicle waiting time3 min
o p i Operation time at point i2 min
t s o p t Operation time at transfer point t3 min
d m a x k Maximum travelling distance300 km
t m a x k Maximum travelling time480 min
Table 7. Operational Results.
Table 7. Operational Results.
123456
Total cost81,418.7276,209.29 121,587.18111,011.5477,595.96 74,534.88
Vehicles in use168159326298153153
Passengers served488488538540396396
c f x 12,04011,52023,52021,48010,86010,860
c t o t a l p e r 166.84 156.17 226.00 205.58 195.95 188.22
c f x p e r 24.67 23.61 43.72 39.78 27.42 27.42
t t r a p e r (min)117.50 108.95 138.56 126.02 141.10133.37
t p w p e r (min)31.01 ± 25.11 31.00 ± 26.44 39.36 ± 27.34 29.58 ± 23.87 40.16 ± 29.92 34.71 ± 27.52
t v w p e r (min)126.78 ± 82.69 106.82 ± 74.60 71.49 ± 44.29 49.21 ± 47.81 109.66 ± 62.58 80.08 ± 60.77
t u p g a p (min)16.58 ± 25.80 15.25 ± 21.78 12.36 ± 11.33 11.82 ± 10.92 9.78 ± 8.49 10.33 ± 10.35
t d o w n g a p (min)17.60 ± 16.51 14.09 ± 13.28 24.10 ± 16.64 14.88 ± 16.76 25.17 ± 19.52 18.16 ± 16.00
t d t p e r (min)4.34 ± 15.31 7.18 ± 19.69 0.86 ± 5.08 7.06 ± 17.830.51 ± 1.32 7.10 ± 21.00
where: c f x : total fixed cost, c t o t a l p e r : total cost per passenger on average, c f x p e r : fixed cost per passenger on average, t u p g a p : difference between the expected picking up time and the actual picking up time of one passenger on average, t d o w n g a p : difference between the expected delivery time and the actual delivery time of one passenger on average, t t r a p e r : Travelling time of one passenger on average, t p w p e r : Waiting time of one passenger on average, t v w p e r : Waiting time of one vehicle on average, t d t p e r : detouring time of one passenger on average.
Table 8. Improvements in quality index after transferring is introduced.
Table 8. Improvements in quality index after transferring is introduced.
IndexCase 1 vs. Case 2Case 3 vs. Case 4Case 5 vs. Case 6
Vehicle in operation5.66%9.40%0.00%
c t o t a l p e r 6.84%9.93%4.11%
c f x p e r 4.51%9.90%0.00%
t t r a p e r 7.84%9.95%5.80%
t p w p e r 0.05%33.05%15.72%
t v w p e r 18.69%45.27%36.94%
Table 9. Details of the optimized solutions.
Table 9. Details of the optimized solutions.
IDVehicle Routing SolutionsTravel Distance (km)
1 r o u t e   1 :   d e 1 , 330 + , 330 , 363 + , 363 , 380 + , 521 + , 380 ,   521 , d e 1
r o u t e   2 :   d e 1 , 389 + , 168 + , 168 , 407 + , 407 , 389 , 440 + , 440 , d e 1

r o u t e   168 :   d e 3 , 338 + , 338 , 272 + , 272 , 299 + , 8 + , 299 , 8 , d e 3
136.10
103.31

139.30
2 r o u t e   1 : d e 1 , 330 + , 373 + , 330 , 380 + , 373 , 380 , 697 + , t s 2 , 650 , 697 , d e 1
r o u t e   2 : d e 1 , 389 + , 392 + , 392 , 407 + , 407 , 389 , t s 2 , 445 , 522 + , 522 , d e 1

r o u t e   159 :   d e 3 , 338 + , 338 , 296 + , 296 , 352 + , 439 + , 439 , 352 , 554 + , t s 2 , d e 3
179.36
102.94

116.18
3 r o u t e   1 : d e 1 , 300 + , 300 , 303 + , 303 , d e 1
r o u t e   2 : d e 1 , 293 + , 293 , 38 + , 38 , d e 1

r o u t e   298 : d e 2 , 61 + , 285 + , 285 , 61 , d e 2
33.93
75.61

64.49
4 r o u t e   1 : d e 1 , 300 + , 300 , t s 3 , 212 , 234 + , 234 ,   d e 1
r o u t e   2 : d e 1 , 293 + , 293 , t s 2 , 295 , 124 + , 124 ,   d e 1

r o u t e   326 : d e 2 , 285 + , 285 , t s 1 , 249 , 149 + , 149 ,   d e 2
83.26
124.94

135.67
5 r o u t e   1 :   d e 1 , 431 + , 432 + , 432 , 461 + , 431 , 461 , d e 1
r o u t e   2 : d e 2 , 388 +   , 388 , 421 + , 421 , 481 + , 481 , d e 2

r o u t e   153 :   d e 3 , 306 + , 306 , 381 + , 401 + , 381 , 401 , 497 + , 497 , d e 3
91.63
115.18

151.15
6 r o u t e   1 :   d e 1 , 431 + , 144 + , 144 , t s 1 , 432 , 461 + , 431 , 461 , d e 1
r o u t e   2 :   d e 2 , 388 + , 388 , 1 + , 481 + , t s 3 , 481 , d e 2

r o u t e   153 :   d e 3 , 306 + , 306 , 381 + , 401 + , 381 , 401 , 185 + , t s 3 , 497 , 185 , d e 3
104.28
149.30

197.87
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheng, H.; Zhang, X.; Chen, J. Study on Customized Shuttle Transit Mode Responding to Spatiotemporal Inhomogeneous Demand in Super-Peak. Information 2021, 12, 429. https://doi.org/10.3390/info12100429

AMA Style

Zheng H, Zhang X, Chen J. Study on Customized Shuttle Transit Mode Responding to Spatiotemporal Inhomogeneous Demand in Super-Peak. Information. 2021; 12(10):429. https://doi.org/10.3390/info12100429

Chicago/Turabian Style

Zheng, Hao, Xingchen Zhang, and Junhua Chen. 2021. "Study on Customized Shuttle Transit Mode Responding to Spatiotemporal Inhomogeneous Demand in Super-Peak" Information 12, no. 10: 429. https://doi.org/10.3390/info12100429

APA Style

Zheng, H., Zhang, X., & Chen, J. (2021). Study on Customized Shuttle Transit Mode Responding to Spatiotemporal Inhomogeneous Demand in Super-Peak. Information, 12(10), 429. https://doi.org/10.3390/info12100429

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