Research Article  Open Access
A Practical Method for Timetable Rescheduling in Subway Networks during the EndofService Period
Abstract
This study proposes a biobjective optimization method for timetable rescheduling during the endofservice period of a subway network, taking all stakeholders’ interests into consideration. We seek to minimize the total transfer waiting time for all transfer passengers, meanwhile minimizing the deviation to the scheduled timetable. The εconstraint method and linearization techniques are utilized to obtain the approximate Pareto optimal solutions within limited seconds, allowing for figuring out the tradeoff between the two objectives. The method is validated by numerical experiments for different delay scenarios based on a realworld case: the Beijing subway network.
1. Introduction
Recently, there are several contributions to the lasttrain timetabling problem of a subway system, which focused only on the last train of each line in the network and expected to generate a more efficiently scheduled timetable for all last trains [1–4]. However, it is very common that the scheduled timetable cannot be implemented because of unavoidable delays that occur at an operational level [5]. As a result, this study is devoted to the timetable rescheduling problem during a specific period: the endofservice period.
Typically, most subway systems will be closed to the public at midnight or thereabouts for maintenance. Owing to the differences in passenger flow characteristics between different lines in a subway network, the operational time frames vary considerably among different lines. To be specific, the endofservice period in this study is defined as a period of time from the scheduled departure time of the earliest last train from its originating station (among all last trains of all lines) to the time when all trains finish their jobs.
Because of unavoidable disturbances in the daily operation, a lot of contributions have been made to the timetable rescheduling problem during other periods (e.g., peak hours) [6]. But there is a considerable difference between the endofservice period and other periods. For example, if a transfer passenger misses a connecting train due to a disturbance before the endofservice period, there is always a next connecting train which the transfer passenger can board. But during the endofservice period, the missed connecting train may be the last train of the connecting line; if so, the transfer passenger cannot finish his/ her trip via the subway system, which will bring a lot of inconvenience to the transfer passenger.
As a result, in order to deal with the disturbances occurring during the endofservice period, the first contribution of this study is that a timetable rescheduling model is proposed from a stakeholderoriented perspective with the consideration of benefits of both passengers and operating agencies. On the one hand, we seek to minimize the total transfer waiting time (TTWT) of all transfer passengers, and a penalty time is adopted if transfer passengers miss their last connecting trains, which benefits improving the level of service (LOS) after disturbances. On the other hand, we try to minimize the deviation between the rescheduled timetable and the scheduled timetable, which also benefits passengers who do not need to transfer.
In addition, in contrast to previous studies that focused only on the last train of each line, there are multiple trains running on each line during the endofservice period, which means that the train connection relationship in transfer stations becomes much more complicated. But, the timetable rescheduling is carried out through realtime adjustment of an existing schedule, with a consequent need for fast computation. In order to solve the practical problem of a largescale and complex network efficiently, we utilize the εconstraint method and some linearization techniques to convert the proposed model into an integer linear programming (ILP) model that can be solved by Cplex speedily, which constitutes the second contribution of this study.
The third contribution of this study is that a realworld case study of the Beijing subway network is presented to validate the effectiveness of the proposed method. Historical automatic fare collection (AFC) data of the Beijing subway system is available to obtain the number of transfer passengers of each connection, as an important input of our model. The approximate Pareto frontier is obtained by calculating the approximate Pareto optimal solutions, which helps us understand the tradeoff of the two objectives.
The rest of this study is organized as follows. Section 2 reviews some recent studies about timetable rescheduling and lasttrain timetabling. The stakeholderoriented model for timetable rescheduling during the endofservice period is proposed in Section 3. Section 4 presents the model conversion and the solution strategy. Some experiments based on a realworld case, the Beijing subway network, are carried out in Section 5. Section 6 draws some conclusions and future directions in brief.
2. Literature Review
The literature review presented in this section focuses on two aspects: timetable rescheduling and lasttrain timetabling. Some recent publications are reviewed below in detail.
2.1. Timetable Rescheduling
There are a lot of studies focusing on timetable rescheduling, which can be classified by disturbance or disruption, microscopic or macroscopic, and passengeroriented or trainoriented [6]. Various approaches have been developed in these previous studies.
From a trainoriented perspective, D’Ariano et al. [7] regarded the timetable rescheduling problem as a huge job shop scheduling problem with nostore constraints and modeled the problem with an alternative graph formulation to minimize the deviation from the scheduled timetable. They proposed a branch and bound algorithm with implication rules enabling the speed up of the computation. Törnquist and Persson [8] presented a MIP model to minimize a cost function based on train delays considering reordering and rerouting of trains. But for certain scenarios, it is difficult to find good solutions within seconds. Therefore, Krasemann [9] designed a greedy algorithm to quickly find a good solution by performing a depthfirst search. Dündar and Șahin [10] developed a genetic algorithm to minimize the total weighted delay. The algorithm could reduce total delay time by around half in comparison to an artificial neural network method developed to mimic the decision behavior of dispatchers.
From a passengeroriented perspective, since Schöbel [11] proposed the first MIP model for the delay management problem to minimize the total delay time of all passengers, the model has been further extended in Schöbel [12], Schachtebeck and Schöbel [13], and Dollevoet et al. [14, 15] by considering limited capacity of tracks, priority decisions, rerouting passengers, and limited capacity of stations. Binder et al. [16] proposed an ILP model with three objectives: the passenger satisfaction, the operational costs, and the deviation from the scheduled timetable. Strategies include canceling, delaying, rerouting the trains, and scheduling emergency trains.
More recently, there are several publications focusing on timetable rescheduling of a subway system. However, methods were mostly proposed at a singleline level. Xu et al. [17] modeled the problem as a discrete event model considering service balance performance of both directions on a doubletrack subway line. The model is expected to minimize the total delay time of all trains rather than all passengers. Gao et al. [18] integrated the skipstop pattern into the rescheduling model for a doubletrack subway line. An iterative algorithm was proposed to solve the model based on the model decomposition. Xu et al. [19] proposed a passengeroriented model for rescheduling on a subway line considering the limited train capacity. The delay time of alighting passengers and the penalty time of stranded passengers constitute the generalized delay time, which is expected to be minimized.
2.2. LastTrain Timetabling
An enormous amount of literatures contribute to the timetabling problem, like Caprara et al. [20], Zhou and Zhong [21], Lee and Chen [22], Cacchiani and Toth [5], and Yang et al. [23]. However, the lasttrain timetabling problem of a subway system is an emerging issue in recent years. However, all related publications focused only on the last train of each line. Zhou et al. [1] developed an optimization model to reduce passengers’ transfer waiting time for last trains and inaccessible passenger volume of all origindestination pairs. Coordinated departure times for all last trains are determined by the model. Kang et al. [2] established a programing model with adjustable running times and dwell times to obtain coordinated arrival and departure times of last trains at transfer stations. A genetic algorithm was designed to solve the model. Kang et al. [3] modeled the problem as a meanvariance model to improve the efficiency of transfer passengers. The model was solved by a genetic simulated annealing algorithm. Kang and Zhu [4] studied the same problem in Kang et al. [2] and designed a new heuristic algorithm outperforming both genetic algorithms and simulated annealing algorithms.
Based on all publications reviewed above, we present the focus of this study here. In case there is a disturbance occurring in a subway network during the endofservice period, this study is working to offer a practical method for timetable rescheduling from the stakeholderoriented perspective, at a macroscopic level. Disturbances (i.e., delays of 3 to 10 minutes) will not make passengers change their predetermined origin and destination stations and paths. To the best of our knowledge, this study is the first attempt on timetable rescheduling during the endofservice period in a subway network and real data from the AFC system is used as the model input. We want to figure out the tradeoff between different objectives and provide a method of decision support to dispatchers.
3. Model Formulation
3.1. Notations
Some necessary parameters are defined as follows:
: the set of subway lines in a network, , where is the total number of lines. Specifically, a doubletrack subway line in practice is considered as two oneway lines.
: the set of all stations on line , , where is the total number of stations on line .
: the set of trains still in operation on line during the endofservice period, , where is the last train of line .
: the minimum headway of line during the endofservice period.
: the minimum running time for trains running from station to on line , including additional time of train starting and braking at stations.
: the minimum dwell time for trains stopping at station on line .
: the planned arrival time for train at station on line .
: the planned departure time for train from station on line .
: the number of transfer passengers on train , who need to transfer from station on line to station on line . In reality, station and station are the same station with different serial numbers on different lines.
: the average time for transfer passengers walking from the platform of station on line to the platform of station on line . It is obvious that different people have different walking speed. For model simplification, the transfer walking time is assumed to be constant for passengers of the same transfer direction in this study.
The decision variables are presented as follows:
: the actual arrival time for train at station on line .
: the actual departure time for train from station on line .
: the waiting time of transfer passengers who are from train and transfer to station on line . This is a period of time from passengers reaching the platform of station to the actual departure time of train from station .
3.2. Train Connection Relationship
During the endofservice period, once a delay occurs in the network, train connections determined in the scheduled timetable might change. To describe connections between trains of different lines, the binary variable is introduced. indicates that transfer passengers on feeder train can transfer from station on line to station on line and board the connecting train successfully. Figures 1 and 2 depict the successful and failed transfer connections, respectively.
As a result, can be determined by the following formula:
3.3. Transfer Waiting Time and Penalty
In normal daytime operation, travelers waiting at a given station may be unable to take the first available train, e.g., if that train has no spare capacity for additional passengers (Schmöcker et al., 2011). However, it seems reasonable to suppose that demand during the endofservice period is generally low enough to permit the assumption that capacity is always available. As a result, all passengers are assumed to board the first arriving train after they reach the platform in this study. When there is a successful transfer connection, can be determined by the following formula:
But for transfer passengers, the last train of the connecting line is the last chance to finish their trips. If the connection to the last connecting train is broken, it will bring a lot of inconvenience to transfer passengers. To avoid this undesired phenomenon as much as possible, a penalty time is introduced in this study. If the missed connecting train is the last train of the connecting line, the transfer waiting time of these transfer passengers (i.e., failed transfer passengers, FTP) equals the penalty time; see the following formula:
In summary, the complete can be determined by the following formula:
3.4. Model Constraints
The rescheduling model is mainly subject to some operational requirements to ensure the safety of the operation and the feasibility of the rescheduled timetable.
3.4.1. Initial Delay
This constraint is to input the delay information (e.g., the delayed train, delay time, and position) to the model; see the following formula:where indicates the delayed train running on line . represents the station, where train is located when the disturbance occurs or the first station that is going to be visited by train after the disturbance occurs. represents the delay time.
3.4.2. Actual Arrival and Departure Time
During the process of rescheduling, the actual arrival and departure times of trains at stations cannot be earlier than the scheduled times; see the following formulas:
3.4.3. Section Running Time
Under the limitations of the traction and brake performance of trains, the length of each section, safety requirements, and the actual running times of trains in sections must be longer than the minimum running times [24]; see the following formula:
3.4.4. Dwell Time
Adjusting the dwell time is an important measure for dispatchers to control subway trains. Similar to the section running time, the actual dwell times of trains at stations must be longer than the minimum dwell times [24]; see the following formula:
3.4.5. Headway
As we mentioned above, there is more than one train still running on each line during the endofservice period. Thus, all trains running on each line should meet the requirements of minimum headway during the endofservice period; see the following formulas:
3.5. Optimization Objectives
We present two objectives to be optimized here. First, we try to minimize the total transfer waiting time (TTWT) for all transfer passengers, which helps to improve the LOS of the system after disturbances; see formula (12). The first objective also benefits increasing the number of successful transfer passengers after disturbances because of the penalty time set for failed transfer passengers (FTP) who miss their last connecting trains. Second, we seek to minimize the deviation between the rescheduled timetable and the scheduled timetable, which is usually dispatchers’ first goal in practice; see formula (13). Passengers who do not need to transfer after disturbances also benefit from the second objective.
The two objectives and constraints (1) and ((4)(11)) above consist of the complete rescheduling model during the endofservice period. The model is called a stakeholderoriented model because interests of both operation and passengers (transfer and nontransfer) are considered in the two objectives. The tradeoff between the two objectives may help dispatchers make decisions.
4. Model Solution
Owing to the huge complexity of the timetable rescheduling problem, especially when solving a realworld case of a largescale and complex network, many heuristic algorithms have been proposed to speed solving this problem. Examples include greedy algorithm [9], particle swarm algorithm [25], and genetic algorithm [10]. However, in this study, we utilize the εconstraint method [26] to convert the model into a singleobjective model, and some linearization techniques are adopted to reformulate all nonlinear constraints into linear constraints. Then, the original model is converted to a singleobjective ILP model, which can be solved by Cplex rapidly. By solving the problem for different values of , the approximate Pareto frontier can be shown to understand the tradeoff between the two objectives.
4.1. The εConstraint Method
The εconstraint method is good at solving multiobjective models to obtain a set of approximate Pareto optimal solutions. During the endofservice period, we should put transfer passengers on the first place because they may miss their last connecting trains due to delayed feeder trains. As a result, the objective function (13) is chosen as the εconstraint and reformulated by the following formula:where indicates the optimal objective value when only objective function (13) is considered in the model. is a coefficient representing dispatchers’ tolerability to the deviation between the rescheduled and scheduled timetables, .
By the εconstraint method, the original model is converted into a singleobjective model with the objective function (12) and constraints (1), ((4)(11)), and (14).
4.2. Linearization
Among all constraints, constraints (1) and (4) are nonlinear constraints. In order to speed the process of model solution, constraints (1) and (4) need to be reformulated into linear constraints.
is introduced to represent a big enough positive integer; then formula (1) can be easily replaced by the following linear formulas:
On the premise of the objective to minimize the total transfer waiting time of all transfer passengers, formula (4) can be relaxed into formula (17). Then formula (17) can be replaced by linear formulas ((18)(21)).
Finally, the singleobjective model obtained in Section 4.1 is converted into a singleobjective ILP model with objective function (12) and linear constraints ((5)(11)), ((14)(16)), and ((18)(21)), which can be solved by Cplex within limited seconds.
5. Case Study
To validate the method proposed in this study, the Beijing subway network is used as a realworld case study. By the end of 2016, the Beijing subway network consisted of 18 doubletrack lines (i.e., 36 oneway lines), 53 transfer stations, and 225 ordinary stations with an average daily ridership of 9.998 million passengers. A sketch map of the Beijing subway network without Airport Express is shown in Figure 3. Words with all letters in uppercase are acronyms of stations’ names.
Owing to the difference in passenger flow characteristics, different lines have different operational time frames. Among all last trains of all lines in Beijing subway network, the earliest one is the last train of Fangshan Line from SZ to GGZ, starting at 22:00. According to the definition in this study, the endofservice period of the Beijing subway case is from 22:00 to the time when all trains finish their jobs, a period of time about 2.5 hours. In addition, the starting time is (22:00) reset to 0 and then all times are changed according to the time lag and the minimum time unit is second.
Table 1 shows a sample of the real AFC data with key information of the Beijing subway system. The number of transfer passengers is the key to deciding whether a connecting train should wait for a delayed feeder train or depart on time. From the real data recorded in the AFC system, we can obtain the approximate number of transfer passengers during the endofservice period by a “PassengertoTrain” assignment method [27]. The AFC data, the data of average transfer walking time in all transfer stations, and the data of the scheduled timetable are all provided by the Beijing Rail Transit Control Center.

5.1. ScenarioBased Experiments
In order to prove that the proposed method is effective, various delay scenarios are generated randomly in terms of delayed train, delay position, and delay time. Numerical experiments based on these delay scenarios are carried out. Detailed information about these delay scenarios is listed in Table 2.

We test these delay scenarios with , minutes, and hour (i.e., 3600 seconds). All corresponding problems are solved within 2 seconds by Cplex 12.6.2 on a laptop computer with Intel Core i77700HQ CPU @ 2.8 GHz, 8 GB RAM. For benchmarking, we also test these delay scenarios with , which is similar to actual behaviors of dispatchers to minimize the deviation. Table 3 reports the solution results in detail.

The TTWT includes the transfer waiting time of all successful transfer passengers and the penalty time of failed transfer passengers. For most scenarios, there is a considerable decrease of about 20% in the TTWT as well as a big decline in the number of FTP, about 40% compared to those of the rescheduled timetable with , which is used to mimic dispatchers’ behaviors. These numerical results show that the LOS of subway systems after disturbances can be improved obviously by our method.
5.2. Delay Time Analysis
During the daily operation, different disturbances may lead to different delay times. In this experiment, we focus on Scenario 1 and set delay time changing from 5 to 10 minutes to test the effect of the method on different delay times. Similarly, all corresponding problems are solved within 2 seconds by Cplex 12.6.2 on the same computer. Table 4 shows the solution results in detail.

With the delay time increasing from 5 to 10 minutes, our method can reduce the TTWT by 17.91% to 23.44% compared with that of the rescheduled timetable with . In other words, the gap in TTWT between the rescheduled timetables by dispatchers and by our method is becoming bigger with the increasing delay time. As a result, with the improving requirements for LOS from passengers, our method is much better than dispatchers when tackling a disturbance which leads to a long delay time.
5.3. TradeOff between Objectives
Our proposed biobjective model for timetable rescheduling during the endofservice period aims to minimize the TTWT for all transfer passengers, meanwhile minimizing the deviation to the scheduled timetable. However, in the actual process of decisionmaking, it is difficult for dispatchers to obtain the optimal solution for multiple criteria. As a result, we are interested in the tradeoff between objectives and adopt the εconstraint method to obtain a set of approximate Pareto optimal solutions, allowing dispatchers to choose one solution simply.
Scenario 1 is still an example to obtain approximate Pareto optimal solutions by changing the value of ε from 0 to 0.5. The numerical results are shown in the left part of Figure 4. With ε rising from 0 to 0.5, the TTWT decreases rapidly; however, the deviation increases steadily. The right part of Figure 4 shows us the approximate Pareto frontier, which can demonstrate the tradeoff between the two objectives. There is an obvious trend between the two objectives: a decrease in the TTWT corresponding to an increase in the deviation. Dispatchers should weigh up interests of all stakeholders and make decisions with the best solution in their minds according to the actual situation.
In addition, for each approximate Pareto optimal solution, we calculate the TTWT, the number of FTP, and the total travel time (TT) of all trains involved in the endofservice period by , as shown in Table 5. According to the numerical results we can conclude that a high tolerability to the deviation only causes a rather small extension in the total TT of all involved trains. For example, among all rescheduled timetables, the biggest increment in total TT is only 132 seconds (), but there is a 27.82% decrement in the TTWT and a 58.13% decrement in the number of FTP compared to those of the rescheduled timetable with .

During the endofservice period, a passenger who does not need to transfer can catch his or her train definitely even if the train is late, but a transfer passenger may miss the last connecting train because of the late feeder train. As a result, based on all our numerical results in Section 5.3, Figure 4, and Table 5 in particular, we strongly suggest that dispatchers should take more interests from transfer passengers into consideration when rescheduling timetable after disturbances during the endofservice period.
6. Conclusions
The timetable rescheduling problem can be optimized by many objectives because of its inherently multicriterion nature. It is difficult to tell which solution is the optimal solution. But in terms of some specific criteria, we can figure out that a solution is better or worse than others. As a result, one of the major contributions of this study is that a biobjective optimization method is proposed from the stakeholderoriented perspective to tackle the timetable rescheduling problem during the endofservice period of a subway network, which allows us to figure out the tradeoff between the LOS (in terms of the TTWT and the number of FTP) and the operation (in terms of the deviation to the scheduled timetable). We utilize the εconstraint method to obtain approximate Pareto optimal solutions within limited seconds for different delay scenarios based on a realworld case, the Beijing subway network, which can help dispatchers to make decisions during the process of rescheduling.
In addition, given the actual characteristics of the endofservice period as well as the fact that a high tolerability to the deviation will not lead to a big extension in the total travel time, we think that dispatchers should put transfer passengers’ interests in the first place when rescheduling during the endofservice period, which will benefit the overall LOS after disturbances.
For future extension, the timetable rescheduling problem during the endofservice period can take passengers rerouting into consideration, especially when the delay time is long. However, train rerouting seems an unreasonable option because in most subway systems, trains belonging to a line cannot run on other lines. Anyway, contributions should be devoted to improving the LOS of subway systems after disturbances.
Data Availability
All related data are included within this article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (51478036) and the Fundamental Research Funds for the Central Universities (2018YJS074).
References
 W. Zhou, L. Deng, M. Xie, and X. Yang, “Coordination optimization of the first and last trains' departure time on urban rail transit network,” Advances in Mechanical Engineering, vol. 2013, Article ID 848292, 12 pages, 2013. View at: Publisher Site  Google Scholar
 L. Kang, J. Wu, H. Sun, X. Zhu, and Z. Gao, “A case study on the coordination of last trains for the Beijing subway network,” Transportation Research Part B: Methodological, vol. 72, pp. 112–127, 2015. View at: Publisher Site  Google Scholar
 L. Kang, X. Zhu, J. Wu, H. Sun, S. Siriya, and T. Kanokvate, “Departure time optimization of last trains in subway networks: meanvariance model and GSA algorithm,” Journal of Computing in Civil Engineering, vol. 29, no. 6, Article ID 04014081, 2014. View at: Publisher Site  Google Scholar
 L. Kang and X. Zhu, “Strategic timetable scheduling for last trains in urban railway transit networks,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 45, pp. 209–225, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 V. Cacchiani and P. Toth, “Nominal and robust train timetabling problems,” European Journal of Operational Research, vol. 219, no. 3, pp. 727–737, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 V. Cacchiani, D. Huisman, M. Kidd et al., “An overview of recovery models and algorithms for realtime railway rescheduling,” Transportation Research Part B: Methodological, vol. 63, pp. 15–37, 2014. View at: Publisher Site  Google Scholar
 A. D'Ariano, D. Pacciarelli, and M. Pranzo, “A branch and bound algorithm for scheduling trains in a railway network,” European Journal of Operational Research, vol. 183, no. 2, pp. 643–657, 2007. View at: Publisher Site  Google Scholar
 J. Törnquist and J. A. Persson, “Ntracked railway traffic rescheduling during disturbances,” Transportation Research Part B: Methodological, vol. 41, no. 3, pp. 342–362, 2007. View at: Publisher Site  Google Scholar
 J. Törnquist Krasemann, “Design of an effective algorithm for fast response to the rescheduling of railway traffic during disturbances,” Transportation Research Part C: Emerging Technologies, vol. 20, no. 1, pp. 62–78, 2012. View at: Publisher Site  Google Scholar
 S. Dündar and I. Şahin, “Train rescheduling with genetic algorithms and artificial neural networks for singletrack railways,” Transportation Research Part C: Emerging Technologies, vol. 27, no. 2, pp. 1–15, 2013. View at: Publisher Site  Google Scholar
 A. Schöbel, “A model for the delay management problem based on mixedintegerprogramming,” Electronic Notes in Theoretical Computer Science, vol. 50, no. 1, pp. 1–10, 2001. View at: Publisher Site  Google Scholar
 A. Schöbel, “Capacity constraints in delay management,” Public Transport, vol. 1, no. 2, pp. 135–154, 2009. View at: Publisher Site  Google Scholar
 M. Schachtebeck and A. Schöbel, “To wait or not to waitand who goes first? Delay management with priority decisions,” Transportation Science, vol. 44, no. 3, pp. 307–321, 2010. View at: Publisher Site  Google Scholar
 T. Dollevoet, D. Huisman, M. Schmidt, and A. Schöbel, “Delay management with rerouting of passengers,” Transportation Science, vol. 46, no. 1, pp. 74–89, 2012. View at: Publisher Site  Google Scholar
 T. Dollevoet, D. Huisman, L. Kroon, M. Schmidt, and A. Schöbel, “Delay management including capacities of stations,” Transportation Science, pp. 1–19, 2014. View at: Publisher Site  Google Scholar
 S. Binder, Y. Maknoon, and M. Bierlaire, “The multiobjective railway timetable rescheduling problem,” Transportation Research Part C: Emerging Technologies, vol. 78, pp. 78–94, 2017. View at: Publisher Site  Google Scholar
 X. Xu, K. Li, and L. Yang, “Rescheduling subway trains by a discrete event model considering service balance performance,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 40, no. 2, pp. 1446–1466, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Gao, L. Kroon, M. Schmidt, and L. Yang, “Rescheduling a metro line in an overcrowded situation after disruptions,” Transportation Research Part B: Methodological, vol. 93, pp. 425–449, 2016. View at: Publisher Site  Google Scholar
 W. Xu, P. Zhao, and L. Ning, “A passengeroriented model for train rescheduling on an urban rail transit line considering train capacity constraint,” Mathematical Problems in Engineering, vol. 2017, Article ID 1010745, 9 pages, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 A. Caprara, M. Fischetti, and P. Toth, “Modeling and solving the train timetabling problem,” Operations Research, vol. 50, no. 5, pp. 851–861, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 X. Zhou and M. Zhong, “Singletrack train timetabling with guaranteed optimality: branchandbound algorithms with enhanced lower bounds,” Transportation Research Part B: Methodological, vol. 41, no. 3, pp. 320–341, 2007. View at: Publisher Site  Google Scholar
 Y. Lee and C. Chen, “A heuristic for the train pathing and timetabling problem,” Transportation Research Part B: Methodological, vol. 43, no. 89, pp. 837–851, 2009. View at: Publisher Site  Google Scholar
 X. Yang, X. Li, Z. Gao, H. Wang, and T. Tang, “A cooperative scheduling model for timetable optimization in subway systems,” IEEE Transactions on Intelligent Transportation Systems, vol. 14, no. 1, pp. 438–447, 2013. View at: Publisher Site  Google Scholar
 R. C. W. Wong, T. W. Y. Yuen, K. W. Fung, and J. M. Y. Leung, “Optimizing timetable synchronization for rail mass transit,” Transportation Science, vol. 42, no. 1, pp. 57–69, 2008. View at: Publisher Site  Google Scholar
 X. Meng, L. Jia, and Y. Qin, “Train timetable optimizing and rescheduling based on improved particle swarm algorithm,” Transportation Research Record, no. 2197, pp. 71–79, 2010. View at: Publisher Site  Google Scholar
 O. J. IbarraRojas, R. Giesen, and Y. A. RiosSolis, “An integrated approach for timetabling and vehicle scheduling problems to analyze the tradeoff between level of service and operating costs of transit networks,” Transportation Research Part B: Methodological, vol. 70, pp. 35–46, 2014. View at: Publisher Site  Google Scholar
 W. Zhu, W. Wang, and Z. Huang, “Estimating Train Choices of Rail Transit Passengers with Real Timetable and Automatic Fare Collection Data,” Journal of Advanced Transportation, vol. 2017, Article ID 5824051, 12 pages, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Wenkai Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.