Next Article in Journal
Arc-Completion of 2-Colored Best Match Graphs to Binary-Explainable Best Match Graphs
Previous Article in Journal
Median Filter Aided CNN Based Image Denoising: An Ensemble Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Optimization of Multi-Objective Regional Public Transportation Scheduling

School of Traffic & Transportation Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
*
Author to whom correspondence should be addressed.
Submission received: 25 February 2021 / Revised: 23 March 2021 / Accepted: 26 March 2021 / Published: 28 March 2021

Abstract

:
The optimization of bus scheduling is a key method to improve bus service. So, the purpose of this paper is to address the regional public transportation dispatching problem, while taking into account the association between the departure time of buses and the waiting time of passengers. A bi-objective optimization model for regional public transportation scheduling is established to minimize the total waiting cost of passengers and to maximize the comprehensive service rate of buses. Moreover, a NSGA-II algorithm with adaptive adjusted model for crossover and mutation probability is designed to obtain the Pareto solution set of this problem, and the entropy weight-TOPSIS method is utilized to make a decision. Then the algorithms are compared with examples, and the results show that the model is feasible, and the proposed algorithms are achievable in solving the regional public transportation scheduling problem.

1. Introduction

With the continuous expansion of urban space scale and population, public transportation has become a travel choice for more and more people. On ordinary weekdays, the number of public transport trips in Beijing is about 6.1 million people/day, accounting for 30% of the city’s permanent resident population, and on weekends, it is about 4.4 million people/day. Among them, the proportion of the direct bus travel is only 43.1% [1]. However, the lack of regional public transportation coordinated dispatching tends to result in a consequent worsening of both the magnitude and variability of the average waiting time. This in turn impacts heavily impacts the level of service. Owing to the inconvenience of transfer, public transportation in China is not attractive to travelers. Some studies have shown that transfer is one of the main factors that affect the travel rate of public transportation [2]. The travel time with transfer is obviously longer than the travel time without transfer in Beijing. The average time of a trip in the public transportation system is about 1 h, and the transfer time is more than 10 min [1,3]. So, it is necessary to study the coordinated dispatching problem of regional public transportation to reduce the transfer waiting time of passengers and improve the service level of public transportation.

2. Literature Review

Compared with single-line public transportation dispatching, regional public transportation dispatching can make public resources more fully utilized [4,5]. Transfer is the main factor of regional public transportation scheduling, and synchronous optimization of regional public transportation timetable can improve the transfer efficiency. Therefore, some researchers have studied the problem of synchronous transfer of regional public transportation scheduling. Liebchen and Stiller [6] considered the bus timetable affected by running late to achieve the optimization of departure interval and fleet size. Liu and Shen [7] established a bilevel programming model between the timetable generation and vehicle scheduling in the bus scheduling system according to the regional public transport scheduling mode. Ceder, Golany and Tal [8] established a collaborative optimization model to maximize the number of buses arriving at the station synchronously in the region. Ibarra-Rojas and Rios-Solis [9] optimized the departure schedule under the determination of the synchronous transfer time window to maximize the synchronous number. Eranki [10] considered the problem of developing synchronized timetables for bus transit systems with fixed routes to have maximum number of simultaneous arrivals. Wei, Chen, and Sun [11] optimized the departure schedule based on the shortest waiting time of passengers and the maximum number of bus berths. Liu et al. [12] evaluated the degree of synchronous bus arrival by defining the synergetic coefficient. Nesheli and Ceder [13] took the total travel time as the target and adopts the synchronous transfer problem of buses under different strategies. Hu, Zhao and Wang [14] realized synchronous transfer between urban bus lines and rural bus lines by optimizing the vehicle arrival time of transfer routes. Wagale et al. [15] presented a model to optimize the bus scheduling by taking into consideration of both bus stop and route segments of the city. Taking into account congestion on the public transport system and works by considering elastic demand, Ibeas et al. [16] proposed an optimization model for designing the intervals and sizes of buses circulating on public transport networks by minimizing the system’s operating and user costs. Wang et al. [17] built an optimal model to schedule the departure time of each bus service with the objective of minimizing the average waiting time. Boyer, Ibarra-Rojas and Ríos-Solís [18] dealt with the flexible vehicle and crew scheduling problem faced by urban bus transport agencies. Kumar, Prasath and Vanajakshi [19] proposed a demand and travel time responsive model to maximize the benefit of the operator by preparing an optimal schedule. Kang, Chen and Meng [20] developed three explicit integer linear programming models to formulate: the bus driver scheduling, bus and driver scheduling, and bus and driver scheduling with mealtime windows, respectively. Burdett et al. [21], Bevrani et al. [22,23], Ghaderi et al. [24] formulated some multi-objective capacity models to deal with transportation planning problem.
A review of the literature on regional public transport scheduling reveals that the existing research mainly considers the constraints on the waiting time for passengers, system’s operating cost or the degree of transfer synchronous. However, the losses of non-synchronized transfer passengers who choose to use alternative modes of transport because of the waiting and the bus capacity are rarely considered. The total waiting cost of passengers should include the waiting cost of passengers who have not achieved easy transfer and those who have to choose other modes of transportation, and the waiting time cost of non-transfer passengers. In terms of operating cost, departure frequency or full load rate are always considered. Effective operating time is also a major component of operating costs. In view of the above, this paper proposes a multi-objective optimization model of regional public transportation to minimize the total waiting cost of passengers and to maximize the comprehensive service rate of buses. According to the characteristics of this model, an improved NSGA-II algorithm is designed to obtain the Pareto solution set of this problem, and the entropy weight-TOPSIS method is used to make a decision.
The remainder of the paper is organized as follows: Some mathematical notations are defined in Section 3.1, and the bus operation process and passenger waiting time are analyzed in Section 3.2. A bi-objective optimization model for regional public transport scheduling is established in Section 3.3. The NSGA-II algorithm with adaptive adjusted model for crossover and mutation probability as the first step of the two-step algorithm is proposed in Section 4.1, and the entropy weight-TOPSIS method as the second step is implemented in Section 4.2. The proposed algorithm is tested via a numerical example in Section 5. Finally, the conclusions are given in Section 6.

3. The Regional Public Transportation Scheduling Model

Compared with traditional single-line public transportation dispatching, regional public transportation dispatching is more complicated. If only considering the number of easy transfer buses in the process of dispatching, it may lead to ignoring the transfer satisfaction of passengers at other transfer stations. So, the number of passengers who have not realized easy transfer at all the transfer stations should be taken into account so as to improve the service level of public transport.
Some assumptions involved in this paper are as follows:
(1)
The buses operating on each line are considered to be full-range buses.
(2)
There are no overtaking, stopovers and U-turns.

3.1. Notations

(1) Parameters
L = { 1 , 2 , , N } is the bus line, N is the number of lines;
s is the station set, s i , q represents station i of line q , and m i is the number of stations on line i ;
S t = { ( s i , q , s j , p ) | s i , q , s j , p S , i , j L } is the set of pairs for transfer stations;
x i X implies the bus x uses line i, where X is the bus set;
V is the passenger capacity of the bus;
p x i s i , q a , p x i s i , q g denote the arrival number and alighting rate of passengers on the bus x i at station s i , q , respectively;
f s i , q s j , p is the proportion of passenger transfer from line i to j at transfer-station s i , q , s j , p ;
Q x i s i , q a and Q x i s i , q g indicate the number of passengers getting on and off the bus x i at station s i , q of line i ;
F x i s i , q is the number of passengers who failed to get on the bus x i of line i at station s i , q ;
C x i s i , q indicates the number of passengers on the bus x i of line i at station s i , q ;
t s i , q 1 s i , q denotes the travel time from station s i , q 1 to station s i , q , when q = 1 , t s i , q 1 s i , q represents the time from the parking lot to the departure station;
U x i s i , q indicates the dwelling time of bus x i of line i at station s i , q ;
δ 1 and δ 2 indicate the time window for easy transfer;
t x i s i , q indicates the time when the bus x i of line i arrives at station s i , q ;
l x i s i , q denotes the time when the bus x i of line i leaves station s i , q ;
T 1 and T 2 are the scheduling simulation time window;
I min and I max indicate the minimum and maximum departure interval during the simulation;
T x j x i s j , p s i , q { 1 , δ 1 t x i s i , q l x j s j , p δ 2 0 , other indicates whether the bus x j of line j at station s j , q can be easily transferred to the bus x i of station s i , q on line i, where δ 1 , δ 2 are the time window for synchronous transfer;
P indicates the proportion of choice other transportation modes for passengers who are failed to get on the bus;
Y x i s i , q indicates the number of passengers who failed to get on the bus at the bus x i of line i at station s i , q and chose other transportation modes;
R x j x i s j , p s i , q = { 1 , l x j s j , p t x i s i , q 0 , other indicates whether the two lines can be transferred.
(2) Decision variables
A x i indicates the departure time of bus x i of line i from the parking lot.

3.2. Operation Analysis

(1) Analysis of bus operation process
The dwelling time of each bus can be expressed as U x i s i , q = a 0 + max ( a 1 Q x i s i , q g , a 2 Q x i s i , q a ) , where a 0 is the time of buses to enter and leave the station, a 1 and a 2 indicate the average boarding and alighting time for each passenger, respectively.
The bus arrival time at the departure station can written as t x i s i , 1 = A x i + t s i , 0 s i , 1 , where t s i , 0 s i , 1 indicates the running time of line i from the parking lot to the departure station, and the bus departure time can be computed as l x i s i , q = t x i s i , q + U x i s i , q .
The total running time of a bus can be calculated by l x i s i , q A x i .
(2) Passenger number analysis
The number of passengers getting off the bus Q x i s i , q g is determined by the number of passengers on the bus, C x i s i , q , and the alighting rate of passengers, p x i s i , q g , which can be calculated using
Q x i s i , q g = C x i s i , q p x i s i , q g
There are two aspects to consider in terms of the number of passengers getting on the bus: the remaining capacity of the bus and the number of passengers the who expect to board the bus which includes the number of people arriving randomly, the number of people who need to transfer, and the number of people who fail to get on the previous bus (do not include those who choose other modes of transportation). So, the number of passengers getting on the bus can be formulated as follows:
Q x i s i , q a = min ( V ( C x i s i , q Q x i s i , q g ) , ( t x i s i , q t x i 1 s i , q ) p x i s i , q a + R x j x i s j , p s i , q Q x j s j , p g f s j , p s i , q + F x i 1 s i , q Y x i 1 s i , q )
Then, the number of passengers on the bus at station, C x i s i , q , is the number of passengers on the vehicle, C x i 1 s i , q , minus the number of passengers getting off the vehicle, Q x i 1 s i , q g , plus the number of passengers getting on the vehicle, which can be obtained:
C x i s i , q = C x i 1 s i , q + Q x i 1 s i , q a Q x i 1 s i , q g
In addition, the number of passengers who have not achieved easy transfer can be expressed as Q x j s j , p g f s j , p s i , q ( 1 T x j x i s j , p s i , q ) .
Thus, the number of passengers who fail to get on the bus is the number of passengers who need to get on the bus in the current period minus the number of passengers who actually get on the vehicle, which can be computed as follows:
(3) Analysis of passenger waiting time
F x i s i , q = ( t x i s i , q l x i 1 s i , q ) p x i s i , q a + R x j x i s j , p s i , q Q x j s j , p g f s j , p s i , q   + F x i 1 s i , q Y x i 1 s i , q Q x i s i , q a .
The waiting time for transfer passengers is the arrival departure interval of two buses with transfer relationship multiplied by the number of transfer people with easy transfer, which can be calculated as follows:
H 1 = x i , x j X i , j L ( s j , p , s i , q ) S t ( t x i s i , q l x j s j , p ) Q x j s j , p g f s j , p s i , q R x j x i s j , p s i , q .
Furthermore, the waiting time for non-transfer passengers includes the average waiting time of random arriving passengers and the extra waiting time of passengers who fail to get on the bus according to Larsen and Sunde [25], which can be calculated by
H 2 = x i X q = 1 m i i L ( ( t x i s i , q t x i 1 s i , q ) 2 p x i s i , q a 2 + ( F x i s i , q Y x i s i , q ) ( t x i + 1 s i , q t x i s i , q ) ) .
The ratio R of the inter-station operating time of all buses to the total operating time is 1 minus the proportion of stop time, which can reflect the effective operation cost and can be expressed as:
R = 1 x i X q = 1 m i i L U x i s i , q x i X i L ( l x i s i , m i A x i )

3.3. The Optimization Model

Based on the above analysis, the first objective of this model is to minimize the total waiting cost of passengers who have not achieved easy transfer and those who have to choose other modes of transportation, and the waiting time cost of non-transfer passengers and transfer passengers. The second objective of this model is to consider the efficiency of public transport operation, including the average full load rate of buses and the ratio of inter-station running time of all buses to the total operation time. Then, the bi-objective optimization model for regional public transport scheduling can be formulated as follows:
min Z 1 = c 1 ( x i , x j X i , j L ( s j , p , s i , q ) S t Q x j s j , p g f s j , p s i , q R x j x i s j , p s i , q ( 1 T x j x i s j , p s i , q ) + x i X q = 1 m i i L Y x i s i , q ) + c 2 ( H 1 + H 2 )
max Z 2 = x i X i L max ( C x i s i , q ) V | X | + R
s.t.
T 1 A x i T 2 , x i X
C x i s i , q V 1 , x i X , s i , q S
I min A x i A x i 1 I max , x i X
Formula (8) indicates the total waiting cost of passengers, c 1 , c 2 represent the waiting cost coefficients of passengers who did not realize easy transfer and choose other modes of transportation, and the time value cost of waiting time. Formula (9) indicates the comprehensive service ratio of buses, | X | is the number of elements of the set X. Constraint (10) ensures that the departure time is within the scheduling simulation time window. Constraint (11) is the bus capacity constraint. Constraint (12) indicates the arrival interval limits.

4. Solution Algorithms

As most of the multi-objective problems have to meet several objectives, it is impossible to obtain a unique solution that satisfies all the objectives simultaneously. Usually, Simulated Annealing (SA) [21], Ant Colony Optimization (ACO) [26], Multi-objective Heuristic Algorithms [27], Discrete Evolutionary Multi-objective Optimization (DEMO) [28], Niched Genetic Algorithm [5], Non-dominated Sorting Genetic Algorithm (NSGA) [29,30] are often chosen to solve multi-objective problems. One of the widely used algorithms proposed to solve multi-objective problems is the Non-dominated Sorting Genetic Algorithm-II (NSGA-II).
The scheduling of different bus lines will affect the transfer efficiency. Moreover, there are many schemes for each bus line and the two objectives conflict with each other. Therefore, it is necessary to design an algorithm that can deal with multi-objective public transportation scheduling quickly. In this section, a two-step algorithm is employed in order to solve this problem. The first step is to obtain the Pareto solutions, and the second step is to obtain a satisfactory scheme from the Pareto solution set on the basis of the entropy weight-TOPSIS method.

4.1. The NSGA-II Algorithm

The NSGA-II algorithm is one of the most commonly used multi-objective optimization algorithms, but it has some drawbacks of an uneven convergence of the population and the ease of falling into the local optimal solution [31]. In order to avoid premature convergence of Genetic algorithm, it is necessary to keep the diversity of population. According to the characteristics of this model, an adaptive adjusted model of the crossover and the mutation probability is introduced, and an adaptive crossover operator for judging the similarity of chromosomes and multipoint mutation operator are designed.
(1) Construction of chromosomes
Chromosome X = ( X 1 , X 2 , , X N ) is designed by a positive integer, where gene X i = ( A 1 , A 2 , , A x i ) of the chromosome represents the set of the departure time for each bus of line i. Figure 1 illustrates the specific structure of a randomized chromosome.
The process of generating a chromosome is as follows:
Step 1: Let x i = 1 and A 1 = r 0 ( I max I min ) + I min + T 1 (where r 0 is a random number obeying 0–1 uniform distribution).
Step 2: x i = x i + 1 , If A x i 1 + r 0 ( I max I min ) + I min T 2 , let A x i = A x i 1 + r 0 ( I max I min ) + I min , return to Step 2. Otherwise, go to Step 3.
Step 3: if i < N , i = i + 1 , return to Step 2. Otherwise, stop.
(2) Select operator
The selection of the parental population is using binary tournament method calculate the performing non-dominated sorting and crowding calculation of the parental population.
(3) Adaptive adjusted model for crossover and mutation probability
In order to improve the search performance of the algorithm and avoid falling into the local optimal solution, an adjusted model is adopted to adjust the crossover probability and mutation probability.
The crossover probability can be adjusted by the following method.
P c = { g 1 g 4 g 1 + P c 1 , 0 g g 1 , g 2 g 4 ( g 2 g 1 ) + P c 2 , g 1 g g 2 , g max g 2 ( g max g 2 ) ( 1 β ) + P c 2 β , g 2 g g max ,
where P c 1 , P c 2 are the crossover probability for different generations, β is the adjustment coefficient of crossover and mutation probability in the later stage of evolution. g is the current evolutionary generation, g 1 , g 2 are the dividing points of the early stage and the middle stage, and g max is the total maximum evolution generation.
Furthermore, the mutation probability can be adjusted by the following method according to the running phase.
P m = { ( P m 1 1 L c ) ( g 1 g ) g 1 + 1 L c , 0 g g 1 , ( 1 L c P m 1 L c ) ( g 2 g ) g 2 g 1 + P m 1 L c , g 1 g g 2 , P m 1 L c ( g max g ) ( g max g 2 ) ( 1 β ) + P m 1 L c β , g 2 g g max .
where L c represents the total number of bus departures for all lines in a chromosome, P m 1 is the maximum mutation probability.
(4) Crossover operator
According to the character of the chromosome structure, the traditional single-point crossover is easy to fall into the local optimal solution. Therefore, an adaptive crossover operator is designed through judging the similarity of chromosomes.
The process of adaptive crossover operator is stated as follows:
Step1: Let i = 1 , k = 0 .
Step2: If i > N , then go to step 4. Otherwise, compare the number of bus departures of a line for the two selected chromosomes, if the departure times are equal, let the departure times of line i be x i and go to Step3, otherwise, let i = i + 1 , return to Step2.
Step3: Count the number Z of equal departure time of lines i of the two chromosomes. If Z x i 2 , then regenerate the departure time of line i for one chromosome and exchange the genes of line i of two chromosomes, and let i = i + 1 , k = k + 1 , go to Step2.
Step4: If k 1 , randomly generate two different lines and exchange the corresponding genes of two chromosomes.
In the case of Z x i 2 , the crossover process is shown in Figure 2. Assume that the first genes of parents P 1 , P 2 are highly similar to the N 1 genes, then randomly select one chromosome from P 1 , P 2 , regenerate the departure time of the high similarity line, and exchange the corresponding genes of P 1 , P 2 (where C x i P q is the regenerated departure time). The crossover process when k 1 is shown in Figure 3.
(5) Mutation operator
Randomly select two different lines from a chromosome and calculate the length of each line gene L g .
Randomly select some places form line gene and regenerate these departure time according to the following conditions ( A x i , A x i 1 and A x i + 1 denote the departure time of bus x i , x i 1 and x i + 1 respectively, and r 0 indicates a random number from uniform distribution u ( 0 , 1 ) ).
Case 1: When x i = 1 , if A x i + 1 2 I min > I max , then A x i = ( I max I min ) r 0 + I min . Otherwise, if A x i + 1 2 I min I max I min , then A x i = ( A x i + 1 2 I min ) r 0 + I min , otherwise, A x i = ( A x i + 1 3 I min ) r 0 + I min .
Case 2: When x i = L g , A x i = A x i 1 + ( I max I min ) r 0 + I min .
Case 3: When 1 < x i < L g , if A x i + 1 A x i 1 2 I min > I max , then A x i = A x i 1 + ( I max I min ) r 0 + I min . Otherwise, if A x i + 1 A x i 1 2 I min I max I min , then A x i = A x i 1 + ( A x i + 1 A x i 1 2 I min ) r 0 + I min , otherwise A x i = A x i 1 + ( A x i + 1 A x i 1 3 I min ) r 0 + I min .
Figure 4 and Figure 5 are used to illustrate the mutation process for Case 1 and Case 3 respectively, where I min = 3 min , I max = 10 min and the scheduling time window is 30 min.
Since Case 2 has only one calculation method, we will give examples of Cases 1 and 3 respectively, the shaded part is the point to be mutated, assuming that the departure interval is [3,10] min and the scheduling time window is 30 min.

4.2. Obtain Satisfactory Scheme from Pareto Solution Set

The second step is to obtain a satisfactory scheme from the Pareto solution set. The entropy weight-TOPSIS method is an improvement over the entropy weight method and the TOPSIS method [32]. Moreover, it is a common evaluation method used in multi-objective decision-making problems. Consequently, the entropy weight-TOPSIS method is used to obtain satisfactory schedule from the Pareto set.
Suppose that the number of Pareto solutions is M. Then, a decision matrix D M × 2 is obtained with M Pareto solutions and 2 objectives. The process is stated as follows:
Step 1: Obtaining the standardization matrix B i j = [ b i j ] M × 2 and normalized matrix C i j = [ c i j ] M × 2 .
Step 2: Calculating the information entropy of each objective j:
K j = ( 1 ln M ) i = 1 M c i j ln c i j , j = 1 , 2   ( if   c i j = 0 ,   then   c i j ln c i j = 0 ) .
Step 3: Computing the weight of objective j:
w j = 1 K j q = 1 j ( 1 K q ) , j = 1 , 2 .
Step 4: Computing the weighted judgment matrix:
Q M × 2 = [ q i j ] M × 2 ,   where   q i j = w j b i j , i = 1 , 2 , M , j = 1 , 2 .
Step 5: Calculating the Euclidean distance and closeness for each solution:
D i + = j = 1 2 ( q i j q j + ) 2 , i = 1 , 2 , M , j = 1 , 2 ,
D i = j = 1 2 ( q i j q j ) 2 , i = 1 , 2 , M , j = 1 , 2 ,
F i = D i D i + D i + , i = 1 , 2 , M .
Moreover, the satisfactory schedule can be obtained according to the closeness F i .

5. Example Analysis

A case study is carried out to test the dispatch approach and algorithm on a regional public transportation network with four bus lines, as shown in Figure 6. The scheduling simulation time window T 1 and T 2 were set to be 7:00, 7:30, respectively. The minimum and maximum departure interval, I min and I max , were set to be 3 min, 10 min, and the time window for easy transfer, δ 1 and δ 2 , were 1 min and 5 min, respectively. Let the waiting cost coefficient c 1 = 1.5 , c 2 = 0.18 . The passenger capacity of the bus V = 100 .
Based on the investigation of bus transfer stations in line 7, line 15 and line 131 of Lanzhou City, it is found that the more the number of passengers who failed to get on the bus, the higher the proportion of passengers who choose other modes of transportation. The corresponding proportion of choice other transportation modes for those passengers is as follows:
P = { 0.09 , 1 F x i s i , q 5 , 0.18 , 6 F x i s i , q 10 , 0.25 , 11 F x i s i , q .
The size of population is 50, the maximum number of iterations is 150, the maximum mutation rate P m 1 = 0.3 , and P c 1 = 0.6 , P c 2 = 0.4 . The main transfer stations and transfer proportions are shown in Table 1, and the average arrival rate, alighting rate of passengers and travel time between stations are listed in Table 2.
The Pareto solutions obtained and the weighted objective value of them are shown in Table 3, and the Euclidean distance between solutions and positive and negative ideal points and the closeness of solutions are listed in Table 4 after using the entropy-weight TOPSIS method. It can be seen from Table 4 that solution 7 is the most satisfactory operation schedule, and its operation schedule is given in Table 5.
By comparing solution 7 with the average value of schedules at a fixed departure interval of 3~10 min (as shown in Table 6), it can be seen that the total cost of loss for passengers is reduced by 31.42%, the comprehensive service ratio of buses is increased by 0.04%, and the number of passengers who not achieved easy transfer decreased by 74.29%.
For comparing the performance of NSGA-II with the fixed crossover probability and mutation probability (FNSGA-II) and NSGA-II with adaptive adjusted model for crossover and mutation probability (ANSGA-II), the two algorithms solve the test experiment for 10 times. The Pareto solutions obtained of two algorithms after running 10 times are shown in Figure 7, and the number of nondominated set with different iterations of each algorithm is shown in Figure 8. The results show that the obtained solution set by ANSGA is superior to that obtained by FNSGA-II.

6. Conclusions

A bi-objective optimization model for regional public transport scheduling is established considering the total cost of passengers and the efficiency of public transport operation. Furthermore, a two-step algorithm is designed to solve this problem. In the first step, the NSGA-II algorithm with adaptive adjusted model for crossover and mutation probability is used to obtain the Pareto solutions. In the second step, the entropy weight-TOPSIS method is implemented to find a satisfactory scheme. At last, a regional public transportation network is collected to test the dispatch approach and algorithm. The results show that the scheduling scheme obtained in this paper is obviously better than the schemes with fixed departure interval, and the obtained solution set by ANSGA is superior to that obtained by FNSGA-II.
In this paper, we assumed that the running time between stations is known. In reality, the running time may be uncertain due to delay. Therefore, future research can consider the followings:
(1)
Considering the uncertainty of the running time to design the regional public transport scheduling model and algorithm.
(2)
Considering other factors, such as multiple types of buses, energy consumption, depots capacities, to the regional bus scheduling problem is studied, to research on the regional public transportation scheduling.

Author Contributions

Conceptualization, X.Y.; methodology, X.Y. and Y.Q.; software, X.Y.; validation, Y.Q.; formal analysis, X.Y.; investigation, Y.Q.; data curation, Y.Q.; writing—original draft preparation, Y.Q.; writing—review and editing, X.Y.; supervision, X.Y.; project administration, X.Y.; funding acquisition, X.Y. Both authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 71761024.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beijing Transport Institute: Analysis on the Characteristics of Beijing Residents’ Public Transport Travel. Available online: https://mp.weixin.qq.com/s/oeK1d8iaqrRN8PS39smfyw (accessed on 5 June 2019).
  2. Fu, J.G. The Key Technologies research to improve transit trip rate. Urban Dev. Stud. 2014, 21, 79–83. [Google Scholar]
  3. Guo, Z.; Wilson, N.H.M. Assessing the cost of transfer inconvenience in public transport systems: A case study of the London underground. Transp. Res. Part A Policy Pract. 2011, 45, 91–104. [Google Scholar] [CrossRef]
  4. Wei, M.; Chen, X.W.; Sun, B. Multi-objective single line transit mixed scheduling model considering express bus Service. J. Transp. Syst. Eng. Inf. Technol. 2015, 15, 169–174. [Google Scholar]
  5. Yang, X.F.; Liu, L.F. A Multi-objective bus rapid transit energy saving dispatching optimization considering multiple types of vehicles. IEEE Access 2020, 8, 79459–79471. [Google Scholar] [CrossRef]
  6. Liebchen, C.; Stiller, S. Delay resistant timetabling. Public Transp. 2009, 1, 55–72. [Google Scholar] [CrossRef]
  7. Liu, Z.G.; Shen, J.S. Regional bus operation bi-level programming model integrating timetabling and vehicle scheduling. Syst. Eng. Theory Pract. 2007, 11, 135–141. [Google Scholar] [CrossRef]
  8. Ceder, A.; Golany, B.; Tal, O. Creating bus timetables with maximal synchronization. Transp. Res. Part A Policy Pract. 2001, 35, 913–928. [Google Scholar] [CrossRef]
  9. Ibarra-Rojas, O.J.; Rios-Solis, Y.A. Synchronization of bus timetabling. Transp. Res. Part B Methods 2012, 46, 599–614. [Google Scholar] [CrossRef]
  10. Eranki, A. A Model to Create Bus Timetables to Attain Maximum Synchronization Considering Waiting Times at Transfer Stops. Master’s Thesis, University of South Florida, Tampa, FL, USA, 2004. [Google Scholar]
  11. Wei, M.; Chen, X.W.; Sun, B. A model and an algorithm of schedule coordination for multi-mode regional bus transit. J. Highw. Transp. Res. Dev. 2015, 32, 136–142. [Google Scholar]
  12. Liu, Z.G.; Shen, J.S.; Wang, H.X.; Yang, W. Regional public transportation timetabling model with synchronization. J. Transp. Syst. Eng. Inf. Technol. 2007, 7, 109–113. [Google Scholar]
  13. Nesheli, M.M.; Ceder, A. Optimal combinations of selected tactics for public-transport transfer synchronization. Transp. Res. Part C Emerg. Technol. 2014, 48, 491–504. [Google Scholar] [CrossRef]
  14. Hu, B.Y.; Zhao, H.; Wang, D.X. Synchronous transfer model of urban public transport and rural passenger transport. J. Highw. Transp. Res. Dev. 2019, 13, 144–150. [Google Scholar]
  15. Wagale, M.; Singh, A.P.; Sarkar, A.K.; Arkatkar, S. Real-time optimal bus scheduling for a city using a DTR model. Procedia Soc. Behav. Sci. 2013, 104, 845–854. [Google Scholar] [CrossRef] [Green Version]
  16. Ibeas, A.; Alonso, B.; Dell’Olio, L.; Moura, J.L. Bus size and headways optimization model considering elastic demand. J. Transp. Eng. 2014, 140, 370–380. [Google Scholar] [CrossRef]
  17. Wang, Y.; Zhang, D.X.; Hu, L.; Yang, Y.; Lee, L.H. A data-driven and optimal bus scheduling model with time-dependent traffic and demand. IEEE Trans. Intell. Transp. Syst. 2017, 18, 2443–2452. [Google Scholar] [CrossRef]
  18. Boyer, V.; Ibarra-Rojas, O.J.; Ríos-Solís, Y.Á. Vehicle and crew scheduling for flexible bus transportation systems. Transp. Res. Part B Meth. 2018, 112, 216–229. [Google Scholar] [CrossRef]
  19. Kumar, B.A.; Prasath, G.H.; Vanajakshi, L. Dynamic bus scheduling based on real-time demand and travel time. Int. J. Civ. Eng. 2019, 17, 1481–1489. [Google Scholar] [CrossRef]
  20. Kang, L.J.; Chen, S.K.; Meng, Q. Bus and driver scheduling with mealtime windows for a single public bus route. Transp. Res. 2019, 101, 145–160. [Google Scholar] [CrossRef]
  21. Burdett, R.L. Multi-objective models and techniques for analysing the absolute capacity of railway networks. Eur. J. Oper. Res. 2015, 245, 489–505. [Google Scholar] [CrossRef] [Green Version]
  22. Bevrani, B.; Burdett, R.; Bhaskar, A.; Yarlagadda, P.K.D.V. A multi-criteria multi-commodity flow model for analysing transportation networks. Oper. Res. Perspect. 2020, 7, 100159. [Google Scholar] [CrossRef]
  23. Bevrani, B.; Burdett, R.L.; Yarlagadda, P.K.D.V. A capacity assessment approach for multi-modal transportation systems. Eur. J. Oper. Res. 2017, 263, 864–878. [Google Scholar] [CrossRef] [Green Version]
  24. Ghaderi, A.; Burdett, R.L. An integrated location and routing approach for transporting hazardous materials in a bi-modal transportation network. Transp. Res. Part E Logist. Transp. Rev. 2019, 127, 49–65. [Google Scholar] [CrossRef]
  25. Larsen, O.I.; Sunde, Ø. Waiting time and the role and value of information in scheduled transport. Res. Transp. Econ. 2008, 23, 41–52. [Google Scholar] [CrossRef]
  26. Xiang, W. A multi-objective ACO for operating room scheduling optimization. Nat. Comput. 2017, 16, 607–617. [Google Scholar] [CrossRef]
  27. Dulebenets, M.A.; Pasha, J.; Kavoosi, M.; Abioye, O.F.; Ozguven, E.E.; Moses, R.; Boot, W.R.; Sando, T. Multiobjective optimization model for emergency evacuation planning in geographical locations with vulnerable population groups. J. Manag. Eng. 2020, 36, 04019043. [Google Scholar] [CrossRef]
  28. Han, Y.Y.; Li, J.Q.; Sang, H.Y.; Liu, Y.P.; Gao, K.Z.; Pan, Q.K. Discrete evolutionary multi-objective optimization for energy-efficient blocking flow shop scheduling with setup time. Appl. Soft. Comput. 2020, 93, 106343. [Google Scholar] [CrossRef]
  29. Li, Y.Y.; Wang, S.Q.; Duan, X.B.; Liu, S.J.; Liu, J.P.; Hu, S. Multi-objective energy management for Atkinson cycle engine and series hybrid electric vehicle based on evolutionary NSGA-II algorithm using digital twins. Energy Convers. Manag. 2021, 230, 113788. [Google Scholar] [CrossRef]
  30. Liu, J.; Ma, B.; Zhao, H.B. Combustion parameters optimization of a diesel/natural gas dual fuel engine using genetic algorithm. Fuel 2020, 260, 116365. [Google Scholar] [CrossRef]
  31. Deb, K.; Pratap, A.; Agarwal, S. A fast and elitist multi-objective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  32. Luo, J.; Chen, S.Y.; Sun, X.; Zhu, Y.Y.; Zeng, J.X.; Chen, G.P. Analysis of city centrality based on entropy weight TOPSIS and population mobility: A case study of cities in the Yangtze River Economic Belt. J. Geogr. Sci. 2020, 30, 515–534. [Google Scholar] [CrossRef]
Figure 1. An example of a chromosome.
Figure 1. An example of a chromosome.
Algorithms 14 00108 g001
Figure 2. The crossover process case Z x i 2 .
Figure 2. The crossover process case Z x i 2 .
Algorithms 14 00108 g002
Figure 3. The crossover process case k 1 .
Figure 3. The crossover process case k 1 .
Algorithms 14 00108 g003
Figure 4. Illustration of the mutation operator for case 1.
Figure 4. Illustration of the mutation operator for case 1.
Algorithms 14 00108 g004
Figure 5. Illustration of the mutation operator for case 3.
Figure 5. Illustration of the mutation operator for case 3.
Algorithms 14 00108 g005
Figure 6. A regional public transportation network.
Figure 6. A regional public transportation network.
Algorithms 14 00108 g006
Figure 7. Comparison of calculation results for two algorithms.
Figure 7. Comparison of calculation results for two algorithms.
Algorithms 14 00108 g007
Figure 8. The variation trend of non-dominated individuals.
Figure 8. The variation trend of non-dominated individuals.
Algorithms 14 00108 g008
Table 1. Transfer station and transfer proportion.
Table 1. Transfer station and transfer proportion.
Transfer Station PairTransfer ProportionTransfer StationTransfer Proportion
s l 1 , 2 , s l 3 , 5 0.35 s l 1 , 3 , s l 2 , 4 0.5
s l 1 , 3 , s l 3 , 4 0.45 s l 1 , 4 , s l 2 , 5 0.25
s l 1 , 4 , s l 4 , 4 0.35 s l 1 , 6 , s l 4 , 6 0.15
s l 1 , 8 , s l 4 , 8 0.05 s l 2 , 2 , s l 3 , 2 0.2
s l 2 , 4 , s l 3 , 4 0.4 s l 2 , 5 , s l 4 , 4 0.3
s l 2 , 6 , s l 4 , 3 0.45 s l 2 , 8 , s l 3 , 8 0.05
s l 3 , 7 , s l 4 , 2 0.35
Table 2. The average arrival rate, alighting rate of passengers and travel time between stations.
Table 2. The average arrival rate, alighting rate of passengers and travel time between stations.
l i s l i , q t s i , q 1 s i , q / min p x i s i , q a p x i s i , q g l i s l i , q t s i , q 1 s i , q / min p x i s i , q a p x i s i , q g l i s l i , q t s i , q 1 s i , q / min p x i s i , q a p x i s i , q g
l 1 s l 1 , 1 210 l 2 s l 2 , 3 320.05 l 3 s l 3 , 7 3.510.35
l 1 s l 1 , 2 3.510.05 l 2 s l 2 , 4 320.25 l 3 s l 3 , 8 220.4
l 1 s l 1 , 3 2.520.05 l 2 s l 2 , 5 2.520.2 l 3 s l 3 , 9 2.501
l 1 s l 1 , 4 3.530.1 l 2 s l 2 , 6 320.25 l 4 s l 4 , 1 1.510
l 1 s l 1 , 5 320.25 l 2 s l 2 , 7 310.25 l 4 s l 4 , 2 210.05
l 1 s l 1 , 6 3.530.2 l 2 s l 2 , 8 4.510.4 l 4 s l 4 , 3 220.35
l 1 s l 1 , 7 420.2 l 2 s l 2 , 9 301 l 4 s l 4 , 4 420.25
l 1 s l 1 , 8 310.25 l 3 s l 3 , 1 210 l 4 s l 4 , 5 330.3
l 1 s l 1 , 9 4.510.35 l 3 s l 3 , 2 220.2 l 4 s l 4 , 6 320.25
l 1 s l 1 , 10 310.25 l 3 s l 3 , 3 3.520.15 l 4 s l 4 , 7 2.510.2
l 1 s l 1 , 11 301 l 3 s l 3 , 4 2.510.25 l 4 s l 4 , 8 310.3
l 2 s l 2 , 1 210 l 3 s l 3 , 5 2.520.2 l 4 s l 4 , 9 210.35
l 2 s l 2 , 2 2.510.05 l 3 s l 3 , 6 330.25 l 4 s l 4 , 10 201
Table 3. Pareto solution set and weighted target value obtained by the entropy weight method.
Table 3. Pareto solution set and weighted target value obtained by the entropy weight method.
Solution No. Z 1 / Yuan Z 2 Weighted Objective ValueSolution No. Z 1 / Yuan Z 2 Weighted Objective Value
13910.78 1.5967 1262749.43 1.3955 0.1732
22578.71 1.3270 0272970.66 1.4790 0.3906
32589.56 1.3571 0.0451283491.44 1.5538 0.741
43497.20 1.5612 0.7535292966.80 1.4757 0.3843
52580.83 1.3287 0.0032303691.16 1.5744 0.8646
63795.25 1.5852 0.9292313437.34 1.5525 0.7131
72641.18 1.3760 0.095323223.69 1.5172 0.5633
82609.90 1.3619 0.0613333756.70 1.5768 0.8994
92955.62 1.4616 0.3602343380.22 1.5417 0.6713
103613.80 1.5638 0.8132352898.54 1.4475 0.314
113091.34 1.5025 0.48363869.82 1.5920 0.9741
122718.75 1.3875 0.1477372686.30 1.3821 0.125
133126.04 1.5083 0.5044382813.10 1.4142 0.2287
143084.49 1.4888 0.4586392764.92 1.3994 0.1859
152934.42 1.4532 0.3388402825.46 1.4202 0.2426
163271.71 1.5290 0.602413436.29 1.5454 0.7032
172868.64 1.4271 0.2726423760.29 1.5786 0.9035
183330.26 1.5399 0.6448433175.89 1.5114 0.5325
193688.07 1.5708 0.8582443044.15 1.4863 0.4357
203282.35 1.5371 0.618452777.11 1.4064 0.2011
213041.99 1.4813 0.4281462853.55 1.4225 0.2591
223616.53 1.5653 0.8166473177.62 1.5141 0.5369
232884.80 1.4407 0.2984482677.38 1.3786 0.1161
243867.57 1.5880 0.9676492793.91 1.4095 0.2132
253225.69 1.5231 0.572503396.20 1.5444 0.6826
Table 4. The Euclidean distance and closeness of solutions.
Table 4. The Euclidean distance and closeness of solutions.
Solution No. D i + D i F i Solution No. D i + D i F i
10.6423 0.3577 0.3577 260.2793 0.5673 0.6701
20.3577 0.6423 0.6423 270.2452 0.4961 0.6693
30.3179 0.6383 0.6675 280.4437 0.3625 0.4496
40.4454 0.3691 0.4532 290.2465 0.4960 0.6680
50.3555 0.6413 0.6433 300.5372 0.3449 0.3910
60.5868 0.3471 0.3717 310.4181 0.3762 0.4736
70.2944 0.6156 0.6765 320.3284 0.4164 0.5591
80.3118 0.6289 0.6686 330.5686 0.3396 0.3739
90.2552 0.4939 0.6593 340.3933 0.3829 0.4933
100.5010 0.3452 0.4079 350.2509 0.5136 0.6718
110.2769 0.4586 0.6235 360.6226 0.3521 0.3613
120.2856 0.5803 0.6702 370.2894 0.5949 0.6728
130.2888 0.4483 0.6083 380.2671 0.5418 0.6697
140.2827 0.4526 0.6155 390.2766 0.5608 0.6697
150.2563 0.4996 0.6610 400.2626 0.5377 0.6718
160.3460 0.4083 0.5413 410.4191 0.3691 0.4683
170.2649 0.5197 0.6624 420.5702 0.3416 0.3746
180.3701 0.3976 0.5179 430.3094 0.4306 0.5819
190.5360 0.3407 0.3886 440.2680 0.4682 0.6360
200.3483 0.4117 0.5417 450.2699 0.5567 0.6735
210.2708 0.4662 0.6326 460.2664 0.5252 0.6635
220.5021 0.3465 0.4083 470.3089 0.4319 0.5830
230.2542 0.5172 0.6705 480.2931 0.5986 0.6713
240.6215 0.3468 0.3581 490.2691 0.5495 0.6712
250.3269 0.4205 0.5626 500.4002 0.3805 0.4874
Table 5. The most satisfactory operation schedule.
Table 5. The most satisfactory operation schedule.
Line No.Departure TimeWaiting Time for Transfer Passengers (min)Number of Passengers Who Not Achieved Easy TransferFull Load RateR
17:04 7:12 7:21238.0460.8670.697
27:07 7:14 7:22 7:26139.0390.6200.714
37:10 7:19 7:25215.0290.6830.684
47:08 7:13 7:17 7:22 7:27117.4110.4930.712
Table 6. Comparison and analysis of schedules.
Table 6. Comparison and analysis of schedules.
Departure IntervalThe Total Waiting Cost of Passengers (Yuan)The Comprehensive Service Ratio of BusesNumber of Passengers Who Not Achieved Easy Transfer
Optimal departure interval2641.18 1.3760 125
3 min3578.25 1.1383 387
4 min3610.00 1.2105 413
5 min3703.40 1.2802 443
6 min3889.94 1.3432 462
7 min3834.57 1.4177 510
8 min3717.22 1.4981 533
9 min4041.29 1.5253 548
10 min4433.40 1.5511 594
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, X.; Qi, Y. Research on Optimization of Multi-Objective Regional Public Transportation Scheduling. Algorithms 2021, 14, 108. https://0-doi-org.brum.beds.ac.uk/10.3390/a14040108

AMA Style

Yang X, Qi Y. Research on Optimization of Multi-Objective Regional Public Transportation Scheduling. Algorithms. 2021; 14(4):108. https://0-doi-org.brum.beds.ac.uk/10.3390/a14040108

Chicago/Turabian Style

Yang, Xinfeng, and Yicheng Qi. 2021. "Research on Optimization of Multi-Objective Regional Public Transportation Scheduling" Algorithms 14, no. 4: 108. https://0-doi-org.brum.beds.ac.uk/10.3390/a14040108

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