Post on 23-Dec-2015
1
Minimizing Waiting Time at Intermediate Nodes for Intercity Public Bus Transportation
Saharidis G.K.D.
Dimitropoulos Ch.
Skordilis E.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
ΠΑΝΕΠΙΣΤΗΜΙΟ ΘΕΣΣΑΛΙΑΣ
2
Introduction
The purpose of this research is to create a suitable timetable for intercity buses, departing from various nodes of the network, such that the total waiting time of passengers at intermediate nodes is minimized.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
3
Description of the problem
• By using the current infrastructure of the public bus companies in Greece (dubbed KTEL), we tried to reform the bus schedule, for better service of passengers who use intermediate nodes to reach their destination.
• In many cases, the waiting time at these nodes is very long. The reason for this is that there are interconnections of various itineraries, which increments the difficulty of a suitable time schedule.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
4
Formulation of the problem
• The problem is formulated as a mixed integer linear program.
• The objective function is the minimization of the sum of waiting times for every intermediate node of the network
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
5
Formulation of the problem
• Assumptions:
1. Steady travel time between nodes2. Steady number of itineraries (routes) between connected
nodes3. Sufficient parking space in bus terminals (no bottleneck)4. Generally, parameters describing the problem are
considered steady
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
6
Decision variables
• The departure time of every bus ()
• The waiting time at an intermediate node ()
• A binary variable that defines the combination of routes to and from an intermediate node as active or inactive ()
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
7
Decision variables
• The time difference between the bus departure and a high priority time (will be used as penalty factor) ()
• A binary variable that defines the number of buses affected by the high priority time ()
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
8
Illustration of the problem
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
Το σχήμα περιλαμβάνει3 παράλληλες γραμμέςπου αναπαριστούν τους κόμβους αναχώρησης, άφιξηςκαι τον ενδιάμεσο. Το χρησιμοποιούμεγια να δείξουμε πως υπολογίζουμετον χρόνο αναμονής στον ενδιάμεσο.
9
Sets, Parameters
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
10
Sets, Parameters
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
11
Constraints • Constraint (1) defines the combination of
routes to and from an intermediate node as active or inactive, based on whether the bus from the intermediate node has already left or not
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
,, , , ,, , , , , ,
,, , , ,, , , , , ,
, ,
( _ ) (2 )
1 ( _ ) (2 )
, , , ,
m ni k i k k j i jk j n i k m i k jk
m ni k i k k j i ji k j k j n i k m k
i k k j
sm LT LT TT WT MIN BigM D D D Y
Y sm LT LT TT WT MIN BigM D D D
i k j m R n R
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E. 12
Constraints
• Constraint (2) calculates the waiting time between the arrival on each intermediate node and the first two available departing routes from that node towards the destination node.
,, , , ,, , , , , ,
,
_ (2 )
, , ,
m ni k i k k j i jk j n i k m i k jk
i k
LT LT TT WT MIN BigM D D D WT
i k j m R
13
Constraints • Constraint (3) assures that a bus departing
from any node will return to this node before and after a certain time
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
, , ,, , , , ,,(1 ) _ , ,i j j i i ji j m j i m i ji jBigM H H LT TT BR MIN LT i j m R
, , j ,, , , , ,,_ (1 ) , ,i j i j ij i m i j m i ji jLT LT TT BR MAX BigM H H i j m R
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E. 14
Constraints • Constraint (4) assures that each bus departs
after another.
• Constraint (5) defines the time window in which a bus must depart
, , , 1 , , ,,(1 ) _ , , (1,..., )i j i j m i j m i ji jBigM D LT NR MIN LT i j m R
,, , , , 1 ,,_ (1 ) , , (1,..., )i ji j m i j m i ji jLT LT NR MAX BigM D i j m R
, ,, , , , ,, ,i j mi j m i j m i jST LT ET i j m R
15
Extensions
Furthermore, an extension regarding high priority times is introduced. By using this extension bus departures are gravitated towards certain times.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E. 16
Extensions• Constraints (7) and (8) calculates the time
difference between the departing time of a route and the time where a high passenger demand occurs.
• Constraint (9) defines the number of buses affected by these high priority times.
, , ,, , , , , , , , ,(2 ) , , ,i j s i ji j m s i j m i j m s i jG BT LT BigM D B i j s m R
, , ,, , , , , , , , ,( ) (2 ) , , ,i j s i ji j m s i j m i j m s i jG BT LT BigM D B i j s m R
, ,, , , ,0
, , ,M
i j si j m s i jm
B C i j s m R
, ,
, ,
, ,
, ,, ,
, mod 0
1, mod 0
i j si j s
i j s
i j si j s
Lif L MC
MCC
Lif L MC
MC
17
Objective function•Minimization of the waiting time for all
available buses departing from the intermediate node
• Using the penalty factor :
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
, ,,,
, ,, , , , , ,0 0 0 0 0 0 0 0 0
min ( 1)
, ,
k j i ji k R RRI K I I I Sm n
i j sk j i k j i j m si k j m n i j m s
R n WT A G
k i k j i j
,,,
, , ,0 0 0 0 0
min ( 1) , ,k ji k RRI K I
m nk j i k j
i k j m n
R n WT k i k j i j
18
Case study
• The formulation was applied to the intercity bus network of the island of Crete.
• Numerical data for Crete network:122 nodes total.13 of these nodes considered intermediate.Number of routes between nodes varies
from 2 to 25.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
19
Case study
We will examine three different cases:
Case 1 : Existing bus timetable (benchmark)
Case 2 : Minimizing waiting time without penalty factor
Case 3 : Minimizing waiting time with penalty factor to better approach the existing bus timetables
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
20
• For the penalty factor we introduced five high priority times, each corresponding to five time zones
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
330 540 750 960 1170 1380Minute of the day:
435 645 855 1065 1275Zone 1 Zone 2 Zone 3 Zone 4 Zone 5
Case study
High Priority Time:
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E. 21
Case study
• For each time zone, the number of buses affected by the high priority time was set by the number of buses departing, based on the existing bus timetable
• This also applies to the weight factor of the penalty cost in the objective function
22
Results
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
Total Waiting Time
(%)* Number of Active Routes
(%)* Average Waiting Time
(%)*
Case 1 1,994,385 - 33051 - 60 -
Case 2 493,924 75.2 19361 41.4 25 58.3
Case 3 925,444 53.5 25565 22.6 36 40.0
*Percentage of decrease over existing timetable (Case 1).
23
Results• Distribution of waiting times :
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
• Distribution of itineraries :
24
Conclusions• For all the examined cases there is a clear
reduction of the waiting times• The use of a penalty factor (Case 3)results in an
increased waiting time compared to not using it (Case 2)• However, it is greatly reduced compare to the
existing bus timetable• Using of a penalty factor (Case 3) leads to
timetable more suitable for realistic cases, giving passengers more choices considering their departure from intermediate nodes
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E.
Saharidis G.K.D., Dimitropoulos Ch., Skordilis E. 25
ReferencesSteer Davies Gleave. Study of passenger transport by coach. Publication TREN/E1/409-2007. European Commission, 2009.
Ceder, A, Golany B, Tal O. (2001) Creating Bus Timetables with Maximal Synchronization. Transportation Research Part A . 35:913-928.
Eranki, A. “A model to create bus timetables to attain maximum synchronization considering waiting times at transfer stops”. Master’s thesis. Department of Industrial and Management Systems Engineering, University of South Florida, 2004.
Ibarra-Rojas, O. J., and Y. A. Rios-Solis. “Synchronization of bus timetabling”. Transportation Research Part B, Vol. 46, 2012, pp. 599-614.
Hall R., Dessouky M. and Lu Q. “Optimal Holding Times at Transfer Stations”. Computer and Industrial Engineering. Vol. 40, 2001, pp. 379-397.
Bussieck M. R., Winter T., and Zimmermann U. T. “Discrete Optimization in public rail transport”. Mathematical Programming. Vol. 79, Issue 1-3, 1997, pp. 415-444.
Goverde R. M. P. “Synchronization Control of Scheduled Train Services to Minimize Passenger Waiting Time”. Transport, Infrastructure and Logistics, Proceedings 4th TRAIL Congress, 1998.
Chen D. and Wu K. “Research on Optimization Model and Algorithm of Initial Schedule of Intercity Passenger Trains based on Fuzzy Sets”. Journal of Software, Vol. 7, No. 1, 2012, pp. 49-54.
Reinhardt L. B., Clausen T., and Pisinger D. “Synchronized dial-a-ride transportation of disabled passengers at airports”. European Journal of Operational Research. Vol. 225, 2013, pp. 106-117.
Wong R. C. W., Yuen T. W. Y., Fung K. W. and Leung J. M. Y. “Optimizing Timetable Synchronization for Rail Mass Transit”, Transportation Science, Vol. 42, No. 1, 2008, pp. 57-69.