Research Article  Open Access
Metro Timetabling for TimeVarying Passenger Demand and Congestion at Stations
Abstract
For the train timetabling problem (TTP) in a metro system, the operatororiented and passengeroriented objectives are both important but partly conflicting. This paper aims to minimize both objectives by considering frequency (in the line planning stage) and train cost (in the vehicle scheduling stage). Timevarying passenger demand and train capacity are considered in a nonsmooth, nonconvex programming model, which is transformed into a mixed integer programming model with a discrete timespace graph (DTSG). A novel dwell time determining process considering congestion at stations is proposed, which turns the dwell times into dependent variables. In the solution approach, we decompose the TTP into a subproblem for optimizing segment travel times (OST) and a subproblem for optimizing departure headways from the shunting yard (OH). Branchandbound and frequency determining algorithms are designed to solve OST. A novel rolling optimization algorithm is designed to solve OH. The numerical experiments include case studies on a short metro line and Beijing Metro Line 4, as well as sensitivity analyses. The results demonstrate the predictive ability of the model, verify the effectiveness and efficiency of the proposed approach, and provide practical insights for different scenarios, which can be used for decisionmaking support in daily operations.
1. Introduction
With increasing concerns about urban congestion and climate change, urban metro rail transportation receives increasing attention due to its high capacity, punctuality and sustainability. Metro timetabling is the problem of assigning precise utilization times for infrastructure resources to every train in the metro system [1]. In large cities, a metro system is often the busiest public transportation system. To meet travel demand and passenger satisfaction, a timetable should be as compact and flexible as possible. Meanwhile, since many metro systems need government subsidies to cover operating expenses, the planners focus more on operating cost than on passenger factors when designing a timetable [2]. This motivates some recent studies to combine these two conflicting objectives in the train timetabling problem [3, 4].
Two types of timetables are used in metro systems: cyclic and noncyclic. In a cyclic timetable, the departure times of trains are scheduled in equally spaced cycles (e.g., a half hour cycle). In a railway system, the cyclic timetables are assumed to be preferred by both operators and passengers because such timetables are easy to operate and remember. However, in a metro system, passengers usually arrive randomly at stations and wait for the next available train, without checking precise timetables beforehand. Besides, as pointed out by Robenek et al. [5], a cyclic timetable provides an inefficient operation as there is a mismatch between the supply (determined by the timetables) and the demand (characterized by the timevarying demand). On the other hand, a noncyclic timetable imposes no special rules on the departure times of trains [6]. The noncyclic timetables are more flexible regarding timevarying passenger demand, especially when the demand is large [2]. Thus, the noncyclic timetabling problem is worth exploring [2, 3, 7].
The planning process in public transportation can be generally split into three stages: line planning, timetabling, and vehicle scheduling [8]. The line planning stage determines the route/station layouts and frequency. Timetabling is based on the output of line planning, after which the vehicle scheduling can be designed. It is pointed out by Schöbel [8] that going through all these stages sequentially (i.e., independently) leads to unsatisfactory solutions. Since they are interrelated, important factors and objectives at the line planning stage and the vehicle scheduling stage should be accounted for in the timetabling problem.
This paper proposes a rolling optimization algorithm to obtain noncyclic timetables considering timevarying passenger demand and the effects of congestion at stations. It integrates the objectives in the line planning (frequency), timetabling (conflicting objectives including passenger wait/invehicle time and energy), and vehicle scheduling (train cost). The main contributions of this paper are as follows.
(1) This paper specifically models the dynamic evolution of passenger loads on trains at each station, by considering passenger arrival rates, limited train capacity and actual passenger alighting/boarding rates associated with congestion. This is an extension of existing studies [9]. Particularly, Niu and Zhou [10] and Niu et al. [11] propose the concept of “effective loading time” that represents the actual time interval within which the arriving passengers can board a train. The detailed boarding and alighting processes are not discussed in their model. Wang et al. [3] set a lower bound (which accounts for details including the number of boarding/alighting passengers and their boarding/alighting rates) for dwell times but consider the dwell times as variables to be optimized, rather than parameters that depend on arrival/departure times and passenger demand. Impacts of gradually increasing passenger demand are analyzed in Robenek et al. [2]; however, in modelling they do not consider the effects of passenger congestion.
(2) This paper integrates objectives from different planning stages in the train timetabling problem, for which solution approaches are scarce [8]. Particularly, Schöbel [8] considers the line planning, timetabling and vehicle scheduling in an integrated way, with passenger and costoriented objectives. However, since his eigenmodel is general and resulting approaches are generic, important objectives described in this paper are not considered in his work, that is, passenger wait time, dwell time, and train capacity. These objectives are included in [2, 5], whereas the frequency and dwell times are fixed in their studies. Besides, energy is not considered in [2, 5].
(3) This paper proposes a novel and effective solution approach, which decomposes the master timetabling model into two subproblems. The first subproblem is solved by branchandbound and frequency determining algorithms, and the second one is solved by a rolling optimization algorithm which optimizes the departure headways from the shunting yard (which are equivalent to arrival times of trains at their first station) and considers interdependent variables in each rolling step. The design of decision variables in this paper is similar to [2, 5], but the models and solution approaches are different. Compared to Wang et al. [3], which addresses the similar timetabling problem with traditional approaches such as sequential quadratic programming and a genetic algorithm, the approach proposed in this paper performs better computationally in terms of run time and solution quality.
The remainder of this paper is as follows. Section 2 reviews relevant studies on the train timetabling problem and summarizes the research gap in the literature. Section 3 first describes the problem, introduces formulation of different cost functions, and constructs the train timetabling model. Then the discrete timespace graph is proposed to transform the model into a mixed integer program (MIP). Section 4 discusses the solution approach for the optimization model. Section 5 presents different numerical experiments and sensitivity analyses on a short metro line and on Beijing Metro Line 4, to demonstrate the predictive ability of the model as well as the effectiveness and efficiency of the proposed approach. Section 6 summarizes the works done in this paper and presents the limitations and potential topics in future studies.
2. Literature Review
The train timetabling problem (TTP) aims to schedule trains to transport passengers (or goods) without conflicts, by specifying train arrival and departure times at stations. It is interrelated with other planning stages, that is, the line planning stage and the vehicle scheduling stage. We organize the related literature from the perspectives of operatororiented objectives, passengeroriented objectives, and their integration.
2.1. OperatorOriented Objectives
In the literature, operatororiented objectives can be found only for the noncyclic version of the TTP [2]. Most studies focus on the operation indicators, for example, energy [17, 18], train delay [19–21], and train travel time [22–24].
The total train delay and train travel time are usually considered in the TTP for railway networks, which aims to find a feasible timetable by minimizing the profit loss resulting from changes to the ideal/planned timetable. Energy is the focus and main cost for metro system operators. Most energy is consumed in train operations [25, 26], and thus obtaining an energyefficient train timetable receives considerable attention [17, 18, 27]. Note that energy is affected by the travel time in the segment, and the segment travel time is determined by the arrival and departure times at adjacent stations, that is, the timetable. Hence, the energyefficient train timetabling problem may be extended to include other objectives. For example, Yang et al. [23] consider both energy and train travel time as optimization objectives. Some recent studies also account for regenerative braking. Yang et al. [28] develop a scheduling approach to coordinate arrivals and departures of trains within the same electricity supply zones, thereby effectively recovering the regenerative energy. Generally, these studies focus solely on the timetabling stage, where frequency (in line planning) is fixed and train cost (invehicle scheduling) is not considered.
The efficient use of railway rolling stock (vehicle scheduling) is an important objective pursued by a railway agency or company because of intensive capital investment in rolling stock [29]. To this end, Lai et al. [29] develop an optimization model to improve the efficiency of rolling stock usage considering necessary regulations and practical constraints, where a hybrid heuristic process is designed to improve solution quality and efficiency. Haahr et al. [27] use CPLEX and a column and row generation approach to assign rolling stock units to timetable services in passenger railways, prepare daily schedules, and check their realtime applicability by testing different disruption scenarios. These vehicle scheduling studies require the trips as input data. Based on Schöbel [8], in a metro system, every trip describes the operation of a train between the start and end time of the line at its first and its last stations (given from the timetable). The objective function then aims at minimizing the number of trains and the costs for their movements.
Particularly, Schmid and Ehmke [30] and Schöbel [8] demonstrate that the integration of timetabling and vehicle scheduling is more beneficial than the sequential planning process. However, the studies on the integration of timetabling and vehicle scheduling usually adopt relatively general and generic models, and the TTP model is often simplified [31]. For example, Cadarso and Marin [32] focus on shunting operations and their timetable is calculated by readjusting frequencies, whereas specific details such as the dwell times and segment travel times are not considered. This motivates us to consider train cost in the timetabling stage, while focusing on the specified TTP formulation.
In addition, the conflictfree timetable is also pursued by the operators [33, 34]. It is worth mentioning that the discrete timespace graph proposed in Caprara et al. [6] is a directed multigraph in which nodes correspond to arrivals/departures at a certain station and at a given time stamp. This graph is widely used to formulate the TTP and derive different integer programming models that correspond to specific objectives. For example, Caprara et al. [6] use the graph to derive an integer linear programming (ILP) model with Lagrangian relaxation, which is embedded within a heuristic algorithm. Cacchiani et al. [12] use the graph to formulate an ILP with linear programming (LP) relaxation. Cacchiani et al. [35] use the timespace graphbased LP relaxation of an ILP to derive a dual bound in the TTP for a set of stations in an urban area interconnected by tracks, thus aiming to resolve the conflicts and evaluate the capacity saturation.
2.2. PassengerOriented Objectives
For a metro system, passengeroriented timetables that consider reliability and reduction of passenger time are most desirable [11]. In some studies, the passenger demand is considered stable [36]. Assuming that passengers prefer easily memorable timetables, such timetables are usually cyclic. They are designed on the basis of a period event scheduling problem (PESP), aiming to minimize passenger travel and wait time [15, 37] and to maximize network stability [38]. However, as mentioned above, metro passengers arrive randomly and do not remember the timetable. Besides, cyclic timetables are less flexible than noncyclic ones in accounting for the passenger demand [2].
Recent studies consider timevarying demands due to the growing concerns for service level and congestion at stations. Such timetables are noncyclic, with objectives of improving passenger satisfaction [3, 7, 10]. Particularly, Barrena et al. [7] propose two nonlinear programming (NP) formulations to generalize noncyclic train timetables on a single line, which are solved by a fast adaptive large neighborhood search metaheuristic. The objective is to minimize passenger wait times at stations. Niu et al. [11] construct mathematically rigorous and algorithmically tractable nonlinear mixed integer programming (MIP) models for both realtime scheduling and mediumterm planning applications, to jointly synchronize effective passenger loading time windows and train arrival and departure times at each station. Their work aims to minimize the total waiting times of passengers at stations. Robenek et al. [5] define a timetable as a set of departure times of every train from its origin station on every line and consider four attributes in passenger satisfaction: the invehicletime, the waiting time at transfers, the number of transfers, and the schedule passenger delay. They describe the hybrid timetables through additional constraints that are imposed on the original passenger centric train timetabling problem [2], so that the passengers would obtain the same level of service as a cyclic timetable, with more flexibility. A specifically defined simulated annealing heuristic is proposed to solve the problem. Luan et al. [16] integrate the TTP and preventive maintenance time slots (PMTSs) planning problem on a general railway network. They propose an MILP model to jointly optimize train routes, orders and passing times at each station, and PMTSs. The objective is to minimize the sum of the absolute arrival time deviations of real trains at destinations between the ideal and actual timetables. A Lagrangianbased label correcting algorithm is designed for solving the timedependent least cost path problem.
In addition, the delay management (DM) problem determines whether trains should wait for a delayed feeder train or should depart on time, while considering passengeroriented objectives [13]. Some recent studies integrate the macroscopic DM with the microscopic TTP [39]. Particularly, Dollevoet et al. [14] formulate an integer programming (IP) model and propose an iterative optimization approach to solve such a bilevel problem that the macroscopic level is the delay management and the microscopic level is the train timetabling. The optimization approach repeats a process that uses DM to find a solution and uses TTP to validate it, until a feasible solution is found. It should be noted that the graphbased models are adopted in their work, which simplify the formulations: the TTP is formulated with the alternative graph [40], and the DM is formulated with an eventactivity network [15].
2.3. Integrated Objectives
Since the operator and passengeroriented objectives are partly conflicting, focusing on only one side would yield undesirable solutions; for example, the best possible service for passengers may also be the most costly alternative for the operator [2]. Thus, many TTP studies integrate these objectives jointly [22, 41].
Wang et al. [3] aim to minimize passenger travel time and energy in their nonsmooth nonconvex programming (NSNCP) model. They propose an eventdriven model that involves arrival and departure events and passenger arrival rates change events and use nonlinear programming approaches and evolutionary algorithms to solve the problem.
Qi et al. [42] focus on a line planning problem with budget constraints and evaluate the station layout with train timetable indicators. The objectives include constructing cost for additional tracks, total travel time, and network capacity. A mixed integer linear programming (MILP) model is formulated to address the TTP, and then a bilevel programming (BP) model is designed to address the integrated problem. Both CPLEX and a local searchbased heuristic are developed to solve the models.
Robenek et al. [2] model the PCTTP as an MILP model with objectives of maximizing the operator’s profit while maintaining passenger satisfaction. The model uses the output of the line planning problem as the input of the train timetabling problem, to address the gap found between them. CPLEX is used to construct a Pareto Frontier of the contradictory objectives.
Li et al. [43] combine dynamic train regulation and passenger flow control to minimize the timetable and the headway deviations for metro lines, thereby reducing the operator profit loss and the passenger delay. The model predictive control and a quadratic programming algorithm are proposed to solve the problem, and stability analysis is conducted to verify the system convergence.
However, the above studies do not consider the impacts of frequency. Schöbel [44] reviews some approaches to model and solve the line planning problem, which demonstrate the impacts of frequency to both operator and passengeroriented objectives. Particularly, Samanta and Jha [45] consider different objectives in the line planning stage which include minimizing user cost, operator cost, and location cost and maximizing the ridership or the service coverage of the line. The objectives are separately formulated with NP models, and a GA is developed to optimize the objective functions subject to constraints on the number and spacing of stations. Fu et al. [46] propose a BP model to determine line frequencies and the stopping patterns for a mix of fast and slow train lines, with the objectives of minimizing passenger travel time and the number of transfers and maximizing train capacity occupancy. The heuristics solve the problem by combining additional conflicting demandsupply factors. Cascetta and Coppola [47] specify both the frequencybased assignment models and the timetablebased models in the demand forecasting, to analyze their differences in modal split estimates and flows on individual trains. The results show that when passenger demand is timevarying or the timetable is noncyclic, the difference is significant. Thus, the integration of frequency and timetabling is worth investigating.
2.4. Research Gap
Table 1 highlights some publications related to the TTP. All the graphbased models use integer programs, including the model in this paper. Although different objectives are considered, very few studies integrate the objectives in all three stages. Particularly, Wang et al. [3] only consider the objectives in one stage, that is, timetabling stage; Robenek et al. [2] integrate passenger satisfaction and operating profit (including train cost) but do not consider frequency and energy; Schöbel [8] takes into account objectives in all stages with an iterating framework; however, crucial factors such as passenger wait time, dwell time, and train capacity are missing. Kroon et al. [15] integrate the objectives in the TTP and vehicle scheduling but omit energy.

Apart from the studies focusing on the delay management, most TTPrelated studies optimize arrival and departure times of trains at each station; that is, the dwell times are optimized in the solution approach. However, the dwell times are heavily affected by the number of boarding and alighting passengers and their boarding and alighting rates, associated with congestion at stations. Thus, they should be considered as parameters that are dependent on other decision variables, that is, the departure times from the first station [2, 5] and the segment travel times. In this paper, we account for these characteristics and introduce specific equations to determine the dwell times in different scenarios. Such a process for determining dwell times has not been found in the literature.
Finally, heuristics and exact approaches are designed to find solutions based on the properties of the models. Particularly, rolling horizon [3] and iterative optimization [8, 14] are used to solve problems that have similar properties to our proposed model. However, the boundary constraints of the rolling horizon limit the solution flexibility of our model; the iterative optimization considers different stages with same importance, while this paper focuses only on the timetabling stage. Besides, the algorithms used in the rolling horizon and iterative optimization are relatively generic [3, 8]. Thus, a more specific algorithm is desirable for better solutions. In this regard, we develop a novel rolling optimization approach to solve the model, which combines the benefits of rolling horizon and iterative optimization.
3. Model Formulation
In this section, we present a mixed integer programming formulation for the TTP. It is a multiobjective optimization problem. The weightedsum scalarization of different objective functions in terms of generalized monetary cost is used to resolve the tradeoffs among them [8].
3.1. Problem Formulation
Table 2 summarizes the parameters and variables used in this paper. Some are also explained further in the text. The model’s simplifying assumptions are listed below and also explained in the problem description.

Assumption 1. Within one period, trains travel through each segment with the same speed profile; that is, the travel times for different trains in a segment are equal.
Assumption 2. We assume that the last train arrives at the first station at the end time of the period, to ensure that no latearriving passengers will be left behind.
Assumption 3. Passenger arrival rate should not exceed total passenger boarding rate . If , the queue of waiting passengers increases even if a train is dwelling at the station. In practice, passenger entrance flow is taken to avoid such a case.
. At most busy metro stations in China, passenger entrance flow control is adopted during peak hours to ensure operation safety, and thus the passenger arrival rate is a restricted number, which is expected not to exceed passenger boarding rate. Assumption 3 is introduced for consistency with practical experience.
Assumption 4. All passengers in train who have reached their destination are assumed to exit the train at its arrival time of station , that is, . All waiting passengers at station able to board train are assumed to board the train at its departure time, that is, . In other words, the number of invehicle passengers changes only at the arrival/departure times of train at stations.
This paper focuses on a bidirectional metro line with a single track, as shown in Figure 1. Properties of such a metro line are described in [48]. The timetabling problem is studied within a given period. The timevarying passenger demand is given as input, and a series of trains (whose frequency is to be optimized) are scheduled to transport the passengers.
The arrival times at the first station are decision variables that must be optimized. The departure times at the first station are determined by the dwell times. As mentioned earlier, this paper obtains the dwell times based on trains’ conditions when they arrive at stations. In this regard, optimizing arrival times at the first station is equivalent to optimizing departure times at the first station, as in [2, 5]. For the second station, the arrival times are determined by the segment travel time in the first segment. Practically, for the metro operators, uniform segment travel times are preferred within one period, because they not only make the system more stable and easier to manage but also simplify synchronized passenger transfers with connecting transit lines. Assumption 1 is introduced to consider this preference. Thus, within one period, the segment travel time in each segment needs to be optimized only once (for all trains within this period). The speed profile optimization is outside the scope of this paper. The departure times at the second station are also determined by the dwell time determining process.
Analogously, for subsequent train arrivals and departures, the arrival times are determined by the previous segment travel time, and the departure times are determined by the dwell times. Thus, within a period, the arrival times at the first station and the segment travel times determine the train timetable.
To represent the TTP, a directed multigraph called discrete timespace graph (DTSG) is proposed [6]. As shown in Figure 2, represents the shunting yard of the metro line and corresponds to the set of time stamps in which trains can arrive at (depart from) station . The nodes in are called arrival (departure) nodes. Unlike the original graph that is designed for a railway network [6, 12], this paper considers a singletrack metro line, and thus the constraints are different and implicitly imposed in the definition of the structure of the DTSG.
The continuous time is discretized into evenly spaced time intervals (1 second in this paper), and the period extends from 0 to . Note that the values of and are different. is the last time stamp of the passenger demand. The period mentioned above refers to the one extending from 0 to . Assumption 2 imposes that the last train departs from the first station at time . However, if the DTSG only considers the period from 0 to , the arrival/departures of the last train at subsequent stations are not considered, and the passengeroriented objectives are hard to formulate. Thus, the DTSG considers the period from 0 to . The relation between and is set as:where is the segment travel time of segment and is the maximum dwell time of station . This setting ensures that even the slowest train can be completely included in the DTSG.
The construction of the DTSG is noncyclic. The objective functions of passenger time, energy, and train cost are associated with the chosen arcs and are described below.
3.2. Passenger Time
3.2.1. Passenger Arriving, Alighting, and Boarding
(1) TimeVarying Demand Information. The passenger demand at each station can be expressed by , which is defined as
Passenger arrival rates at stations are , . In this paper, demands are regarded as deterministic values, which are given by . Thus, the number of arriving passengers at station with destination during time interval is
The number of arriving passengers at station within time interval is the sum of arriving passengers with different destinations:
(2) Determination of Dwell Arcs. Letting be the number of doors per train, the average number of alighting (boarding) passengers per door when train arrives at (departs from) station is . In addition, we denote the average number of boarding passengers per second (boarding rate) as and alighting rate as , which vary with passenger density in trains. For simplicity, we consider that boarding passengers wait until all alighting passengers get off the train; that is, the boarding process starts at time .
The dwell time is the sum of alighting and boarding times. The calculation of alighting time is solely based on the number of alighting passengers:
The calculation of boarding time is slightly more complex. Figure 3 shows the relations of train remaining capacity, waiting passengers, boarding rate, and max/min dwell times. The horizontal axis represents time, and the vertical axis represents the number of passengers. These four factors determine the boarding time for train at station , and in Figure 3 they are represented by remaining capacity of train at station : , a horizontal line; number of waiting passengers at station at time , a curve ; maximum allowable number of boarding passengers, a diagonal function with the slope of ; minimum and maximum dwell times at station , that is, and . In Figure 3, the diagonal intersects the horizontal line at time and intersects the curve at time . The curve intersects the horizontal line at time . Note that the time is discrete, and Assumption 3 ensures that the slope of always exceeds the slope of . Thus, we have
Practically, in an oversaturated scenario, a train departs from a station as soon as the number of intrain passengers reaches its capacity. Based on this strategy, (6) shows that the boarding passengers who arrive after exceed train capacity, and they will be left behind; (7) shows that before there are passengers waiting at station, and after all passengers board the train; (8) points that the train is full at , and passengers arriving after should wait for the next train.
The departure time of train at station is determined based on the relations among , , , , and . Based on Assumption 3, we consider two cases: and . Note that and should be considered in both cases. In case , we can further derive that (which is the solid curve in Figure 3), where the number of boarding passengers per second is determined by passenger arrival rate. Thus, the departure time is determined as . In case , we can further derive that (which is the dashed curve in Figure 3), where the number of boarding passengers per second is determined by the actual passenger boarding rate. Thus, the departure time is determined as . The departure time of train at station is expressed as
Finally, we consider the relations between alighting/boarding rates and congestion at stations. Here the congestion is described by the number of invehicle passengers. Practically, the greater the passenger density in a train is, the slower the passengers can get off/on the train. Here, we design a monotonically decreasing function with respect to the number of invehicle passengers, to approximate actual alighting/boarding rates according to the passenger density in trains. For simplicity, in this paper, the value of during the dwell process is fixed and is determined by the number of invehicle passengers at its arrival time. That is, and remain constant during the dwell process:where and are the ideal alighting and boarding rates in the uncongested condition.
(3) Passenger Alighting and Boarding. Once the dwell arcs are determined, the alighting and boarding times for train at station are determined. Based on the results of dwell arcs, three cases should be considered: ; and ; and . In case , waiting passengers board the train at rate . In case , actual passenger boarding rate is equal to the passenger arrival rate. In case , the actual boarding rate is . Thus, the number of boarding passengers for train at station is
Wang et al. [3] point out that, in a busy metro line, passengers with different destinations are well mixed at each station. Thus, the ratio between the numbers of boarding passengers with destination and of total boarding passengers is equivalent to the ratio between and :
The number of alighting passengers from train at station consists of all previously boarding passengers in train with destination :
There are no alighting passengers at the first station in each direction (station 1 and ).
3.2.2. Passenger Wait and Travel Time
(1) Passenger Wait Time. Total passenger wait time at a station can be calculated by summing up the products of the cumulative number of passengers and their time. Let be the index of last train that has departed from station before time , that is, . The cumulative number of passengers at station at time is the difference between cumulative numbers of arrival passengers and boarding passengers within the period :
Since the continuous time is discretized into 1 second intervals, the passenger wait time between times and is . Thus, the total passenger wait time is
(2) Passenger Travel Time. The number of intrain passengers in train at time is the difference between cumulative numbers of boarding passengers and alighting passengers. Let be the index of the last station that train has visited before time . We consider two cases: train is dwelling at station , that is, ; train is in the segment between station and , that is, . Based on Assumption 4, the boarding passengers at station are not considered in case but are included in case . Thus, the number of invehicle passengers in train at time is
Analogously, the total passenger travel time is
Finally, the total cost of passenger time is formulated as
The cost unit is yuan .
3.3. Energy
In general, for each segment , with given segment travel time and segment length , the speed profile determines energy use [49]. Thus, given a speed profile generating method, the energy per unit mass of different trains in segment is the same. We denote the energy per unit mass for one train in segment as . Its calculation usually considers the following train dynamic equations [48]:
The first two equations are train motion equations, where and , respectively, represent traveling speed and distance. The third one is the Davis resistance function where Davis parameters are , , and . The last one requires that every train stops at each station. An optimal train control strategy can be obtained based on (19), in which energy per unit mass is
However, in (20), the relation between energy and decision variable is unclear. It is also difficult to integrate the optimal train control with other objective functions [18]. Since the speed profile optimization is not considered, for simplicity, we use the method in [50] to generate practical speed profiles. Based on these speed profiles, we then adopt a linear piecewise curve to approximate the energy per unit mass of train in segment [48]:
We introduce here an estimation procedure for the values of and , which is based on mean squared error (MSE). First, we set several feasible values (e.g., ) of and . Then, we obtain optimal speed profiles based on (19) and (20), thereby obtaining the optimal energy with respect to and . Finally, we define the optimal energy per unit mass as , and the estimated energy per unit mass obtained with (21) as . The MSE is formulated as
The optimal values of and can be obtained when reaches its minimum and are found as follows:
With and obtained, the total energy can be calculated as
Thus, the total cost of energy consumption is
Note that when all tracks are straight and level, and suffice for estimating energy. However, if grades and curves are considered, we may have , where and are specified parameters to describe the characteristics in segment . The parameter estimation procedure is the same.
3.4. Train Cost
In this paper, the train cost is estimated with empirical equations with respect to the interest rate, train price, and average round trip time, where longterm requirements (e.g., the consideration of the highest peak during a year and demand evolution over the planning horizon) are not considered.
Train cost consists of train capital cost and operating cost (excluding energy). Train capital cost is related to the price and economic life of a train, as well as interest rate. Let be the annual cost per train, be the current price per train, be the interest rate, and be the lifespan of a train. The relations among these four factors can be expressed as
The capital cost per train hour can be then calculated by dividing the annual cost by the product of average working hours per day and working days per year
In addition, a train can be reassigned after it finishes a round trip and returns to the shunting yard, and the round trip time can be obtained as
Thus, the actual period for train running is . The estimated fleet size within period [] depends on the round trip time and the departure headway from the shunting yard (which is equivalent to the arrival time of train at the first station ):where the term represents the average round trip time and the term represents the average headway. The total cost of trains is formulated as
The first term represents the capital cost of trains for the actual running period . The second term represents the operating cost with frequency (trains) on a line with the length of the sum of segments (Km). Both decision variables (segment travel times and arrival times at the first station) affect the train cost.
3.5. Mixed Integer Programming Formulation
The optimization model aims to minimize the weightedsum of different cost functions, so that the conflicting objectives are generalized in terms of monetary cost. In this paper, the weights of different functions are reflected by their average costs, that is, , , , , and . For example, an operatororiented timetable (neglecting passenger satisfaction) can be obtained by setting and ; a passengeroriented timetable can be obtained by setting . The weights for integrated objectives can be set similarly, with larger cost coefficients for larger weights. The model is formulated as follows:
Minimize
Subject to
This model is a nonsmooth, nonconvex programming model, where the nonsmoothness is caused by the limited train capacity, and the nonconvexity is caused by the determination of dwell times [3]. Constraint (32) is equivalent to Assumption 2. Constraint (33) considers the turnaround from station to as traveling in segment , with fixed travel time . Constraint (34) ensures that the segment travel time should be optimized within a reasonable range . Constraint (35) imposes the minimum headway constraint for both departures and arrivals (based on Assumption 1). Constraint (36) requires that the number of intrain passengers must not exceed train capacity.
In Caprara et al. [6], the timespace graph is used to derive an integer linear programming model, thereby allowing a considerably faster solution. Analogously, this paper transforms the nonsmooth, nonconvex programming model into a mixed integer programming (MIP) model based on the DTSG. First, we construct two matrixes and , respectively, representing the selection of arrival and departure nodes in the DTSG:where if a train arrives at (departs from) station at time . As mentioned earlier, the value of is determined by , for , . The value of is determined by and the dwell arcs, that is, , for , , , . The value of is determined by and the travel time in segment , that is, , for , , and , . Thus, the decision variables determine and .
Based on the structures of and , as well as Assumption 4, we introduce two matrixes and , respectively, representing passenger alighting and boarding information (number and time):
The cumulative numbers of waiting passengers and invehicle passengers can be constructed based on and :
In (39), represents the total waiting passengers at station , where is the number of waiting passengers at the start of this period. According to the definition of cumulative passengers, we have . Thus, the total passenger wait time can be expressed by since . In (40), the initial number of invehicle passengers is zero, that is, . The number of invehicle passengers changes only at arrival/departure times, which is consistent with Assumption 4. Analogously, the total passenger travel time can be expressed as . Thus, the total cost for passenger time can be rewritten as
For energy cost, the energy per unit mass in segment can be rewritten as
The total energy can be rewritten aswhere if ; otherwise . The total energy cost can be rewritten correspondingly.
For train cost, the round trip time can be rewritten as
The departure headway from the shunting yard can be rewritten as
The total cost of trains can be rewritten correspondingly.
Finally, the last departing train constraint (32) is equivalent to
The turnaround constraint (33) is equivalent to
The segment travel time window constraint (34) is equivalent to
The minimum headway constraint (35) is equivalent to
The train capacity constraint (36) is equivalent to
The selection of arcs determines the arrival nodes and departure nodes , thereby determining every element in the optimization model. Thus, the model is transformed into an MIP. The next section discusses how its solution can be found.
4. Solution Approach
In this section, we decompose the TTP into two subproblems. The first subproblem optimizes the segment travel times (OST), aiming at minimizing the costs for energy and passenger segment travel times. A branchandbound algorithm and a frequency determining algorithm are introduced to solve the OST. The second subproblem optimizes departure headways from the shunting yard (OH), aiming to minimize passenger wait time, dwell time and train cost. A novel rolling optimization algorithm is designed to solve the OH.
4.1. Decomposition
The rationale of the decomposition is explained as follows. The objectives in the optimization model include passenger wait time (at their original stations), passenger travel time (segment travel times and station dwell times), energy, and train cost. Based on (43), the total energy is related to the energy per unit mass and the number of passengers in each segment. The energy per unit mass is solely determined by segment travel times . Since passengers board at their origins and alight at their destinations regardless of which train they take, the sum of invehicle passengers in one segment is determined by passenger demand. In this regard, the total energy can be expressed as
Therefore, the total energy cost depends only on segment travel times . The passenger segment travel times also depend only on segment travel times.
On the other hand, passenger wait time, dwell time, and train cost depend on the train arrival and departure times at stations. Based on Assumption 1, train arrival/departure times at stations are affected by the departure headways from the shunting yard. Thus, it is reasonable to use segment travel times as input and consider the costs for passenger wait time, dwell time, and trains as functions of the departure headways.
In general, the OST subproblem is formulated as follows:
Minimize
Subject to constraints (47) and (48).
It should be noted that since segment travel times are given inputs, minimizing the cost of passenger dwell time is equivalent to minimizing the cost of total travel time. Since we already obtain the intrain passenger matrix , there is no need to further derive dwell times at stations. Thus, the objectives in the OH subproblem include the total passenger time cost (both and ) instead of costs of passenger wait and dwell times. The OH subproblem is formulated as follows (with segment travel times as input):
Minimize
Subject to constraints (46), (49) and (50).
4.2. Solution Approach for OST
The solution approach for OST consists of branchandbound (B&B) and frequency determining algorithms. On one hand, the B&B algorithm needs frequency as input. On the other hand, the frequency determining algorithm needs segment travel times obtained by the B&B to evaluate the performance. Thus, we iteratively evaluate the performance with different frequencies (each frequency is used as input for B&B to find a solution of segment travel times) and choose the one with the best performance. The segment travel times are then determined correspondingly.
4.2.1. BranchandBound (B&B) Algorithm
The OST subproblem is an integer linear programming problem with decision variables (, for ). Here we use a generic B&B, which is widely used to solve linear programming problems [7]. Generally, a B&B needs an input model formulated as follows:where is the coefficient vector and is the decision variable vector. Both and have elements. and are coefficient matrixes, and and are constant vectors representing the bounds. Here is a matrix, and is a matrix. Correspondingly, has elements, and has one element. We define a vector with elements as
The decision variable vector can be represented as for , where is the th element in the vector . We then define for as the total number of passengers that travel through segment . Based on (52), can be represented as