Research Article  Open Access
Shuo Yang, Kai Yang, Ziyou Gao, Lixing Yang, Jungang Shi, "LastTrain Timetabling under Transfer Demand Uncertainty: MeanVariance Model and Heuristic Solution", Journal of Advanced Transportation, vol. 2017, Article ID 5095021, 13 pages, 2017. https://doi.org/10.1155/2017/5095021
LastTrain Timetabling under Transfer Demand Uncertainty: MeanVariance Model and Heuristic Solution
Abstract
Traditional models of timetable generation for last trains do not account for the fact that decisionmaker (DM) often incorporates transfer demand variability within his/her decisionmaking process. This study aims to develop such a model with particular consideration of the decisionmakers’ risk preferences in subway systems under uncertainty. First, we formulate an optimization model for lasttrain timetabling based on meanvariance (MV) theory that explicitly considers two significant factors including the number of successful transfer passengers and the running time of last trains. Then, we add the meanvariance risk measure into the model to generate timetables by adjusting the last trains’ departure times and running times for each line. Furthermore, we normalize two heterogeneous terms of the risk measure to provide assistance in getting reasonable results. Due to the complexity of MV model, we design a tabu search (TS) algorithm with specifically designed operators to solve the proposed timetabling problem. Through computational experiments involving the Beijing subway system, we demonstrate the computational efficiency of the proposed MV model and the heuristic approach.
1. Introduction
The subway is characterized by a high level of safety, wide accessibility, and high energy performance and offers reliable service with sufficient punctuality [1]. Generally, the subway does not provide passengers with 24hour service, which means that train operations stop at some point each day (typically late at night or in the very early morning hours) for track or fleet maintenance. The last train naturally becomes the last chance for the passengers transported by the subway to reach their destinations. In addition, the subway usually does not provide a pointtopoint transportation service, and the passengers are usually required to transfer one or several times within the network. Statistics show that the lasttrain transfer passengers accounted for more than 65% of the total lasttrain passengers in the entire Beijing subway network in 2015. From the perspective of passengers, they strongly desire welldesigned lasttrain timetables, with good coordination between last feeder trains and last connecting trains so that they can transfer smoothly. Therefore, lasttrain timetabling should be determined to satisfy the transfer demands.
In many practical situations, the transfer demands of last trains are stochastic. There are several reasons for this, such as inherent stochasticity of realworld traffic flow evolution, multiple sources of data, and fluctuating nature of parameter values. To tackle such type of uncertainty, a samplebased representation may be an appropriate research tool. Based on the historical data, this representation collects the dayspecific data, used as stochastic samples with specified probabilities, to capture the characteristics of randomness associated with the transfer demands in subway networks. Furthermore, the scenariospecific samples enable a flexible way for decisionmaker (DM) to formulate a model for lasttrain timetabling problem under uncertainty.
The uncertainty of transfer demands can also bring risk to the DM in designing lasttrain timetables process. The formulation of meanvariance theory is a basic optimization tool to model risk concept with regard to a quadratic utility function for the DM’s risk preference. In other words, the DM would want to generate lasttrain timetables with a higher expected value and a lower variance of positive utilities, such as the number of successful transfers in a risky situation. With these concerns, we are particularly interested in developing a systematic approach that synchronously considers stochastic transfer demands characteristics and management of risks in a lasttrain timetabling problem.
By combining these aforementioned aspects, we intend to address the lasttrain timetabling problem in uncertain conditions of transfer demands in this paper. For modeling uncertainty, we first employ the meanvariance method to formulate an optimization model for lasttrain timetabling. Subsequently, we embed the meanvariance risk measure into the model to get timetables for last trains, which is the main novelty of this study compared to previous works. Specifically, this paper aims to make the following contributions to the study of lasttrain timetable problems.(i)We adopt the representation of samplebased transfer demands to capture the characteristics of randomness in a subway network. By doing so, the measured transfer demands from realworld applications can be directly embedded in the proposed modeling framework. Furthermore, scenariospecific samples enable a flexible way to measure the uncertainty in decision process, compared to probability function distributions that typically only represent an aggregated version of transfer demands randomness.(ii)We formulate the meanvariance utility criterion to fully capture the risk preference under uncertainty. With the data structure, we rigorously propose a meanvariance model for lasttrain timetabling problem that increases the number of successful transfer passengers and at the same time reduces the total running time of last trains. Compared with existing researches, the proposed lasttrain timetabling framework simultaneously considers uncertainty in the transfer demands and risk modeling.(iii)We provide a highly efficient method framework to solve the proposed model. In order to eliminate the effect of the heterogeneity between the variance and the expected value in the objective function, we firstly suggest a normalization method by normalizing two heterogeneous terms of the risk measure. Then, we further design a TS algorithm with specifically designed operators to solve the proposed model. The proposed approaches provide a promising algorithmic foundation for solving other closely related lasttrain timetabling problems in reallife subway networks.(iv)We conduct a case study based on the Beijing subway instances to demonstrate the effectiveness of the proposed model and solution method. Compared to its counterpart with expected utility criteria, the computational results reported in the numerical experiments demonstrate the superiority of the proposed approach. Further analysis shows that the DM can select distinct values of riskaversion coefficients depending on their preferences in our optimization method.
The rest of this paper is organized as follows. In Section 2, we briefly review related literature. In Section 3, we give a detailed description for the problem and state the formulation of lasttrain timetabling problem based on the meanvariance theory. Then, we propose the TS heuristic algorithm to solve the developed model in Section 4. In Section 5, we design a case study based on the realworld operation data of Beijing subway, to demonstrate the effectiveness of the proposed approaches. Finally, some conclusions and further studies are presented in Section 6.
2. Literature Review
The literature relevant to this paper comes from two streams: the train timetabling literature and the risk analysis and management literature.
The first stream studies the train timetabling problems. Peeters and Kroon [2] presented a mixedinteger nonlinear programming formulation for cyclic railway timetabling, where the integer variables corresponded to cycles in the graph induced by the constraints. Albrecht and Oettich [3] proposed an algorithm for the dynamic modification of train running times to increase the probability of making connections to other means of public transport. Liebchen [4] formulated a timetable optimization model that optimized the arrival and departure times at the transfer stations of the metro network in Berlin so that the passenger transfer waiting times were minimized. Stanojević et al. [5] developed an integer programming approach for determining an optimal timetable for a single, oneway track linking two major stations. Yin et al. [6] established linear mixedinteger programming models for timetable optimization to decrease the energy consumption and simultaneously to shorten the passengers’ waiting time. Gao et al. [7] investigated a rescheduling problem on a metro line in an overcrowded situation after disruptions. Huang et al. [8] studied a joint train scheduling optimization with service quality and energy efficiency in urban rail transit networks, which was solved by a heuristic.
Addressing the lasttrain timetabling problem in subway networks is a new branch in the train timetabling problem domain. Zhou et al. [9] built coordination optimization models of the last trains’ departure time to reduce passengers’ transfer waiting time for the last trains and inaccessible passenger volume of all origindestination (OD) and improve passengers’ accessible reliability for the last trains. Kang et al. [10] proposed a rescheduling model for last trains with the consideration of train delays caused by incidents that occurred in train operations. Kang et al. [11] established a lasttrain network transfer model to maximize passenger transfer connection headway, which reflected lasttrain connections and transfer waiting time. Li et al. [12] studied the lasttrain coordination problem for metro networks based on actual passenger demand obtained from AFC system data, aiming to maximize the total number of passengers who can reach their destinations by metro prior to the end of operation.
In addition, there is other research concerned with the train timetabling problem under uncertainty. Yang et al. [13] formulated an expected value goalprogramming model for the train timetable problem with fuzzy passenger demand on a singleline railway, in which two objectives, that is, fuzzy total passengers’ time and total delay time, were considered. Barrena et al. [14] studied the design and optimization of train timetables for a rail rapid transit line to minimize the average passenger waiting time in a dynamic demand environment, in which an adaptive large neighborhood search metaheuristic was proposed to solve large instances of the problem. Yin et al. [15] developed a stochastic programming model for metro train rescheduling problem in order to jointly reduce the time delay of affected passengers, their total traveling time, and operational costs of trains. Yin et al. [16] converted the train operation problem into a Markov decision process with nondeterministic state transition probabilities to minimize the cost for both the total time delay and energy consumption in a subway line.
The second stream of literature related to this paper is on risk analysis and management. Choi et al. [17] presented a risk assessment methodology for underground construction projects with an uncertainty model based on fuzzy concept. Hensher et al. [18] used a nonlinear probability weighting function to accommodate choice made under risk in the context of the expected utility framework. Woodward et al. [19] designed the flood risk strategies to allow for flexible adaptive measures by capturing the concepts of real options and multiobjective optimization to evaluate potential flood risk management opportunities. Zio [20] examined a framework of vulnerability and risk analysis in relation to its application for the protection and resilience of critical infrastructure. Liu et al. [21] presented a risk analysis of derailment rates on Class I railroad mainlines based on data from the US Federal Railroad Administration and the major freight railroads.
In regard to meanvariance risk measure, Markowitz [22] firstly introduced the meanvariance model into financial portfolio management. Choi et al. [23] used the meanvariance theory to quantify the risk and the payoff of an investment decision by the variance of return and the expected return, respectively. Fu et al. [24] investigated a continuoustime meanvariance portfolio selection problem under the constraint of a borrowing rate which is higher than the lending rate. Wu et al. [25] presented a meanvariance optimization model for solving the track allocation problem that minimizes the occupation time costs in groups of turnouts at station bottlenecks. Ahmadi et al. [26] developed a systematic approach to make robust decisions for the single locationallocation hub median problem based on meanvariance theory. Li et al. [27] presented a systematic approach integrated the meanvariance and scheduling to allow for schedule delay and trip time variability under uncertainty.
To the best of our knowledge, there has been no significant attempt to study the lasttrain timetabling problems in stochastic environments considering reliability of decisions; even stochastic models of this problem without considering risk are quiet rare in the literature. This paper takes the first initiative to develop a systematic risk management approach to make decisions according to decisionmakers’ preferences for a lasttrain timetabling problem. In this study, more comprehensive and realistic factors are taken into account; these features make this study different from the literature in lasttrain timetabling area.
3. Problem Description and Formulations
Lasttrain timetabling problem aims at coordinating the arrival and departure times of feeder trains with connecting trains at transfer stations so that passengers can enjoy smooth transfers in subway networks. The practical transfer situations based on the arrival and departure times of the last feeder train and the last connecting train is illustrated in Figure 1. Case A depicts that the passengers can smoothly transfer from train to train since the departure time of the train is later than the passengers’ arrival time at the platform of train . Conversely, the transfer passengers cannot do that before the train departs, as shown in the case B. That is, the passengers fail to transfer in this case. Therefore, lasttrain timetabling should be determined to satisfy the transfer demands. With these concerns, we generate timetables to pick up as many transfers as possible by adjusting the last trains’ departure times and running times for each line.
Next, we first present the basic setting upon which our modeling framework is built. We also present several underlying assumptions and most of the notation considered throughout the paper. Then, we propose a modeling framework that synchronously considers uncertainty in the transfer demand and risk modeling, which will be discussed in the following subsections.
3.1. Basic Model
In this section, we introduce the basic setting for our analysis. We consider a simple base situation in order to get a more focused demonstration of the methodology we propose.
To provide a precise statement of this problem, we define the following:
Sets : the set of last trains in the network, . : the set of all transfer stations on train , . : the set of transfer stations up to station on train . : the set of transfer stations before station on train .
Parameters : the passengers’ average walking time from train to at station . : the transfer demand from trains to at station . : the dwell time of the last train at station . : the minimum running time on train between station and the previous transfer station, including additional time caused by train start and brake at stations. : the maximum running time on train between station and the previous transfer station, including additional time caused by train start and brake at stations. : the earliest permitted departure time for last trains. : the last permitted departure time for last trains.
Decision Variables : the departure time of the last train from the starting station. : the arrival time of the last train at transfer station . : the departure time of the last train from transfer station . : the segment running time of the last train between station and the previous transfer station, including additional time caused by train start and brake at station. : the delay on the platform of for passengers who transfer from to at station . : a binary variable which takes a value of 1 if passengers successfully transfer from train to at station and a value of 0 otherwise.
Assumptions 1. To facilitate the presentation of the essential ideas without loss of generality, the following basic assumptions are made in this study:(A1) The passengers entering the system during the lasttrain period proceed directly to trains and do not linger or engage in nontravel activities (purchasing coffee, etc.).(A2) The passengers’ average walking time between two lines at a station is fixed. With fixed passenger’s average walking time, the model can be simplified significantly.(A3) The last trains’ dwelling time at transfer stations is fixed with the same value. The simplicity makes the analysis of processing data more convenient.(A4) The capacity of each last train is sufficient to receive all passengers. The assumption is reasonable because the passenger volume in last trains is smaller than the train capacity.
Objective Functions. This part analyzes two objectives in the problem of interest, that is, the number of successful transfer passengers and the running time of all the last trains. Besides, we are interested in formulating the objectives into linear formulas that could make the model easier to be handled and also provide a basis for subsequent heuristic solution algorithm.
(i) Total of Running Time. One of the most significant indicators in lasttrain timetabling is the total of running time for last trains in the view of the DM. The sum of last trains’ running time is represented as which is introduced to save the operational time and the passengers’ traveling time, thus avoiding operational costs and energy wastage.
(ii) Number of Successful Transfer Passengers. The other significant indicator is the number of successful transfer passengers by comparing the arrival time of feeder trains and the departure time of connecting trains. The number of successful transfer passengers can be calculated as where indicates that the passengers cannot make the change from to .
To obtain a tradeoff solution, we give two weight factors, that is, and , to, respectively, indicate the importance of each objective. The lasttrain timetabling problem can be expressed by the following mathematical program, which we call the basic model:
Constraints (4) and (5) are operational constraints. , constraint (4) is given to tracks the arrival time, , of the last train at station , where denotes the total segment running time of the last train from the starting station to station , and means the total dwell time along the way to station . Constraint (5) states that the departure time of train at station , , can be obtained by adding the dwell time at station to , in which the set contains transfer stations passed through by the last train including the current station , and the set contains members of except station , since the arrival time at a certain station has nothing to do with how long it will stay here later.
Constraints (6) and (7) are qualitative constraints. In constraint (6), represents the difference between the departure time of the connecting train and the arrival time of transfer passengers from the feeder train at the platform of train . It is determined by the departure time of the connecting train , , the arrival time of the feeder train , , and the passengers’ walking time from platform of to the one of , . In constraint (7), a binary variable is introduced to indicate whether passengers can successfully transfer from to at station . Here, is a sufficiently large positive number. That is, if the passengers succeed in transferring from to at station ; , otherwise.
Constraints (8) and (9) are timerange constraints. There are upper and lower bounds on segment running time for each interval in the subway network [28], and the departure times of last trains at starting stations are also supposed to be in a limited range for all trains , stations .
3.2. MeanVariance Model
3.2.1. Uncertainty in the Transfer Demands
As it is mentioned in the introductory section, stochastic variability (randomness) always exists in the lasttrain timetabling problems. In this study, we use the samplebased representation to characterize the randomness associated with the transfer demands. Specifically, we collect a set of day transfer demands during the lasttrain period (between 22:00 and 24:00 daily) acquired from the automated fare collection (AFC) system, used as stochastic samples. Each day corresponds to one sample with the even probability . Alternatively, if we use the notation sets to represent the considered days and their probabilities, then the probability of each dayspecified data set is deduced by .
Note that, we only consider discrete type of uncertainty where uncertainty is represented by a finite set of samples. When this is not the situation, one can resort to an approximate sampling technique [29, 30] and reduce the problem to the finitescenario setting that we address. Hereafter, we are intended to use the samplebased representation method proposed above to generate the corresponding random variables and then present an optimization model for lasttrain timetabling problem based on the meanvariance theory.
3.2.2. MeanVariance Theory
In a risky situation, the DM exhibits three different behaviors in the forms of riskseeking, riskneutral, and riskaverse that originated from concavity, linearity, and convexity of their utility function [26], respectively, as shown in Figure 2. Apparently, a rational DM prefers to choose a policy with a higher expected value and a lower variance of positive utilities. And in general, the DM has a riskaversion behavior to deal with a risky situation. Therefore, in this study, we assume that the DM is riskaverse and trying to make decisions with less risk.
The formulation of meanvariance theory is a basic optimization tool to model risk concept with regard to a quadratic utility function for the DM’s preference. In this situation, a tradeoff between risk and utility function can be formulated by using a riskaversion coefficient , given as follows:in which riskseeking, riskneutral, and riskaverse behavior is adjusted based on the negative, zero, and positive value for riskaversion coefficient , respectively. Therefore, in the stochastic situation with random variable , the quadratic utility function (positive utility) with a riskaverse preference can be derived aswhere and represent the expected value and variance of the function , respectively.
3.2.3. MeanVariance Risk Model
Due to the stochastic trait of transfer demands, it is necessary to have an especial tool to analyze the lasttrain timetabling problem and generate timetables. Moreover, the adopted timetable should be able to meet the subway system requirements for a long time, since changing the implemented initial timetable would be very costly. This clarifies the importance of obtaining solutions with higher reliability or lower risk in the lasttrain timetabling. Based on this consideration, we present a novel approach, considering risk management in the lasttrain timetabling problem under uncertainty. In addition to the parameters introduced in the previous section, we now define the following: : the index of scenarios, . : the number of passengers transferring from trains to at station under scenario . : the probability of scenario .
As stated above, we can reformulate the objective function with stochastic transfer demands as follows:where
To obtain the expected value and variance of the above function , the basic terms are introduced as
By applying the main concept of the first and second moment and substituting (14), the expected value and variance of are derived as (15) and (16), respectively.where and are the expected value and variance mathematical operators.
Therefore, the quadratic utility function based on meanvariance theory could be written as follows:where is the riskaversion coefficient for the meanvariance function that takes different values depending on the DM’s preference in dealing with risky conditions.
The heterogeneity between the variance and the expected value of the function leads to a difficulty in computing the quadratic utility function. It can be seen that greater values of variance compared to the expected value will lead to loss of sensitivity of the model in different scenarios. In order to reduce this deficiency, the objective function is normalized as suggested in multiobjective optimization approaches. In this research, the two functions of and are individually optimized in both minimum and maximum states and then the risk measure is defined as the objective function of risk modeling, which can be minimized as follows:
On the basis of aforementioned analysis, we can formally present the corresponding meanvariance (MV) model for lasttrain timetabling problem:
It is clear that, the MV model (19) is quite complicated because it involves a variety of constraints and contains many variables. IbarraRojas and RiosSolis [31] have shown that the train timetabling problem belongs to the NPhard class problem set because of the large number of variables and constraints. In order to obtain a satisfactory lasttrain timetable within an acceptable time, we design a TS algorithm with specifically designed operators to solve the MV model (19) in this paper.
4. Tabu Search Algorithm
For train timetabling problems, researchers usually use a very large number of different techniques, including exact and approximation algorithms. Among these solution approaches, the TS algorithm developed by Glover [32] has been very successful in finding close to optimal solutions for engineering problems [33–37]. The TS algorithm starts searching with an initial solution, based on which a set of feasible solutions is constructed as the neighborhood. Then the procedure will select the best local candidate among the neighborhood (although it has worse fitness than the current best) as the next current solution in order to escape the local optimal. The history record, that is, the tabu list, is modified after it. The process is repeated until the userspecified stopping criteria are satisfied. Finally, the best solution during the search process will be returned.
First, we introduce some notations. Let Initial_s be an initial solution and formed as an array composed by two vectors representing the departure time of last trains from origin stations and the segment running times, as shown in Figure 3. In this paper, we input original values of the two vectors as an initial solution. Given a solution Present_s, we define the set to be the set of neighboring solutions. To obtain the neighboring solutions, we applied the procedure illustrated in Figure 4.
Step 1 (selection). Randomly select positions in an array to make changes based on Present_s.
Step 2 (change). Produce values at the picked positions based on the Solution Generation Algorithm [38], while meeting constraints (8) and (9).
Step 3 (substitute). A new solution is generated by replacing the corresponding values with the newly generated ones in Present_s.
Based on the aforementioned analyses, the procedure of this heuristic algorithm is summarized by Algorithm 1.

5. Case Study
In this section, we verify the effectiveness of the TS algorithm and the application of the MV model using real cases from the Beijing subway system. The proposed TS algorithm has been coded in MATLAB 2013a. All numerical tests are carried out on a Windows Server 2008 R2 Enterprise workstation with Intel®Xeon® CPU E52640v3 @ 2.60 GHz and 64 G RAM to improve the solving efficiency.
Next, we perform this case study on the practical Beijing subway network. Since this paper focuses only on transfer between different lines, we simplify the original Beijing Subway Map for computational convenience. In detail, we give a simplified structure of the considered Beijing subway network in Figure 5, which includes 12 lines, 45 transfer stations, and 81 segments. Note that symbol different last trains with directions suggested by the arrows next to them.
5.1. Data Collection
Before conducting the experiments, we need to precollect some basic data in the solution process. More specifically, we adopt tenday transfer demands (i.e., 10 samples) during the lasttrain period (between 22:00 and 24:00 daily) as the experimental data. In addition, we set the probability of each sample as 0.1, which is described in Section 3.2.1. Tables 1 and 2 demonstrate the present timetable of last trains in Beijing subway. Note that 22:00 is set as the zero point for simplicity in the process and 50 means the 22:50 in 24 hour clock, to take an example.


5.2. Computational Results
After implementing the proposed model and solution approach, we examine the effect of varying riskaversion coefficients on the expected value of and its standard deviation, as summarized in Figures 6 and 7, respectively. As anticipated, the expected value of objective function increases as the riskaversion coefficient increases, while the standard deviation of objective function takes an opposite tendency with respective to riskaversion coefficient. Furthermore, it is obvious that there are unique decisions for each level of riskaversion coefficient.
Figure 8 demonstrates the threedimensional characteristics of the relationship among , , and . This curve contains all the different preferences of decisionmakers. In other words, there will be a riskaversion coefficient for every preference of the DM () and each riskaversion coefficient corresponds to a certain preference and then the curve consists of all the riskaversion DMs’ preferences. The value of riskaversion coefficient is chosen based on the DM’s preferences and his/her risk attitude, and it may be different for each DM. Therefore, the DM should be able to select a riskaversion coefficient on the curve based on their preference to obtain a reliable decision for implementation.
5.3. Comparison Results
To further validate the advantages of the proposed MV model, more computational results were obtained to solve a similar problem with expected utility criteria. When , the MV model is degenerated to expectation utility model. Table 3 summarizes the comparison of solution results on different riskaversion coefficient . It can be observed that the variance of objective function in expectation utility model is higher than that in MV model in terms of optimal solutions. This result suggests that the MV model obtains better results compared to its counterpart with expected utility criteria in a risky situation. Furthermore, the MV model can also provide comprehensive strategies to generate the lasttrain timetables according to the DM’ preferences.

6. Conclusions
Focusing on generating timetables for last trains, this paper proposed a modeling framework in subway systems with uncertainty. Specifically, the samplebased transfer demands were adopted to represent the randomness in the subway network. Considering the data structure, the meanvariance utility based model for lasttrain timetabling problems was formulated to fully capture the risk preference under uncertainty. Then, the meanvariance risk measure was added to the proposed model. One of the main problems associated with this measure is heterogeneity between the variance and the expected value in the objective function. In order to eliminate the effect of this heterogeneity, a normalization method was employed, which is of huge assistance in getting reasonable results. To solve the models more efficiently, a TS based algorithmic framework was proposed that can obtain a good solution in an acceptable computational time. Finally, a practical case study of the Beijing subway system demonstrated the effectiveness of the proposed model and optimization method.
Future research can focus on the following several aspects. (1) The improvements to the algorithm efficiency, designing hybrid algorithms, or developing an adaptive method can be considered in the further study. (2) The influence of feeder train delay needs to be studied in advance. (3) Except for the meanvariation criterion when coping with riskaversion preferences, we should carry out some other methods and investigate which can better signify risk levels in our future research.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (nos. 71422002, 71621001, 71631007, and 71701012), the Research Foundation of the State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University (no. RCS2017ZZ01), and the Fundamental Research Funds of Beijing Jiaotong University (no. 2015RC075).
References
 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
 L. Peeters and L. Kroon, “A Cycle Based Optimization Model for the Cyclic Railway Timetabling Problem,” in ComputerAided Scheduling of Public Transport, vol. 505 of Lecture Notes in Economics and Mathematical Systems, pp. 275–296, Springer Berlin Heidelberg, Berlin, Heidelberg, 2001. View at: Publisher Site  Google Scholar
 T. Albrecht and S. Oettich, “A new integrated approach to dynamic schedule synchronization and energysaving train control,” in Computers in Railways VIII, pp. 847–856, WIT Press, Southampton, UK, 2002. View at: Google Scholar
 C. Liebchen, “The first optimized railway timetable in practice,” Transportation Science, vol. 42, no. 4, pp. 420–435, 2008. View at: Publisher Site  Google Scholar
 P. Stanojević, M. Marić, J. Kratica, N. Bojović, and M. Milenković, “Mathematical optimization for the train timetabling problem,” Mathematica Balkanica New Series, vol. 24, no. 3, pp. 303–312, 2010. View at: Google Scholar
 J. Yin, L. Yang, T. Tang, Z. Gao, and B. Ran, “Dynamic passenger demand oriented metro train scheduling with energyefficiency and waiting time minimization: Mixedinteger linear programming approaches,” Transportation Research Part B: Methodological, vol. 97, pp. 182–213, 2017. View at: Publisher Site  Google Scholar
 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
 Y. Huang, L. Yang, T. Tang, Z. Gao, and F. Cao, “Joint train scheduling optimization with service quality and energy efficiency in urban rail transit networks,” Energy, vol. 138, pp. 1124–1147, 2017. View at: Publisher Site  Google Scholar
 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. 5, no. 3, pp. 848292–848292, 2013. View at: Publisher Site  Google Scholar
 L. Kang, J. Wu, H. Sun, X. Zhu, and B. Wang, “A practical model for last train rescheduling with train delay in urban railway transit networks,” Omega, vol. 50, pp. 29–42, 2015. 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
 W. Li, R. Xu, Q. Luo, and S. Jones, “Coordination of last train transfers using automated fare collection (AFC) system data,” Journal of Advanced Transportation, vol. 50, no. 8, pp. 2209–2225, 2016. View at: Publisher Site  Google Scholar
 L. X. Yang, K. P. Li, and Z. Y. Gao, “Train timetable problem on a singleline railway with fuzzy passenger demand,” IEEE Transactions on Fuzzy Systems, vol. 17, no. 3, pp. 617–629, 2009. View at: Publisher Site  Google Scholar
 E. Barrena, D. Canca, L. C. Coelho, and G. Laporte, “Singleline rail rapid transit timetabling under dynamic passenger demand,” Transportation Research Part B: Methodological, vol. 70, pp. 134–150, 2014. View at: Publisher Site  Google Scholar
 J. Yin, T. Tang, L. Yang, Z. Gao, and B. Ran, “Energyefficient metro train rescheduling with uncertain timevariant passenger demands: An approximate dynamic programming approach,” Transportation Research Part B: Methodological, vol. 91, pp. 178–210, 2016. View at: Publisher Site  Google Scholar
 J. Yin, D. Chen, L. Yang, T. Tang, and B. Ran, “Efficient RealTime Train Operation Algorithms with Uncertain Passenger Demands,” IEEE Transactions on Intelligent Transportation Systems, vol. 17, no. 9, pp. 1–13, 2016. View at: Publisher Site  Google Scholar
 H.H. Choi, H.N. Cho, and J. W. Seo, “Risk assessment methodology for underground construction projects,” Journal of Construction Engineering and Management, vol. 130, no. 2, pp. 258–272, 2004. View at: Publisher Site  Google Scholar
 D. A. Hensher, Z. Li, and J. M. Rose, “Accommodating risk in the valuation of expected travel time savings,” Journal of Advanced Transportation, vol. 47, no. 2, pp. 206–224, 2013. View at: Publisher Site  Google Scholar
 M. Woodward, Z. Kapelan, and B. Gouldby, “Adaptive flood risk management under climate change uncertainty using real options and optimization,” Risk Analysis, vol. 34, no. 1, pp. 75–92, 2014. View at: Publisher Site  Google Scholar
 E. Zio, “Challenges in the vulnerability and risk analysis of critical infrastructures,” Reliability Engineering and System Safety, vol. 152, pp. 137–150, 2016. View at: Publisher Site  Google Scholar
 X. Liu, M. Rapik Saat, and C. P. L. Barkan, “Freighttrain derailment rates for railroad safety and risk analysis,” Accident Analysis and Prevention, vol. 98, pp. 1–9, 2017. View at: Publisher Site  Google Scholar
 H. Markowitz, “Portfolio selection,” The Journal of Finance, vol. 7, no. 1, pp. 77–91, 1952. View at: Publisher Site  Google Scholar
 T.M. Choi, D. Li, H. Yan, and C.H. Chiu, “Channel coordination in supply chains with agents having meanvariance objectives,” Omega, vol. 36, no. 4, pp. 565–576, 2008. View at: Publisher Site  Google Scholar
 C. Fu, A. LariLavassani, and X. Li, “Dynamic meanvariance portfolio selection with borrowing constraint,” European Journal of Operational Research, vol. 200, no. 1, pp. 312–319, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 J. J. Wu, L. J. Kang, H. J. Sun, and X. L. Jia, “Track allocation optimization in railway station: meanvariance model and case study,” Journal of Transportation Engineering, vol. 139, no. 5, pp. 540–547, 2013. View at: Publisher Site  Google Scholar
 T. Ahmadi, H. Karimi, H. Davoudpour, and S. A. Hosseinijou, “A robust decisionmaking approach for phub median location problems based on twostage stochastic programming and meanvariance theory: a real case study,” International Journal of Advanced Manufacturing Technology, vol. 77, no. 912, pp. 1943–1953, 2015. View at: Publisher Site  Google Scholar
 H. Li, H. Tu, and D. A. Hensher, “Integrating the meanvariance and scheduling approaches to allow for schedule delay and trip time variability under uncertainty,” Transportation Research Part A: Policy and Practice, vol. 89, pp. 151–163, 2016. 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
 J. R. Birge and F. Louveaux, Introduction to Stochastic Programming, Springer, 1997. View at: MathSciNet
 A. J. Kleywegt, A. Shapiro, and T. HomemdeMello, “The sample average approximation method for stochastic discrete optimization,” SIAM Journal on Optimization, vol. 12, no. 2, pp. 479–502, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 O. J. IbarraRojas and Y. A. RiosSolis, “Synchronization of bus timetabling,” Transportation Research Part B: Methodological, vol. 46, no. 5, pp. 599–614, 2012. View at: Publisher Site  Google Scholar
 F. Glover, “Future paths for integer programming and links to artificial intelligence,” Computers & Operations Research, vol. 13, no. 5, pp. 533–549, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 F. Glover, T. Ye, A. P. Punnen, and G. Kochenberger, “Integrating tabu search and VLSN search to develop enhanced algorithms: a case study using bipartite boolean quadratic programs,” European Journal of Operational Research, vol. 241, no. 3, pp. 697–707, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Xue, C. Zhang, W.H. Lin, L. Miao, and P. Yang, “A tabu search heuristic for the local container drayage problem under a new operation mode,” Transportation Research Part E: Logistics and Transportation Review, vol. 62, pp. 136–150, 2014. View at: Publisher Site  Google Scholar
 J. Paquette, J.F. Cordeau, G. Laporte, and M. M. B. Pascoal, “Combining multicriteria analysis and tabu search for dialaride problems,” Transportation Research Part B: Methodological, vol. 52, pp. 1–16, 2013. View at: Publisher Site  Google Scholar
 F. Ruisanchez, L. Dell'Olio, and A. Ibeas, “Design of a tabu search algorithm for assigning optimal bus sizes and frequencies in urban transport services,” Journal of Advanced Transportation, vol. 46, no. 4, pp. 366–377, 2012. View at: Publisher Site  Google Scholar
 O. Shahvari and R. Logendran, “An enhanced tabu search algorithm to minimize a bicriteria objective in batching and scheduling problems on unrelatedparallel machines with desired lower bounds on batch sizes,” Computers & Operations Research, vol. 77, pp. 154–176, 2017. View at: Publisher Site  Google Scholar
 L. Kang and X. Zhu, “A simulated annealing algorithm for first train transfer problem in urban railway networks,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 40, no. 1, pp. 419–435, 2016. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2017 Shuo Yang 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.